RNA 5-Methylcytosine Analysis by Bisulfite Sequencing

The bisulfite reaction with cytosine and genomic DNA structure

The bisulfite reaction with native DNA has been extensively employed in the detection of non-B DNA structures that can form spontaneously in DNA. These sequences are dynamic in that they can adopt both normal Watson-Crick paired B-DNA or unusual structures like the Triplex, G-Quadruplex, i-motif and Cruciform or Hairpin. Considerable evidence now suggests that these dynamic sequences play roles in both epigenetics and mutagenesis. The bisulfite reaction with native DNA offers a key approach to their detection. In this application whole cells, isolated nuclei or isolated DNA are treated with bisulfite under non-denaturing conditions in order to detect bisulfite accessible regions DNA that are associated with these structures. Here I review the stereochemistry of the bisulfite reaction, the electronic structure of its DNA cytosine substrates and its application in the detection of unusual structures in native DNA.

Dynamic RNA methylation modifications and their regulatory role in mammalian development and diseases

Among over 170 different types of chemical modifications on RNA nucleobases identified so far, RNA methylation is the major type of epitranscriptomic modifications existing on almost all types of RNAs, and has been demonstrated to participate in the entire process of RNA metabolism, including transcription, pre-mRNA alternative splicing and maturation, mRNA nucleus export, mRNA degradation and stabilization, mRNA translation. Attributing to the development of high-throughput detection technologies and the identification of both dynamic regulators and recognition proteins, mechanisms of RNA methylation modification in regulating the normal development of the organism as well as various disease occurrence and developmental abnormalities upon RNA methylation dysregulation have become increasingly clear. Here, we particularly focus on three types of RNA methylations: N6-methylcytosine (m6A), 5-methylcytosine (m5C), and N7-methyladenosine (m7G). We summarize the elements related to their dynamic installment and removal, specific binding proteins, and the development of high-throughput detection technologies. Then, for a comprehensive understanding of their biological significance, we also overview the latest knowledge on the underlying mechanisms and key roles of these three mRNA methylation modifications in gametogenesis, embryonic development, immune system development, as well as disease and tumor progression.

Distinct characteristics of putative archaeal 5-methylcytosine RNA methyltransferases unveil their substrate specificities and evolutionary ancestries

5-Methylcytosine methyltransferases (m5C MTases) are known to be involved in the modification of RNA. Although these enzymes have been relatively well characterized in bacteria and eukarya, a complete understanding of the archaeal counterparts is lacking. In this study, the identification and characterization of archaeal RNA m5C MTases were performed. As a case study, a hyperthermophilic archaeon, Pyrococcus horikoshii OT3, which possesses five putative RNA m5C MTases, was chosen. Among the five putative RNA m5C MTases, two proteins (PH0851 and PH1991) have been characterized as homologs of a bacterial rRNA MTase (RsmB) and eukaryal tRNA MTase (NSUN6), respectively. The in-depth characterization of the remaining three putative RNA m5C MTases (PH1078, PH1374, and PH1537) in this study suggests the presence of the signature architecture and catalytic residues plausibly involved in the binding of their cognate RNA substrates. Additionally, the results also suggest the existence of two RsmB-like proteins (PH0851 and PH1078) belonging to the same subfamily IV of m5C RNA MTase. However, the proteins PH1374 and PH1537 belong to the same subfamily V but bind to different substrates, rRNA and tRNA, respectively. The findings further indicate that archaeal RNA m5C MTases link those from bacteria and eukarya.Communicated by Ramaswamy H. Sarma.

General Principles and Limitations for Detection of RNA Modifications by Sequencing

Among the many analytical methods applied to RNA modifications, a particularly pronounced surge has occurred in the past decade in the field of modification mapping. The occurrence of modifications such as m6A in mRNA, albeit known since the 1980s, became amenable to transcriptome-wide analyses through the advent of next-generation sequencing techniques in a rather sudden manner. The term "mapping" here refers to detection of RNA modifications in a sequence context, which has a dramatic impact on the interpretation of biological functions. As a consequence, an impressive number of mapping techniques were published, most in the perspective of what now has become known as "epitranscriptomics". While more and more different modifications were reported to occur in mRNA, conflicting reports and controversial results pointed to a number of technical and theoretical problems rooted in analytics, statistics, and reagents. Rather than finding the proverbial needle in a haystack, the tasks were to determine how many needles of what color in what size of a haystack one was looking at.As the authors of this Account, we think it important to outline the limitations of different mapping methods since many life scientists freshly entering the field confuse the accuracy and precision of modification mapping with that of normal sequencing, which already features numerous caveats by itself. Indeed, we propose here to qualify a specific mapping method by the size of the transcriptome that can be meaningfully analyzed with it.We here focus on high throughput sequencing by Illumina technology, referred to as RNA-Seq. We noted with interest the development of methods for modification detection by other high throughput sequencing platforms that act directly on RNA, e.g., PacBio SMRT and nanopore sequencing, but those are not considered here.In contrast to approaches relying on direct RNA sequencing, current Illumina RNA-Seq protocols require prior conversion of RNA into DNA. This conversion relies on reverse transcription (RT) to create cDNA; thereafter, the cDNA undergoes a sequencing-by-synthesis type of analysis. Thus, a particular behavior of RNA modified nucleotides during the RT-step is a prerequisite for their detection (and quantification) by deep sequencing, and RT properties have great influence on the detection efficiency and reliability. Moreover, the RT-step requires annealing of a synthetic primer, a prerequisite with a crucial impact on library preparation. Thus, all RNA-Seq protocols must feature steps for the introduction of primers, primer landing sites, or adapters on both the RNA 3'- and 5'-ends.

Regulation and functions of non-m6A mRNA modifications

Nucleobase modifications are prevalent in eukaryotic mRNA and their discovery has resulted in the emergence of epitranscriptomics as a research field. The most abundant internal (non-cap) mRNA modification is N6-methyladenosine (m6A), the study of which has revolutionized our understanding of post-transcriptional gene regulation. In addition, numerous other mRNA modifications are gaining great attention because of their major roles in RNA metabolism, immunity, development and disease. In this Review, we focus on the regulation and function of non-m6A modifications in eukaryotic mRNA, including pseudouridine (Ψ), N6,2'-O-dimethyladenosine (m6Am), N1-methyladenosine (m1A), inosine, 5-methylcytidine (m5C), N4-acetylcytidine (ac4C), 2'-O-methylated nucleotide (Nm) and internal N7-methylguanosine (m7G). We highlight their regulation, distribution, stoichiometry and known roles in mRNA metabolism, such as mRNA stability, translation, splicing and export. We also discuss their biological consequences in physiological and pathological processes. In addition, we cover research techniques to further study the non-m6A mRNA modifications and discuss their potential future applications.

Benchmarking RNA Editing Detection Tools

RNA, like DNA and proteins, can undergo modifications. To date, over 170 RNA modifications have been identified, leading to the emergence of a new research area known as epitranscriptomics. RNA editing is the most frequent RNA modification in mammalian transcriptomes, and two types have been identified: (1) the most frequent, adenosine to inosine (A-to-I); and (2) the less frequent, cysteine to uracil (C-to-U) RNA editing. Unlike other epitranscriptomic marks, RNA editing can be readily detected from RNA sequencing (RNA-seq) data without any chemical conversions of RNA before sequencing library preparation. Furthermore, analyzing RNA editing patterns from transcriptomic data provides an additional layer of information about the epitranscriptome. As the significance of epitranscriptomics, particularly RNA editing, gains recognition in various fields of biology and medicine, there is a growing interest in detecting RNA editing sites (RES) by analyzing RNA-seq data. To cope with this increased interest, several bioinformatic tools are available. However, each tool has its advantages and disadvantages, which makes the choice of the most appropriate tool for bench scientists and clinicians difficult. Here, we have benchmarked bioinformatic tools to detect RES from RNA-seq data. We provide a comprehensive view of each tool and its performance using previously published RNA-seq data to suggest recommendations on the most appropriate for utilization in future studies.

Selective Functionalisation of 5-Methylcytosine by Organic Photoredox Catalysis

The epigenetic modification 5-methylcytosine plays a vital role in development, cell specific gene expression and disease states. The selective chemical modification of the 5-methylcytosine methyl group is challenging. Currently, no such chemistry exists. Direct functionalisation of 5-methylcytosine would improve the detection and study of this epigenetic feature. We report a xanthone-photosensitised process that introduces a 4-pyridine modification at a C(sp3)-H bond in the methyl group of 5-methylcytosine. We propose a reaction mechanism for this type of reaction based on density functional calculations and apply transition state analysis to rationalise differences in observed reaction efficiencies between cyanopyridine derivatives. The reaction is initiated by single electron oxidation of 5-methylcytosine followed by deprotonation to generate the methyl group radical. Cross coupling of the methyl radical with 4-cyanopyridine installs a 4-pyridine label at 5-methylcytosine. We demonstrate use of the pyridination reaction to enrich 5-methylcytosine-containing ribonucleic acid.

Selective Functionalisation of 5-Methylcytosine by Organic Photoredox Catalysis

The epigenetic modification 5-methylcytosine plays a vital role in development, cell specific gene expression and disease states. The selective chemical modification of the 5-methylcytosine methyl group is challenging. Currently, no such chemistry exists. Direct functionalisation of 5-methylcytosine would improve the detection and study of this epigenetic feature. We report a xanthone-photosensitised process that introduces a 4-pyridine modification at a C(sp3 )-H bond in the methyl group of 5-methylcytosine. We propose a reaction mechanism for this type of reaction based on density functional calculations and apply transition state analysis to rationalise differences in observed reaction efficiencies between cyanopyridine derivatives. The reaction is initiated by single electron oxidation of 5-methylcytosine followed by deprotonation to generate the methyl group radical. Cross coupling of the methyl radical with 4-cyanopyridine installs a 4-pyridine label at 5-methylcytosine. We demonstrate use of the pyridination reaction to enrich 5-methylcytosine-containing ribonucleic acid.

Small RNA modifications: regulatory molecules and potential applications

Small RNAs (also referred to as small noncoding RNAs, sncRNA) are defined as polymeric ribonucleic acid molecules that are less than 200 nucleotides in length and serve a variety of essential functions within cells. Small RNA species include microRNA (miRNA), PIWI-interacting RNA (piRNA), small interfering RNA (siRNA), tRNA-derived small RNA (tsRNA), etc. Current evidence suggest that small RNAs can also have diverse modifications to their nucleotide composition that affect their stability as well as their capacity for nuclear export, and these modifications are relevant to their capacity to drive molecular signaling processes relevant to biogenesis, cell proliferation and differentiation. In this review, we highlight the molecular characteristics and cellular functions of small RNA and their modifications, as well as current techniques for their reliable detection. We also discuss how small RNA modifications may be relevant to the clinical applications for the diagnosis and treatment of human health conditions such as cancer.

Navigating the pitfalls of mapping DNA and RNA modifications

Chemical modifications to nucleic acids occur across the kingdoms of life and carry important regulatory information. Reliable high-resolution mapping of these modifications is the foundation of functional and mechanistic studies, and recent methodological advances based on next-generation sequencing and long-read sequencing platforms are critical to achieving this aim. However, mapping technologies may have limitations that sometimes lead to inconsistent results. Some of these limitations are technical in nature and specific to certain types of technology. Here, however, we focus on common (yet not always widely recognized) pitfalls that are shared among frequently used mapping technologies and discuss strategies to help technology developers and users mitigate their effects. Although the emphasis is primarily on DNA modifications, RNA modifications are also discussed.

Quantification of substoichiometric modification reveals global tsRNA hypomodification, preferences for angiogenin-mediated tRNA cleavage, and idiosyncratic epitranscriptomes of human neuronal cell-lines

Modification of tRNA is an integral part of the epitranscriptome with a particularly pronounced potential to generate diversity in RNA expression. Eukaryotic tRNA contains modifications in up to 20% of their nucleotides, but not all sites are always fully modified. Combinations and permutations of partially modified sites in tRNAs can generate a plethora of tRNA isoforms, termed modivariants. Here, we investigate the stoichiometry of incompletely modified sites in tRNAs from human cell lines for their information content. Using a panel of RNA modification mapping methods, we assess the stoichiometry of sites that contain the modifications 5-methylcytidine (m⁵C), 2'-O-ribose methylation (Nm), 3-methylcytidine (m³C), 7-methylguanosine (m⁷G), and Dihydrouridine (D). We discovered that up to 75% of sites can be incompletely modified and that the differential modification status of a cellular tRNA population holds information that allows to discriminate e.g. different cell lines. As a further aspect, we investigated potential causal connectivity between tRNA modification and its processing into tRNA fragments (tiRNAs and tRFs). Upon exposure of cultured living cells to cell-penetrating angiogenin, the modification patterns of the corresponding RNA populations was changed. Importantly, we also found that tsRNAs were significantly less modified than their parent tRNAs at numerous sites, suggesting that tsRNAs might derive chiefly from hypomodified tRNAs.

The Transition from Cancer "omics" to "epi-omics" through Next- and Third-Generation Sequencing

Deciphering cancer etiopathogenesis has proven to be an especially challenging task since the mechanisms that drive tumor development and progression are far from simple. An astonishing amount of research has revealed a wide spectrum of defects, including genomic abnormalities, epigenomic alterations, disturbance of gene transcription, as well as post-translational protein modifications, which cooperatively promote carcinogenesis. These findings suggest that the adoption of a multidimensional approach can provide a much more precise and comprehensive picture of the tumor landscape, hence serving as a powerful tool in cancer research and precision oncology. The introduction of next- and third-generation sequencing technologies paved the way for the decoding of genetic information and the elucidation of cancer-related cellular compounds and mechanisms. In the present review, we discuss the current and emerging applications of both generations of sequencing technologies, also referred to as massive parallel sequencing (MPS), in the fields of cancer genomics, transcriptomics and proteomics, as well as in the progressing realms of epi-omics. Finally, we provide a brief insight into the expanding scope of sequencing applications in personalized cancer medicine and pharmacogenomics.

RNA modifications in aging-associated cardiovascular diseases

Cardiovascular disease (CVD) is a leading cause of morbidity and mortality worldwide that bears an enormous healthcare burden and aging is a major contributing factor to CVDs. Functional gene expression network during aging is regulated by mRNAs transcriptionally and by non-coding RNAs epi-transcriptionally. RNA modifications alter the stability and function of both mRNAs and non-coding RNAs and are involved in differentiation, development, and diseases. Here we review major chemical RNA modifications on mRNAs and non-coding RNAs, including N6-adenosine methylation, N1-adenosine methylation, 5-methylcytidine, pseudouridylation, 2′ -O-ribose-methylation, and N7-methylguanosine, in the aging process with an emphasis on cardiovascular aging. We also summarize the currently available methods to detect RNA modifications and the bioinformatic tools to study RNA modifications. More importantly, we discussed the specific implication of the RNA modifications on mRNAs and non-coding RNAs in the pathogenesis of aging-associated CVDs, including atherosclerosis, hypertension, coronary heart diseases, congestive heart failure, atrial fibrillation, peripheral artery disease, venous insufficiency, and stroke.

Systematic evaluation of parameters in RNA bisulfite sequencing data generation and analysis

The presence of 5-methylcytosine (m5C) in RNA molecules has been known for decades and its importance in regulating RNA metabolism has gradually become appreciated. Despite recent advances made in the functional and mechanistic understanding of RNA m5C modifications, the detection and quantification of methylated RNA remains a challenge. In this study, we compared four library construction procedures for RNA bisulfite sequencing and implemented an analytical pipeline to assess the key parameters in the process of m5C calling. We found that RNA fragmentation after bisulfite conversion increased the yield significantly, and an additional high temperature treatment improved bisulfite conversion efficiency especially for sequence reads mapped to the mitochondrial transcriptome. Using Unique Molecular Identifiers (UMIs), we observed that PCR favors the amplification of unmethylated templates. The low sequencing quality of bisulfite-converted bases is a major contributor to the methylation artifacts. In addition, we found that mitochondrial transcripts are frequently resistant to bisulfite conversion and no p-m5C sites with high confidence could be identified on mitochondrial mRNAs. Taken together, this study reveals the various sources of artifacts in RNA bisulfite sequencing data and provides an improved experimental procedure together with analytical methodology.

Cross-Talk Between m6A- and m5C-Related lncRNAs to Construct a Novel Signature and Predict the Immune Landscape of Colorectal Cancer Patients

Background

N6-methyladenosine (m⁶A) and 5-methylcytosine (m⁵C) can modify long non-coding RNAs (lncRNAs), thereby affecting tumorigenesis and tumor progression. However, there is a lack of knowledge regarding the potential roles and cross-talk of m⁶A- and m⁵C-related lncRNAs in the tumor microenvironment (TME) and their effect on prognosis.

Methods

We systematically evaluated the expression patterns of m⁶A- and m⁵C-related lncRNAs in 1358 colorectal cancer (CRC) samples from four datasets. Consensus clustering was conducted to identify molecular subtypes of CRC, and the clinical significance, TME, tumor-infiltrating immune cells (TIICs), and immune checkpoints in the different molecular subtypes were analyzed. Finally, we established a m⁶A- and m⁵C-related lncRNA signature and a prognostic nomogram.

Results

We identified 141 m⁶A- and m⁵C-related lncRNAs by co-expression analysis, among which 23 lncRNAs were significantly associated with the overall survival (OS) of CRC patients. Two distinct molecular subtypes (cluster A and cluster B) were identified, and these two distinct molecular subtypes could predict clinicopathological features, prognosis, TME stromal activity, TIICs, immune checkpoints. Next, a m⁶A- and m⁵C-related lncRNA signature for predicting OS was constructed, and its predictive capability in CRC patients was validated. We then constructed a highly accurate nomogram for improving the clinical applicability of the signature. Analyses of clinicopathological features, prognosis, TIICs, cancer stem cell (CSC), and drug response revealed significant differences between two risk groups. In addition, we found that patients with a low-risk score exhibited enhanced response to anti-PD-1/L1 immunotherapy. Functional enrichment analysis showed that these lncRNAs related to the high-risk group were involved in the development and progression of CRC.

Conclusions

We conducted a comprehensive analysis of m⁶A- and m⁵C-related lncRNAs in CRC and revealed their potential functions in predicting tumor-immune-stromal microenvironment, clinicopathological features, and prognosis, and determined their role in immunotherapy. These findings may improve our understanding of the cross-talk between m⁶A- and m⁵C-related lncRNAs in CRC and pave a new road for prognosis assessment and more effective immunotherapy strategies.

Analysis of RNA Modifications by Second- and Third-Generation Deep Sequencing: 2020 Update

The precise mapping and quantification of the numerous RNA modifications that are present in tRNAs, rRNAs, ncRNAs/miRNAs, and mRNAs remain a major challenge and a top priority of the epitranscriptomics field. After the keystone discoveries of massive m⁶A methylation in mRNAs, dozens of deep sequencing-based methods and protocols were proposed for the analysis of various RNA modifications, allowing us to considerably extend the list of detectable modified residues. Many of the currently used methods rely on the particular reverse transcription signatures left by RNA modifications in cDNA; these signatures may be naturally present or induced by an appropriate enzymatic or chemical treatment. The newest approaches also include labeling at RNA abasic sites that result from the selective removal of RNA modification or the enhanced cleavage of the RNA ribose-phosphate chain (perhaps also protection from cleavage), followed by specific adapter ligation. Classical affinity/immunoprecipitation-based protocols use either antibodies against modified RNA bases or proteins/enzymes, recognizing RNA modifications. In this survey, we review the most recent achievements in this highly dynamic field, including promising attempts to map RNA modifications by the direct single-molecule sequencing of RNA by nanopores.

Deciphering RNA modifications at base resolution: from chemistry to biology

Nearly 200 distinct chemical modifications of RNAs have been discovered to date. Their analysis via direct methods has been possible in abundant RNA species, such as ribosomal, transfer or viral RNA, since several decades. However, their analysis in less abundant RNAs species, especially cellular messenger RNAs, was rendered possible only recently with the advent of high throughput sequencing techniques. Given the growing biomedical interest of the proteins that write, erase and read RNA modifications, ingenious new methods to enrich and identify RNA modifications at base resolution have been implemented, and more efforts are underway to render them more quantitative. Here, we review several crucial modification-specific (bio)chemical approaches and discuss their advantages and shortcomings for exploring the epitranscriptome.

Dynamic transcriptomic m5 C and its regulatory role in RNA processing

RNA 5-methylcytosine (m⁵ C) is a prevalent RNA modification in multiple RNA species, including messenger RNAs (mRNAs), transfer RNAs (tRNAs), ribosomal RNAs (rRNAs), and noncoding RNAs (ncRNAs), and broadly distributed from archaea, prokaryotes to eukaryotes. The multiple detecting techniques of m⁵ C have been developed, such as m⁵ C-RIP-seq, miCLIP-seq, AZA-IP-seq, RNA-BisSeq, TAWO-seq, and Nanopore sequencing. These high-throughput techniques, combined with corresponding analysis pipeline, provide a precise m⁵ C landscape contributing to the deciphering of its biological functions. The m⁵ C modification is distributed along with mRNA and enriched around 5'UTR and 3'UTR, and conserved in tRNAs and rRNAs. It is dynamically regulated by its related enzymes, including methyltransferases (NSUN, DNMT, and TRDMT family members), demethylases (TET families and ALKBH1), and binding proteins (ALYREF and YBX1). So far, accumulative studies have revealed that m⁵ C participates in a variety of RNA metabolism, including mRNA export, RNA stability, and translation. Depletion of m⁵ C modification in the organism could cause dysfunction of mitochondria, drawback of stress response, frustration of gametogenesis and embryogenesis, abnormality of neuro and brain development, and has been implicated in cell migration and tumorigenesis. In this review, we provide a comprehensive summary of dynamic regulatory elements of RNA m⁵ C, including methyltransferases (writers), demethylases (erasers), and binding proteins (readers). We also summarized the related detecting technologies and biological functions of the RNA 5-methylcytosine, and provided future perspectives in m⁵ C research. This article is categorized under: RNA Processing > RNA Editing and Modification.

Manganese Ions Individually Alter the Reverse Transcription Signature of Modified Ribonucleosides

Reverse transcription of RNA templates containing modified ribonucleosides transfers modification-related information as misincorporations, arrest or nucleotide skipping events to the newly synthesized cDNA strand. The frequency and proportion of these events, merged from all sequenced cDNAs, yield a so-called RT signature, characteristic for the respective RNA modification and reverse transcriptase (RT). While known for DNA polymerases in so-called error-prone PCR, testing of four different RTs by replacing Mg²⁺ with Mn²⁺ in reaction buffer revealed the immense influence of manganese chloride on derived RT signatures, with arrest rates on m¹A positions dropping from 82% down to 24%. Additionally, we observed a vast increase in nucleotide skipping events, with single positions rising from 4% to 49%, thus implying an enhanced read-through capability as an effect of Mn²⁺ on the reverse transcriptase, by promoting nucleotide skipping over synthesis abortion. While modifications such as m¹A, m²₂G, m¹G and m³C showed a clear influence of manganese ions on their RT signature, this effect was individual to the polymerase used. In summary, the results imply a supporting effect of Mn²⁺ on reverse transcription, thus overcoming blockades in the Watson-Crick face of modified ribonucleosides and improving both read-through rate and signal intensity in RT signature analysis.

The critical function of the plastid rRNA methyltransferase, CMAL, in ribosome biogenesis and plant development

Methylation of nucleotides in ribosomal RNAs (rRNAs) is a ubiquitous feature that occurs in all living organisms. The formation of methylated nucleotides is performed by a variety of RNA-methyltransferases. Chloroplasts of plant cells result from an endosymbiotic event and possess their own genome and ribosomes. However, enzymes responsible for rRNA methylation and the function of modified nucleotides in chloroplasts remain to be determined. Here, we identified an rRNA methyltransferase, CMAL (Chloroplast MraW-Like), in the Arabidopsis chloroplast and investigated its function. CMAL is the Arabidopsis ortholog of bacterial MraW/ RsmH proteins and accounts to the N4-methylation of C1352 in chloroplast 16S rRNA, indicating that CMAL orthologs and this methyl-modification nucleotide is conserved between bacteria and the endosymbiont-derived eukaryotic organelle. The knockout of CMAL in Arabidopsis impairs the chloroplast ribosome accumulation and accordingly reduced the efficiency of mRNA translation. Interestingly, the loss of CMAL leads not only to defects in chloroplast function, but also to abnormal leaf and root development and overall plant morphology. Further investigation showed that CMAL is involved in the plant development probably by modulating auxin derived signaling pathways. This study uncovered the important role of 16S rRNA methylation mediated by CMAL in chloroplast ribosome biogenesis and plant development.

Epitranscriptomic 5-Methylcytosine Profile in PM2.5-induced Mouse Pulmonary Fibrosis

Exposure of airborne particulate matter (PM) with an aerodynamic diameter less than 2.5 μm (PM2.5) is epidemiologically associated with lung dysfunction and respiratory symptoms, including pulmonary fibrosis. However, whether epigenetic mechanisms are involved in PM2.5-induced pulmonary fibrosis is currently poorly understood. Herein, using a PM2.5-induced pulmonary fibrosis mouse model, we found that PM2.5 exposure leads to aberrant mRNA 5-methylcytosine (m5C) gain and loss in fibrotic lung tissues. Moreover, we showed the m5C-mediated regulatory map of gene functions in pulmonary fibrosis after PM2.5 exposure. Several genes act as m5C gain-upregulated factors, probably critical for the development of PM2.5-induced fibrosis in mouse lungs. These genes, including Lcn2, Mmp9, Chi3l1, Adipoq, Atp5j2, Atp5l, Atpif1, Ndufb6, Fgr, Slc11a1, and Tyrobp, are highly related to oxidative stress response, inflammatory responses, and immune system processes. Our study illustrates the first epitranscriptomic RNA m5C profile in PM2.5-induced pulmonary fibrosis and will be valuable in identifying biomarkers for PM2.5 exposure-related lung pathogenesis with translational potential.

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.

Dynamic methylome of internal mRNA N7-methylguanosine and its regulatory role in translation

Over 150 types of RNA modifications are identified in RNA molecules. Transcriptome profiling is one of the key steps in decoding the epitranscriptomic panorama of these chemical modifications and their potential functions. N⁷-methylguanosine (m⁷G) is one of the most abundant modifications present in tRNA, rRNA and mRNA 5'cap, and has critical roles in regulating RNA processing, metabolism and function. Besides its presence at the cap position in mRNAs, m⁷G is also identified in internal mRNA regions. However, its transcriptome-wide distribution and dynamic regulation within internal mRNA regions remain unknown. Here, we have established m⁷G individual-nucleotide-resolution cross-linking and immunoprecipitation with sequencing (m⁷G miCLIP-seq) to specifically detect internal mRNA m⁷G modification. Using this approach, we revealed that m⁷G is enriched at the 5'UTR region and AG-rich contexts, a feature that is well-conserved across different human/mouse cell lines and mouse tissues. Strikingly, the internal m⁷G modification is dynamically regulated under both H₂O₂ and heat shock treatments, with remarkable accumulations in the CDS and 3'UTR regions, and functions in promoting mRNA translation efficiency. Consistently, a PCNA 3'UTR minigene reporter harboring the native m⁷G modification site displays both enriched m⁷G modification and increased mRNA translation upon H₂O₂ treatment compared to the m⁷G site-mutated minigene reporter (G to A). Taken together, our findings unravel the dynamic profiles of internal mRNA m⁷G methylome and highlight m⁷G as a novel epitranscriptomic marker with regulatory roles in translation.

METTL1 Promotes let-7 MicroRNA Processing via m7G Methylation

7-methylguanosine (m7G) is present at mRNA caps and at defined internal positions within tRNAs and rRNAs. However, its detection within low-abundance mRNAs and microRNAs (miRNAs) has been hampered by a lack of sensitive detection strategies. Here, we adapt a chemical reactivity assay to detect internal m7G in miRNAs. Using this technique (Borohydride Reduction sequencing [BoRed-seq]) alongside RNA immunoprecipitation, we identify m7G within a subset of miRNAs that inhibit cell migration. We show that the METTL1 methyltransferase mediates m7G methylation within miRNAs and that this enzyme regulates cell migration via its catalytic activity. Using refined mass spectrometry methods, we map m7G to a single guanosine within the let-7e-5p miRNA. We show that METTL1-mediated methylation augments let-7 miRNA processing by disrupting an inhibitory secondary structure within the primary miRNA transcript (pri-miRNA). These results identify METTL1-dependent N7-methylation of guanosine as a new RNA modification pathway that regulates miRNA structure, biogenesis, and cell migration.

5-Methylcytosine Analysis by RNA-BisSeq

5-Methylcytosine (m⁵C) is a posttranscriptional RNA modification identified in both stable and highly abundant tRNAs and rRNAs, and in mRNAs. Many known or novel m⁵C sites have been validated by using advanced high-throughput techniques combined with next-generation sequencing (NGS), especially RNA bisulfite sequencing (RNA-BisSeq). Here we introduce an optimized RNA-BisSeq method by using ACT random hexamers to prime the reverse transcription of bisulfite-treated RNA samples to detect the m⁵C sites.

Transcriptome-wide profiling of multiple RNA modifications simultaneously at single-base resolution

The field of RNA modification would be significantly advanced by the development of sensitive, accurate, single-base resolution methods for profiling multiple common RNA modifications in the same RNA molecule. Our work provides several advances toward that goal, including (i) quantitative methods for profiling Ψ sites at true base-pair resolution transcriptome-wide, (ii) a chemical understanding of our observed Ψ-dependent deletion signature, (iii) improved methods for profiling m⁵C and m¹A, and (iv) a coupling of these methods for the simultaneous detection of all three modifications in the same RNA. Together, the combinatorial ability and relative ease of execution provided by this procedure should greatly forward epitranscriptome studies involving these three very common RNA modifications.

The breadth and importance of RNA modifications are growing rapidly as modified ribonucleotides can impact the sequence, structure, function, stability, and fate of RNAs and their interactions with other molecules. Therefore, knowing cellular RNA modifications at single-base resolution could provide important information regarding cell status and fate. A current major limitation is the lack of methods that allow the reproducible profiling of multiple modifications simultaneously, transcriptome-wide and at single-base resolution. Here we developed RBS-Seq, a modification of RNA bisulfite sequencing that enables the sensitive and simultaneous detection of m⁵C, Ψ, and m¹A at single-base resolution transcriptome-wide. With RBS-Seq, m⁵C and m¹A are accurately detected based on known signature base mismatches and are detected here simultaneously along with Ψ sites that show a 1–2 base deletion. Structural analyses revealed the mechanism underlying the deletion signature, which involves Ψ-monobisulfite adduction, heat-induced ribose ring opening, and Mg²⁺-assisted reorientation, causing base-skipping during cDNA synthesis. Detection of each of these modifications through a unique chemistry allows high-precision mapping of all three modifications within the same RNA molecule, enabling covariation studies. Application of RBS-Seq on HeLa RNA revealed almost all known m⁵C, m¹A, and ψ sites in tRNAs and rRNAs and provided hundreds of new m⁵C and Ψ sites in noncoding RNAs and mRNAs. However, our results diverge greatly from earlier work, suggesting ∼10-fold fewer m⁵C sites in noncoding and coding RNAs and the absence of substantial m¹A in mRNAs. Taken together, the approaches and refined datasets in this work will greatly enable future epitranscriptome studies.

Nucleotide resolution sequencing of N4-acetylcytidine in RNA

Posttranscriptional modifications of RNA represent an emerging class of regulatory elements in human biology. Improved methods for studying how these elements are controlled and where they occur has the potential to transform our understanding of gene expression in development and disease. Here we describe a chemical method for nucleotide resolution sequencing of N4-acetylcytidine (ac4C), a highly conserved modified nucleobase whose formation is catalyzed by the essential cytidine acetyltransferase enzyme NAT10. This approach enables the sensitive, PCR-amplifiable detection of individual ac4C sites from nanograms of unfractionated cellular RNA. The sensitive and quantitative nature of this assay provides a powerful tool to understand how cytidine acetylation is targeted, profile RNA acetyltransferase dynamics, and validate the sites and stoichiometry of ac4C in novel RNA species.

Epitranscriptomics: RNA Modifications in Bacteria and Archaea

The increasingly complex functionality of RNA is contrasted by its simple chemical composition. RNA is generally built from only four different nucleotides (adenine, guanine, cytosine, and uracil). To date, >160 chemical modifications are known to decorate RNA molecules and thereby alter their function or stability. Many RNA modifications are conserved throughout bacteria, archaea, and eukaryotes, while some are unique to each branch of life. Most known modifications occur at internal positions, while there is limited diversity at the termini. The dynamic nature of RNA modifications and newly discovered regulatory functions of some of these RNA modifications gave birth to a new field, now often referred to as "epitranscriptomics." This review highlights the major developments in this field and summarizes detection principles for internal as well as 5'-terminal mRNA modifications in prokaryotes and archaea to investigate their biological significance.

Clustering Count-based RNA Methylation Data Using a Nonparametric Generative Model

Background:

RNA methylome has been discovered as an important layer of gene regulation and can be profiled directly with count-based measurements from high-throughput sequencing data. Although the detailed regulatory circuit of the epitranscriptome remains uncharted, clustering effect in methylation status among different RNA methylation sites can be identified from transcriptome-wide RNA methylation profiles and may reflect the epitranscriptomic regulation. Count-based RNA methylation sequencing data has unique features, such as low reads coverage, which calls for novel clustering approaches. <P><P> Objective: Besides the low reads coverage, it is also necessary to keep the integer property to approach clustering analysis of count-based RNA methylation sequencing data. <P><P> Method: We proposed a nonparametric generative model together with its Gibbs sampling solution for clustering analysis. The proposed approach implements a beta-binomial mixture model to capture the clustering effect in methylation level with the original count-based measurements rather than an estimated continuous methylation level. Besides, it adopts a nonparametric Dirichlet process to automatically determine an optimal number of clusters so as to avoid the common model selection problem in clustering analysis. <P><P> Results: When tested on the simulated system, the method demonstrated improved clustering performance over hierarchical clustering, K-means, MClust, NMF and EMclust. It also revealed on real dataset two novel RNA N6-methyladenosine (m6A) co-methylation patterns that may be induced directly by METTL14 and WTAP, which are two known regulatory components of the RNA m6A methyltransferase complex. <P><P> Conclusion: Our proposed DPBBM method not only properly handles the count-based measurements of RNA methylation data from sites of very low reads coverage, but also learns an optimal number of clusters adaptively from the data analyzed. <P><P> Availability: The source code and documents of DPBBM R package are freely available through the Comprehensive R Archive Network (CRAN): https://cran.r-project.org/web/packages/DPBBM/.

The dynamic RNA modification 5-methylcytosine and its emerging role as an epitranscriptomic mark

It is a well‐known fact that RNA is the target of a plethora of modifications which currently amount to over a hundred. The vast majority of these modifications was observed in the two most abundant classes of RNA, rRNA and tRNA. With the recent advance in mapping technologies, modifications have been discovered also in mRNA and in less abundant non‐coding RNA species. These developments have sparked renewed interest in elucidating the nature and functions of those “epitransciptomic” modifications in RNA. N6‐methyladenosine (m⁶A) is the best understood and most frequent mark of mRNA with demonstrated functions ranging from pre‐mRNA processing, translation, miRNA biogenesis to mRNA decay. By contrast, much less research has been conducted on 5‐methylcytosine (m5C), which was detected in tRNAs and rRNAs and more recently in poly(A)RNAs. In this review, we discuss recent developments in the discovery of m5C RNA methylomes, the functions of m5C as well as the proteins installing, translating and manipulating this modification. Although our knowledge about m5C in RNA transcripts is just beginning to consolidate, it has become clear that cytosine methylation represents a powerful mechanistic strategy to regulate cellular processes on an epitranscriptomic level.

This article is categorized under:

RNA Processing > RNA Editing and Modification

RNA Interactions with Proteins and Other Molecules > Protein–RNA Interactions: Functional Implications

RNA Processing > tRNA Processing

RNA Turnover and Surveillance > Regulation of RNA Stability

Topological Characterization of Human and Mouse m5C Epitranscriptome Revealed by Bisulfite Sequencing

Background

Compared with the well-studied 5-methylcytosine (m⁵C) in DNA, the role and topology of epitranscriptome m⁵C remain insufficiently characterized.

Results

Through analyzing transcriptome-wide m⁵C distribution in human and mouse, we show that the m⁵C modification is significantly enriched at 5′ untranslated regions (5′UTRs) of mRNA in human and mouse. With a comparative analysis of the mRNA and DNA methylome, we demonstrate that, like DNA methylation, transcriptome m⁵C methylation exhibits a strong clustering effect. Surprisingly, an inverse correlation between mRNA and DNA m⁵C methylation is observed at CpG sites. Further analysis reveals that RNA m⁵C methylation level is positively correlated with both RNA expression and RNA half-life. We also observed that the methylation level of mitochondrial RNAs is significantly higher than RNAs transcribed from the nuclear genome.

Conclusions

This study provides an in-depth topological characterization of transcriptome-wide m⁵C modification by associating RNA m⁵C methylation patterns with transcriptional expression, DNA methylations, RNA stabilities, and mitochondrial genome.

Emerging role of dynamic RNA modifications during animal development

The central dogma of molecular biology statically says that the information flows from DNA to messenger RNA to protein. But the recent advances in mass spectrometry and high throughput technology have helped the scientists to view RNA as little more than a courier of genetic information encoded in the DNA. The dynamics of RNA modifications in coding and non-coding RNAs are just emerging as a carrier of non-genetic information, uncovering a new layer of complexity in the regulation of gene expression and protein translation. In this review, we summarize about the current knowledge of N6-methyladenosine (m6A), N1-methyladenosine (m1A), 5-methylcytosine (m5C) and pseudouridine (Ψ) modifications in RNA, and described how these RNA modifications are implicated in early animal development and in several human diseases.

Understanding RNA modifications: the promises and technological bottlenecks of the 'epitranscriptome'

The discovery of mechanisms that alter genetic information via RNA editing or introducing covalent RNA modifications points towards a complexity in gene expression that challenges long-standing concepts. Understanding the biology of RNA modifications represents one of the next frontiers in molecular biology. To this date, over 130 different RNA modifications have been identified, and improved mass spectrometry approaches are still adding to this list. However, only recently has it been possible to map selected RNA modifications at single-nucleotide resolution, which has created a number of exciting hypotheses about the biological function of RNA modifications, culminating in the proposition of the ‘epitranscriptome’. Here, we review some of the technological advances in this rapidly developing field, identify the conceptual challenges and discuss approaches that are needed to rigorously test the biological function of specific RNA modifications.

Statistical Methods for Transcriptome-Wide Analysis of RNA Methylation by Bisulfite Sequencing

For the transcriptome-wide detection and quantification of the 5-methylcytosine (m⁵C) methylation modification of RNA, one experimental approach is via bisulfite conversion. In this chapter we discuss statistical methods, and a corresponding computational pipeline, to perform transcriptome-wide differential m⁵C methylation analysis between RNA samples, specialized for this assay.

5-methylcytosine promotes mRNA export - NSUN2 as the methyltransferase and ALYREF as an m5C reader

5-methylcytosine (m⁵C) is a post-transcriptional RNA modification identified in both stable and highly abundant tRNAs and rRNAs, and in mRNAs. However, its regulatory role in mRNA metabolism is still largely unknown. Here, we reveal that m⁵C modification is enriched in CG-rich regions and in regions immediately downstream of translation initiation sites and has conserved, tissue-specific and dynamic features across mammalian transcriptomes. Moreover, m⁵C formation in mRNAs is mainly catalyzed by the RNA methyltransferase NSUN2, and m⁵C is specifically recognized by the mRNA export adaptor ALYREF as shown by in vitro and in vivo studies. NSUN2 modulates ALYREF's nuclear-cytoplasmic shuttling, RNA-binding affinity and associated mRNA export. Dysregulation of ALYREF-mediated mRNA export upon NSUN2 depletion could be restored by reconstitution of wild-type but not methyltransferase-defective NSUN2. Our study provides comprehensive m⁵C profiles of mammalian transcriptomes and suggests an essential role for m⁵C modification in mRNA export and post-transcriptional regulation.

High-throughput identification of C/D box snoRNA targets with CLIP and RiboMeth-seq

High-throughput sequencing has greatly facilitated the discovery of long and short non-coding RNAs (ncRNAs), which frequently guide ribonucleoprotein complexes to RNA targets, to modulate their metabolism and expression. However, for many ncRNAs, the targets remain to be discovered. In this study, we developed computational methods to map C/D box snoRNA target sites using data from core small nucleolar ribonucleoprotein crosslinking and immunoprecipitation and from transcriptome-wide mapping of 2΄-O-ribose methylation sites. We thereby assigned the snoRNA guide to a known methylation site in the 18S rRNA, we uncovered a novel partially methylated site in the 28S ribosomal RNA, and we captured a site in the 28S rRNA in interaction with multiple snoRNAs. Although we also captured mRNAs in interaction with snoRNAs, we did not detect 2΄-O-methylation of these targets. Our study provides an integrated approach to the comprehensive characterization of 2΄-O-methylation targets of snoRNAs in species beyond those in which these interactions have been traditionally studied and contributes to the rapidly developing field of 'epitranscriptomics'.

Analysis of High-Throughput RNA Bisulfite Sequencing Data

Methylation of the 5-cytosine (m⁵C) is a common but not well-understood RNA modification, which can be detected by sequencing of bisulfite-treated transcripts (RNA-BSseq). In this Chapter, we discuss computational RNA-BSseq data analysis methods for transcriptome-wide identification and quantification of m⁵C.

Rapid and dynamic transcriptome regulation by RNA editing and RNA modifications

Advances in next-generation sequencing and mass spectrometry have revealed widespread messenger RNA modifications and RNA editing, with dramatic effects on mammalian transcriptomes. Factors introducing, deleting, or interpreting specific modifications have been identified, and analogous with epigenetic terminology, have been designated "writers," "erasers," and "readers." Such modifications in the transcriptome are referred to as epitranscriptomic changes and represent a fascinating new layer of gene expression regulation that has only recently been appreciated. Here, we outline how RNA editing and RNA modification can rapidly affect gene expression, making both processes as well suited to respond to cellular stress and to regulate the transcriptome during development or circadian periods.