Comparative genomics of trypanosomatid parasitic protozoa

Innovation through instability? Genome (dis)organisation in Trypanosoma cruzi

Chagas disease affects millions globally and is caused by Trypanosoma cruzi, a parasite with a remarkable ability to adapt and persist in complex environments. Emerging evidence suggests that this adaptability may be driven by its extensive genome instability yet the underlying mechanisms remain poorly understood. Drawing insights from related trypanosomatids, we explore how processes such as aneuploidy, copy number variations and genetic rearrangements can be strategically exploited for immune evasion, host adaptation and drug resistance, and consider whether similar processes shape the adaptive potential of T. cruzi. Unravelling the fundamental biology governing genome instability in this parasite will undoubtedly advance our understanding of its evolution, pathogenesis and resilience - crucial steps towards reducing the burden of Chagas disease.

PCR-based diagnosis of Surra using a newly identified conserved region of the variant surface glycoprotein (VSG) gene

Trypanosoma evansi, the causative agent of Surra, is the most widespread pathogenic trypanosome that parasitizes the widest variety of mammals worldwide; however, rapid and easily accessible diagnostics still need to be improved. Proteomic research identified the VSG (Variant Surface Glycoprotein) TevSTIB805.3.100 gene as a potential biomarker for T. evansi diagnosis. The aim of this study was to design primers (called Tev3.100) for the gene encoding this VSG, testing the specificity and sensitivity of these primers on genomic DNA (gDNA) from different species and on biological samples. The specificity of primers was tested against gDNA from T. evansi, T. brucei, T. equiperdum, T. rangeli, T. cruzi, T. vivax, Babesia bovis, B. bigemina, and Anaplasma sp. Seventy-one biological samples from Lageana Creole cattle DNA were used, testing the sensitivity, specificity and concordance ​​in relation to RoTat 1.2 primers. The Tev3.100 primers were able to produce amplicons with a single band of approximately 1800 bp for gDNA from T. evansi, but showed cross-reactions with T. brucei, and T. equiperdum, diverging from the in silico predictions. These primers indicated high sensitivity (98.28 %) and specificity (84.62 %) in the detection of biological samples from Lageano Creole cattle, in addition to high concordance values (κ: 0.854; SE: 0.082; 95 % CI: 0.695-1.000) in relation to RoTat primers results. The Tev3.100 primers are a new molecular tool with good sensitivity and specificity for Surra infections, but the cross-reactions with T. equiperdum, diverging from the databases, indicate that new genomic studies should be carried out for these species in Latin America.

Exploring a series of multifunctional Mn(I) tricarbonyls as prospective agents against trypanosomatid parasites: a comparative study with the Re(I) analogues

Diseases caused by trypanosomatid parasites are among the most pressing neglected illnesses. Chagas disease, caused by Trypanosoma cruzi, and visceral Leishmaniasis, caused by Leishmania infantum, have a severe health impact in developing countries. Searching for prospective metal-based drugs against these diseases, five multifunctional fac-[Mn(CO)3(CTZ)(NN)](PF6) compounds, including four new derivatives, were synthesized and thoroughly characterized, featuring NN polypyridyl derivatives and Clotrimazole (CTZ) as bioactive ligands. The biological behavior was compared with that previously reported for the Re analogues. Mn compounds showed EC50 values in the low micromolar range against the infective trypomastigote form of Trypanosoma cruzi and the promastigote form of Leishmania infantum and moderate selectivity indexes. While their potency against T. cruzi was comparable to the Re analogues, their selectivity was lower. Key physicochemical properties relevant to drug development were assessed: Mn(I) compounds showed lower stability in relevant tested media compared with their Re(I) counterparts and higher lipophilicity than the free ligands and the Re analogues. To gain insight into the potential mechanisms of action, the interaction with DNA and the effects on ergosterol biosynthesis in T. cruzi and L. infantum were investigated. Minimal DNA association (<1%) and moderate interaction with this target discarded DNA binding as the primary mechanism of action. In contrast, inhibition of lanosterol 14-α-demethylase (CYP51), key enzyme involved in the parasites' ergosterol biosynthetic pathway, was experimentally confirmed. Metallomic study revealed an uptake by T. cruzi of the most promising compound, fac-[Mn(CO)3(CTZ)(tmp)](PF6), more than twice that of the Re(I) analogue and preferential association to soluble proteins. Proteomic analysis of T. cruzi epimastigotes treated with the Mn(I) and Re(I) analogues showed no change in CYP51 abundance, suggesting that reduced ergosterol levels may arise from post-translational modifications of the enzyme. Raman confocal microscopy allowed us to detect effects of the most promising Mn compound in treated T. cruzi. Furthermore, the photoinduced CO release properties of both Mn and Re analogues were examined, searching for an additional and yet non-studied potential mechanism of action of metal-tricarbonyls in these trypanosomatid parasites. Collectively, the results highlight the potential of Mn(I) tricarbonyls as promising candidates for further drug development.

An Overview of Trypanosoma cruzi Biology Through the Lens of Proteomics: A Review

The protozoan parasite Trypanosoma cruzi, causative agent of Chagas disease, affects millions of people in endemic Latin American countries and beyond. In Latin America, Chagas disease is an important cause of death and disability, for which vaccines are lacking and improved treatment options are required. Additionally, the factors governing the development of a variety of clinical manifestations during Chagas disease, ranging from complete lack of symptoms to severe irreversible chronic organ damage (mainly cardiac or digestive), remain largely unknown. Much remains to be learned regarding the biology of T. cruzi in order to enhance our understanding of these lines of inquiry. In this context, proteomic methods have been leveraged to investigate different parasite strains, life-cycle forms, subcellular compartments, macromolecular complexes, signaling events and secreted molecules. The factors driving morphological transformation during the life cycle, the composition and functions of the parasite's organelles and secreted molecules as well as the determinants of pathogenicity have been explored via proteomic methods, yielding insights into the fundamental processes behind the parasite biology and informing drug design and vaccine development. Importantly, the correlation between the wide genetic and phenotypic variability displayed by T. cruzi has been examined through proteomic methods as well. Here, we review the literature on T. cruzi proteomics and discuss it in the light of its limitations and in the context of the parasite's genetic diversity.

A limitation lifted: A conditional knockdown system reveals essential roles for Polo-like kinase and Aurora kinase 1 in Trypanosoma cruzi cell division

While advances in genome editing technologies have simplified gene disruption in many organisms, the study of essential genes requires development of conditional disruption or knockdown systems that are not available in most organisms. Such is the case for Trypanosoma cruzi, a parasite that causes Chagas disease, a severely neglected tropical disease endemic to Latin America that is often fatal. Our knowledge of the identity of essential genes and their functions in T. cruzi has been severely constrained by historical challenges in very basic genetic manipulation and the absence of RNA interference machinery. Here, we describe the development and use of self-cleaving RNA sequences to conditionally regulate essential gene expression in T. cruzi. Using these tools, we identified essential roles for Polo-like and Aurora kinases in T. cruzi cell division, mirroring their functions in Trypanosoma brucei. Importantly, we demonstrate conditional knockdown of essential genes in intracellular amastigotes, the disease-causing stage of the parasite in its human host. This conditional knockdown system enables the efficient and scalable functional characterization of essential genes in T. cruzi and provides a framework for the development of conditional gene knockdown systems for other nonmodel organisms.

R-loops acted on by RNase H1 influence DNA replication timing and genome stability in Leishmania

Genomes in eukaryotes normally undergo DNA replication in a choreographed temporal order, resulting in early and late replicating chromosome compartments. Leishmania, a human protozoan parasite, displays an unconventional DNA replication program in which the timing of DNA replication completion is chromosome size-dependent: larger chromosomes complete replication later then smaller ones. Here we show that both R-loops and RNase H1, a ribonuclease that resolves RNA-DNA hybrids, accumulate in Leishmania major chromosomes in a pattern that reflects their replication timing. Furthermore, we demonstrate that such differential organisation of R-loops, RNase H1 and DNA replication timing across the parasite's chromosomes correlates with size-dependent differences in chromatin accessibility, G quadruplex distribution and sequence content. Using conditional gene excision, we show that loss of RNase H1 leads to transient growth perturbation and permanently abrogates the differences in DNA replication timing across chromosomes, as well as altering levels of aneuploidy and increasing chromosome instability in a size-dependent manner. This work provides a link between R-loop homeostasis and DNA replication timing in a eukaryotic parasite and demonstrates that orchestration of DNA replication dictates levels of genome plasticity in Leishmania.

Trypanosoma cruzi: Genomic Diversity and Structure

Trypanosoma cruzi is the causative agent of Chagas disease, a neglected tropical disease, and one of the most important parasitic diseases worldwide. The first genome of T. cruzi was sequenced in 2005, and its complexity made assembly and annotation challenging. Nowadays, new sequencing methods have improved some strains' genome sequence and annotation, revealing this parasite's extensive genetic diversity and complexity. In this review, we examine the genetic diversity, the genomic structure, and the principal multi-gene families involved in the pathogenicity of T. cruzi. The T. cruzi genome sequence is divided into two compartments: the core (conserved) and the disruptive (variable in length and multicopy gene families among strains). The disruptive region has also been described as genome plasticity and plays a key role in the parasite survival and infection process. This region comprises several multi-gene families, including trans-sialidases, mucins, and mucin-associated surface proteins (MASPs). Trans-sialidases are the most prevalent genes in the genome with a key role in the infection process, while mucins and MASPs are also significant glycosylated proteins expressed on the parasite surface, essential for its biological functions, as host-parasite interaction, host cell invasion or protection against the host immune system, in both insect and mammalian stages. Collectively, in this review, some of the most recent advances in the structure and composition of the T. cruzi genome are reviewed.

Characterization of Tau95 led to the identification of a four-subunit TFIIIC complex in trypanosomatid parasites

RNA polymerase III (RNAP III) synthetizes small essential non-coding RNA molecules such as tRNAs and 5S rRNA. In yeast and vertebrates, RNAP III needs general transcription factors TFIIIA, TFIIIB, and TFIIIC to initiate transcription. TFIIIC, composed of six subunits, binds to internal promoter elements in RNAP III-dependent genes. Limited information is available about RNAP III transcription in the trypanosomatid protozoa Trypanosoma brucei and Leishmania major, which diverged early from the eukaryotic lineage. Analyses of the first published draft of the trypanosomatid genome sequences failed to recognize orthologs of any of the TFIIIC subunits, suggesting that this transcription factor is absent in these parasites. However, a putative TFIIIC subunit was recently annotated in the databases. Here we characterize this subunit in T. brucei and L. major and demonstrate that it corresponds to Tau95. In silico analyses showed that both proteins possess the typical Tau95 sequences: the DNA binding region and the dimerization domain. As anticipated for a transcription factor, Tau95 localized to the nucleus in insect forms of both parasites. Chromatin immunoprecipitation (ChIP) assays demonstrated that Tau95 binds to tRNA and U2 snRNA genes in T. brucei. Remarkably, by performing tandem affinity purifications we identified orthologs of TFIIIC subunits Tau55, Tau131, and Tau138 in T. brucei and L. major. Thus, contrary to what was assumed, trypanosomatid parasites do possess a TFIIIC complex. Other putative interacting partners of Tau95 were identified in T. brucei and L. major. KEY POINTS: • A four-subunit TFIIIC complex is present in T. brucei and L. major • TbTau95 associates with tRNA and U2 snRNA genes • Putative interacting partners of Tau95 might include some RNAP II regulators.

Analyses of the essential C82 subunit uncovered some differences in RNA polymerase III transcription between Trypanosoma brucei and Leishmania major

The 17-subunit RNA polymerase III (RNAP III) synthesizes essential untranslated RNAs such as tRNAs and 5S rRNA. In yeast and vertebrates, subunit C82 forms a stable subcomplex with C34 and C31 that is necessary for promoter-specific transcription initiation. Little is known about RNAP III transcription in trypanosomatid parasites. To narrow this knowledge gap, we characterized the C82 subunit in Trypanosoma brucei and Leishmania major. Bioinformatic analyses showed that the 4 distinctive extended winged-helix (eWH) domains and the coiled-coil motif are present in C82 in these microorganisms. Nevertheless, C82 in trypanosomatids presents certain unique traits, including an exclusive loop within the eWH1 domain. We found that C82 localizes to the nucleus and binds to RNAP III-dependent genes in the insect stages of both parasites. Knock-down of C82 by RNA interference significantly reduced the levels of tRNAs and 5S rRNA and led to the death of procyclic forms of T. brucei. Tandem affinity purifications with both parasites allowed the identification of several C82-interacting partners, including C34 and some genus-specific putative regulators of transcription. However, the orthologue of C31 was not found in trypanosomatids. Interestingly, our data suggest a strong association of C82 with TFIIIC subunits in T. brucei, but not in L. major.

Insights into Peptidyl-Prolyl cis-trans Isomerases from Clinically Important Protozoans: From Structure to Potential Biotechnological Applications

Peptidyl-prolyl cis/trans isomerases (PPIases) are present in a wide variety of microorganisms, including protozoan parasites such as Trypanosoma cruzi, Trypanosoma brucei, Trichomonas vaginalis, Leishmania major, Leishmania donovani, Plasmodium falciparum, Plasmodium vivax, Entamoeba histolytica, Giardia intestinalis, Cryptosporidium parvum, and Cryptosporidium hominis, all of which cause important neglected diseases. PPIases are classified as cyclophilins, FKBPs, or parvulins and play crucial roles in catalyzing the cis-trans isomerization of the peptide bond preceding a proline residue. This activity assists in correct protein folding. However, experimentally, the biological structure-function characterization of PPIases from these protozoan parasites has been poorly addressed. The recombinant production of these enzymes is highly relevant for this ongoing research. Thus, this review explores the structural diversity, functions, recombinant production, activity, and inhibition of protozoan PPIases. We also highlight their potential as biotechnological tools for the in vitro refolding of other recombinant proteins from these parasites. These applications are invaluable for the development of diagnostic and therapeutic tools.

Protein phosphatase PP1 regulation of RNA polymerase II transcription termination and allelic exclusion of VSG genes in trypanosomes

The genomes of Leishmania and trypanosomes are organized into polycistronic transcription units flanked by a modified DNA base J involved in promoting RNA polymerase II (Pol II) termination. We recently characterized a Leishmania complex containing a J-binding protein, PP1 protein phosphatase 1, and PP1 regulatory protein (PNUTS) that controls transcription termination potentially via dephosphorylation of Pol II by PP1. While T. brucei contains eight PP1 isoforms, none purified with the PNUTS complex, complicating the analysis of PP1 function in termination. We now demonstrate that the PP1-binding motif of TbPNUTS is required for function in termination in vivo and that TbPP1-1 modulates Pol II termination in T. brucei and dephosphorylation of the large subunit of Pol II. PP1-1 knock-down results in increased cellular levels of phosphorylated RPB1 accompanied by readthrough transcription and aberrant transcription of the chromosome by Pol II, including Pol I transcribed loci that are typically silent, such as telomeric VSG expression sites involved in antigenic variation. These results provide important insights into the mechanism underlying Pol II transcription termination in primitive eukaryotes that rely on polycistronic transcription and maintain allelic exclusion of VSG genes.

Integrated genomics provides insights into the evolution of the polyphosphate accumulation trait of Ca. Accumulibacter

Candidatus Accumulibacter, a prominent polyphosphate-accumulating organism (PAO) in wastewater treatment, plays a crucial role in enhanced biological phosphorus removal (EBPR). The genetic underpinnings of its polyphosphate accumulation capabilities, however, remain largely unknown. Here, we conducted a comprehensive genomic analysis of Ca. Accumulibacter-PAOs and their relatives within the Rhodocyclaceae family, identifying 124 core genes acquired via horizontal gene transfer (HGT) at its least common ancestor. Metatranscriptomic analysis of an enrichment culture of Ca. Accumulibacter revealed active transcription of 44 of these genes during an EBPR cycle, notably including the polyphosphate kinase 2 (PPK2) gene instead of the commonly recognized polyphosphate kinase 1 (PPK1) gene. Intriguingly, the phosphate regulon (Pho) genes showed minimal transcriptions, pointing to a distinctive fact of Pho dysregulation, where PhoU, the phosphate signaling complex protein, was not regulating the high-affinity phosphate transport (Pst) system, resulting in continuous phosphate uptake. To prevent phosphate toxicity, Ca. Accumulibacter utilized the laterally acquired PPK2 to condense phosphate into polyphosphate, resulting in the polyphosphate-accumulating feature. This study provides novel insights into the evolutionary emergence of the polyphosphate-accumulating trait in Ca. Accumulibacter, offering potential advancements in understanding the PAO phenotype in the EBPR process.

Multifunctional Organometallic Compounds Active against Infective Trypanosomes: Ru(II) Ferrocenyl Derivatives with Two Different Bioactive Ligands

Human African trypanosomiasis (HAT, sleeping sickness) and American trypanosomiasis (Chagas disease) are endemic zoonotic diseases caused by genomically related trypanosomatid protozoan parasites (Trypanosoma brucei and Trypanosoma cruzi, respectively). Just a few old drugs are available for their treatment, with most of them sharing poor safety, efficacy, and pharmacokinetic profiles. Only fexinidazole has been recently incorporated into the arsenal for the treatment of HAT. In this work, new multifunctional Ru(II) ferrocenyl compounds were rationally designed as potential agents against these pathogens by including in a single molecule 1,1'-bis(diphenylphosphino)ferrocene (dppf) and two bioactive bidentate ligands: pyridine-2-thiolato-1-oxide ligand (mpo) and polypyridyl ligands (NN). Three [Ru(mpo)(dppf)(NN)](PF6) compounds and their derivatives with chloride as a counterion were synthesized and fully characterized in solid state and solution. They showed in vitro activity on bloodstream T. brucei (EC50 = 31-160 nM) and on T. cruzi trypomastigotes (EC50 = 190-410 nM). Compounds showed the lowest EC50 values on T. brucei when compared to the whole set of metal-based compounds previously developed by us. In addition, several of the Ru compounds showed good selectivity toward the parasites, particularly against the highly proliferative bloodstream form of T. brucei. Interaction with DNA and generation of reactive oxygen species (ROS) were ruled out as potential targets and modes of action of the Ru compounds. Biochemical assays and in silico analysis led to the insight that they are able to inhibit the NADH-dependent fumarate reductase from T. cruzi. One representative hit induced a mild oxidation of low molecular weight thiols in T. brucei. The compounds were stable for at least 72 h in two different media and more lipophilic than both bioactive ligands, mpo and NN. An initial assessment of the therapeutic efficacy of one of the most potent and selective candidates, [Ru(mpo)(dppf)(bipy)]Cl, was performed using a murine infection model of acute African trypanosomiasis. This hit compound lacks acute toxicity when applied to animals in the dose/regimen described, but was unable to control parasite proliferation in vivo, probably because of its rapid clearance or low biodistribution in the extracellular fluids. Future studies should investigate the pharmacokinetics of this compound in vivo and involve further research to gain deeper insight into the mechanism of action of the compounds.

Genome-wide distribution of 5-hydroxymethyluracil and chromatin accessibility in the Breviolum minutum genome

BACKGROUND

In dinoflagellates, a unique and extremely divergent genomic and nuclear organization has evolved. The highly unusual features of dinoflagellate nuclei and genomes include permanently condensed liquid crystalline chromosomes, primarily packaged by proteins other than histones, genes organized in very long unidirectional gene arrays, a general absence of transcriptional regulation, high abundance of the otherwise very rare DNA modification 5-hydroxymethyluracil (5-hmU), and many others. While most of these fascinating properties are originally identified in the 1970s and 1980s, they have not yet been investigated using modern genomic tools.

RESULTS

In this work, we address some of the outstanding questions regarding dinoflagellate genome organization by mapping the genome-wide distribution of 5-hmU (using both immunoprecipitation-based and basepair-resolution chemical mapping approaches) and of chromatin accessibility in the genome of the Symbiodiniaceae dinoflagellate Breviolum minutum. We find that the 5-hmU modification is preferentially enriched over certain classes of repetitive elements, often coincides with the boundaries between gene arrays, and is generally correlated with decreased chromatin accessibility, the latter otherwise being largely uniform along the genome. We discuss the potential roles of 5-hmU in the functional organization of dinoflagellate genomes and its relationship to the transcriptional landscape of gene arrays.

CONCLUSIONS

Our results provide the first window into the 5-hmU and chromatin accessibility landscapes in dinoflagellates.

Molecular tools to regulate gene expression in Trypanosoma cruzi

Developing molecular strategies to manipulate gene expression in trypanosomatids is challenging, particularly with respect to the unique gene expression mechanisms adopted by these unicellular parasites, such as polycistronic mRNA transcription and multi-gene families. In the case of Trypanosoma cruzi (T. cruzi), the causative agent of Chagas Disease, the lack of RNA interference machinery further complicated functional genetic studies important for understanding parasitic biology and developing biomarkers and potential therapeutic targets. Therefore, alternative methods of performing knockout and/or endogenous labelling experiments were developed to identify and understand the function of proteins for survival and interaction with the host. In this review, we present the main tools for the genetic manipulation of T. cruzi, focusing on the Clustered Regularly Interspaced Short Palindromic Repeats Cas9-associated system technique widely used in this organism. Moreover, we highlight the importance of using these tools to elucidate the function of uncharacterized and glycosylated proteins. Further developments of these technologies will allow the identification of new biomarkers, therapeutic targets and potential vaccines against Chagas disease with greater efficiency and speed.

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.

Unwrap RAP1's Mystery at Kinetoplastid Telomeres

Although located at the chromosome end, telomeres are an essential chromosome component that helps maintain genome integrity and chromosome stability from protozoa to mammals. The role of telomere proteins in chromosome end protection is conserved, where they suppress various DNA damage response machineries and block nucleolytic degradation of the natural chromosome ends, although the detailed underlying mechanisms are not identical. In addition, the specialized telomere structure exerts a repressive epigenetic effect on expression of genes located at subtelomeres in a number of eukaryotic organisms. This so-called telomeric silencing also affects virulence of a number of microbial pathogens that undergo antigenic variation/phenotypic switching. Telomere proteins, particularly the RAP1 homologs, have been shown to be a key player for telomeric silencing. RAP1 homologs also suppress the expression of Telomere Repeat-containing RNA (TERRA), which is linked to their roles in telomere stability maintenance. The functions of RAP1s in suppressing telomere recombination are largely conserved from kinetoplastids to mammals. However, the underlying mechanisms of RAP1-mediated telomeric silencing have many species-specific features. In this review, I will focus on Trypanosoma brucei RAP1's functions in suppressing telomeric/subtelomeric DNA recombination and in the regulation of monoallelic expression of subtelomere-located major surface antigen genes. Common and unique mechanisms will be compared among RAP1 homologs, and their implications will be discussed.

Fifteen Years after the Definition of Trypanosoma cruzi DTUs: What Have We Learned?

Trypanosoma cruzi, the protozoan causative of Chagas disease (ChD), exhibits striking genetic and phenotypic intraspecific diversity, along with ecoepidemiological complexity. Human-pathogen interactions lead to distinct clinical presentations of ChD. In 2009, an international consensus classified T. cruzi strains into six discrete typing units (DTUs), TcI to TcVI, later including TcBat, and proposed reproducible genotyping schemes for DTU identification. This article aims to review the impact of classifying T. cruzi strains into DTUs on our understanding of biological, ecoepidemiological, and pathogenic aspects of T. cruzi. We will explore the likely origin of DTUs and the intrinsic characteristics of each group of strains concerning genome organization, genomics, and susceptibility to drugs used in ChD treatment. We will also provide an overview of the association of DTUs with mammalian reservoirs, and summarize the geographic distribution, and the clinical implications, of prevalent specific DTUs in ChD patients. Throughout this review, we will emphasize the crucial roles of both parasite and human genetics in defining ChD pathogenesis and chemotherapy outcome.

Leishmania PNUTS discriminates between PP1 catalytic subunits through an RVxF-ΦΦ-F motif and polymorphisms in the PP1 C-tail and catalytic domain

Phosphoprotein phosphatase 1 (PP1) associates with specific regulatory subunits to achieve, among other functions, substrate selectivity. Among the eight PP1 isotypes in Leishmania, PP1-8e associates with the regulatory protein PNUTS along with the structural factors JBP3 and Wdr82 in the PJW/PP1 complex that modulates RNA polymerase II (pol II) phosphorylation and transcription termination. Little is known regarding interactions involved in PJW/PP1 complex formation, including how PP1-8e is the selective isotype associated with PNUTS. Here, we show that PNUTS uses an established RVxF-ΦΦ-F motif to bind the PP1 catalytic domain with similar interfacial interactions as mammalian PP1-PNUTS and noncanonical motifs. These atypical interactions involve residues within the PP1-8e catalytic domain and N and C terminus for isoform-specific regulator binding. This work advances our understanding of PP1 isoform selectivity and reveals key roles of PP1 residues in regulator binding. We also explore the role of PNUTS as a scaffold protein for the complex by identifying the C-terminal region involved in binding JBP3 and Wdr82 and impact of PNUTS on the stability of complex components and function in pol II transcription in vivo. Taken together, these studies provide a potential mechanism where multiple motifs within PNUTS are used combinatorially to tune binding affinity to PP1, and the C terminus for JBP3 and Wdr82 association, in the Leishmania PJW/PP1 complex. Overall, our data provide insights in the formation of the PJW/PP1 complex involved in regulating pol II transcription in divergent protozoans where little is understood.

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.

Morphogenesis in Trypanosoma cruzi epimastigotes proceeds via a highly asymmetric cell division

Trypanosoma cruzi is a protist parasite that is the causative agent of Chagas disease, a neglected tropical disease endemic to the Americas. T. cruzi cells are highly polarized and undergo morphological changes as they cycle within their insect and mammalian hosts. Work on related trypanosomatids has described cell division mechanisms in several life-cycle stages and identified a set of essential morphogenic proteins that serve as markers for key events during trypanosomatid division. Here, we use Cas9-based tagging of morphogenic genes, live-cell imaging, and expansion microscopy to study the cell division mechanism of the insect-resident epimastigote form of T. cruzi, which represents an understudied trypanosomatid morphotype. We find that T. cruzi epimastigote cell division is highly asymmetric, producing one daughter cell that is significantly smaller than the other. Daughter cell division rates differ by 4.9 h, which may be a consequence of this size disparity. Many of the morphogenic proteins identified in T. brucei have altered localization patterns in T. cruzi epimastigotes, which may reflect fundamental differences in the cell division mechanism of this life cycle stage, which widens and shortens the cell body to accommodate the duplicated organelles and cleavage furrow rather than elongating the cell body along the long axis of the cell, as is the case in life-cycle stages that have been studied in T. brucei. This work provides a foundation for further investigations of T. cruzi cell division and shows that subtle differences in trypanosomatid cell morphology can alter how these parasites divide.

Protein Phosphatase PP1 Regulation of Pol II Phosphorylation is Linked to Transcription Termination and Allelic Exclusion of VSG Genes and TERRA in Trypanosomes

The genomes of Leishmania and trypanosomes are organized into polycistronic transcription units flanked by a modified DNA base J involved in promoting RNA polymerase II (Pol II) termination. We recently characterized a Leishmania complex containing a J-binding protein, PP1 protein phosphatase 1, and PP1 regulatory protein (PNUTS) that controls transcription termination potentially via dephosphorylation of Pol II by PP1. While T. brucei contains eight PP1 isoforms, none purified with the PNUTS complex, suggesting a unique PP1-independent mechanism of termination. We now demonstrate that the PP1-binding motif of TbPNUTS is required for function in termination in vivo and that TbPP1-1 modulates Pol II termination in T. brucei involving dephosphorylation of the C-terminal domain of the large subunit of Pol II. PP1-1 knock-down results in increased cellular levels of phosphorylated large subunit of Pol II accompanied by readthrough transcription and pervasive transcription of the entire genome by Pol II, including Pol I transcribed loci that are typically silent, such as telomeric VSG expression sites involved in antigenic variation and production of TERRA RNA. These results provide important insights into the mechanism underlying Pol II transcription termination in primitive eukaryotes that rely on polycistronic transcription and maintain allelic exclusion of VSG genes.

Genome-wide distribution of 5-hydroxymethyluracil and chromatin accessibility in the Breviolum minutum genome

In dinoflagellates, a unique and extremely divergent genomic and nuclear organization has evolved. The highly unusual features of dinoflagellate nuclei and genomes include permanently condensed liquid crystalline chromosomes, primarily packaged by proteins other than histones, genes organized in very long unidirectional gene arrays, a general absence of transcriptional regulation, high abundance of the otherwise very rare DNA modification 5-hydroxymethyluracil (5-hmU), and many others. While most of these fascinating properties were originally identified in the 1970s and 1980s, they have not yet been investigated using modern genomic tools. In this work, we address some of the outstanding questions regarding dinoflagellate genome organization by mapping the genome-wide distribution of 5-hmU (using both immunoprecipitation-based and basepair-resolution chemical mapping approaches) and of chromatin accessibility in the genome of the Symbiodiniaceae dinoflagellate Breviolum minutum. We find that the 5-hmU modification is preferentially enriched over certain classes of repetitive elements, often coincides with the boundaries between gene arrays, and is generally correlated with decreased chromatin accessibility, the latter otherwise being largely uniform along the genome. We discuss the potential roles of 5-hmU in the functional organization of dinoflagellate genomes and its relationship to the transcriptional landscape of gene arrays.

Leishmania PNUTS discriminates between PP1 catalytic subunits through a RVxF-ΦΦ-F motif and polymorphisms in the PP1 C-tail and catalytic domain

PP1 phosphatases lack substrate specificity and associate with specific regulatory subunits to achieve selectivity. Among the eight PP1 isotypes in Leishmania, PP1-8e associates with the regulatory protein PNUTS along with the structural factors JBP3 and Wdr82 in the PJW/PP1 complex that modulates RNA polymerase II (Pol II) phosphorylation and transcription termination. Little is known regarding interactions involved in PJW/PP1 complex formation, including how PP1-8e is the selective isotype associated with PNUTS. Here, we show that PNUTS uses an established RVxF-ΦΦ-F motif to bind the PP1 catalytic domain with similar interfacial interactions as mammalian PP1- PNUTS and non-canonical motifs. These atypical interactions involve residues within the PP1-8e catalytic domain and N- and C-terminus for isoform specific regulator binding. This work advances our understanding of PP1 isoform selectivity and reveals key roles of PP1 residues in regulator binding. We also explore the role of PNUTS as a scaffold protein for the complex by identifying the C-terminal region involved in binding JBP3 and Wdr82, and impact of PNUTS on the stability of complex components and function in Pol II transcription in vivo . Taken together, these studies provide a potential mechanism where multiple motifs within PNUTS are used combinatorially to tune binding affinity to PP1, and the C-termini for independent binding of JBP3 and Wdr82, in the Leishmania PJW/PP1 complex. Overall, our data provide insights in the formation of the PJW/PP1 complex involved in regulating Pol II transcription in divergent protozoans where little is understood.

Transporter Proteins as Ecological Assets and Features of Microbial Eukaryotic Pangenomes

Here we review two connected themes in evolutionary microbiology: (a) the nature of gene repertoire variation within species groups (pangenomes) and (b) the concept of metabolite transporters as accessory proteins capable of providing niche-defining "bolt-on" phenotypes. We discuss the need for improved sampling and understanding of pangenome variation in eukaryotic microbes. We then review the factors that shape the repertoire of accessory genes within pangenomes. As part of this discussion, we outline how gene duplication is a key factor in both eukaryotic pangenome variation and transporter gene family evolution. We go on to outline how, through functional characterization of transporter-encoding genes, in combination with analyses of how transporter genes are gained and lost from accessory genomes, we can reveal much about the niche range, the ecology, and the evolution of virulence of microbes. We advocate for the coordinated systematic study of eukaryotic pangenomes through genome sequencing and the functional analysis of genes found within the accessory gene repertoire.

mt-LAF3 is a pseudouridine synthase ortholog required for mitochondrial rRNA and mRNA gene expression in Trypanosoma brucei

Trypanosoma brucei and related kinetoplastid parasites possess unique RNA processing pathways, including in their mitochondria, that regulate metabolism and development. Altering RNA composition or conformation through nucleotide modifications is one such pathway, and modifications including pseudouridine regulate RNA fate and function in many organisms. We surveyed pseudouridine synthase (PUS) orthologs in trypanosomatids, with a particular interest in mitochondrial enzymes due to their potential importance for mitochondrial function and metabolism. Trypanosoma brucei mitochondrial (mt)-LAF3 is an ortholog of human and yeast mitochondrial PUS enzymes, and a mitoribosome assembly factor, but structural studies differ in their conclusion as to whether it has PUS catalytic activity. Here, we generated T. brucei cells that are conditionally null (CN) for mt-LAF3 expression and showed that mt-LAF3 loss is lethal and disrupts mitochondrial membrane potential (ΔΨm). Addition of a mutant gamma ATP synthase allele to the CN cells permitted ΔΨm maintenance and cell survival, allowing us to assess primary effects on mitochondrial RNAs. As expected, these studies showed that loss of mt-LAF3 dramatically decreases levels of mitochondrial 12S and 9S rRNAs. Notably, we also observed decreases in mitochondrial mRNA levels, including differential effects on edited vs. pre-edited mRNAs, indicating that mt-LAF3 is required for mitochondrial rRNA and mRNA processing, including of edited transcripts. To assess the importance of PUS catalytic activity in mt-LAF3 we mutated a conserved aspartate that is necessary for catalysis in other PUS enzymes and showed it is not essential for cell growth, or maintenance of ΔΨm and mitochondrial RNA levels. Together, these results indicate that mt-LAF3 is required for normal expression of mitochondrial mRNAs in addition to rRNAs, but that PUS catalytic activity is not required for these functions. Instead, our work, combined with previous structural studies, suggests that T. brucei mt-LAF3 acts as a mitochondrial RNA-stabilizing scaffold.

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.

Knockout of protein phosphatase 1 in Leishmania major reveals its role during RNA polymerase II transcription termination

The genomes of kinetoplastids are organized into polycistronic transcription units that are flanked by a modified DNA base (base J, beta-D-glucosyl-hydroxymethyluracil). Previous work established a role of base J in promoting RNA polymerase II (Pol II) termination in Leishmania major and Trypanosoma brucei. We recently identified a PJW/PP1 complex in Leishmania containing a J-binding protein (JBP3), PP1 phosphatase 1, PP1 interactive-regulatory protein (PNUTS) and Wdr82. Analyses suggested the complex regulates transcription termination by recruitment to termination sites via JBP3-base J interactions and dephosphorylation of proteins, including Pol II, by PP1. However, we never addressed the role of PP1, the sole catalytic component, in Pol II transcription termination. We now demonstrate that deletion of the PP1 component of the PJW/PP1 complex in L. major, PP1-8e, leads to readthrough transcription at the 3'-end of polycistronic gene arrays. We show PP1-8e has in vitro phosphatase activity that is lost upon mutation of a key catalytic residue and associates with PNUTS via the conserved RVxF motif. Additionally, purified PJW complex with associated PP1-8e, but not complex lacking PP1-8e, led to dephosphorylation of Pol II, suggesting a direct role of PNUTS/PP1 holoenzymes in regulating transcription termination via dephosphorylating Pol II in the nucleus.

State-of-the-Art in the Drug Discovery Pathway for Chagas Disease: A Framework for Drug Development and Target Validation

Chagas disease is the most important protozoan infection in the Americas, and constitutes a significant public health concern throughout the world. Development of new medications against its etiologic agent, Trypanosoma cruzi, has been traditionally slow and difficult, lagging in comparison with diseases caused by other kinetoplastid parasites. Among the factors that explain this are the incompletely understood mechanisms of pathogenesis of T. cruzi infection and its complex set of interactions with the host in the chronic stage of the disease. These demand the performance of a variety of in vitro and in vivo assays as part of any drug development effort. In this review, we discuss recent breakthroughs in the understanding of the parasite's life cycle and their implications in the search for new chemotherapeutics. For this, we present a framework to guide drug discovery efforts against Chagas disease, considering state-of-the-art preclinical models and recently developed tools for the identification and validation of molecular targets.

Enhanced apoptotic index in hepatocytes, Kupffer cells, and inflammatory infiltrate showed positive correlation with hepatic lesion intensity, parasite load, and clinical status in naturally Leishmania-infected dogs

It is unknown if Leishmania amastigote infections affect hepatocytes and Kupffer cell apoptosis, and the role played by apoptosis in liver lesions in leishmaniasis is still unclear. Clinically affected and subclinically infected dogs with leishmaniosis and uninfected controls were assessed. Parasite load, biochemical markers for evaluation of liver damage, morphometry (area, perimeter, number of inflammatory focus, major and minor diameters), apoptosis in hepatic tissue (hepatocytes, Kupffer cells, and inflammatory infiltrates) and cellularity in inflammatory foci were quantified. The parasite load in clinically affected dogs proved to be higher than in the other groups. All morphometric parameters (area, perimeter, number of inflammatory focus, major and minor diameters) from clinically affected were higher than the values found in the subclinically infected and uninfected control dogs. Only clinically affected dogs presented high levels of ALT, FA, GGT and cholesterol in serum. Strong positive correlation was observed between biochemical markers for evaluation of liver damage (ALT, FA, GGT and cholesterol) and hepatic apoptosis (hepatocytes, Kupffer cells, and inflammation). Clinically affected dogs showed a more intense hepatic lesion. Hepatocytes showed a higher rate of apoptosis in Leishmania-infected dogs than in uninfected control dogs. The Kupffer cell apoptotic index and apoptosis within the inflammatory infiltrates were higher in clinically affected dogs. The apoptotic index evaluated in hepatocytes, Kupffer cells, and inflammatory infiltrates showed a positive correlation with the intensity of the hepatic lesion, parasite load, and clinical status. Apoptotic cells also showed positive immunostaining for TUNEL, Bcl2, and Bax. Our data showed that hepatic apoptosis was related to the severity of liver damage, the progression of infection, and the parasite load in leishmaniasis. Apoptotic regulated cell recruitment modulated the inflammatory response and favored the survival and dissemination of parasites, depending on the clinical status of the Leishmania-infected dogs.

Morphogenesis in Trypanosoma cruzi epimastigotes proceeds via a highly asymmetric cell division

Trypanosoma cruzi is a protist parasite that is the causative agent of Chagas' disease, a neglected tropical disease endemic to the Americas. T. cruzi cells are highly polarized and undergo morphological changes as they cycle within their insect and mammalian hosts. Work on related trypanosomatids has described cell division mechanisms in several life-cycle stages and identified a set of essential morphogenic proteins that serve as markers for key events during trypanosomatid division. Here, we use Cas9-based tagging of morphogenic genes, live-cell imaging, and expansion microscopy to study the cell division mechanism of the insect-resident epimastigote form of T. cruzi, which represents an understudied trypanosomatid morphotype. We find that T. cruzi epimastigote cell division is highly asymmetric, producing one daughter cell that is significantly smaller than the other. Daughter cell division rates differ by 4.9 h, which may be a consequence of this size disparity. Many of the morphogenic proteins identified in T. brucei have altered localization patterns in T. cruzi epimastigoes, which may reflect fundamental differences in the cell division mechanism of this life cycle stage, which widens and shortens the cell body to accommodate the duplicated organelles and cleavage furrow rather than elongating the cell body along the long axis of the cell, as is the case in life-cycle stages that have been studied in T. brucei. This work provides a foundation for further investigations of T. cruzi cell division and shows that subtle differences in trypansomatid cell morphology can alter how these parasites divide.

Genome-Wide Libraries for Protozoan Pathogen Drug Target Screening Using Yeast Surface Display

The lack of genetic tools to manipulate protozoan pathogens has limited the use of genome-wide approaches to identify drug or vaccine targets and understand these organisms' biology. We have developed an efficient method to construct genome-wide libraries for yeast surface display (YSD) and developed a YSD fitness screen (YSD-FS) to identify drug targets. We show the efficacy of our method by generating genome-wide libraries for Trypanosoma brucei, Trypanosoma cruzi, and Giardia lamblia parasites. Each library has a diversity of ∼105 to 106 clones, representing ∼6- to 30-fold of the parasite's genome. Nanopore sequencing confirmed the libraries' genome coverage with multiple clones for each parasite gene. Western blot and imaging analysis confirmed surface expression of the G. lamblia library proteins in yeast. Using the YSD-FS assay, we identified bonafide interactors of metronidazole, a drug used to treat protozoan and bacterial infections. We also found enrichment in nucleotide-binding domain sequences associated with yeast increased fitness to metronidazole, indicating that this drug might target multiple enzymes containing nucleotide-binding domains. The libraries are valuable biological resources for discovering drug or vaccine targets, ligand receptors, protein-protein interactions, and pathogen-host interactions. The library assembly approach can be applied to other organisms or expression systems, and the YSD-FS assay might help identify new drug targets in protozoan pathogens.

Omics data integration facilitates target selection for new antiparasitic drugs against TriTryp infections

Introduction:Trypanosoma cruzi, Trypanosoma brucei, and Leishmania spp., commonly referred to as TriTryps, are a group of protozoan parasites that cause important human diseases affecting millions of people belonging to the most vulnerable populations worldwide. Current treatments have limited efficiencies and can cause serious side effects, so there is an urgent need to develop new control strategies. Presently, the identification and prioritization of appropriate targets can be aided by integrative genomic and computational approaches.

Methods: In this work, we conducted a genome-wide multidimensional data integration strategy to prioritize drug targets. We included genomic, transcriptomic, metabolic, and protein structural data sources, to delineate candidate proteins with relevant features for target selection in drug development.

Results and Discussion: Our final ranked list includes proteins shared by TriTryps and covers a range of biological functions including essential proteins for parasite survival or growth, oxidative stress-related enzymes, virulence factors, and proteins that are exclusive to these parasites. Our strategy found previously described candidates, which validates our approach as well as new proteins that can be attractive targets to consider during the initial steps of drug discovery.

Advances in Protozoan Epigenetic Targets and Their Inhibitors for the Development of New Potential Drugs

Protozoan parasite diseases cause significant mortality and morbidity worldwide. Factors such as climate change, extreme poverty, migration, and a lack of life opportunities lead to the propagation of diseases classified as tropical or non-endemic. Although there are several drugs to combat parasitic diseases, strains resistant to routinely used drugs have been reported. In addition, many first-line drugs have adverse effects ranging from mild to severe, including potential carcinogenic effects. Therefore, new lead compounds are needed to combat these parasites. Although little has been studied regarding the epigenetic mechanisms in lower eukaryotes, it is believed that epigenetics plays an essential role in vital aspects of the organism, from controlling the life cycle to the expression of genes involved in pathogenicity. Therefore, using epigenetic targets to combat these parasites is foreseen as an area with great potential for development. This review summarizes the main known epigenetic mechanisms and their potential as therapeutics for a group of medically important protozoal parasites. Different epigenetic mechanisms are discussed, highlighting those that can be used for drug repositioning, such as histone post-translational modifications (HPTMs). Exclusive parasite targets are also emphasized, including the base J and DNA 6 mA. These two categories have the greatest potential for developing drugs to treat or eradicate these diseases.

Ligand-based virtual screening, molecular dynamics, and biological evaluation of repurposed drugs as inhibitors of Trypanosoma cruzi proteasome

Chagas disease is a well-known Neglected Tropical Disease, mostly endemic in continental Latin America, but that has spread to North America and Europe. Unfortunately, current treatments against such disease are ineffective and produce known and undesirable side effects. To find novel effective drug candidates to treat Chagas disease, we uniquely explore the Trypanosoma cruzi proteasome as a recent biological target and, also, apply drug repurposing through different computational methodologies. For this, we initially applied protein homology modeling to build a robust model of proteasome β4/β5 subunits, since there is no crystallographic structure of this target. Then, we used it on a drug repurposing via a virtual screening campaign starting with more than 8,000 drugs and including the methodologies: ligand-based similarity, toxicity predictions, and molecular docking. Three drugs were selected concerning their favorable interactions at the protein binding site and subsequently submitted to molecular dynamics simulations, which allowed us to elucidate their behavior and compare such theoretical results with experimental ones, obtained in biological assays also described in this paper.Communicated by Ramaswamy H. Sarma.

mt-LAF3 is a pseudouridine synthase ortholog required for mitochondrial rRNA and mRNA gene expression in Trypanosoma brucei

Trypanosoma brucei and related kinetoplastid parasites possess unique RNA processing pathways, including in their mitochondria, that regulate metabolism and development. Altering RNA composition or conformation through nucleotide modifications is one such pathway, and modifications including pseudouridine regulate RNA fate and function in many organisms. We surveyed pseudouridine synthase (PUS) orthologs in Trypanosomatids, with a particular interest in mitochondrial enzymes due to their potential importance for mitochondrial function and metabolism. T. brucei mt-LAF3 is an ortholog of human and yeast mitochondrial PUS enzymes, and a mitoribosome assembly factor, but structural studies differ in their conclusion as to whether it has PUS catalytic activity. Here, we generated T. brucei cells that are conditionally null for mt-LAF3 and showed that mt-LAF3 loss is lethal and disrupts mitochondrial membrane potential (ΔΨm). Addition of a mutant gamma-ATP synthase allele to the conditionally null cells permitted ΔΨm maintenance and cell survival, allowing us to assess primary effects on mitochondrial RNAs. As expected, these studies showed that loss of mt-LAF3 dramatically decreases levels of mitochondrial 12S and 9S rRNAs. Notably, we also observed decreases in mitochondrial mRNA levels, including differential effects on edited vs. pre-edited mRNAs, indicating that mt-LAF3 is required for mitochondrial rRNA and mRNA processing, including of edited transcripts. To assess the importance of PUS catalytic activity in mt-LAF3 we mutated a conserved aspartate that is necessary for catalysis in other PUS enzymes and showed it is not essential for cell growth, or maintenance of ΔΨm and mitochondrial RNA levels. Together, these results indicate that mt-LAF3 is required for normal expression of mitochondrial mRNAs in addition to rRNAs, but that PUS catalytic activity is not required for these functions. Instead, our work, combined with previous structural studies, suggests that T. brucei mt-LAF3 acts as a mitochondrial RNA-stabilizing scaffold.

The paradigm of intracellular parasite survival and drug resistance in leishmanial parasite through genome plasticity and epigenetics: Perception and future perspective

Leishmania is an intracellular, zoonotic, kinetoplastid eukaryote with more than 1.2 million cases all over the world. The leishmanial chromosomes are divided into polymorphic chromosomal ends, conserved central domains, and antigen-encoding genes found in telomere-proximal regions. The genome flexibility of chromosomal ends of the leishmanial parasite is known to cause drug resistance and intracellular survival through the evasion of host defense mechanisms. Therefore, in this review, we discuss the plasticity of Leishmania genome organization which is the primary cause of drug resistance and parasite survival. Moreover, we have not only elucidated the causes of such genome plasticity which includes aneuploidy, epigenetic factors, copy number variation (CNV), and post-translation modification (PTM) but also highlighted their impact on drug resistance and parasite survival.

Revisiting epidemiology of leishmaniasis in central Asia: lessons learnt

In this work we reviewed historical and recent data on Leishmania spp. infection combining data collected in Turkmenistan, Uzbekistan, Kazakhstan, Kyrgyzstan, Iran, China and Mongolia. We specifically focused on a complex of co-existing species (Leishmania major, Leishmania turanica and Leishmania gerbilli) sharing the same animal reservoirs and vectors. In addition, we analysed the presence of dsRNA viruses in these species and discussed future research directions to identify species-specific traits, which may determine susceptibility of different Leishmania spp. to viral infection.

Accessing the Variability of Multicopy Genes in Complex Genomes using Unassembled Next-Generation Sequencing Reads: The Case of Trypanosoma cruzi Multigene Families

Repetitive elements cause assembly fragmentation in complex eukaryotic genomes, limiting the study of their variability. The genome of Trypanosoma cruzi, the parasite that causes Chagas disease, has a high repetitive content, including multigene families. Although many T. cruzi multigene families encode surface proteins that play pivotal roles in host-parasite interactions, their variability is currently underestimated, as their high repetitive content results in collapsed gene variants. To estimate sequence variability and copy number variation of multigene families, we developed a read-based approach that is independent of gene-specific read mapping and de novo assembly. This methodology was used to estimate the copy number and variability of MASP, TcMUC, and Trans-Sialidase (TS), the three largest T. cruzi multigene families, in 36 strains, including members of all six parasite discrete typing units (DTUs). We found that these three families present a specific pattern of variability and copy number among the distinct parasite DTUs. Inter-DTU hybrid strains presented a higher variability of these families, suggesting that maintaining a larger content of their members could be advantageous. In addition, in a chronic murine model and chronic Chagasic human patients, the immune response was focused on TS antigens, suggesting that targeting TS conserved sequences could be a potential avenue to improve diagnosis and vaccine design against Chagas disease. Finally, the proposed approach can be applied to study multicopy genes in any organism, opening new avenues to access sequence variability in complex genomes. IMPORTANCE Sequences that have several copies in a genome, such as multicopy-gene families, mobile elements, and microsatellites, are among the most challenging genomic segments to study. They are frequently underestimated in genome assemblies, hampering the correct assessment of these important players in genome evolution and adaptation. Here, we developed a new methodology to estimate variability and copy numbers of repetitive genomic regions and employed it to characterize the T. cruzi multigene families MASP, TcMUC, and transsialidase (TS), which are important virulence factors in this parasite. We showed that multigene families vary in sequence and content among the parasite's lineages, whereas hybrid strains have a higher sequence variability that could be advantageous to the parasite's survivability. By identifying conserved sequences within multigene families, we showed that the mammalian host immune response toward these multigene families is usually focused on the TS multigene family. These TS conserved and immunogenic peptides can be explored in future works as diagnostic targets or vaccine candidates for Chagas disease. Finally, this methodology can be easily applied to any organism of interest, which will aid in our understanding of complex genomic regions.

New multifunctional Ru(II) organometallic compounds show activity against Trypanosoma brucei and Leishmania infantum

Human African trypanosomiasis (sleeping sickness) and leishmaniasis are prevalent zoonotic diseases caused by genomically related trypanosomatid protozoan parasites (Trypanosoma brucei and Leishmania spp). Additionally, both are co-endemic in certain regions of the world. Only a small number of old drugs exist for their treatment, with most of them sharing poor safety, efficacy, and pharmacokinetic profiles. In this work, new multifunctional Ru(II) ferrocenyl compounds were rationally designed as potential agents against these trypanosomatid parasites by including in a single molecule 1,1'-bis(diphenylphosphino)ferrocene (dppf) and two bioactive bidentate ligands: 8-hydroxyquinoline derivatives (8HQs) and polypyridyl ligands (NN). Three [Ru(8HQs)(dppf)(NN)](PF₆) compounds were synthesized and fully characterized. They showed in vitro activity on bloodstream Trypanosoma brucei (IC₅₀ 140-310 nM) and on Leishmania infantum promastigotes (IC₅₀ 3.0-4.8 μM). The compounds showed good selectivity towards T. brucei in respect to J774 murine macrophages as mammalian cell model (SI 15-38). Changing hexafluorophosphate counterion by chloride led to a three-fold increase in activity on both parasites and to a two to three-fold increase in selectivity towards the pathogens. The compounds affect in vitro at least the targets of the individual bioactive moieties included in the new chemical entities: DNA and generation of ROS. The compounds are stable in solution and are more lipophilic than the free bioactive ligands. No clear correlation between lipophilicity, interaction with DNA or generation of ROS and activity was detected, which agrees with their overall similar anti-trypanosoma potency and selectivity. These compounds are promising candidates for further drug development.

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.

In silico analysis of the HSP90 chaperone system from the African trypanosome, Trypanosoma brucei

African trypanosomiasis is a neglected tropical disease caused by Trypanosoma brucei (T. brucei) and spread by the tsetse fly in sub-Saharan Africa. The trypanosome relies on heat shock proteins for survival in the insect vector and mammalian host. Heat shock protein 90 (HSP90) plays a crucial role in the stress response at the cellular level. Inhibition of its interactions with chaperones and co-chaperones is being explored as a potential therapeutic target for numerous diseases. This study provides an in silico overview of HSP90 and its co-chaperones in both T. brucei brucei and T. brucei gambiense in relation to human and other trypanosomal species, including non-parasitic Bodo saltans and the insect infecting Crithidia fasciculata. A structural analysis of T. brucei HSP90 revealed differences in the orientation of the linker and C-terminal domain in comparison to human HSP90. Phylogenetic analysis displayed the T. brucei HSP90 proteins clustering into three distinct groups based on subcellular localizations, namely, cytosol, mitochondria, and endoplasmic reticulum. Syntenic analysis of cytosolic HSP90 genes revealed that T. b. brucei encoded for 10 tandem copies, while T. b. gambiense encoded for three tandem copies; Leishmania major (L. major) had the highest gene copy number with 17 tandem copies. The updated information on HSP90 from recently published proteomics on T. brucei was examined for different life cycle stages and subcellular localizations. The results show a difference between T. b. brucei and T. b. gambiense with T. b. brucei encoding a total of twelve putative HSP90 genes, while T. b. gambiense encodes five HSP90 genes. Eighteen putative co-chaperones were identified with one notable absence being cell division cycle 37 (Cdc37). These results provide an updated framework on approaching HSP90 and its interactions as drug targets in the African trypanosome.

Delineating transitions during the evolution of specialised peroxisomes: Glycosome formation in kinetoplastid and diplonemid protists

One peculiarity of protists belonging to classes Kinetoplastea and Diplonemea within the phylum Euglenozoa is compartmentalisation of most glycolytic enzymes within peroxisomes that are hence called glycosomes. This pathway is not sequestered in peroxisomes of the third Euglenozoan class, Euglenida. Previous analysis of well-studied kinetoplastids, the ‘TriTryps’ parasites Trypanosoma brucei, Trypanosoma cruzi and Leishmania spp., identified within glycosomes other metabolic processes usually not present in peroxisomes. In addition, trypanosomatid peroxins, i.e. proteins involved in biogenesis of these organelles, are divergent from human and yeast orthologues. In recent years, genomes, transcriptomes and proteomes for a variety of euglenozoans have become available. Here, we track the possible evolution of glycosomes by querying these databases, as well as the genome of Naegleria gruberi, a non-euglenozoan, which belongs to the same protist supergroup Discoba. We searched for orthologues of TriTryps proteins involved in glycosomal metabolism and biogenesis. Predicted cellular location(s) of each metabolic enzyme identified was inferred from presence or absence of peroxisomal-targeting signals. Combined with a survey of relevant literature, we refine extensively our previously postulated hypothesis about glycosome evolution. The data agree glycolysis was compartmentalised in a common ancestor of the kinetoplastids and diplonemids, yet additionally indicates most other processes found in glycosomes of extant trypanosomatids, but not in peroxisomes of other eukaryotes were either sequestered in this ancestor or shortly after separation of the two lineages. In contrast, peroxin divergence is evident in all euglenozoans. Following their gain of pathway complexity, subsequent evolution of peroxisome/glycosome function is complex. We hypothesize compartmentalisation in glycosomes of glycolytic enzymes, their cofactors and subsequently other metabolic enzymes provided selective advantage to kinetoplastids and diplonemids during their evolution in changing marine environments. We contend two specific properties derived from the ancestral peroxisomes were key: existence of nonselective pores for small solutes and the possibility of high turnover by pexophagy. Critically, such pores and pexophagy are characterised in extant trypanosomatids. Increasing amenability of free-living kinetoplastids and recently isolated diplonemids to experimental study means our hypothesis and interpretation of bioinformatic data are suited to experimental interrogation.

Genome-scale RNA interference profiling of Trypanosoma brucei cell cycle progression defects

Trypanosomatids, which include major pathogens of humans and livestock, are flagellated protozoa for which cell cycle controls and the underlying mechanisms are not completely understood. Here, we describe a genome-wide RNA-interference library screen for cell cycle defects in Trypanosoma brucei. We induced massive parallel knockdown, sorted the perturbed population using high-throughput flow cytometry, deep-sequenced RNAi-targets from each stage and digitally reconstructed cell cycle profiles at a genomic scale; also enabling data visualisation using an online tool ( https://tryp-cycle.pages.dev/ ). Analysis of several hundred genes that impact cell cycle progression reveals >100 flagellar component knockdowns linked to genome endoreduplication, evidence for metabolic control of the G₁-S transition, surface antigen regulatory mRNA-binding protein knockdowns linked to G₂M accumulation, and a putative nucleoredoxin required for both mitochondrial genome segregation and for mitosis. The outputs provide comprehensive functional genomic evidence for the known and novel machineries, pathways and regulators that coordinate trypanosome cell cycle progression.

Light chain 2 is a Tctex-type related axonemal dynein light chain that regulates directional ciliary motility in Trypanosoma brucei

Flagellar motility is essential for the cell morphology, viability, and virulence of pathogenic kinetoplastids. Trypanosoma brucei flagella beat with a bending wave that propagates from the flagellum's tip to its base, rather than base-to-tip as in other eukaryotes. Thousands of dynein motor proteins coordinate their activity to drive ciliary bending wave propagation. Dynein-associated light and intermediate chains regulate the biophysical mechanisms of axonemal dynein. Tctex-type outer arm dynein light chain 2 (LC2) regulates flagellar bending wave propagation direction, amplitude, and frequency in Chlamydomonas reinhardtii. However, the role of Tctex-type light chains in regulating T. brucei motility is unknown. Here, we used a combination of bioinformatics, in-situ molecular tagging, and immunofluorescence microscopy to identify a Tctex-type light chain in the procyclic form of T. brucei (TbLC2). We knocked down TbLC2 expression using RNAi in both wild-type and FLAM3, a flagellar attachment zone protein, knockdown cells and quantified TbLC2's effects on trypanosome cell biology and biophysics. We found that TbLC2 knockdown reduced the directional persistence of trypanosome cell swimming, induced an asymmetric ciliary bending waveform, modulated the bias between the base-to-tip and tip-to-base beating modes, and increased the beating frequency. Together, our findings are consistent with a model of TbLC2 as a down-regulator of axonemal dynein activity that stabilizes the forward tip-to-base beating ciliary waveform characteristic of trypanosome cells. Our work sheds light on axonemal dynein regulation mechanisms that contribute to pathogenic kinetoplastids' unique tip-to-base ciliary beating nature and how those mechanisms underlie dynein-driven ciliary motility more generally.

Nucleosome Structures Built from Highly Divergent Histones: Parasites and Giant DNA Viruses

In eukaryotes, genomic DNA is bound with histone proteins and packaged into chromatin. The nucleosome, a fundamental unit of chromatin, regulates the accessibility of DNA to enzymes involved in gene regulation. During the past few years, structural analyses of chromatin architectures have been limited to evolutionarily related organisms. The amino acid sequences of histone proteins are highly conserved from humans to yeasts, but are divergent in the deeply branching protozoan groups, including human parasites that are directly related to human health. Certain large DNA viruses, as well as archaeal organisms, contain distant homologs of eukaryotic histone proteins. The divergent sequences give rise to unique and distinct nucleosome architectures, although the fundamental principles of histone folding and DNA contact are highly conserved. In this article, we review the structures and biophysical properties of nucleosomes containing histones from the human parasites Giardia lamblia and Leishmania major, and histone-like proteins from the Marseilleviridae amoeba virus family. The presented data confirm the sharing of the overall DNA compaction system among evolutionally distant species and clarify the deviations from the species-specific nature of the nucleosome.

Current Status of Regulatory Non-Coding RNAs Research in the Tritryp

Trypanosomatids are protozoan parasites that cause devastating vector-borne human diseases. Gene expression regulation of these organisms depends on post-transcriptional control in responding to diverse environments while going through multiple developmental stages of their complex life cycles. In this scenario, non-coding RNAs (ncRNAs) are excellent candidates for a very efficient, quick, and economic strategy to regulate gene expression. The advent of high throughput RNA sequencing technologies show the presence and deregulation of small RNA fragments derived from canonical ncRNAs. This review seeks to depict the ncRNA landscape in trypanosomatids, focusing on the small RNA fragments derived from functional RNA molecules observed in RNA sequencing studies. Small RNA fragments derived from canonical ncRNAs (tsRNAs, snsRNAs, sdRNAs, and sdrRNAs) were identified in trypanosomatids. Some of these RNAs display changes in their levels associated with different environments and developmental stages, demanding further studies to determine their functional characterization and potential roles. Nevertheless, a comprehensive and detailed ncRNA annotation for most trypanosomatid genomes is still needed, allowing better and more extensive comparative and functional studies.

A novel SNF2 ATPase complex in Trypanosoma brucei with a role in H2A.Z-mediated chromatin remodelling

A cascade of histone acetylation events with subsequent incorporation of a histone H2A variant plays an essential part in transcription regulation in various model organisms. A key player in this cascade is the chromatin remodelling complex SWR1, which replaces the canonical histone H2A with its variant H2A.Z. Transcriptional regulation of polycistronic transcription units in the unicellular parasite Trypanosoma brucei has been shown to be highly dependent on acetylation of H2A.Z, which is mediated by the histone-acetyltransferase HAT2. The chromatin remodelling complex which mediates H2A.Z incorporation is not known and an SWR1 orthologue in trypanosomes has not yet been reported. In this study, we identified and characterised an SWR1-like remodeller complex in T. brucei that is responsible for Pol II-dependent transcriptional regulation. Bioinformatic analysis of potential SNF2 DEAD/Box helicases, the key component of SWR1 complexes, identified a 1211 amino acids-long protein that exhibits key structural characteristics of the SWR1 subfamily. Systematic protein-protein interaction analysis revealed the existence of a novel complex exhibiting key features of an SWR1-like chromatin remodeller. RNAi-mediated depletion of the ATPase subunit of this complex resulted in a significant reduction of H2A.Z incorporation at transcription start sites and a subsequent decrease of steady-state mRNA levels. Furthermore, depletion of SWR1 and RNA-polymerase II (Pol II) caused massive chromatin condensation. The potential function of several proteins associated with the SWR1-like complex and with HAT2, the key factor of H2A.Z incorporation, is discussed.

Trypanosoma brucei is the causative agent of African trypanosomiasis (sleeping sickness) in humans and nagana in cattle. Its unusual genomic organisation featuring large polycistronic units requires a general mechanism of transcription initiation, because individual gene promoters are mostly absent. Despite the fact that the histone variant H2A.Z has previously been identified as a key player of transcription regulation, the complex responsible for correct H2A.Z incorporation at transcription start sites (TSS) remains elusive. In other eukaryotes, SWR1, a SNF2 ATPase-associated chromatin remodelling complex, is responsible for correct incorporation of this histone variant. This study identified a SWR1-like complex in T. brucei. Depletion of the SNF2 ATPase resulted in a reduction of H2A.Z incorporation at the TSS and decreased steady-state mRNA levels accompanied by chromatin condensation. In addition to the SWR1-like complex, we also identified a trypanosome-specific HAT2 complex that includes the histone acetyltransferases HAT2, a key player in the H2A.Z incorporation process. This complex has a trypanosome-specific composition that is different from the NuA4/TIP60 complex in Saccharomyces cerevisiae.