The evolution of amastin surface glycoproteins in trypanosomatid parasites

Genomic and epidemiological evidence for the emergence of a L. infantum/L. donovani hybrid with unusual epidemiology in northern Italy

Leishmania (L.) infantum is one of the main causative agents of animal and human leishmaniasis across many endemic areas in South America, Europe, North Africa, and Asia. Despite its clinical significance, little is known about the genetic diversity of L. infantum circulating in a given endemic area. Here, we investigate this important open question by applying a comparative genomics approach to seven L. infantum isolates from different hosts and Italian regions, including the northern part of the country (Emilia-Romagna, RER), Sicily, and Sardinia, as an initial attempt to explore the breadth of parasite genetic heterogeneity in Italy. Additionally, microsatellite analysis was carried out to compare the isolates from RER with other 70 L. infantum strains from the same region as well as 65 strains belonging to the L. donovani complex from other countries. We revealed important karyotypic instability and identified strain-specific changes in gene dosage, which affected important virulence factors such as amastins and surface antigen-like proteins. Single nucleotide polymorphism-based clustering analysis of these genomes together with over 80 publicly available L. infantum and L. donovani genomes placed the Italian isolates into three geographically distinct clusters within the Mediterranean basin and uncovered three isolates clustering with putative L. infantum/L. donovani hybrids isolated in Cyprus. As judged by microsatellite profiling, these hybrid isolates are representative of a sub-population of parasites circulating in northern Italy that preferentially infect humans but not dogs. Our results place Italy at the crossroads of L. infantum infection in the Mediterranean and call attention to the public health risk represented by the introduction of non-European Leishmania species.IMPORTANCEThis study closes important knowledge gaps with respect to Leishmania (L.) infantum genetic heterogeneity in a given endemic country, as exemplified here for Italy, and reveals genetic hybridization as a main cause for re-emerging human leishmaniasis in northern Italy. The observed high diversity of Leishmania parasites on the Italian peninsula suggests different geographical origins, with genomic adaptation to various ecologies affecting both pathogenicity and transmission potential. This is documented by the discovery of a putative L. infantum/L. donovani hybrid strain, which has been shown to preferentially infect humans but not dogs. Our results provide important information to health authorities, which need to consider the public health risk represented by the introduction of new Leishmania species into EU countries due to population displacement or travel from countries where exotic/allochthonous parasite species are endemic.

Immunogenicity and protective efficacy of tuzin protein as a vaccine candidate in Leishmania donovani-infected BALB/c mice

Visceral leishmaniasis (VL) is referred to as the most severe and fatal type of leishmaniasis basically caused by Leishmania donovani and L. infantum. The most effective method for preventing the spread of the disease is vaccination. Till today, there is no promising licensed vaccination for human VL. Hence, investigation for vaccines is necessary to enrich the therapeutic repertoire against leishmaniasis. Tuzin is a rare trans-membrane protein that has been reported in Trypanosoma cruzi with unknown function. However, tuzin is not characterized in Leishmania parasites. In this study, we for the first time demonstrated that tuzin protein was expressed in both stages (promastigote and amastigote) of L. donovani parasites. In-silico studies revealed that tuzin has potent antigenic properties. Therefore, we analyzed the immunogenicity of tuzin protein and immune response in BALB/c mice challenged with the L. donovani parasite. We observed that tuzin-vaccinated mice have significantly reduced parasite burden in the spleen and liver compared with the control. The number of granulomas in the liver was also significantly decreased compared with the control groups. We further measured the IgG2a antibody level, a marker of Th1 immune response in VL, which was significantly higher in the serum of immunized mice when compared with the control. Splenocytes stimulated with soluble Leishmania antigen (SLA) displayed a significant increase in NO and ROS levels compared with the control groups. Tuzin-immunized and parasite-challenged mice exhibit a notable rise in the IFN-γ/IL-10 ratio by significantly suppressing IL-10 expression level, an immunosuppressive cytokine that inhibits leishmanicidal immune function and encourages disease progression. In conclusion, tuzin immunizations substantially increase the protective immune response in L. donovani-challenged mice groups compared with control.

Genomes of Endotrypanum monterogeii from Panama and Zelonia costaricensis from Brazil: Expansion of Multigene Families in Leishmaniinae Parasites That Are Close Relatives of Leishmania spp

The Leishmaniinae subfamily of the Trypanosomatidae contains both genus Zelonia (monoxenous) and Endotrypanum (dixenous). They are amongst the nearest known relatives of Leishmania, which comprises many human pathogens widespread in the developing world. These closely related lineages are models for the genomic biology of monoxenous and dixenous parasites. Herein, we used comparative genomics to identify the orthologous groups (OGs) shared among 26 Leishmaniinae species to investigate gene family expansion/contraction and applied two phylogenomic approaches to confirm relationships within the subfamily. The Endotrypanum monterogeii and Zelonia costaricensis genomes were assembled, with sizes of 29.9 Mb and 38.0 Mb and 9.711 and 12.201 predicted protein-coding genes, respectively. The genome of E. monterogeii displayed a higher number of multicopy cell surface protein families, including glycoprotein 63 and glycoprotein 46, compared to Leishmania spp. The genome of Z. costaricensis presents expansions of BT1 and amino acid transporters and proteins containing leucine-rich repeat domains, as well as a loss of ABC-type transporters. In total, 415 and 85 lineage-specific OGs were identified in Z. costaricensis and E. monterogeii. The evolutionary relationships within the subfamily were confirmed using the supermatrix (3384 protein-coding genes) and supertree methods. Overall, this study showed new expansions of multigene families in monoxenous and dixenous parasites of the subfamily Leishmaniinae.

Modulation of Virulence Factors during Trypanosoma cruzi Differentiation

Chagas disease is a neglected tropical disease caused by Trypanosoma cruzi. This protozoan developed several mechanisms to infect, propagate, and survive in different hosts. The specific expression of proteins is responsible for morphological and metabolic changes in different parasite stages along the parasite life cycle. The virulence strategies at the cellular and molecular levels consist of molecules responsible for mediating resistance mechanisms to oxidative damage, cellular invasion, and immune evasion, performed mainly by surface proteins. Since parasite surface coat remodeling is crucial to invasion and infectivity, surface proteins are essential virulence elements. Understanding the factors involved in these processes improves the knowledge of parasite pathogenesis. Genome sequencing has opened the door to high-throughput technologies, allowing us to obtain a deeper understanding of gene reprogramming along the parasite life cycle and identify critical molecules for survival. This review therefore focuses on proteins regulated during differentiation into infective forms considered virulence factors and addresses the current known mechanisms acting in the modulation of gene expression, emphasizing mRNA signals, regulatory factors, and protein complexes.

Quantitative monitoring of experimental and human leishmaniasis employing amastigote-specific genes

The gold standard for diagnosis of leishmaniasis is the microscopic detection of amastigotes/Leishman Donovan (LD) bodies, but its moderate sensitivity necessitates the development of molecular approaches. This study aimed to quantify in experimental animal models and human leishmaniasis the expression of amastigote-specific virulence genes, A2 and amastin by droplet digital polymerase chain reaction (ddPCR). Total RNA was isolated from L. donovani-infected hamsters or murine peritoneal macrophages and lesional biopsies from patients with post kala-azar dermal leishmaniasis (PKDL). Following cDNA conversion, EvaGreen-based ddPCR was performed using specific primers for A2 or amastin and parasite load expressed in copies per μL. Assay was optimized and the specificity of amastigote-specific A2 and amastin was confirmed. In hepatic and splenic tissues of L. donovani-infected hamsters and peritoneal macrophages, ddPCR demonstrated a greater abundance of A2 than amastin. Treatment of L. donovani-infected peritoneal macrophages with conventional anti-leishmanials, miltefosine and amphotericin B translated into a dose-dependent reduction in copies per μL of A2 and amastin, and the extrapolated IC50 was comparable with results obtained by counting LD bodies in Giemsa-stained macrophages. Similarly, in dermal biopsies of patients with PKDL, A2 and amastin were detected. Overall, monitoring of A2 by ddPCR can be an objective measure of parasite burden and potentially adaptable into a high throughput approach necessary for drug development and monitoring disease progression when the causative species is L. donovani.

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.

Genomics of Trypanosomatidae: Where We Stand and What Needs to Be Done?

Trypanosomatids are easy to cultivate and they are (in many cases) amenable to genetic manipulation. Genome sequencing has become a standard tool routinely used in the study of these flagellates. In this review, we summarize the current state of the field and our vision of what needs to be done in order to achieve a more comprehensive picture of trypanosomatid evolution. This will also help to illuminate the lineage-specific proteins and pathways, which can be used as potential targets in treating diseases caused by these parasites.

AQPX-cluster aquaporins and aquaglyceroporins are asymmetrically distributed in trypanosomes

Major Intrinsic Proteins (MIPs) are membrane channels that permeate water and other small solutes. Some trypanosomatid MIPs mediate the uptake of antiparasitic compounds, placing them as potential drug targets. However, a thorough study of the diversity of these channels is still missing. Here we place trypanosomatid channels in the sequence-function space of the large MIP superfamily through a sequence similarity network. This analysis exposes that trypanosomatid aquaporins integrate a distant cluster from the currently defined MIP families, here named aquaporin X (AQPX). Our phylogenetic analyses reveal that trypanosomatid MIPs distribute exclusively between aquaglyceroporin (GLP) and AQPX, being the AQPX family expanded in the Metakinetoplastina common ancestor before the origin of the parasitic order Trypanosomatida. Synteny analysis shows how African trypanosomes specifically lost AQPXs, whereas American trypanosomes specifically lost GLPs. AQPXs diverge from already described MIPs on crucial residues. Together, our results expose the diversity of trypanosomatid MIPs and will aid further functional, structural, and physiological research needed to face the potentiality of the AQPXs as gateways for trypanocidal drugs.

Genome Analysis of Endotrypanum and Porcisia spp., Closest Phylogenetic Relatives of Leishmania, Highlights the Role of Amastins in Shaping Pathogenicity

While numerous genomes of Leishmania spp. have been sequenced and analyzed, an understanding of the evolutionary history of these organisms remains limited due to the unavailability of the sequence data for their closest known relatives, Endotrypanum and Porcisia spp., infecting sloths and porcupines. We have sequenced and analyzed genomes of three members of this clade in order to fill this gap. Their comparative analyses revealed only minute differences from Leishmaniamajor genome in terms of metabolic capacities. We also documented that the number of genes under positive selection on the Endotrypanum/Porcisia branch is rather small, with the flagellum-related group of genes being over-represented. Most significantly, the analysis of gene family evolution revealed a substantially reduced repertoire of surface proteins, such as amastins and biopterin transporters BT1 in the Endotrypanum/Porcisia species when compared to amastigote-dwelling Leishmania. This reduction was especially pronounced for δ-amastins, a subfamily of cell surface proteins crucial in the propagation of Leishmania amastigotes inside vertebrate macrophages and, apparently, dispensable for Endotrypanum/Porcisia, which do not infect such cells.

VLP-Based Vaccines as a Suitable Technology to Target Trypanosomatid Diseases

Research on vaccines against trypanosomatids, a family of protozoa that cause neglected tropical diseases, such as Chagas disease, leishmaniasis, and sleeping sickness, is a current need. Today, according to modern vaccinology, virus-like particle (VLP) technology is involved in many vaccines, including those undergoing studies related to COVID-19. The potential use of VLPs as vaccine adjuvants opens an opportunity for the use of protozoan antigens for the development of vaccines against diseases caused by Trypanosoma cruzi, Leishmania spp., and Trypanosoma brucei. In this context, it is important to consider the evasion mechanisms of these protozoa in the host and the antigens involved in the mechanisms of the parasite–host interaction. Thus, the immunostimulatory properties of VLPs can be part of an important strategy for the development and evaluation of new vaccines. This work aims to highlight the potential of VLPs as vaccine adjuvants for the development of immunity in complex diseases, specifically in the context of tropical diseases caused by trypanosomatids.

Genomic and Transcriptomic Analysis for Identification of Genes and Interlinked Pathways Mediating Artemisinin Resistance in Leishmania donovani

Current therapy for visceral leishmaniasis (VL), compromised by drug resistance, toxicity, and high cost, demands for more effective, safer, and low-cost drugs. Artemisinin has been found to be an effectual drug alternative in experimental models of leishmaniasis. Comparative genome and transcriptome analysis of in vitro-adapted artesunate-resistant (K133AS-R) and -sensitive wild-type (K133WT) Leishmania donovani parasites was carried out using next-generation sequencing and single-color DNA microarray technology, respectively, to identify genes and interlinked pathways contributing to drug resistance. Whole-genome sequence analysis of K133WT vs. K133AS-R parasites revealed substantial variation among the two and identified 240 single nucleotide polymorphisms (SNPs), 237 insertion deletions (InDels), 616 copy number variations (CNVs) (377 deletions and 239 duplications), and trisomy of chromosome 12 in K133AS-R parasites. Transcriptome analysis revealed differential expression of 208 genes (fold change ≥ 2) in K133AS-R parasites. Functional categorization and analysis of modulated genes of interlinked pathways pointed out plausible adaptations in K133AS-R parasites, such as (i) a dependency on lipid and amino acid metabolism for generating energy, (ii) reduced DNA and protein synthesis leading to parasites in the quiescence state, and (iii) active drug efflux. The upregulated expression of cathepsin-L like protease, amastin-like surface protein, and amino acid transporter and downregulated expression of the gene encoding ABCG2, pteridine receptor, adenylatecyclase-type receptor, phosphoaceylglucosamine mutase, and certain hypothetical proteins are concordant with genomic alterations suggesting their potential role in drug resistance. The study provided an understanding of the molecular basis linked to artemisinin resistance in Leishmania parasites, which may be advantageous for safeguarding this drug for future use.

Reductionist Pathways for Parasitism in Euglenozoans? Expanded Datasets Provide New Insights

The unicellular trypanosomatids belong to the phylum Euglenozoa and all known species are obligate parasites. Distinct lineages infect plants, invertebrates, and vertebrates, including humans. Genome data for marine diplonemids, together with freshwater euglenids and free-living kinetoplastids, the closest known nonparasitic relatives to trypanosomatids, recently became available. Robust phylogenetic reconstructions across Euglenozoa are now possible and place the results of parasite-focused studies into an evolutionary context. Here we discuss recent advances in identifying the factors shaping the evolution of Euglenozoa, focusing on ancestral features generally considered parasite-specific. Remarkably, most of these predate the transition(s) to parasitism, suggesting that the presence of certain preconditions makes a significant lifestyle change more likely.

Dual transcriptome analysis reveals differential gene expression modulation influenced by Leishmania arginase and host genetic background

The outcome of Leishmania infection is strongly influenced by the host's genetic background. BALB/c mice are susceptible to Leishmania infection, while C57BL/6 mice show discrete resistance. Central to the fate of the infection is the availability of l-arginine and the related metabolic processes in the host and parasite. Depending on l-arginine availability, nitric oxide synthase 2 (NOS2) of the host cell produces nitric oxide (NO) controlling the parasite growth. On the other hand, Leishmania can also use host l-arginine for the production of polyamines through its own arginase activity, thus favouring parasite replication. Considering RNA-seq data, we analysed the dual modulation of host and parasite gene expression of BALB/c or C57BL/6 mouse bone marrow-derived macrophages (BMDMs) after 4 h of infection with Leishmania amazonensis wild-type (La-WT) or L. amazonensis arginase knockout (La-arg-). We identified 12 641 host transcripts and 8282 parasite transcripts by alignment analysis with the respective Mus musculus and L. mexicana genomes. The comparison of BALB/c_La-arg-versus BALB/c_La-WT revealed 233 modulated transcripts, with most related to the immune response and some related to the amino acid transporters and l-arginine metabolism. In contrast, the comparison of C57BL/6_La-arg-vs. C57BL/6_La-WT revealed only 30 modulated transcripts, including some related to the immune response but none related to amino acid transport or l-arginine metabolism. The transcriptome profiles of the intracellular amastigote revealed 94 modulated transcripts in the comparison of La-arg-_BALB/c vs. La-WT_BALB/c and 45 modulated transcripts in the comparison of La-arg-_C57BL/6 vs. La-WT_C57BL/6. Taken together, our data present new insights into the impact of parasite arginase activity on the orchestration of the host gene expression modulation, including in the immune response and amino acid transport and metabolism, mainly in susceptible BALB/c-infected macrophages. Moreover, we show how parasite arginase activity affects parasite gene expression modulation, including amino acid uptake and amastin expression.

Frequent Recombination Events in Leishmania donovani: Mining Population Data

The Leishmania donovani species complex consists of all L. donovani and L. infantum strains mainly responsible for visceral leishmaniasis (VL). It was suggested that genome rearrangements in Leishmania spp. occur very often, thus enabling parasites to adapt to the different environmental conditions. Some of these rearrangements may be directly linked to the virulence or explain the reduced efficacy of antimonial drugs in some isolates. In the current study, we focused on a large-scale analysis of putative gene conversion events using publicly available datasets. Previous population study of L. donovani suggested that population variability of L. donovani is relatively low, however the authors used masking procedures and strict read selection criteria. We decided to re-analyze DNA-seq data without masking sequences, because we were interested in the most dynamic fraction of the genome. The majority of samples have an excess of putative gene conversion/recombination events in the noncoding regions, however we found an overall excess of putative intrachromosomal gene conversion/recombination in the protein coding genes, compared to putative interchromosomal gene conversion/recombination events.

Trypanosomatid parasite dynamically changes the transcriptome during infection and modifies honey bee physiology

It is still not understood how honey bee parasite changes the gene expression to adapt to the host environment and how the host simultaneously responds to the parasite infection by modifying its own gene expression. To address this question, we studied a trypanosomatid, Lotmaria passim, which can be cultured in medium and inhabit the honey bee hindgut. We found that L. passim decreases mRNAs associated with protein translation, glycolysis, detoxification of radical oxygen species, and kinetoplast respiratory chain to adapt to the anaerobic and nutritionally poor honey bee hindgut during the infection. After the long term infection, the host appears to be in poor nutritional status, indicated by the increase and decrease of take-out and vitellogenin mRNAs, respectively. Simultaneous gene expression profiling of L. passim and honey bee during infection by dual RNA-seq provided insight into how both parasite and host modify their gene expressions.

Comparative genomics of Leishmania (Mundinia)

Background

Trypanosomatids of the genus Leishmania are parasites of mammals or reptiles transmitted by bloodsucking dipterans. Many species of these flagellates cause important human diseases with clinical symptoms ranging from skin sores to life-threatening damage of visceral organs. The genus Leishmania contains four subgenera: Leishmania, Sauroleishmania, Viannia, and Mundinia. The last subgenus has been established recently and remains understudied, although Mundinia contains human-infecting species. In addition, it is interesting from the evolutionary viewpoint, representing the earliest branch within the genus and possibly with a different type of vector. Here we analyzed the genomes of L. (M.) martiniquensis, L. (M.) enriettii and L. (M.) macropodum to better understand the biology and evolution of these parasites.

Results

All three genomes analyzed were approximately of the same size (~ 30 Mb) and similar to that of L. (Sauroleishmania) tarentolae, but smaller than those of the members of subgenera Leishmania and Viannia, or the genus Endotrypanum (~ 32 Mb). This difference was explained by domination of gene losses over gains and contractions over expansions at the Mundinia node, although only a few of these genes could be identified. The analysis predicts significant changes in the Mundinia cell surface architecture, with the most important ones relating to losses of LPG-modifying side chain galactosyltransferases and arabinosyltransferases, as well as β-amastins. Among other important changes were gene family contractions for the oxygen-sensing adenylate cyclases and FYVE zinc finger-containing proteins.

Conclusions

We suggest that adaptation of Mundinia to different vectors and hosts has led to alternative host-parasite relationships and, thereby, made some proteins redundant. Thus, the evolution of genomes in the genus Leishmania and, in particular, in the subgenus Mundinia was mainly shaped by host (or vector) switches.

Transcriptional responses of Leishmania (Leishmania) amazonensis in the presence of trivalent sodium stibogluconate

Background

In the last decade, resistance to antimonials has become a serious problem due to the emergence of drug-resistant strains. Therefore, understanding the mechanisms used by Leishmania parasites to survive under drug pressure is essential, particularly for species of medical-veterinary importance such as L. amazonensis.

Methods

Here, we used RNA-seq technology to analyse transcriptome profiles and identify global changes in gene expression between antimony-resistant and -sensitive L. amazonensis promastigotes.

Results

A total of 723 differentially expressed genes were identified between resistant and sensitive lines. Comparative transcriptomic analysis revealed that genes encoding proteins involved in metabolism (fatty acids) and stress response, as well as those associated with antimony resistance in other Leishmania species, were upregulated in the antimony-resistant line. Most importantly, we observed upregulation of genes encoding autophagy proteins, suggesting that in the presence of trivalent stibogluconate (Sb^III) L. amazonensis can activate these genes either as a survival strategy or to induce cell death, as has been observed in other parasites.

Conclusions

This work identified global transcriptomic changes in an in vitro-adapted strain in response to Sb^III. Our results provide relevant information to continue understanding the mechanism used by parasites of the subgenus Leishmania (L. amazonensis) to generate an antimony-resistant phenotype.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3603-8) contains supplementary material, which is available to authorized users.

Analysis by RNA-seq of transcriptomic changes elicited by heat shock in Leishmania major

Besides their medical relevance, Leishmania is an adequate model for studying post-transcriptional mechanisms of gene expression. In this microorganism, mRNA degradation/stabilization mechanisms together with translational control and post-translational modifications of proteins are the major drivers of gene expression. Leishmania parasites develop as promastigotes in sandflies and as amastigotes in mammalians, and during host transmission, the parasite experiences a sudden temperature increase. Here, changes in the transcriptome of Leishmania major promastigotes after a moderate heat shock were analysed by RNA-seq. Several of the up-regulated transcripts code for heat shock proteins, other for proteins previously reported to be amastigote-specific and many for hypothetical proteins. Many of the transcripts experiencing a decrease in their steady-state levels code for transporters, proteins involved in RNA metabolism or translational factors. In addition, putative long noncoding RNAs were identified among the differentially expressed transcripts. Finally, temperature-dependent changes in the selection of the spliced leader addition sites were inferred from the RNA-seq data, and particular cases were further validated by RT-PCR and Northern blotting. This study provides new insights into the post-transcriptional mechanisms by which Leishmania modulate gene expression.

Diagnostic evaluation of the amastin protein from Leishmania infantum in canine and human visceral leishmaniasis and immunogenicity in human cells derived from patients and healthy controls

The diagnosis of visceral leishmaniasis (VL) presents problems due to the toxicity and/or high cost of drugs. In addition, no vaccine exists to protect against human disease. In this study, the antigenicity and immunogenicity of amastin protein were evaluated in L. infantum-infected dogs and humans. For the diagnosis, besides the recombinant protein, 1 linear B-cell epitope was synthetized and evaluated in serological assays. Results showed high sensitivity and specificity values to detect the disease when both antigens were employed against a canine and human serological panel. By contrast, when using rA2 and a soluble Leishmania antigenic preparation, sensitivity and specificity values proved to be lower. A preliminary immunogenicity study showed that the amastin protein induced high IFN-γ and low IL-10 production in stimulated PBMC derived from treated VL patients and healthy subjects, thus suggesting a potential use of this protein as an immunogen to protect against human disease.

Bioinformatics analysis of four proteins of Leishmania donovani to guide epitopes vaccine design and drug targets selection

Visceral leishmaniasis (VL) is a serious and widespread parasitic disease caused by Leishmania donovani complex. The threat of this fatal disease continues due to the lack of ideal drugs or vaccines. In this study, we selected Amastin, CaNA2, Kmp-11 and PDI proteins of Leishmania donovani for study, which are VL vaccine candidates or possible drug targets. Eleven bioinformatics tools were used to analyze different aspects of these proteins, including amino acid composition, topology, signal peptide, secondary structure, surface properties, phosphorylation sites and kinases, protein binding sites, 3D homology modeling, B cell epitopes, MHC class Ⅰ and Ⅱ epitopes and protein-protein interactions. Finally, the functionally related amino acid sites and dominant epitopes of these proteins were founded. Some possible relationships between protein structure, phosphorylation sites, protein binding sites and epitopes were also discovered. High flexibility and random coils regions of protein have a tendency to be phosphorylated, bind proteins and present epitopes. Since some phosphorylation sites and their kinases are involved in Leishmania invasion and survival in host cells, they may be potential drug targets. Bioinformatics analysis helps us better understand protein function and find dominant epitopes to guide drug design and vaccine development.

Genomic insights into virulence mechanisms of Leishmania donovani: evidence from an atypical strain

Background

Leishmaniasis is a neglected tropical disease with diverse clinical phenotypes, determined by parasite, host and vector interactions. Despite the advances in molecular biology and the availability of more Leishmania genome references in recent years, the association between parasite species and distinct clinical phenotypes remains poorly understood. We present a genomic comparison of an atypical variant of Leishmania donovani from a South Asian focus, where it mostly causes cutaneous form of leishmaniasis.

Results

Clinical isolates from six cutaneous leishmaniasis patients (CL-SL); 2 of whom were poor responders to antimony (CL-PR), and two visceral leishmaniasis patients (VL-SL) were sequenced on an Illumina MiSeq platform. Chromosome aneuploidy was observed in both groups but was more frequent in CL-SL. 248 genes differed by 2 fold or more in copy number among the two groups. Genes involved in amino acid use (LdBPK_271940) and energy metabolism (LdBPK_271950), predominated the VL-SL group with the same distribution pattern reflected in gene tandem arrays. Genes encoding amastins were present in higher copy numbers in VL-SL and CL-PR as well as being among predicted pseudogenes in CL-SL. Both chromosome and SNP profiles showed CL-SL and VL-SL to form two distinct groups. While expected heterozygosity was much higher in VL-SL, SNP allele frequency patterns did not suggest potential recent recombination breakpoints. The SNP/indel profile obtained using the more recently generated PacBio sequence did not vary markedly from that based on the standard LdBPK282A1 reference. Several genes previously associated with resistance to antimonials were observed in higher copy numbers in the analysis of CL-PR. H-locus amplification was seen in one cutaneous isolate which however did not belong to the CL-PR group.

Conclusions

The data presented suggests that intra species variations at chromosome and gene level are more likely to influence differences in tropism as well as response to treatment, and contributes to greater understanding of parasite molecular mechanisms underpinning these differences. These findings should be substantiated with a larger sample number and expression/functional studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5271-z) contains supplementary material, which is available to authorized users.

Genomic comparison of Trypanosoma conorhini and Trypanosoma rangeli to Trypanosoma cruzi strains of high and low virulence

Background

Trypanosoma conorhini and Trypanosoma rangeli, like Trypanosoma cruzi, are kinetoplastid protist parasites of mammals displaying divergent hosts, geographic ranges and lifestyles. Largely nonpathogenic T. rangeli and T. conorhini represent clades that are phylogenetically closely related to the T. cruzi and T. cruzi-like taxa and provide insights into the evolution of pathogenicity in those parasites. T. rangeli, like T. cruzi is endemic in many Latin American countries, whereas T. conorhini is tropicopolitan. T. rangeli and T. conorhini are exclusively extracellular, while T. cruzi has an intracellular stage in the mammalian host.

Results

Here we provide the first comprehensive sequence analysis of T. rangeli AM80 and T. conorhini 025E, and provide a comparison of their genomes to those of T. cruzi G and T. cruzi CL, respectively members of T. cruzi lineages TcI and TcVI. We report de novo assembled genome sequences of the low-virulent T. cruzi G, T. rangeli AM80, and T. conorhini 025E ranging from ~ 21–25 Mbp, with ~ 10,000 to 13,000 genes, and for the highly virulent and hybrid T. cruzi CL we present a ~ 65 Mbp in-house assembled haplotyped genome with ~ 12,500 genes per haplotype. Single copy orthologs of the two T. cruzi strains exhibited ~ 97% amino acid identity, and ~ 78% identity to proteins of T. rangeli or T. conorhini. Proteins of the latter two organisms exhibited ~ 84% identity. T. cruzi CL exhibited the highest heterozygosity. T. rangeli and T. conorhini displayed greater metabolic capabilities for utilization of complex carbohydrates, and contained fewer retrotransposons and multigene family copies, i.e. trans-sialidases, mucins, DGF-1, and MASP, compared to T. cruzi.

Conclusions

Our analyses of the T. rangeli and T. conorhini genomes closely reflected their phylogenetic proximity to the T. cruzi clade, and were largely consistent with their divergent life cycles. Our results provide a greater context for understanding the life cycles, host range expansion, immunity evasion, and pathogenesis of these trypanosomatids.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5112-0) contains supplementary material, which is available to authorized users.

Conserved epitopes in variants of amastin protein of Trypanosoma cruzi for vaccine design: A bioinformatics approach

Chagas disease caused by protozoan parasite Trypanosoma cruzi is endemic disease in South and Central American countries but due to migrating human populations it has shown emergence in Europe, North America and Australia. With only two drugs, benznidazole and nifurtimox for its treatment there is need for newer therapies. In the current study, we have tried to analyse the potential of amastin, a major surface protein as a vaccine target using bioinformatics tools. Using 282 variants of this protein available in NCBI protein database we have found out five conserved potential T_c cell and two T_H cell epitopes. These epitopes are conserved in more than 90% of the cohort of sequences used in the study. The epitopes showed binding to the peptide binding cleft of HLA-A02 and HLA-DR molecules. With coverage of pan world populations and being non-toxic and non-allergic these epitopes could be used for future vaccine applications.

A systematic review of the Trypanosoma cruzi genetic heterogeneity, host immune response and genetic factors as plausible drivers of chronic chagasic cardiomyopathy

Chagas disease is a complex tropical pathology caused by the kinetoplastid Trypanosoma cruzi. This parasite displays massive genetic diversity and has been classified by international consensus in at least six Discrete Typing Units (DTUs) that are broadly distributed in the American continent. The main clinical manifestation of the disease is the chronic chagasic cardiomyopathy (CCC) that is lethal in the infected individuals. However, one intriguing feature is that only 30-40% of the infected individuals will develop CCC. Some authors have suggested that the immune response, host genetic factors, virulence factors and even the massive genetic heterogeneity of T. cruzi are responsible of this clinical pattern. To date, no conclusive data support the reason why a few percentages of the infected individuals will develop CCC. Therefore, we decided to conduct a systematic review analysing the host genetic factors, immune response, cytokine production, virulence factors and the plausible association of the parasite DTUs and CCC. The epidemiological and clinical implications are herein discussed.

The cotton stainer's gut microbiota suppresses infection of a cotransmitted trypanosomatid parasite

The evolutionary and ecological success of many insects is attributed to mutualistic partnerships with bacteria that confer hosts with novel traits including food digestion, nutrient supplementation, detoxification of harmful compounds and defence against natural enemies. Dysdercus fasciatus firebugs (Hemiptera: Pyrrhocoridae), commonly known as cotton stainers, possess a simple but distinctive gut bacterial community including B vitamin-supplementing Coriobacteriaceae symbionts. In addition, their guts are often infested with the intestinal trypanosomatid parasite Leptomonas pyrrhocoris (Kinetoplastida: Trypanosomatidae). In this study, using experimental bioassays and fluorescence in situ hybridization (FISH), we report on the protective role of the D. fasciatus gut bacteria against L. pyrrhocoris. We artificially infected 2nd instars of dysbiotic and symbiotic insects with a parasite culture and measured parasite titres, developmental time and survival rates. Our results show that L. pyrrhocoris infection increases developmental time and slightly modifies the quantitative composition of the gut microbiota. More importantly, we found significantly higher parasite titres and a tendency towards lower survival rates in parasite-infected dysbiotic insects compared to symbiotic controls, indicating that the gut bacteria successfully interfere with the establishment or proliferation of L. pyrrhocoris. The colonization of symbiotic bacteria on the peritrophic matrix along the gut wall, as revealed by FISH, likely acts as a barrier blocking parasite attachment or entry into the hemolymph. Our findings show that in addition to being nutritionally important, D. fasciatus' gut bacteria complement the host's immune system in preventing parasite invasions and that a stable gut microbial community is integral for the host's health.

An Alternative Strategy for Trypanosome Survival in the Mammalian Bloodstream Revealed through Genome and Transcriptome Analysis of the Ubiquitous Bovine Parasite Trypanosoma (Megatrypanum) theileri

There are hundreds of Trypanosoma species that live in the blood and tissue spaces of their vertebrate hosts. The vast majority of these do not have the ornate system of antigenic variation that has evolved in the small number of African trypanosome species, but can still maintain long-term infections in the face of the vertebrate adaptive immune system. Trypanosoma theileri is a typical example, has a restricted host range of cattle and other Bovinae, and is only occasionally reported to cause patent disease although no systematic survey of the effect of infection on agricultural productivity has been performed. Here, a detailed genome sequence and a transcriptome analysis of gene expression in bloodstream form T. theileri have been performed. Analysis of the genome sequence and expression showed that T. theileri has a typical kinetoplastid genome structure and allowed a prediction that it is capable of meiotic exchange, gene silencing via RNA interference and, potentially, density-dependent growth control. In particular, the transcriptome analysis has allowed a comparison of two distinct trypanosome cell surfaces, T. brucei and T. theileri, that have each evolved to enable the maintenance of a long-term extracellular infection in cattle. The T. theileri cell surface can be modeled to contain a mixture of proteins encoded by four novel large and divergent gene families and by members of a major surface protease gene family. This surface composition is distinct from the uniform variant surface glycoprotein coat on African trypanosomes providing an insight into a second mechanism used by trypanosome species that proliferate in an extracellular milieu in vertebrate hosts to avoid the adaptive immune response.

Genomic and transcriptomic alterations in Leishmania donovani lines experimentally resistant to antileishmanial drugs

Leishmaniasis is a serious medical issue in many countries around the World, but it remains largely neglected in terms of research investment for developing new control and treatment measures. No vaccines exist for human use, and the chemotherapeutic agents currently used are scanty. Furthermore, for some drugs, resistance and treatment failure are increasing to alarming levels. The aim of this work was to identify genomic and trancriptomic alterations associated with experimental resistance against the common drugs used against VL: trivalent antimony (Sb^III, S line), amphotericin B (AmB, A line), miltefosine (MIL, M line) and paromomycin (PMM, P line). A total of 1006 differentially expressed transcripts were identified in the S line, 379 in the A line, 146 in the M line, and 129 in the P line. Also, changes in ploidy of chromosomes and amplification/deletion of particular regions were observed in the resistant lines regarding the parental one. A series of genes were identified as possible drivers of the resistance phenotype and were validated in both promastigotes and amastigotes from Leishmania donovani, Leishmania infantum and Leishmania major species. Remarkably, a deletion of the gene LinJ.36.2510 (coding for 24-sterol methyltransferase, SMT) was found to be associated with AmB-resistance in the A line. In the P line, a dramatic overexpression of the transcripts LinJ.27.T1940 and LinJ.27.T1950 that results from a massive amplification of the collinear genes was suggested as one of the mechanisms of PMM resistance. This conclusion was reinforced after transfection experiments in which significant PMM-resistance was generated in WT parasites over-expressing either gene LinJ.27.1940 (coding for a D-lactate dehydrogenase-like protein, D-LDH) or gene LinJ.27.1950 (coding for an aminotransferase of branched-chain amino acids, BCAT). This work allowed to identify new drivers, like SMT, the deletion of which being associated with resistance to AmB, and the tandem D-LDH-BCAT, the amplification of which being related to PMM resistance.

Leishmania naiffi and Leishmania guyanensis reference genomes highlight genome structure and gene evolution in the Viannia subgenus

The unicellular protozoan parasite Leishmania causes the neglected tropical disease leishmaniasis, affecting 12 million people in 98 countries. In South America, where the Viannia subgenus predominates, so far only L. (Viannia) braziliensis and L. (V.) panamensis have been sequenced, assembled and annotated as reference genomes. Addressing this deficit in molecular information can inform species typing, epidemiological monitoring and clinical treatment. Here, L. (V.) naiffi and L. (V.) guyanensis genomic DNA was sequenced to assemble these two genomes as draft references from short sequence reads. The methods used were tested using short sequence reads for L. braziliensis M2904 against its published reference as a comparison. This assembly and annotation pipeline identified 70 additional genes not annotated on the original M2904 reference. Phylogenetic and evolutionary comparisons of L. guyanensis and L. naiffi with 10 other Viannia genomes revealed four traits common to all Viannia: aneuploidy, 22 orthologous groups of genes absent in other Leishmania subgenera, elevated TATE transposon copies and a high NADH-dependent fumarate reductase gene copy number. Within the Viannia, there were limited structural changes in genome architecture specific to individual species: a 45 Kb amplification on chromosome 34 was present in all bar L. lainsoni, L. naiffi had a higher copy number of the virulence factor leishmanolysin, and laboratory isolate L. shawi M8408 had a possible minichromosome derived from the 3' end of chromosome 34. This combination of genome assembly, phylogenetics and comparative analysis across an extended panel of diverse Viannia has uncovered new insights into the origin and evolution of this subgenus and can help improve diagnostics for leishmaniasis surveillance.

RNA-seq analysis reveals differences in transcript abundance between cultured and sand fly-derived Leishmania infantum promastigotes

Leishmania infantum is responsible for human and canine leishmaniasis in the Mediterranean basin, where the major vector is Phlebotomus perniciosus. Because isolation of sufficient parasites from the sand fly gut is technically challenging, axenic cultivation of promastigotes is routinely used to obtain material for biochemical and genetic analyses. Here, we report the use of Spliced Leader RNA-seq (SL-seq) to compare transcript abundance in cultured promastigotes and those obtained from the whole midgut of the sand fly 5 days after infection. SL-seq allows for amplification of RNA from the parasite avoiding contamination with RNA from the gut of the insect. The study has been performed by means of a single technical replicate comparing pools of samples obtained from sand fly-derived (sfPro) and axenic culture promastigotes (acPro). Although there was a moderate correlation (R² = 0.83) in gene expression, 793 genes showed significantly different (≥2-fold, p <0.05) mRNA levels in sand fly-derived promastigotes and in culture, of which 31 were up-regulated ≥8-fold (p < 10^-8 in most cases). These included several genes that are typically up-regulated during metacyclogenesis, suggesting that sand fly-derived promastigotes contain a substantial number of metacyclics, and/or that their differentiation status as metacyclics is more advanced in these populations. Infection experiments and studies evaluating the proportion of metacyclic promastigotes in culture and within the sand fly gut, previously reported by us, support the last hypothesis.

Trypanosomatids Are Much More than Just Trypanosomes: Clues from the Expanded Family Tree

Trypanosomes and leishmanias are widely known parasites of humans. However, they are just two out of several phylogenetic lineages that constitute the family Trypanosomatidae. Although dixeny - the ability to infect two hosts - is a derived trait of vertebrate-infecting parasites, the majority of trypanosomatids are monoxenous. Like their common ancestor, the monoxenous Trypanosomatidae are mostly parasites or commensals of insects. This review covers recent advances in the study of insect trypanosomatids, highlighting their diversity as well as genetic, morphological and biochemical complexity, which, until recently, was underappreciated. The investigation of insect trypanosomatids is providing an important foundation for understanding the origin and evolution of parasitism, including colonization of vertebrates and the appearance of human pathogens.

Proteomic Profiling of Leishmania donovani Promastigote Subcellular Organelles

To facilitate a greater understanding of the biological processes in the medically important Leishmania donovani parasite, a combination of differential and density-gradient ultracentrifugation techniques were used to achieve a comprehensive subcellular fractionation of the promastigote stage. An in-depth label-free proteomic LC-MS/MS analysis of the density gradients resulted in the identification of ∼50% of the Leishmania proteome (3883 proteins detected), which included ∼645 integral membrane proteins and 1737 uncharacterized proteins. Clustering and subcellular localization of proteins was based on a subset of training Leishmania proteins with known subcellular localizations that had been determined using biochemical, confocal microscopy, or immunoelectron microscopy approaches. This subcellular map will be a valuable resource that will help dissect the cell biology and metabolic processes associated with specific organelles of Leishmania and related kinetoplastids.

Gene and Chromosomal Copy Number Variations as an Adaptive Mechanism Towards a Parasitic Lifestyle in Trypanosomatids

Trypanosomatids are a group of kinetoplastid parasites including some of great public health importance, causing debilitating and life-long lasting diseases that affect more than 24 million people worldwide. Among the trypanosomatids, Trypanosoma cruzi, Trypanosoma brucei and species from the Leishmania genus are the most well studied parasites, due to their high prevalence in human infections. These parasites have an extreme genomic and phenotypic variability, with a massive expansion in the copy number of species-specific multigene families enrolled in host-parasite interactions that mediate cellular invasion and immune evasion processes. As most trypanosomatids are heteroxenous, and therefore their lifecycles involve the transition between different hosts, these parasites have developed several strategies to ensure a rapid adaptation to changing environments. Among these strategies, a rapid shift in the repertoire of expressed genes, genetic variability and genome plasticity are key mechanisms. Trypanosomatid genomes are organized into large directional gene clusters that are transcribed polycistronically, where genes derived from the same polycistron may have very distinct mRNA levels. This particular mode of transcription implies that the control of gene expression operates mainly at post-transcriptional level. In this sense, gene duplications/losses were already associated with changes in mRNA levels in these parasites. Gene duplications also allow the generation of sequence variability, as the newly formed copy can diverge without loss of function of the original copy. Recently, aneuploidies have been shown to occur in several Leishmania species and T. cruzi strains. Although aneuploidies are usually associated with debilitating phenotypes in superior eukaryotes, recent data shows that it could also provide increased fitness in stress conditions and generate drug resistance in unicellular eukaryotes. In this review, we will focus on gene and chromosomal copy number variations and their relevance to the evolution of trypanosomatid parasites.

The genomes of Crithidia bombi and C. expoeki, common parasites of bumblebees

Trypanosomatids (Trypanosomatidae, Kinetoplastida) are flagellated protozoa containing many parasites of medical or agricultural importance. Among those, Crithidia bombi and C. expoeki, are common parasites in bumble bees around the world, and phylogenetically close to Leishmania and Leptomonas. They have a simple and direct life cycle with one host, and partially castrate the founding queens greatly reducing their fitness. Here, we report the nuclear genome sequences of one clone of each species, extracted from a field-collected infection. Using a combination of Roche 454 FLX Titanium, Pacific Biosciences PacBio RS, and Illumina GA2 instruments for C. bombi, and PacBio for C. expoeki, we could produce high-quality and well resolved sequences. We find that these genomes are around 32 and 34 MB, with 7,808 and 7,851 annotated genes for C. bombi and C. expoeki, respectively-which is somewhat less than reported from other trypanosomatids, with few introns, and organized in polycistronic units. A large fraction of genes received plausible functional support in comparison primarily with Leishmania and Trypanosoma. Comparing the annotated genes of the two species with those of six other trypanosomatids (C. fasciculata, L. pyrrhocoris, L. seymouri, B. ayalai, L. major, and T. brucei) shows similar gene repertoires and many orthologs. Similar to other trypanosomatids, we also find signs of concerted evolution in genes putatively involved in the interaction with the host, a high degree of synteny between C. bombi and C. expoeki, and considerable overlap with several other species in the set. A total of 86 orthologous gene groups show signatures of positive selection in the branch leading to the two Crithidia under study, mostly of unknown function. As an example, we examined the initiating glycosylation pathway of surface components in C. bombi, finding it deviates from most other eukaryotes and also from other kinetoplastids, which may indicate rapid evolution in the extracellular matrix that is involved in interactions with the host. Bumble bees are important pollinators and Crithidia-infections are suspected to cause substantial selection pressure on their host populations. These newly sequenced genomes provide tools that should help better understand host-parasite interactions in these pollinator pathogens.

Structural and evolutionary analysis of Leishmania Alba proteins

The Alba superfamily proteins share a common RNA-binding domain. These proteins participate in a variety of regulatory pathways by controlling developmental gene expression. They also interact with ribosomal subunits, translation factors, and other RNA-binding proteins. The Leishmania infantum genome encodes two Alba-domain proteins, LiAlba1 and LiAlba3. In this work, we used homology modeling, protein-protein docking, and molecular dynamics (MD) simulations to explore the details of the Alba1-Alba3-RNA complex from Leishmania infantum at the molecular level. In addition, we compared the structure of LiAlba3 with the human ribonuclease P component, Rpp20. We also mapped the ligand-binding residues on the Alba3 surface to analyze its druggability and performed mutational analyses in Alba3 using alanine scanning to identify residues involved in its function and structural stability. These results suggest that the RGG-box motif of LiAlba1 is important for protein function and stability. Finally, we discuss the function of Alba proteins in the context of pathogen adaptation to host cells. The data provided herein will facilitate further translational research regarding Alba structure and function.

A Brief View of the Surface Membrane Proteins from Trypanosoma cruzi

Trypanosoma cruzi is the causal agent of Chagas' disease which affects millions of people around the world mostly in Central and South America. T. cruzi expresses a wide variety of proteins on its surface membrane which has an important role in the biology of these parasites. Surface molecules of the parasites are the result of the environment to which the parasites are exposed during their life cycle. Hence, T. cruzi displays several modifications when they move from one host to another. Due to the complexity of this parasite's cell surface, this review presents some membrane proteins organized as large families, as they are the most abundant and/or relevant throughout the T. cruzi membrane.

Modulation of Aneuploidy in Leishmania donovani during Adaptation to Different In Vitro and In Vivo Environments and Its Impact on Gene Expression

Aneuploidy is usually deleterious in multicellular organisms but appears to be tolerated and potentially beneficial in unicellular organisms, including pathogens. Leishmania, a major protozoan parasite, is emerging as a new model for aneuploidy, since in vitro-cultivated strains are highly aneuploid, with interstrain diversity and intrastrain mosaicism. The alternation of two life stages in different environments (extracellular promastigotes and intracellular amastigotes) offers a unique opportunity to study the impact of environment on aneuploidy and gene expression. We sequenced the whole genomes and transcriptomes of Leishmania donovani strains throughout their adaptation to in vivo conditions mimicking natural vertebrate and invertebrate host environments. The nucleotide sequences were almost unchanged within a strain, in contrast to highly variable aneuploidy. Although high in promastigotes in vitro, aneuploidy dropped significantly in hamster amastigotes, in a progressive and strain-specific manner, accompanied by the emergence of new polysomies. After a passage through a sand fly, smaller yet consistent karyotype changes were detected. Changes in chromosome copy numbers were correlated with the corresponding transcript levels, but additional aneuploidy-independent regulation of gene expression was observed. This affected stage-specific gene expression, downregulation of the entire chromosome 31, and upregulation of gene arrays on chromosomes 5 and 8. Aneuploidy changes in Leishmania are probably adaptive and exploited to modulate the dosage and expression of specific genes; they are well tolerated, but additional mechanisms may exist to regulate the transcript levels of other genes located on aneuploid chromosomes. Our model should allow studies of the impact of aneuploidy on molecular adaptations and cellular fitness.IMPORTANCE Aneuploidy is usually detrimental in multicellular organisms, but in several microorganisms, it can be tolerated and even beneficial. Leishmania-a protozoan parasite that kills more than 30,000 people each year-is emerging as a new model for aneuploidy studies, as unexpectedly high levels of aneuploidy are found in clinical isolates. Leishmania lacks classical regulation of transcription at initiation through promoters, so aneuploidy could represent a major adaptive strategy of this parasite to modulate gene dosage in response to stressful environments. For the first time, we document the dynamics of aneuploidy throughout the life cycle of the parasite, in vitro and in vivo We show its adaptive impact on transcription and its interaction with regulation. Besides offering a new model for aneuploidy studies, we show that further genomic studies should be done directly in clinical samples without parasite isolation and that adequate methods should be developed for this.

Involvement of an RNA binding protein containing Alba domain in the stage-specific regulation of beta-amastin expression in Trypanosoma cruzi

Amastins are surface glycoproteins, first identified in amastigotes of T. cruzi but later found to be expressed in several Leishmania species, as well as in T. cruzi epimastigotes. Amastins are encoded by a diverse gene family that can be grouped into four subfamilies named α, β, γ, and δ amastins. Differential expression of amastin genes results from regulatory mechanisms involving changes in mRNA stability and/or translational control. Although distinct regulatory elements were identified in the 3' UTR of T. cruzi and Leishmania amastin mRNAs, RNA binding proteins involved with amastin gene regulation have only being characterized in L. infantum where an Alba-domain protein (LiAlba20) able to bind to the 3' UTR of a δ-amastin mRNA was identified. Here we investigated the role of TcAlba30, the LiAlba20 homologue in T. cruzi, in the post transcriptional regulation of amastin genes. TcAlba30 transcripts are present in all stages of the T. cruzi life cycle. RNA immunoprecipitation assays using a transfected cell line expressing a cMyc tagged TcAlba30 revealed that TcAlba30 can interact with β-amastin mRNA. In addition, over-expression of TcAlba30 in epimastigotes resulted in 50% decreased levels of β-amastin mRNAs compared to wild type parasites. Since luciferase assays indicated the presence of regulatory elements in the 3' UTR of β-amastin mRNA and reduced levels of luciferase mRNA were found in parasites over expressing TcAlba30, we conclude that TcAlba30 acts as a T. cruzi RNA binding protein involved in the negative control of β-amastin expression through interactions with its 3'UTR.

Developmental differentiation in Leishmania lifecycle progression: post-transcriptional control conducts the orchestra

The successful progression of Leishmania spp. through their lifecycle entails a series of differentiation processes; the proliferative procyclic promastigote forms become quiescent, human-infective metacyclic promastigotes during metacyclogenesis in the sandfly vector, which then differentiate into amastigotes during amastigogenesis in the mammalian host. The progression to these infective forms requires two components: environmental cues and a coordinated cellular response. Recent studies have shown that the Leishmania cellular transformation into mammalian-infective stages is triggered by broad changes in the absolute and relative RNA and protein levels. In this review, we will discuss the implications of Leishmania transcriptomic and proteomic fluctuations, which adapt the parasitic cell for survival.

Development of a toolbox to dissect host-endosymbiont interactions and protein trafficking in the trypanosomatid Angomonas deanei

Background

Bacterial endosymbionts are found across the eukaryotic kingdom and profoundly impacted eukaryote evolution. In many endosymbiotic associations with vertically inherited symbionts, highly complementary metabolic functions encoded by host and endosymbiont genomes indicate integration of metabolic processes between the partner organisms. While endosymbionts were initially expected to exchange only metabolites with their hosts, recent evidence has demonstrated that also host-encoded proteins can be targeted to the bacterial symbionts in various endosymbiotic systems. These proteins seem to participate in regulating symbiont growth and physiology. However, mechanisms required for protein targeting and the specific endosymbiont targets of these trafficked proteins are currently unexplored owing to a lack of molecular tools that enable functional studies of endosymbiotic systems.

Results

Here we show that the trypanosomatid Angomonas deanei, which harbors a β-proteobacterial endosymbiont, is readily amenable to genetic manipulation. Its rapid growth, availability of full genome and transcriptome sequences, ease of transfection, and high frequency of homologous recombination have allowed us to stably integrate transgenes into the A. deanei nuclear genome, efficiently generate null mutants, and elucidate protein localization by heterologous expression of a fluorescent protein fused to various putative targeting signals. Combining these novel tools with proteomic analysis was key for demonstrating the routing of a host-encoded protein to the endosymbiont, suggesting the existence of a specific endosymbiont-sorting machinery in A. deanei.

Conclusions

After previous reports from plants, insects, and a cercozoan amoeba we found here that also in A. deanei, i.e. a member of a fourth eukaryotic supergroup, host-encoded proteins can be routed to the bacterial endosymbiont. This finding adds further evidence to our view that the targeting of host proteins is a general strategy of eukaryotes to gain control over and interact with a bacterial endosymbiont. The molecular resources reported here establish A. deanei as a time and cost efficient reference system that allows for a rigorous dissection of host-symbiont interactions that have been, and are still being shaped over evolutionary time. We expect this system to greatly enhance our understanding of the biology of endosymbiosis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0820-z) contains supplementary material, which is available to authorized users.

Symbiont modulates expression of specific gene categories in Angomonas deanei

Trypanosomatids are parasites that cause disease in humans, animals, and plants. Most are non-pathogenic and some harbor a symbiotic bacterium. Endosymbiosis is part of the evolutionary process of vital cell functions such as respiration and photosynthesis. Angomonas deanei is an example of a symbiont-containing trypanosomatid. In this paper, we sought to investigate how symbionts influence host cells by characterising and comparing the transcriptomes of the symbiont-containing A. deanei (wild type) and the symbiont-free aposymbiotic strains. The comparison revealed that the presence of the symbiont modulates several differentially expressed genes. Empirical analysis of differential gene expression showed that 216 of the 7625 modulated genes were significantly changed. Finally, gene set enrichment analysis revealed that the largest categories of genes that downregulated in the absence of the symbiont were those involved in oxidation-reduction process, ATP hydrolysis coupled proton transport and glycolysis. In contrast, among the upregulated gene categories were those involved in proteolysis, microtubule-based movement, and cellular metabolic process. Our results provide valuable information for dissecting the mechanism of endosymbiosis in A. deanei.

Autonomy and integration in complex parasite life cycles

Complex life cycles are common in free-living and parasitic organisms alike. The adaptive decoupling hypothesis postulates that separate life cycle stages have a degree of developmental and genetic autonomy, allowing them to be independently optimized for dissimilar, competing tasks. That is, complex life cycles evolved to facilitate functional specialization. Here, I review the connections between the different stages in parasite life cycles. I first examine evolutionary connections between life stages, such as the genetic coupling of parasite performance in consecutive hosts, the interspecific correlations between traits expressed in different hosts, and the developmental and functional obstacles to stage loss. Then, I evaluate how environmental factors link life stages through carryover effects, where stressful larval conditions impact parasites even after transmission to a new host. There is evidence for both autonomy and integration across stages, so the relevant question becomes how integrated are parasite life cycles and through what mechanisms? By highlighting how genetics, development, selection and the environment can lead to interdependencies among successive life stages, I wish to promote a holistic approach to studying complex life cycle parasites and emphasize that what happens in one stage is potentially highly relevant for later stages.

Genomic Appraisal of the Multifactorial Basis for In Vitro Acquisition of Miltefosine Resistance in Leishmania donovani

Protozoan parasites of the Leishmania donovani complex are the causative agents of visceral leishmaniasis (VL), the most severe form of leishmaniasis, with high rates of mortality if left untreated. Leishmania parasites are transmitted to humans through the bite of infected female sandflies (Diptera: Phlebotominae), and approximately 500,000 new cases of VL are reported each year. In the absence of a safe human vaccine, chemotherapy, along with vector control, is the sole tool with which to fight the disease. Miltefosine (hexadecylphosphatidylcholine [HePC]), an antitumoral drug, is the only successful oral treatment for VL. In the current study, we describe the phenotypic traits of L. donovani clonal lines that have acquired resistance to HePC. We performed whole-genome and RNA sequencing of these resistant lines to provide an inclusive overview of the multifactorial acquisition of experimental HePC resistance, circumventing the challenge of identifying changes in membrane-bound proteins faced by proteomics. This analysis was complemented by assessment of the in vitro infectivity of HePC-resistant parasites. Our work underscores the importance of complementary "omics" to acquire the most comprehensive insight for multifaceted processes, such as HePC resistance.

Genome of Leptomonas pyrrhocoris: a high-quality reference for monoxenous trypanosomatids and new insights into evolution of Leishmania

Many high-quality genomes are available for dixenous (two hosts) trypanosomatid species of the genera Trypanosoma, Leishmania, and Phytomonas, but only fragmentary information is available for monoxenous (single-host) trypanosomatids. In trypanosomatids, monoxeny is ancestral to dixeny, thus it is anticipated that the genome sequences of the key monoxenous parasites will be instrumental for both understanding the origin of parasitism and the evolution of dixeny. Here, we present a high-quality genome for Leptomonas pyrrhocoris, which is closely related to the dixenous genus Leishmania. The L. pyrrhocoris genome (30.4 Mbp in 60 scaffolds) encodes 10,148 genes. Using the L. pyrrhocoris genome, we pinpointed genes gained in Leishmania. Among those genes, 20 genes with unknown function had expression patterns in the Leishmania mexicana life cycle suggesting their involvement in virulence. By combining differential expression data for L. mexicana, L. major and Leptomonas seymouri, we have identified several additional proteins potentially involved in virulence, including SpoU methylase and U3 small nucleolar ribonucleoprotein IMP3. The population genetics of L. pyrrhocoris was also addressed by sequencing thirteen strains of different geographic origin, allowing the identification of 1,318 genes under positive selection. This set of genes was significantly enriched in components of the cytoskeleton and the flagellum.

Serum Removal from Culture Induces Growth Arrest, Ploidy Alteration, Decrease in Infectivity and Differential Expression of Crucial Genes in Leishmania infantum Promastigotes

Leishmania infantum is one of the species responsible for visceral leishmaniasis. This species is distributed basically in the Mediterranean basin. A recent outbreak in humans has been reported in Spain. Axenic cultures are performed for most procedures with Leishmania spp. promastigotes. This model is stable and reproducible and mimics the conditions of the gut of the sand fly host, which is the natural environment of promastigote development. Culture media are undefined because they contain mammalian serum, which is a rich source of complex lipids and proteins. Serum deprivation slows down the growth kinetics and therefore, yield in biomass. In fact, we have confirmed that the growth rate decreases, as well as infectivity. Ploidy is also affected. Regarding the transcriptome, a high-throughput approach has revealed a low differential expression rate but important differentially regulated genes. The most remarkable profiles are: up-regulation of the GINS Psf3, the fatty acyl-CoA synthase (FAS1), the glyoxylase I (GLO1), the hydrophilic surface protein B (HASPB), the methylmalonyl-CoA epimerase (MMCE) and an amastin gene; and down-regulation of the gPEPCK and the arginase. Implications for metabolic adaptations, differentiation and infectivity are discussed herein.

Kinetoplastid Phylogenomics Reveals the Evolutionary Innovations Associated with the Origins of Parasitism

The evolution of parasitism is a recurrent event in the history of life and a core problem in evolutionary biology. Trypanosomatids are important parasites and include the human pathogens Trypanosoma brucei, Trypanosoma cruzi, and Leishmania spp., which in humans cause African trypanosomiasis, Chagas disease, and leishmaniasis, respectively. Genome comparison between trypanosomatids reveals that these parasites have evolved specialized cell-surface protein families, overlaid on a well-conserved cell template. Understanding how these features evolved and which ones are specifically associated with parasitism requires comparison with related non-parasites. We have produced genome sequences for Bodo saltans, the closest known non-parasitic relative of trypanosomatids, and a second bodonid, Trypanoplasma borreli. Here we show how genomic reduction and innovation contributed to the character of trypanosomatid genomes. We show that gene loss has “streamlined” trypanosomatid genomes, particularly with respect to macromolecular degradation and ion transport, but consistent with a widespread loss of functional redundancy, while adaptive radiations of gene families involved in membrane function provide the principal innovations in trypanosomatid evolution. Gene gain and loss continued during trypanosomatid diversification, resulting in the asymmetric assortment of ancestral characters such as peptidases between Trypanosoma and Leishmania, genomic differences that were subsequently amplified by lineage-specific innovations after divergence. Finally, we show how species-specific, cell-surface gene families (DGF-1 and PSA) with no apparent structural similarity are independent derivations of a common ancestral form, which we call “bodonin.” This new evidence defines the parasitic innovations of trypanosomatid genomes, revealing how a free-living phagotroph became adapted to exploiting hostile host environments.

•

The Bodo saltans genome reveals evolutionary changes at the origin of parasitism

•

Parasite genomes are streamlined, consistent with a loss of functional redundancy

•

Expanded parasite transporter genes reflect a reorientation of membrane function

•

Non-homologous, parasite cell-surface proteins evolved from a common ancestor

Amastin Knockdown in Leishmania braziliensis Affects Parasite-Macrophage Interaction and Results in Impaired Viability of Intracellular Amastigotes

Leishmaniasis, a human parasitic disease with manifestations ranging from cutaneous ulcerations to fatal visceral infection, is caused by several Leishmania species. These protozoan parasites replicate as extracellular, flagellated promastigotes in the gut of a sandfly vector and as amastigotes inside the parasitophorous vacuole of vertebrate host macrophages. Amastins are surface glycoproteins encoded by large gene families present in the genomes of several trypanosomatids and highly expressed in the intracellular amastigote stages of Trypanosoma cruzi and Leishmania spp. Here, we showed that the genome of L. braziliensis contains 52 amastin genes belonging to all four previously described amastin subfamilies and that the expression of members of all subfamilies is upregulated in L. braziliensis amastigotes. Although primary sequence alignments showed no homology to any known protein sequence, homology searches based on secondary structure predictions indicate that amastins are related to claudins, a group of proteins that are components of eukaryotic tight junction complexes. By knocking-down the expression of δ-amastins in L. braziliensis, their essential role during infection became evident. δ-amastin knockdown parasites showed impaired growth after in vitro infection of mouse macrophages and completely failed to produce infection when inoculated in BALB/c mice, an attenuated phenotype that was reverted by the re-expression of an RNAi-resistant amastin gene. Further highlighting their essential role in host-parasite interactions, electron microscopy analyses of macrophages infected with amastin knockdown parasites showed significant alterations in the tight contact that is normally observed between the surface of wild type amastigotes and the membrane of the parasitophorous vacuole.

Leishmaniasis is a parasitic disease caused by more than 20 species of the genus Leishmania that affects about 12 million people throughout the world and for which there is not an effective vaccine. Depending on the Leishmania species, clinical manifestation of the disease varies from self-resolving skin lesions to life-threatening visceralizing diseases. In addition to the toxicity of currently available drugs, their long treatment course, and limited efficacy, a major concern is the development of drug resistant parasite and more virulent variants. Together with the urgent need to develop new drugs that are more effective against this parasite as well as a vaccine to prevent new infections, it is also imperative to develop a better understanding of the factors that determine Leishmania virulence. Here, we describe the characterization of a gene family encoding surface proteins preferentially expressed in the mammalian stage of Leishmania that may be directly involved with the close interaction that is established between the intracellular parasite and host cell membranes. By inhibiting amastin gene expression in L. braziliensis in a mouse model of infection, we showed that these proteins are essential for intracellular parasite survival.

The Leishmania metaphylome: a comprehensive survey of Leishmania protein phylogenetic relationships

Background

Leishmaniasis is a neglected parasitic disease with diverse clinical manifestations and a complex epidemiology. It has been shown that its parasite-related traits vary between species and that they modulate infectivity, pathogenicity, and virulence. However, understanding of the species-specific adaptations responsible for these features and their evolutionary background is limited. To improve our knowledge regarding the parasite biology and adaptation mechanisms of different Leishmania species, we conducted a proteome-wide phylogenomic analysis to gain insights into Leishmania evolution.

Results

The analysis of the reconstructed phylomes (totaling 45,918 phylogenies) allowed us to detect genes that are shared in pathogenic Leishmania species, such as calpain-like cysteine peptidases and 3'a2rel-related proteins, or genes that could be associated with visceral or cutaneous development. This analysis also established the phylogenetic relationship of several hypothetical proteins whose roles remain to be characterized. Our findings demonstrated that gene duplication constitutes an important evolutionary force in Leishmania, acting on protein families that mediate host-parasite interactions, such as amastins, GP63 metallopeptidases, cathepsin L-like proteases, and our methods permitted a deeper analysis of their phylogenetic relationships.

Conclusions

Our results highlight the importance of proteome wide phylogenetic analyses to detect adaptation and evolutionary processes in different organisms and underscore the need to characterize the role of expanded and species-specific proteins in the context of Leishmania evolution by providing a framework for the phylogenetic relationships of Leishmania proteins.

Phylogenomic data are publicly available for use through PhylomeDB (http://www.phylomedb.org).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2091-2) contains supplementary material, which is available to authorized users.