Dataset Information

Leishmania braziliensis causing human disease in Northeast Brazil presents loci with genotypes in long-term equilibrium.

ABSTRACT:

Background

Leishmaniases are neglected tropical diseases that inflict great burden to poor areas of the globe. Intense research has aimed to identify parasite genetic signatures predictive of infection outcomes. Consistency of diagnostic tools based on these markers would greatly benefit from accurate understanding of Leishmania spp. population genetics. We explored two chromosomal loci to characterize a population of L. braziliensis causing human disease in Northeast Brazil.

Methodology/principal findings

Two temporally distinct samples of L. braziliensis were obtained from patients attending the leishmaniasis clinic at the village of Corte de Pedra: (2008-2011) primary sample, N = 120; (1999-2001) validation sample, N = 35. Parasites were genotyped by Sanger's sequencing of two 600 base pairs loci starting at nucleotide positions 3,074 and 425,451 of chromosomes 24 and 28, respectively. Genotypes based on haplotypes of biallelic positions in each locus were tested for several population genetic parameters as well as for geographic clustering within the region. Ample geographic overlap of genotypes at the two loci was observed as indicated by non-significant Cusick and Edward's comparisons. No linkage disequilibrium was detected among combinations of haplotypes for both parasite samples. Homozygous and heterozygous genotypes displayed Hardy-Weinberg equilibrium (HWE) at both loci in the two samples when straight observed and expected counts were compared by Chi-square (p>0.5). However, Bayesian statistics using one million Monte-Carlo randomizations disclosed a less robust HWE for chromosome 24 genotypes, particularly in the primary sample (p = 0.04). Fixation indices (Fst) were consistently lower than 0.05 among individuals of the two samples at both tested loci, and no intra-populational structuralization could be detected using STRUCTURE software.

Conclusions/significance

These findings suggest that L. braziliensis can maintain stable populations in foci of human leishmaniasis and are capable of robust genetic recombination possibly due to events of sexual reproduction during the parasite's lifecycle.

SUBMITTER: Silva JA

PROVIDER: S-EPMC9239440 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundThe Leishmania (Viannia) braziliensis complex is responsible for most cases of New World tegumentary leishmaniasis. This complex includes two closely related species but with different geographic distribution and disease phenotypes, L. (V.) peruviana and L. (V.) braziliensis. However, the genetic basis of these differences is not well understood and the status of L. (V.) peruviana as distinct species has been questioned by some. Here we sequenced the genomes of two L. (V.) peruviana isolates (LEM1537 and PAB-4377) using Illumina high throughput sequencing and performed comparative analyses against the L. (V.) braziliensis M2904 reference genome. Comparisons were focused on the detection of Single Nucleotide Polymorphisms (SNPs), insertions and deletions (INDELs), aneuploidy and gene copy number variations.ResultsWe found 94,070 variants shared by both L. (V.) peruviana isolates (144,079 in PAB-4377 and 136,946 in LEM1537) against the L. (V.) braziliensis M2904 reference genome while only 26,853 variants separated both L. (V.) peruviana genomes. Analysis in coding sequences detected 26,750 SNPs and 1,513 indels shared by both L. (V.) peruviana isolates against L. (V.) braziliensis M2904 and revealed two L. (V.) braziliensis pseudogenes that are likely to have coding potential in L. (V.) peruviana. Chromosomal read density and allele frequency profiling showed a heterogeneous pattern of aneuploidy with an overall disomic tendency in both L. (V.) peruviana isolates, in contrast with a trisomic pattern in the L. (V.) braziliensis M2904 reference. Read depth analysis allowed us to detect more than 368 gene expansions and 14 expanded gene arrays in L. (V.) peruviana, and the likely absence of expanded amastin gene arrays.ConclusionsThe greater numbers of interspecific SNP/indel differences between L. (V.) peruviana and L. (V.) braziliensis and the presence of different gene and chromosome copy number variations support the classification of both organisms as closely related but distinct species. The extensive nucleotide polymorphisms and differences in gene and chromosome copy numbers in L. (V.) peruviana suggests the possibility that these may contribute to some of the unique features of its biology, including a lower pathology and lack of mucosal development.