The DEAD-box RNA helicase PfDOZI imposes opposing actions on RNA metabolism in Plasmodium falciparum.

In malaria parasites, the regulation of mRNA translation, storage and degradation during development and life-stage transitions remains largely unknown. Here, we functionally characterized the DEAD-box RNA helicase PfDOZI in P. falciparum. Disruption of pfdozi enhanced asexual proliferation but reduced sexual commitment and impaired gametocyte development. By quantitative transcriptomics, we show that PfDOZI is involved in the regulation of invasion-related genes and sexual stage-specific genes during different developmental stages. PfDOZI predominantly participates in processing body-like mRNPs in schizonts but germ cell granule-like mRNPs in gametocytes to impose opposing actions of degradation and protection on different mRNA targets. We further show the formation of stress granule-like mRNPs during nutritional deprivation, highlighting an essential role of PfDOZI-associated mRNPs in stress response. We demonstrate that PfDOZI participates in distinct mRNPs to maintain mRNA homeostasis in response to life-stage transition and environmental changes by differentially executing post-transcriptional regulation on the target mRNAs.

Transcript tinkering: RNA modifications in protozoan parasites.

Apicomplexan and trypanosomatid parasites have evolved a wide range of post-transcriptional processes that allow them to replicate, differentiate, and transmit within and among multiple different tissue, host, and vector environments. In this review, we highlight the recent advances that point toward the regulatory potential of RNA modifications in mediating these processes on the coding and noncoding transcriptome throughout the life cycle of protozoan parasites. We discuss the recent technical advancements enabling the study of the 'epitranscriptome' and how parasites evolved RNA modification-mediated mechanisms adapted to their unique lifestyles.

A role for a Trypanosoma brucei cytosine RNA methyltransferase homolog in ribosomal RNA processing.

In Trypanosoma brucei, gene expression is primarily regulated posttranscriptionally making RNA metabolism critical. T. brucei has an epitranscriptome containing modified RNA bases. Yet, the identity of the enzymes catalyzing modified RNA base addition and the functions of the enzymes and modifications remain unclear. Homology searches indicate the presence of numerous T. brucei cytosine RNA methyltransferase homologs. One such homolog, TbNop2 was studied in detail. TbNop2 contains the six highly conserved motifs found in cytosine RNA methyltransferases and is evolutionarily related to the Nop2 protein family required for rRNA modification and processing. RNAi experiments targeting TbNop2 resulted in reduced levels of TbNop2 RNA and protein, and a cessation of parasite growth. Next generation sequencing of bisulfite-treated RNA (BS-seq) detected the presence of two methylation sites in the large rRNA; yet TbNop2 RNAi did not result in a significant reduction of methylation. However, TbNop2 RNAi resulted in the retention of 28S internal transcribed spacer RNAs, indicating a role for TbNop2 in rRNA processing.

Unconventional features in the transcription and processing of spliceosomal small nuclear RNAs in the protozoan parasite Trichomonas vaginalis.

Trichomonas vaginalis is a medically important protozoan parasite, and a deep-branching, evolutionarily divergent unicellular eukaryote that has conserved several key features of eukaryotic gene expression. Trichomonas vaginalis possesses a metazoan/plant-like capping apparatus, mRNAs with a cap 1 structure and spliceosomes containing the five small nuclear RNAs (snRNAs). However, in contrast to metazoan and plant snRNAs, the structurally conserved T. vaginalis snRNAs were initially identified as lacking the canonical guanosine cap nucleotide. To explain this unusual condition, we sought to investigate transcriptional and processing features of the spliceosomal snRNAs in this protist. Here, we show that T. vaginalis spliceosomal snRNA genes mostly lack typical eukaryotic promoters. In contrast to other eukaryotes, the putative TATA box in the T. vaginalis U6 snRNA gene was found to be dispensable for transcription or RNA polymerase selectivity. Moreover, U6 transcription in T. vaginalis was virtually insensitive to tagetitoxin compared with other cellular transcripts produced by the same RNA polymerase III. Most important and unexpected, snRNA transcription in T. vaginalis appears to bypass capping as we show that these transcripts retain their original 5'-triphosphate groups. In conclusion, transcription and processing of spliceosomal snRNAs in T. vaginalis deviate considerably from the conventional rules of other eukaryotes.

From a bimodal to a multi-stage view on trypanosomes' differential RNA editing.

Significant variations in the abundance of mitochondrial RNA processing proteins and their target RNAs across trypanosome life stages present an opportunity to explore the regulatory mechanisms that drive these changes. Utilizing omics approaches can uncover unconventional targets, aiding our understanding of the parasites' adaptation and enabling targeted interventions for differentiation.

RNA editing catalytic complexes edit multiple mRNA sites non-processively in Trypanosoma brucei.

RNA editing generates mature mitochondrial mRNAs in T. brucei by extensive uridine insertion and deletion at numerous editing sites (ESs) as specified by guide RNAs (gRNAs). The editing is performed by three RNA Editing Catalytic Complexes (RECCs) which each have a different endonuclease in addition to 12 proteins in common resulting in RECC1 that is specific for deletion ESs and RECC2 and RECC3 that are specific for insertion ESs. Thus, different RECCs are required for editing of mRNA sequence regions where single gRNAs specify a combination of insertion and deletion ESs. We investigated how the three different RECCs might edit combinations of insertion and deletion ESs that are specified by single gRNAs by testing whether their endonuclease compositions are stable or dynamic during editing. We analyzed in vivo BirA* proximity labeling and found that the endonucleases remain associated with their set of common RECC proteins during editing when expressed at normal physiological levels. We also found that overexpression of endonuclease components resulted in minor effects on RECCs but did not affect growth. Thus, the protein stoichiometries that exist within each RECC can be altered by perturbations of RECC expression levels. These results indicate that editing of consecutive insertion and deletion ESs occurs by successive engagement and disengagement of RECCs, i.e., is non-processive, which is likely the case for consecutive pairs of insertion or deletion ESs. This clarifies the nature of the complex patterns of partially edited mRNAs that occur in vivo.

Multiple domains of the integral KREPA3 protein are critical for the structure and precise functions of RNA editing catalytic complexes in Trypanosoma brucei.

The gRNA directed U-insertion and deletion editing of mitochondrial mRNAs that is essential in different life-cycle stages for the protozoan parasite Trypanosoma brucei is performed by three similar multiprotein catalytic complexes (CCs) that contain the requisite enzymes. These CCs also contain a common set of eight proteins that have no apparent direct catalytic function, including six that have an OB-fold domain. We show here that one of these OB-fold proteins, KREPA3 (A3), has structural homology to other editing proteins, is essential for editing, and is multifunctional. We investigated A3 function by analyzing the effects of single amino acid loss of function mutations, most of which were identified by screening bloodstream form (BF) parasites for loss of growth following random mutagenesis. Mutations in the zinc fingers (ZFs), an intrinsically disordered region (IDR), and several within or near the carboxy-terminal OB-fold domain variably impacted CC structural integrity and editing. Some mutations resulted in almost complete loss of CCs and its proteins and editing, whereas others retained CCs but had aberrant editing. All but a mutation which is near the OB-fold affected growth and editing in BF but not procyclic form (PF) parasites. These data indicate that multiple positions within A3 have essential functions that contribute to the structural integrity of CCs, the precision of editing and the developmental differences in editing between BF and PF stages.

Trypanosome RNA helicase KREH2 differentially controls non-canonical editing and putative repressive structure via a novel proposed 'bifunctional' gRNA in mRNA A6.

U-insertion/deletion (U-indel) RNA editing in trypanosome mitochondria is directed by guide RNAs (gRNAs). This editing may developmentally control respiration in bloodstream forms (BSF) and insect procyclic forms (PCF). Holo-editosomes include the accessory RNA Editing Substrate Binding Complex (RESC) and RNA Editing Helicase 2 Complex (REH2C), but the specific proteins controlling differential editing remain unknown. Also, RNA editing appears highly error prone because most U-indels do not match the canonical pattern. However, despite extensive non-canonical editing of unknown functions, accurate canonical editing is required for normal cell growth. In PCF, REH2C controls editing fidelity in RESC-bound mRNAs. Here, we report that KREH2, a REH2C-associated helicase, developmentally controls programmed non-canonical editing, including an abundant 3' element in ATPase subunit 6 (A6) mRNA. The 3' element sequence is directed by a proposed novel regulatory gRNA. In PCF, KREH2 RNAi-knockdown up-regulates the 3' element, which establishes a stable structure hindering element removal by canonical initiator-gRNA-directed editing. In BSF, KREH2-knockdown does not up-regulate the 3' element but reduces its high abundance. Thus, KREH2 differentially controls extensive non-canonical editing and associated RNA structure via a novel regulatory gRNA, potentially hijacking factors as a 'molecular sponge'. Furthermore, this gRNA is bifunctional, serving in canonical CR4 mRNA editing whilst installing a structural element in A6 mRNA.

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

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

Structural basis for guide RNA selection by the RESC1-RESC2 complex.

Kinetoplastid parasites, such as trypanosomes or leishmania, rely on RNA-templated RNA editing to mature mitochondrial cryptic pre-mRNAs into functional protein-coding transcripts. Processive pan-editing of multiple editing blocks within a single transcript is dependent on the 20-subunit RNA editing substrate binding complex (RESC) that serves as a platform to orchestrate the interactions between pre-mRNA, guide RNAs (gRNAs), the catalytic RNA editing complex (RECC), and a set of RNA helicases. Due to the lack of molecular structures and biochemical studies with purified components, neither the spacio-temporal interplay of these factors nor the selection mechanism for the different RNA components is understood. Here we report the cryo-EM structure of Trypanosoma brucei RESC1-RESC2, a central hub module of the RESC complex. The structure reveals that RESC1 and RESC2 form an obligatory domain-swapped dimer. Although the tertiary structures of both subunits closely resemble each other, only RESC2 selectively binds 5'-triphosphate-nucleosides, a defining characteristic of gRNAs. We therefore propose RESC2 as the protective 5'-end binding site for gRNAs within the RESC complex. Overall, our structure provides a starting point for the study of the assembly and function of larger RNA-bound kinetoplast RNA editing modules and might aid in the design of anti-parasite drugs.

KRGG1 function in RNA editing in Trypanosoma brucei.

Mitochondrial gene expression in trypanosomes requires numerous multiprotein complexes that are unique to kinetoplastids. Among these, the most well characterized are RNA editing catalytic complexes (RECCs) that catalyze the guide RNA (gRNA)-specified insertion and deletion of uridines during mitochondrial mRNA maturation. This post-transcriptional resequencing of mitochondrial mRNAs can be extensive, involving dozens of different gRNAs and hundreds of editing sites with most of the mature mRNA sequences resulting from the editing process. Proper coordination of the editing with the cognate gRNAs is attributed to RNA editing substrate-binding complexes (RESCs), which are also required for RNA editing. Although the precise mechanism of RESC function is less well understood, their affinity for binding both editing substrates and products suggests that these complexes may provide a scaffold for RECC catalytic processing. KRGG1 has been shown to bind RNAs, and although affinity purification co-isolates RESC complexes, its role in RNA editing remains uncertain. We show here that KRGG1 is essential in BF parasites and required for normal editing. KRGG1 repression results in reduced amounts of edited A6 mRNA and increased amounts of edited ND8 mRNA. Sequence and structure analysis of KRGG1 identified a region of homology with RESC6, and both proteins have predicted tandem helical repeats that resemble ARM/HEAT motifs. The ARM/HEAT-like region is critical for function as exclusive expression of mutated KRGG1 results in growth inhibition and disruption of KRGG1 association with RESCs. These results indicate that KRGG1 is critical for RNA editing and its specific function is associated with RESC activity.

Conserved and transcript-specific functions of the RESC factors, RESC13 and RESC14, in kinetoplastid RNA editing.

Uridine insertion/deletion RNA editing is an extensive post-transcriptional modification of mitochondrial mRNAs in kinetoplastid organisms, including Trypanosoma brucei This process is carried out using trans-acting gRNAs and complex protein machinery. The essential RNA editing substrate binding complex (RESC) serves as the scaffold that modulates protein and RNA interactions during editing, and contains the guide RNA binding complex (GRBC), the RNA editing mediator complexes (REMCs), and organizer proteins. Despite the importance of RESC in editing, the functions of each protein comprising this complex are not completely understood. Here, we further define the roles of a REMC protein, RESC13, and a RESC organizer, RESC14, using high-throughput sequencing on two large pan-edited mRNAs, A6 and COIII. When comparing our analyses to that of a previously published small pan-edited mRNA, RPS12, we find that RESC13 has conserved functions across the three transcripts with regard to editing initiation, gRNA utilization, gRNA exchange, and restricting the formation of long misedited junctions that likely arise from its ability to modulate RNA structure. However, RESC13 does have transcript-specific effects on the types of long junctions whose formation it restricts. RESC14 has a conserved effect on gRNA utilization across the three transcripts analyzed, but has transcript-specific effects on editing initiation, gRNA exchange, and junction formation. Our data suggest that transcript-specific effects of both proteins are due to differences in transcript length and sequences as well as transcript-specific protein interactions. These findings highlight the importance of studying multiple transcripts to determine the function of editing factors.

Domain function and predicted structure of three heterodimeric endonuclease subunits of RNA editing catalytic complexes in Trypanosoma brucei.

Each of the three similar RNA Editing Catalytic Complexes (RECCs) that perform gRNA-directed uridine insertion and deletion during Trypanosoma brucei mitochondrial (mt) mRNA editing has a distinct endonuclease activity that requires two related RNase III proteins, with only one competent for catalysis. We identified multiple loss-of-function mutations in the RNase III and other motifs of the non-catalytic KREPB6, KREPB7, and KREPB8 components by random mutagenesis and screening. These mutations had various effects on growth, editing, and both the abundances and RECC associations of these RNase III protein pairs in bloodstream form (BF) and procyclic form (PF) cells. Protein structure modelling predicted that the Zinc Finger (ZnF) of each paired RNase III protein contacts RNA positioned at the heterodimeric active site which is flanked by helices of a novel RNase III-Associated Motif (RAM). The results indicate that the protein domains of the non-catalytic subunits function together in RECC integrity, substrate binding, and editing site recognition during the multistep RNA editing process. Additionally, several mutants display distinct functional consequences in different life cycle stages. These results highlight the complementary roles of protein pairs and three RECCs within the complicated T. brucei mRNA editing machinery that matures mt mRNAs differentially between developmental stages.

Fuzzy RNA recognition by the Trypanosoma brucei editosome.

The assembly of high molecular mass ribonucleoprotein complexes typically relies on the binary interaction of defined RNA sequences or precisely folded RNA motifs with dedicated RNA-binding domains on the protein side. Here we describe a new molecular recognition principle of RNA molecules by a high molecular mass protein complex. By chemically probing the solvent accessibility of mitochondrial pre-mRNAs when bound to the Trypanosoma brucei editosome, we identified multiple similar but non-identical RNA motifs as editosome contact sites. However, by treating the different motifs as mathematical graph objects we demonstrate that they fit a consensus 2D-graph consisting of 4 vertices (V) and 3 edges (E) with a Laplacian eigenvalue of 0.5477 (λ2). We establish that synthetic 4V(3E)-RNAs are sufficient to compete for the editosomal pre-mRNA binding site and that they inhibit RNA editing in vitro. Furthermore, we demonstrate that only two topological indices are necessary to predict the binding of any RNA motif to the editosome with a high level of confidence. Our analysis corroborates that the editosome has adapted to the structural multiplicity of the mitochondrial mRNA folding space by recognizing a fuzzy continuum of RNA folds that fit a consensus graph descriptor.

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

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

Transcriptome-Wide Identification of Coding and Noncoding RNA-Binding Proteins Defines the Comprehensive RNA Interactome of Leishmania mexicana.

Proteomic profiling of RNA-binding proteins in Leishmania is currently limited to polyadenylated mRNA-binding proteins, leaving proteins that interact with nonadenylated RNAs, including noncoding RNAs and pre-mRNAs, unidentified. Using a combination of unbiased orthogonal organic phase separation methodology and tandem mass tag-labeling-based high resolution quantitative proteomic mass spectrometry, we robustly identified 2,417 RNA-binding proteins, including 1289 putative novel non-poly(A)-RNA-binding proteins across the two main Leishmania life cycle stages. Eight out of 20 Leishmania deubiquitinases, including the recently characterized L. mexicana DUB2 with an elaborate RNA-binding protein interactome were exclusively identified in the non-poly(A)-RNA-interactome. Additionally, an increased representation of WD40 repeat domains were observed in the Leishmania non-poly(A)-RNA-interactome, thus uncovering potential involvement of this protein domain in RNA-protein interactions in Leishmania. We also characterize the protein-bound RNAs using RNA-sequencing and show that in addition to protein coding transcripts ncRNAs are also enriched in the protein-RNA interactome. Differential gene expression analysis revealed enrichment of 142 out of 195 total L. mexicana protein kinase genes in the protein-RNA-interactome, suggesting important role of protein-RNA interactions in the regulation of the Leishmania protein kinome. Additionally, we characterize the quantitative changes in RNA-protein interactions in hundreds of Leishmania proteins following inhibition of heat shock protein 90 (Hsp90). Our results show that the Hsp90 inhibition in Leishmania causes widespread disruption of RNA-protein interactions in ribosomal proteins, proteasomal proteins and translation factors in both life cycle stages, suggesting downstream effect of the inhibition on protein synthesis and degradation pathways in Leishmania. This study defines the comprehensive RNA interactome of Leishmania and provides in-depth insight into the widespread involvement of RNA-protein interactions in Leishmania biology. IMPORTANCE Advances in proteomics and mass spectrometry have revealed the mRNA-binding proteins in many eukaryotic organisms, including the protozoan parasites Leishmania spp., the causative agents of leishmaniasis, a major infectious disease in over 90 tropical and subtropical countries. However, in addition to mRNAs, which constitute only 2 to 5% of the total transcripts, many types of non-coding RNAs participate in crucial biological processes. In Leishmania, RNA-binding proteins serve as primary gene regulators. Therefore, transcriptome-wide identification of RNA-binding proteins is necessary for deciphering the distinctive posttranscriptional mechanisms of gene regulation in Leishmania. Using a combination of highly efficient orthogonal organic phase separation method and tandem mass tag-labeling-based quantitative proteomic mass spectrometry, we provide unprecedented comprehensive molecular definition of the total RNA interactome across the two main Leishmania life cycle stages. In addition, we characterize for the first time the quantitative changes in RNA-protein interactions in Leishmania following inhibition of heat shock protein 90, shedding light into hitherto unknown large-scale downstream molecular effect of the protein inhibition in the parasite. This work provides insight into the importance of total RNA-protein interactions in Leishmania, thus significantly expanding our knowledge of the emergence of RNA-protein interactions in Leishmania biology.

In vivo architecture of the telomerase RNA catalytic core in Trypanosoma brucei.

Telomerase is a unique ribonucleoprotein (RNP) reverse transcriptase that utilizes its cognate RNA molecule as a template for telomere DNA repeat synthesis. Telomerase contains the reverse transcriptase protein, TERT and the template RNA, TR, as its core components. The 5'-half of TR forms a highly conserved catalytic core comprising of the template region and adjacent domains necessary for telomere synthesis. However, how telomerase RNA folding takes place in vivo has not been fully understood due to low abundance of the native RNP. Here, using unicellular pathogen Trypanosoma brucei as a model, we reveal important regional folding information of the native telomerase RNA core domains, i.e. TR template, template boundary element, template proximal helix and Helix IV (eCR4-CR5) domain. For this purpose, we uniquely combined in-cell probing with targeted high-throughput RNA sequencing and mutational mapping under three conditions: in vivo (in WT and TERT-/- cells), in an immunopurified catalytically active telomerase RNP complex and ex vivo (deproteinized). We discover that TR forms at least two different conformers with distinct folding topologies in the insect and mammalian developmental stages of T. brucei. Also, TERT does not significantly affect the RNA folding in vivo, suggesting that the telomerase RNA in T. brucei exists in a conformationally preorganized stable structure. Our observed differences in RNA (TR) folding at two distinct developmental stages of T. brucei suggest that important conformational changes are a key component of T. brucei development.

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.

G-quadruplex RNA motifs influence gene expression in the malaria parasite Plasmodium falciparum.

G-quadruplexes are non-helical secondary structures that can fold in vivo in both DNA and RNA. In human cells, they can influence replication, transcription and telomere maintenance in DNA, or translation, transcript processing and stability of RNA. We have previously showed that G-quadruplexes are detectable in the DNA of the malaria parasite Plasmodium falciparum, despite a very highly A/T-biased genome with unusually few guanine-rich sequences. Here, we show that RNA G-quadruplexes can also form in P. falciparum RNA, using rG4-seq for transcriptome-wide structure-specific RNA probing. Many of the motifs, detected here via the rG4seeker pipeline, have non-canonical forms and would not be predicted by standard in silico algorithms. However, in vitro biophysical assays verified formation of non-canonical motifs. The G-quadruplexes in the P. falciparum transcriptome are frequently clustered in certain genes and associated with regions encoding low-complexity peptide repeats. They are overrepresented in particular classes of genes, notably those that encode PfEMP1 virulence factors, stress response genes and DNA binding proteins. In vitro translation experiments and in vivo measures of translation efficiency showed that G-quadruplexes can influence the translation of P. falciparum mRNAs. Thus, the G-quadruplex is a novel player in post-transcriptional regulation of gene expression in this major human pathogen.

Malaria Parasite Stress Tolerance Is Regulated by DNMT2-Mediated tRNA Cytosine Methylation.

Malaria parasites need to cope with changing environmental conditions that require strong countermeasures to ensure pathogen survival in the human and mosquito hosts. The molecular mechanisms that protect Plasmodium falciparum homeostasis during the complex life cycle remain unknown. Here, we identify cytosine methylation of tRNAAsp (GTC) as being critical to maintain stable protein synthesis. Using conditional knockout (KO) of a member of the DNA methyltransferase family, called Pf-DNMT2, RNA bisulfite sequencing demonstrated the selective cytosine methylation of this enzyme of tRNAAsp (GTC) at position C38. Although no growth defect on parasite proliferation was observed, Pf-DNMT2KO parasites showed a selective downregulation of proteins with a GAC codon bias. This resulted in a significant shift in parasite metabolism, priming KO parasites for being more sensitive to various types of stress. Importantly, nutritional stress made tRNAAsp (GTC) sensitive to cleavage by an unknown nuclease and increased gametocyte production (>6-fold). Our study uncovers an epitranscriptomic mechanism that safeguards protein translation and homeostasis of sexual commitment in malaria parasites. IMPORTANCE P. falciparum is the most virulent malaria parasite species, accounting for the majority of the disease mortality and morbidity. Understanding how this pathogen is able to adapt to different cellular and environmental stressors during its complex life cycle is crucial in order to develop new strategies to tackle the disease. In this study, we identified the writer of a specific tRNA cytosine methylation site as a new layer of epitranscriptomic regulation in malaria parasites that regulates the translation of a subset of parasite proteins (>400) involved in different metabolic pathways. Our findings give insight into a novel molecular mechanism that regulates P. falciparum response to drug treatment and sexual commitment.