Recent advances in trypanosomatid research: genome organization, expression, metabolism, taxonomy and evolution

Blastocrithidia nonstop mitochondrial genome and its expression are remarkably insulated from nuclear codon reassignment

The canonical stop codons of the nuclear genome of the trypanosomatid Blastocrithidia nonstop are recoded. Here, we investigated the effect of this recoding on the mitochondrial genome and gene expression. Trypanosomatids possess a single mitochondrion and protein-coding transcripts of this genome require RNA editing in order to generate open reading frames of many transcripts encoded as 'cryptogenes'. Small RNAs that can number in the hundreds direct editing and produce a mitochondrial transcriptome of unusual complexity. We find B. nonstop to have a typical trypanosomatid mitochondrial genetic code, which presumably requires the mitochondrion to disable utilization of the two nucleus-encoded suppressor tRNAs, which appear to be imported into the organelle. Alterations of the protein factors responsible for mRNA editing were also documented, but they have likely originated from sources other than B. nonstop nuclear genome recoding. The population of guide RNAs directing editing is minimal, yet virtually all genes for the plethora of known editing factors are still present. Most intriguingly, despite lacking complex I cryptogene guide RNAs, these cryptogene transcripts are stochastically edited to high levels.

Cultivation of monoxenous trypanosomatids: A minireview

Trypanosomatids are obligatory parasites, some of which are responsible for important human and animal diseases, but the vast majority of trypanosomatids are restricted to invertebrate hosts. Isolation and in vitro cultivation of trypanosomatids from insect hosts enable their description, characterization, and subsequently genetic and genomic studies. However, exact nutritional requirements are still unknown for most trypanosomatids and thus very few defined media are available. This mini review provides information about the role of different ingredients, recommendations and advice on essential supplements and important physicochemical parameters of culture media with the aim of facilitating first attempts to cultivate insect-infesting trypanosomatids, with a focus on monoxenous trypanosomatids.

In silico prediction of the metabolism of Blastocrithidia nonstop, a trypanosomatid with non-canonical genetic code

BACKGROUND

Almost all extant organisms use the same, so-called canonical, genetic code with departures from it being very rare. Even more exceptional are the instances when a eukaryote with non-canonical code can be easily cultivated and has its whole genome and transcriptome sequenced. This is the case of Blastocrithidia nonstop, a trypanosomatid flagellate that reassigned all three stop codons to encode amino acids.

RESULTS

We in silico predicted the metabolism of B. nonstop and compared it with that of the well-studied human parasites Trypanosoma brucei and Leishmania major. The mapped mitochondrial, glycosomal and cytosolic metabolism contains all typical features of these diverse and important parasites. We also provided experimental validation for some of the predicted observations, concerning, specifically presence of glycosomes, cellular respiration, and assembly of the respiratory complexes.

CONCLUSIONS

In an unusual comparison of metabolism between a parasitic protist with a massively altered genetic code and its close relatives that rely on a canonical code we showed that the dramatic differences on the level of nucleic acids do not seem to be reflected in the metabolisms. Moreover, although the genome of B. nonstop is extremely AT-rich, we could not find any alterations of its pyrimidine synthesis pathway when compared to other trypanosomatids. Hence, we conclude that the dramatic alteration of the genetic code of B. nonstop has no significant repercussions on the metabolism of this flagellate.

Phylogenetic framework to explore trait evolution in Trypanosomatidae

The number of sequenced trypanosomatid genomes has reached a critical point so that they are now available for almost all genera and subgenera. Based on this, we inferred a phylogenomic tree and propose it as a framework to study trait evolution together with some examples of how to do it.

Linear motifs regulating protein secretion, sorting and autophagy in Leishmania parasites are diverged with respect to their host equivalents

The pathogenic, tropical Leishmania flagellates belong to an early-branching eukaryotic lineage (Kinetoplastida) with several unique features. Unfortunately, they are poorly understood from a molecular biology perspective, making development of mechanistically novel and selective drugs difficult. Here, we explore three functionally critical targeting short linear motif systems as well as their receptors in depth, using a combination of structural modeling, evolutionary sequence divergence and deep learning. Secretory signal peptides, endoplasmic reticulum (ER) retention motifs (KDEL motifs), and autophagy signals (motifs interacting with ATG8 family members) are ancient and essential components of cellular life. Although expected to be conserved amongst the kinetoplastids, we observe that all three systems show a varying degree of divergence from their better studied equivalents in animals, plants, or fungi. We not only describe their behaviour, but also build models that allow the prediction of localization and potential functions for several uncharacterized Leishmania proteins. The unusually Ala/Val-rich secretory signal peptides, endoplasmic reticulum resident proteins ending in Asp-Leu-COOH and atypical ATG8-like proteins are all unique molecular features of kinetoplastid parasites. Several of their critical protein-protein interactions could serve as targets of selective antimicrobial agents against Leishmaniasis due to their systematic divergence from the host.

Deletion of the lipid droplet protein kinase gene affects lipid droplets biogenesis, parasite infectivity, and resistance to trivalent antimony in Leishmania infantum

The Lipid Droplet Protein Kinase (LDK) facilitates lipid droplet (LD) biogenesis, organelles involved in various metabolic and signaling functions in trypanosomatids. As LDK's function has not been previously explored in Leishmania spp., we utilized CRISPR/Cas9 technology to generate LDK-knockout lines of Leishmania infantum to investigate its role in this parasite. Our findings demonstrate that LDK is not an essential gene for L. infantum, as its deletion did not impede parasite survival. Furthermore, removing LDK did not impact the growth of promastigote forms of L. infantum lacking LDK. However, a noticeable reduction in LDs occurred during the stationary phase of parasite growth following LDK deletion. In the presence of myriocin, a LD inducer, LDK-knockout parasites displayed reduced LD abundance during both logarithmic and stationary growth phases compared to control parasites. Moreover, an infection analysis involving THP-1 cells revealed that 72 h post-infection, LDK-knockout L. infantum lines exhibited fewer infected macrophages and intracellular amastigotes than control parasites. LDK-knockout L. infantum lines also displayed 1.7 to 1.8 -fold greater resistance to trivalent antimony than control parasites. There were no observed alterations in susceptibility to amphotericin B, miltefosine, or menadione in LDK-knockout L. infantum lines. Our results suggest that LDK plays a crucial role in the biogenesis and/or maintenance of LDs in L. infantum, as well as in parasite infectivity and resistance to trivalent antimony.

A sticky situation: When trypanosomatids attach to insect tissues

Transmission of trypanosomatids to their mammalian hosts requires a complex series of developmental transitions in their insect vectors, including stable attachment to an insect tissue. While there are many ultrastructural descriptions of attached cells, we know little about the signaling events and molecular mechanisms involved in this process. Each trypanosomatid species attaches to a specific tissue in the insect at a particular stage of its life cycle. Attachment is mediated by the flagellum, which is modified to accommodate a filament-rich plaque within an expanded region of the flagellar membrane. Attachment immediately precedes differentiation to the mammal-infectious stage and in some cases a direct mechanistic link has been demonstrated. In this review, we summarize the current state of knowledge of trypanosomatid attachment in insects, including structure, function, signaling, candidate molecules, and changes in gene expression. We also highlight remaining questions about this process and how the field is poised to address them through modern approaches.

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.

Wide reference databases for typing Trypanosoma cruzi based on amplicon sequencing of the minicircle hypervariable region

BACKGROUND

Trypanosoma cruzi, the etiological agent of Chagas Disease, exhibits remarkable genetic diversity and is classified into different Discrete Typing Units (DTUs). Strain typing techniques are crucial for studying T. cruzi, because their DTUs have significant biological differences from one another. However, there is currently no methodological strategy for the direct typing of biological materials that has sufficient sensitivity, specificity, and reproducibility. The high diversity and copy number of the minicircle hypervariable regions (mHVRs) makes it a viable target for typing.

METHODOLOGY/PRINCIPAL FINDINGS

Approximately 24 million reads obtained by amplicon sequencing of the mHVR were analyzed for 62 strains belonging to the six main T. cruzi DTUs. To build reference databases of mHVR diversity for each DTU and to evaluate this target as a typing tool. Strains of the same DTU shared more mHVR clusters than strains of different DTUs, and clustered together. Different identity thresholds were used to build the reference sets of the mHVR sequences (85% and 95%, respectively). The 95% set had a higher specificity and was more suited for detecting co-infections, whereas the 85% set was excellent for identifying the primary DTU of a sample. The workflow's capacity for typing samples obtained from cultures, a set of whole-genome data, under various simulated PCR settings, in the presence of co-infecting lineages and for blood samples was also assessed.

CONCLUSIONS/SIGNIFICANCE

We present reference databases of mHVR sequences and an optimized typing workflow for T. cruzi including a simple online tool for deep amplicon sequencing analysis (https://ntomasini.github.io/cruzityping/). The results show that the workflow displays an equivalent resolution to that of the other typing methods. Owing to its specificity, sensitivity, relatively low cost, and simplicity, the proposed workflow could be an alternative for screening different types of samples.

Diversity of RNA viruses in the cosmopolitan monoxenous trypanosomatid Leptomonas pyrrhocoris

BACKGROUND

Trypanosomatids are parasitic flagellates well known because of some representatives infecting humans, domestic animals, and cultural plants. Many trypanosomatid species bear RNA viruses, which, in the case of human pathogens Leishmania spp., influence the course of the disease. One of the close relatives of leishmaniae, Leptomonas pyrrhocoris, has been previously shown to harbor viruses of the groups not documented in other trypanosomatids. At the same time, this species has a worldwide distribution and high prevalence in the natural populations of its cosmopolitan firebug host. It therefore represents an attractive model to study the diversity of RNA viruses.

RESULTS

We surveyed 106 axenic cultures of L. pyrrhocoris and found that 64 (60%) of these displayed 2-12 double-stranded RNA fragments. The analysis of next-generation sequencing data revealed four viral groups with seven species, of which up to five were simultaneously detected in a single trypanosomatid isolate. Only two of these species, a tombus-like virus and an Ostravirus, were earlier documented in L. pyrrhocoris. In addition, there were four new species of Leishbuviridae, the family encompassing trypanosomatid-specific viruses, and a new species of Qinviridae, the family previously known only from metatranscriptomes of invertebrates. Currently, this is the only qinvirus with an unambiguously determined host. Our phylogenetic inferences suggest reassortment in the tombus-like virus owing to the interaction of different trypanosomatid strains. Two of the new Leishbuviridae members branch early on the phylogenetic tree of this family and display intermediate stages of genomic segment reduction between insect Phenuiviridae and crown Leishbuviridae.

CONCLUSIONS

The unprecedented wide range of viruses in one protist species and the simultaneous presence of up to five viral species in a single Leptomonas pyrrhocoris isolate indicate the uniqueness of this flagellate. This is likely determined by the peculiarity of its firebug host, a highly abundant cosmopolitan species with several habits ensuring wide distribution and profuseness of L. pyrrhocoris, as well as its exposure to a wider spectrum of viruses compared to other trypanosomatids combined with a limited ability to transmit these viruses to its relatives. Thus, L. pyrrhocoris represents a suitable model to study the adoption of new viruses and their relationships with a protist host.

Shining the spotlight on the neglected: new high-quality genome assemblies as a gateway to understanding the evolution of Trypanosomatidae

BACKGROUND

Protists of the family Trypanosomatidae (phylum Euglenozoa) have gained notoriety as parasites affecting humans, domestic animals, and agricultural plants. However, the true extent of the group's diversity spreads far beyond the medically and veterinary relevant species. We address several knowledge gaps in trypanosomatid research by undertaking sequencing, assembly, and analysis of genomes from previously overlooked representatives of this protistan group.

RESULTS

We assembled genomes for twenty-one trypanosomatid species, with a primary focus on insect parasites and Trypanosoma spp. parasitizing non-human hosts. The assemblies exhibit sizes consistent with previously sequenced trypanosomatid genomes, ranging from approximately 18 Mb for Obscuromonas modryi to 35 Mb for Crithidia brevicula and Zelonia costaricensis. Despite being the smallest, the genome of O. modryi has the highest content of repetitive elements, contributing nearly half of its total size. Conversely, the highest proportion of unique DNA is found in the genomes of Wallacemonas spp., with repeats accounting for less than 8% of the assembly length. The majority of examined species exhibit varying degrees of aneuploidy, with trisomy being the most frequently observed condition after disomy.

CONCLUSIONS

The genome of Obscuromonas modryi represents a very unusual, if not unique, example of evolution driven by two antidromous forces: i) increasing dependence on the host leading to genomic shrinkage and ii) expansion of repeats causing genome enlargement. The observed variation in somy within and between trypanosomatid genera suggests that these flagellates are largely predisposed to aneuploidy and, apparently, exploit it to gain a fitness advantage. High heterogeneity in the genome size, repeat content, and variation in chromosome copy numbers in the newly-sequenced species highlight the remarkable genome plasticity exhibited by trypanosomatid flagellates. These new genome assemblies are a robust foundation for future research on the genetic basis of life cycle changes and adaptation to different hosts in the family Trypanosomatidae.

Kinetoplast Genome of Leishmania spp. Is under Strong Purifying Selection

Instability is an intriguing characteristic of many protist genomes, and trypanosomatids are not an exception in this respect. Some regions of trypanosomatid genomes evolve fast. For instance, the trypanosomatid mitochondrial (kinetoplast) genome consists of fairly conserved maxicircle and minicircle molecules that can, nevertheless, possess high nucleotide substitution rates between closely related strains. Recent experiments have demonstrated that rapid laboratory evolution can result in the non-functionality of multiple genes of kinetoplast genomes due to the accumulation of mutations or loss of critical genomic components. An example of a loss of critical components is the reported loss of entire minicircle classes in Leishmania tarentolae during laboratory cultivation, which results in an inability to generate some correctly encoded genes. In the current work, we estimated the evolutionary rates of mitochondrial and nuclear genome regions of multiple natural Leishmania spp. We analyzed synonymous and non-synonymous substitutions and, rather unexpectedly, found that the coding regions of kinetoplast maxicircles are among the most variable regions of both genomes. In addition, we demonstrate that synonymous substitutions greatly predominate among maxicircle coding regions and that most maxicircle genes show signs of purifying selection. These results imply that maxicircles in natural Leishmania populations remain functional despite their high mutation rate.

KREH1 RNA helicase activity promotes utilization of initiator gRNAs across multiple mRNAs in trypanosome RNA editing

Mitochondrial U-indel RNA editing in kinetoplastid protozoa is directed by trans-acting gRNAs and mediated by a holoenzyme with associated factors. Here, we examine the function of the holoenzyme-associated KREH1 RNA helicase in U-indel editing. We show that KREH1 knockout (KO) impairs editing of a small subset of mRNAs. Overexpression of helicase-dead mutants results in expanded impairment of editing across multiple transcripts, suggesting the existence of enzymes that can compensate for KREH1 in KO cells. In depth analysis of editing defects using quantitative RT-PCR and high-throughput sequencing reveals compromised editing initiation and progression in both KREH1-KO and mutant-expressing cells. In addition, these cells exhibit a distinct defect in the earliest stages of editing in which the initiator gRNA is bypassed, and a small number of editing events takes place just outside this region. Wild type KREH1 and a helicase-dead KREH1 mutant interact similarly with RNA and holoenzyme, and overexpression of both similarly disorders holoenzyme homeostasis. Thus, our data support a model in which KREH1 RNA helicase activity facilitates remodeling of initiator gRNA-mRNA duplexes to permit accurate utilization of initiating gRNAs on multiple transcripts.

Contrasting outcomes of genome reduction in mikrocytids and microsporidians

BACKGROUND

Intracellular symbionts often undergo genome reduction, losing both coding and non-coding DNA in a process that ultimately produces small, gene-dense genomes with few genes. Among eukaryotes, an extreme example is found in microsporidians, which are anaerobic, obligate intracellular parasites related to fungi that have the smallest nuclear genomes known (except for the relic nucleomorphs of some secondary plastids). Mikrocytids are superficially similar to microsporidians: they are also small, reduced, obligate parasites; however, as they belong to a very different branch of the tree of eukaryotes, the rhizarians, such similarities must have evolved in parallel. Since little genomic data are available from mikrocytids, we assembled a draft genome of the type species, Mikrocytos mackini, and compared the genomic architecture and content of microsporidians and mikrocytids to identify common characteristics of reduction and possible convergent evolution.

RESULTS

At the coarsest level, the genome of M. mackini does not exhibit signs of extreme genome reduction; at 49.7 Mbp with 14,372 genes, the assembly is much larger and gene-rich than those of microsporidians. However, much of the genomic sequence and most (8075) of the protein-coding genes code for transposons, and may not contribute much of functional relevance to the parasite. Indeed, the energy and carbon metabolism of M. mackini share several similarities with those of microsporidians. Overall, the predicted proteome involved in cellular functions is quite reduced and gene sequences are extremely divergent. Microsporidians and mikrocytids also share highly reduced spliceosomes that have retained a strikingly similar subset of proteins despite having reduced independently. In contrast, the spliceosomal introns in mikrocytids are very different from those of microsporidians in that they are numerous, conserved in sequence, and constrained to an exceptionally narrow size range (all 16 or 17 nucleotides long) at the shortest extreme of known intron lengths.

CONCLUSIONS

Nuclear genome reduction has taken place many times and has proceeded along different routes in different lineages. Mikrocytids show a mix of similarities and differences with other extreme cases, including uncoupling the actual size of a genome with its functional reduction.

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.

Development of a TaqMan qPCR assay for trypanosomatid multi-species detection and quantification in insects

BACKGROUND

Trypanosomatid parasites are widely distributed in nature and can have a monoxenous or dixenous life-cycle. These parasites thrive in a wide number of insect orders, some of which have an important economic and environmental value, such as bees. The objective of this study was to develop a robust and sensitive real-time quantitative PCR (qPCR) assay for detecting trypanosomatid parasites in any type of parasitized insect sample.

METHODS

A TaqMan qPCR assay based on a trypanosomatid-conserved region of the α-tubulin gene was standardized and evaluated. The limits of detection, sensitivity and versatility of the α-tubulin TaqMan assay were tested and validated using field samples of honeybee workers, wild bees, bumblebees and grasshoppers, as well as in the human infective trypanosomatid Leishmania major.

RESULTS

The assay showed a detection limit of 1 parasite equivalent/µl and successfully detected trypanosomatids in 10 different hosts belonging to the insect orders Hymenoptera and Orthoptera. The methodology was also tested using honeybee samples from four apiaries (n = 224 worker honeybees) located in the Alpujarra region (Granada, Spain). Trypanosomatids were detected in 2.7% of the honeybees, with an intra-colony prevalence of 0% to 13%. Parasite loads in the four different classes of insects ranged from 40.6 up to 1.1 × 10⁸ cell equivalents per host.

CONCLUSIONS

These results show that the α-tubulin TaqMan qPCR assay described here is a versatile diagnostic tool for the accurate detection and quantification of trypanosomatids in a wide range of environmental settings.

Short tRNA anticodon stem and mutant eRF1 allow stop codon reassignment

Cognate tRNAs deliver specific amino acids to translating ribosomes according to the standard genetic code, and three codons with no cognate tRNAs serve as stop codons. Some protists have reassigned all stop codons as sense codons, neglecting this fundamental principle^1-4. Here we analyse the in-frame stop codons in 7,259 predicted protein-coding genes of a previously undescribed trypanosomatid, Blastocrithidia nonstop. We reveal that in this species in-frame stop codons are underrepresented in genes expressed at high levels and that UAA serves as the only termination codon. Whereas new tRNAs^Glu fully cognate to UAG and UAA evolved to reassign these stop codons, the UGA reassignment followed a different path through shortening the anticodon stem of tRNA^Trp_CCA from five to four base pairs (bp). The canonical 5-bp tRNA^Trp recognizes UGG as dictated by the genetic code, whereas its shortened 4-bp variant incorporates tryptophan also into in-frame UGA. Mimicking this evolutionary twist by engineering both variants from B. nonstop, Trypanosoma brucei and Saccharomyces cerevisiae and expressing them in the last two species, we recorded a significantly higher readthrough for all 4-bp variants. Furthermore, a gene encoding B. nonstop release factor 1 acquired a mutation that specifically restricts UGA recognition, robustly potentiating the UGA reassignment. Virtually the same strategy has been adopted by the ciliate Condylostoma magnum. Hence, we describe a previously unknown, universal mechanism that has been exploited in unrelated eukaryotes with reassigned stop codons.

The Roles of Mosquitoes in the Circulation of Monoxenous Trypanosomatids in Temperate Climates

Monoxenous (insect-restricted) trypanosomatids are highly diverse and abundant in nature. While many papers focus on the taxonomy and distribution of these parasites, studies on their biology are still scarce. In particular, this concerns trypanosomatids inhabiting the ubiquitous mosquitoes. To shed light on the circulation of monoxenous trypanosomatids with the participation of mosquitoes, we performed a multifaceted study combining the examination of naturally- and experimentally-infected insects using light and electron microscopy and molecular identification of parasites. Our examination of overwintering mosquitoes (genera Culex and Culiseta) revealed that their guts contained living trypanosomatids, which can be spread during the next season. Experimental infections with Crithidia spp. demonstrated that imagines represent permissive hosts, while larvae are resistant to these parasites. We argue that for the parasites with wide specificity, mosquitoes act as facultative hosts. Other trypanosomatids may have specific adaptations for vertical transmission in these insects at the expense of their potential to infect a wider range of hosts and, consequently, abundance in nature.

Kinetoplastid Species Maintained by a Small Mammal Community in the Pantanal Biome

Kinetoplastids include species economically important in agriculture, livestock, and human health. We evaluated the richness of kinetoplastids that infect small mammals in patches of unflooded forests in the Pantanal biome, an area where we hypothesize that its diversity is higher than currently recognized. Hemocultures (HC) and Next Generation Sequencing (NGS) targeting the 18S rDNA gene were employed for the detection of kinetoplastids. We grouped the positive samples into pools for each small mammal species (Monodelphis domestica, Thylamys macrurus, Oecomys mamorae, Thrichomys fosteri, Clyomys laticeps, and Holochilus chacarius). Eight parasite species were identified: Leishmania amazonensis, L. infantum; Trypanosoma cascavelli (HC + NGS), T. cruzi, T. lainsoni, T. rangeli (HC + NGS), Trypanosoma sp. DID, and Neobodo sp. The use of a tool as sensitive as NGS has increased our awareness of the diversity of kinetoplastids, as well as their host range, with emphasis on the species O. mamorae (seven kinetoplastid species, excepting T. cascavelli in a pool of nine individuals) and T. macrurus (four kinetoplastid species in a single individual). Furthermore, L. infantum and L. amazonensis infections were described in small mammals from this region for the first time. These findings make it mandatory to revisit the kinetoplastids/host associations proposed so far.

Uncovering Pseudogenes and Intergenic Protein-coding Sequences in TriTryps' Genomes

Trypanosomatids belong to a remarkable group of unicellular, parasitic organisms of the order Kinetoplastida, an early diverging branch of the phylogenetic tree of eukaryotes, exhibiting intriguing biological characteristics affecting gene expression (intronless polycistronic transcription, trans-splicing, and RNA editing), metabolism, surface molecules, and organelles (compartmentalization of glycolysis, variation of the surface molecules, and unique mitochondrial DNA), cell biology and life cycle (phagocytic vacuoles evasion and intricate patterns of cell morphogenesis). With numerous genomic-scale data of several trypanosomatids becoming available since 2005 (genomes, transcriptomes, and proteomes), the scientific community can further investigate the mechanisms underlying these unusual features and address other unexplored phenomena possibly revealing biological aspects of the early evolution of eukaryotes. One fundamental aspect comprises the processes and mechanisms involved in the acquisition and loss of genes throughout the evolutionary history of these primitive microorganisms. Here, we present a comprehensive in silico analysis of pseudogenes in three major representatives of this group: Leishmania major, Trypanosoma brucei, and Trypanosoma cruzi. Pseudogenes, DNA segments originating from altered genes that lost their original function, are genomic relics that can offer an essential record of the evolutionary history of functional genes, as well as clues about the dynamics and evolution of hosting genomes. Scanning these genomes with functional proteins as proxies to reveal intergenic regions with protein-coding features, relying on a customized threshold to distinguish statistically and biologically significant sequence similarities, and reassembling remnant sequences from their debris, we found thousands of pseudogenes and hundreds of open reading frames, with particular characteristics in each trypanosomatid: mutation profile, number, content, density, codon bias, average size, single- or multi-copy gene origin, number and type of mutations, putative primitive function, and transcriptional activity. These features suggest a common process of pseudogene formation, different patterns of pseudogene evolution and extant biological functions, and/or distinct genome organization undertaken by those parasites during evolution, as well as different evolutionary and/or selective pressures acting on distinct lineages.

Distinct mRNA and protein interactomes highlight functional differentiation of major eIF4F-like complexes from Trypanosoma brucei

Gene expression in pathogenic protozoans of the family Trypanosomatidae has several novel features, including multiple eIF4F-like complexes involved in protein synthesis. The eukaryotic eIF4F complex, formed mainly by eIF4E and eIF4G subunits, is responsible for the canonical selection of mRNAs required for the initiation of mRNA translation. The best-known complexes implicated in translation in trypanosomatids are based on two related pairs of eIF4E and eIF4G subunits (EIF4E3/EIF4G4 and EIF4E4/EIF4G3), whose functional distinctions remain to be fully described. Here, to define interactomes associated with both complexes in Trypanosoma brucei procyclic forms, we performed parallel immunoprecipitation experiments followed by identification of proteins co-precipitated with the four tagged eIF4E and eIF4G subunits. A number of different protein partners, including RNA binding proteins and helicases, specifically co-precipitate with each complex. Highlights with the EIF4E4/EIF4G3 pair include RBP23, PABP1, EIF4AI and the CRK1 kinase. Co-precipitated partners with the EIF4E3/EIF4G4 pair are more diverse and include DRBD2, PABP2 and different zinc-finger proteins and RNA helicases. EIF4E3/EIF4G4 are essential for viability and to better define their role, we further investigated their phenotypes after knockdown. Depletion of either EIF4E3/EIF4G4 mRNAs lead to aberrant morphology with a more direct impact on events associated with cytokinesis. We also sought to identify those mRNAs differentially associated with each complex through CLIP-seq with the two eIF4E subunits. Predominant among EIF4E4-bound transcripts are those encoding ribosomal proteins, absent from those found with EIF4E3, which are generally more diverse. RNAi mediated depletion of EIF4E4, which does not affect proliferation, does not lead to changes in mRNAs or proteins associated with EIF4E3, confirming a lack of redundancy and distinct roles for the two complexes.

Sending the message: specialized RNA export mechanisms in trypanosomes

Export of RNA from the nucleus is essential for all eukaryotic cells and has emerged as a major step in the control of gene expression. mRNA molecules are required to complete a complex series of processing events and pass a quality control system to protect the cytoplasm from the translation of aberrant proteins. Many of these events are highly conserved across eukaryotes, reflecting their ancient origin, but significant deviation from a canonical pathway as described from animals and fungi has emerged in the trypanosomatids. With significant implications for the mechanisms that control gene expression and hence differentiation, responses to altered environments and fitness as a parasite, these deviations may also reveal additional, previously unsuspected, mRNA export pathways.

Bee Trypanosomatids: First Steps in the Analysis of the Genetic Variation and Population Structure of Lotmaria passim, Crithidia bombi and Crithidia mellificae

Trypanosomatids are among the most prevalent parasites in bees but, despite the fact that their impact on the colonies can be quite important and that their infectivity may potentially depend on their genotypes, little is known about the population diversity of these pathogens. Here we cloned and sequenced three non-repetitive single copy loci (DNA topoisomerase II, glyceraldehyde-3-phosphate dehydrogenase and RNA polymerase II large subunit, RPB1) to produce new genetic data from Crithidia bombi, C. mellificae and Lotmaria passim isolated from honeybees and bumblebees. These were analysed by applying population genetic tools in order to quantify and compare their variability within and between species, and to obtain information on their demography and population structure. The general pattern for the three species was that (1) they were subject to the action of purifying selection on nonsynonymous variants, (2) the levels of within species diversity were similar irrespective of the host, (3) there was evidence of recombination among haplotypes and (4) they showed no haplotype structuring according to the host. C. bombi exhibited the lowest levels of synonymous variation (π_S= 0.06 ± 0.04 %) - and a mutation frequency distribution compatible with a population expansion after a bottleneck - that contrasted with the extensive polymorphism displayed by C. mellificae (π_S= 2.24 ± 1.00 %), which likely has a more ancient origin. L. passim showed intermediate values (π_S= 0.40 ± 0.28 %) and an excess of variants a low frequencies probably linked to the spread of this species to new geographical areas.

Differences in Charge Distribution in Leishmania tarentolae Leishmanolysin Result in a Reduced Enzymatic Activity

Leishmania tarentolae is a non-pathogenic trypanosomatid isolated from lizards widely used for heterologous protein expression and extensively studied to understand the pathogenic mechanisms of leishmaniasis. The repertoire of leishmanolysin genes was reported to be expanded in L. tarentolae genome, but no proteolytic activity was detected. Here, we analyzed L. tarentolae leishmanolysin proteins from the genome to the structural levels and evaluated the enzymatic activity of the wild-type and overexpressing mutants of leishmanolysin. A total of 61 leishmanolysin sequences were retrieved from the L. tarentolae genome. Five of them were selected for phylogenetic analysis, and for three of them, we built 3D models based on the crystallographic structure of L. major ortholog. Molecular dynamics simulations of these models disclosed a less negative electrostatic potential compared to the template. Subsequently, L. major LmjF.10.0460 and L. tarentolae LtaP10.0650 leishmanolysins were cloned in a pLEXSY expression system into L. tarentolae. Proteins from the wild-type and the overexpressing parasites were submitted to enzymatic analysis. Our results revealed that L. tarentolae leishmanolysins harbor a weak enzymatic activity about three times less abundant than L. major leishmanolysin. Our findings strongly suggest that the less negative electrostatic potential of L. tarentolae leishmanolysin can be the reason for the reduced proteolytic activity detected in this parasite.

Trypanosoma cruzi strain and starvation-driven mitochondrial RNA editing and transcriptome variability

Trypanosoma cruzi is a unicellular protistan parasitic species that is comprised of strains and isolates exhibiting high levels of genetic and metabolic variability. In the insect vector, it is known to be highly responsive to starvation, a signal for progression to a life stage in which it can infect mammalian cells. Most mRNAs encoded in its mitochondrion require the targeted insertion and deletion of uridines to become translatable transcripts. This study defined differences in uridine-insertion/deletion RNA editing among three strains and established the mechanism whereby abundances of edited (and, thus, translatable) mitochondrial gene products increase during starvation. Our approach utilized our custom T-Aligner toolkit to describe transcriptome-wide editing events and reconstruct editing products from high-throughput sequencing data. We found that the relative abundance of mitochondrial transcripts and the proportion of mRNAs that are edited varies greatly between analyzed strains, a characteristic that could potentially impact metabolic capacity. Starvation typically led to an increase in overall editing activity rather than affecting a specific step in the process. We also determined that transcripts CR3, CR4, and ND3 produce multiple open reading frames that, if translated, would generate different proteins. Finally, we quantitated the inherent flexibility of editing in T. cruzi and found it to be higher relative to that in a related trypanosomatid lineage. Over time, new editing domains or patterns could prove advantageous to the organism and become more widespread within individual transcriptomes or among strains.

Trypanosomatid Richness Among Rats, Opossums, and Dogs in the Caatinga Biome, Northeast Brazil, a Former Endemic Area of Chagas Disease

Parasites are important components of the immense n-dimensional trophic network that connects all living beings because they, among others, forge biodiversity and deeply influence ecological evolution and host behavior. In this sense, the influence of Trypanosomatidae remains unknown. The aim of this study was to determine trypanosomatid infection and richness in rats, opossums, and dogs in the semiarid Caatinga biome. We submitted DNA samples from trypanosomatids obtained through axenic cultures of the blood of these mammals to mini exon multiplex-PCR, Sanger, and next-generation sequencing targeting the 18S rDNA gene. Phylogenetic analyses were performed to identify genetic diversity in the Trypanosomatidae family. Shannon, Simpson, equability, and beta-diversity indices were calculated per location and per mammalian host. Dogs were surveyed for trypanosomatid infection through hemocultures and serological assays. The examined mammal species of this area of the Caatinga biome exhibited an enormous trypanosomatid species/genotypes richness. Ten denoised Operational Taxonomic Units (ZOTUs), including three species (Trypanosoma cruzi, Trypanosoma rangeli and Crithidia mellificae) and one Trypanosoma sp. five genotypes/lineages (T. cruzi DTU TcI, TcII, and TcIV; T. rangeli A and B) and four DTU TcI haplotypes (ZOTU1, ZOTU2, ZOTU5, and ZOTU10 merged), as well as 13 Amplicon Sequence Variants (ASVs), including five species (T. cruzi, T. rangeli, C. mellificae, Trypanosoma dionisii, and Trypanosoma lainsoni), five genotypes/lineages (same as the ZOTUs) and six DTU TcI haplotypes (ASV, ASV1, ASV2, ASV3, ASV5 and ASV13), were identified in single and mixed infections. We observed that trypanosomatids present a broad host spectrum given that species related to a single host are found in other mammals from different taxa. Concomitant infections between trypanosomatids and new host-parasite relationships have been reported, and this immense diversity in mammals raised questions, such as how this can influence the course of the infection in these animals and its transmissibility. Dogs demonstrated a high infection rate by T. cruzi as observed by positive serological results (92% in 2005 and 76% in 2007). The absence of positive parasitological tests confirmed their poor infectivity potential but their importance as sentinel hosts of T. cruzi transmission.

Leishmania guyanensis M4147 as a new LRV1-bearing model parasite: Phosphatidate phosphatase 2-like protein controls cell cycle progression and intracellular lipid content

Leishmaniasis is a parasitic vector-borne disease caused by the protistan flagellates of the genus Leishmania. Leishmania (Viannia) guyanensis is one of the most common causative agents of the American tegumentary leishmaniasis. It has previously been shown that L. guyanensis strains that carry the endosymbiotic Leishmania RNA virus 1 (LRV1) cause more severe form of the disease in a mouse model than those that do not. The presence of the virus was implicated into the parasite’s replication and spreading. In this respect, studying the molecular mechanisms of cellular control of viral infection is of great medical importance. Here, we report ~30.5 Mb high-quality genome assembly of the LRV1-positive L. guyanensis M4147. This strain was turned into a model by establishing the CRISPR-Cas9 system and ablating the gene encoding phosphatidate phosphatase 2-like (PAP2L) protein. The orthologue of this gene is conspicuously absent from the genome of an unusual member of the family Trypanosomatidae, Vickermania ingenoplastis, a species with mostly bi-flagellated cells. Our analysis of the PAP2L-null L. guyanensis showed an increase in the number of cells strikingly resembling the bi-flagellated V. ingenoplastis, likely as a result of the disruption of the cell cycle, significant accumulation of phosphatidic acid, and increased virulence compared to the wild type cells.

Worldwide, over one million people are getting infected by the parasitic flagellates of the genus Leishmania annually leading to ~30,000 deaths. Notably, there is still no approved vaccine against human leishmaniases. A range of methods of forward and reverse genetics has recently been developed for several model Leishmania species. Unfortunately, these methods are often not transferrable to non-model species, which may be of even greater medical importance. Leishmania guyanensis is one of such cases. It frequently carries a symbiotic RNA virus that contributes to the development of a more aggressive form of leishmaniasis in an experimental murine model. In order to establish and optimize the system for genetic manipulations in L. guyanensis, we sequenced and annotated its genome. Next, we applied the CRISPR-Cas9 technology to target a gene of interest. This approach was validated by ablating a gene encoding a protein involved in lipid metabolism. In this work, we document that deletion of this gene leads to the disturbance of cell cycle and affects the ratio of critical intracellular lipids. We believe that our study will facilitate research into more effective treatment of leishmaniases.

G Protein-Coupled Receptors as Potential Intercellular Communication Mediators in Trypanosomatidae

Detection and transduction of environmental signals, constitute a prerequisite for successful parasite invasion; i.e., Leishmania transmission, survival, pathogenesis and disease manifestation and dissemination, with diverse molecules functioning as inter-cellular signaling ligands. Receptors [i.e., G protein-coupled receptors (GPCRs)] and their associated transduction mechanisms, well conserved through evolution, specialize in this function. However, canonical GPCR-related signal transduction systems have not been described in Leishmania, although orthologs, with reduced domains and function, have been identified in Trypanosomatidae. These inter-cellular communication means seem to be essential for multicellular and unicellular organism’s survival. GPCRs are flexible in their molecular architecture and may interact with the so-called receptor activity-modifying proteins (RAMPs), which modulate their function, changing GPCRs pharmacology, acting as chaperones and regulating signaling and/or trafficking in a receptor-dependent manner. In the skin, vasoactive- and neuro- peptides released in response to the noxious stimuli represented by the insect bite may trigger parasite physiological responses, for example, chemotaxis. For instance, in Leishmania (V.) braziliensis, sensory [Substance P, SP, chemoattractant] and autonomic [Vasoactive Intestinal Peptide, VIP, and Neuropeptide Y, NPY, chemorepellent] neuropeptides at physiological levels stimulate in vitro effects on parasite taxis. VIP and NPY chemotactic effects are impaired by their corresponding receptor antagonists, suggesting that the stimulated responses might be mediated by putative GPCRs (with essential conserved receptor domains); the effect of SP is blocked by [(D-Pro 2, D-Trp7,9]-Substance P (10^-6 M)] suggesting that it might be mediated by neurokinin-1 transmembrane receptors. Additionally, vasoactive molecules like Calcitonin Gene-Related Peptide [CGRP] and Adrenomedullin [AM], exert a chemorepellent effect and increase the expression of a 24 kDa band recognized in western blot analysis by (human-)-RAMP-2 antibodies. In-silico search oriented towards GPCRs-like receptors and signaling cascades detected a RAMP-2-aligned sequence corresponding to Leishmania folylpolyglutamate synthase and a RAMP-3 aligned protein, a hypothetical Leishmania protein with yet unknown function, suggesting that in Leishmania, CGRP and AM activities may be modulated by RAMP- (-2) and (-3) homologs. The possible presence of proteins and molecules potentially involved in GPCRs cascades, i.e., RAMPs, signpost conservation of ancient signaling systems associated with responses, fundamental for cell survival, (i.e., taxis and migration) and may constitute an open field for description of pharmacophores against Leishmania parasites.

The genomic basis of host and vector specificity in non-pathogenic trypanosomatids

Trypanosoma theileri, a non-pathogenic parasite of bovines, has a predicted surface protein architecture that likely aids survival in its mammalian host. Their surface proteins are encoded by genes which account for ∼10% of their genome. A non-pathogenic parasite of sheep, Trypanosoma melophagium, is transmitted by the sheep ked and is closely related to T. theileri. To explore host and vector specificity between these species, we sequenced the T. melophagium genome and transcriptome and an annotated draft genome was assembled. T. melophagium was compared to 43 kinetoplastid genomes, including T. theileri. T. melophagium and T. theileri have an AT biased genome, the greatest bias of publicly available trypanosomatids. This trend may result from selection acting to decrease the genomic nucleotide cost. The T. melophagium genome is 6.3Mb smaller than T. theileri and large families of proteins, characteristic of the predicted surface of T. theileri, were found to be absent or greatly reduced in T. melophagium. Instead, T. melophagium has modestly expanded protein families associated with the avoidance of complement-mediated lysis. We propose that the contrasting genomic features of these species is linked to their mode of transmission from their insect vector to their mammalian host. This article has an associated First Person interview with the first author of the paper.

Trypanosoma Species in Small Nonflying Mammals in an Area With a Single Previous Chagas Disease Case

Trypanosomatids are hemoflagellate parasites that even though they have been increasingly studied, many aspects of their biology and taxonomy remain unknown. The aim of this study was to investigate the Trypanosoma sp. transmission cycle in nonflying small mammals in an area where a case of acute Chagas disease occurred in Mangaratiba municipality, Rio de Janeiro state. Three expeditions were conducted in the area: the first in 2012, soon after the human case, and two others in 2015. Sylvatic mammals were captured and submitted to blood collection for trypanosomatid parasitological and serological exams. Dogs from the surrounding areas where the sylvatic mammals were captured were also tested for T. cruzi infection. DNA samples were extracted from blood clots and positive hemocultures, submitted to polymerase chain reaction targeting SSU rDNA and gGAPDH genes, sequenced and phylogenetic analysed. Twenty-one wild mammals were captured in 2012, mainly rodents, and 17 mammals, mainly marsupials, were captured in the two expeditions conducted in 2015. Only four rodents demonstrated borderline serological T. cruzi test (IFAT), two in 2012 and two in 2015. Trypanosoma janseni was the main Trypanosoma species identified, and isolates were obtained solely from Didelphis aurita. In addition to biological differences, molecular differences are suggestive of genetic diversity in this flagellate species. Trypanosoma sp. DID was identified in blood clots from D. aurita in single and mixed infections with T. janseni. Concerning dogs, 12 presented mostly borderline serological titers for T. cruzi and no positive hemoculture. In blood clots from 11 dogs, T. cruzi DNA was detected and characterized as TcI (n = 9) or TcII (n = 2). Infections by Trypanosoma rangeli lineage E (n = 2) and, for the first time, Trypanosoma caninum, Trypanosoma dionisii, and Crithidia mellificae (n = 1 each) were also detected in dogs. We concluded that despite the low mammalian species richness and degraded environment, a high Trypanosoma species richness species was being transmitted with the predominance of T. janseni and not T. cruzi, as would be expected in a locality of an acute case of Chagas disease.

Mutation Characteristics and Phylogenetic Analysis of Five Leishmania Clinical Isolates

Simple Summary

Leishmaniasis, a neglected tropical disease, is caused by infection with the Leishmania species, threatening millions of people in approximately 100 endemic countries. The emergence of antimony-resistant Leishmania strains have brought difficulties to the treatment and elimination of leishmaniasis. This study performed genome-wide resequencing and phylogenetic analysis of five isolates from the Leishmania donovani complex, focusing on finding mutations related to antimony resistance and virulence of the newly isolated Leishmania strain L_HCZ in 2016. By combining whole-genome sequencing and whole-genome phylogenetic analysis, Leishmania isolates L_801, L_9044 and L_Liu were identified as Leishmania donovani, and L_HCZ as Leishmania infantum. By discovering genome-wide single-nucleotide polymorphisms and structural variations, we identified mutations of drug resistance-related genes in the antimony-resistant Leishmania isolate L_HCZ. The new Leishmania isolate L_HCZ has strong virulence and strong drug resistance, which should be taken seriously by the relevant health departments and scientific researchers.

Abstract

Leishmaniasis is a neglected tropical disease threatening millions of people worldwide. The emergence of antimony-resistant Leishmania strains have brought difficulties to the treatment and elimination of leishmaniasis. This study performed genome sequencing, phylogenetic analysis and mutation analysis of five Leishmania clinical isolates, especially the Leishmania strain L_HCZ isolated in 2016, which shows strong virulence and antimony resistance. By phylogenetic analysis, four isolates (L_DD8, L_801, L_Liu and L_9044) were identified as Leishmania donovani, the isolate L_HCZ was identified as Leishmania infantum and the isolate L_DD8 as a standard strain of L.donovani. Genome-wide mutation analysis was applied to identify mutations related to the drug resistance and virulence of the newly isolated L_HCZ. Compared with the other four Leishmania isolates, L_HCZ had the most mutations in genes associated with antimony resistance, including the ABC transporter, ascorbate-dependent peroxidase, gamma–glutamylcysteine synthetase, glucose-6-phosphate 1-dehydrogenase, ATP-binding cassette protein subfamily A and multi-drug resistance protein-like genes. Among the genes associated with virulence, L_HCZ had the most mutations in cysteine peptidase A, cysteine peptidase B, cysteine peptidase C, heat-shock protein 70, gp63, acid phosphatase, kinesin k39, kinesin, phosphoglycan beta 1, amastin-like surface protein and amastin-like proteins. The mutations in L_HCZ might possibly contribute to its antimony resistance and strong virulence in clinical patients. Whole-genome resequencing has exhibited broad application prospects and may be put into clinical use in the future for parasite identifying and epidemiological investigations.

Leishmania and its relationships with bacteria

Leishmaniasis is a zoonotic and neglected disease, which represents an important public health problem worldwide. Different species of Leishmania are associated with different manifestations, and a practical problem that can worsen the condition of hosts infected with Leishmania is the secondary infection caused by bacteria. This review aims to examine the importance and prevalence of bacteria co-infection during leishmaniasis and the nature of this ecological relationship. In the cases discussed in this review, the facilitation phenomenon, defined as any interaction where the action of one organism has a beneficial effect on an organism of another species, was considered in the Leishmania-bacteria interaction, as well as the effects on one another and their consequences for the host.

Comparative Analysis of Three Trypanosomatid Catalases of Different Origin

Most trypanosomatid flagellates do not have catalase. In the evolution of this group, the gene encoding catalase has been independently acquired at least three times from three different bacterial groups. Here, we demonstrate that the catalase of Vickermania was obtained by horizontal gene transfer from Gammaproteobacteria, extending the list of known bacterial sources of this gene. Comparative biochemical analyses revealed that the enzymes of V. ingenoplastis, Leptomonas pyrrhocoris, and Blastocrithidia sp., representing the three independent catalase-bearing trypanosomatid lineages, have similar properties, except for the unique cyanide resistance in the catalase of the latter species.

A new species of mammalian trypanosome, Trypanosoma (Megatrypanum) bubalisi sp. nov., found in the freshwater leech Hirudinaria manillensis

Leeches have long been considered potential vectors for the aquatic lineage of trypanosomes, while bloodsucking insects are generally considered as the vectors for the terrestrial lineage of trypanosomes. The freshwater leech, Hirudinaria manillensis, is a widely distributed species in southern China and could potentially act as the vector for trypanosomes. Prior to this study, no trypanosomes had been reported from this leech. However, in this study, leeches were collected from three different places in Guangdong province, China, and a large number of flagellates were isolated and successfully cultured in vitro. Based on morphology, these flagellates looked like a typical trypanosome species. Analysis was carried out on the molecular sequences of the 18S rRNA gene and the glycosomal glyceraldehyde-3-phosphate dehydrogenase (gGAPDH) gene. To our surprise, these flagellates were identified as likely to be a mammalian trypanosome belonging to the clade containing Trypanosoma (Megatrypanum) theileri but they are significantly different from the typical TthI and TthII stocks. Analyses of blood composition indicated that the source of the blood meal in these leeches was from the water buffalo (Bubalus bubalis). To further test if this flagellate from the freshwater leech was indeed a mammalian trypanosome, we transferred the trypanosomes cultured at 27-37 °C and they were able to successfully adapt to this mammalian body temperature, providing further supporting evidence. Due to the significant genetic differences from other related trypanosomes in the subgenus Megatrypanum, we propose that this flagellate, isolated from H. manillensis, is a new species and have named it Trypanosoma bubalisi. Our results indicate that freshwater leeches may be a potential vector of this new mammalian trypanosome.

Ku80 is involved in telomere maintenance but dispensable for genomic stability in Leishmania mexicana

BACKGROUND

Telomeres are indispensable for genome stability maintenance. They are maintained by the telomere-associated protein complex, which include Ku proteins and a telomerase among others. Here, we investigated a role of Ku80 in Leishmania mexicana. Leishmania is a genus of parasitic protists of the family Trypanosomatidae causing a vector-born disease called leishmaniasis.

METHODOLOGY/PRINCIPAL FINDINGS

We used the previously established CRISPR/Cas9 system to mediate ablation of Ku80- and Ku70-encoding genes in L. mexicana. Complete knock-outs of both genes were confirmed by Southern blotting, whole-genome Illumina sequencing, and RT-qPCR. Resulting telomeric phenotypes were subsequently investigated using Southern blotting detection of terminal restriction fragments. The genome integrity in the Ku80- deficient cells was further investigated by whole-genome sequencing. Our work revealed that telomeres in the ΔKu80 L. mexicana are elongated compared to those of the wild type. This is a surprising finding considering that in another model trypanosomatid, Trypanosoma brucei, they are shortened upon ablation of the same gene. A telomere elongation phenotype has been documented in other species and associated with a presence of telomerase-independent alternative telomere lengthening pathway. Our results also showed that Ku80 appears to be not involved in genome stability maintenance in L. mexicana.

CONCLUSION/SIGNIFICANCE

Ablation of the Ku proteins in L. mexicana triggers telomere elongation, but does not have an adverse impact on genome integrity.

Catalase impairs Leishmania mexicana development and virulence

Catalase is one of the most abundant enzymes on Earth. It decomposes hydrogen peroxide, thus protecting cells from dangerous reactive oxygen species. The catalase-encoding gene is conspicuously absent from the genome of most representatives of the family Trypanosomatidae. Here, we expressed this protein from the Leishmania mexicana Β-TUBULIN locus using a novel bicistronic expression system, which relies on the 2A peptide of Teschovirus A. We demonstrated that catalase-expressing parasites are severely compromised in their ability to develop in insects, to be transmitted and to infect mice, and to cause clinical manifestation in their mammalian host. Taken together, our data support the hypothesis that the presence of catalase is not compatible with the dixenous life cycle of Leishmania, resulting in loss of this gene from the genome during the evolution of these parasites.

Highly flexible metabolism of the marine euglenozoan protist Diplonema papillatum

BACKGROUND

The phylum Euglenozoa is a group of flagellated protists comprising the diplonemids, euglenids, symbiontids, and kinetoplastids. The diplonemids are highly abundant and speciose, and recent tools have rendered the best studied representative, Diplonema papillatum, genetically tractable. However, despite the high diversity of diplonemids, their lifestyles, ecological functions, and even primary energy source are mostly unknown.

RESULTS

We designed a metabolic map of D. papillatum cellular bioenergetic pathways based on the alterations of transcriptomic, proteomic, and metabolomic profiles obtained from cells grown under different conditions. Comparative analysis in the nutrient-rich and nutrient-poor media, as well as the absence and presence of oxygen, revealed its capacity for extensive metabolic reprogramming that occurs predominantly on the proteomic rather than the transcriptomic level. D. papillatum is equipped with fundamental metabolic routes such as glycolysis, gluconeogenesis, TCA cycle, pentose phosphate pathway, respiratory complexes, β-oxidation, and synthesis of fatty acids. Gluconeogenesis is uniquely dominant over glycolysis under all surveyed conditions, while the TCA cycle represents an eclectic combination of standard and unusual enzymes.

CONCLUSIONS

The identification of conventional anaerobic enzymes reflects the ability of this protist to survive in low-oxygen environments. Furthermore, its metabolism quickly reacts to restricted carbon availability, suggesting a high metabolic flexibility of diplonemids, which is further reflected in cell morphology and motility, correlating well with their extreme ecological valence.

A resource for improved predictions of Trypanosoma and Leishmania protein three-dimensional structure

AlphaFold2 and RoseTTAfold represent a transformative advance for predicting protein structure. They are able to make very high-quality predictions given a high-quality alignment of the protein sequence with related proteins. These predictions are now readily available via the AlphaFold database of predicted structures and AlphaFold or RoseTTAfold Colaboratory notebooks for custom predictions. However, predictions for some species tend to be lower confidence than model organisms. Problematic species include Trypanosoma cruzi and Leishmania infantum: important unicellular eukaryotic human parasites in an early-branching eukaryotic lineage. The cause appears to be due to poor sampling of this branch of life (Discoba) in the protein sequences databases used for the AlphaFold database and ColabFold. Here, by comprehensively gathering openly available protein sequence data for Discoba species, significant improvements to AlphaFold2 protein structure prediction over the AlphaFold database and ColabFold are demonstrated. This is made available as an easy-to-use tool for the parasitology community in the form of Colaboratory notebooks for generating multiple sequence alignments and AlphaFold2 predictions of protein structure for Trypanosoma, Leishmania and related species.

Trypanosomatid Richness in Wild and Synanthropic Small Mammals from a Biological Station in Rio de Janeiro, Brazil

Trypanosomatids are diverse and can infect several host species, including small mammals (rodents and marsupials). Between 2012 and 2014, 91 small mammals were surveyed for trypanosomatid infection in the Estação Biológica FIOCRUZ Mata Atlântica (EFMA), an Atlantic Forest area in Rio de Janeiro that presents different levels of conserved and degraded areas. Blood, skin, liver, and spleen samples were submitted to parasitological, serological, and molecular assays to detect the infection and determine the taxonomic status of their parasites. Sixty-eight individuals (74.7%; n = 91) were infected by trypanosomatids, including fourteen mixed infected by different trypanosomatid parasites. These hosts were infected by: T. cruzi DTU TcI (n = 12), T. cruzi DTU TcIV (n = 2), T. janseni (n = 15), T. dionisii (n = 1), and T. rangeli A (n = 1) detected in blood or tissue cultures, in addition to T. cruzi DTU TcI (n = 9) and Leishmania sp. (n = 1) only by the molecular diagnosis. Serological diagnosis was positive in 38 (71.6%) individuals for T. cruzi, the same amount for Leishmania spp., and 23 (43.3%) individuals were mixed infected. These data indicate a remarkable richness of trypanosomatid species/genotypes infecting small mammals, even in a disturbed area with low mammal species diversity—as is the case of the EFMA—reinforcing the generalist aspect of these parasites.

Diversity and Epidemiology of Bat Trypanosomes: A One Health Perspective

Bats (order Chiroptera) have been increasingly recognised as important reservoir hosts for human and animal pathogens worldwide. In this context, molecular and microscopy-based investigations to date have revealed remarkably high diversity of Trypanosoma spp. harboured by bats, including species of recognised medical and veterinary importance such as Trypanosoma cruzi and Trypanosoma evansi (aetiological agents of Chagas disease and Surra, respectively). This review synthesises current knowledge on the diversity, taxonomy, evolution and epidemiology of bat trypanosomes based on both molecular studies and morphological records. In addition, we use a One Health approach to discuss the significance of bats as reservoirs (and putative vectors) of T. cruzi, with a focus on the complex associations between intra-specific genetic diversity and eco-epidemiology of T. cruzi in sylvatic and domestic ecosystems. This article also highlights current knowledge gaps on the biological implications of trypanosome co-infections in a single host, as well as the prevalence, vectors, life-cycle, host-range and clinical impact of most bat trypanosomes recorded to date. Continuous research efforts involving molecular surveillance of bat trypanosomes are required for improved disease prevention and control, mitigation of biosecurity risks and potential spill-over events, ultimately ensuring the health of humans, domestic animals and wildlife globally.

A New Model Trypanosomatid, Novymonas esmeraldas: Genomic Perception of Its "Candidatus Pandoraea novymonadis" Endosymbiont

The closest relative of human pathogen Leishmania, the trypanosomatid Novymonas esmeraldas, harbors a bacterial endosymbiont “Candidatus Pandoraea novymonadis.” Based on genomic data, we performed a detailed characterization of the metabolic interactions of both partners. While in many respects the metabolism of N. esmeraldas resembles that of other Leishmaniinae, the endosymbiont provides the trypanosomatid with heme, essential amino acids, purines, some coenzymes, and vitamins. In return, N. esmeraldas shares with the bacterium several nonessential amino acids and phospholipids. Moreover, it complements its carbohydrate metabolism and urea cycle with enzymes missing from the “Ca. Pandoraea novymonadis” genome. The removal of the endosymbiont from N. esmeraldas results in a significant reduction of the overall translation rate, reduced expression of genes involved in lipid metabolism and mitochondrial respiratory activity, and downregulation of several aminoacyl-tRNA synthetases, enzymes involved in the synthesis of some amino acids, as well as proteins associated with autophagy. At the same time, the genes responsible for protection against reactive oxygen species and DNA repair become significantly upregulated in the aposymbiotic strain of this trypanosomatid. By knocking out a component of its flagellum, we turned N. esmeraldas into a new model trypanosomatid that is amenable to genetic manipulation using both conventional and CRISPR-Cas9-mediated approaches.

First description of Lotmaria passim and Crithidia mellificae haptomonad stages in the honeybee hindgut

The remodelling of flagella into attachment structures is a common and important event in the trypanosomatid life cycle. Lotmaria passim and Crithidia mellificae can parasitize Apis mellifera, and as a result they might have a significant impact on honeybee health. However, there are details of their life cycle and the mechanisms underlying their pathogenicity in this host that remain unclear. Here we show that both L. passim promastigotes and C. mellificae choanomastigotes differentiate into haptomonad stages covering the ileum and rectum of honeybees. These haptomonad cells remain attached to the host surface via zonular hemidesmosome-like structures, as revealed by transmission electron microscopy. This work describes for the first known time the haptomonad morphotype of these species and their hemidesmosome-like attachments in A. mellifera, a key trait used by other trypanosomatid species to proliferate in the insect host hindgut.

Experimental Selection of Paromomycin Resistance in Leishmania donovani Amastigotes Induces Variable Genomic Polymorphisms

The relatively high post-treatment relapse rates of paromomycin (PMM) in visceral leishmaniasis treatment and the swift emergence of experimental drug resistance challenge its broad application and urge for rational use and monitoring of resistance. However, no causal molecular mechanisms to Leishmania PMM resistance have been identified so far. To gain insights into potential resistance mechanisms, twelve experimentally selected Leishmania donovani clonal lines and the non-cloned preselection population, with variable degrees of PMM resistance, were subjected to whole genome sequencing. To identify genomic variations potentially associated with resistance, SNPs, Indels, chromosomal somy and gene copy number variations were compared between the different parasite lines. A total of 11 short nucleotide variations and the copy number alterations in 39 genes were correlated to PMM resistance. Some of the identified genes are involved in transcription, translation and protein turn-over (transcription elongation factor-like protein, RNA-binding protein, ribosomal protein L1a, 60S ribosomal protein L6, eukaryotic translation initiation factor 4E-1, proteasome regulatory non-ATP-ase subunit 3), virulence (major surface protease gp63, protein-tyrosine phosphatase 1-like protein), mitochondrial function (ADP/ATP mitochondrial carrier-like protein), signaling (phosphatidylinositol 3-related kinase, protein kinase putative and protein-tyrosine phosphatase 1-like protein) and vesicular trafficking (ras-related protein RAB1). These results indicate that, in Leishmania, the aminoglycoside PMM affects protein translational processes and underlines the complex and probably multifactorial origin of resistance.

Basic Biology of Trypanosoma cruzi

The present review addresses basic aspects of the biology of the pathogenic protozoa Trypanosoma cruzi and some comparative information of Trypanosoma brucei. Like eukaryotic cells, their cellular organization is similar to that of mammalian hosts. However, these parasites present structural particularities. That is why the following topics are emphasized in this paper: developmental stages of the life cycle in the vertebrate and invertebrate hosts; the cytoskeleton of the protozoa, especially the sub-pellicular microtubules; the flagellum and its attachment to the protozoan body through specialized junctions; the kinetoplast-mitochondrion complex, including its structural organization and DNA replication; glycosome and its role in the metabolism of the cell; acidocalcisome, describing its morphology, biochemistry, and functional role; cytostome and the endocytic pathway; the organization of the endoplasmic reticulum and Golgi complex; the nucleus, describing its structural organization during interphase and division; and the process of interaction of the parasite with host cells. The unique characteristics of these structures also make them interesting chemotherapeutic targets. Therefore, further understanding of cell biology aspects contributes to the development of drugs for chemotherapy.

The Complete Mitochondrial DNA of Trypanosoma cruzi: Maxicircles and Minicircles

The mitochondrial DNA of Trypanosomatids, known as the kinetoplast DNA or kDNA or mtDNA, consists of a few maxicircles and thousands of minicircles concatenated together into a huge complex network. These structures present species-specific sizes, from 20 to 40 Kb in maxicircles and from 0.5 to 10 Kb in minicircles. Maxicircles are equivalent to other eukaryotic mitochondrial DNAs, while minicircles contain coding guide RNAs involved in U-insertion/deletion editing processes exclusive of Trypanosomatids that produce the maturation of the maxicircle-encoded transcripts. The knowledge about this mitochondrial genome is especially relevant since the expression of nuclear and mitochondrial genes involved in oxidative phosphorylation must be coordinated. In Trypanosoma cruzi (T. cruzi), the mtDNA has a dual relevance; the production of energy, and its use as a phylogenetic marker due to its high conservation among strains. Therefore, this study aimed to assemble, annotate, and analyze the complete repertoire of maxicircle and minicircle sequences of different T. cruzi strains by using DNA sequencing. We assembled and annotated the complete maxicircle sequence of the Y and Bug2148 strains. For Bug2148, our results confirm that the maxicircle sequence is the longest assembled to date, and is composed of 21 genes, most of them conserved among Trypanosomatid species. In agreement with previous results, T. cruzi minicircles show a conserved structure around 1.4 Kb, with four highly conserved regions and other four hypervariable regions interspersed between them. However, our results suggest that the parasite minicircles display several sizes and numbers of conserved and hypervariable regions, contrary to those previous studies. Besides, this heterogeneity is also reflected in the three conserved sequence blocks of the conserved regions that play a key role in the minicircle replication. Our results using sequencing technologies of second and third-generation indicate that the different consensus sequences of the maxicircles and minicircles seem to be more complex than previously described indicating at least four different groups in T. cruzi minicircles.

Endosymbiont Capture, a Repeated Process of Endosymbiont Transfer with Replacement in Trypanosomatids Angomonas spp

Trypanosomatids of the subfamily Strigomonadinae bear permanent intracellular bacterial symbionts acquired by the common ancestor of these flagellates. However, the cospeciation pattern inherent to such relationships was revealed to be broken upon the description of Angomonas ambiguus, which is sister to A. desouzai, but bears an endosymbiont genetically close to that of A. deanei. Based on phylogenetic inferences, it was proposed that the bacterium from A. deanei had been horizontally transferred to A. ambiguus. Here, we sequenced the bacterial genomes from two A. ambiguus isolates, including a new one from Papua New Guinea, and compared them with the published genome of the A. deanei endosymbiont, revealing differences below the interspecific level. Our phylogenetic analyses confirmed that the endosymbionts of A. ambiguus were obtained from A. deanei and, in addition, demonstrated that this occurred more than once. We propose that coinfection of the same blowfly host and the phylogenetic relatedness of the trypanosomatids facilitate such transitions, whereas the drastic difference in the occurrence of the two trypanosomatid species determines the observed direction of this process. This phenomenon is analogous to organelle (mitochondrion/plastid) capture described in multicellular organisms and, thereafter, we name it endosymbiont capture.

Mitochondrial RNA quality control in trypanosomes

Unicellular parasites Trypanosoma brucei spp. cause African human and animal trypanosomiasis, a spectrum of diseases that jeopardize public health and afflict the economy in sub-Saharan Africa. These hemoflagellates are distinguished by a single mitochondrion, which contains a kinetoplast nucleoid composed of DNA and histone-like proteins. Kinetoplast DNA (kDNA) represents a densely packed network of interlinked relaxed circular molecules: a few ~23-kb maxicircles encoding ribosomal RNAs (rRNAs) and proteins, and approximately 5,000 1-kb minicircles bearing guide RNA (gRNA) genes. The transcription start site defines the mRNA's 5' terminus while the primary RNA is remodeled into a monocistronic messenger by 3'-5' exonucleolytic trimming, 5' and 3' end modifications, and, in most cases, by internal U-insertion/deletion editing. Ribosomal and guide RNA precursors are also trimmed, and the processed molecules are uridylated. For 35 years, mRNA editing has attracted a major effort, but more recently the essential pre- and postediting processing and turnover events have been discovered and the key effectors have been identified. Among these, pentatricopeptide repeat (PPR) RNA binding proteins emerged as conduits coupling modifications of mRNA termini with internal sequence changes introduced by editing. Among 39 annotated PPRs, 20 belong to ribosomal subunits or assembly intermediates, four function as polyadenylation factors, a single factor directs 5' mRNA modification, and one protein is found in F1-ATPase. Nuclear and mitochondrial RNases P consist of a single PPR polypeptide, PRORP1 and PROP2, respectively. Here, we review PPR-mediated mitochondrial processes and discuss their potential roles in mRNA maturation, quality control, translational activation, and decay. This article is categorized under: RNA Processing > Capping and 5' End Modifications RNA Processing > 3' End Processing RNA Processing > RNA Editing and Modification.

Crithidia mellificae infection in different mammalian species in Brazil

Crithidia mellificae, a monoxenous trypanosomatid considered restricted to insects, was recently reported to infect a bat. Herein, C. mellificae has been demonstrated to have a wider range of vertebrate hosts and distribution in Brazilian biomes than once thought. Parasites isolated from haemocultures were characterized using V7V8 SSU rDNA and glyceraldehyde 3-phosphate dehydrogenase genes. Coatis (Nasua nasua) in the Cerrado; marmosets (Callithrix sp.) and bats (Carollia perspicillata, Myotis lavali, M. izecksohni, Artibeus lituratus) in the Atlantic Forest; crab-eating foxes (Cerdocyon thous) and ocelot (Leopardus pardalis) in the Pantanal biomes were infected by trypanosomatids that displayed choanomastigote forms in haemoculture in Giemsa-stained slide smears. Molecular characterization and phylogenetic inference confirmed the infection of C. mellificae in these animals. Moreover, slight differences in C. mellificae sequences were observed. Crithidia mellificae growth curves were counted at 27°C, 36°C and 37°C, and the morphotypes were able to grow and survive for up to 16 days. Serological titers for C. mellificae were observed in nonhuman primates, demonstrating that this parasite is able to induce a humoral immune response in an infected mammal. These results showed that host specificity in trypanosomatids is complex and far from understood.