Transfer RNA modifications and genes for modifying enzymes in Arabidopsis thaliana

RNA modification: A contemporary review of pseudouridine (Ψ) and its role in functional plant biology

Pseudouridine (Ψ) is a modified nucleoside present in diverse RNA species, including mRNA (messenger RNA), snRNA (small nuclear RNA), rRNA (ribosomal RNA) and tRNA (transfer RNA). In plants, Ψ serves a critical function in RNA modification, supporting the stability, structural integrity, and functionality of RNA molecules. This review provides the various roles that Ψ fulfils in the modification of plant RNA biology, encompassing effects on biosynthesis pathways, regulatory mechanisms, stability, and translation efficiency. Additionally, we discuss recent advancements in the dynamic regulation of Ψ deposition in response to environmental stimuli and stressors. Elucidating Ψ's roles contributes to the comprehension of plant biology and may facilitate developments in biotechnology and crop improvement.

Unveiling transfer RNA modifications of oil palm and their dynamic changes during fruit ripening

BACKGROUND

The oil palm (Elaeis guineensis) is a crucial agricultural commodity, yielding the highest oil output among oil-bearing crops. Despite its significance, productivity challenges persist due to genetic and environmental factors. This study breaks new ground by mapping tRNA modifications in oil palm, exploring their roles during fruit ripening, an area not extensively studied in non-model crops.

RESULTS

Utilizing advanced RNA mass spectrometry techniques, we identified 48 distinct tRNA modifications across 88 sites, alongside 164 genes associated with tRNA modifying enzymes. This comprehensive mapping reveals the decreasing nature of most tRNA modifications during fruit development, except for adenosine 2'-O methylation (Am). It hints at a gradual weakening of protein translation quality control and highlights a unique role for Am during fruit ripening. Additionally, lipidomic analysis tracked 674 lipids in oil palm fruits, indicating a correlation between tRNA modifications and the accumulation of specific lipids.

CONCLUSIONS

This study mapped tRNA modifications in oil palm for the first time and showed the diversity of dynamic changes in tRNA modifications as the fruits develop.

AtTRM11 as a tRNA 2-methylguanosine methyltransferase modulates flowering and bacterial resistance via translational regulation

2-methylguanosine is an eukaryote-specific modified nucleoside in transfer RNAs, and m2G10 is catalyzed by Trm11-Trm112 protein complex in eukaryotic tRNAs. Here, we show that loss-of-function mutation of the Arabidopsis Trm11 homolog AtTRM11 resulted in m2G deficiency associated with disturbed ribosome assembly and overall transcriptome changes, including genes involved in flowering regulation and plant-pathogen interaction. The attrm11 mutant showed phenotypes of enlarged rosette leaves and early flowering, as well as enhanced resistance to Pseudomonas bacterial infection. AtTRM11 could partially rescue the m2G nucleoside level in yeast trm11 mutant, and AtTRM11 protein mostly resided in cytosol and physically interacted with AtTRM112b in planta. AtTRM11 was mostly expressed in shoot apex, root tip, and distal end of rosette leaves. KEGG enrichment analysis of differentially expressed genes between trm11 mutant and wild type indicated changes in pathways including phenopropanoid biosynthesis, plant-pathogen interaction, plant hormone signal transduction and MAPK signaling, suggesting that the pleiotropic phenotypes of the attrm11 mutant can be ascribed to translational and transcriptional changes.

RNA modifications in plant adaptation to abiotic stresses

Epitranscriptomic chemical modifications of RNAs have emerged as potent regulatory mechanisms in the process of plant stress adaptation. Currently, over 170 distinct chemical modifications have been identified in mRNAs, tRNAs, rRNAs, microRNAs (miRNAs), and long noncoding RNAs (lncRNAs). Genetic and molecular studies have identified the genes responsible for addition and removal of chemical modifications from RNA molecules, which are known as "writers" and "erasers," respectively. N6-methyladenosine (m6A) is the most prevalent chemical modification identified in eukaryotic mRNAs. Recent studies have identified m6A writers and erasers across different plant species, including Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), cotton (Gossypium hirsutum), and tomato (Solanum lycopersicum). Accumulating discoveries have improved our understanding of the functions of RNA modifications in plant stress responses. This review highlights the latest research on RNA modification, emphasizing the biological and cellular roles of diverse chemical modifications of mRNAs, tRNAs, rRNAs, miRNAs, and lncRNAs in plant responses to environmental and hormonal signals. We also propose and discuss critical questions and future challenges for enhancing our understanding of the cellular and mechanistic roles of RNA modifications in plant stress responses. Integrating molecular insights into the regulatory roles of RNA modifications in stress responses with novel genome- and RNA-editing technologies will facilitate the breeding of stress-tolerant crops through precise engineering of RNA modifications.

The m6A-YTH regulatory system in plants: A status

Plants use mRNA methylation to regulate gene expression. As in other eukaryotes, the only abundant methylated nucleotide in plant mRNA bodies is N6-methyladenosine (m6A). The conserved core components of m6A-based genetic control are a multi-subunit nuclear methyltransferase, and a set of nuclear and cytoplasmic RNA-binding proteins consisting of an m6A recognition module, the YT521-B homology (YTH) domain, and long intrinsically disordered regions (IDRs). In plants, this system is essential for growth during embryonic and post-embryonic development, but emerging evidence also points to key functions in plant-virus interactions and stimulus-dependent gene regulation. Cytoplasmic YTH-domain proteins are particularly important for these functions, and recent progress has identified two elements of the underlying molecular mechanisms: IDR-mediated phase separation and conserved short linear motifs mediating interactions with other key mRNA-binding proteins.

Molecular basis of A. thaliana KEOPS complex in biosynthesizing tRNA t6A

In archaea and eukaryotes, the evolutionarily conserved KEOPS is composed of four core subunits-Kae1, Bud32, Cgi121 and Pcc1, and a fifth Gon7/Pcc2 that is found in fungi and metazoa. KEOPS cooperates with Sua5/YRDC to catalyze the biosynthesis of tRNA N6-threonylcarbamoyladenosine (t6A), an essential modification needed for fitness of cellular organisms. Biochemical and structural characterizations of KEOPSs from archaea, yeast and humans have determined a t6A-catalytic role for Kae1 and auxiliary roles for other subunits. However, the precise molecular workings of KEOPSs still remain poorly understood. Here, we investigated the biochemical functions of A. thaliana KEOPS and determined a cryo-EM structure of A. thaliana KEOPS dimer. We show that A. thaliana KEOPS is composed of KAE1, BUD32, CGI121 and PCC1, which adopts a conserved overall arrangement. PCC1 dimerization leads to a KEOPS dimer that is needed for an active t6A-catalytic KEOPS-tRNA assembly. BUD32 participates in direct binding of tRNA to KEOPS and modulates the t6A-catalytic activity of KEOPS via its C-terminal tail and ATP to ADP hydrolysis. CGI121 promotes the binding of tRNA to KEOPS and potentiates the t6A-catalytic activity of KEOPS. These data and findings provide insights into mechanistic understanding of KEOPS machineries.

C2-methyladenosine in tRNA promotes protein translation by facilitating the decoding of tandem m2A-tRNA-dependent codons

RNA modification C2-methyladenosine (m2A) exists in both rRNA and tRNA of Escherichia coli (E. coli), installed by the methyltransferase RlmN using a radical-S-adenosylmethionine (SAM) mechanism. However, the precise function of m2A in tRNA and its ubiquity in plants have remained unclear. Here we discover the presence of m2A in chloroplast rRNA and tRNA, as well as cytosolic tRNA, in multiple plant species. We identify six m2A-modified chloroplast tRNAs and two m2A-modified cytosolic tRNAs across different plants. Furthermore, we characterize three Arabidopsis m2A methyltransferases-RLMNL1, RLMNL2, and RLMNL3-which methylate chloroplast rRNA, chloroplast tRNA, and cytosolic tRNA, respectively. Our findings demonstrate that m2A37 promotes a relaxed conformation of tRNA, enhancing translation efficiency in chloroplast and cytosol by facilitating decoding of tandem m2A-tRNA-dependent codons. This study provides insights into the molecular function and biological significance of m2A, uncovering a layer of translation regulation in plants.

Glycosylated queuosines in tRNAs optimize translational rate and post-embryonic growth

Transfer RNA (tRNA) modifications are critical for protein synthesis. Queuosine (Q), a 7-deaza-guanosine derivative, is present in tRNA anticodons. In vertebrate tRNAs for Tyr and Asp, Q is further glycosylated with galactose and mannose to generate galQ and manQ, respectively. However, biogenesis and physiological relevance of Q-glycosylation remain poorly understood. Here, we biochemically identified two RNA glycosylases, QTGAL and QTMAN, and successfully reconstituted Q-glycosylation of tRNAs using nucleotide diphosphate sugars. Ribosome profiling of knockout cells revealed that Q-glycosylation slowed down elongation at cognate codons, UAC and GAC (GAU), respectively. We also found that galactosylation of Q suppresses stop codon readthrough. Moreover, protein aggregates increased in cells lacking Q-glycosylation, indicating that Q-glycosylation contributes to proteostasis. Cryo-EM of human ribosome-tRNA complex revealed the molecular basis of codon recognition regulated by Q-glycosylations. Furthermore, zebrafish qtgal and qtman knockout lines displayed shortened body length, implying that Q-glycosylation is required for post-embryonic growth in vertebrates.

Phosphate Limitation Responses in Marine Green Algae Are Linked to Reprogramming of the tRNA Epitranscriptome and Codon Usage Bias

Marine algae are central to global carbon fixation, and their productivity is dictated largely by resource availability. Reduced nutrient availability is predicted for vast oceanic regions as an outcome of climate change; however, there is much to learn regarding response mechanisms of the tiny picoplankton that thrive in these environments, especially eukaryotic phytoplankton. Here, we investigate responses of the picoeukaryote Micromonas commoda, a green alga found throughout subtropical and tropical oceans. Under shifting phosphate availability scenarios, transcriptomic analyses revealed altered expression of transfer RNA modification enzymes and biased codon usage of transcripts more abundant during phosphate-limiting versus phosphate-replete conditions, consistent with the role of transfer RNA modifications in regulating codon recognition. To associate the observed shift in the expression of the transfer RNA modification enzyme complement with the transfer RNAs encoded by M. commoda, we also determined the transfer RNA repertoire of this alga revealing potential targets of the modification enzymes. Codon usage bias was particularly pronounced in transcripts encoding proteins with direct roles in managing phosphate limitation and photosystem-associated proteins that have ill-characterized putative functions in "light stress." The observed codon usage bias corresponds to a proposed stress response mechanism in which the interplay between stress-induced changes in transfer RNA modifications and skewed codon usage in certain essential response genes drives preferential translation of the encoded proteins. Collectively, we expose a potential underlying mechanism for achieving growth under enhanced nutrient limitation that extends beyond the catalog of up- or downregulated protein-encoding genes to the cell biological controls that underpin acclimation to changing environmental conditions.

RNA pseudouridine modification in plants

Pseudouridine is one of the well-known chemical modifications in various RNA species. Current advances to detect pseudouridine show that the pseudouridine landscape is dynamic and affects multiple cellular processes. Although our understanding of this post-transcriptional modification mainly depends on yeast and human models, the recent findings provide strong evidence for the critical role of pseudouridine in plants. Here, we review the current knowledge of pseudouridine in plant RNAs, including its synthesis, degradation, regulatory mechanisms, and functions. Moreover, we propose future areas of research on pseudouridine modification in plants.

Detection of queuosine and queuosine precursors in tRNAs by direct RNA sequencing

Queuosine (Q) is a complex tRNA modification found in bacteria and eukaryotes at position 34 of four tRNAs with a GUN anticodon, and it regulates the translational efficiency and fidelity of the respective codons that differ at the Wobble position. In bacteria, the biosynthesis of Q involves two precursors, preQ0 and preQ1, whereas eukaryotes directly obtain Q from bacterial sources. The study of queuosine has been challenging due to the limited availability of high-throughput methods for its detection and analysis. Here, we have employed direct RNA sequencing using nanopore technology to detect the modification of tRNAs with Q and Q precursors. These modifications were detected with high accuracy on synthetic tRNAs as well as on tRNAs extracted from Schizosaccharomyces pombe and Escherichia coli by comparing unmodified to modified tRNAs using the tool JACUSA2. Furthermore, we present an improved protocol for the alignment of raw sequence reads that gives high specificity and recall for tRNAs ex cellulo that, by nature, carry multiple modifications. Altogether, our results show that 7-deazaguanine-derivatives such as queuosine are readily detectable using direct RNA sequencing. This advancement opens up new possibilities for investigating these modifications in native tRNAs, furthering our understanding of their biological function.

The Small RNA Component of Arabidopsis thaliana Phloem Sap and Its Response to Iron Deficiency

In order to discover sRNA that might function during iron deficiency stress, RNA was prepared from phloem exudates of Arabidopsis thaliana, and used for RNA-seq. Bioanalyzer results indicate that abundant RNA from phloem is small in size-less than 200 nt. Moreover, typical rRNA bands were not observed. Sequencing of eight independent phloem RNA samples indicated that tRNA-derived fragments, specifically 5' tRFs and 5' tRNA halves, are highly abundant in phloem sap, comprising about 46% of all reads. In addition, a set of miRNAs that are present in phloem sap was defined, and several miRNAs and sRNAs were identified that are differentially expressed during iron deficiency.

High-throughput profiling of RNA modifications by ultra-performance liquid chromatography coupled to complementary mass spectrometry: Methods, quality control, and applications

Although next-generation sequencing technology has been used to delineate RNA modifications in recent years, the paucity of appropriate converting reactions or specific antibodies impedes the accurate characterization and quantification of numerous RNA modifications, especially when these modifications demonstrate wide variations across developmental stages and cell types. In this study, we developed a high-throughput analytical platform coupling ultra-performance liquid chromatograph (UPLC) with complementary mass spectrometry (MS) to identify and quantify RNA modifications in both synthetic and biological samples. Sixty-four types of RNA modifications, including positional isomers and hypermodified ribonucleosides, were successfully monitored within a 16-min single run of UPLC-MS. Two independent methods to cross-validate the purity of RNA extracted from Caenorhabditis elegans (C. elegans) were developed using the coexisting C. elegans and Escherichia coli (E. coli) as a surveillance system. To test the validity of the method, we investigated the RNA modification landscape of three model organisms, C. elegans, E. coli, and Arabidopsis thaliana (A. thaliana). Both the identity and molarity of modified ribonucleosides markedly varied among the species. Moreover, our platform is not only useful for exploring the dynamics of RNA modifications in response to environmental cues (e.g., cold shock) but can also help with the identification of RNA-modifying enzymes in genetic studies. Cumulatively, our method presents a novel platform for the comprehensive analysis of RNA modifications, which will be of benefit to both analytical chemists involved in biomarker discovery and biologists conducting functional studies of RNA modifications.

tRNA derived small RNAs—Small players with big roles

In the past 2 decades, small non-coding RNAs derived from tRNA (tsRNAs or tRNA derived fragments; tRFs) have emerged as new powerful players in the field of small RNA mediated regulation of gene expression, translation, and epigenetic control. tRFs have been identified from evolutionarily divergent organisms from Archaea, the higher plants, to humans. Recent studies have confirmed their roles in cancers and other metabolic disorders in humans and experimental models. They have been implicated in biotic and abiotic stress responses in plants as well. In this review, we summarize the current knowledge on tRFs including types of tRFs, their biogenesis, and mechanisms of action. The review also highlights recent studies involving differential expression profiling of tRFs and elucidation of specific functions of individual tRFs from various species. We also discuss potential considerations while designing experiments involving tRFs identification and characterization and list the available bioinformatics tools for this purpose.

RNA Sequencing Unveils Very Small RNAs With Potential Regulatory Functions in Bacteria

RNA sequencing (RNA-seq) is the gold standard for the discovery of small non-coding RNAs. Following a long-standing approach, reads shorter than 16 nucleotides (nt) are removed from the small RNA sequencing libraries or datasets. The serendipitous discovery of an eukaryotic 12 nt-long RNA species capable of modulating the microRNA from which they derive prompted us to challenge this dogma and, by expanding the window of RNA sizes down to 8 nt, to confirm the existence of functional very small RNAs (vsRNAs <16 nt). Here we report the detailed profiling of vsRNAs in Escherichia coli, E. coli-derived outer membrane vesicles (OMVs) and five other bacterial strains (Pseudomonas aeruginosa PA7, P. aeruginosa PAO1, Salmonella enterica serovar Typhimurium 14028S, Legionella pneumophila JR32 Philadelphia-1 and Staphylococcus aureus HG001). vsRNAs of 8–15 nt in length [RNAs (8-15 nt)] were found to be more abundant than RNAs of 16–30 nt in length [RNAs (16–30 nt)]. vsRNA biotypes were distinct and varied within and across bacterial species and accounted for one third of reads identified in the 8–30 nt window. The tRNA-derived fragments (tRFs) have appeared as a major biotype among the vsRNAs, notably Ile-tRF and Ala-tRF, and were selectively loaded in OMVs. tRF-derived vsRNAs appear to be thermodynamically stable with at least 2 G-C basepairs and stem-loop structure. The analyzed tRF-derived vsRNAs are predicted to target several human host mRNAs with diverse functions. Bacterial vsRNAs and OMV-derived vsRNAs could be novel players likely modulating the intricate relationship between pathogens and their hosts.

The Dihydrouridine landscape from tRNA to mRNA: a perspective on synthesis, structural impact and function

The universal dihydrouridine (D) epitranscriptomic mark results from a reduction of uridine by the Dus family of NADPH-dependent reductases and is typically found within the eponym D-loop of tRNAs. Despite its apparent simplicity, D is structurally unique, with the potential to deeply affect the RNA backbone and many, if not all, RNA-connected processes. The first landscape of its occupancy within the tRNAome was reported 20 years ago. Its potential biological significance was highlighted by observations ranging from a strong bias in its ecological distribution to the predictive nature of Dus enzymes overexpression for worse cancer patient outcomes. The exquisite specificity of the Dus enzymes revealed by a structure-function analyses and accumulating clues that the D distribution may expand beyond tRNAs recently led to the development of new high-resolution mapping methods, including Rho-seq that established the presence of D within mRNAs and led to the demonstration of its critical physiological relevance.

Data-Independent Acquisition for the Detection of Mononucleoside RNA Modifications by Mass Spectrometry

RNA is dynamically modified in cells by a plethora of chemical moieties to modulate molecular functions and processes. Over 140 modifications have been identified across species and RNA types, with the highest density and diversity of modifications found in tRNA (tRNA). The methods used to identify and quantify these modifications have developed over recent years and continue to advance, primarily in the fields of next-generation sequencing (NGS) and mass spectrometry (MS). Most current NGS methods are limited to antibody-recognized or chemically derivatized modifications and have limitations in identifying multiple modifications simultaneously. Mass spectrometry can overcome both of these issues, accurately identifying a large number of modifications in a single run. Here, we present advances in MS data acquisition for the purpose of RNA modification identification and quantitation. Using this approach, we identified multiple tRNA wobble position modifications in Arabidopsis thaliana that are upregulated in salt-stressed growth conditions and may stabilize translation of salt stress induced proteins. This work presents improvements in methods for studying RNA modifications and introduces a possible regulatory role of wobble position modifications in A. thaliana translation.

Epitranscriptomics: An Additional Regulatory Layer in Plants' Development and Stress Response

Epitranscriptomics has added a new layer of regulatory machinery to eukaryotes, and the advancement of sequencing technology has revealed more than 170 post-transcriptional modifications in various types of RNAs, including messenger RNA (mRNA), transfer RNA (tRNA), ribosomal RNA (rRNA), and long non-coding RNA (lncRNA). Among these, N6-methyladenosine (m6A) and N5-methylcytidine (m5C) are the most prevalent internal mRNA modifications. These regulate various aspects of RNA metabolism, mainly mRNA degradation and translation. Recent advances have shown that regulation of RNA fate mediated by these epitranscriptomic marks has pervasive effects on a plant's development and responses to various biotic and abiotic stresses. Recently, it was demonstrated that the removal of human-FTO-mediated m6A from transcripts in transgenic rice and potatoes caused a dramatic increase in their yield, and that the m6A reader protein mediates stress responses in wheat and apple, indicating that regulation of m6A levels could be an efficient strategy for crop improvement. However, changing the overall m6A levels might have unpredictable effects; therefore, the identification of precise m6A levels at a single-base resolution is essential. In this review, we emphasize the roles of epitranscriptomic modifications in modulating molecular, physiological, and stress responses in plants, and provide an outlook on epitranscriptome engineering as a promising tool to ensure food security by editing specific m6A and m5C sites through robust genome-editing technology.

TRM61 is essential for Arabidopsis embryo and endosperm development

Post-transcriptional modifications of tRNA molecules play crucial roles in gene expression and protein biosynthesis. Across the genera, methylation of tRNAs at N¹ of adenosine 58 (A58) by AtTRM61/AtTRM6 complex plays a critical role in maintaining the stability of initiator methionyl-tRNA (tRNA_i^Met). Recently, it was shown that mutation in AtTRM61 or AtTRM6 leads to seed abortion. However, a detailed study about the AtTRM61/AtTRM6 function in plants remains vague. Here, we found that AtTRM61 has a conserved functional structure and possesses conserved binding motifs for cofactor S-adenosyl-L-methionine (AdoMet). Mutations of the complex subunits AtTRM61/AtTRM6 result in embryo and endosperm developmental defects. The endosperm and embryo developmental defects were conditionally complemented by Attrm61-1/ + FIS2pro::AtTRM61 and Attrm61-1/ + ABI3pro::AtTRM61 indicating that AtTRM61 is required for early embryo and endosperm development. Besides, the rescue of the fertility defects in trm61/ + by overexpression of initiator tRNA suggests that AtTRM61 mutation could diminish tRNA_i^Met stability. Moreover, using yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays, we showed that AtMPK4 physically interacts with AtTRM61. The data presented here suggest that AtTRM61 has a conserved structure and function in Arabidopsis. Also, AtTRM61 may be required for tRNA_i^Met stability, embryo and endosperm development.

Overexpression of GhKTI12 Enhances Seed Yield and Biomass Production in Nicotiana Tabacum

Crop molecular breeding primarily focuses on increasing the trait of plant yield. An elongator-associated protein, KTI12, is closely associated with plant biomass and yield. KTI12 is involved in developmental processes of most organs, including the leaf, root, flower, and seed, through regulating cell division and differentiation. Previous work has shown that in upland cotton (Gossypium hirsutum), GhKTI12 regulates plant height, flowering, and tolerance to salt and drought stress. However, little is known about the molecular regulation mechanism of GhKTI12 in plant developmental processes. In this study, we identified the main GhKTI12 (Gh_D02G144400) gene and transformed it into tobacco (Nicotonia tabacum cv NC89). From seven transgenic lines, we obtained three (OE5, OE6 and OE8) with high expression of GhKTI12; compared with wild type plants, these three lines exhibited larger plant size, later flowering, and higher seed yield. Microscopic observation revealed that the number of leaf epidermal cells and stem parenchyma cells was increased by ~55%. Biochemical analysis showed that chlorophyll content and starch accumulation were significantly increased in younger leaves at the top canopy of transgenic plants, which may contribute to improved photosynthetic rate and, in turn, increased seed yield. To understand the molecular mechanism of GhKTI12 in transgenic plants development, two lines (OE6 and OE8) with higher expression levels of GhKTI12 were used as representative plants to conduct RNA-seq analysis. Through transcriptome analysis of the plant's shoot apical meristematic tissue of these two lines, we identified 518 upregulated genes and 406 downregulated genes common to both overexpression lines. A large number of cellular component genes associated with cell division and differentiation, such as RD21, TET8, KTN80, AOX1, AOX2, CP1, and KIC, were found to be upregulated, and genes showing the most downregulation included MADS-box genes related to flowering time, such as MADS6, AP1, AP3, AGL8, AGL6, SEP1, and SEP2. Downregulation of these genes caused delayed flowering time and longer vegetative stage during development. Combined with the upregulation of the yield-related gene RD21, the GhKTI12 transgenic plants could produce a higher seed yield. We here show that the overexpression of GhKTI12 could positively improve key agronomic traits in tobacco by regulating cell proliferation, photosynthesis, and organ development, and suggest that homologs of GhKTI12 may also be important in the genetic improvement of other crop plants.

Variations in transfer and ribosomal RNA epitranscriptomic status can adapt eukaryote translation to changing physiological and environmental conditions

The timely reprogramming of gene expression in response to internal and external cues is essential to eukaryote development and acclimation to changing environments. Chemically modifying molecular receptors and transducers of these signals is one way to efficiently induce proper physiological responses. Post-translation modifications, regulating protein biological activities, are central to many well-known signal-responding pathways. Recently, messenger RNA (mRNA) chemical (i.e. epitranscriptomic) modifications were also shown to play a key role in these processes. In contrast, transfer RNA (tRNA) and ribosomal RNA (rRNA) chemical modifications, although critical for optimal function of the translation apparatus, and much more diverse and quantitatively important compared to mRNA modifications, were until recently considered as mainly static chemical decorations. We present here recent observations that are challenging this view and supporting the hypothesis that tRNA and rRNA modifications dynamically respond to various cell and environmental conditions and contribute to adapt translation to these conditions.

RsmD, a Chloroplast rRNA m2G Methyltransferase, Plays a Role in Cold Stress Tolerance by Possibly Affecting Chloroplast Translation in Arabidopsis

Ribosomal RNA (rRNA) methylation is a pivotal process in the assembly and activity of ribosomes, which in turn play vital roles in the growth, development and stress responses of plants. Although few methyltransferases responsible for rRNA methylation have been identified in plant chloroplasts, the nature and function of these enzymes in chloroplasts remain largely unknown. In this study, we characterized ArabidopsisRsmD (At3g28460), an ortholog of the methyltransferase responsible for N2-methylguanosine (m2G) modification of 16S rRNA in Escherichia coli. Confocal microscopic analysis of an RsmD- green fluorescent protein fusion protein revealed that RsmD is localized to chloroplasts. Primer extension analysis indicated that RsmD is responsible for m2G methylation at position 915 in the 16S rRNA of Arabidopsis chloroplasts. Under cold stress, rsmd mutant plants exhibited retarded growth, i.e. had shorter roots, lower fresh weight and pale-green leaves, compared with wild-type (WT) plants. However, these phenotypes were not detected in response to drought or salt stress. Notably, the rsmd mutant was hypersensitive to erythromycin or lincomycin and accumulated fewer chloroplast proteins compared with the WT, suggesting that RsmD influences translation in chloroplasts. Complementation lines expressing RsmD in the rsmd mutant background recovered WT phenotypes. Importantly, RsmD harbored RNA methyltransferase activity. Collectively, the findings of this study indicate that RsmD is a chloroplast 16S rRNA methyltransferase responsible for m2G915 modification that plays a role in the adaptation of Arabidopsisto cold stress.

Deciphering Epitranscriptome: Modification of mRNA Bases Provides a New Perspective for Post-transcriptional Regulation of Gene Expression

Gene regulation depends on dynamic and reversibly modifiable biological and chemical information in the epigenome/epitranscriptome. Accumulating evidence suggests that messenger RNAs (mRNAs) are generated in flashing bursts in the cells in a precisely regulated manner. However, the different aspects of the underlying mechanisms are not fully understood. Cellular RNAs are post-transcriptionally modified at the base level, which alters the metabolism of mRNA. The current understanding of epitranscriptome in the animal system is far ahead of that in plants. The accumulating evidence indicates that the epitranscriptomic changes play vital roles in developmental processes and stress responses. Besides being non-genetically encoded, they can be of reversible nature and involved in fine-tuning the expression of gene. However, different aspects of base modifications in mRNAs are far from adequate to assign the molecular basis/functions to the epitranscriptomic changes. Advances in the chemogenetic RNA-labeling and high-throughput next-generation sequencing techniques are enabling functional analysis of the epitranscriptomic modifications to reveal their roles in mRNA biology. Mapping of the common mRNA modifications, including N 6-methyladenosine (m6A), and 5-methylcytidine (m5C), have enabled the identification of other types of modifications, such as N 1-methyladenosine. Methylation of bases in a transcript dynamically regulates the processing, cellular export, translation, and stability of the mRNA; thereby influence the important biological and physiological processes. Here, we summarize the findings in the field of mRNA base modifications with special emphasis on m6A, m5C, and their roles in growth, development, and stress tolerance, which provide a new perspective for the regulation of gene expression through post-transcriptional modification. This review also addresses some of the scientific and technical issues in epitranscriptomic study, put forward the viewpoints to resolve the issues, and discusses the future perspectives of the research in this area.

PtKTI12 genes influence wobble uridine modifications and drought stress tolerance in hybrid poplar

The multisubunit Elongator complex plays key roles in transcription by interacting with RNA polymerase II and chromatin modeling. Kti proteins have been identified as the auxiliary protein for the Elongator complex. However, our knowledge of Kti proteins in woody plants remains limited. In this study, in total 16 KTI gene homologs were identified in Populus trichocarpa. Among them, the two KTI12 candidates were named PtKTI12A and PtKTI12B. Although PtKTI12A and PtKTI12B were largely different in gene expression level and tissue specificity, both genes were induced by heat and drought stresses. PtKTI12A and PtKTI12B RNAi transgenic poplar plants showed reduced levels of modified nucleosides, in particular 5-carbamoylmethyluridine and 5-methoxycarbonylmethyl-2-thiouridine. Meanwhile, their tolerance to drought was improved when subjected to withdrawal of watering. Also, the protein products of PtKTI12A and PtKTI12B had similar subcellular localization and predicted tertiary structure. The results suggest that Kti12 proteins are involved in tRNA wobble uridine modification, stress response and drought stress tolerance in hybrid poplar.

Plant Elongator-Protein Complex of Diverse Activities Regulates Growth, Development, and Immune Responses

Contrary to the conserved Elongator composition in yeast, animals, and plants, molecular functions and catalytic activities of the complex remain controversial. Elongator was identified as a component of elongating RNA polymerase II holoenzyme in yeast, animals, and plants. Furthermore, it was suggested that Elonagtor facilitates elongation of transcription via histone acetyl transferase activity. Accordingly, phenotypes of Arabidopsis elo mutants, which show development, growth, or immune response defects, correlate with transcriptional downregulation and the decreased histone acetylation in the coding regions of crucial genes. Plant Elongator was also implicated in other processes: transcription and processing of miRNA, regulation of DNA replication by histone acetylation, and acetylation of alpha-tubulin. Moreover, tRNA modification, discovered first in yeast and confirmed in plants, was claimed as the main activity of Elongator, leading to specificity in translation that might also result indirectly in a deficiency in transcription. Heterologous overexpression of individual Arabidopsis Elongator subunits and their respective phenotypes suggest that single Elongator subunits might also have another function next to being a part of the complex. In this review, we shall present the experimental evidence of all molecular mechanisms and catalytic activities performed by Elongator in nucleus and cytoplasm of plant cells, which might explain how Elongator regulates growth, development, and immune responses.

The Arabidopsis TRM61/TRM6 complex is a bona fide tRNA N1-methyladenosine methyltransferase

tRNA molecules, which contain the most abundant post-transcriptional modifications, are crucial for proper gene expression and protein biosynthesis. Methylation at N1 of adenosine 58 (A58) is critical for maintaining the stability of initiator methionyl-tRNA (tRNAiMet) in bacterial, archaeal, and eukaryotic tRNAs. However, although research has been conducted in yeast and mammals, it remains unclear how A58 in plant tRNAs is modified and involved in development. In this study, we identify the nucleus-localized complex AtTRM61/AtTRM6 in Arabidopsis as tRNA m1A58 methyltransferase. Deficiency or a lack of either AtTRM61 or AtTRM6 leads to embryo arrest and seed abortion. The tRNA m1A level decreases in conditionally complemented Attrm61/LEC1pro::AtTRM61 plants and this is accompanied by reduced levels of tRNAiMet, indicating the importance of the tRNA m1A modification for tRNAiMet stability. Taken together, our results demonstrate that tRNA m1A58 modification is necessary for tRNAiMet stability and is required for embryo development in Arabidopsis.

Multi-faceted nature of the tRNA isopentenyltransferase

Transfer RNA isopentenylation an adenine 37 position (A37) is a universal modification known in prokaryotes and eukaryotes. A set of highly homologous enzymes catalyse a series of reactions, leading to tRNA modifications, aimed to increase adaptation to environmental condition through the control of translation efficiency and reading frame maintenance. Transfer RNA-isopentenylation-related (TI-related) functions are well studied in bacteria, mitochondria of yeast and human, but completely unexplored in plants. Transfer RNA-isopentenylation-unrelated (TI-unrelated) functions participate in adaptation to environmental stresses via the regulation of sterol metabolism, gene silencing/suppression and amyloid fibrils formation. TI-unrelated functions are mostly studied in yeast. Finally, the degradation of A37-modified tRNA releases a set of bioactive compounds known as cis-cytokinins. Although all organisms are able to produce cis-cytokinins, its physiological role is still a matter of debates. For several species of bacteria and fungi, cis-cytokinins are known to play a crucial role in pathogenesis. In mammalian and human models cis-cytokinins have tumour-suppressing and anti-inflammation effects. This review aims to summarise current knowledge of the TI-related and TI-unrelated functions and main bioactive by-products of isopentenylated tRNA degradation.

LC-MS/MS Profiling of Post-Transcriptional Modifications in Ginseng tRNA Purified by a Polysaccharase-Aided Extraction Method

Transfer RNAs (tRNAs) are the most heavily modified RNA species in life entities. Post-transcriptional modifications severely impact the structure and function of tRNAs. To date, hundreds of modifications have been identified in tRNAs, mainly from microorganisms and animals. However, tRNAs in plant roots or tubers that have been widely used for food and medical purpose for centuries are rarely studied because isolation of RNA from plants still remains a challenge. In this paper, a polysaccharase-aided RNA isolation (PARI) method for extraction of high-quality RNA from plants containing large quantities of polysaccharides is developed. This method presents a new strategy of “digesting” polysaccharides that is completely different from the conventional method of “dissolving” the contaminants. By using this method, RNA of high integrity and purity were successfully extracted from ginseng roots because polysaccharide contaminations were removed efficiently with α-amylase digestion. Ginseng tRNAs were first sequenced by NGS and a total of 41 iso acceptors were identified. ChloroplastictRNA^Gly(GCC) in ginseng root was purified and four modified nucleosides, including m⁷G, D, T, and Ψ, were identified by LC-MS/MS. The results also revealed that the m⁷G occurs at a novel position 18, which may be related to the deformation of D-loop. PARI is the first enzyme-assisted technique for RNA isolation from plants, which could fundamentally solve the problem of polysaccharide contaminations. By using the PARI method, more individual tRNAs could be isolated easily from polysaccharide-rich plant tissues, which would have a positive impact on the feasibility of research on structure and function of tRNA in plants.

Transcriptome Analysis of ppdnmt2 and Identification of Superoxide Dismutase as a Novel Interactor of DNMT2 in the Moss Physcomitrella patens

DNMT2 is a DNA/tRNA cytosine methyltransferase that is highly conserved in structure and function in eukaryotes. In plants however, limited information is available on the function of this methyltransferase. We have previously reported that in the moss Physcomitrella patens, DNMT2 plays a crucial role in stress recovery and tRNA^Asp transcription/stability under salt stress. To further investigate the role of PpDNMT2 at genome level, in this study we have performed RNA sequencing of ppdnmt2. Transcriptome analysis reveals a number of genes and pathways to function differentially and suggests a close link between PpDNMT2 function and osmotic and ionic stress tolerance. We propose PpDNMT2 to play a pivotal role in regulating salt tolerance by affecting molecular networks involved in stress perception and signal transduction that underlie maintenance of ion homeostasis in cells. We also examined interactome of PpDNMT2 using affinity purification (AP) coupled to mass spectrometry (AP-MS). Quantitative proteomic analysis reveals several chloroplast proteins involved in light reactions and carbon assimilation and proteins involved in stress response and some not implicated in stress to co-immunoprecipitate with PpDNMT2. Comparison between transcriptome and interactome datasets has revealed novel association between PpDNMT2 activity and the antioxidant enzyme Superoxide dismutase (SOD), protein turnover mediated by the Ubiquitin-proteasome system and epigenetic gene regulation. PpDNMT2 possibly exists in complex with CuZn-SODs in vivo and the two proteins also directly interact in the yeast nucleus as observed by yeast two-hybrid assay. Taken together, the work presented in this study sheds light on diverse roles of PpDNMT2 in maintaining molecular and physiological homeostasis in P. patens. This is a first report describing transcriptome and interactome of DNMT2 in any land plant.

AtTrm5a catalyses 1-methylguanosine and 1-methylinosine formation on tRNAs and is important for vegetative and reproductive growth in Arabidopsis thaliana

Modified nucleosides on tRNA are critical for decoding processes and protein translation. tRNAs can be modified through 1-methylguanosine (m¹G) on position 37; a function mediated by Trm5 homologs. We show that AtTRM5a (At3g56120) is a Trm5 ortholog in Arabidopsis thaliana. AtTrm5a is localized to the nucleus and its function for m¹G and m¹I methylation was confirmed by mutant analysis, yeast complementation, m¹G nucleoside level on single tRNA, and tRNA in vitro methylation. Arabidopsis attrm5a mutants were dwarfed and had short filaments, which led to reduced seed setting. Proteomics data indicated differences in the abundance of proteins involved in photosynthesis, ribosome biogenesis, oxidative phosphorylation and calcium signalling. Levels of phytohormone auxin and jasmonate were reduced in attrm5a mutant, as well as expression levels of genes involved in flowering, shoot apex cell fate determination, and hormone synthesis and signalling. Taken together, loss-of-function of AtTrm5a impaired m¹G and m¹I methylation and led to aberrant protein translation, disturbed hormone homeostasis and developmental defects in Arabidopsis plants.

Translational Regulation of Plant Response to High Temperature by a Dual-Function tRNAHis Guanylyltransferase in Rice

As sessile organisms, plants have evolved numerous strategies to acclimate to changes in environmental temperature. However, the molecular basis of this acclimation remains largely unclear. In this study we identified a tRNA^His guanylyltransferase, AET1, which contributes to the modification of pre-tRNA^His and is required for normal growth under high-temperature conditions in rice. Interestingly, AET1 possibly interacts with both RACK1A and eIF3h in the endoplasmic reticulum. Notably, AET1 can directly bind to OsARF mRNAs including the uORFs of OsARF19 and OsARF23, indicating that AET1 is associated with translation regulation. Furthermore, polysome profiling assays suggest that the translational status remains unaffected in the aet1 mutant, but that the translational efficiency of OsARF19 and OsARF23 is reduced; moreover, OsARF23 protein levels are obviously decreased in the aet1 mutant under high temperature, implying that AET1 regulates auxin signaling in response to high temperature. Our findings provide new insights into the molecular mechanisms whereby AET1 regulates the environmental temperature response in rice by playing a dual role in tRNA modification and translational control.

CyMIRA: The Cytonuclear Molecular Interactions Reference for Arabidopsis

The function and evolution of eukaryotic cells depend upon direct molecular interactions between gene products encoded in nuclear and cytoplasmic genomes. Understanding how these cytonuclear interactions drive molecular evolution and generate genetic incompatibilities between isolated populations and species is of central importance to eukaryotic biology. Plants are an outstanding system to investigate such effects because of their two different genomic compartments present in the cytoplasm (mitochondria and plastids) and the extensive resources detailing subcellular targeting of nuclear-encoded proteins. However, the field lacks a consistent classification scheme for mitochondrial- and plastid-targeted proteins based on their molecular interactions with cytoplasmic genomes and gene products, which hinders efforts to standardize and compare results across studies. Here, we take advantage of detailed knowledge about the model angiosperm Arabidopsis thaliana to provide a curated database of plant cytonuclear interactions at the molecular level. CyMIRA (Cytonuclear Molecular Interactions Reference for Arabidopsis) is available at http://cymira.colostate.edu/ and https://github.com/dbsloan/cymira and will serve as a resource to aid researchers in partitioning evolutionary genomic data into functional gene classes based on organelle targeting and direct molecular interaction with cytoplasmic genomes and gene products. It includes 11 categories (and 27 subcategories) of different cytonuclear complexes and types of molecular interactions, and it reports residue-level information for cytonuclear contact sites. We hope that this framework will make it easier to standardize, interpret, and compare studies testing the functional and evolutionary consequences of cytonuclear interactions.

NEP-TC a rRNA Methyltransferase Involved on Somatic Embryogenesis of Tamarillo (Solanum betaceum Cav.)

Somatic embryogenesis (SE) is an important biotechnological tool for large-scale clonal propagation and for embryogenesis research. Moreover, genetic transformation and cryopreservation procedures in many species rely on efficient SE protocols. We have been studying different aspects related to SE induction and somatic embryo development in tamarillo (Solanum betaceum Cav.), a small tree from the Solanaceae family. Previous proteomic analyses identified a protein (NEP-TC, 26.5 kDa) consistently present in non-embryogenic calluses of tamarillo, but absent in the embryogenic ones. In this work, the role of NEP-TC during SE was assessed by gene expression analysis and immunolocalization. The results obtained demonstrated that NEP-TC is a putative member of the SpoU rRNA methylase family. This protein, present in the cytoplasm and nucleus, is expressed in non-embryogenic cells and not expressed in embryogenic cells. Slightly enhanced SE induction levels in tamarillo plants with NEP-TC down-regulated levels also supports the role of this protein on SE induction. Heterologous expression was used to confirm NEP-TC rRNA methyltransferase activity, with enhanced activity levels when rRNA was used as a substrate. These data relate a putative member of the SpoU methylase family with plant morphogenesis, in particular with SE induction.

Messenger RNA Modifications in Plants

Over 160 distinct RNA modifications are known and collectively termed the epitranscriptome. Some of these modifications have been discovered in mRNA, uncovering a new layer of gene regulation. Transcriptome-wide mapping of epitranscriptomic codes and the discovery of their writers, erasers, and readers that dynamically install, remove, and interpret RNA modifications, respectively, are fundamental to understanding the epitranscriptome. Recent technologies have enabled the transcriptome-wide profiling of several mRNA modifications in Arabidopsis thaliana, providing key insights into regulating these modifications and their effects on plant development. Here we review technological innovations and recent progress in epitranscriptomics, with specific focus on N⁶-methyladenosine (m⁶A), 5-methylcytosine (m⁵C), uridylation, and their roles in multiple aspects of plant development.

Structural studies on dihydrouridine synthase A (DusA) from Pseudomonas aeruginosa

Dihydrouridination is one of the abundant modifications in tRNA editing. The presence of dihydrouridine is attributed to tRNA stability desired for the efficient gene translation process. The conversion of uridine to dihydrouridine is catalyzed by flavine containing enzyme called dihydrouridine synthase (Dus). We report first ever information about DusA enzyme from Pseudomonas aeruginosa in form of structural and functional studies. The gene coding for DusA from P. aeruginosa (PADusA) was cloned, expressed and purified, using recombinant DNA technology methods. Thermal and chemical stability of PADusA was determined with respect to temperature and urea-induced equilibrium unfolding experiments, with monitoring the change of ellipticity at 200-260 nm by Circular Dichroism (CD) spectroscopy. Unfolding studies revealed that PADusA has acquired a stable tertiary structure fold with a T_m value of 46.2 °C and C_m of 2.7 M for urea. The enzyme contains 43% α-helices and 16% β-strands. The three dimensional structure of PADusA was modeled using insilico methods. In order to understand the mechanism of substrate recognition and catalysis, tRNA and puromycin were docked on PADusA structure and their binding was analyzed. The structural features suggested that PADusA may also form a novel target for structure based drug design of antimicrobial agents.

A 2'-O-Methyltransferase Responsible for Transfer RNA Anticodon Modification Is Pivotal for Resistance to Pseudomonas syringae DC3000 in Arabidopsis

Transfer RNA (tRNA) is the most highly modified class of RNA species in all living organisms. Recent discoveries have revealed unprecedented complexity in the tRNA chemical structures, modification patterns, regulation, and function, suggesting that each modified nucleoside in tRNA may have its own specific function. However, in plants, our knowledge of the role of individual tRNA modifications and how they are regulated is very limited. In a genetic screen designed to identify factors regulating disease resistance in Arabidopsis, we identified SUPPRESSOR OF CSB3 9 (SCS9). Our results reveal SCS9 encodes a tRNA methyltransferase that mediates the 2'-O-ribose methylation of selected tRNA species in the anticodon loop. These SCS9-mediated tRNA modifications enhance susceptibility during infection with the virulent bacterial pathogen Pseudomonas syringae DC3000. Lack of such tRNA modification, as observed in scs9 mutants, specifically dampens plant resistance against DC3000 without compromising the activation of the salicylic acid signaling pathway or the resistance to other biotrophic pathogens. Our results support a model that gives importance to the control of certain tRNA modifications for mounting an effective disease resistance in Arabidopsis toward DC3000 and, therefore, expands the repertoire of molecular components essential for an efficient disease resistance response.

Fine-Tuning of Gene Expression by tRNA-Derived Fragments during Abiotic Stress Signal Transduction

When plants are subjected to unfavorable environmental conditions, overall gene expression in stressed cells is altered from a programmed pattern for normal development to an adaptive pattern for survival. Rapid changes in plant gene expression include production of stress responsive proteins for protection as well as reduction of irrelevant proteins to minimize energy consumption during growth. In addition to the many established mechanisms known to modulate gene expression in eukaryotes, a novel strategy involving tRNA-derived fragments (tRFs) was recently reported to control gene expression. In animals, tRFs are shown to play a certain role in infected or cancer cells. However, tRFs are expected to function in the regulation of gene expression against abiotic stress conditions in plants. Moreover, the underlying mechanism linking up-regulation of tRFs under stress conditions with the stress tolerant response remains unknown. In this review, the biogenesis and putative function of diverse tRFs in abiotic stress signaling are discussed with a focus on tRFs as a transcriptional/post-transcriptional/translational regulator.

Identification of tRNA nucleoside modification genes critical for stress response and development in rice and Arabidopsis

BACKGROUND

Modification of nucleosides on transfer RNA (tRNA) is important either for correct mRNA decoding process or for tRNA structural stabilization. Nucleoside methylations catalyzed by MTase (methyltransferase) are the most common type among all tRNA nucleoside modifications. Although tRNA modified nucleosides and modification enzymes have been extensively studied in prokaryotic systems, similar research remains preliminary in higher plants, especially in crop species, such as rice (Oryza sativa). Rice is a monocot model plant as well as an important cereal crop, and stress tolerance and yield are of great importance for rice breeding.

RESULTS

In this study, we investigated how the composition and abundance of tRNA modified nucleosides could change in response to drought, salt and cold stress, as well as in different tissues during the whole growth season in two model plants-O. sativa and Arabidopsis thaliana. Twenty two and 20 MTase candidate genes were identified in rice and Arabidopsis, respectively, by protein sequence homology and conserved domain analysis. Four methylated nucleosides, Am, Cm, m¹A and m⁷G, were found to be very important in stress response both in rice and Arabidopsis. Additionally, three nucleosides,Gm, m⁵U and m⁵C, were involved in plant development. Hierarchical clustering analysis revealed consistency on Am, Cm, m¹A and m⁷G MTase candidate genes, and the abundance of the corresponding nucleoside under stress conditions. The same is true for Gm, m⁵U and m⁵C modifications and corresponding methylation genes in different tissues during different developmental stages.

CONCLUSIONS

We identified candidate genes for various tRNA modified nucleosides in rice and Arabidopsis, especially on MTases for methylated nucleosides. Based on bioinformatics analysis, nucleoside abundance assessments and gene expression profiling, we propose four methylated nucleosides (Am, Cm, m¹A and m⁷G) that are critical for stress response in rice and Arabidopsis, and three methylated nucleosides (Gm, m⁵U and m⁵C) that might be important during development.

5-Methylcytosine RNA Methylation in Arabidopsis Thaliana

5-Methylcytosine (m⁵C) is a well-characterized DNA modification, and is also predominantly reported in abundant non-coding RNAs in both prokaryotes and eukaryotes. However, the distribution and biological functions of m⁵C in plant mRNAs remain largely unknown. Here, we report transcriptome-wide profiling of RNA m⁵C in Arabidopsis thaliana by applying m⁵C RNA immunoprecipitation followed by a deep-sequencing approach (m⁵C-RIP-seq). LC-MS/MS and dot blot analyses reveal a dynamic pattern of m⁵C mRNA modification in various tissues and at different developmental stages. m⁵C-RIP-seq analysis identified 6045 m⁵C peaks in 4465 expressed genes in young seedlings. We found that m⁵C is enriched in coding sequences with two peaks located immediately after start codons and before stop codons, and is associated with mRNAs with low translation activity. We further demonstrated that an RNA (cytosine-5)-methyltransferase, tRNA-specific methyltransferase 4B (TRM4B), exhibits m⁵C RNA methyltransferase activity. Mutations in TRM4B display defects in root development and decreased m⁵C peaks. TRM4B affects the transcript levels of the genes involved in root development, which is positively correlated with their mRNA stability and m⁵C levels. Our results suggest that m⁵C in mRNA is a new epitranscriptome marker inArabidopsis, and that regulation of this modification is an integral part of gene regulatory networks underlying plant development.

The nuclear and organellar tRNA-derived RNA fragment population in Arabidopsis thaliana is highly dynamic

In the expanding repertoire of small noncoding RNAs (ncRNAs), tRNA-derived RNA fragments (tRFs) have been identified in all domains of life. Their existence in plants has been already proven but no detailed analysis has been performed. Here, short tRFs of 19–26 nucleotides were retrieved from Arabidopsis thaliana small RNA libraries obtained from various tissues, plants submitted to abiotic stress or fractions immunoprecipitated with ARGONAUTE 1 (AGO1). Large differences in the tRF populations of each extract were observed. Depending on the tRNA, either tRF-5D (due to a cleavage in the D region) or tRF-3T (via a cleavage in the T region) were found and hot spots of tRNA cleavages have been identified. Interestingly, up to 25% of the tRFs originate from plastid tRNAs and we provide evidence that mitochondrial tRNAs can also be a source of tRFs. Very specific tRF-5D deriving not only from nucleus-encoded but also from plastid-encoded tRNAs are strongly enriched in AGO1 immunoprecipitates. We demonstrate that the organellar tRFs are not found within chloroplasts or mitochondria but rather accumulate outside the organelles. These observations suggest that some organellar tRFs could play regulatory functions within the plant cell and may be part of a signaling pathway.

Transcriptome-Wide Mapping of RNA 5-Methylcytosine in Arabidopsis mRNAs and Noncoding RNAs

Posttranscriptional methylation of RNA cytosine residues to 5-methylcytosine (m5C) is an important modification with diverse roles, such as regulating stress responses, stem cell proliferation, and RNA metabolism. Here, we used RNA bisulfite sequencing for transcriptome-wide quantitative mapping of m5C in the model plant Arabidopsis thaliana We discovered more than a thousand m5C sites in Arabidopsis mRNAs, long noncoding RNAs, and other noncoding RNAs across three tissue types (siliques, seedling shoots, and roots) and validated a number of these sites. Quantitative differences in methylated sites between these three tissues suggest tissue-specific regulation of m5C. Perturbing the RNA m5C methyltransferase TRM4B resulted in the loss of m5C sites on mRNAs and noncoding RNAs and reduced the stability of tRNAAsp(GTC) We also demonstrate the importance of m5C in plant development, as trm4b mutants have shorter primary roots than the wild type due to reduced cell division in the root apical meristem. In addition, trm4b mutants show increased sensitivity to oxidative stress. Finally, we provide insights into the targeting mechanism of TRM4B by demonstrating that a 50-nucleotide sequence flanking m5C C3349 in MAIGO5 mRNA is sufficient to confer methylation of a transgene reporter in Nicotiana benthamiana.

Transcriptome-Wide Mapping of RNA 5-Methylcytosine in Arabidopsis mRNAs and Noncoding RNAs

Posttranscriptional methylation of RNA cytosine residues to 5-methylcytosine (m⁵C) is an important modification with diverse roles, such as regulating stress responses, stem cell proliferation, and RNA metabolism. Here, we used RNA bisulfite sequencing for transcriptome-wide quantitative mapping of m⁵C in the model plant Arabidopsis thaliana We discovered more than a thousand m⁵C sites in Arabidopsis mRNAs, long noncoding RNAs, and other noncoding RNAs across three tissue types (siliques, seedling shoots, and roots) and validated a number of these sites. Quantitative differences in methylated sites between these three tissues suggest tissue-specific regulation of m⁵C. Perturbing the RNA m⁵C methyltransferase TRM4B resulted in the loss of m⁵C sites on mRNAs and noncoding RNAs and reduced the stability of tRNA^Asp(GTC) We also demonstrate the importance of m⁵C in plant development, as trm4b mutants have shorter primary roots than the wild type due to reduced cell division in the root apical meristem. In addition, trm4b mutants show increased sensitivity to oxidative stress. Finally, we provide insights into the targeting mechanism of TRM4B by demonstrating that a 50-nucleotide sequence flanking m⁵C C3349 in MAIGO5 mRNA is sufficient to confer methylation of a transgene reporter in Nicotiana benthamiana.

Plant DNA Methyltransferase Genes: Multiplicity, Expression, Methylation Patterns

Expression and methylation patterns of genes encoding DNA methyltransferases and their functionally related proteins were studied in organs of Arabidopsis thaliana plants. Genes coding for the major maintenance-type DNA methyltransferases, MET1 and CMT3, and the major de novo-type DNA methyltransferase, DRM2, are actively expressed in all organs. Similar constitutively active expression was observed for genes encoding their functionally related proteins, a histone H3K9 methyltransferase KYP and a catalytically non-active protein DRM3. Expression of the MET1 and CMT3 genes is significantly lower in developing endosperm compared with embryo. Vice versa, expression of the MET2a, MET2b, MET3, and CMT2 genes in endosperm is much more active compared with embryo. A special maintenance DNA methylation system seems to operate in endosperm. The DNMT2 and N6AMT genes encoding putative methyltransferases are constitutively expressed at low levels. CMT1 and DRM1 genes are expressed rather weakly in all investigated organs. Most of the studied genes have methylation patterns conforming to the "body-methylated gene" prototype. A peculiar feature of the MET family genes is methylation at all three possible site types (CG, CHG, and CHH). The most weakly expressed among genes of their respective families, CMT1 and DRM1, are practically unmethylated. The MET3 and N6AMT genes have unusual methylation patterns, promoter region, and most of the gene body devoid of any methylation, and the 3'-end proximal part of the gene body is highly methylated.

Deciphering the epitranscriptome: A green perspective

The advent of high-throughput sequencing technologies coupled with new detection methods of RNA modifications has enabled investigation of a new layer of gene regulation - the epitranscriptome. With over 100 known RNA modifications, understanding the repertoire of RNA modifications is a huge undertaking. This review summarizes what is known about RNA modifications with an emphasis on discoveries in plants. RNA ribose modifications, base methylations and pseudouridylation are required for normal development in Arabidopsis, as mutations in the enzymes modifying them have diverse effects on plant development and stress responses. These modifications can regulate RNA structure, turnover and translation. Transfer RNA and ribosomal RNA modifications have been mapped extensively and their functions investigated in many organisms, including plants. Recent work exploring the locations, functions and targeting of N⁶ -methyladenosine (m⁶ A), 5-methylcytosine (m⁵ C), pseudouridine (Ψ), and additional modifications in mRNAs and ncRNAs are highlighted, as well as those previously known on tRNAs and rRNAs. Many questions remain as to the exact mechanisms of targeting and functions of specific modified sites and whether these modifications have distinct functions in the different classes of RNAs.