Resequencing of the Leishmania infantum (strain JPCM5) genome and de novo assembly into 36 contigs

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.

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.

Refinement of Leishmania donovani Genome Annotations in the Light of Ribosome-Protected mRNAs Fragments (Ribo-Seq Data)

Advances in next-generation sequencing methodologies have facilitated the assembly of an ever-increasing number of genomes. Gene annotations are typically conducted via specialized software, but the most accurate results require additional manual curation that incorporates insights derived from functional and bioinformatic analyses (e.g., transcriptomics, proteomics, and phylogenetics). In this study, we improved the annotation of the Leishmania donovani (strain HU3) genome using publicly available data from the deep sequencing of ribosome-protected mRNA fragments (Ribo-Seq). As a result of this analysis, we uncovered 70 previously non-annotated protein-coding genes and improved the annotation of around 600 genes. Additionally, we present evidence for small upstream open reading frames (uORFs) in a significant number of transcripts, indicating their potential role in the translational regulation of gene expression. The bioinformatics pipelines developed for these analyses can be used to improve the genome annotations of other organisms for which Ribo-Seq data are available. The improvements provided by these studies will bring us closer to the ultimate goal of a complete and accurately annotated L. donovani genome and will enhance future transcriptomics, proteomics, and genetics studies.

Leishmania infantum (JPCM5) Transcriptome, Gene Models and Resources for an Active Curation of Gene Annotations

Leishmania infantum is one of the causative agents of visceral leishmaniases, the most severe form of leishmaniasis. An improved assembly for the L. infantum genome was published five years ago, yet delineation of its transcriptome remained to be accomplished. In this work, the transcriptome annotation was attained by a combination of both short and long RNA-seq reads. The good agreement between the results derived from both methodologies confirmed that transcript assembly based on Illumina RNA-seq and further delimitation according to the positions of spliced leader (SAS) and poly-A (PAS) addition sites is an adequate strategy to annotate the transcriptomes of Leishmania, a procedure previously used for transcriptome annotation in other Leishmania species and related trypanosomatids. These analyses also confirmed that the Leishmania transcripts boundaries are relatively slippery, showing extensive heterogeneity at the 5′- and 3′-ends. However, the use of RNA-seq reads derived from the PacBio technology (referred to as Iso-Seq) allowed the authors to uncover some complex transcription patterns occurring at particular loci that would be unnoticed by the use of short RNA-seq reads alone. Thus, Iso-Seq analysis provided evidence that transcript processing at particular loci would be more dynamic than expected. Another noticeable finding was the observation of a case of allelic heterozygosity based on the existence of chimeric Iso-Seq reads that might be generated by an event of intrachromosomal recombination. In addition, we are providing the L. infantum gene models, including both UTRs and CDS regions, that would be helpful for undertaking whole-genome expression studies. Moreover, we have built the foundations of a communal database for the active curation of both gene/transcript models and functional annotations for genes and proteins.

Targeting trypanosomes: how chemogenomics and artificial intelligence can guide drug discovery

Trypanosomatids are protozoan parasites that cause human and animal neglected diseases. Despite global efforts, effective treatments are still much needed. Phenotypic screens have provided several chemical leads for drug discovery, but the mechanism of action for many of these chemicals is currently unknown. Recently, chemogenomic screens assessing the susceptibility or resistance of parasites carrying genome-wide modifications started to define the mechanism of action of drugs at large scale. In this review, we discuss how genomics is being used for drug discovery in trypanosomatids, how integration of chemical and genomics data from these and other organisms has guided prioritisations of candidate therapeutic targets and additional chemical starting points, and how these data can fuel the expansion of drug discovery pipelines into the era of artificial intelligence.

An insight into differential protein abundance throughout Leishmania donovani promastigote growth and differentiation

Leishmania donovani causes anthroponotic visceral leishmaniasis, responsible for about 50,000 annual deaths worldwide. Current therapies have considerable side effects. Drug resistance has been reported and no vaccine is available nowadays. The development of undifferentiated promastigotes in the sand fly vector's gut leads to the promastigote form that is highly infective to the mammalian host. Fully differentiated promastigotes play a crucial role in the initial stages of mammalian host infection before internalization in the host phagocytic cell. Therefore, the study of protein levels in the promastigote stage is relevant for disease control, and proteomics analysis is an ideal source of vaccine candidate discovery. This study aims to get insight into the protein levels during the differentiation process of promastigotes by 2DE-MALDI-TOF/TOF. This partial proteome analysis has led to the identification of 75 proteins increased in at least one of the L. donovani promastigote differentiation and growth phases. This study has revealed the differential abundance of said proteins during growth and differentiation. According to previous studies, some are directly involved in parasite survival or are immunostimulatory. The parasite survival-related proteins are ascorbate peroxidase; cystathionine β synthase; an elongation factor 1β paralog; elongation factor 2; endoribonuclease L-PSP; an iron superoxide dismutase paralog; GDP-mannose pyrophosphorylase; several heat shock proteins-HSP70, HSP83-17, mHSP70-rel, HSP110; methylthioadenosine phosphorylase; two thiol-dependent reductase 1 paralogs; transitional endoplasmic reticulum ATPase; and the AhpC thioredoxin paralog. The confirmed immunostimulatory proteins are the heat shock proteins, enolase, and protein kinase C receptor analog. The potential immunostimulatory molecules according to findings in patogenic bacteria are fructose-1,6-diphophate aldolase, dihydrolipoamide acetyltransferase, isocitrate dehydrogenase, pyruvate dehydrogenase E1α and E1β subunits, and triosephosphate isomerase. These proteins may become disease control candidates through future intra-vector control methods or vaccines.

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.

Identification of Leishmania infantum Epidemiology, Drug Resistance and Pathogenicity Biomarkers with Nanopore Sequencing

The emergence of drug-resistant strains of the parasite Leishmania infantum infecting dogs and humans represents an increasing threat. L. infantum genomes are complex and unstable with extensive structural variations, ranging from aneuploidies to multiple copy number variations (CNVs). These CNVs have recently been validated as biomarkers of Leishmania concerning virulence, tissue tropism, and drug resistance. As a proof-of-concept to develop a novel diagnosis platform (LeishGenApp), four L. infantum samples from humans and dogs were nanopore sequenced. Samples were epidemiologically typed within the Mediterranean L. infantum group, identifying members of the JCP5 and non-JCP5 subgroups, using the conserved region (CR) of the maxicircle kinetoplast. Aneuploidies were frequent and heterogenous between samples, yet only chromosome 31 tetrasomy was common between all the samples. A high frequency of aneuploidies was observed for samples with long passage history (MHOM/TN/80/IPT-1), whereas fewer were detected for samples maintained in vivo (MCRI/ES/2006/CATB033). Twenty-two genes were studied to generate a genetic pharmacoresistance profile against miltefosine, allopurinol, trivalent antimonials, amphotericin, and paromomycin. MHOM/TN/80/IPT-1 and MCRI/ES/2006/CATB033 displayed a genetic profile with potential resistance against miltefosine and allopurinol. Meanwhile, MHOM/ES/2016/CATB101 and LCAN/ES/2020/CATB102 were identified as potentially resistant against paromomycin. All four samples displayed a genetic profile for resistance against trivalent antimonials. Overall, this proof-of-concept revealed the potential of nanopore sequencing and LeishGenApp for the determination of epidemiological, drug resistance, and pathogenicity biomarkers in L. infantum.

Transcriptome Analysis of Intracellular Amastigotes of Clinical Leishmania infantum Lines from Therapeutic Failure Patients after Infection of Human Macrophages

Leishmaniasis is considered to be one of the most neglected tropical diseases affecting humans and animals around the world. Due to the absence of an effective vaccine, current treatment is based on chemotherapy. However, the continuous appearance of drug resistance and therapeutic failure (TF) lead to an early obsolescence of treatments. Identification of the factors that contribute to TF and drug resistance in leishmaniasis will constitute a useful tool for establishing future strategies to control this disease. In this manuscript, we evaluated the transcriptomic changes in the intracellular amastigotes of the Leishmania infantum parasites isolated from patients with leishmaniasis and TF at 96 h post-infection of THP-1 cells. The adaptation of the parasites to their new environment leads to expression alterations in the genes involved mainly in the transport through cell membranes, energy and redox metabolism, and detoxification. Specifically, the gene that codes for the prostaglandin f2α synthase seems to be relevant in the pathogenicity and TF since it appears substantially upregulated in all the L. infantum lines. Overall, our results show that at the late infection timepoint, the transcriptome of the parasites undergoes significant changes that probably improve the survival of the Leishmania lines in the host cells, contributing to the TF phenotype as well as drug therapy evasion.

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

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

The Astonishing Large Family of HSP40/DnaJ Proteins Existing in Leishmania

Abrupt environmental changes are faced by Leishmania parasites during transmission from a poikilothermic insect vector to a warm-blooded host. Adaptation to harsh environmental conditions, such as nutrient deprivation, hypoxia, oxidative stress and heat shock needs to be accomplished by rapid reconfiguration of gene expression and remodeling of protein interaction networks. Chaperones play a central role in the maintenance of cellular homeostasis, and they are responsible for crucial tasks such as correct folding of nascent proteins, protein translocation across different subcellular compartments, avoiding protein aggregates and elimination of damaged proteins. Nearly one percent of the gene content in the Leishmania genome corresponds to members of the HSP40 family, a group of proteins that assist HSP70s in a variety of cellular functions. Despite their expected relevance in the parasite biology and infectivity, little is known about their functions or partnership with the different Leishmania HSP70s. Here, we summarize the structural features of the 72 HSP40 proteins encoded in the Leishmania infantum genome and their classification into four categories. A review of proteomic data, together with orthology analyses, allow us to postulate cellular locations and possible functional roles for some of them. A detailed study of the members of this family would provide valuable information and opportunities for drug discovery and improvement of current treatments against leishmaniasis.

In-Depth Quantitative Proteomics Characterization of In Vitro Selected Miltefosine Resistance in Leishmania infantum

Visceral leishmaniasis (VL) is a neglected disease caused by Leishmania parasites. Although significant morbidity and mortality in tropical and subtropical regions of the world are associated with VL, the low investment for developing new treatment measures is chronic. Moreover, resistance and treatment failure are increasing for the main medications, but the emergence of resistance phenotypes is poorly understood at the protein level. Here, we analyzed the development of resistance to miltefosine upon experimental selection in a L. infantum strain. Time to miltefosine resistance emergence was ~six months and label-free quantitative mass-spectrometry-based proteomics analyses revealed that this process involves a remodeling of components of the membrane and mitochondrion, with significant increase in oxidative phosphorylation complexes, particularly on complex IV and ATP synthase, accompanied by increased energy metabolism mainly dependent on β-oxidation of fatty acids. Proteins canonically involved in ROS detoxification did not contribute to the resistant process whereas sterol biosynthesis enzymes could have a role in this development. Furthermore, changes in the abundance of proteins known to be involved in miltefosine resistance such as ABC transporters and phospholipid transport ATPase were detected. Together, our data show a more complete picture of the elements that make up the miltefosine resistance phenotype in L. infantum.

ARAMIS: From systematic errors of NGS long reads to accurate assemblies

NGS long-reads sequencing technologies (or third generation) such as Pacific BioSciences (PacBio) have revolutionized the sequencing field over the last decade improving multiple genomic applications like de novo genome assemblies. However, their error rate, mostly involving insertions and deletions (indels), is currently an important concern that requires special attention to be solved. Multiple algorithms are available to fix these sequencing errors using short reads (such as Illumina), although they require long processing times and some errors may persist. Here, we present Accurate long-Reads Assembly correction Method for Indel errorS (ARAMIS), the first NGS long-reads indels correction pipeline that combines several correction software in just one step using accurate short reads. As a proof OF concept, six organisms were selected based on their different GC content, size and genome complexity, and their PacBio-assembled genomes were corrected thoroughly by this pipeline. We found that the presence of systematic sequencing errors in long-reads PacBio sequences affecting homopolymeric regions, and that the type of indel error introduced during PacBio sequencing are related to the GC content of the organism. The lack of knowledge of this fact leads to the existence of numerous published studies where such errors have been found and should be resolved since they may contain incorrect biological information. ARAMIS yields better results with less computational resources needed than other correction tools and gives the possibility of detecting the nature of the found indel errors found and its distribution along the genome. The source code of ARAMIS is available at https://github.com/genomics-ngsCBMSO/ARAMIS.git.

Revisiting the heterogeneous global genomic population structure of Leishmania infantum

Leishmania infantum is the main causative agent responsible for visceral leishmaniasis (VL), a disease with global distribution. The genomic structure and genetic variation of this species have been widely studied in different parts of the world. However, in some countries, this information is still yet unknown, as is the genomic behaviour of the main antigens used in VL diagnosis (rK39 and rK28), which have demonstrated variable sensitivity and specificity in a manner dependent on the geographic region analysed. The objective of this study was to explore the genomic architecture and diversity of four Colombian L. infantum isolates obtained in this study and to compare these results with the genetic analysis of 183 L. infantum isolates from across the world (obtained from public databases), as well as to analyse the whole rK39 and rK28 antigen sequences in our dataset. The results showed that, at the global level, L. infantum has high genetic homogeneity and extensive aneuploidy. Furthermore, we demonstrated that there are distinct populations of L. infantum circulating in various countries throughout the globe and that populations of distant countries have close genomic relationships. Additionally, this study demonstrated the high genetic variability of the rK28 antigen worldwide. In conclusion, our study allowed us to (i) expand our knowledge of the genomic structure of global L. infantum; (ii) describe the intra-specific genomic variability of this species; and (iii) understand the genomic characteristics of the main antigens used in the diagnosis of VL. Additionally, this is the first study to report whole-genome sequences of Colombian L. infantum isolates.

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.

Gene Annotation and Transcriptome Delineation on a De Novo Genome Assembly for the Reference Leishmania major Friedlin Strain

Leishmania major is the main causative agent of cutaneous leishmaniasis in humans. The Friedlin strain of this species (LmjF) was chosen when a multi-laboratory consortium undertook the objective of deciphering the first genome sequence for a parasite of the genus Leishmania. The objective was successfully attained in 2005, and this represented a milestone for Leishmania molecular biology studies around the world. Although the LmjF genome sequence was done following a shotgun strategy and using classical Sanger sequencing, the results were excellent, and this genome assembly served as the reference for subsequent genome assemblies in other Leishmania species. Here, we present a new assembly for the genome of this strain (named LMJFC for clarity), generated by the combination of two high throughput sequencing platforms, Illumina short-read sequencing and PacBio Single Molecular Real-Time (SMRT) sequencing, which provides long-read sequences. Apart from resolving uncertain nucleotide positions, several genomic regions were reorganized and a more precise composition of tandemly repeated gene loci was attained. Additionally, the genome annotation was improved by adding 542 genes and more accurate coding-sequences defined for around two hundred genes, based on the transcriptome delimitation also carried out in this work. As a result, we are providing gene models (including untranslated regions and introns) for 11,238 genes. Genomic information ultimately determines the biology of every organism; therefore, our understanding of molecular mechanisms will depend on the availability of precise genome sequences and accurate gene annotations. In this regard, this work is providing an improved genome sequence and updated transcriptome annotations for the reference L. major Friedlin strain.

Proteins interacting with Leishmania major PUF1: A proteomic dataset

Leishmania parasites must deal with stressful environmental conditions (thermal, nutritional and oxidative) along their digenetic life cycles. This requires drastic changes in gene expression, which in this parasite occurs mainly through post-transcriptional mechanisms involving RNA binding proteins (RBPs). PUF proteins, a class of RBPs existing in most eukaryotic organisms, might play too an essential role modulating the fate of mRNAs and regulating gene expression in Leishmania parasites. A proteomic approach to identify putative protein partners (interactome) of the Leishmania major PUF1 protein was performed. The PUF1 interactome was characterized by co-immunoprecipitation using L. major cellular extracts and an anti-LiPUF1 antibody, and a subsequent analysis of the co-immunoprecipitated proteins by mass-spectrometry, identifying 90 LmPUF1 candidate partners. Remarkably, many of the identified proteins are other RBPs and/or putative P bodies and mRNA-exporting machinery components. Raw mass spectrometry data are available via ProteomeXchange with identifier PXD022581.

The In Silico Identification of Potential Members of the Ded1/DDX3 Subfamily of DEAD-Box RNA Helicases from the Protozoan Parasite Leishmania infantum and Their Analyses in Yeast

DEAD-box RNA helicases are ubiquitous proteins found in all kingdoms of life and that are associated with all processes involving RNA. Their central roles in biology make these proteins potential targets for therapeutic or prophylactic drugs. The Ded1/DDX3 subfamily of DEAD-box proteins is of particular interest because of their important role(s) in translation. In this paper, we identified and aligned the protein sequences of 28 different DEAD-box proteins from the kinetoplast-protozoan parasite Leishmania infantum, which is the cause of the visceral form of leishmaniasis that is often lethal if left untreated, and compared them with the consensus sequence derived from DEAD-box proteins in general, and from the Ded1/DDX3 subfamily in particular, from a wide variety of other organisms. We identified three potential homologs of the Ded1/DDX3 subfamily and the equivalent proteins from the related protozoan parasite Trypanosoma brucei, which is the causative agent of sleeping sickness. We subsequently tested these proteins for their ability to complement a yeast strain deleted for the essential DED1 gene. We found that the DEAD-box proteins from Trypanosomatids are highly divergent from other eukaryotes, and consequently they are suitable targets for protein-specific drugs.

Comparative transcriptomic analysis of antimony resistant and susceptible Leishmania infantum lines

Background

One of the major challenges to leishmaniasis treatment is the emergence of parasites resistant to antimony. To study differentially expressed genes associated with drug resistance, we performed a comparative transcriptomic analysis between wild-type and potassium antimonyl tartrate (Sb^III)-resistant Leishmania infantum lines using high-throughput RNA sequencing.

Methods

All the cDNA libraries were constructed from promastigote forms of each line, sequenced and analyzed using STAR for mapping the reads against the reference genome (L. infantum JPCM5) and DESeq2 for differential expression statistical analyses. All the genes were functionally annotated using sequence similarity search.

Results

The analytical pipeline considering an adjusted p-value < 0.05 and fold change > 2.0 identified 933 transcripts differentially expressed (DE) between wild-type and Sb^III-resistant L. infantum lines. Out of 933 DE transcripts, 504 presented functional annotation and 429 were assigned as hypothetical proteins. A total of 837 transcripts were upregulated and 96 were downregulated in the Sb^III-resistant L. infantum line. Using this DE dataset, the proteins were further grouped in functional classes according to the gene ontology database. The functional enrichment analysis for biological processes showed that the upregulated transcripts in the Sb^III-resistant line are associated with protein phosphorylation, microtubule-based movement, ubiquitination, host–parasite interaction, cellular process and other categories. The downregulated transcripts in the Sb^III-resistant line are assigned in the GO categories: ribonucleoprotein complex, ribosome biogenesis, rRNA processing, nucleosome assembly and translation.

Conclusions

The transcriptomic profile of L. infantum showed a robust set of genes from different metabolic pathways associated with the antimony resistance phenotype in this parasite. Our results address the complex and multifactorial antimony resistance mechanisms in Leishmania, identifying several candidate genes that may be further evaluated as molecular targets for chemotherapy of leishmaniasis.

The Experimental Proteome of Leishmania infantum Promastigote and Its Usefulness for Improving Gene Annotations

Leishmania infantum causes visceral leishmaniasis (kala-azar), the most severe form of leishmaniasis, which is lethal if untreated. A few years ago, the re-sequencing and de novo assembling of the L. infantum (JPCM5 strain) genome was accomplished, and now we aimed to describe and characterize the experimental proteome of this species. In this work, we performed a proteomic analysis from axenic cultured promastigotes and carried out a detailed comparison with other Leishmania experimental proteomes published to date. We identified 2352 proteins based on a search of mass spectrometry data against a database built from the six-frame translated genome sequence of L. infantum. We detected many proteins belonging to organelles such as glycosomes, mitochondria, or flagellum, as well as many metabolic enzymes and many putative RNA binding proteins and molecular chaperones. Moreover, we listed some proteins presenting post-translational modifications, such as phosphorylations, acetylations, and methylations. On the other hand, the identification of peptides mapping to genomic regions previously annotated as non-coding allowed for the correction of annotations, leading to the N-terminal extension of protein sequences and the uncovering of eight novel protein-coding genes. The alliance of proteomics, genomics, and transcriptomics has resulted in a powerful combination for improving the annotation of the L. infantum reference genome.

LPG2 Gene Duplication in Leishmania infantum: A Case for CRISPR-Cas9 Gene Editing

On the surface of the Leishmania promastigote, phosphoglycans (PG) such as lipophosphoglycan (LPG), proteophosphoglycan (PPG), free phosphoglycan polymers (PGs), and acid phosphatases (sAP), are dominant and contribute to the invasion and survival of Leishmania within the host cell by modulating macrophage signaling and intracellular trafficking. Phosphoglycan synthesis depends on the Golgi GDP-mannose transporter encoded by the LPG2 gene. Aiming to investigate the role of PG-containing molecules in Leishmania infantum infection process, herein we describe the generation and characterization of L. infantum LPG2-deficient parasites. This gene was unexpectedly identified as duplicated in the L. infantum genome, which impaired gene targeting using the conventional homologous recombination approach. This limitation was circumvented by the use of CRISPR/Cas9 technology. Knockout parasites were selected by agglutination assays using CA7AE antibodies followed by a lectin (RCA 120). Five clones were isolated and molecularly characterized, all revealing the expected edited genome, as well as the complete absence of LPG and PG-containing molecule expression. Finally, the deletion of LPG2 was found to impair the outcome of infection in human neutrophils, as demonstrated by a pronounced reduction (~83%) in intracellular load compared to wild-type parasite infection. The results obtained herein reinforce the importance of LPG and other PGs as virulence factors in host-parasite interactions.

Application of next generation sequencing (NGS) for descriptive analysis of 30 genomes of Leishmania infantum isolates in Middle-North Brazil

Visceral leishmaniasis (VL) is a life-threatening disease caused by the protozoa Leishmania donovani and L. infantum. Likely, L. infantum was introduced in the New World by the Iberic colonizers. Due to recent introduction, the genetic diversity is low. Access to genomic information through the sequencing of Leishmania isolates allows the characterization of populations through the identification and analysis of variations. Population structure information may reveal important data on disease dynamics. Aiming to describe the genetic diversity of L. infantum from the Middle-North, Brazil, next generation sequencing of 30 Leishmania isolates obtained in the city of Teresina, from where the disease dispersed, was performed. The variations were categorized accordingly to the genome region and impact and provided the basis for chromosomal ploidy and population structure analysis. The results showed low diversity between the isolates and the Iberic reference genome JPCM5. Most variations were seen in non-coding regions, with modifying impact. The ploidy number analysis showed aneuploid profile. The population structure analysis revealed the presence of two L. infantum populations identified in Teresina. Further population genetics studies with a larger number of isolates should be performed in order to identify the genetic background associated with virulence and parasite ecology.

A review of the systematics, species identification and diagnostics of the Trypanosomatidae using the maxicircle kinetoplast DNA: from past to present

The Trypanosomatid family are a diverse and widespread group of protozoan parasites that belong to the higher order class Kinetoplastida. Containing predominantly monoxenous species (i.e. those having only a single host) that are confined to invertebrate hosts, this class is primarily known for its pathogenic dixenous species (i.e. those that have two hosts), serving as the aetiological agents of the important neglected tropical diseases including leishmaniasis, American trypanosomiasis (Chagas disease) and human African trypanosomiasis. Over the past few decades, a multitude of studies have investigated the diversity, classification and evolutionary history of the trypanosomatid family using different approaches and molecular targets. The mitochondrial-like DNA of the trypanosomatid parasites, also known as the kinetoplast, has emerged as a unique taxonomic and diagnostic target for exploring the evolution of this diverse group of parasitic eukaryotes. This review discusses recent advancements and important developments that have made a significant impact in the field of trypanosomatid systematics and diagnostics in recent years.

Seroprevalence, Clinical, and Pathological Characteristics of Canine Leishmaniasis in a Central Region of Colombia

Introduction

Leishmaniasis is a zoonotic disease which is caused by protozoan parasites of the genus Leishmania. Canids are the most important reservoir of the parasites; however, limited data are available on the species of Leishmania prevalent in these animals and their impact on human health. The objective of this study was to estimate the seroprevalence of leishmaniasis in dogs from an inter-Andean region of Colombia during July 2016–July 2017, and to describe the clinical and histopathological features of the disease.

Material and Methods

A total of 155 dogs were subjected to clinical examination and a serological test for detection of antibodies against Leishmania. Necropsy was carried out on positive animals and tissue samples were processed by routine histopathology.

Results

Altogether 19 dogs were positive in the serological test, establishing a 12% seroprevalence of Leishmania. Clinical examination and necropsy revealed exfoliative and ulcerative dermatitis with haemorrhagic borders on the ears, head, nose, and legs. Histopathology revealed severe multifocal dermatitis with abundant Leishmania amastigotes within the cytoplasm of phagocytic cells, depletion of lymphocytes in lymphoid tissues, interstitial pneumonia, and interstitial nephritis. Tissue samples were positive for Leishmania by PCR.

Conclusion

The macro- and microscopic changes correlated with the presence of Leishmania as established by serological test and PCR.

Proteins involved in the biosynthesis of lipophosphoglycan in Leishmania: a comparative genomic and evolutionary analysis

BACKGROUND

Leishmania spp. are digenetic parasites capable of infecting humans and causing a range of diseases collectively known as leishmaniasis. The main mechanisms involved in the development and permanence of this pathology are linked to evasion of the immune response. Crosstalk between the immune system and particularities of each pathogenic species is associated with diverse disease manifestations. Lipophosphoglycan (LPG), one of the most important molecules present on the surface of Leishmania parasites, is divided into four regions with high molecular variability. Although LPG plays an important role in host-pathogen and vector-parasite interactions, the distribution and phylogenetic relatedness of the genes responsible for its synthesis remain poorly explored. The recent availability of full genomes and transcriptomes of Leishmania parasites offers an opportunity to leverage insight on how LPG-related genes are distributed and expressed by these pathogens.

RESULTS

Using a phylogenomics-based framework, we identified a catalog of genes involved in LPG biosynthesis across 22 species of Leishmania from the subgenera Viannia and Leishmania, as well as 5 non-Leishmania trypanosomatids. The evolutionary relationships of these genes across species were also evaluated. Nine genes related to the production of the glycosylphosphatidylinositol (GPI)-anchor were highly conserved among compared species, whereas 22 genes related to the synthesis of the repeat unit presented variable conservation. Extensive gain/loss events were verified, particularly in genes SCG1-4 and SCA1-2. These genes act, respectively, on the synthesis of the side chain attached to phosphoglycans and in the transfer of arabinose residues. Phylogenetic analyses disclosed evolutionary patterns reflective of differences in host specialization, geographic origin and disease manifestation.

CONCLUSIONS

The multiple gene gain/loss events identified by genomic data mining help to explain some of the observed intra- and interspecies variation in LPG structure. Collectively, our results provide a comprehensive catalog that details how LPG-related genes evolved in the Leishmania parasite specialization process.

Detecting sequence variants in clinically important protozoan parasites

Second and third generation sequencing methods are crucial for population genetic studies, and variant detection is a popular approach for exploiting this sequence data. While mini- and microsatellites are historically useful markers for studying important Protozoa such as Toxoplasma and Plasmodium spp., detecting non-repetitive variants such as those found in genes can be fundamental to investigating a pathogen's biology. These variants, namely single nucleotide polymorphisms and insertions and deletions, can help elucidate the genetic basis of an organism's pathogenicity, identify selective pressures, and resolve phylogenetic relationships. They also have the added benefit of possessing a comparatively low mutation rate, which contributes to their stability. However, there is a plethora of variant analysis tools with nuanced pipelines and conflicting recommendations for best practise, which can be confounding. This lack of standardisation means that variant analysis requires careful parameter optimisation, an understanding of its limitations, and the availability of high quality data. This review explores the value of variant detection when applied to non-model organisms such as clinically important protozoan pathogens. The limitations of current methods are discussed, including special considerations that require the end-users' attention to ensure that the results generated are reproducible, and the biological conclusions drawn are valid.

Leishmania Mitochondrial Genomes: Maxicircle Structure and Heterogeneity of Minicircles

The mitochondrial DNA (mtDNA), which is present in almost all eukaryotic organisms, is a useful marker for phylogenetic studies due to its relative high conservation and its inheritance manner. In Leishmania and other trypanosomatids, the mtDNA (also referred to as kinetoplast DNA or kDNA) is composed of thousands of minicircles and a few maxicircles, catenated together into a complex network. Maxicircles are functionally similar to other eukaryotic mtDNAs, whereas minicircles are involved in RNA editing of some maxicircle-encoded transcripts. Next-generation sequencing (NGS) is increasingly used for assembling nuclear genomes and, currently, a large number of genomic sequences are available. However, most of the time, the mitochondrial genome is ignored in the genome assembly processes. The aim of this study was to develop a pipeline to assemble Leishmania minicircles and maxicircle DNA molecules, exploiting the raw data generated in the NGS projects. As a result, the maxicircle molecules and the plethora of minicircle classes for Leishmania major, Leishmania infantum and Leishmania braziliensis have been characterized. We have observed that whereas the heterogeneity of minicircle sequences existing in a single cell hampers their use for Leishmania typing and classification, maxicircles emerge as an extremely robust genetic marker for taxonomic studies within the clade of kinetoplastids.

Single-Strand Annealing Plays a Major Role in Double-Strand DNA Break Repair following CRISPR-Cas9 Cleavage in Leishmania

CRISPR-Cas9 genome editing relies on an efficient double-strand DNA break (DSB) and repair. Contrary to mammalian cells, the protozoan parasite Leishmania lacks the most efficient nonhomologous end-joining pathway and uses microhomology-mediated end joining (MMEJ) and, occasionally, homology-directed repair to repair DSBs. Here, we reveal that Leishmania predominantly uses single-strand annealing (SSA) (>90%) instead of MMEJ (<10%) for DSB repair (DSBR) following CRISPR targeting of the miltefosine transporter gene, resulting in 9-, 18-, 20-, and 29-kb sequence deletions and multiple gene codeletions. Strikingly, when targeting the Leishmania donovani LdBPK_241510 gene, SSA even occurred by using direct repeats 77 kb apart, resulting in the codeletion of 15 Leishmania genes, though with a reduced frequency. These data strongly indicate that DSBR is not efficient in Leishmania, which explains why more than half of DSBs led to cell death and why the CRISPR gene-targeting efficiency is low compared with that in other organisms. Since direct repeat sequences are widely distributed in the Leishmania genome, we predict that many DSBs created by CRISPR are repaired by SSA. It is also revealed that DNA polymerase theta is involved in both MMEJ and SSA in Leishmania Collectively, this study establishes that DSBR mechanisms and their competence in an organism play an important role in determining the outcome and efficacy of CRISPR gene targeting. These observations emphasize the use of donor DNA templates to improve gene editing specificity and efficiency in Leishmania In addition, we developed a novel Staphylococcus aureus Cas9 constitutive expression vector (pLdSaCN) for gene targeting in Leishmania IMPORTANCE Due to differences in double-strand DNA break (DSB) repair mechanisms, CRISPR-Cas9 gene editing efficiency can vary greatly in different organisms. In contrast to mammalian cells, the protozoan parasite Leishmania uses microhomology-mediated end joining (MMEJ) and, occasionally, homology-directed repair (HDR) to repair DSBs but lacks the nonhomologous end-joining pathway. Here, we show that Leishmania predominantly uses single-strand annealing (SSA) instead of MMEJ for DSB repairs (DSBR), resulting in large deletions that can include multiple genes. This strongly indicates that the overall DSBR in Leishmania is inefficient and therefore can influence the outcome of CRISPR-Cas9 gene editing, highlighting the importance of using a donor DNA to improve gene editing fidelity and efficiency in Leishmania.

Genome and transcriptome analyses of Leishmania spp.: opening Pandora's box

In the last 30 years, significant advances in genetic manipulation tools along with complete genome and transcriptome sequencing have advanced our understanding of the biology of Leishmania parasites and their interplay with the sand fly and mammalian hosts. High-throughput sequencing in association with CRISPR/Cas9 have prepared the ground for significant advances. Given the richness of the progress made over the last decade, in this article, we focused on the most recent contributions of genome-wide and transcriptome analyses of Leishmania spp., which permit the comparison of life cycle stages, the evaluation of different strains and species in their natural niches and in the field and the simultaneously comparison of the gene expression profiles of parasites and hosts.

Complete assembly of the Leishmania donovani (HU3 strain) genome and transcriptome annotation

Leishmania donovani is a unicellular parasite that causes visceral leishmaniasis, a fatal disease in humans. In this study, a complete assembly of the genome of L. donovani is provided. Apart from being the first published genome of this strain (HU3), this constitutes the best assembly for an L. donovani genome attained to date. The use of a combination of sequencing platforms enabled to assemble, without any sequence gap, the 36 chromosomes for this species. Additionally, based on this assembly and using RNA-seq reads derived from poly-A + RNA, the transcriptome for this species, not yet available, was delineated. Alternative SL addition sites and heterogeneity in the poly-A addition sites were commonly observed for most of the genes. After a complete annotation of the transcriptome, 2,410 novel transcripts were defined. Additionally, the relative expression for all transcripts present in the promastigote stage was determined. Events of cis-splicing have been documented to occur during the maturation of the transcripts derived from genes LDHU3_07.0430 and LDHU3_29.3990. The complete genome assembly and the availability of the gene models (including annotation of untranslated regions) are important pieces to understand how differential gene expression occurs in this pathogen, and to decipher phenotypic peculiarities like tissue tropism, clinical disease, and drug susceptibility.

The Tritryps Comparative Repeatome: Insights on Repetitive Element Evolution in Trypanosomatid Pathogens

The major human pathogens Trypanosoma cruzi, Trypanosoma brucei, and Leishmania major are collectively known as the Tritryps. The initial comparative analysis of their genomes has uncovered that Tritryps share a great number of genes, but repetitive DNA seems to be extremely variable between them. However, the in-depth characterization of repetitive DNA in these pathogens has been in part neglected, mainly due to the well-known technical challenges of studying repetitive sequences from de novo assemblies using short reads. Here, we compared the repetitive DNA repertories between the Tritryps genomes using genome-wide, low-coverage Illumina sequencing coupled to RepeatExplorer analysis. Our work demonstrates that this extensively implemented approach for studying higher eukaryote repeatomes is also useful for protozoan parasites like trypanosomatids, as we recovered previously observed differences in the presence and amount of repetitive DNA families. Additionally, our estimations of repetitive DNA abundance were comparable to those obtained from enhanced-quality assemblies using longer reads. Importantly, our methodology allowed us to describe a previously undescribed transposable element in Leishmania major (TATE element), highlighting its potential to accurately recover distinctive features from poorly characterized repeatomes. Together, our results support the application of this low-cost, low-coverage sequencing approach for the extensive characterization of repetitive DNA evolutionary dynamics in trypanosomatid and other protozoan genomes.

In silico characterization of multiple genes encoding the GP63 virulence protein from Leishmania braziliensis: identification of sources of variation and putative roles in immune evasion

Background

The leishmaniasis are parasitic diseases caused by protozoans of the genus Leishmania, highly divergent eukaryotes, characterized by unique biological features. To survive in both the mammalian hosts and insect vectors, these pathogens make use of a number of mechanisms, many of which are associated with parasite specific proteases. The metalloprotease GP63, the major Leishmania surface antigen, has been found to have multiple functions required for the parasite’s survival. GP63 is encoded by multiple genes and their copy numbers vary considerably between different species and are increased in those from the subgenus Viannia, including L. braziliensis.

Results

By comparing multiple sequences from Leishmania and related organisms this study sought to characterize paralogs in silico, evaluating their differences and similarities and the implications for the GP63 function. The Leishmania GP63 genes are encoded on chromosomes 10, 28 and 31, with the genes from the latter two chromosomes more related to genes found in insect or plant parasites. Those from chromosome 10 have experienced independent expansions in numbers in Leishmania, especially in L. braziliensis. These could be clustered in three groups associated with different mRNA 3′ untranslated regions as well as distinct C-terminal ends for the encoded proteins, with presumably distinct expression patterns and subcellular localizations. Sequence variations between the chromosome 10 genes were linked to intragenic recombination events, mapped to the external surface of the proteins and predicted to be immunogenic, implying a role against the host immune response.

Conclusions

Our results suggest a greater role for the sequence variation found among the chromosome 10 GP63 genes, possibly related to the pathogenesis of L. braziliensis and closely related species within the mammalian host. They also indicate different functions associated to genes mapped to different chromosomes. For the chromosome 10 genes, variable subcellular localizations were found to be most likely associated with multiple functions and target substrates for this versatile protease.

Electronic supplementary material

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

A Guide to Next Generation Sequence Analysis of Leishmania Genomes

Next generation sequencing (NGS) technology transformed Leishmania genome studies and became an indispensable tool for Leishmania researchers. Recent Leishmania genomics analyses facilitated the discovery of various genetic diversities including single nucleotide polymorphisms (SNPs), copy number variations (CNVs), somy variations, and structural variations in detail and provided valuable insights into the complexity of the genome and gene regulation. Many aspects of Leishmania NGS analyses are similar to those of related pathogens like trypanosomes. However, the analyses of Leishmania genomes face a unique challenge because of the presence of frequent aneuploidy. This makes characterization and interpretation of read depth and somy a key part of Leishmania NGS analyses because read depth affects the accuracy of detection of all genetic variations. However, there are no general guidelines on how to explore and interpret the impact of aneuploidy, and this has made it difficult for biologists and bioinformaticians, especially for beginners, to perform their own analyses and interpret results across different analyses. In this guide we discuss a wide range of topics essential for Leishmania NGS analyses, ranging from how to set up a computational environment for genome analyses, to how to characterize genetic variations among Leishmania samples, and we will particularly focus on chromosomal copy number variation and its impact on genome analyses.

Complete and de novo assembly of the Leishmania braziliensis (M2904) genome

Leishmania braziliensis is the etiological agent of American mucosal leishmaniasis, one of the most severe clinical forms of leishmaniasis. Here, we report the assembly of the L. braziliensis (M2904) genome into 35 continuous chromosomes. Also, the annotation of 8395 genes is provided. The public availability of this information will contribute to a better knowledge of this pathogen and help in the search for vaccines and novel drug targets aimed to control the disease caused by this Leishmania species.

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.

A complete Leishmania donovani reference genome identifies novel genetic variations associated with virulence

Leishmania donovani is responsible for visceral leishmaniasis, a neglected and lethal parasitic disease with limited treatment options and no vaccine. The study of L. donovani has been hindered by the lack of a high-quality reference genome and this can impact experimental outcomes including the identification of virulence genes, drug targets and vaccine development. We therefore generated a complete genome assembly by deep sequencing using a combination of second generation (Illumina) and third generation (PacBio) sequencing technologies. Compared to the current L. donovani assembly, the genome assembly reported within resulted in the closure over 2,000 gaps, the extension of several chromosomes up to telomeric repeats and the re-annotation of close to 15% of protein coding genes and the annotation of hundreds of non-coding RNA genes. It was possible to correctly assemble the highly repetitive A2 and Amastin virulence gene clusters. A comparative sequence analysis using the improved reference genome confirmed 70 published and identified 15 novel genomic differences between closely related visceral and atypical cutaneous disease-causing L. donovani strains providing a more complete map of genes associated with virulence and visceral organ tropism. Bioinformatic tools including protein variation effect analyzer and basic local alignment search tool were used to prioritize a list of potential virulence genes based on mutation severity, gene conservation and function. This complete genome assembly and novel information on virulence factors will support the identification of new drug targets and the development of a vaccine for L. donovani.

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.