Trans-splicing in trypanosomes: machinery and its impact on the parasite transcriptome

Identification of diverse RNA viruses in Obscuromonas flagellates (Euglenozoa: Trypanosomatidae: Blastocrithidiinae)

Trypanosomatids (Euglenozoa) are a diverse group of unicellular flagellates predominately infecting insects (monoxenous species) or circulating between insects and vertebrates or plants (dixenous species). Monoxenous trypanosomatids harbor a wide range of RNA viruses belonging to the families Narnaviridae, Totiviridae, Qinviridae, Leishbuviridae, and a putative group of tombus-like viruses. Here, we focus on the subfamily Blastocrithidiinae, a previously unexplored divergent group of monoxenous trypanosomatids comprising two related genera: Obscuromonas and Blastocrithidia. Members of the genus Blastocrithidia employ a unique genetic code, in which all three stop codons are repurposed to encode amino acids, with TAA also used to terminate translation. Obscuromonas isolates studied here bear viruses of three families: Narnaviridae, Qinviridae, and Mitoviridae. The latter viral group is documented in trypanosomatid flagellates for the first time. While other known mitoviruses replicate in the mitochondria, those of trypanosomatids appear to reside in the cytoplasm. Although no RNA viruses were detected in Blastocrithidia spp., we identified an endogenous viral element in the genome of B. triatomae indicating its past encounter(s) with tombus-like viruses.

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.

Emerging roles of the epitranscriptome in parasitic protozoan biology and pathogenesis

RNA modifications (epitranscriptome) - such as N6-methyladenosine (m6A), 5-methylcytosine (m5C), and pseudouridine (Ψ) - modulate RNA processing, stability, interaction, and translation, thereby playing critical roles in the development, replication, virulence, metabolism, and life cycle adaptations of parasitic protozoa. Here, we summarize potential homologs of the major human RNA modification regulatory factors in parasites, outline current knowledge on how RNA modifications affect parasitic protozoa, highlight the regulation of RNA modifications and their crosstalk, and discuss current progress in exploring RNA modifications as potential drug targets. This review contributes to our understanding of epitranscriptomic regulation of parasitic protozoa biology and pathogenesis and provides new perspectives for the treatment of parasitic diseases.

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.

A single pseudouridine on rRNA regulates ribosome structure and function in the mammalian parasite Trypanosoma brucei

Trypanosomes are protozoan parasites that cycle between insect and mammalian hosts and are the causative agent of sleeping sickness. Here, we describe the changes of pseudouridine (Ψ) modification on rRNA in the two life stages of the parasite using four different genome-wide approaches. CRISPR-Cas9 knock-outs of all four snoRNAs guiding Ψ on helix 69 (H69) of the large rRNA subunit were lethal. A single knock-out of a snoRNA guiding Ψ530 on H69 altered the composition of the 80S monosome. These changes specifically affected the translation of only a subset of proteins. This study correlates a single site Ψ modification with changes in ribosomal protein stoichiometry, supported by a high-resolution cryo-EM structure. We propose that alteration in rRNA modifications could generate ribosomes preferentially translating state-beneficial proteins.

How much (ATP) does it cost to build a trypanosome? A theoretical study on the quantity of ATP needed to maintain and duplicate a bloodstream-form Trypanosoma brucei cell

ATP hydrolysis is required for the synthesis, transport and polymerization of monomers for macromolecules as well as for the assembly of the latter into cellular structures. Other cellular processes not directly related to synthesis of biomass, such as maintenance of membrane potential and cellular shape, also require ATP. The unicellular flagellated parasite Trypanosoma brucei has a complex digenetic life cycle. The primary energy source for this parasite in its bloodstream form (BSF) is glucose, which is abundant in the host's bloodstream. Here, we made a detailed estimation of the energy budget during the BSF cell cycle. As glycolysis is the source of most produced ATP, we calculated that a single parasite produces 6.0 x 1011 molecules of ATP/cell cycle. Total biomass production (which involves biomass maintenance and duplication) accounts for ~63% of the total energy budget, while the total biomass duplication accounts for the remaining ~37% of the ATP consumption, with in both cases translation being the most expensive process. These values allowed us to estimate a theoretical YATP of 10.1 (g biomass)/mole ATP and a theoretical [Formula: see text] of 28.6 (g biomass)/mole ATP. Flagellar motility, variant surface glycoprotein recycling, transport and maintenance of transmembrane potential account for less than 30% of the consumed ATP. Finally, there is still ~5.5% available in the budget that is being used for other cellular processes of as yet unknown cost. These data put a new perspective on the assumptions about the relative energetic weight of the processes a BSF trypanosome undergoes during its cell cycle.

Roles of alternative splicing in infectious diseases: from hosts, pathogens to their interactions

ABSTRACT

Alternative splicing (AS) is an evolutionarily conserved mechanism that removes introns and ligates exons to generate mature messenger RNAs (mRNAs), extremely improving the richness of transcriptome and proteome. Both mammal hosts and pathogens require AS to maintain their life activities, and inherent physiological heterogeneity between mammals and pathogens makes them adopt different ways to perform AS. Mammals and fungi conduct a two-step transesterification reaction by spliceosomes to splice each individual mRNA (named cis -splicing). Parasites also use spliceosomes to splice, but this splicing can occur among different mRNAs (named trans -splicing). Bacteria and viruses directly hijack the host's splicing machinery to accomplish this process. Infection-related changes are reflected in the spliceosome behaviors and the characteristics of various splicing regulators (abundance, modification, distribution, movement speed, and conformation), which further radiate to alterations in the global splicing profiles. Genes with splicing changes are enriched in immune-, growth-, or metabolism-related pathways, highlighting approaches through which hosts crosstalk with pathogens. Based on these infection-specific regulators or AS events, several targeted agents have been developed to fight against pathogens. Here, we summarized recent findings in the field of infection-related splicing, including splicing mechanisms of pathogens and hosts, splicing regulation and aberrant AS events, as well as emerging targeted drugs. We aimed to systemically decode host-pathogen interactions from a perspective of splicing. We further discussed the current strategies of drug development, detection methods, analysis algorithms, and database construction, facilitating the annotation of infection-related splicing and the integration of AS with disease phenotype.

Nanopore-Based Direct RNA Sequencing of the Trypanosoma brucei Transcriptome Identifies Novel lncRNAs

Trypanosomatids are single-cell eukaryotic parasites. Unlike higher eukaryotes, they control gene expression post-transcriptionally and not at the level of transcription initiation. This involves all known cellular RNA circuits, from mRNA processing to mRNA decay, to translation, in addition to a large panel of RNA-interacting proteins that modulate mRNA abundance. However, other forms of gene regulation, for example by lncRNAs, cannot be excluded. LncRNAs are poorly studied in trypanosomatids, with only a single lncRNA characterized to date. Furthermore, it is not clear whether the complete inventory of trypanosomatid lncRNAs is known, because of the inherent cDNA-recoding and DNA-amplification limitations of short-read RNA sequencing. Here, we overcome these limitations by using long-read direct RNA sequencing (DRS) on nanopore arrays. We analyze the native RNA pool of the two main lifecycle stages of the African trypanosome Trypanosoma brucei, with a special emphasis on the inventory of lncRNAs. We identify 207 previously unknown lncRNAs, 32 of which are stage-specifically expressed. We also present insights into the complexity of the T. brucei transcriptome, including alternative transcriptional start and stop sites and potential transcript isoforms, to provide a bias-free understanding of the intricate RNA landscape in T. brucei.

Construction and Application of Orthogonal T7 Expression System in Eukaryote: An Overview

The T7 system is an orthogonal transcription-system, which is characterized by simplicity, higher efficiency, and higher processivity, and it is used for protein or mRNA synthesis in various biological-systems. In comparison with prokaryotes, the construction of the T7 expression system is still on-going in eukaryotes, but it shows greatly applicable prospects. In the present paper, development of T7 expression system construction in eukaryotes is reviewed, including its construction in animal (mammalian cells, trypanosomatid protozoa, Xenopus oocytes, zebrafish), plant, and microorganism and its application in vaccine production and gene therapy. In addition, the innate challenges of T7 expression system construction in eukaryote and its potential application in vaccine production and gene therapy are discussed.

slORFfinder: a tool to detect open reading frames resulting from trans-splicing of spliced leader sequences

Trans-splicing of a spliced leader (SL) to the 5' ends of mRNAs is used to produce mature mRNAs in several phyla of great importance to human health and the marine ecosystem. One of the consequences of the addition of SL sequences is the change or disruption of the open reading frames (ORFs) in the recipient transcripts. Given that most SL sequences have one or more of the trinucleotide NUG, including AUG in flatworms, trans-splicing of SL sequences can potentially supply a start codon to create new ORFs, which we refer to as slORFs, in the recipient mRNAs. Due to the lack of a tool to precisely detect them, slORFs were usually neglected in previous studies. In this work, we present the tool slORFfinder, which automatically links the SL sequences to the recipient mRNAs at the trans-splicing sites identified from SL-containing reads of RNA-Seq and predicts slORFs according to the distribution of ribosome-protected footprints (RPFs) on the trans-spliced transcripts. By applying this tool to the analyses of nematodes, ascidians and euglena, whose RPFs are publicly available, we find wide existence of slORFs in these taxa. Furthermore, we find that slORFs are generally translated at higher levels than the annotated ORFs in the genomes, suggesting they might have important functions. Overall, this study provides a tool, slORFfinder (https://github.com/songbo446/slORFfinder), to identify slORFs, which can enhance our understanding of ORFs in taxa with SL machinery.

Modulation of Virulence Factors during Trypanosoma cruzi Differentiation

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

Clinically relevant benzoxaboroles inhibit mRNA processing in Trypanosoma brucei

OBJECTIVE

The cleavage and polyadenylation endonuclease CPSF73 is thought to be the target of the anti-trypanosomal benzoxaboroles AN7973, acoziborole and AN11736. We previously showed that AN7973 inhibits mRNA processing. We here investigated whether the drug candidates acoziborole (for human sleeping sickness) and AN11736 (for nagana in cattle) have the same effect. We also affinity purified tagged CPSF73 from parasites without, or after, AN7973 treatment, and analysed differentially co-purified proteins by mass spectrometry.

RESULTS

AN11736 and acoziborole both inhibited mRNA processing, as demonstrated by decreased levels of spliced mRNAs and accumulation of di- and tri-cistronic mRNAs from the alpha-beta tubulin locus. Treating the cells with AN7973 for 30 min. did not significantly affect the proteins that copurified with CPSF73.

The endoplasmic reticulum of trypanosomatids: An unrevealed road for chemotherapy

The endoplasmic reticulum (ER) of higher eukaryotic cells forms an intricate membranous network that serves as the main processing facility for folding and assembling of secreted and membrane proteins. The ER is a highly dynamic organelle that interacts with other intracellular structures, as well as endosymbiotic pathogenic and non-pathogenic microorganisms. A strict ER quality control (ERQC) must work to ensure that proteins entering the ER are folded and processed correctly. Unfolded or misfolded proteins are usually identified, selected, and addressed to Endoplasmic Reticulum-Associated Degradation (ERAD) complex. Conversely, when there is a large demand for secreted proteins or ER imbalance, the accumulation of unfolded or misfolded proteins activates the Unfold Protein Response (UPR) to restore the ER homeostasis or, in the case of persistent ER stress, induces the cell death. Pathogenic trypanosomatids, such as Trypanosoma cruzi, Trypanosoma brucei and Leishmania spp are the etiological agents of important neglected diseases. These protozoans have a complex life cycle alternating between vertebrate and invertebrate hosts. The ER of trypanosomatids, like those found in higher eukaryotes, is also specialized for secretion, and depends on the ERAD and non-canonical UPR to deal with the ER stress. Here, we reviewed the basic aspects of ER biology, organization, and quality control in trypanosomatids. We also focused on the unusual way by which T. cruzi, T. brucei, and Leishmania spp. respond to ER stress, emphasizing how these parasites' ER-unrevealed roads might be an attractive target for chemotherapy.

Sequences and proteins that influence mRNA processing in Trypanosoma brucei: Evolutionary conservation of SR-domain and PTB protein functions

BACKGROUND

Spliced leader trans splicing is the addition of a short, capped sequence to the 5' end of mRNAs. It is widespread in eukaryotic evolution, but factors that influence trans splicing acceptor site choice have been little investigated. In Kinetoplastids, all protein-coding mRNAs are 5' trans spliced. A polypyrimidine tract is usually found upstream of the AG splice acceptor, but there is no branch point consensus; moreover, splicing dictates polyadenylation of the preceding mRNA, which is a validated drug target.

METHODOLOGY AND PRINCIPAL FINDINGS

We here describe a trans splicing reporter system that can be used for studies and screens concerning the roles of sequences and proteins in processing site choice and efficiency. Splicing was poor with poly(U) tracts less than 9 nt long, and was influenced by an intergenic region secondary structure. A screen for signals resulted in selection of sequences that were on average 45% U and 35% C. Tethering of either the splicing factor SF1, or the cleavage and polyadenylation factor CPSF3 within the intron stimulated processing in the correct positions, while tethering of two possible homologues of Opisthokont PTB inhibited processing. In contrast, tethering of SR-domain proteins RBSR1, RBSR2, or TSR1 or its interaction partner TSR1IP, promoted use of alternative signals upstream of the tethering sites. RBSR1 interacts predominantly with proteins implicated in splicing, whereas the interactome of RBSR2 is more diverse.

CONCLUSIONS

Our selectable constructs are suitable for screens of both sequences, and proteins that affect mRNA processing in T. brucei. Our results suggest that the functions of PTB and SR-domain proteins in splice site definition may already have been present in the last eukaryotic common ancestor.

The Power of Gene Technologies: 1001 Ways to Create a Cell Model

Modern society faces many biomedical challenges that require urgent solutions. Two of the most important include the elucidation of mechanisms of socially significant diseases and the development of prospective drug treatments for these diseases. Experimental cell models are a convenient tool for addressing many of these problems. The power of cell models is further enhanced when combined with gene technologies, which allows the examination of even more subtle changes within the structure of the genome and permits testing of proteins in a native environment. The list and possibilities of these recently emerging technologies are truly colossal, which requires a rethink of a number of approaches for obtaining experimental cell models. In this review, we analyze the possibilities and limitations of promising gene technologies for obtaining cell models, and also give recommendations on the development and creation of relevant models. In our opinion, this review will be useful for novice cell biologists, as it provides some reference points in the rapidly growing universe of gene and cell technologies.

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.

The Trypanosoma brucei RNA-binding protein DRBD18 ensures correct mRNA trans splicing and polyadenylation patterns

The parasite Trypanosoma brucei grows as bloodstream forms in mammals, and as procyclic forms in tsetse flies. Transcription is polycistronic, all mRNAs are trans spliced, and polyadenylation sites are defined by downstream splicing signals. Expression regulation therefore depends heavily on post-transcriptional mechanisms. The RNA-binding protein DRBD18 was previously implicated in the export of some mRNAs from the nucleus in procyclic forms. It copurifies the outer ring of the nuclear pore, mRNA export factors and exon-junction-complex proteins. We show that for more than 200 mRNAs, DRBD18 depletion caused preferential accumulation of versions with shortened 3'-untranslated regions, arising from use of polyadenylation sites that were either undetectable or rarely seen in nondepleted cells. The shortened mRNAs were often, but not always, more abundant in depleted cells than the corresponding longer versions in normal cells. Their appearance was linked to the appearance of trans-spliced, polyadenylated RNAs containing only downstream 3'-untranslated region-derived sequences. Experiments with one mRNA suggested that nuclear retention alone, through depletion of MEX67, did not affect mRNA length, suggesting a specific effect of DRBD18 on processing. DRBD18-bound mRNAs were enriched in polypyrimidine tract motifs, and DRBD18 was found in both the nucleus and the cytoplasm. We therefore suggest that in the nucleus, DRBD18 might bind to polypyrimidine tracts in 3'-UTRs of mRNA precursors. Such binding might both prevent recognition of mRNA-internal polypyrimidine tracts by splicing factors, and promote export of the processed bound mRNAs to the cytosol.

Finding Correlations Between mRNA and Protein Levels in Leishmania Development: Is There a Discrepancy?

Multiple genes and proteins have been identified as differentially expressed in the stages of the Leishmania life cycle. The differentiation processes are implicated in specific transcriptional and proteomic adjustments driven by gene expression regulation mechanisms. Leishmania parasites lack gene-specific transcriptional control, and gene expression regulation mostly depends on posttranscriptional mechanisms. Due to the lack of transcriptional regulation, criticism regarding the relevance of transcript quantification as a possible and efficient prediction of protein levels is recurrent in studies that use transcriptomic information. The advent of high-throughput technologies has improved the analysis of genomes, transcriptomes and proteomes for different organisms under several conditions. Nevertheless, defining the correlation between transcriptional and proteomic profiles requires arduous and expensive work and remains a challenge in Leishmania. In this review, we analyze transcriptomic and proteomic data for several Leishmania species in two different stages of the parasite life cycle: metacyclogenesis and amastigogenesis (amastigote differentiation). We found a correlation between mRNA and protein levels of 60.9% and 69.8% for metacyclogenesis and amastigogenesis, respectively; showing that majority mRNA and protein levels increase or decrease concomitantly. Among the analyzed genes that did not present correlation indicate that transcriptomic data should be carefully interpreted as protein expression. We also discuss possible explanations and mechanisms involved for this lack of correlation.

Identification and functional implications of pseudouridine RNA modification on small noncoding RNAs in the mammalian pathogen Trypanosoma brucei

Trypanosoma brucei, the parasite that causes sleeping sickness, cycles between an insect and a mammalian host. However, the effect of RNA modifications such as pseudouridinylation on its ability to survive in these two different host environments is unclear. Here, two genome-wide approaches were applied for mapping pseudouridinylation sites (Ψs) on small nucleolar RNA (snoRNA), 7SL RNA, vault RNA, and tRNAs from T. brucei. We show using HydraPsiSeq and RiboMeth-seq that the Ψ on C/D snoRNA guiding 2'-O-methylation increased the efficiency of the guided modification on its target, rRNA. We found differential levels of Ψs on these noncoding RNAs in the two life stages (insect host and mammalian host) of the parasite. Furthermore, tRNA isoform abundance and Ψ modifications were characterized in these two life stages demonstrating stage-specific regulation. We conclude that the differential Ψ modifications identified here may contribute to modulating the function of noncoding RNAs involved in rRNA processing, rRNA modification, protein synthesis, and protein translocation during cycling of the parasite between its two hosts.

Curvicollide D Isolated from the Fungus Amesia sp. Kills African Trypanosomes by Inhibiting Transcription

Sleeping sickness or African trypanosomiasis is a serious health concern with an added socio-economic impact in sub-Saharan Africa due to direct infection in both humans and their domestic livestock. There is no vaccine available against African trypanosomes and its treatment relies only on chemotherapy. Although the current drugs are effective, most of them are far from the modern concept of a drug in terms of toxicity, specificity and therapeutic regime. In a search for new molecules with trypanocidal activity, a high throughput screening of 2000 microbial extracts was performed. Fractionation of one of these extracts, belonging to a culture of the fungus Amesia sp., yielded a new member of the curvicollide family that has been designated as curvicollide D. The new compound showed an inhibitory concentration 50 (IC₅₀) 16-fold lower in Trypanosoma brucei than in human cells. Moreover, it induced cell cycle arrest and disruption of the nucleolar structure. Finally, we showed that curvicollide D binds to DNA and inhibits transcription in African trypanosomes, resulting in cell death. These results constitute the first report on the activity and mode of action of a member of the curvicollide family in T. brucei.

Several different sequences are implicated in bloodstream-form-specific gene expression in Trypanosoma brucei

The parasite Trypanosoma brucei grows as bloodstream forms in mammalian hosts, and as procyclic forms in tsetse flies. In trypanosomes, gene expression regulation depends heavily on post-transcriptional mechanisms. Both the RNA-binding protein RBP10 and glycosomal phosphoglycerate kinase PGKC are expressed only in mammalian-infective forms. RBP10 targets procyclic-specific mRNAs for destruction, while PGKC is required for bloodstream-form glycolysis. Developmental regulation of both is essential: expression of either RBP10 or PGKC in procyclic forms inhibits their proliferation. We show that the 3’-untranslated region of the RBP10 mRNA is extraordinarily long—7.3kb—and were able to identify six different sequences, scattered across the untranslated region, which can independently cause bloodstream-form-specific expression. The 3’-untranslated region of the PGKC mRNA, although much shorter, still contains two different regions, of 125 and 153nt, that independently gave developmental regulation. No short consensus sequences were identified that were enriched either within these regulatory regions, or when compared with other mRNAs with similar regulation, suggesting that more than one regulatory RNA-binding protein is important for repression of mRNAs in procyclic forms. We also identified regions, including an AU repeat, that increased expression in bloodstream forms, or suppressed it in both forms. Trypanosome mRNAs that encode RNA-binding proteins often have extremely extended 3’-untranslated regions. We suggest that one function of this might be to act as a fail-safe mechanism to ensure correct regulation even if mRNA processing or expression of trans regulators is defective.

The parasite Trypanosoma brucei causes sleeping sickness in humans, and nagana in cattle, and is transmitted by Tsetse flies. It grows in the bloodstream and tissue fluids of mammalian hosts, as "bloodstream forms", and as "procyclic forms" in the midgut of tsetse flies. Several hundred proteins are expressed in a stage-specific fashion, and this is essential for parasite survival in the different environments. RBP10 is an RNA-binding protein that is expressed only in bloodstream forms. It binds to procyclic-specific mRNAs, and causes their destruction. PGKC is an enzyme that is also specifically expressed in bloodstream forms. Developmental regulation of both is essential: expression of either RBP10 or PGKC in procyclic forms prevents their growth. The mRNAs encoding both proteins are very unstable in procyclic forms, and the sequences responsible are in an "untranslated region" of the mRNA—sequences that follow the part that codes for protein. We here show that the mRNA encoding PGKC has two regions that independently cause developmental regulation, and that the very long untranslated region of the RBP10 mRNA has no fewer than six regulatory regions, but there were no obvious similarities between them. We suggest that the presence of several different regulatory sequences in trypanosome mRNAs might be a fail-safe mechanism to ensure correct regulation.

Experimental evolution links post-transcriptional regulation to Leishmania fitness gain

The protozoan parasite Leishmania donovani causes fatal human visceral leishmaniasis in absence of treatment. Genome instability has been recognized as a driver in Leishmania fitness gain in response to environmental change or chemotherapy. How genome instability generates beneficial phenotypes despite potential deleterious gene dosage effects is unknown. Here we address this important open question applying experimental evolution and integrative systems approaches on parasites adapting to in vitro culture. Phenotypic analyses of parasites from early and late stages of culture adaptation revealed an important fitness tradeoff, with selection for accelerated growth in promastigote culture (fitness gain) impairing infectivity (fitness costs). Comparative genomics, transcriptomics and proteomics analyses revealed a complex regulatory network associated with parasite fitness gain, with genome instability causing highly reproducible, gene dosage-independent and -dependent changes. Reduction of flagellar transcripts and increase in coding and non-coding RNAs implicated in ribosomal biogenesis and protein translation were not correlated to dosage changes of the corresponding genes, revealing a gene dosage-independent, post-transcriptional mechanism of regulation. In contrast, abundance of gene products implicated in post-transcriptional regulation itself correlated to corresponding gene dosage changes. Thus, RNA abundance during parasite adaptation is controled by direct and indirect gene dosage changes. We correlated differential expression of small nucleolar RNAs (snoRNAs) with changes in rRNA modification, providing first evidence that Leishmania fitness gain in culture may be controlled by post-transcriptional and epitranscriptomic regulation. Our findings propose a novel model for Leishmania fitness gain in culture, where differential regulation of mRNA stability and the generation of modified ribosomes may potentially filter deleterious from beneficial gene dosage effects and provide proteomic robustness to genetically heterogenous, adapting parasite populations. This model challenges the current, genome-centric approach to Leishmania epidemiology and identifies the Leishmania transcriptome and non-coding small RNome as potential novel sources for the discovery of biomarkers that may be associated with parasite phenotypic adaptation in clinical settings.

Genome instability plays a central yet poorly understood role in human disease. Gene amplifications and deletions drive cancer development, microbial infection and therapeutic failure. The molecular mechanisms that harness the deleterious effects of genome instability to generate beneficial phenotypes in pathogenic systems are unknown. Here we study this important open question in the protozoan parasite Leishmania that causes devastating human diseases termed leishmaniases. Leishmania parasites lack transcriptional control and instead exploit genome instability to adapt to their host environment. Analyzing in vitro adaptation of hamster-derived parasites via gene copy number (genomic level) and gene expression changes (transcriptomic and proteomic levels), we show that these parasites likely exploit small nucleolar RNAs (snoRNAs) to mitigate toxic effects of genome instability by post-transcriptional regulation and the establishment of modified ribosomes. Our findings propose non-coding RNAs as potential novel biomarkers with diagnostic and prognostic value that may be linked to changes in parasite tissue tropism or drug susceptibility. This novel insight into Leishmania adaptation will be likely applicable to other fast evolving eukaryotic systems with unstable genomes, such as fungi or cancer cells.

The SPARC complex defines RNAPII promoters in Trypanosoma brucei

Kinetoplastids are a highly divergent lineage of eukaryotes with unusual mechanisms for regulating gene expression. We previously surveyed 65 putative chromatin factors in the kinetoplastid Trypanosoma brucei. Our analyses revealed that the predicted histone methyltransferase SET27 and the Chromodomain protein CRD1 are tightly concentrated at RNAPII transcription start regions (TSRs). Here, we report that SET27 and CRD1, together with four previously uncharacterized constituents, form the SET27 promoter-associated regulatory complex (SPARC), which is specifically enriched at TSRs. SET27 loss leads to aberrant RNAPII recruitment to promoter sites, accumulation of polyadenylated transcripts upstream of normal transcription start sites, and conversion of some normally unidirectional promoters to bidirectional promoters. Transcriptome analysis in the absence of SET27 revealed upregulated mRNA expression in the vicinity of SPARC peaks within the main body of chromosomes in addition to derepression of genes encoding variant surface glycoproteins (VSGs) located in subtelomeric regions. These analyses uncover a novel chromatin-associated complex required to establish accurate promoter position and directionality.

High genome plasticity and frequent genetic exchange in Leishmania tropica isolates from Afghanistan, Iran and Syria

Background

The kinetoplastid protozoan Leishmania tropica mainly causes cutaneous leishmaniasis in humans in the Middle East, and relapse or treatment failure after treatment are common in this area. L. tropica’s digenic life cycle includes distinct stages in the vector sandfly and the mammalian host. Sexual reproduction and genetic exchange appear to occur more frequently than in other Leishmania species. Understanding these processes is complicated by chromosome instability during cell division that yields aneuploidy, recombination and heterozygosity. This combination of rare recombination and aneuploid permits may reveal signs of hypothetical parasexual mating, where diploid cells fuse to form a transient tetraploid that undergoes chromosomal recombination and gradual chromosomal loss.

Methodology/principal findings

The genome-wide SNP diversity from 22 L. tropica isolates showed chromosome-specific runs of patchy heterozygosity and extensive chromosome copy number variation. All these isolates were collected during 2007–2017 in Sweden from patients infected in the Middle East and included isolates from a patient possessing two genetically distinct leishmaniasis infections three years apart with no evidence of re-infection. We found differing ancestries on the same chromosome (chr36) across multiple samples: matching the reference genome with few derived alleles, followed by blocks of heterozygous SNPs, and then by clusters of homozygous SNPs with specific recombination breakpoints at an inferred origin of replication. Other chromosomes had similar marked changes in heterozygosity at strand-switch regions separating polycistronic transcriptional units.

Conclusion/significance

These large-scale intra- and inter-chromosomal changes in diversity driven by recombination and aneuploidy suggest multiple mechanisms of cell reproduction and diversification in L. tropica, including mitotic, meiotic and parasexual processes. It underpins the need for more genomic surveillance of Leishmania, to detect emerging hybrids that could spread more widely and to better understand the association between genetic variation and treatment outcome. Furthering our understanding of Leishmania genome evolution and ancestry will aid better diagnostics and treatment for cutaneous leishmaniasis caused by L.tropica in the Middle East.

Cutaneous leishmaniasis is mainly caused by Leishmania tropica in the Middle East, where it is known for treatment failure and a need for prolonged and/or multiple treatments. Several factors affect the clinical presentation and treatment outcome, such as host genetic variability and specific immune response, as well as environmental factors and the vector species. Little is known about the parasite genome and its influence on treatment response. By analysing the genome of 22 isolates of L. tropica, we have revealed extensive genomic variation and a complex population structure with evidence of genetic exchange within and among the isolates, indicating a possible presence of sexual or parasexual mechanisms. Understanding the Leishmania genome better may improve future treatment and better understanding of treatment failure and relapse.

The Spliced Leader RNA Silencing (SLS) Pathway in Trypanosoma brucei Is Induced by Perturbations of Endoplasmic Reticulum, Golgi Complex, or Mitochondrial Protein Factors: Functional Analysis of SLS-Inducing Kinase PK3

In the parasite Trypanosoma brucei, the causative agent of human African sleeping sickness, all mRNAs are trans-spliced to generate a common 5′ exon derived from the spliced leader (SL) RNA. Perturbations of protein translocation across the endoplasmic reticulum (ER) induce the spliced leader RNA silencing (SLS) pathway. SLS activation is mediated by a serine-threonine kinase, PK3, which translocates from the cytosolic face of the ER to the nucleus, where it phosphorylates the TATA-binding protein TRF4, leading to the shutoff of SL RNA transcription, followed by induction of programmed cell death. Here, we demonstrate that SLS is also induced by depletion of the essential ER-resident chaperones BiP and calreticulin, ER oxidoreductin 1 (ERO1), and the Golgi complex-localized quiescin sulfhydryl oxidase (QSOX). Most strikingly, silencing of Rhomboid-like 1 (TIMRHOM1), involved in mitochondrial protein import, also induces SLS. The PK3 kinase, which integrates SLS signals, is modified by phosphorylation on multiple sites. To determine which of the phosphorylation events activate PK3, several individual mutations or their combination were generated. These mutations failed to completely eliminate the phosphorylation or translocation of the kinase to the nucleus. The structures of PK3 kinase and its ATP binding domain were therefore modeled. A conserved phenylalanine at position 771 was proposed to interact with ATP, and the PK3^F771L mutation completely eliminated phosphorylation under SLS, suggesting that the activation involves most if not all of the phosphorylation sites. The study suggests that the SLS occurs broadly in response to failures in protein sorting, folding, or modification across multiple compartments.

Genome instability drives epistatic adaptation in the human pathogen Leishmania

Chromosome and gene copy number variations often correlate with the evolution of microbial and cancer drug resistance, thus causing important human mortality. How genome instability is harnessed to generate beneficial phenotypes and how deleterious gene dosage effects are compensated remain open questions. The protist pathogen Leishmania exploits genome instability to regulate expression via gene dosage changes. Using these parasites as a unique model system, we uncover complex epistatic interactions between gene copy number variations and compensatory transcriptomic responses as key processes that harness genome instability for adaptive evolution in Leishmania. Our results propose a model of eukaryotic fitness gain that may be broadly applicable to pathogenic fungi or tumor cells known to exploit genome instability for adaptation.

How genome instability is harnessed for fitness gain despite its potential deleterious effects is largely elusive. An ideal system to address this important open question is provided by the protozoan pathogen Leishmania, which exploits frequent variations in chromosome and gene copy number to regulate expression levels. Using ecological genomics and experimental evolution approaches, we provide evidence that Leishmania adaptation relies on epistatic interactions between functionally associated gene copy number variations in pathways driving fitness gain in a given environment. We further uncover posttranscriptional regulation as a key mechanism that compensates for deleterious gene dosage effects and provides phenotypic robustness to genetically heterogenous parasite populations. Finally, we correlate dynamic variations in small nucleolar RNA (snoRNA) gene dosage with changes in ribosomal RNA 2′-O-methylation and pseudouridylation, suggesting translational control as an additional layer of parasite adaptation. Leishmania genome instability is thus harnessed for fitness gain by genome-dependent variations in gene expression and genome-independent compensatory mechanisms. This allows for polyclonal adaptation and maintenance of genetic heterogeneity despite strong selective pressure. The epistatic adaptation described here needs to be considered in Leishmania epidemiology and biomarker discovery and may be relevant to other fast-evolving eukaryotic cells that exploit genome instability for adaptation, such as fungal pathogens or cancer.

The economics of organellar gene loss and endosymbiotic gene transfer

BACKGROUND

The endosymbiosis of the bacterial progenitors of the mitochondrion and the chloroplast are landmark events in the evolution of life on Earth. While both organelles have retained substantial proteomic and biochemical complexity, this complexity is not reflected in the content of their genomes. Instead, the organellar genomes encode fewer than 5% of the genes found in living relatives of their ancestors. While many of the 95% of missing organellar genes have been discarded, others have been transferred to the host nuclear genome through a process known as endosymbiotic gene transfer.

RESULTS

Here, we demonstrate that the difference in the per-cell copy number of the organellar and nuclear genomes presents an energetic incentive to the cell to either delete organellar genes or transfer them to the nuclear genome. We show that, for the majority of transferred organellar genes, the energy saved by nuclear transfer exceeds the costs incurred from importing the encoded protein into the organelle where it can provide its function. Finally, we show that the net energy saved by endosymbiotic gene transfer can constitute an appreciable proportion of total cellular energy budgets and is therefore sufficient to impart a selectable advantage to the cell.

CONCLUSION

Thus, reduced cellular cost and improved energy efficiency likely played a role in the reductive evolution of mitochondrial and chloroplast genomes and the transfer of organellar genes to the nuclear genome.

A distinct complex of PRP19-related and trypanosomatid-specific proteins is required for pre-mRNA splicing in trypanosomes

The pre-mRNA splicing factor PRP19 is recruited into the spliceosome after forming the PRP19/CDC5L complex in humans and the Nineteen complex in yeast. Additionally, ‘PRP19-related’ proteins enter the spliceosome individually or in pre-assemblies that differ in these systems. The protistan family Trypanosomatidae, which harbors parasites such as Trypanosoma brucei, diverged early during evolution from opisthokonts. While introns are rare in these organisms, spliced leader trans splicing is an obligatory step in mRNA maturation. So far, ∼70 proteins have been identified as homologs of human and yeast splicing factors. Moreover, few proteins of unknown function have recurrently co-purified with splicing proteins. Here we silenced the gene of one of these proteins, termed PRC5, and found it to be essential for cell viability and pre-mRNA splicing. Purification of PRC5 combined with sucrose gradient sedimentation revealed a complex of PRC5 with a second trypanosomatid-specific protein, PRC3, and PRP19-related proteins SYF1, SYF3 and ISY1, which we named PRP19-related complex (PRC). Importantly, PRC and the previously described PRP19 complex are distinct from each other because PRC, unlike PRP19, co-precipitates U4 snRNA, which indicates that PRC enters the spliceosome prior to PRP19 and uncovers a unique pre-organization of these proteins in trypanosomes.

Protein Lactylation Critically Regulates Energy Metabolism in the Protozoan Parasite Trypanosoma brucei

Lysine lactylation has been recognized as a novel post-translational modification occurring on histones. However, lactylation in non-histone proteins, especially in proteins of early branching organisms, is not well understood. Energy metabolism and the histone repertoire in the early diverging protozoan parasite Trypanosoma brucei, the causative agent of African trypanosomiasis, markedly diverge from those of conventional eukaryotes. Here, we present the first exhaustive proteome-wide investigation of lactylated sites in T. brucei. We identified 387 lysine-lactylated sites in 257 proteins of various cellular localizations and biological functions. Further, we revealed that glucose metabolism critically regulates protein lactylation in T. brucei although the parasite lacks lactate dehydrogenase. However, unlike mammals, increasing the glucose concentration reduced the level of lactate, and protein lactylation decreased in T. brucei via a unique lactate production pathway. In addition to providing a valuable resource, these foregoing data reveal the regulatory roles of protein lactylation of trypanosomes in energy metabolism and gene expression.

Identification of novel proteins and mRNAs differentially bound to the Leishmania Poly(A) Binding Proteins reveals a direct association between PABP1, the RNA-binding protein RBP23 and mRNAs encoding ribosomal proteins

Poly(A) Binding Proteins (PABPs) are major eukaryotic RNA-binding proteins (RBPs) with multiple roles associated with mRNA stability and translation and characterized mainly from multicellular organisms and yeasts. A variable number of PABP homologues are seen in different organisms however the biological reasons for multiple PABPs are generally not well understood. In the unicellular Leishmania, dependent on post-transcriptional mechanisms for the control of its gene expression, three distinct PABPs are found, with yet undefined functional distinctions. Here, using RNA-immunoprecipitation sequencing analysis we show that the Leishmania PABP1 preferentially associates with mRNAs encoding ribosomal proteins, while PABP2 and PABP3 bind to an overlapping set of mRNAs distinct to those enriched in PABP1. Immunoprecipitation studies combined to mass-spectrometry analysis identified RBPs differentially associated with PABP1 or PABP2, including RBP23 and DRBD2, respectively, that were investigated further. Both RBP23 and DRBD2 bind directly to the three PABPs in vitro, but reciprocal experiments confirmed preferential co-immunoprecipitation of PABP1, as well as the EIF4E4/EIF4G3 based translation initiation complex, with RBP23. Other RBP23 binding partners also imply a direct role in translation. DRBD2, in contrast, co-immunoprecipitated with PABP2, PABP3 and with RBPs unrelated to translation. Over 90% of the RBP23-bound mRNAs code for ribosomal proteins, mainly absent from the transcripts co-precipitated with DRBD2. These experiments suggest a novel and specific route for translation of the ribosomal protein mRNAs, mediated by RBP23, PABP1 and the associated EIF4E4/EIF4G3 complex. They also highlight the unique roles that different PABP homologues may have in eukaryotic cells associated with mRNA translation.

A nuclear enterprise: zooming in on nuclear organization and gene expression control in the African trypanosome

African trypanosomes are early divergent protozoan parasites responsible for high mortality and morbidity as well as a great economic burden among the world's poorest populations. Trypanosomes undergo antigenic variation in their mammalian hosts, a highly sophisticated immune evasion mechanism. Their nuclear organization and mechanisms for gene expression control present several conventional features but also a number of striking differences to the mammalian counterparts. Some of these unorthodox characteristics, such as lack of controlled transcription initiation or enhancer sequences, render their monogenic antigen transcription, which is critical for successful antigenic variation, even more enigmatic. Recent technological developments have advanced our understanding of nuclear organization and gene expression control in trypanosomes, opening novel research avenues. This review is focused on Trypanosoma brucei nuclear organization and how it impacts gene expression, with an emphasis on antigen expression. It highlights several dedicated sub-nuclear bodies that compartmentalize specific functions, whilst outlining similarities and differences to more complex eukaryotes. Notably, understanding the mechanisms underpinning antigen as well as general gene expression control is of great importance, as it might help designing effective control strategies against these organisms.

Comparative analysis of the transcriptional responses of five Leishmania species to trivalent antimony

BACKGROUND

Leishmaniasis is a neglected tropical disease caused by several species of Leishmania. The resistance phenotype of these parasites depends on the characteristics of each species, which contributes to increased therapeutic failures. Understanding the mechanism used by the parasite to survive under treatment pressure in order to identify potential common and specific therapeutic targets is essential for the control of leishmaniasis. The aim of this study was to investigate the expression profiles and potential shared and specific resistance markers of the main Leishmania species of medical importance [subgenus L. (Leishmania): L. donovani, L. infantum and L. amazonensis; subgenus L. (Viannia): L. panamensis and L. braziliensis)] resistant and sensitive to trivalent stibogluconate (SbIII).

METHODS

We conducted comparative analysis of the transcriptomic profiles (only coding sequences) of lines with experimentally induced resistance to SbIII from biological replicates of five Leishmania species available in the databases of four articles based on ortholog attribution. Simultaneously, we carried out functional analysis of ontology and reconstruction of metabolic pathways of the resulting differentially expressed genes (DEGs).

RESULTS

Resistant lines for each species had differential responses in metabolic processes, compound binding, and membrane components concerning their sensitive counterpart. One hundred and thirty-nine metabolic pathways were found, with the three main pathways comprising cysteine and methionine metabolism, glycolysis, and the ribosome. Differentially expressed orthologous genes assigned to species-specific responses predominated, with 899 self-genes. No differentially expressed genes were found in common among the five species. Two common upregulated orthologous genes were found among four species (L. donovani, L. braziliensis, L. amazonensis, and L. panamensis) related to an RNA-binding protein and the NAD(P)H cytochrome-B5-oxidoreductase complex, associated with transcriptional control and de novo synthesis of linoleic acid, critical mechanisms in resistance to antimonials.

CONCLUSION

Herein, we identified potential species-specific genes related to resistance to SbIII. Therefore, we suggest that future studies consider a treatment scheme that is species-specific. Despite the limitations of our study, this is the first approach toward unraveling the pan-genus genetic mechanisms of resistance in leishmaniasis.

The enigmatic RNase MRP of kinetoplastids

The ribonucleoprotein RNase MRP is responsible for the processing of ribosomal RNA precursors. It is found in virtually all eukaryotes that have been examined. In the Euglenozoa, including the genera Euglena, Diplonema and kinetoplastids, MRP RNA and protein subunits have so far escaped detection using bioinformatic methods. However, we now demonstrate that the RNA component is widespread among the Euglenozoa and that these RNAs have secondary structures that conform to the structure of all other phylogenetic groups. In Euglena, we identified the same set of P/MRP protein subunits as in many other protists. However, we failed to identify any of these proteins in the kinetoplastids. This finding poses interesting questions regarding the structure and function of RNase MRP in these species.

Cell-to-Cell Heterogeneity in Trypanosomes

Recent developments in single-cell and single-molecule techniques have revealed surprising levels of heterogeneity among isogenic cells. These advances have transformed the study of cell-to-cell heterogeneity into a major area of biomedical research, revealing that it can confer essential advantages, such as priming populations of unicellular organisms for future environmental stresses. Protozoan parasites, such as trypanosomes, face multiple and often hostile environments, and to survive, they undergo multiple changes, including changes in morphology, gene expression, and metabolism. But why does only a subset of proliferative cells differentiate to the next life cycle stage? Why do only some bloodstream parasites undergo antigenic switching while others stably express one variant surface glycoprotein? And why do some parasites invade an organ while others remain in the bloodstream? Building on extensive research performed in bacteria, here we suggest that biological noise can contribute to the fitness of eukaryotic pathogens and discuss the importance of cell-to-cell heterogeneity in trypanosome infections. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Pseudouridines on Trypanosoma brucei mRNAs are developmentally regulated: Implications to mRNA stability and protein binding

The parasite Trypanosoma brucei cycles between an insect and a mammalian host and is the causative agent of sleeping sickness. Here, we performed high-throughput mapping of pseudouridines (Ψs) on mRNA from two life stages of the parasite. The analysis revealed ~273 Ψs, including developmentally regulated Ψs that are guided by homologs of pseudouridine synthases (PUS1, 3, 5, and 7). Mutating the U that undergoes pseudouridylation in the 3' UTR of valyl-tRNA synthetase destabilized the mRNA level. To investigate the mechanism by which Ψ affects the stability of this mRNA, proteins that bind to the 3' UTR were identified, including the RNA binding protein RBSR1. The binding of RBSR1 protein to the 3' UTR was stronger when lacking Ψ compared to transcripts carrying the modification, suggesting that Ψ can inhibit the binding of proteins to their target and thus affect the stability of mRNAs. Consequently, Ψ modification on mRNA adds an additional level of regulation to the dominant post-transcriptional control in these parasites.

Selective nuclear export of mRNAs is promoted by DRBD18 in Trypanosoma brucei

Kinetoplastids, including Trypanosoma brucei, control gene expression primarily at the posttranscriptional level. Nuclear mRNA export is an important, but understudied, step in this process. The general heterodimeric export factors, Mex67/Mtr2, function in the export of mRNAs and tRNAs in T. brucei, but RNA binding proteins (RBPs) that regulate export processes by controlling the dynamics of Mex67/Mtr2 ribonucleoprotein formation or transport have not been identified. Here, we report that DRBD18, an essential and abundant T. brucei RBP, associates with Mex67/Mtr2 in vivo, likely through its direct interaction with Mtr2. DRBD18 downregulation results in partial accumulation of poly(A)⁺ mRNA in the nucleus, but has no effect on the localization of intron-containing or mature tRNAs. Comprehensive analysis of transcriptomes from whole-cell and cytosol in DRBD18 knockdown parasites demonstrates that depletion of DRBD18 leads to impairment of nuclear export of a subset of mRNAs. CLIP experiments reveal the association of DRBD18 with several of these mRNAs. Moreover, DRBD18 knockdown leads to a partial accumulation of the Mex67/Mtr2 export receptors in the nucleus. Taken together, the current study supports a model in which DRBD18 regulates the selective nuclear export of mRNAs by promoting the mobilization of export competent mRNPs to the cytosol through the nuclear pore complex.

An RNA-binding protein complex regulates the purine-dependent expression of a nucleobase transporter in trypanosomes

Post-transcriptional regulation of gene expression is particularly important in trypanosomatid protozoa. RNA-binding proteins (RBPs) regulate mRNA stability and translation, yet information about how RBPs are able to link environmental cues to post-transcriptional control is scarce. In Trypanosoma brucei, we have previously characterized a short RNA stem-loop cis-element (PuRE, Purine Responsive Element) within the 3'-UTR of the NT8 nucleobase transporter mRNA that is necessary and sufficient to confer a strong repression of gene expression in response to purines. In this study, we have identified a protein complex composed of two RNA-binding proteins (PuREBP1 and PuREBP2) that binds to the PuRE in vitro and to NT8 mRNA in vivo. Depletion of PuREBP1 by RNA interference results in the upregulation of just NT8 and the mRNAs encoding the amino acid transporter AATP6 paralogues. Moreover, we found that the PuREBP1/2 complex is associated with only a handful of mRNAs, and that it is responsible for the observed purine-dependent regulation of NT8 expression.

Computational study on the allosteric mechanism of Leishmania major IF4E-1 by 4E-interacting protein-1: Unravelling the determinants of m7GTP cap recognition

•

Atomistic details on perturbations induced by Lm4E-IP1 binding are described.

•

The modulation of LmIF4E-1 affinity for the cap is confirmed by energetic analyses.

•

Signaling paths between the allosteric and othosteric sites of LmIF4E-1 are predicted.

•

Lm4E-IP1 binding increases the side-chain entropy of W83 and R172 of LmIF4E-1.

•

A mechanism of dynamic allostery is proposed for the regulation mediated by Lm4E-IP1.

During their life cycle, Leishmania parasites display a fine-tuned regulation of the mRNA translation through the differential expression of isoforms of eukaryotic translation initiation factor 4E (LeishIF4Es). The interaction between allosteric modulators such as 4E-interacting proteins (4E-IPs) and LeishIF4E affects the affinity of this initiation factor for the mRNA cap. Here, several computational approaches were employed to elucidate the molecular bases of the previously-reported allosteric modulation in L. major exerted by 4E-IP1 (Lm4E-IP1) on eukaryotic translation initiation factor 4E 1 (LmIF4E-1). Molecular dynamics (MD) simulations and accurate binding free energy calculations (ΔG_bind) were combined with network-based modeling of residue-residue correlations. We also describe the differences in internal motions of LmIF4E-1 apo form, cap-bound, and Lm4E-IP1-bound systems. Through community network calculations, the differences in the allosteric pathways of allosterically-inhibited and active forms of LmIF4E-1 were revealed. The ΔG_bind values show significant differences between the active and inhibited systems, which are in agreement with the available experimental data. Our study thoroughly describes the dynamical perturbations of LmIF4E-1 cap-binding site triggered by Lm4E-IP1. These findings are not only essential for the understanding of a critical process of trypanosomatids’ gene expression but also for gaining insight into the allostery of eukaryotic IF4Es, which could be useful for structure-based design of drugs against this protein family.

Diverse telomeres in trypanosomatids

Telomeres are the ends of linear eukaryotic chromosomes facilitating the resolution of the ‘end replication and protection’ problems, associated with linearity. At the nucleotide level, telomeres typically represent stretches of tandemly arranged telomeric repeats, which vary in length and sequence among different groups of organisms. Recently, a composition of the telomere-associated protein complex has been scrutinized in Trypanosoma brucei. In this work, we subjected proteins from that list to a more detailed bioinformatic analysis and delineated a core set of 20 conserved proteins putatively associated with telomeres in trypanosomatids. Out of these, two proteins (Ku70 and Ku80) are conspicuously missing in representatives of the genus Blastocrithidia, yet telomeres in these species do not appear to be affected. In this work, based on the analysis of a large set of trypanosomatids widely different in their phylogenetic position and life strategies, we demonstrated that telomeres of trypanosomatids are diverse in length, even within groups of closely related species. Our analysis showed that the expression of two proteins predicted to be associated with telomeres (those encoding telomerase and telomere-associated hypothetical protein orthologous to Tb927.6.4330) may directly affect and account for the differences in telomere length within the species of the Leishmania mexicana complex.

Nuclear mRNA maturation and mRNA export control: from trypanosomes to opisthokonts

The passage of mRNAs through the nuclear pores into the cytoplasm is essential in all eukaryotes. For regulation, mRNA export is tightly connected to the full machinery of nuclear mRNA processing, starting at transcription. Export competence of pre-mRNAs gradually increases by both transient and permanent interactions with multiple RNA processing and export factors. mRNA export is best understood in opisthokonts, with limited knowledge in plants and protozoa. Here, I review and compare nuclear mRNA processing and export between opisthokonts and Trypanosoma brucei. The parasite has many unusual features in nuclear mRNA processing, such as polycistronic transcription and trans-splicing. It lacks several nuclear complexes and nuclear-pore-associated proteins that in opisthokonts play major roles in mRNA export. As a consequence, trypanosome mRNA export control is not tight and export can even start co-transcriptionally. Whether trypanosomes regulate mRNA export at all, or whether leakage of immature mRNA to the cytoplasm is kept to a low level by a fast kinetics of mRNA processing remains to be investigated. mRNA export had to be present in the last common ancestor of eukaryotes. Trypanosomes are evolutionary very distant from opisthokonts and a comparison helps understanding the evolution of mRNA export.

Capsid Structure of Leishmania RNA Virus 1

Twelve million people worldwide suffer from leishmaniasis, resulting in more than 30 thousand deaths annually. The disease has several variants that differ in their symptoms.

Leishmania parasites cause a variety of symptoms, including mucocutaneous leishmaniasis, which results in the destruction of the mucous membranes of the nose, mouth, and throat. The species of Leishmania carrying Leishmania RNA virus 1 (LRV1), from the family Totiviridae, are more likely to cause severe disease and are less sensitive to treatment than those that do not contain the virus. Although the importance of LRV1 for the severity of leishmaniasis was discovered a long time ago, the structure of the virus remained unknown. Here, we present a cryo-electron microscopy reconstruction of the virus-like particle of LRV1 determined to a resolution of 3.65 Å. The capsid has icosahedral symmetry and is formed by 120 copies of a capsid protein assembled in asymmetric dimers. RNA genomes of viruses from the family Totiviridae are synthetized, but not capped at the 5′ end, by virus RNA polymerases. To protect viral RNAs from degradation, capsid proteins of the L-A totivirus cleave the 5′ caps of host mRNAs, creating decoys to overload the cellular RNA quality control system. Capsid proteins of LRV1 form positively charged clefts, which may be the cleavage sites for the 5′ cap of Leishmania mRNAs. The putative RNA binding site of LRV1 is distinct from that of the related L-A virus. The structure of the LRV1 capsid enables the rational design of compounds targeting the putative decapping site. Such inhibitors may be developed into a treatment for mucocutaneous leishmaniasis caused by LRV1-positive species of Leishmania.

IMPORTANCE Twelve million people worldwide suffer from leishmaniasis, resulting in more than 30 thousand deaths annually. The disease has several variants that differ in their symptoms. The mucocutaneous form, which leads to disintegration of the nasal septum, lips, and palate, is caused predominantly by Leishmania parasites carrying Leishmania RNA virus 1 (LRV1). Here, we present the structure of the LRV1 capsid determined using cryo-electron microscopy. Capsid proteins of a related totivirus, L-A virus, protect viral RNAs from degradation by cleaving the 5′ caps of host mRNAs. Capsid proteins of LRV1 may have the same function. We show that the LRV1 capsid contains positively charged clefts that may be sites for the cleavage of mRNAs of Leishmania cells. The structure of the LRV1 capsid enables the rational design of compounds targeting the putative mRNA cleavage site. Such inhibitors may be used as treatments for mucocutaneous leishmaniasis.

Spatial integration of transcription and splicing in a dedicated compartment sustains monogenic antigen expression in African trypanosomes

Highly selective gene expression is a key requirement for antigenic variation in several pathogens, allowing evasion of host immune responses and maintenance of persistent infections ¹. African trypanosomes, parasites that cause lethal diseases in humans and livestock, employ an antigenic variation mechanism that involves monogenic antigen expression from a pool of >2600 antigen-coding genes ². In other eukaryotes, the expression of individual genes can be enhanced by mechanisms involving the juxtaposition of otherwise distal chromosomal loci in the three-dimensional nuclear space ^3–5. However, trypanosomes lack classical enhancer sequences or regulated transcription initiation ^6,7. In this context, it has remained unclear how genome architecture contributes to monogenic transcription elongation and transcript processing. Here, we show that the single expressed antigen coding gene displays a specific inter-chromosomal interaction with a major mRNA splicing locus. Chromosome conformation capture (Hi-C) revealed a dynamic reconfiguration of this inter-chromosomal interaction upon activation of another antigen. Super-resolution microscopy showed the interaction to be heritable and splicing dependent. We find a specific association of the two genomic loci with the antigen exclusion complex, whereby VEX1 occupied the splicing locus and VEX2 the antigen coding locus. Following VEX2 depletion, loss of monogenic antigen expression was accompanied by increased interactions between previously silent antigen genes and the splicing locus. Our results reveal a mechanism to ensure monogenic expression, where antigen transcription and mRNA splicing occur in a specific nuclear compartment. These findings suggest a new means of post-transcriptional gene regulation.

Nucleosome landscape reflects phenotypic differences in Trypanosoma cruzi life forms

Trypanosoma cruzi alternates between replicative and nonreplicative life forms, accompanied by a shift in global transcription levels and by changes in the nuclear architecture, the chromatin proteome and histone posttranslational modifications. To gain further insights into the epigenetic regulation that accompanies life form changes, we performed genome-wide high-resolution nucleosome mapping using two T. cruzi life forms (epimastigotes and cellular trypomastigotes). By combining a powerful pipeline that allowed us to faithfully compare nucleosome positioning and occupancy, more than 125 thousand nucleosomes were mapped, and approximately 20% of them differed between replicative and nonreplicative forms. The nonreplicative forms have less dynamic nucleosomes, possibly reflecting their lower global transcription levels and DNA replication arrest. However, dynamic nucleosomes are enriched at nonreplicative regulatory transcription initiation regions and at multigenic family members, which are associated with infective-stage and virulence factors. Strikingly, dynamic nucleosome regions are associated with GO terms related to nuclear division, translation, gene regulation and metabolism and, notably, associated with transcripts with different expression levels among life forms. Finally, the nucleosome landscape reflects the steady-state transcription expression: more abundant genes have a more deeply nucleosome-depleted region at putative 5' splice sites, likely associated with trans-splicing efficiency. Taken together, our results indicate that chromatin architecture, defined primarily by nucleosome positioning and occupancy, reflects the phenotypic differences found among T. cruzi life forms despite the lack of a canonical transcriptional control context.

Structural Differences in Translation Initiation between Pathogenic Trypanosomatids and Their Mammalian Hosts

Canonical mRNA translation in eukaryotes begins with the formation of the 43S pre-initiation complex (PIC). Its assembly requires binding of initiator Met-tRNA_i^Met and several eukaryotic initiation factors (eIFs) to the small ribosomal subunit (40S). Compared to their mammalian hosts, trypanosomatids present significant structural differences in their 40S, suggesting substantial variability in translation initiation. Here, we determine the structure of the 43S PIC from Trypanosoma cruzi, the parasite causing Chagas disease. Our structure shows numerous specific features, such as the variant eIF3 structure and its unique interactions with the large rRNA expansion segments (ESs) 9^S, 7^S, and 6^S, and the association of a kinetoplastid-specific DDX60-like helicase. It also reveals the 40S-binding site of the eIF5 C-terminal domain and structures of key terminal tails of several conserved eIFs underlying their activities within the PIC. Our results are corroborated by glutathione S-transferase (GST) pull-down assays in both human and T. cruzi and mass spectrometry data.

•

Structure of the 43S pre-initiation complex from Trypanosoma cruzi is solved at 3.33 Å

•

The kinetoplastids’ eIF3 core is a septamer that binds mainly the unique, extended ES7^s

•

A kinetoplastid-specific DDX60-like helicase binds to the 43S PIC entry pore

•

The 40S positions of eIF5-CTD and key tails of several eIFs are determined

Developmentally Regulated Novel Non-coding Anti-sense Regulators of mRNA Translation in Trypanosoma b rucei

The parasite Trypanosoma brucei is the causative agent of sleeping sickness and cycles between insect and mammalian hosts. The parasite appears to lack conventional transcriptional regulation of protein coding genes, and mRNAs are processed from polycistronic transcripts by the concerted action of trans-splicing and polyadenylation. Regulation of mRNA function is mediated mainly by RNA binding proteins affecting mRNA stability and translation. In this study, we describe the identification of 62 non-coding (nc) RNAs that are developmentally regulated and/or respond to stress. We characterized two novel anti-sense RNA regulators (TBsRNA-33 and 37) that originate from the rRNA loci, associate with ribosomes and polyribosomes, and interact in vivo with distinct mRNA species to regulate translation. Thus, this study suggests for the first-time anti-sense RNA regulators as an additional layer for controlling gene expression in these parasites.

•

Trypanosome non-coding RNAs (ncRNAs) are developmentally regulated during cycling between two hosts

•

ncRNAs originate from rRNA locus and associate with the ribosome en route to cytoplasm

•

In vivo cross-linking enable identification of target RNA species regulated by ncRNAs

•

Trypanosomes possess anti-sense ncRNAs that regulate translation

Veterinary trypanocidal benzoxaboroles are peptidase-activated prodrugs

Livestock diseases caused by Trypanosoma congolense, T. vivax and T. brucei, collectively known as nagana, are responsible for billions of dollars in lost food production annually. There is an urgent need for novel therapeutics. Encouragingly, promising antitrypanosomal benzoxaboroles are under veterinary development. Here, we show that the most efficacious subclass of these compounds are prodrugs activated by trypanosome serine carboxypeptidases (CBPs). Drug-resistance to a development candidate, AN11736, emerged readily in T. brucei, due to partial deletion within the locus containing three tandem copies of the CBP genes. T. congolense parasites, which possess a larger array of related CBPs, also developed resistance to AN11736 through deletion within the locus. A genome-scale screen in T. brucei confirmed CBP loss-of-function as the primary mechanism of resistance and CRISPR-Cas9 editing proved that partial deletion within the locus was sufficient to confer resistance. CBP re-expression in either T. brucei or T. congolense AN11736-resistant lines restored drug-susceptibility. CBPs act by cleaving the benzoxaborole AN11736 to a carboxylic acid derivative, revealing a prodrug activation mechanism. Loss of CBP activity results in massive reduction in net uptake of AN11736, indicating that entry is facilitated by the concentration gradient created by prodrug metabolism.

Identification of Novel Interspersed DNA Repetitive Elements in the Trypanosoma cruzi Genome Associated with the 3'UTRs of Surface Multigenic Families

Trypanosoma cruzi is the etiological agent of Chagas disease, which affects millions of people in Latin America. No transcriptional control of gene expression has been demonstrated in this organism, and 50% of its genome consists of repetitive elements and members of multigenic families. In this study, we applied a novel bioinformatics approach to predict new repetitive elements in the genome sequence of T. cruzi. A new repetitive sequence measuring 241 nt was identified and found to be interspersed along the genome sequence from strains of different DTUs. This new repeat was mostly on intergenic regions, and upstream and downstream regions of the 241 nt repeat were enriched in surface protein genes. RNAseq analysis revealed that the repeat was part of processed mRNAs and was predominantly found in the 3′ untranslated regions (UTRs) of genes of multigenic families encoding surface proteins. Moreover, we detected a correlation between the presence of the repeat in the 3′UTR of multigenic family genes and the level of differential expression of these genes when comparing epimastigote and trypomastigote transcriptomes. These data suggest that this sequence plays a role in the posttranscriptional regulation of the expression of multigenic families.

A newly identified Leishmania IF4E-interacting protein, Leish4E-IP2, modulates the activity of cap-binding protein paralogs

Translation of most cellular mRNAs in eukaryotes proceeds through a cap-dependent pathway, whereby the cap-binding complex, eIF4F, anchors the preinitiation complex at the 5′ end of mRNAs and regulates translation initiation. The requirement of Leishmania to survive in changing environments can explain why they encode multiple eIF4E (LeishIF4Es) and eIF4G (LeishIF4Gs) paralogs, as each could be assigned a discrete role during their life cycle. Here we show that the expression and activity of different LeishIF4Es change during the growth of cultured promastigotes, urging a search for regulatory proteins. We describe a novel LeishIF4E-interacting protein, Leish4E-IP2, which contains a conserved Y(X)₄LΦ IF4E-binding-motif. Despite its capacity to bind several LeishIF4Es, Leish4E-IP2 was not detected in m⁷GTP-eluted cap-binding complexes, suggesting that it could inhibit the cap-binding activity of LeishIF4Es. Using a functional assay, we show that a recombinant form of Leish4E-IP2 inhibits the cap-binding activity of LeishIF4E-1 and LeishIF4E-3. Furthermore, we show that transgenic parasites expressing a tagged version of Leish4E-IP2 also display reduced cap-binding activities of tested LeishIF4Es, and decreased global translation. Given its ability to bind more than a single LeishIF4E, we suggest that Leish4E-IP2 could serve as a broad-range repressor of Leishmania protein synthesis.

Proteomic analysis of Trypanosoma cruzi spliceosome complex

The unicellular protists of the group Kinetoplastida include the genera Leishmania and Trypanosoma, which are pathogens of invertebrate and vertebrate animals. Despite their medical and economical importance, critical aspects of their biology such as specific molecular characteristics of gene expression regulation are just beginning to be deciphered. Gene expression regulation also depends on post-transcriptional processing steps, such as the trans-splicing process. Despite being widely used in trypanosomes, trans-splicing is a rare event in other eukaryotes. We sought to describe the protein composition of spliceosomes in epimastigotes of T. cruzi, the etiological agent of Chagas disease. We used two TAP-tagged proteins to affinity purify spliceosomes and analyzed their composition by mass spectrometry. Among the 115 identified proteins we detected conserved spliceosome components, as Sm and LSm proteins, RNA helicases, U2- and U5-snRNP specific proteins. Importantly, by comparing our data with proteomic data of human and T. brucei spliceosome complexes, we observed a core group of proteins common to all spliceosomes. By using amino acid sequence comparisons, we identified RNA-associated proteins that might be involved with splicing regulation in T. cruzi, namely the orthologous of WDR33, PABPCL1 and three different HNRNPs. Data are available via ProteomeXchange with identifier PXD018776.

Inducible aging in Hydra oligactis implicates sexual reproduction, loss of stem cells, and genome maintenance as major pathways

Freshwater polyps of the genus Hydra do not age. However, temperature stress induces aging and a shift from reproduction by asexual budding to sexual gamete production in a cold-sensitive (CS) strain of H. oligactis. We sequenced the transcriptome of a male CS strain before and after this life history shift and compared changes in gene expression relative to those seen in a cold-resistant (CR) strain that does not undergo a life history shift in response to altered temperature. We found that the switch from non-aging asexual reproduction to aging and sexual reproduction involves upregulation of genes not only involved in gametogenesis but also genes involved in cellular senescence, apoptosis, and DNA repair accompanied by a downregulation of genes involved in stem cell maintenance. These results suggest that aging is a byproduct of sexual reproduction-associated cellular reprogramming and underscore the power of these H. oligactis strains to identify intrinsic mechanisms of aging.