Candidates for Balancing Selection in Leishmania donovani Complex Parasites

Ancestral aneuploidy and stable chromosomal duplication resulting in differential genome structure and gene expression control in trypanosomatid parasites

Aneuploidy is widely observed in both unicellular and multicellular eukaryotes, usually associated with adaptation to stress conditions. Chromosomal duplication stability is a tradeoff between the fitness cost of having unbalanced gene copies and the potential fitness gained from increased dosage of specific advantageous genes. Trypanosomatids, a family of protozoans that include species that cause neglected tropical diseases, are a relevant group to study aneuploidies. Their life cycle has several stressors that could select for different patterns of chromosomal duplications and/or losses, and their nearly universal use of polycistronic transcription increases their reliance on gene expansion/contraction, as well as post-transcriptional control as mechanisms for gene expression regulation. By evaluating the data from 866 isolates covering seven trypanosomatid genera, we have revealed that aneuploidy tolerance is an ancestral characteristic of trypanosomatids but has a reduced occurrence in a specific monophyletic clade that has undergone large genomic reorganization and chromosomal fusions. We have also identified an ancient chromosomal duplication that was maintained across these parasite's speciation, named collectively as the trypanosomatid ancestral supernumerary chromosome (TASC). TASC has most genes in the same coding strand, is expressed as a disomic chromosome (even having four copies), and has increased potential for functional variation, but it purges highly deleterious mutations more efficiently than other chromosomes. The evidence of stringent control over gene expression in this chromosome suggests that these parasites have adapted to mitigate the fitness cost associated with this ancient chromosomal duplication.

Parasite Genotype Is a Major Predictor of Mortality from Visceral Leishmaniasis

Visceral leishmaniasis (VL) is a potentially fatal disease caused mainly by Leishmania infantum in South America and Leishmania donovani in Asia and Africa. Disease outcomes have been associated with patient genotype, nutrition, age, sex, comorbidities, and coinfections. In this study, we examine the effects of parasite genetic variation on VL disease severity in Brazil. We collected and sequenced the genomes of 109 L. infantum isolates from patients in northeastern Brazil and retrieved matching patient clinical data from medical records, including mortality, sex, HIV coinfection, and laboratory data (creatinine, hemoglobin, and leukocyte and platelet counts). We identified genetic differences between parasite isolates, including single nucleotide polymorphisms (SNPs), small insertions/deletions (indels), and variations in genic, intergenic, and chromosome copy numbers (copy number variants [CNVs]). To describe associations between the parasite genotypes and clinical outcomes, we applied quantitative genetics methods of heritability and genome-wide association studies (GWAS), treating clinical outcomes as traits that may be influenced by parasite genotype. Multiple aspects of the genetic analysis indicate that parasite genotype affects clinical outcomes. We estimate that parasite genotype explains 83% chance of mortality (narrow-sense heritability [h²] = 0.83 ± 0.17) and has a significant relationship with patient sex (h² = 0.60 ± 0.27). Impacts of parasite genotype on other clinical traits are lower (h² ≤ 0.34). GWAS analysis identified multiple parasite genetic loci that were significantly associated with clinical outcomes; 17 CNVs were significantly associated with mortality, two with creatinine, and one with bacterial coinfection, jaundice, and HIV coinfection, and two SNPs/indels and six CNVs were associated with age, jaundice, HIV and bacterial coinfections, creatinine, and/or bleeding sites. Parasite genotype is an important factor in VL disease severity in Brazil. Our analysis indicates that specific genetic differences between parasites act as virulence factors, enhancing risks of severe disease and mortality. More detailed understanding of these virulence factors could be exploited for novel therapies. IMPORTANCE Multiple factors contribute to the risk of mortality from visceral leishmaniasis (VL), including, patient genotype, comorbidities, and nutrition. Many of these factors are influenced by socioeconomic biases. Our work suggests that the virulence of the infecting parasite is an important risk factor for mortality. We pinpoint some specific genomic markers that are associated with mortality, which can lead to a greater understanding of the molecular mechanisms that cause severe VL disease, to the identification of genetic markers for virulent parasites, and to the development of drug and vaccine therapies.

Assembly of a Large Collection of Maxicircle Sequences and Their Usefulness for Leishmania Taxonomy and Strain Typing

Parasites of medical importance, such as Leishmania and Trypanosoma, are characterized by the presence of thousands of circular DNA molecules forming a structure known as kinetoplast, within the mitochondria. The maxicircles, which are equivalent to the mitochondrial genome in other eukaryotes, have been proposed as a promising phylogenetic marker. Using whole-DNA sequencing data, it is also possible to assemble maxicircle sequences as shown here and in previous works. In this study, based on data available in public databases and using a bioinformatics workflow previously reported by our group, we assembled the complete coding region of the maxicircles for 26 prototypical strains of trypanosomatid species. Phylogenetic analysis based on this dataset resulted in a robust tree showing an accurate taxonomy of kinetoplastids, which was also able to discern between closely related Leishmania species that are usually difficult to discriminate by classical methodologies. In addition, we provide a dataset of the maxicircle sequences of 60 Leishmania infantum field isolates from America, Western Europe, North Africa, and Eastern Europe. In agreement with previous studies, our data indicate that L. infantum parasites from Brazil are highly homogeneous and closely related to European strains, which were transferred there during the discovery of America. However, this study showed the existence of different L. infantum populations/clades within the Mediterranean region. A maxicircle signature for each clade has been established. Interestingly, two L. infantum clades were found coexisting in the same region of Spain, one similar to the American strains, represented by the Spanish JPCM5 reference strain, and the other, named “non-JPC like”, may be related to an important leishmaniasis outbreak that occurred in Madrid a few years ago. In conclusion, the maxicircle sequence emerges as a robust molecular marker for phylogenetic analysis and species typing within the kinetoplastids, which also has the potential to discriminate intraspecific variability.

Tissue Specific Dual RNA-Seq Defines Host-Parasite Interplay in Murine Visceral Leishmaniasis Caused by Leishmania donovani and Leishmania infantum

Visceral leishmaniasis is associated with hepato-splenomegaly and altered immune and hematological parameters in both preclinical animal models and humans. We studied mouse experimental visceral leishmaniasis caused by Leishmania infantum and Leishmania donovani in BALB/c mice using dual RNA-seq to investigate the transcriptional response of host and parasite in liver and spleen. We identified only 4 species-specific parasite expressed genes (SSPEGs; log2FC >1, FDR <0.05) in the infected spleen, and none in the infected liver. For the host transcriptome, we found 789 differentially expressed genes (DEGs; log2FC >1, FDR <0.05) in the spleen that were common to both infections, with IFNγ signaling and complement and coagulation cascade pathways highly enriched, and an additional 286 and 186 DEGs that were selective to L. donovani and L. infantum infection, respectively. Among those, there were network interactions between genes of amino acid metabolism and PPAR signaling in L. donovani infection and increased IL1β and positive regulation of fatty acid transport in L. infantum infection, although no pathway enrichment was observed. In the liver, there were 1,939 DEGs in mice infected with either L. infantum or L. donovani in comparison to uninfected mice, and the most enriched pathways were IFNγ signaling, neutrophil mediated immunity, complement and coagulation, cytokine-chemokine responses, and hemostasis. Additionally, 221 DEGs were selective in L. donovani and 429 DEGs in L. infantum infections. These data show that the host response for these two visceral leishmaniasis infection models is broadly similar, and ∼10% of host DEGs vary in infections with either parasite species. IMPORTANCE Visceral leishmaniasis (VL) is caused by two species of Leishmania parasites, L. donovani in the Old World and L. infantum in the New World and countries bordering the Mediterranean. Although cardinal features such as hepato-splenomegaly and alterations in blood and immune function are evident, clinical presentation may vary by geography, with for example severe bleeding often associated with VL in Brazil. Although animal models of both L. donovani and L. infantum have been widely used to study disease pathogenesis, a direct side-by-side comparison of how these parasites species impact the infected host and/or how they might respond to the stresses of mammalian infection has not been previously reported. Identifying common and distinct pathways to pathogenesis will be important to ensure that new therapeutic or prophylactic approaches will be applicable across all forms of VL.