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 drives epistatic adaptation in the human pathogen Leishmania.

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.

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.

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.

Small nucleolar RNAs controlling rRNA processing in Trypanosoma brucei.

In trypanosomes, in contrast to most eukaryotes, the large subunit (LSU) ribosomal RNA is fragmented into two large and four small ribosomal RNAs (srRNAs) pieces, and this additional processing likely requires trypanosome-specific factors. Here, we examined the role of 10 abundant small nucleolar RNAs (snoRNAs) involved in rRNA processing. We show that each snoRNA involved in LSU processing associates with factors engaged in either early or late biogenesis steps. Five of these snoRNAs interact with the intervening sequences of rRNA precursor, whereas the others only guide rRNA modifications. The function of the snoRNAs was explored by silencing snoRNAs. The data suggest that the LSU rRNA processing events do not correspond to the order of rRNA transcription, and that srRNAs 2, 4 and 6 which are part of LSU are processed before srRNA1. Interestingly, the 6 snoRNAs that affect srRNA1 processing guide modifications on rRNA positions that span locations from the protein exit tunnel to the srRNA1, suggesting that these modifications may serve as check-points preceding the liberation of srRNA1. This study identifies the highest number of snoRNAs so far described that are involved in rRNA processing and/or rRNA folding and highlights their function in the unique trypanosome rRNA maturation events.

Pseudouridines on Trypanosoma brucei spliceosomal small nuclear RNAs and their implication for RNA and protein interactions.

The parasite Trypanosoma brucei, the causative agent of sleeping sickness, cycles between an insect and a mammalian host. Here, we investigated the presence of pseudouridines (Ψs) on the spliceosomal small nuclear RNAs (snRNAs), which may enable growth at the very different temperatures characterizing the two hosts. To this end, we performed the first high-throughput mapping of spliceosomal snRNA Ψs by small RNA Ψ-seq. The analysis revealed 42 Ψs on T. brucei snRNAs, which is the highest number reported so far. We show that a trypanosome protein analogous to human protein WDR79, is essential for guiding Ψ on snRNAs but not on rRNAs. snoRNA species implicated in snRNA pseudouridylation were identified by a genome-wide approach based on ligation of RNAs following in vivo UV cross-linking. snRNA Ψs are guided by single hairpin snoRNAs, also implicated in rRNA modification. Depletion of such guiding snoRNA by RNAi compromised the guided modification on snRNA and reduced parasite growth at elevated temperatures. We further demonstrate that Ψ strengthens U4/U6 RNA-RNA and U2B"/U2A' proteins-U2 snRNA interaction at elevated temperatures. The existence of single hairpin RNAs that modify both the spliceosome and ribosome RNAs is unique for these parasites, and may be related to their ability to cycle between their two hosts that differ in temperature.

The large repertoire of 2'-O-methylation guided by C/D snoRNAs on Trypanosoma brucei rRNA.

The parasite Trypanosoma brucei cycles between insect and mammalian hosts, and is the causative agent of sleeping sickness. Here, we performed genome-wide mapping of 2'-O-methylations (Nms) on trypanosome rRNA using three high-throughput sequencing methods; RibOxi-seq, RiboMeth-seq and 2'-OMe-seq. This is the first study using three genome-wide mapping approaches on rRNA from the same species showing the discrepancy among the methods. RibOxi-seq detects all the sites, but RiboMeth-seq is the only method to evaluate the level of a single Nm site. The sequencing revealed at least ninety-nine Nms guided by eighty-five snoRNAs among these thirty-eight Nms are trypanosome specific sites. We present the sequence and target of the C/D snoRNAs guiding on rRNA. This is the highest number of Nms detected to date on rRNA of a single cell parasite. Based on RiboMeth-seq, several Nm sites were found to be differentially regulated at the two stages of the parasite life cycle, the insect procyclic form (PCF) versus the bloodstream form (BSF) in the mammalian host.

Exosome secretion affects social motility in Trypanosoma brucei.

Extracellular vesicles (EV) secreted by pathogens function in a variety of biological processes. Here, we demonstrate that in the protozoan parasite Trypanosoma brucei, exosome secretion is induced by stress that affects trans-splicing. Following perturbations in biogenesis of spliced leader RNA, which donates its spliced leader (SL) exon to all mRNAs, or after heat-shock, the SL RNA is exported to the cytoplasm and forms distinct granules, which are then secreted by exosomes. The exosomes are formed in multivesicular bodies (MVB) utilizing the endosomal sorting complexes required for transport (ESCRT), through a mechanism similar to microRNA secretion in mammalian cells. Silencing of the ESCRT factor, Vps36, compromised exosome secretion but not the secretion of vesicles derived from nanotubes. The exosomes enter recipient trypanosome cells. Time-lapse microscopy demonstrated that cells secreting exosomes or purified intact exosomes affect social motility (SoMo). This study demonstrates that exosomes are delivered to trypanosome cells and can change their migration. Exosomes are used to transmit stress signals for communication between parasites.

The hnRNP F/H homologue of Trypanosoma brucei is differentially expressed in the two life cycle stages of the parasite and regulates splicing and mRNA stability.

Trypanosomes are protozoan parasites that cycle between a mammalian host (bloodstream form) and an insect host, the Tsetse fly (procyclic stage). In trypanosomes, all mRNAs are trans-spliced as part of their maturation. Genome-wide analysis of trans-splicing indicates the existence of alternative trans-splicing, but little is known regarding RNA-binding proteins that participate in such regulation. In this study, we performed functional analysis of the Trypanosoma brucei heterogeneous nuclear ribonucleoproteins (hnRNP) F/H homologue, a protein known to regulate alternative splicing in metazoa. The hnRNP F/H is highly expressed in the bloodstream form of the parasite, but is also functional in the procyclic form. Transcriptome analyses of RNAi-silenced cells were used to deduce the RNA motif recognized by this protein. A purine rich motif, AAGAA, was enriched in both the regulatory regions flanking the 3' splice site and poly (A) sites of the regulated genes. The motif was further validated using mini-genes carrying wild-type and mutated sequences in the 3' and 5' UTRs, demonstrating the role of hnRNP F/H in mRNA stability and splicing. Biochemical studies confirmed the binding of the protein to this proposed site. The differential expression of the protein and its inverse effects on mRNA level in the two lifecycle stages demonstrate the role of hnRNP F/H in developmental regulation.

Structural and mechanistic insights into the function of Leishmania ribosome lacking a single pseudouridine modification.

Leishmania is the causative agent of cutaneous and visceral diseases affecting millions of individuals worldwide. Pseudouridine (Ψ), the most abundant modification on rRNA, changes during the parasite life cycle. Alterations in the level of a specific Ψ in helix 69 (H69) affected ribosome function. To decipher the molecular mechanism of this phenotype, we determine the structure of ribosomes lacking the single Ψ and its parental strain at ∼2.4-3 Å resolution using cryo-EM. Our findings demonstrate the significance of a single Ψ on H69 to its structure and the importance for its interactions with helix 44 and specific tRNAs. Our study suggests that rRNA modification affects translation of mRNAs carrying codon bias due to selective accommodation of tRNAs by the ribosome. Based on the high-resolution structures, we propose a mechanism explaining how the ribosome selects specific tRNAs.

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.

Lytic Reactivation of the Kaposi's Sarcoma-Associated Herpesvirus (KSHV) Is Accompanied by Major Nucleolar Alterations.

The nucleolus is a subnuclear compartment whose primary function is the biogenesis of ribosomal subunits. Certain viral infections affect the morphology and composition of the nucleolar compartment and influence ribosomal RNA (rRNA) transcription and maturation. However, no description of nucleolar morphology and function during infection with Kaposi's sarcoma-associated herpesvirus (KSHV) is available to date. Using immunofluorescence microscopy, we documented extensive destruction of the nuclear and nucleolar architecture during the lytic reactivation of KSHV. This was manifested by the redistribution of key nucleolar proteins, including the rRNA transcription factor UBF. Distinct delocalization patterns were evident; certain nucleolar proteins remained together whereas others dissociated, implying that nucleolar proteins undergo nonrandom programmed dispersion. Significantly, the redistribution of UBF was dependent on viral DNA replication or late viral gene expression. No significant changes in pre-rRNA levels and no accumulation of pre-rRNA intermediates were found by RT-qPCR and Northern blot analysis. Furthermore, fluorescent in situ hybridization (FISH), combined with immunofluorescence, revealed an overlap between Fibrillarin and internal transcribed spacer 1 (ITS1), which represents the primary product of the pre-rRNA, suggesting that the processing of rRNA proceeds during lytic reactivation. Finally, small changes in the levels of pseudouridylation (Ψ) and 2'-O-methylation (Nm) were documented across the rRNA; however, none were localized to the functional domain. Taken together, our results suggest that despite dramatic changes in the nucleolar organization, rRNA transcription and processing persist during lytic reactivation of KSHV. Whether the observed nucleolar alterations favor productive infection or signify cellular anti-viral responses remains to be determined.

A long noncoding RNA promotes parasite differentiation in African trypanosomes.

The parasite Trypanosoma brucei causes African sleeping sickness that is fatal to patients if untreated. Parasite differentiation from a replicative slender form into a quiescent stumpy form promotes host survival and parasite transmission. Long noncoding RNAs (lncRNAs) are known to regulate cell differentiation in other eukaryotes. To determine whether lncRNAs are also involved in parasite differentiation, we used RNA sequencing to survey the T. brucei genome, identifying 1428 previously uncharacterized lncRNA genes. We find that grumpy lncRNA is a key regulator that promotes parasite differentiation into the quiescent stumpy form. This function is promoted by a small nucleolar RNA encoded within the grumpy lncRNA. snoGRUMPY binds to messenger RNAs of at least two stumpy regulatory genes, promoting their expression. grumpy overexpression reduces parasitemia in infected mice. Our analyses suggest that T. brucei lncRNAs modulate parasite-host interactions and provide a mechanism by which grumpy regulates cell differentiation in trypanosomes.

Dynamic refolding of IFN-gamma mRNA enables it to function as PKR activator and translation template.

Interferon-gamma mRNA activates the RNA-dependent protein kinase PKR, which in turn strongly attenuates translation of interferon-gamma mRNA. Unlike riboswitches restricted to noncoding regions, the interferon-gamma RNA domain that activates PKR comprises the 5' UTR and 26 translated codons. Extensive interferon-gamma coding sequence is thus dedicated to activating PKR and blocking interferon-gamma synthesis. This implies that the PKR activator is disrupted by ribosomes during translation initiation and must refold promptly to restore PKR activation. The activator structure harbors an essential kink-turn, probably to allow formation of a pseudoknot that is critical for PKR activation. Three indispensable short helices, bordered by orientation-sensitive base pairs, align with the pseudoknot stem, generating RNA helix of sufficient length to activate PKR. Through gain-of-function mutations, we show that the RNA activator can adopt alternative conformations that activate PKR. This flexibility promotes efficient refolding of interferon-gamma mRNA, which is necessary for its dual function as translation template and activator of PKR, and which thus prevents overexpression of this inflammatory cytokine.

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.

The effect of reverse transcription enzymes and conditions on high throughput amplicon sequencing of the 16S rRNA.

It is assumed that the sequencing of ribosomes better reflects the active microbial community than the sequencing of the ribosomal RNA encoding genes. Yet, many studies exploring microbial communities in various environments, ranging from the human gut to deep oceans, questioned the validity of this paradigm due to the discrepancies between the DNA and RNA based communities. Here, we focus on an often neglected key step in the analysis, the reverse transcription (RT) reaction. Previous studies showed that RT may introduce biases when expressed genes and ribosmal rRNA are quantified, yet its effect on microbial diversity and community composition was never tested. High throughput sequencing of ribosomal RNA is a valuable tool to understand microbial communities as it better describes the active population than DNA analysis. However, the necessary step of RT may introduce biases that have so far been poorly described. In this manuscript, we compare three RT enzymes, commonly used in soil microbiology, in two temperature modes to determine a potential source of bias due to non-standardized RT conditions. In our comparisons, we have observed up to six fold differences in bacterial class abundance. A temperature induced bias can be partially explained by G-C content of the affected bacterial groups, thus pointing toward a need for higher reaction temperatures. However, another source of bias was due to enzyme processivity differences. This bias is potentially hard to overcome and thus mitigating it might require the use of one enzyme for the sake of cross-study comparison.

Control of mRNA Splicing by Intragenic RNA Activators of Stress Signaling: Potential Implications for Human Disease.

A critical step in the cellular stress response is transient activation of the RNA-dependent protein kinase PKR by double-helical RNA, resulting in down-regulation of protein synthesis through phosphorylation of the α chain of translation initiation factor eIF2, a major PKR substrate. However, intragenic elements of 100-200 nucleotides in length within primary transcripts of cellular genes, exemplified by the tumor necrosis factor (TNF)-α gene and fetal and adult globin genes, are capable of forming RNA structures that potently activate PKR and thereby strongly enhance mRNA splicing efficiency. By inducing nuclear eIF2α phosphorylation, these PKR activator elements enable highly efficient early spliceosome assembly yet do not impair translation of the mature spliced mRNA. The TNF-α RNA activator of PKR folds into a compact pseudoknot that is highly conserved within the phylogeny. Upon excision of ß-globin first intron, the RNA activator of PKR, located in exon 1, is silenced through strand displacement by a short sequence within exon 2, restricting thereby the ability to activate PKR to the splicing process without impeding subsequent synthesis of ß-globin essential for survival. This activator/silencer mechanism likewise controls splicing of α-globin pre-mRNA, but the exonic locations of PKR activator and silencer sequences are reversed, demonstrating evolutionary flexibility. Impaired splicing efficiency may underlie numerous human ß-thalassemia mutations that map to the ß-globin RNA activator of PKR or its silencer. Even where such mutations change the encoded amino acid sequence during subsequent translation, they carry the potential of first impairing PKR-dependent mRNA splicing or shutoff of PKR activation needed for optimal translation.

PKR activation and eIF2α phosphorylation mediate human globin mRNA splicing at spliceosome assembly.

Short elements in mammalian mRNA can control gene expression by activating the RNA-dependent protein kinase PKR that attenuates translation by phosphorylating cytoplasmic eukaryotic initiation factor 2α (eIF2α). We demonstrate a novel, positive role for PKR activation and eIF2α phosphorylation in human globin mRNA splicing. PKR localizes in splicing complexes and associates with splicing factor SC35. Splicing and early-stage spliceosome assembly on ß-globin pre-mRNA depend strictly on activation of PKR by a codon-containing RNA fragment within exon 1 and on phosphorylation of nuclear eIF2α on Serine 51. Nonphosphorylatable mutant eIF2αS51A blocked ß-globin mRNA splicing in cells and nuclear extract. Mutations of the ß-globin RNA activator abrogated PKR activation and profoundly affected mRNA splicing efficiency. PKR depletion abrogated splicing and spliceosome assembly; recombinant PKR effectively restored splicing. Excision of the first intron of ß-globin induces strand displacement within the RNA activator of PKR by a sequence from exon 2, a structural rearrangement that silences the ability of spliced ß-globin mRNA to activate PKR. Thus, the ability to activate PKR is transient, serving solely to enable splicing. α-Globin pre-mRNA splicing is controlled likewise but positions of PKR activator and silencer are reversed, demonstrating evolutionary flexibility in how PKR activation regulates globin mRNA splicing through eIF2α phosphorylation.

The Canonical Poly (A) Polymerase PAP1 Polyadenylates Non-Coding RNAs and Is Essential for snoRNA Biogenesis in Trypanosoma brucei.

The parasite Trypanosoma brucei is the causative agent of African sleeping sickness and is known for its unique RNA processing mechanisms that are common to all the kinetoplastidea including Leishmania and Trypanosoma cruzi. Trypanosomes possess two canonical RNA poly (A) polymerases (PAPs) termed PAP1 and PAP2. PAP1 is encoded by one of the only two genes harboring cis-spliced introns in this organism, and its function is currently unknown. In trypanosomes, all mRNAs, and non-coding RNAs such as small nucleolar RNAs (snoRNAs) and long non-coding RNAs (lncRNAs), undergo trans-splicing and polyadenylation. Here, we show that the function of PAP1, which is located in the nucleus, is to polyadenylate non-coding RNAs, which undergo trans-splicing and polyadenylation. Major substrates of PAP1 are the snoRNAs and lncRNAs. Under the silencing of either PAP1 or PAP2, the level of snoRNAs is reduced. The dual polyadenylation of snoRNA intermediates is carried out by both PAP2 and PAP1 and requires the factors essential for the polyadenylation of mRNAs. The dual polyadenylation of the precursor snoRNAs by PAPs may function to recruit the machinery essential for snoRNA processing.

Thiamine pyrophosphate riboswitches are targets for the antimicrobial compound pyrithiamine.

Thiamine metabolism genes are regulated in numerous bacteria by a riboswitch class that binds the coenzyme thiamine pyrophosphate (TPP). We demonstrate that the antimicrobial action of the thiamine analog pyrithiamine (PT) is mediated by interaction with TPP riboswitches in bacteria and fungi. For example, pyrithiamine pyrophosphate (PTPP) binds the TPP riboswitch controlling the tenA operon in Bacillus subtilis. Expression of a TPP riboswitch-regulated reporter gene is reduced in transgenic B. subtilis or Escherichia coli when grown in the presence of thiamine or PT, while mutant riboswitches in these organisms are unresponsive to these ligands. Bacteria selected for PT resistance bear specific mutations that disrupt ligand binding to TPP riboswitches and derepress certain TPP metabolic genes. Our findings demonstrate that riboswitches can serve as antimicrobial drug targets and expand our understanding of thiamine metabolism in bacteria.