Dynamics of epitranscriptomes uncover translational reprogramming directed by ac4C in rice during pathogen infection

Messenger RNA modifications play pivotal roles in RNA biology, but comprehensive landscape changes of epitranscriptomes remain largely unknown in plant immune response. Here we report translational reprogramming directed by ac4C mRNA modification upon pathogen challenge. We first investigate the dynamics of translatomes and epitranscriptomes and uncover that the change in ac4C at single-base resolution promotes translational reprogramming upon Magnaporthe oryzae infection. Then by characterizing the specific distributions of m1A, 2'O-Nm, ac4C, m5C, m6A and m7G, we find that ac4Cs, unlike other modifications, are enriched at the 3rd position of codons, which stabilizes the Watson-Crick base pairing. Importantly, we demonstrate that upon pathogen infection, the increased expression of the ac4C writer OsNAT10/OsACYR (N-ACETYLTRANSFERASE FOR CYTIDINE IN RNA) promotes translation to facilitate rapid activation of immune responses, including the enhancement of jasmonic acid biosynthesis. Our study provides an atlas of mRNA modifications and insights into ac4C function in plant immunity.

Splice_sim: a nucleotide conversion-enabled RNA-seq simulation and evaluation framework

Nucleotide conversion RNA sequencing techniques interrogate chemical RNA modifications in cellular transcripts, resulting in mismatch-containing reads. Biases in mapping the resulting reads to reference genomes remain poorly understood. We present splice_sim, a splice-aware RNA-seq simulation and evaluation pipeline that introduces user-defined nucleotide conversions at set frequencies, creates mixture models of converted and unconverted reads, and calculates mapping accuracies per genomic annotation. By simulating nucleotide conversion RNA-seq datasets under realistic experimental conditions, including metabolic RNA labeling and RNA bisulfite sequencing, we measure mapping accuracies of state-of-the-art spliced-read mappers for mouse and human transcripts and derive strategies to prevent biases in the data interpretation.

Prediction of m6A and m5C at single-molecule resolution reveals a transcriptome-wide co-occurrence of RNA modifications

The epitranscriptome embodies many new and largely unexplored functions of RNA. A significant roadblock hindering progress in epitranscriptomics is the identification of more than one modification in individual transcript molecules. We address this with CHEUI (CH3 (methylation) Estimation Using Ionic current). CHEUI predicts N6-methyladenosine (m6A) and 5-methylcytosine (m5C) in individual molecules from the same sample, the stoichiometry at transcript reference sites, and differential methylation between any two conditions. CHEUI processes observed and expected nanopore direct RNA sequencing signals to achieve high single-molecule, transcript-site, and stoichiometry accuracies in multiple tests using synthetic RNA standards and cell line data. CHEUI's capability to identify two modification types in the same sample reveals a co-occurrence of m6A and m5C in individual mRNAs in cell line and tissue transcriptomes. CHEUI provides new avenues to discover and study the function of the epitranscriptome.

RNA m5C modification upregulates E2F1 expression in a manner dependent on YBX1 phase separation and promotes tumor progression in ovarian cancer

5-Methylcytosine (m5C) is a common RNA modification that modulates gene expression at the posttranscriptional level, but the crosstalk between m5C RNA modification and biomolecule condensation, as well as transcription factor-mediated transcriptional regulation, in ovarian cancer, is poorly understood. In this study, we revealed that the RNA methyltransferase NSUN2 facilitates mRNA m5C modification and forms a positive feedback regulatory loop with the transcription factor E2F1 in ovarian cancer. Specifically, NSUN2 promotes m5C modification of E2F1 mRNA and increases its stability, and E2F1 binds to the NSUN2 promoter, subsequently reciprocally activating NSUN2 transcription. The RNA binding protein YBX1 functions as the m5C reader and is involved in NSUN2-mediated E2F1 regulation. m5C modification promotes YBX1 phase separation, which upregulates E2F1 expression. In ovarian cancer, NSUN2 and YBX1 are amplified and upregulated, and higher expression of NSUN2 and YBX1 predicts a worse prognosis for ovarian cancer patients. Moreover, E2F1 transcriptionally regulates the expression of the oncogenes MYBL2 and RAD54L, driving ovarian cancer progression. Thus, our study delineates a NSUN2-E2F1-NSUN2 loop regulated by m5C modification in a manner dependent on YBX1 phase separation, and this previously unidentified pathway could be a promising target for ovarian cancer treatment.

SRSF2 plays an unexpected role as reader of m5C on mRNA, linking epitranscriptomics to cancer

A common mRNA modification is 5-methylcytosine (m5C), whose role in gene-transcript processing and cancer remains unclear. Here, we identify serine/arginine-rich splicing factor 2 (SRSF2) as a reader of m5C and impaired SRSF2 m5C binding as a potential contributor to leukemogenesis. Structurally, we identify residues involved in m5C recognition and the impact of the prevalent leukemia-associated mutation SRSF2P95H. We show that SRSF2 binding and m5C colocalize within transcripts. Furthermore, knocking down the m5C writer NSUN2 decreases mRNA m5C, reduces SRSF2 binding, and alters RNA splicing. We also show that the SRSF2P95H mutation impairs the ability of the protein to read m5C-marked mRNA, notably reducing its binding to key leukemia-related transcripts in leukemic cells. In leukemia patients, low NSUN2 expression leads to mRNA m5C hypomethylation and, combined with SRSF2P95H, predicts poor outcomes. Altogether, we highlight an unrecognized mechanistic link between epitranscriptomics and a key oncogenesis driver.

TET2-mediated mRNA demethylation regulates leukemia stem cell homing and self-renewal

TET2 is recurrently mutated in acute myeloid leukemia (AML) and its deficiency promotes leukemogenesis (driven by aggressive oncogenic mutations) and enhances leukemia stem cell (LSC) self-renewal. However, the underlying cellular/molecular mechanisms have yet to be fully understood. Here, we show that Tet2 deficiency significantly facilitates leukemogenesis in various AML models (mediated by aggressive or less aggressive mutations) through promoting homing of LSCs into bone marrow (BM) niche to increase their self-renewal/proliferation. TET2 deficiency in AML blast cells increases expression of Tetraspanin 13 (TSPAN13) and thereby activates the CXCR4/CXCL12 signaling, leading to increased homing/migration of LSCs into BM niche. Mechanistically, TET2 deficiency results in the accumulation of methyl-5-cytosine (m5C) modification in TSPAN13 mRNA; YBX1 specifically recognizes the m5C modification and increases the stability and expression of TSPAN13 transcripts. Collectively, our studies reveal the functional importance of TET2 in leukemogenesis, leukemic blast cell migration/homing, and LSC self-renewal as an mRNA m5C demethylase.

Dynamic DNA 5-hydroxylmethylcytosine and RNA 5-methycytosine Reprogramming During Early Human Development

After implantation, complex and highly specialized molecular events render functionally distinct organ formation, whereas how the epigenome shapes organ-specific development remains to be fully elucidated. Here, nano-hmC-Seal, RNA bisulfite sequencing (RNA-BisSeq), and RNA sequencing (RNA-Seq) were performed, and the first multilayer landscapes of DNA 5-hydroxymethylcytosine (5hmC) and RNA 5-methylcytosine (m5C) epigenomes were obtained in the heart, kidney, liver, and lung of the human foetuses at 13-28 weeks with 123 samples in total. We identified 70,091 and 503 organ- and stage-specific differentially hydroxymethylated regions (DhMRs) and m5C-modified mRNAs, respectively. The key transcription factors (TFs), T-box transcription factor 20 (TBX20), paired box 8 (PAX8), krueppel-like factor 1 (KLF1), transcription factor 21 (TCF21), and CCAAT enhancer binding protein beta (CEBPB), specifically contribute to the formation of distinct organs at different stages. Additionally, 5hmC-enriched Alu elements may participate in the regulation of expression of TF-targeted genes. Our integrated studies reveal a putative essential link between DNA modification and RNA methylation, and illustrate the epigenetic maps during human foetal organogenesis, which provide a foundation for for an in-depth understanding of the epigenetic mechanisms underlying early development and birth defects.

NSUN2-mediated mRNA m5C Modification Regulates the Progression of Hepatocellular Carcinoma

RNA modifications affect many biological processes and physiological diseases. The 5-methylcytosine (m5C) modification regulates the progression of multiple tumors. However, its characteristics and functions in hepatocellular carcinoma (HCC) remain largely unknown. Here, we found that HCC tissues had a higher m5C methylation level than the adjacent normal tissues. Transcriptome analysis revealed that the hypermethylated genes mainly participated in the phosphokinase signaling pathways, such as the Ras and PI3K-Akt pathways. The m5C methyltransferase NSUN2 was highly expressed in HCC tissues. Interestingly, the expression of many genes was positively correlated with the expression of NSUN2, including GRB2, RNF115, AATF, ADAM15, RTN3, and HDGF. Real-time PCR assays further revealed that the expression of the mRNAs of GRB2, RNF115, and AATF decreased significantly with the down-regulation of NSUN2 expression in HCC cells. Furthermore, NSUN2 could regulate the cellular sensitivity of HCC cells to sorafenib via modulating the Ras signaling pathway. Moreover, knocking down NSUN2 caused cell cycle arrest. Taken together, our study demonstrates the vital role of NSUN2 in the progression of HCC.

No evidence for epitranscriptomic m5C modification of SARS-CoV-2, HIV and MLV viral RNA

The addition of chemical groups to cellular RNA to modulate RNA fate and/or function is summarized under the term epitranscriptomic modification. More than 170 different modifications have been identified on cellular RNA, such as tRNA, rRNA and, to a lesser extent, on other RNA types. Recently, epitranscriptomic modification of viral RNA has received considerable attention as a possible additional mechanism regulating virus infection and replication. N6-methyladenosine (m6A) and C5-methylcytosine (m5C) have been most broadly studied in different RNA viruses. Various studies, however, reported varying results with regard to number and extent of the modification. Here we investigated the m5C methylome of SARS-CoV-2, and we reexamined reported m5C sites in HIV and MLV. Using a rigorous bisulfite-sequencing protocol and stringent data analysis, we found no evidence for the presence of m5C in these viruses. The data emphasize the necessity for optimizing experimental conditions and bioinformatic data analysis.

Dynamics of RNA m5C modification during brain development

Post-transcriptional RNA modifications have been recognized as key regulators of neuronal differentiation and synapse development in the mammalian brain. While distinct sets of 5-methylcytosine (m⁵C) modified mRNAs have been detected in neuronal cells and brain tissues, no study has been performed to characterize methylated mRNA profiles in the developing brain. Here, together with regular RNA-seq, we performed transcriptome-wide bisulfite sequencing to compare RNA cytosine methylation patterns in neural stem cells (NSCs), cortical neuronal cultures, and brain tissues at three postnatal stages. Among 501 m⁵C sites identified, approximately 6% are consistently methylated across all five conditions. Compared to m⁵C sites identified in NSCs, 96% of them were hypermethylated in neurons and enriched for genes involved in positive transcriptional regulation and axon extension. In addition, brains at the early postnatal stage demonstrated substantial changes in both RNA cytosine methylation and gene expression of RNA cytosine methylation readers, writers, and erasers. Furthermore, differentially methylated transcripts were significantly enriched for genes regulating synaptic plasticity. Altogether, this study provides a brain epitranscriptomic dataset as a new resource and lays the foundation for further investigations into the role of RNA cytosine methylation during brain development.

Folate regulates RNA m5C modification and translation in neural stem cells

BACKGROUND

Folate is an essential B-group vitamin and a key methyl donor with important biological functions including DNA methylation regulation. Normal neurodevelopment and physiology are sensitive to the cellular folate levels. Either deficiency or excess of folate may lead to neurological disorders. Recently, folate has been linked to tRNA cytosine-5 methylation (m5C) and translation in mammalian mitochondria. However, the influence of folate intake on neuronal mRNA m5C modification and translation remains largely unknown. Here, we provide transcriptome-wide landscapes of m5C modification in poly(A)-enriched RNAs together with mRNA transcription and translation profiles for mouse neural stem cells (NSCs) cultured in three different concentrations of folate.

RESULTS

NSCs cultured in three different concentrations of folate showed distinct mRNA methylation profiles. Despite uncovering only a few differentially expressed genes, hundreds of differentially translated genes were identified in NSCs with folate deficiency or supplementation. The differentially translated genes induced by low folate are associated with cytoplasmic translation and mitochondrial function, while the differentially translated genes induced by high folate are associated with increased neural stem cell proliferation. Interestingly, compared to total mRNAs, polysome mRNAs contained high levels of m5C. Furthermore, an integrative analysis indicated a transcript-specific relationship between RNA m5C methylation and mRNA translation efficiency.

CONCLUSIONS

Altogether, our study reports a transcriptome-wide influence of folate on mRNA m5C methylation and translation in NSCs and reveals a potential link between mRNA m5C methylation and mRNA translation.

Neuronal Depolarization Induced RNA m5C Methylation Changes in Mouse Cortical Neurons

Neuronal activity is accomplished via substantial changes in gene expression, which may be accompanied by post-transcriptional modifications including RNA cytosine-5 methylation (m5C). Despite several reports on the transcriptome profiling of activated neurons, the dynamics of neuronal mRNA m5C modification in response to environmental stimuli has not been explored. Here, we provide transcriptome-wide maps of m5C modification, together with gene expression profiles, for mouse cortical neurons at 0 h, 2 h, and 6 h upon membrane depolarization. Thousands of differentially expressed genes (DEGs) were identified during the neuronal depolarization process. In stimulated neurons, the majority of early response genes were found to serve as expression regulators of late response genes, which are involved in signaling pathways and diverse synaptic functions. With RNA bisulfite sequencing data, a union set of 439 m5C sites was identified with high confidence, and approximately 30% of them were shared by neurons at all three time points. Interestingly, over 41% of the m5C sites showed increased methylation upon neuronal activation and were enriched in transcripts coding for proteins with synaptic functions. In addition, a modest negative correlation was observed between RNA expression and methylation. In summary, our study provided dynamic transcriptome-wide landscapes of RNA m5C methylation in neurons, and revealed that mRNA m5C methylation is associated with the regulation of gene expression.

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.

RNA bisulfite sequencing reveals NSUN2-mediated suppression of epithelial differentiation in pancreatic cancer

Posttranscriptional modifications in RNA have been considered to contribute to disease pathogenesis and tumor progression. NOL1/NOP2/Sun domain family member 2 (NSUN2) is an RNA methyltransferase that promotes tumor progression in several cancers. Pancreatic cancer relapse inevitably occurs even in cases where primary tumors have been successfully treated. Associations of cancer progression due to reprogramming of the cancer methyl-metabolome and the cancer genome have been noted, but the effect of base modifications, namely 5-methylcytosine (m5C), in the transcriptome remains unclear. Aberrant regulation of 5-methylcytosine turnover in cancer may affect posttranscriptional modifications in coding and noncoding RNAs in disease pathogenesis. Mutations in NSUN2 have been reported as drivers of neurodevelopmental disorders in mice, and upregulated expression of NSUN2 in tumors of the breast, bladder, and pancreas has been reported. In this study, we conducted mRNA whole transcriptomic bisulfite sequencing to categorize NSUN2 target sites in the mRNA of human pancreatic cancer cells. We identified a total of 2829 frequent m5C sites in mRNA from pancreatic cancer cells. A total of 90.9% (2572/2829) of these m5C sites were mapped to annotated genes in autosomes and sex chromosomes X and Y. Immunohistochemistry staining confirmed that the NSUN2 expression was significantly upregulated in cancer lesions in the LSL-KrasG12D/+;Trp53fl/fl;Pdx1-Cre (KPC) spontaneous pancreatic cancer mouse model induced by Pdx1-driven Cre/lox system expressing mutant KrasG12D and p53 deletion. The in vitro phenotypic analysis of NSUN2 knockdown showed mild effects on pancreatic cancer cell 2D/3D growth, morphology and gemcitabine sensitivity in the early phase of tumorigenesis, but cumulative changes after multiple cell doubling passages over time were required for these mutations to accumulate. Syngeneic transplantation of NSUN2-knockdown KPC cells via subcutaneous injection showed decreased stromal fibrosis and restored differentiation of ductal epithelium in vivo. SIGNIFICANCE: Transcriptome-wide mRNA bisulfite sequencing identified candidate m5C sites of mRNAs in human pancreatic cancer cells. NSUN2-mediated m5C mRNA metabolism was observed in a mouse model of pancreatic cancer. NSUN2 regulates cancer progression and epithelial differentiation via mRNA methylation.

5-Methylcytosine profiles in mouse transcriptomes suggest the randomness of m5C formation catalyzed by RNA methyltransferase

OBJECTIVE

5-Methylcytosine (m⁵C) is a type of chemical modification on the nucleotides and is widespread in both DNA and RNA. Although the DNA m⁵C has been extensively studied over the past years, the distribution and biological function of RNA m⁵C still remain to be elucidated. Here, I explored the profiles of RNA m⁵C in four mouse tissues by applying a RNA cytosine methylation data analysis tool to public mouse RNA m⁵C data.

RESULTS

I found that the methylation rates of cytosine were the same with the averages of methylation level at single-nucleotide level. Furthermore, I gave a mathematical formula to describe the observed relationship and analyzed it deeply. The sufficient necessary condition for the given formula suggests that the methylation levels at most m⁵C sites are the same in four mouse tissues. Therefore, I proposed a hypothesis that the m⁵C formation catalyzed by RNA methyltransferase is random and with the same probability at most m⁵C sites, which is the methylation rate of cytosine. My hypothesis can be used to explain the observed profiles of RNA m⁵C in four mouse tissues and will be benefit to future studies of the distribution and biological function of RNA m⁵C in mammals.

m5C-Atlas: a comprehensive database for decoding and annotating the 5-methylcytosine (m5C) epitranscriptome

5-Methylcytosine (m5C) is one of the most prevalent covalent modifications on RNA. It is known to regulate a broad variety of RNA functions, including nuclear export, RNA stability and translation. Here, we present m5C-Atlas, a database for comprehensive collection and annotation of RNA 5-methylcytosine. The database contains 166 540 m5C sites in 13 species identified from 5 base-resolution epitranscriptome profiling technologies. Moreover, condition-specific methylation levels are quantified from 351 RNA bisulfite sequencing samples gathered from 22 different studies via an integrative pipeline. The database also presents several novel features, such as the evolutionary conservation of a m5C locus, its association with SNPs, and any relevance to RNA secondary structure. All m5C-atlas data are accessible through a user-friendly interface, in which the m5C epitranscriptomes can be freely explored, shared, and annotated with putative post-transcriptional mechanisms (e.g. RBP intermolecular interaction with RNA, microRNA interaction and splicing sites). Together, these resources offer unprecedented opportunities for exploring m5C epitranscriptomes. The m5C-Atlas database is freely accessible at https://www.xjtlu.edu.cn/biologicalsciences/m5c-atlas.

RNA 5-methylcytosine regulates YBX2-dependent liquid-liquid phase separation

5-Methylcytosine (m5C) is one of the most prevalent internal modifications of messenger RNA (mRNA) in higher eukaryotes. Here we report that Y box protein 2 (YBX2) serves as a novel mammalian m5C binding protein to undergo liquid-liquid phase separation (LLPS) both in vivo and in vitro, and this YBX2-dependent LLPS is enhanced by m5C marked RNA. Furthermore, the crystal structure assay revealed that W100, as a distinct m5C binding site of YBX2, is critical in mediating YBX2 phase separation. Our study resolved the relationship between RNA m5C and phase separation, providing a clue for a new regulatory layer of epigenetics.

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

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

HIV Modifies the m6A and m5C Epitranscriptomic Landscape of the Host Cell

The study of RNA modifications, today known as epitranscriptomics, is of growing interest. The N6-methyladenosine (m6A) and 5-methylcytosine (m5C) RNA modifications are abundantly present on mRNA molecules, and impact RNA interactions with other proteins or molecules, thereby affecting cellular processes, such as RNA splicing, export, stability, and translation. Recently m6A and m5C marks were found to be present on human immunodeficiency (HIV) transcripts as well and affect viral replication. Therefore, the discovery of RNA methylation provides a new layer of regulation of HIV expression and replication, and thus offers novel array of opportunities to inhibit replication. However, no study has been performed to date to investigate the impact of HIV replication on the transcript methylation level in the infected cell. We used a productive HIV infection model, consisting of the CD4+ SupT1 T cell line infected with a VSV-G pseudotyped HIVeGFP-based vector, to explore the temporal landscape of m6A and m5C epitranscriptomic marks upon HIV infection, and to compare it to mock-treated cells. Cells were collected at 12, 24, and 36 h post-infection for mRNA extraction and FACS analysis. M6A RNA modifications were investigated by methylated RNA immunoprecipitation followed by high-throughput sequencing (MeRIP-Seq). M5C RNA modifications were investigated using a bisulfite conversion approach followed by high-throughput sequencing (BS-Seq). Our data suggest that HIV infection impacted the methylation landscape of HIV-infected cells, inducing mostly increased methylation of cellular transcripts upon infection. Indeed, differential methylation (DM) analysis identified 59 m6A hypermethylated and only 2 hypomethylated transcripts and 14 m5C hypermethylated transcripts and 7 hypomethylated ones. All data and analyses are also freely accessible on an interactive web resource (http://sib-pc17.unil.ch/HIVmain.html). Furthermore, both m6A and m5C methylations were detected on viral transcripts and viral particle RNA genomes, as previously described, but additional patterns were identified. This work used differential epitranscriptomic analysis to identify novel players involved in HIV life cycle, thereby providing innovative opportunities for HIV regulation.

Alteration of mRNA 5-Methylcytosine Modification in Neurons After OGD/R and Potential Roles in Cell Stress Response and Apoptosis

Epigenetic modifications play an important role in central nervous system disorders. As a widespread posttranscriptional RNA modification, the role of the m5C modification in cerebral ischemia-reperfusion injury (IRI) remains poorly defined. Here, we successfully constructed a neuronal oxygen-glucose deprivation/reoxygenation (OGD/R) model and obtained an overview of the transcriptome-wide m5C profiles using RNA-BS-seq. We discovered that the distribution of neuronal m5C modifications was highly conserved, significantly enriched in CG-rich regions and concentrated in the mRNA translation initiation regions. After OGD/R, modification level of m5C increased, whereas the number of methylated mRNA genes decreased. The amount of overlap of m5C sites with the binding sites of most RNA-binding proteins increased significantly, except for that of the RBM3-binding protein. Moreover, hypermethylated genes in neurons were significantly enriched in pathological processes, and the hub hypermethylated genes RPL8 and RPS9 identified by the protein-protein interaction network were significantly related to cerebral injury. Furthermore, the upregulated transcripts with hypermethylated modification were enriched in the processes involved in response to stress and regulation of apoptosis, and these processes were not identified in hypomethylated transcripts. In final, we verified that OGD/R induced neuronal apoptosis in vitro using TUNEL and western blot assays. Our study identified novel m5C mRNAs associated with ischemia-reperfusion in neurons, providing valuable perspectives for future studies on the role of the RNA methylation in cerebral IRI.

Quantitative and Single-Nucleotide Resolution Profiling of RNA 5-Methylcytosine

RNA has coevolved with numerous posttranscriptional modifications to sculpt interactions with proteins and other molecules. One of these modifications is 5-methylcytosine (m⁵C) and mapping the position and quantifying the level in different types of cellular RNAs and tissues is an important objective in the field of epitranscriptomics. Both in plants and animals bisulfite conversion has long been the gold standard for detection of m⁵C in DNA but it can also be applied to RNA. Here, we detail methods for highly reproducible bisulfite treatment of RNA, efficient locus-specific PCR amplification, detection of candidate sites by sequencing on the Illumina MiSeq platform, and bioinformatic calling of non-converted sites.

Episo: quantitative estimation of RNA 5-methylcytosine at isoform level by high-throughput sequencing of RNA treated with bisulfite

Motivation

RNA 5-methylcytosine (m⁵C) is a type of post-transcriptional modification that may be involved in numerous biological processes and tumorigenesis. RNA m⁵C can be profiled at single-nucleotide resolution by high-throughput sequencing of RNA treated with bisulfite (RNA-BisSeq). However, the exploration of transcriptome-wide profile and potential function of m⁵C in splicing remains to be elucidated due to lack of isoform level m⁵C quantification tool.

Results

We developed a computational package to quantify Epitranscriptomal RNA m⁵C at the transcript isoform level (named Episo). Episo consists of three tools: mapper, quant and Bisulfitefq, for mapping, quantifying and simulating RNA-BisSeq data, respectively. The high accuracy of Episo was validated using an improved m⁵C-specific methylated RNA immunoprecipitation (meRIP) protocol, as well as a set of in silico experiments. By applying Episo to public human and mouse RNA-BisSeq data, we found that the RNA m⁵C is not evenly distributed among the transcript isoforms, implying the m⁵C may subject to be regulated at isoform level.

Availability and implementation

Episo is released under the GNU GPLv3+ license. The resource code Episo is freely accessible from https://github.com/liujunfengtop/Episo (with Tophat/cufflink) and https://github.com/liujunfengtop/Episo/tree/master/Episo_Kallisto (with Kallisto).

Supplementary information

Supplementary data are available at Bioinformatics online.

Bioinformatics approaches for deciphering the epitranscriptome: Recent progress and emerging topics

Post-transcriptional RNA modification occurs on all types of RNA and plays a vital role in regulating every aspect of RNA function. Thanks to the development of high-throughput sequencing technologies, transcriptome-wide profiling of RNA modifications has been made possible. With the accumulation of a large number of high-throughput datasets, bioinformatics approaches have become increasing critical for unraveling the epitranscriptome. We review here the recent progress in bioinformatics approaches for deciphering the epitranscriptomes, including epitranscriptome data analysis techniques, RNA modification databases, disease-association inference, general functional annotation, and studies on RNA modification site prediction. We also discuss the limitations of existing approaches and offer some future perspectives.

Disease Activity-Associated Alteration of mRNA m5 C Methylation in CD4+ T Cells of Systemic Lupus Erythematosus

Epigenetic processes including RNA methylation, post-translational modifications, and non-coding RNA expression have been associated with the heritable risks of systemic lupus erythematosus (SLE). In this study, we aimed to explore the dysregulated expression of 5-methylcytosine (m⁵C) in CD4⁺ T cells from patients with SLE and the potential function of affected mRNAs in SLE pathogenesis. mRNA methylation profiles were ascertained through chromatography-coupled triple quadrupole mass spectrometry in CD4⁺ T cells from two pools of patients with SLE exhibiting stable activity, two pools with moderate-to-major activity, and two pools of healthy controls (HCs). Simultaneously, mRNA methylation profiles and expression profiling were performed using RNA-Bis-Seq and RNA-Seq, respectively. Integrated mRNA methylation and mRNA expression bioinformatics analysis was comprehensively performed. mRNA methyltransferase NSUN2 expression was validated in CD4⁺ T cells from 27 patients with SLE and 28 HCs using real-time polymerase chain reaction and western blot analyses. Hypomethylated-mRNA profiles of NSUN2-knockdown HeLa cells and of CD4⁺ T cells of patients with SLE were jointly analyzed using bioinformatics. Eleven methylation modifications (including elevated Am, 3′OMeA, m¹A, and m⁶A and decreased Ψ, m³C, m¹G, m⁵U, and t⁶A levels) were detected in CD4⁺ T cells of patients with SLE. Additionally, decreased m⁵C levels, albeit increased number of m⁵C-containing mRNAs, were observed in CD4⁺ T cells of patients with SLE compared with that in CD4⁺ T cells of HCs. m⁵C site distribution in mRNA transcripts was highly conserved and enriched in mRNA translation initiation sites. In particular, hypermethylated m⁵C or/and significantly up-regulated genes in SLE were significantly involved in immune-related and inflammatory pathways, including immune system, cytokine signaling pathway, and interferon signaling. Compared to that in HCs, NSUN2 expression was significantly lower in SLE CD4⁺ T cells. Notably, hypomethylated m⁵C genes in SLE and in NSUN2-knockdown HeLa cells revealed linkage between eukaryotic translation elongation and termination, and mRNA metabolism. Our study identified novel aberrant m⁵C mRNAs relevant to critical immune pathways in CD4⁺ T cells from patients with SLE. These data provide valuable perspectives for future studies of the multifunctionality and post-transcriptional significance of mRNA m⁵C modification in SLE.

The human mitochondrial 12S rRNA m4C methyltransferase METTL15 is required for mitochondrial function

Mitochondrial DNA gene expression is coordinately regulated both pre- and post-transcriptionally, and its perturbation can lead to human pathologies. Mitochondrial rRNAs (mt-rRNAs) undergo a series of nucleotide modifications after release from polycistronic mitochondrial RNA precursors, which is essential for mitochondrial ribosomal biogenesis. Cytosine N ⁴-methylation (m⁴C) at position 839 (m⁴C839) of the 12S small subunit mt-rRNA was identified decades ago; however, its biogenesis and function have not been elucidated in detail. Here, using several approaches, including immunofluorescence, RNA immunoprecipitation and methylation assays, and bisulfite mapping, we demonstrate that human methyltransferase-like 15 (METTL15), encoded by a nuclear gene, is responsible for 12S mt-rRNA methylation at m⁴C839 both in vivo and in vitro We tracked the evolutionary history of RNA m⁴C methyltransferases and identified a difference in substrate preference between METTL15 and its bacterial ortholog rsmH. Additionally, unlike the very modest impact of a loss of m⁴C methylation in bacterial small subunit rRNA on the ribosome, we found that METTL15 depletion results in impaired translation of mitochondrial protein-coding mRNAs and decreases mitochondrial respiration capacity. Our findings reveal that human METTL15 is required for mitochondrial function, delineate the evolution of methyltransferase substrate specificities and modification patterns in rRNA, and highlight a differential impact of m⁴C methylation on prokaryotic ribosomes and eukaryotic mitochondrial ribosomes.

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.

m7GHub: deciphering the location, regulation and pathogenesis of internal mRNA N7-methylguanosine (m7G) sites in human

Motivation

Recent progress in N7-methylguanosine (m7G) RNA methylation studies has focused on its internal (rather than capped) presence within mRNAs. Tens of thousands of internal mRNA m7G sites have been identified within mammalian transcriptomes, and a single resource to best share, annotate and analyze the massive m7G data generated recently are sorely needed.

Results

We report here m7GHub, a comprehensive online platform for deciphering the location, regulation and pathogenesis of internal mRNA m7G. The m7GHub consists of four main components, including: the first internal mRNA m7G database containing 44 058 experimentally validated internal mRNA m7G sites, a sequence-based high-accuracy predictor, the first web server for assessing the impact of mutations on m7G status, and the first database recording 1218 disease-associated genetic mutations that may function through regulation of m7G methylation. Together, m7GHub will serve as a useful resource for research on internal mRNA m7G modification.

Availability and implementation

m7GHub is freely accessible online at www.xjtlu.edu.cn/biologicalsciences/m7ghub.

Contact

kunqi.chen@liverpool.ac.uk

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

Advances in methods and software for RNA cytosine methylation analysis

Our understanding of RNA modifications has been growing rapidly over the last decade. Epitranscriptomics has recently emerged as an exciting, new field for understanding the fundamental mechanisms underlying RNA modifications and their impact on gene expression. Among the over one hundred different kinds of RNA modifications, cytosine methylation in mRNA (5-mrC) is now recognized as an important epigenetic mark that modulates mRNA transportation, translation, and stability at the post-transcriptional level. Across plant and animal species, recent studies have revealed the roles of mRNA cytosine methylation in several fundamental biological processes. In mammals, genome-wide profiling has determined thousands of mRNA transcripts carrying the 5-mrC modification in a tissue specific manner. Here, we summarize the experimental techniques that were exploited to determine 5-mrC in mRNA and the computational procedures implemented for RNA bisulfite sequencing data analysis.

Depletion of TRDMT1 affects 5-methylcytosine modification of mRNA and inhibits HEK293 cell proliferation and migration

Human TRDMT1 is a transfer RNA (tRNA) methyltransferase for cytosine-5 methylation and has been suggested to be involved in the regulation of numerous developmental processes. However, little is known about the molecular mechanisms or their biological significance. In this study, we investigated the effects of CRISPR-based TRDMT1 knockdown on phenotypes, mRNA m5C modifications and gene expression changes in HEK293 cells. We found that knockdown of TRDMT1 significantly inhibited cell proliferation and migration but had no effect on clonogenic potential. The inhibitory effects could be attenuated by re-expression of TRDMT1 in HEK293 cells. RNA sequencing (RNA-Seq) and RNA bisulfite sequencing (RNA-BisSeq) were performed in TRDMT1 knockdown and wild-type HEK293 cells. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses indicated that the differentially expressed genes were associated with the cell cycle, RNA transport, and RNA degradation and were enriched in cancer and Notch signaling pathways. We also found that TRDMT1 knockdown could change mRNA methylation levels. For the first time, these findings clarify the role of TRDMT1 in regulating mRNA methylation and inhibiting the proliferation and migration of HEK293 cells. These results provide new insights into a new function of TRDMT1 and elucidate the molecular mechanisms of aberrant RNA m5C during tumorigenesis.

RNA 5-Methylcytosine Facilitates the Maternal-to-Zygotic Transition by Preventing Maternal mRNA Decay

The maternal-to-zygotic transition (MZT) is a conserved and fundamental process during which the maternal environment is converted to an environment of embryonic-driven development through dramatic reprogramming. However, how maternally supplied transcripts are dynamically regulated during MZT remains largely unknown. Herein, through genome-wide profiling of RNA 5-methylcytosine (m⁵C) modification in zebrafish early embryos, we found that m⁵C-modified maternal mRNAs display higher stability than non-m⁵C-modified mRNAs during MZT. We discovered that Y-box binding protein 1 (Ybx1) preferentially recognizes m⁵C-modified mRNAs through π-π interactions with a key residue, Trp45, in Ybx1's cold shock domain (CSD), which plays essential roles in maternal mRNA stability and early embryogenesis of zebrafish. Together with the mRNA stabilizer Pabpc1a, Ybx1 promotes the stability of its target mRNAs in an m⁵C-dependent manner. Our study demonstrates an unexpected mechanism of RNA m⁵C-regulated maternal mRNA stabilization during zebrafish MZT, highlighting the critical role of m⁵C mRNA modification in early development.

5-methylcytosine promotes pathogenesis of bladder cancer through stabilizing mRNAs

Although 5-methylcytosine (m⁵C) is a widespread modification in RNAs, its regulation and biological role in pathological conditions (such as cancer) remain unknown. Here, we provide the single-nucleotide resolution landscape of messenger RNA m⁵C modifications in human urothelial carcinoma of the bladder (UCB). We identify numerous oncogene RNAs with hypermethylated m⁵C sites causally linked to their upregulation in UCBs and further demonstrate YBX1 as an m⁵C 'reader' recognizing m⁵C-modified mRNAs through the indole ring of W65 in its cold-shock domain. YBX1 maintains the stability of its target mRNA by recruiting ELAVL1. Moreover, NSUN2 and YBX1 are demonstrated to drive UCB pathogenesis by targeting the m⁵C methylation site in the HDGF 3' untranslated region. Clinically, a high coexpression of NUSN2, YBX1 and HDGF predicts the poorest survival. Our findings reveal an unprecedented mechanism of RNA m⁵C-regulated oncogene activation, providing a potential therapeutic strategy for UCB.

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.

Bisulfite Sequencing of RNA for Transcriptome-Wide Detection of 5-Methylcytosine

A powerful method to determine the methylation status of specific cytosine residues within RNA is bisulfite sequencing. In combination with high-throughput sequencing methods cytosine methylation can be determined at nucleotide resolution on a transcriptome-wide level. Nevertheless, several critical aspects need to be considered before starting such a project. Below we describe a detailed step-by-step protocol for planning and performing a transcriptome-wide bisulfite sequencing experiment and subsequent data analysis to determine methyl-cytosine in poly(A)RNA from cells and tissues.

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.

Statistically robust methylation calling for whole-transcriptome bisulfite sequencing reveals distinct methylation patterns for mouse RNAs

Cytosine-5 RNA methylation plays an important role in several biologically and pathologically relevant processes. However, owing to methodological limitations, the transcriptome-wide distribution of this mark has remained largely unknown. We previously established RNA bisulfite sequencing as a method for the analysis of RNA cytosine-5 methylation patterns at single-base resolution. More recently, next-generation sequencing has provided opportunities to establish transcriptome-wide maps of this modification. Here, we present a computational approach that integrates tailored filtering and data-driven statistical modeling to eliminate many of the artifacts that are known to be associated with bisulfite sequencing. By using RNAs from mouse embryonic stem cells, we performed a comprehensive methylation analysis of mouse tRNAs, rRNAs, and mRNAs. Our approach identified all known methylation marks in tRNA and two previously unknown but evolutionary conserved marks in 28S rRNA. In addition, mRNAs were found to be very sparsely methylated or not methylated at all. Finally, the tRNA-specific activity of the DNMT2 methyltransferase could be resolved at single-base resolution, which provided important further validation. Our approach can be used to profile cytosine-5 RNA methylation patterns in many experimental contexts and will be important for understanding the function of cytosine-5 RNA methylation in RNA biology and in human disease.

Distinct 5-methylcytosine profiles in poly(A) RNA from mouse embryonic stem cells and brain

BACKGROUND

Recent work has identified and mapped a range of posttranscriptional modifications in mRNA, including methylation of the N6 and N1 positions in adenine, pseudouridylation, and methylation of carbon 5 in cytosine (m5C). However, knowledge about the prevalence and transcriptome-wide distribution of m5C is still extremely limited; thus, studies in different cell types, tissues, and organisms are needed to gain insight into possible functions of this modification and implications for other regulatory processes.

RESULTS

We have carried out an unbiased global analysis of m5C in total and nuclear poly(A) RNA of mouse embryonic stem cells and murine brain. We show that there are intriguing differences in these samples and cell compartments with respect to the degree of methylation, functional classification of methylated transcripts, and position bias within the transcript. Specifically, we observe a pronounced accumulation of m5C sites in the vicinity of the translational start codon, depletion in coding sequences, and mixed patterns of enrichment in the 3' UTR. Degree and pattern of methylation distinguish transcripts modified in both embryonic stem cells and brain from those methylated in either one of the samples. We also analyze potential correlations between m5C and micro RNA target sites, binding sites of RNA binding proteins, and N6-methyladenosine.

CONCLUSION

Our study presents the first comprehensive picture of cytosine methylation in the epitranscriptome of pluripotent and differentiated stages in the mouse. These data provide an invaluable resource for future studies of function and biological significance of m5C in mRNA in mammals.

Transcriptome-Wide Detection of 5-Methylcytosine by Bisulfite Sequencing

While low-throughput RNA bisulfite sequencing is the method of choice to assess the methylation status of specific cytosines in candidate RNAs, the combination of bisulfite treatment of RNA with today's high-throughput sequencing techniques opens the door to methylation studies at nucleotide resolution on a transcriptome-wide scale. Below we describe a protocol for the transcriptome-wide analysis of total or fractionated poly(A)RNA in cells and tissues. Although the nature of the bisulfite sequencing protocol makes it comparably easy to translate from a low to a high-throughput approach, several critical points require attention before starting such a project. We describe a step-by-step protocol for planning and performing the experiment and analyzing the data.

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.