Persistent T cell unresponsiveness associated with chronic visceral leishmaniasis in HIV-coinfected patients.

A large proportion of HIV-coinfected visceral leishmaniasis (VL-HIV) patients exhibit chronic disease with frequent VL recurrence. However, knowledge on immunological determinants underlying the disease course is scarce. We longitudinally profiled the circulatory cellular immunity of an Ethiopian HIV cohort that included VL developers. We show that chronic VL-HIV patients exhibit high and persistent levels of TIGIT and PD-1 on CD8+/CD8- T cells, in addition to a lower frequency of IFN-γ+ TIGIT- CD8+/CD8- T cells, suggestive of impaired T cell functionality. At single T cell transcriptome and clonal resolution, the patients show CD4+ T cell anergy, characterised by a lack of T cell activation and lymphoproliferative response. These findings suggest that PD-1 and TIGIT play a pivotal role in VL-HIV chronicity, and may be further explored for patient risk stratification. Our findings provide a strong rationale for adjunctive immunotherapy for the treatment of chronic VL-HIV patients to break the recurrent disease cycle.

Diversity and dissemination of viruses in pathogenic protozoa.

Viruses are the most abundant biological entities on Earth and play a significant role in the evolution of many organisms and ecosystems. In pathogenic protozoa, the presence of viruses has been linked to an increased risk of treatment failure and severe clinical outcome. Here, we studied the molecular epidemiology of the zoonotic disease cutaneous leishmaniasis in Peru and Bolivia through a joint evolutionary analysis of Leishmania braziliensis and their dsRNA Leishmania virus 1. We show that parasite populations circulate in tropical rainforests and are associated with single viral lineages that appear in low prevalence. In contrast, groups of hybrid parasites are geographically and ecologically more dispersed and associated with an increased prevalence, diversity and spread of viruses. Our results suggest that parasite gene flow and hybridization increased the frequency of parasite-virus symbioses, a process that may change the epidemiology of leishmaniasis in the region.

The adaptive roles of aneuploidy and polyclonality in Leishmania in response to environmental stress.

Aneuploidy is generally considered harmful, but in some microorganisms, it can act as an adaptive mechanism against environmental stress. Here, we use Leishmania-a protozoan parasite with remarkable genome plasticity-to study the early steps of aneuploidy evolution under high drug pressure (using antimony or miltefosine as stressors). By combining single-cell genomics, lineage tracing with cellular barcodes, and longitudinal genome characterization, we reveal that aneuploidy changes under antimony pressure result from polyclonal selection of pre-existing karyotypes, complemented by further and rapid de novo alterations in chromosome copy number along evolution. In the case of miltefosine, early parasite adaptation is associated with independent point mutations in a miltefosine transporter gene, while aneuploidy changes only emerge later, upon exposure to increased drug levels. Therefore, polyclonality and genome plasticity are hallmarks of parasite adaptation, but the scenario of aneuploidy dynamics depends on the nature and strength of the environmental stress as well as on the existence of other pre-adaptive mechanisms.

Genome diversity of Leishmania aethiopica.

Leishmania aethiopica is a zoonotic Old World parasite transmitted by Phlebotomine sand flies and causing cutaneous leishmaniasis in Ethiopia and Kenya. Despite a range of clinical manifestations and a high prevalence of treatment failure, L. aethiopica is one of the most neglected species of the Leishmania genus in terms of scientific attention. Here, we explored the genome diversity of L. aethiopica by analyzing the genomes of twenty isolates from Ethiopia. Phylogenomic analyses identified two strains as interspecific hybrids involving L. aethiopica as one parent and L. donovani and L. tropica respectively as the other parent. High levels of genome-wide heterozygosity suggest that these two hybrids are equivalent to F1 progeny that propagated mitotically since the initial hybridization event. Analyses of allelic read depths further revealed that the L. aethiopica - L. tropica hybrid was diploid and the L. aethiopica - L. donovani hybrid was triploid, as has been described for other interspecific Leishmania hybrids. When focusing on L. aethiopica, we show that this species is genetically highly diverse and consists of both asexually evolving strains and groups of recombining parasites. A remarkable observation is that some L. aethiopica strains showed an extensive loss of heterozygosity across large regions of the nuclear genome, which likely arose from gene conversion/mitotic recombination. Hence, our prospection of L. aethiopica genomics revealed new insights into the genomic consequences of both meiotic and mitotic recombination in Leishmania.

Genome Analysis of Triploid Hybrid Leishmania Parasite from the Neotropics.

We discovered a hybrid Leishmania parasite in Costa Rica that is genetically similar to hybrids from Panama. Genome analyses demonstrated the hybrid is triploid and identified L. braziliensis and L. guyanensis-related strains as parents. Our findings highlight the existence of poorly sampled Leishmania (Viannia) variants infectious to humans.

Parasite hybridization promotes spreading of endosymbiotic viruses.

Viruses are the most abundant biological entities on Earth and play a significant role in the evolution of many organisms and ecosystems. In pathogenic protozoa, the presence of endosymbiotic viruses has been linked to an increased risk of treatment failure and severe clinical outcome. Here, we studied the molecular epidemiology of the zoonotic disease cutaneous leishmaniasis in Peru and Bolivia through a joint evolutionary analysis of Leishmania braziliensis parasites and their endosymbiotic Leishmania RNA virus. We show that parasite populations circulate in isolated pockets of suitable habitat and are associated with single viral lineages that appear in low prevalence. In contrast, groups of hybrid parasites were geographically and ecologically dispersed, and commonly infected from a pool of genetically diverse viruses. Our results suggest that parasite hybridization, likely due to increased human migration and ecological perturbations, increased the frequency of endosymbiotic interactions known to play a key role in disease severity.

Transcriptional Shift and Metabolic Adaptations during Leishmania Quiescence Using Stationary Phase and Drug Pressure as Models.

Microorganisms can adopt a quiescent physiological condition which acts as a survival strategy under unfavorable conditions. Quiescent cells are characterized by slow or non-proliferation and a deep downregulation of processes related to biosynthesis. Although quiescence has been described mostly in bacteria, this survival skill is widespread, including in eukaryotic microorganisms. In Leishmania, a digenetic parasitic protozoan that causes a major infectious disease, quiescence has been demonstrated, but the molecular and metabolic features enabling its maintenance are unknown. Here, we quantified the transcriptome and metabolome of Leishmania promastigotes and amastigotes where quiescence was induced in vitro either, through drug pressure or by stationary phase. Quiescent cells have a global and coordinated reduction in overall transcription, with levels dropping to as low as 0.4% of those in proliferating cells. However, a subset of transcripts did not follow this trend and were relatively upregulated in quiescent populations, including those encoding membrane components, such as amastins and GP63, or processes like autophagy. The metabolome followed a similar trend of overall downregulation albeit to a lesser magnitude than the transcriptome. It is noteworthy that among the commonly upregulated metabolites were those involved in carbon sources as an alternative to glucose. This first integrated two omics layers afford novel insight into cell regulation and show commonly modulated features across stimuli and stages.

Genomic and Phenotypic Characterization of Experimentally Selected Resistant Leishmania donovani Reveals a Role for Dynamin-1-Like Protein in the Mechanism of Resistance to a Novel Antileishmanial Compound.

The implementation of prospective drug resistance (DR) studies in the research-and-development (R&D) pipeline is a common practice for many infectious diseases but not for neglected tropical diseases (NTDs). Here, we explored and demonstrated the importance of this approach using as paradigms Leishmania donovani, the etiological agent of visceral leishmaniasis (VL), and TCMDC-143345, a promising compound of the GlaxoSmithKline (GSK) "Leishbox" to treat VL. We experimentally selected resistance to TCMDC-143345 in vitro and characterized resistant parasites at the genomic and phenotypic levels. We found that it took more time to develop resistance to TCMDC-143345 than to other drugs in clinical use and that there was no cross-resistance to these drugs, suggesting a new and unique mechanism. By whole-genome sequencing, we found two mutations in the gene encoding the L. donovani dynamin-1-like protein (LdoDLP1) that were fixed at the highest drug pressure. Through phylogenetic analysis, we identified LdoDLP1 as a family member of the dynamin-related proteins, a group of proteins that impacts the shapes of biological membranes by mediating fusion and fission events, with a putative role in mitochondrial fission. We found that L. donovani lines genetically engineered to harbor the two identified LdoDLP1 mutations were resistant to TCMDC-143345 and displayed altered mitochondrial properties. By homology modeling, we showed how the two LdoDLP1 mutations may influence protein structure and function. Taken together, our data reveal a clear involvement of LdoDLP1 in the adaptation/reduced susceptibility of L. donovani to TCMDC-143345. IMPORTANCE Humans and their pathogens are continuously locked in a molecular arms race during which the eventual emergence of pathogen drug resistance (DR) seems inevitable. For neglected tropical diseases (NTDs), DR is generally studied retrospectively once it has already been established in clinical settings. We previously recommended to keep one step ahead in the host-pathogen arms race and implement prospective DR studies in the R&D pipeline, a common practice for many infectious diseases but not for NTDs. Here, using Leishmania donovani, the etiological agent of visceral leishmaniasis (VL), and TCMDC-143345, a promising compound of the GSK Leishbox to treat VL, as paradigms, we experimentally selected resistance to the compound and proceeded to genomic and phenotypic characterization of DR parasites. The results gathered in the present study suggest a new DR mechanism involving the L. donovani dynamin-1-like protein (LdoDLP1) and demonstrate the practical relevance of prospective DR studies.

High throughput single-cell genome sequencing gives insights into the generation and evolution of mosaic aneuploidy in Leishmania donovani.

Leishmania, a unicellular eukaryotic parasite, is a unique model for aneuploidy and cellular heterogeneity, along with their potential role in adaptation to environmental stresses. Somy variation within clonal populations was previously explored in a small subset of chromosomes using fluorescence hybridization methods. This phenomenon, termed mosaic aneuploidy (MA), might have important evolutionary and functional implications but remains under-explored due to technological limitations. Here, we applied and validated a high throughput single-cell genome sequencing method to study for the first time the extent and dynamics of whole karyotype heterogeneity in two clonal populations of Leishmania promastigotes representing different stages of MA evolution in vitro. We found that drastic changes in karyotypes quickly emerge in a population stemming from an almost euploid founder cell. This possibly involves polyploidization/hybridization at an early stage of population expansion, followed by assorted ploidy reduction. During further stages of expansion, MA increases by moderate and gradual karyotypic alterations, affecting a defined subset of chromosomes. Our data provide the first complete characterization of MA in Leishmania and pave the way for further functional studies.

Evaluation of whole genome amplification and bioinformatic methods for the characterization of Leishmania genomes at a single cell level.

Here, we report a pilot study paving the way for further single cell genomics studies in Leishmania. First, the performances of two commercially available kits for Whole Genome Amplification (WGA), PicoPLEX and RepliG were compared on small amounts of Leishmania donovani DNA, testing their ability to preserve specific genetic variations, including aneuploidy levels and SNPs. We show here that the choice of WGA method should be determined by the planned downstream genetic analysis, PicoPLEX and RepliG performing better for aneuploidy and SNP calling, respectively. This comparison allowed us to evaluate and optimize corresponding bio-informatic methods. As PicoPLEX was shown to be the preferred method for studying single cell aneuploidy, this method was applied in a second step, on single cells of L. braziliensis, which were sorted by fluorescence activated cell sorting (FACS). Even sequencing depth was achieved in 28 single cells, allowing accurate somy estimation. A dominant karyotype with three aneuploid chromosomes was observed in 25 cells, while two different minor karyotypes were observed in the other cells. Our method thus allowed the detection of aneuploidy mosaicism, and provides a solid basis which can be further refined to concur with higher-throughput single cell genomic methods.

Ecological divergence and hybridization of Neotropical Leishmania parasites.

The tropical Andes are an important natural laboratory to understand speciation in many taxa. Here we examined the evolutionary history of parasites of the Leishmania braziliensis species complex based on whole-genome sequencing of 67 isolates from 47 localities in Peru. We first show the origin of Andean Leishmania as a clade of near-clonal lineages that diverged from admixed Amazonian ancestors, accompanied by a significant reduction in genome diversity and large structural variations implicated in host-parasite interactions. Within the Andean species, patterns of population structure were strongly associated with biogeographical origin. Molecular clock and ecological niche modeling suggested that the history of diversification of the Andean lineages is limited to the Late Pleistocene and intimately associated with habitat contractions driven by climate change. These results suggest that changes in forestation over the past 150,000 y have influenced speciation and diversity of these Neotropical parasites. Second, genome-scale analyses provided evidence of meiotic-like recombination between Andean and Amazonian Leishmania species, resulting in full-genome hybrids. The mitochondrial genome of these hybrids consisted of homogeneous uniparental maxicircles, but minicircles originated from both parental species. We further show that mitochondrial minicircles-but not maxicircles-show a similar evolutionary pattern to the nuclear genome, suggesting that compatibility between nuclear-encoded mitochondrial genes and minicircle-encoded guide RNA genes is essential to maintain efficient respiration. By comparing full nuclear and mitochondrial genome ancestries, our data expand our appreciation on the genetic consequences of diversification and hybridization in parasitic protozoa.

Genomes of Leishmania parasites directly sequenced from patients with visceral leishmaniasis in the Indian subcontinent.

Whole genome sequencing (WGS) is increasingly used for molecular diagnosis and epidemiology of infectious diseases. Current Leishmania genomic studies rely on DNA extracted from cultured parasites, which might introduce sampling and biological biases into the subsequent analyses. Up to now, direct analysis of Leishmania genome in clinical samples is hampered by high levels of human DNA and large variation in parasite load in clinical samples. Here, we present a method, based on target enrichment of Leishmania donovani DNA with Agilent SureSelect technology, that allows the analysis of Leishmania genomes directly in clinical samples. We validated our protocol with a set of artificially mixed samples, followed by the analysis of 63 clinical samples (bone marrow or spleen aspirates) from visceral leishmaniasis patients in Nepal. We were able to identify genotypes using a set of diagnostic SNPs in almost all of these samples (97%) and access comprehensive genome-wide information in most (83%). This allowed us to perform phylogenomic analysis, assess chromosome copy number and identify large copy number variants (CNVs). Pairwise comparisons between the parasite genomes in clinical samples and derived in vitro cultured promastigotes showed a lower aneuploidy in amastigotes as well as genomic differences, suggesting polyclonal infections in patients. Altogether our results underline the need for sequencing parasite genomes directly in the host samples.

Tracking of quiescence in Leishmania by quantifying the expression of GFP in the ribosomal DNA locus.

Under stressful conditions some microorganisms adopt a quiescent stage characterized by a reversible non or slow proliferative condition that allows their survival. This adaptation was only recently discovered in Leishmania. We developed an in vitro model and a biosensor to track quiescence at population and single cell levels. The biosensor is a GFP reporter gene integrated within the 18S rDNA locus, which allows monitoring the expression of 18S rRNA (rGFP expression). We showed that rGFP expression decreased significantly and rapidly during the transition from extracellular promastigotes to intracellular amastigotes and that it was coupled in vitro with a decrease in replication as measured by BrdU incorporation. rGFP expression was useful to track the reversibility of quiescence in live cells and showed for the first time the heterogeneity of physiological stages among the population of amastigotes in which shallow and deep quiescent stages may coexist. We also validated the use of rGFP expression as a biosensor in animal models of latent infection. Our models and biosensor should allow further characterization of quiescence at metabolic and molecular level.

Integrated genomic and metabolomic profiling of ISC1, an emerging Leishmania donovani population in the Indian subcontinent.

Leishmania donovani is the responsible agent for visceral leishmaniasis (VL) in the Indian subcontinent (ISC). The disease is lethal without treatment and causes 0.2 to 0.4 million cases each year. Recently, reports of VL in Nepalese hilly districts have increased as well as VL cases caused by L. donovani from the ISC1 genetic group, a new and emerging genotype. In this study, we perform for the first time an integrated, untargeted genomics and metabolomics approach to characterize ISC1, in comparison with the Core Group (CG), main population that drove the most recent outbreak of VL in the ISC. We show that the ISC1 population is very different from the CG, both at genome and metabolome levels. The genomic differences include SNPs, CNV and small indels in genes coding for known virulence factors, immunogens and surface proteins. Both genomic and metabolic approaches highlighted dissimilarities related to membrane lipids, the nucleotide salvage pathway and the urea cycle in ISC1 versus CG. Many of these pathways and molecules are important for the interaction with the host/extracellular environment. Altogether, our data predict major functional differences in ISC1 versus CG parasites, including virulence. Therefore, particular attention is required to monitor the fate of this emerging ISC1 population in the ISC, especially in a post-VL elimination context.

Macromolecular biosynthetic parameters and metabolic profile in different life stages of Leishmania braziliensis: Amastigotes as a functionally less active stage.

It was recently hypothesized that Leishmania amastigotes could constitute a semi-quiescent stage characterized by low replication and reduced metabolic activity. This concept developed with Leishmania (Leishmania) mexicana and Leishmania (Leishmania) major models might explain numerous clinical and sub-clinical features of Leishmania (Viannia) braziliensis infections, like reactivation of the disease, non-response to chemotherapy or asymptomatic infections. We compared here in vitro the proliferative capability of L. (V.) braziliensis amastigotes and promastigotes, assessed the expression of key molecular parameters and performed metabolomic analysis. We found that contrary to the highly proliferative promastigotes, amastigotes (axenic and intracellular) do not show evidence of extensive proliferation. In parallel, amastigotes showed a significant decrease of (i) the kDNA mini-circle abundance, (ii) the intracellular ATP level, (iii) the ribosomal components: rRNA subunits 18S and 28S α and ribosomal proteins RPS15 and RPL19, (iv) total RNA and protein levels. An untargeted metabolomic study identified clear differences between the different life stages: in comparison to logarithmic promastigotes, axenic amastigotes showed (a) a strong decrease of 14 essential and non-essential amino acids and eight metabolites involved in polyamine synthesis, (b) extensive changes in the phospholipids composition and (c) increased levels of several endogenous and exogenous sterols. Altogether, our results show that L. (V.) braziliensis amastigotes can show a phenotype with negligible rate of proliferation, a lower capacity of biosynthesis, a reduced bio-energetic level and a strongly altered metabolism. Our results pave the way for further exploration of quiescence among amastigotes of this species.

Phylogenetic analysis of the Trypanosoma genus based on the heat-shock protein 70 gene.

Trypanosome evolution was so far essentially studied on the basis of phylogenetic analyses of small subunit ribosomal RNA (SSU-rRNA) and glycosomal glyceraldehyde-3-phosphate dehydrogenase (gGAPDH) genes. We used for the first time the 70kDa heat-shock protein gene (hsp70) to investigate the phylogenetic relationships among 11 Trypanosoma species on the basis of 1380 nucleotides from 76 sequences corresponding to 65 strains. We also constructed a phylogeny based on combined datasets of SSU-rDNA, gGAPDH and hsp70 sequences. The obtained clusters can be correlated with the sections and subgenus classifications of mammal-infecting trypanosomes except for Trypanosoma theileri and Trypanosoma rangeli. Our analysis supports the classification of Trypanosoma species into clades rather than in sections and subgenera, some of which being polyphyletic. Nine clades were recognized: Trypanosoma carassi, Trypanosoma congolense, Trypanosoma cruzi, Trypanosoma grayi, Trypanosoma lewisi, T. rangeli, T. theileri, Trypanosoma vivax and Trypanozoon. These results are consistent with existing knowledge of the genus' phylogeny. Within the T. cruzi clade, three groups of T. cruzi discrete typing units could be clearly distinguished, corresponding to TcI, TcIII, and TcII+V+VI, while support for TcIV was lacking. Phylogenetic analyses based on hsp70 demonstrated that this molecular marker can be applied for discriminating most of the Trypanosoma species and clades.

Evolutionary genomics of epidemic visceral leishmaniasis in the Indian subcontinent.

Leishmania donovani causes visceral leishmaniasis (VL), the second most deadly vector-borne parasitic disease. A recent epidemic in the Indian subcontinent (ISC) caused up to 80% of global VL and over 30,000 deaths per year. Resistance against antimonial drugs has probably been a contributing factor in the persistence of this epidemic. Here we use whole genome sequences from 204 clinical isolates to track the evolution and epidemiology of L. donovani from the ISC. We identify independent radiations that have emerged since a bottleneck coincident with 1960s DDT spraying campaigns. A genetically distinct population frequently resistant to antimonials has a two base-pair insertion in the aquaglyceroporin gene LdAQP1 that prevents the transport of trivalent antimonials. We find evidence of genetic exchange between ISC populations, and show that the mutation in LdAQP1 has spread by recombination. Our results reveal the complexity of L. donovani evolution in the ISC in response to drug treatment.

Association of the Endobiont Double-Stranded RNA Virus LRV1 With Treatment Failure for Human Leishmaniasis Caused by Leishmania braziliensis in Peru and Bolivia.

Cutaneous and mucosal leishmaniasis, caused in South America by Leishmania braziliensis, is difficult to cure by chemotherapy (primarily pentavalent antimonials [Sb(V)]). Treatment failure does not correlate well with resistance in vitro, and the factors responsible for treatment failure in patients are not well understood. Many isolates of L. braziliensis (>25%) contain a double-stranded RNA virus named Leishmaniavirus 1 (LRV1), which has also been reported in Leishmania guyanensis, for which an association with increased pathology, metastasis, and parasite replication was found in murine models. Here we probed the relationship of LRV1 to drug treatment success and disease in 97 L. braziliensis-infected patients from Peru and Bolivia. In vitro cultures were established, parasites were typed as L. braziliensis, and the presence of LRV1 was determined by reverse transcription-polymerase chain reaction, followed by sequence analysis. LRV1 was associated significantly with an increased risk of treatment failure (odds ratio, 3.99; P = .04). There was no significant association with intrinsic Sb(V) resistance among parasites, suggesting that treatment failure arises from LRV1-mediated effects on host metabolism and/or parasite survival. The association of LRV1 with clinical drug treatment failure could serve to guide more-effective treatment of tegumentary disease caused by L. braziliensis.

Differentiation between Trypanosoma cruzi and Trypanosoma rangeli using heat-shock protein 70 polymorphisms.

OBJECTIVE: Differential diagnosis of infection with Trypanosoma cruzi or T. rangeli is relevant for epidemiological studies and clinical practice as both species infect humans, but only T. cruzi causes Chagas' disease. Their common antigen determinants complicate the distinction between both species, while current PCR assays used for differentiation show some drawbacks. We developed and validated a generic PCR discriminating the species by restriction fragment length polymorphism (RFLP) analysis and a duplex PCR specifically amplifying a differently sized fragment of both species. METHODS: The assays are based upon a partial region of the heat-shock protein 70 gene (hsp70). The analytical sensitivity and specificity were determined for both PCRs. The assays were analytically evaluated on a panel of six T. cruzi, one T. cruzi marinkellei and four T. rangeli strains, various other infectious pathogens, a panel of spiked samples of T. cruzi, and artificially mixed infections of T. cruzi and T. rangeli. Finally, the tools were applied on 36 additional isolates of Trypanosoma species. RESULTS: The detection limit of the PCRs was between 0.05 and 0.5 parasite genomes, and 1-10 parasites spiked in 200 µl blood. In artificial mixtures, PCR-RFLP picked up both species in ratios up to 10(2) and duplex PCR up to 10(4) . In the 36 isolates tested, both single and mixed infections were identified. All assays were shown to be specific. CONCLUSION: Our PCRs show high potential for the differential diagnosis of T. cruzi and T. rangeli, which in view of their sensitivity can aid in the confirmation of infection with these parasites in vectors, reservoirs and clinical samples.

Metabolic adaptations of Leishmania donovani in relation to differentiation, drug resistance, and drug pressure.

Antimonial (sodium stibogluconate, SSG) resistance and differentiation have been shown to be closely linked in Leishmania donovani, with SSG-resistant strains showing an increased capacity to generate infectious (metacyclic) forms. This is the first untargeted LC-MS metabolomics study which integrated both phenomena in one experimental design and provided insights into metabolic differences between three clinical L. donovani strains with a similar genetic background but different SSG-susceptibilities. We performed this analysis at different stages during promastigote growth and in the absence or presence of drug pressure. When comparing SSG-resistant and SSG-sensitive strains, a number of metabolic changes appeared to be constitutively present in all growth stages, pointing towards a clear link with SSG-resistance, whereas most metabolic changes were only detected in the stationary stage. These changes reflect the close intertwinement between SSG-resistance and an increased metacyclogenesis in resistant parasites. The metabolic changes suggest that SSG-resistant parasites have (i) an increased capacity for protection against oxidative stress; (ii) a higher fluidity of the plasma membrane; and (iii) a metabolic survival kit to better endure infection. These changes were even more pronounced in a resistant strain kept under Sb(III) drug pressure.