The Arabidopsis epitranscriptome

N6-methyladenosine (m6A) modification: Emerging regulators in plant-virus interactions

N6-methyladenosine (m6A), a reversible epigenetic modification, is widely present on both cellular and viral RNAs. This modification undergoes catalysis by methyltransferases (writers), removal by demethylases (erasers), and recognition by m6A-binding proteins (readers), ultimately influencing the fate and function of modified RNA molecules. With recent advances in sequencing technologies, the genome-wide mapping of m6A has become possible, enabling a deeper exploration of its roles during viral infections. So far, while the significance of m6A in regulating virus-host interactions has been well-established in animal viruses, research on its involvement in plant viruses remains in its early stages. In this review, we summarize the current knowledge regarding the functions and molecular mechanisms of m6A in plant-virus interactions. A better understanding of these complex interactions may provide valuable insights for developing novel antiviral strategies, potentially leading to more effective control of plant viral diseases in the field.

ALBA proteins facilitate cytoplasmic YTHDF-mediated reading of m6A in Arabidopsis

N6-methyladenosine (m6A) exerts many of its regulatory effects on eukaryotic mRNAs by recruiting cytoplasmic YT521-B homology-domain family (YTHDF) proteins. Here, we show that in Arabidopsis thaliana, the interaction between m6A and the major YTHDF protein ECT2 also involves the mRNA-binding ALBA protein family. ALBA and YTHDF proteins physically associate via a deeply conserved short linear motif in the intrinsically disordered region of YTHDF proteins and their mRNA target sets overlap, with ALBA4 binding sites being juxtaposed to m6A sites. These binding sites correspond to pyrimidine-rich elements previously found to be important for m6A binding to ECT2. Accordingly, both the biological functions of ECT2, and its binding to m6A targets in vivo, require ALBA association. Our results introduce the YTHDF-ALBA complex as the functional cytoplasmic m6A-reader in Arabidopsis, and define a molecular foundation for the concept of facilitated m6A reading, which increases the potential for combinatorial control of biological m6A effects.

The m6A reader SlYTH2 negatively regulates tomato fruit aroma by impeding the translation process

N6-methyladenosine (m6A) is a fundamentally important RNA modification for gene regulation, whose function is achieved through m6A readers. However, whether and how m6A readers play regulatory roles during fruit ripening and quality formation remains unclear. Here, we characterized SlYTH2 as a tomato m6A reader protein and profiled the binding sites of SlYTH2 at the transcriptome-wide level. SlYTH2 undergoes liquid-liquid phase separation and promotes RNA-protein condensate formation. The target mRNAs of SlYTH2, namely m6A-modified SlHPL and SlCCD1B associated with volatile synthesis, are enriched in SlYTH2-induced condensates. Through polysome profiling assays and proteomic analysis, we demonstrate that knockout of SlYTH2 expedites the translation process of SlHPL and SlCCD1B, resulting in augmented production of aroma-associated volatiles. This aroma enrichment significantly increased consumer preferences for CRISPR-edited fruit over wild type. These findings shed light on the underlying mechanisms of m6A in plant RNA metabolism and provided a promising strategy to generate fruits that are more attractive to consumers.

MeRIP-Seq for Identifying Stress-Responsive Transcriptome-Wide m6A Profiles in Plants

Recent advancements in detection and mapping methods have enabled researchers to uncover the biological importance of RNA chemical modifications, which play a vital role in post-transcriptional gene regulation. Although numerous types of RNA modifications have been identified in higher eukaryotes, only a few have been extensively studied for their biological functions. Of these, N6-methyladenosine (m6A) is the most prevalent and important mRNA modification that influences various aspects of RNA metabolism, including mRNA stability, degradation, splicing, alternative polyadenylation, export, and localization, as well as translation. Thus, they have implications for a variety of biological processes, including growth, development, and stress responses. The m6A deposition or removal on transcripts is dynamic and is altered in response to internal and external cues. Because this mark can alter gene expression under stress conditions, it is essential to identify the transcripts that can acquire or lose this epitranscriptomic mark upon exposure to stress conditions. Here we describe a step-by-step protocol for identifying stress-responsive transcriptome-wide m6A changes using RNA immunoprecipitation followed by high-throughput sequencing (MeRIP-seq).

A YTHDF-PABP interaction is required for m6 A-mediated organogenesis in plants

N6-methyladenosine (m6 A) in mRNA is key to eukaryotic gene regulation. Many m6 A functions involve RNA-binding proteins that recognize m6 A via a YT521-B Homology (YTH) domain. YTH domain proteins contain long intrinsically disordered regions (IDRs) that may mediate phase separation and interaction with protein partners, but whose precise biochemical functions remain largely unknown. The Arabidopsis thaliana YTH domain proteins ECT2, ECT3, and ECT4 accelerate organogenesis through stimulation of cell division in organ primordia. Here, we use ECT2 to reveal molecular underpinnings of this function. We show that stimulation of leaf formation requires the long N-terminal IDR, and we identify two short IDR elements required for ECT2-mediated organogenesis. Of these two, a 19-amino acid region containing a tyrosine-rich motif conserved in both plant and metazoan YTHDF proteins is necessary for binding to the major cytoplasmic poly(A)-binding proteins PAB2, PAB4, and PAB8. Remarkably, overexpression of PAB4 in leaf primordia partially rescues the delayed leaf formation in ect2 ect3 ect4 mutants, suggesting that the ECT2-PAB2/4/8 interaction on target mRNAs of organogenesis-related genes may overcome limiting PAB concentrations in primordial cells.

Insights into the conservation and diversification of the molecular functions of YTHDF proteins

YT521-B homology (YTH) domain proteins act as readers of N6-methyladenosine (m6A) in mRNA. Members of the YTHDF clade determine properties of m6A-containing mRNAs in the cytoplasm. Vertebrates encode three YTHDF proteins whose possible functional specialization is debated. In land plants, the YTHDF clade has expanded from one member in basal lineages to eleven so-called EVOLUTIONARILY CONSERVED C-TERMINAL REGION1-11 (ECT1-11) proteins in Arabidopsis thaliana, named after the conserved YTH domain placed behind a long N-terminal intrinsically disordered region (IDR). ECT2, ECT3 and ECT4 show genetic redundancy in stimulation of primed stem cell division, but the origin and implications of YTHDF expansion in higher plants are unknown, as it is unclear whether it involves acquisition of fundamentally different molecular properties, in particular of their divergent IDRs. Here, we use functional complementation of ect2/ect3/ect4 mutants to test whether different YTHDF proteins can perform the same function when similarly expressed in leaf primordia. We show that stimulation of primordial cell division relies on an ancestral molecular function of the m6A-YTHDF axis in land plants that is present in bryophytes and is conserved over YTHDF diversification, as it appears in all major clades of YTHDF proteins in flowering plants. Importantly, although our results indicate that the YTH domains of all arabidopsis ECT proteins have m6A-binding capacity, lineage-specific neo-functionalization of ECT1, ECT9 and ECT11 happened after late duplication events, and involves altered properties of both the YTH domains, and, especially, of the IDRs. We also identify two biophysical properties recurrent in IDRs of YTHDF proteins able to complement ect2 ect3 ect4 mutants, a clear phase separation propensity and a charge distribution that creates electric dipoles. Human and fly YTHDFs do not have IDRs with this combination of properties and cannot replace ECT2/3/4 function in arabidopsis, perhaps suggesting different molecular activities of YTHDF proteins between major taxa.

Plant YTHDF proteins are direct effectors of antiviral immunity against an N6-methyladenosine-containing RNA virus

In virus-host interactions, nucleic acid-directed first lines of defense that allow viral clearance without compromising growth are of paramount importance. Plants use the RNA interference pathway as a basal antiviral immune system, but additional RNA-based mechanisms of defense also exist. The infectivity of a plant positive-strand RNA virus, alfalfa mosaic virus (AMV), relies on the demethylation of viral RNA by the recruitment of the cellular N6-methyladenosine (m6 A) demethylase ALKBH9B, but how demethylation of viral RNA promotes AMV infection remains unknown. Here, we show that inactivation of the Arabidopsis cytoplasmic YT521-B homology domain (YTH)-containing m6 A-binding proteins ECT2, ECT3, and ECT5 is sufficient to restore AMV infectivity in partially resistant alkbh9b mutants. We further show that the antiviral function of ECT2 is distinct from its previously demonstrated function in the promotion of primordial cell proliferation: an ect2 mutant carrying a small deletion in its intrinsically disordered region is partially compromised for antiviral defense but not for developmental functions. These results indicate that the m6 A-YTHDF axis constitutes a novel branch of basal antiviral immunity in plants.

The covalent nucleotide modifications within plant mRNAs: What we know, how we find them, and what should be done in the future

Although covalent nucleotide modifications were first identified on the bases of transfer RNAs (tRNAs) and ribosomal RNAs (rRNAs), a number of these epitranscriptome marks have also been found to occur on the bases of messenger RNAs (mRNAs). These covalent mRNA features have been demonstrated to have various and significant effects on the processing (e.g. splicing, polyadenylation, etc.) and functionality (e.g. translation, transport, etc.) of these protein-encoding molecules. Here, we focus our attention on the current understanding of the collection of covalent nucleotide modifications known to occur on mRNAs in plants, how they are detected and studied, and the most outstanding future questions of each of these important epitranscriptomic regulatory signals.

Clinical and molecular significance of the RNA m6A methyltransferase complex in prostate cancer

N⁶-methyladenosine (m⁶A) is the most abundant internal mRNA modification and is dynamically regulated through distinct protein complexes that methylate, demethylate, and/or interpret the m⁶A modification. These proteins, and the m⁶A modification, are involved in the regulation of gene expression, RNA stability, splicing and translation. Given its role in these crucial processes, m⁶A has been implicated in many diseases, including in cancer development and progression. Prostate cancer (PCa) is the most commonly diagnosed non-cutaneous cancer in men and recent studies support a role for m⁶A in PCa. Despite this, the literature currently lacks an integrated analysis of the expression of key components of the m⁶A RNA methyltransferase complex, both in PCa patients and in well-established cell line models. For this reason, this study used immunohistochemistry and functional studies to investigate the mechanistic and clinical significance of the METTL3, METTL14, WTAP and CBLL1 components of the m⁶A methyltransferase complex in PCa specimens and cell lines. Expression of METTL3 and CBLL1, but not METTL14 and WTAP, was associated with poorer PCa patient outcomes. Expression of METTL3, METTL14, WTAP and CBLL1 was higher in PCa cells compared with non-malignant prostate cells, with the highest expression seen in castrate-sensitive, androgen-responsive PCa cells. Moreover, in PCa cell lines, expression of METTL3 and WTAP was found to be androgen-regulated. To investigate the mechanistic role(s) of the m⁶A methyltransferase complex in PCa cells, short hairpin RNA (shRNA)-mediated knockdown coupled with next generation sequencing was used to determine the transcriptome-wide roles of METTL3, the catalytic subunit of the m⁶A methyltransferase complex. Functional depletion of METTL3 resulted in upregulation of the androgen receptor (AR), together with 134 AR-regulated genes. METTL3 knockdown also resulted in altered splicing, and enrichment of cell cycle, DNA repair and metabolic pathways. Collectively, this study identified the functional and clinical significance of four essential m⁶A complex components in PCa patient specimens and cell lines for the first time. Further studies are now warranted to determine the potential therapeutic relevance of METTL3 inhibitors in development to treat leukaemia to benefit patients with PCa.

Epitranscriptome changes triggered by ammonium nutrition regulate the proteome response of maritime pine roots

Epitranscriptome constitutes a gene expression checkpoint in all living organisms. Nitrogen is an essential element for plant growth and development that influences gene expression at different levels such as epigenome, transcriptome, proteome, and metabolome. Therefore, our hypothesis is that changes in the epitranscriptome may regulate nitrogen metabolism. In this study, epitranscriptomic modifications caused by ammonium nutrition were monitored in maritime pine roots using Oxford Nanopore Technology. Transcriptomic responses mainly affected transcripts involved in nitrogen and carbon metabolism, defense, hormone synthesis/signaling, and translation. Global detection of epitranscriptomic marks was performed to evaluate this posttranscriptional mechanism in un/treated seedlings. Increased N⁶-methyladenosine (m⁶A) deposition in the 3'-UTR was observed in response to ammonium, which seems to be correlated with poly(A) lengths and changes in the relative abundance of the corresponding proteins. The results showed that m⁶A deposition and its dynamics seem to be important regulators of translation under ammonium nutrition. These findings suggest that protein translation is finely regulated through epitranscriptomic marks likely by changes in mRNA poly(A) length, transcript abundance and ribosome protein composition. An integration of multiomics data suggests that the epitranscriptome modulates responses to nutritional, developmental and environmental changes through buffering, filtering, and focusing the final products of gene expression.

Tomato SlYTH1 encoding a putative RNA m6A reader affects plant growth and fruit shape

N⁶-methyladenosine (m⁶A), the most abundant and common modification on eukaryotic mRNA, plays crucial roles in multiple biological processes through controlling endogenous gene activity in organisms. The m⁶A reader specifically recognizes the m⁶A mark to mediate the regulation of m⁶A on mRNA, and determines the fate of its target mRNA. In plants, the currently confirmed m⁶A readers are YTH (YT521B homology) domain-containing proteins. We previously reported that tomato contains 9 YTH genes, of which SlYTH1 has the strongest expression. The present study reports the functional characterization of SlYTH1 in tomato. SlYTH1 mutants generated via CRISPR/Cas9 technology exhibited pleiotropic phenotypes, including low seed germination rate, shortened seedling root, retarded plant growth and development during vegetative development, and elongated and longitudinally flattened fruit with reduced the locule number. SlYTH1 knockout reduced GA₃ content and downregulated the expression of related genes in gibberellin biosynthesis pathway. Moreover, exogenous GA₃ application could partially restore the phenotypic defects caused by SlYTH1 mutations. SlYTH1 knockout could alleviate the inhibition of seedling root elongation by exogenous GA₃ application at relatively low concentration. These facts indicated SlYTH1 is involved in regulating gibberellin biosynthesis and plays important roles in multiple physiological processes in tomato.

Covalent RNA modifications and their budding crosstalk with plant epigenetic processes

Our recent cognizance of diverse RNA classes undergoing dynamic covalent chemical modifications (or epitranscriptomic marks) in plants has provided fresh insight into the underlying molecular mechanisms of gene expression regulation. Comparatively, epigenetic marks comprising heritable modifications of DNA and histones have been extensively studied in plants and their impact on plant gene expression is quite established. Based on our growing knowledge of the plant epitranscriptome and epigenome, it is logical to explore how the two regulatory layers intermingle to intricately determine gene expression levels underlying key biological processes such as development and response to stress. Herein, we focus on the emerging evidence of crosstalk between the plant epitranscriptome with epigenetic regulation involving DNA modification, histone modification, and non-coding RNAs.

Transcriptome-wide N6-methyladenosine (m6A) methylation in soybean under Meloidogyne incognita infection

N6-methyladenosine (m6A) is a reversible epigenetic modification of mRNA and other RNAs that plays a significant role in regulating gene expression and biological processes. However, m6A abundance, dynamics, and transcriptional regulatory mechanisms remain unexplored in the context of soybean resistance to Meloidogyne incognita. In this study, we performed a comparative analysis of transcriptome-wide m6A and metabolome profiles of soybean root tissues with and without M. incognita infection. Global m6A hypermethylation was widely induced in response to M. incognita infection and was enriched around the 3' end of coding sequences and in 3' UTR regions. There were 2069 significantly modified m6A sites, 594 differentially expressed genes, and 103 differentially accumulated metabolites between infected and uninfected roots, including coumestrol, psoralidin, and 2-hydroxyethylphosphonate. Among 101 m6A-modified DEGs, 34 genes were hypomethylated and upregulated, and 39 genes were hypermethylated and downregulated, indicating a highly negative correlation between m6A methylation and gene transcript abundance. A number of these m6A-modified DEGs, including WRKY70, ERF60, POD47 and LRR receptor-like serine/threonine-protein kinases, were involved in plant defense responses. Our study provides new insights into the critical role of m6A modification in early soybean responses to M. incognita.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-022-00077-2.

Transcriptional expression of m6A and m5C RNA methyltransferase genes in the brain and fat body of honey bee adult workers

Honey bee (Apis mellifera) adult workers change behaviors and nutrition according to age progression. Young workers, such as nurses, perform in-hive tasks and consume protein-rich pollen, while older workers (foragers) leave the colony to search for food, and consume carbohydrate-rich nectar. These environmentally stimulated events involve transcriptional and DNA epigenetic marks alterations in worker tissues. However, post-transcriptional RNA modifications (epitranscriptomics) are still poorly explored in bees. We investigated the transcriptional profiles of m6A and m5C RNA methyltransferases in the brain and fat body of adult workers of 1) different ages and performing different tasks [nurses of 8 days-old (N-8D) and foragers of 29 days-old (F-29D), sampled from wild-type colonies], and 2) same-aged young workers caged in an incubator and treated with a pollen-rich [PR] or a pollen-deprived [PD] diet for 8 days. In the brain, METTL3, DNMT2, NOP2, NSUN2, NSUN5, and NSUN7 genes increased expression during adulthood (from N-8D to F-29D), while the opposite pattern was observed in the fat body for METTL3, DNMT2, and NSUN2 genes. Regarding diet treatments, high expression levels were observed in the brains of the pollen-deprived group (DNMT2, NOP2, and NSUN2 genes) and the fat bodies of the pollen-rich group (NOP2, NSUN4, and NSUN5 genes) compared to the brains of the PR group and the fat bodies of the PD group, respectively. Our data indicate that RNA epigenetics may be an important regulatory layer in the development of adult workers, presenting tissue-specific signatures of RNA methyltransferases expression in response to age, behavior, and diet content.

mRNA N6 -methyladenosine is critical for cold tolerance in Arabidopsis

Plants respond to low temperatures by altering the mRNA abundance of thousands of genes contributing to numerous physiological and metabolic processes that allow them to adapt. At the post-transcriptional level, these cold stress-responsive transcripts undergo alternative splicing, microRNA-mediated regulation and alternative polyadenylation, amongst others. Recently, m⁶ A, m⁵ C and other mRNA modifications that can affect the regulation and stability of RNA were discovered, thus revealing another layer of post-transcriptional regulation that plays an important role in modulating gene expression. The importance of m⁶ A in plant growth and development has been appreciated, although its significance under stress conditions is still underexplored. To assess the role of m⁶ A modifications during cold stress responses, methylated RNA immunoprecipitation sequencing was performed in Arabidopsis seedlings esposed to low temperature stress (4°C) for 24 h. This transcriptome-wide m⁶ A analysis revealed large-scale shifts in this modification in response to low temperature stress. Because m⁶ A is known to affect transcript stability/degradation and translation, we investigated these possibilities. Interestingly, we found that cold-enriched m⁶ A-containing transcripts demonstrated the largest increases in transcript abundance coupled with increased ribosome occupancy under cold stress. The significance of the m⁶ A epitranscriptome on plant cold tolerance was further assessed using the mta mutant in which the major m⁶ A methyltransferase gene was mutated. Compared to the wild-type, along with the differences in CBFs and COR gene expression levels, the mta mutant exhibited hypersensitivity to cold treatment as determined by primary root growth, biomass, and reactive oxygen species accumulation. Furthermore, and most importantly, both non-acclimated and cold-acclimated mta mutant demonstrated hypersensitivity to freezing tolerance. Taken together, these findings suggest a critical role for the epitranscriptome in cold tolerance of Arabidopsis.

Active genic machinery for epigenetic RNA modifications in bees

Epitranscriptomics is an emerging field of investigation dedicated to the study of post-transcriptional RNA modifications. RNA methylations regulate RNA metabolism and processing, including changes in response to environmental cues. Although RNA modifications are conserved from bacteria to eukaryotes, there is little evidence of an epitranscriptomic pathway in insects. Here we identified genes related to RNA m⁶ A (N6-methyladenine) and m⁵ C (5-methylcytosine) methylation machinery in seven bee genomes (Apis mellifera, Melipona quadrifasciata, Frieseomelitta varia, Eufriesea mexicana, Bombus terrestris, Megachile rotundata and Dufourea novaeangliae). In A. mellifera, we validated the expression of methyltransferase genes and found that the global levels of m⁶ A and m⁵ C measured in the fat body and brain of adult workers differ significantly. Also, m⁶ A levels were differed significantly mainly between the fourth larval instar of queens and workers. Moreover, we found a conserved m⁵ C site in the honeybee 28S rRNA. Taken together, we confirm the existence of epitranscriptomic machinery acting in bees and open avenues for future investigations on RNA epigenetics in a wide spectrum of hymenopteran species.

Genome-wide identification and characterization of YTH domain-containing genes, encoding the m6A readers, and their expression in tomato

9 YTH genes in tomato were identified and cloned, and their expression patterns were comprehensively analyzed, which reveal potential multiple roles in development and fruit ripening. N⁶-methyladenosine (m⁶A) is an abundant and pervasive post-transcriptional modification in eukaryotic mRNAs. The YTH domain-containing proteins act as m⁶A readers to read m⁶A marks and transduce their downstream regulatory effects by altering m⁶A-mRNA metabolism processes. Identification of YTH proteins is essential for understanding the regulatory mechanisms of m⁶A in physiological processes, but little is known about YTH proteins in tomato, a model system for fruit development. Here, we report that tomato genomes contain a total of 9 SlYTH genes. While YTH proteins of both tomato and Arabidopsis can be classified into two subfamilies, the member distributions in subfamilies are very different between the two species. Homology modeling exhibited the similar three-dimensional structures of SlYTH proteins to human YTHDF1 or YTHDC1. Multiple hormone-response elements locating on the promoters of SlYTH genes indicate that they are involved in the physiological processes related to phytohormone. SlYTH genes are ubiquitous and spatiotemporal dynamic expression in tomato. Eight SlTYH genes have the strongest expression in stamens among the parts of flowers. Throughout fruit ontogeny, most of the SlYTH genes display obvious high mRNA levels during the developmental phases (4 dpa to mature green); moreover, SlYTH1 and SlYTH2 have absolute predominant expressions demonstrated by RNA-seq. The results lay a foundation for future characterizations on the functions of YTH proteins and m⁶A regulatory mechanism in tomato.

R-loop resolution promotes co-transcriptional chromatin silencing

RNA-mediated chromatin silencing is central to genome regulation in many organisms. However, how nascent non-coding transcripts regulate chromatin is poorly understood. Here, through analysis of Arabidopsis FLC, we show that resolution of a nascent-transcript-induced R-loop promotes chromatin silencing. Stabilization of an antisense-induced R-loop at the 3' end of FLC enables an RNA binding protein FCA, with its direct partner FY/WDR33 and other 3'-end processing factors, to polyadenylate the nascent antisense transcript. This clears the R-loop and recruits the chromatin modifiers demethylating H3K4me1. FCA immunoprecipitates with components of the m⁶A writer complex, and m⁶A modification affects dynamics of FCA nuclear condensates, and promotes FLC chromatin silencing. This mechanism also targets other loci in the Arabidopsis genome, and consistent with this fca and fy are hypersensitive to a DNA damage-inducing drug. These results show how modulation of R-loop stability by co-transcriptional RNA processing can trigger chromatin silencing.

M6A RNA Methylation Regulator HNRNPC Contributes to Tumorigenesis and Predicts Prognosis in Glioblastoma Multiforme

Glioblastoma multiforme (GBM) is the most malignant glioma with a high death rate. N6-methyladenosine (m6A) RNA methylation plays an increasingly important role in tumors. The current study aimed to determine the function of the regulators of m6A RNA methylation in GBM. We evaluated the difference, interaction, and correlation of these regulators with TCGA database. HNRNPC, WTAP, YTHDF2 and, YTHDF1 were significantly upregulated in GBM. To explore the expression characteristics of regulators in GBM, we defined two subgroups through consensus cluster. HNRNPC, WTAP, and YTHDF2 were significantly upregulated in the cluster2 which had a good overall survival (OS). To investigate the prognostic value of regulators, we used lasso cox regression algorithm to screen an independent prognostic risk characteristic based on the expression of HNRNPC, ZC3H13, and YTHDF2. The prognostic feature between the low and high-risk groups was significantly different (P < 0.05), which could predict significance of prognosis (area under the curve (AUC) = 0.819). Moreover, we used western blot, RT-PCR, and immunohistochemical staining to verify the expression of HNRNPC was associated with malignancy and development of gliomas. Similarly, the high expression of HNRNPC had a good prognosis. In conclusion, HNRNPC is a vital participant in the malignant progression of GBM and might be valuable for prognosis.

Recurrent requirement for the m6A-ECT2/ECT3/ECT4 axis in the control of cell proliferation during plant organogenesis

mRNA methylation at the N6-position of adenosine (m6A) enables multiple layers of post-transcriptional gene control, often via RNA-binding proteins that use a YT521-B homology (YTH) domain for specific m6A recognition. In Arabidopsis, normal leaf morphogenesis and rate of leaf formation require m6A and the YTH-domain proteins ECT2, ECT3 and ECT4. In this study, we show that ect2/ect3 and ect2/ect3/ect4 mutants also exhibit slow root and stem growth, slow flower formation, defective directionality of root growth, and aberrant flower and fruit morphology. In all cases, the m6A-binding site of ECT proteins is required for in vivo function. We also demonstrate that both m6A methyltransferase mutants and ect2/ect3/ect4 exhibit aberrant floral phyllotaxis. Consistent with the delayed organogenesis phenotypes, we observe particularly high expression of ECT2, ECT3 and ECT4 in rapidly dividing cells of organ primordia. Accordingly, ect2/ect3/ect4 mutants exhibit decreased rates of cell division in leaf and vascular primordia. Thus, the m6A-ECT2/ECT3/ECT4 axis is employed as a recurrent module to stimulate plant organogenesis, at least in part by enabling rapid cellular proliferation.

Marking RNA: m6A writers, readers, and functions in Arabidopsis

N6-methyladenosine (m6A) emerges as an important modification in eukaryotic mRNAs. m6A has first been reported in 1974, and its functional significance in mammalian gene regulation and importance for proper development have been well established. An arsenal of writer, eraser, and reader proteins accomplish deposition, removal, and interpretation of the m6A mark, resulting in dynamic function. This led to the concept of an epitranscriptome, the compendium of RNA species with chemical modification of the nucleobases in the cell, in analogy to the epigenome. While m6A has long been known to also exist in plant mRNAs, proteins involved in m6A metabolism have only recently been detected by mutant analysis, homology search, and mRNA interactome capture in the reference plant Arabidopsis thaliana. Dysregulation of the m6A modification causes severe developmental abnormalities of leaves and roots and altered timing of reproductive development. Furthermore, m6A modification affects viral infection. Here, we discuss recent progress in identifying m6A sites transcriptome-wide, in identifying the molecular players involved in writing, removing, and reading the mark, and in assigning functions to this RNA modification in A. thaliana. We highlight similarities and differences to m6A modification in mammals and provide an outlook on important questions that remain to be addressed.

Reactive Oxygen Species (ROS) and Nucleic Acid Modifications During Seed Dormancy

The seed is the propagule of higher plants and allows its dissemination and the survival of the species. Seed dormancy prevents premature germination under favourable conditions. Dormant seeds are only able to germinate in a narrow range of conditions. During after-ripening (AR), a mechanism of dormancy release, seeds gradually lose dormancy through a period of dry storage. This review is mainly focused on how chemical modifications of mRNA and genomic DNA, such as oxidation and methylation, affect gene expression during late stages of seed development, especially during dormancy. The oxidation of specific nucleotides produced by reactive oxygen species (ROS) alters the stability of the seed stored mRNAs, being finally degraded or translated into non-functional proteins. DNA methylation is a well-known epigenetic mechanism of controlling gene expression. In Arabidopsis thaliana, while there is a global increase in CHH-context methylation through embryogenesis, global DNA methylation levels remain stable during seed dormancy, decreasing when germination occurs. The biological significance of nucleic acid oxidation and methylation upon seed development is discussed.

Multiple links between 5-methylcytosine content of mRNA and translation

BACKGROUND

5-Methylcytosine (m5C) is a prevalent base modification in tRNA and rRNA but it also occurs more broadly in the transcriptome, including in mRNA, where it serves incompletely understood molecular functions. In pursuit of potential links of m5C with mRNA translation, we performed polysome profiling of human HeLa cell lysates and subjected RNA from resultant fractions to efficient bisulfite conversion followed by RNA sequencing (bsRNA-seq). Bioinformatic filters for rigorous site calling were devised to reduce technical noise.

RESULTS

We obtained ~ 1000 candidate m5C sites in the wider transcriptome, most of which were found in mRNA. Multiple novel sites were validated by amplicon-specific bsRNA-seq in independent samples of either human HeLa, LNCaP and PrEC cells. Furthermore, RNAi-mediated depletion of either the NSUN2 or TRDMT1 m5C:RNA methyltransferases showed a clear dependence on NSUN2 for the majority of tested sites in both mRNAs and noncoding RNAs. Candidate m5C sites in mRNAs are enriched in 5'UTRs and near start codons and are embedded in a local context reminiscent of the NSUN2-dependent m5C sites found in the variable loop of tRNA. Analysing mRNA sites across the polysome profile revealed that modification levels, at bulk and for many individual sites, were inversely correlated with ribosome association.

CONCLUSIONS

Our findings emphasise the major role of NSUN2 in placing the m5C mark transcriptome-wide. We further present evidence that substantiates a functional interdependence of cytosine methylation level with mRNA translation. Additionally, we identify several compelling candidate sites for future mechanistic analysis.

Epigenetic Modifications of mRNA and DNA in Plants

Advances in the detection and mapping of messenger RNA (mRNA) N⁶-methyladenosine (m⁶A) and 5-methylcytosine (m⁵C), and DNA N⁶-methyldeoxyadenosine (6mA) redefined our understanding of these modifications as additional tiers of epigenetic regulation. In plants, the most prevalent internal mRNA modifications, m⁶A and m⁵C, play crucial and dynamic roles in many processes, including embryo development, stem cell fate determination, trichome branching, leaf morphogenesis, floral transition, stress responses, fruit ripening, and root development. The newly identified and widespread epigenetic marker 6mA DNA methylation is associated with gene expression, plant development, and stress responses. Here, we review the latest research progress on mRNA and DNA epigenetic modifications, including the detection, dynamics, distribution, functions, regulatory proteins, and evolution, with a focus on m⁶A, m⁵C, and 6mA. We also provide some perspectives on future research of the newly identified and unknown epigenetic modifications of mRNA and DNA in plants.

Epigenetic Methylations on N6-Adenine and N6-Adenosine with the same Input but Different Output

Epigenetic modifications on individual bases in DNA and RNA can encode inheritable genetic information beyond the canonical bases. Among the nucleic acid modifications, DNA N6-methadenine (6mA) and RNA N6-methyladenosine (m6A) have recently been well-studied due to the technological development of detection strategies and the functional identification of modification enzymes. The current findings demonstrate a wide spectrum of 6mA and m6A distributions from prokaryotes to eukaryotes and critical roles in multiple cellular processes. It is interesting that the processes of modification in which the methyl group is added to adenine and adenosine are the same, but the outcomes of these modifications in terms of their physiological impacts in organisms are quite different. In this review, we summarize the latest progress in the study of enzymes involved in the 6mA and m6A methylation machinery, including methyltransferases and demethylases, and their functions in various biological pathways. In particular, we focus on the mechanisms by which 6mA and m6A regulate the expression of target genes, and we highlight the future challenges in epigenetic regulation.

Analysis of Transcriptome and Epitranscriptome in Plants Using PacBio Iso-Seq and Nanopore-Based Direct RNA Sequencing

Nanopore sequencing from Oxford Nanopore Technologies (ONT) and Pacific BioSciences (PacBio) single-molecule real-time (SMRT) long-read isoform sequencing (Iso-Seq) are revolutionizing the way transcriptomes are analyzed. These methods offer many advantages over most widely used high-throughput short-read RNA sequencing (RNA-Seq) approaches and allow a comprehensive analysis of transcriptomes in identifying full-length splice isoforms and several other post-transcriptional events. In addition, direct RNA-Seq provides valuable information about RNA modifications, which are lost during the PCR amplification step in other methods. Here, we present a comprehensive summary of important applications of these technologies in plants, including identification of complex alternative splicing (AS), full-length splice variants, fusion transcripts, and alternative polyadenylation (APA) events. Furthermore, we discuss the impact of the newly developed nanopore direct RNA-Seq in advancing epitranscriptome research in plants. Additionally, we summarize computational tools for identifying and quantifying full-length isoforms and other co/post-transcriptional events and discussed some of the limitations with these methods. Sequencing of transcriptomes using these new single-molecule long-read methods will unravel many aspects of transcriptome complexity in unprecedented ways as compared to previous short-read sequencing approaches. Analysis of plant transcriptomes with these new powerful methods that require minimum sample processing is likely to become the norm and is expected to uncover novel co/post-transcriptional gene regulatory mechanisms that control biological outcomes during plant development and in response to various stresses.

Chemical Modifications in the Life of an mRNA Transcript

Investigations over the past eight years of chemical modifications on messenger RNA (mRNA) have revealed a new level of posttranscriptional gene regulation in eukaryotes. Rapid progress in our understanding of these modifications, particularly, N6-methyladenosine (m6A), has revealed their roles throughout the life cycle of an mRNA transcript. m6A methylation provides a rapid mechanism for coordinated transcriptome processing and turnover that is important in embryonic development and cell differentiation. In response to cellular signals, m6A can also regulate the translation of specific pools of transcripts. These mechanisms can be hijacked in human diseases, including numerous cancers and viral infection. Beyond m6A, many other mRNA modifications have been mapped in the transcriptome, but much less is known about their biological functions. As methods continue to be developed, we will be able to study these modifications both more broadly and in greater depth, which will likely reveal a wealth of new RNA biology.

Identification of factors required for m6 A mRNA methylation in Arabidopsis reveals a role for the conserved E3 ubiquitin ligase HAKAI

N6-adenosine methylation (m⁶ A) of mRNA is an essential process in most eukaryotes, but its role and the status of factors accompanying this modification are still poorly understood. Using combined methods of genetics, proteomics and RNA biochemistry, we identified a core set of mRNA m⁶ A writer proteins in Arabidopsis thaliana. The components required for m⁶ A in Arabidopsis included MTA, MTB, FIP37, VIRILIZER and the E3 ubiquitin ligase HAKAI. Downregulation of these proteins led to reduced relative m⁶ A levels and shared pleiotropic phenotypes, which included aberrant vascular formation in the root, indicating that correct m⁶ A methylation plays a role in developmental decisions during pattern formation. The conservation of these proteins amongst eukaryotes and the demonstration of a role in writing m⁶ A for the E3 ubiquitin ligase HAKAI is likely to be of considerable relevance beyond the plant sciences.

RNA-Binding Proteins Revisited - The Emerging Arabidopsis mRNA Interactome

RNA-protein interaction is an important checkpoint to tune gene expression at the RNA level. Global identification of proteins binding in vivo to mRNA has been possible through interactome capture - where proteins are fixed to target RNAs by UV crosslinking and purified through affinity capture of polyadenylated RNA. In Arabidopsis over 500 RNA-binding proteins (RBPs) enriched in UV-crosslinked samples have been identified. As in mammals and yeast, the mRNA interactomes came with a few surprises. For example, a plethora of the proteins caught on RNA had not previously been linked to RNA-mediated processes, for example proteins of intermediary metabolism. Thus, the studies provide unprecedented insights into the composition of the mRNA interactome, highlighting the complexity of RNA-mediated processes.

Auxin regulates functional gene groups in a fold-change-specific manner in Arabidopsis thaliana roots

Auxin plays a pivotal role in virtually every aspect of plant morphogenesis. It simultaneously orchestrates a diverse variety of processes such as cell wall biogenesis, transition through the cell cycle, or metabolism of a wide range of chemical substances. The coordination principles for such a complex orchestration are poorly understood at the systems level. Here, we perform an RNA-seq experiment to study the transcriptional response to auxin treatment  within gene groups of different biological processes, molecular functions, or cell components in a quantitative fold-change-specific manner. We find for Arabidopsis thaliana roots treated with auxin for 6 h that (i) there are functional groups within which genes respond to auxin with a surprisingly similar fold changes and that (ii) these fold changes vary from one group to another. These findings make it tempting to conjecture the existence of some transcriptional logic orchestrating the coordinated expression of genes within functional groups in a fold-change-specific manner. To obtain some initial insight about this coordinated expression, we performed a motif enrichment analysis and found cis-regulatory elements TBX1-3, SBX, REG, and TCP/site2 as the candidates conferring fold-change-specific responses to auxin in Arabidopsis thaliana.

Gene expression atlas of pigeonpea and its application to gain insights into genes associated with pollen fertility implicated in seed formation

Pigeonpea (Cajanus cajan) is an important grain legume of the semi-arid tropics, mainly used for its protein rich seeds. To link the genome sequence information with agronomic traits resulting from specific developmental processes, a Cajanus cajan gene expression atlas (CcGEA) was developed using the Asha genotype. Thirty tissues/organs representing developmental stages from germination to senescence were used to generate 590.84 million paired-end RNA-Seq data. The CcGEA revealed a compendium of 28 793 genes with differential, specific, spatio-temporal and constitutive expression during various stages of development in different tissues. As an example to demonstrate the application of the CcGEA, a network of 28 flower-related genes analysed for cis-regulatory elements and splicing variants has been identified. In addition, expression analysis of these candidate genes in male sterile and male fertile genotypes suggested their critical role in normal pollen development leading to seed formation. Gene network analysis also identified two regulatory genes, a pollen-specific SF3 and a sucrose-proton symporter, that could have implications for improvement of agronomic traits such as seed production and yield. In conclusion, the CcGEA provides a valuable resource for pigeonpea to identify candidate genes involved in specific developmental processes and to understand the well-orchestrated growth and developmental process in this resilient crop.

Epitranscriptomic regulation of viral replication

RNA plays central roles in biology and novel functions and regulation mechanisms are constantly emerging. To accomplish some of their functions within the cell, RNA molecules undergo hundreds of chemical modifications from which N6-methyladenosine (m⁶A), inosine (I), pseudouridine (ψ) and 5-methylcytosine (5mC) have been described in eukaryotic mRNA. Interestingly, the m⁶A modification was shown to be reversible, adding novel layers of regulation of gene expression through what is now recognized as epitranscriptomics. The development of molecular mapping strategies coupled to next generation sequencing allowed the identification of thousand of modified transcripts in different tissues and under different physiological conditions such as viral infections. As intracellular parasites, viruses are confronted to cellular RNA modifying enzymes and, as a consequence, viral RNA can be chemically modified at some stages of the replication cycle. This review focuses on the chemical modifications of viral RNA and the impact that these modifications have on viral gene expression and the output of infection. A special emphasis is given to m⁶A, which was recently shown to play important yet controversial roles in different steps of the HIV-1, HCV and ZIKV replication cycles.

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.

N(6)-Methyladenosine RNA Modification Regulates Shoot Stem Cell Fate in Arabidopsis

N(6)-Methyladenosine (m(6)A) represents the most prevalent internal modification on mRNA and requires a multicomponent m(6)A methyltransferase complex in mammals. How their plant counterparts determine the global m(6)A modification landscape and its molecular link to plant development remain unknown. Here we show that FKBP12 INTERACTING PROTEIN 37 KD (FIP37) is a core component of the m(6)A methyltransferase complex, which underlies control of shoot stem cell fate in Arabidopsis. The mutants lacking FIP37 exhibit massive overproliferation of shoot meristems and a transcriptome-wide loss of m(6)A RNA modifications. We further demonstrate that FIP37 mediates m(6)A RNA modification on key shoot meristem genes inversely correlated with their mRNA stability, thus confining their transcript levels to prevent shoot meristem overproliferation. Our results suggest an indispensable role of FIP37 in mediating m(6)A mRNA modification, which is required for maintaining the shoot meristem as a renewable source for continuously producing all aerial organs in plants.