Nucleotide-Level Profiling of m⁵C RNA Methylation

Comparative RNA Genomics

Over the last quarter of a century it has become clear that RNA is much more than just a boring intermediate in protein expression. Ancient RNAs still appear in the core information metabolism and comprise a surprisingly large component in bacterial gene regulation. A common theme with these types of mostly small RNAs is their reliance of conserved secondary structures. Large-scale sequencing projects, on the other hand, have profoundly changed our understanding of eukaryotic genomes. Pervasively transcribed, they give rise to a plethora of large and evolutionarily extremely flexible non-coding RNAs that exert a vastly diverse array of molecule functions. In this chapter we provide a-necessarily incomplete-overview of the current state of comparative analysis of non-coding RNAs, emphasizing computational approaches as a means to gain a global picture of the modern RNA world.

m5C Regulator-Mediated Methylation Modification Patterns and Tumor Microenvironment Infiltration Characterization in Papillary Thyroid Carcinoma

Recently, immune response modulation at the epigenetic level is illustrated in studies, but the possible function of RNA 5-methylcytosine (m⁵C) modification in cell infiltration within the tumor microenvironment (TME) is still unclear. Three different m⁵C modification patterns were identified, and high differentiation degree was observed in the cell infiltration features within TME under the above three identified patterns. A low m⁵C-score, which was reflected in the activated immunity, predicted the relatively favorable prognostic outcome. A small amount of effective immune infiltration was seen in the high m⁵C-score subtype, indicating the dismal patient survival. Our study constructed a diagnostic model using the 10 signature genes highly related to the m⁵C-score, discovered that the model exhibited high diagnostic accuracy for PTC, and screened out five potential drugs for PTC based on this m5C-score model. m⁵C modification exerts an important part in forming the TME complexity and diversity. It is valuable to evaluate the m⁵C modification patterns in single tumors, so as to enhance our understanding towards the infiltration characterization in TME.

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.

Novel long noncoding RNA NMR promotes tumor progression via NSUN2 and BPTF in esophageal squamous cell carcinoma

Long noncoding RNAs (lncRNA) have been implicated in cancer but most of them remain largely unstudied. Here, we identified a novel NSUN2 methylated lncRNA (NMR), which was significantly upregulated in esophageal squamous cell carcinoma (ESCC), functioned as a key regulator of ESCC tumor metastasis and drug resistance. Upregulation of NMR correlated with tumor metastasis and indicated poor overall survival in ESCC patients. Functionally, NMR could promote tumor cell migration and invasion, inhibit cisplatin-induced apoptosis and increase drug resistance in ESCC cells. Mechanistically, transcription of NMR could be upregulated by NF-κB activation after IL-1β and TNF-α treatment. NMR was methylated by NSUN2 and might competitively inhibit methylation of potential mRNAs. NMR could directly bind to chromatin regulator BPTF, and potentially promote MMP3 and MMP10 expression by ERK1/2 pathway through recruiting BPTF to chromatin. Taken together, NMR functions as an oncogenic gene and may serve as new biomarker and therapeutic target in ESCC.

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.

Comparative RNA Genomics

Over the last two decades it has become clear that RNA is much more than just a boring intermediate in protein expression. Ancient RNAs still appear in the core information metabolism and comprise a surprisingly large component in bacterial gene regulation. A common theme with these types of mostly small RNAs is their reliance of conserved secondary structures. Large scale sequencing projects, on the other hand, have profoundly changed our understanding of eukaryotic genomes. Pervasively transcribed, they give rise to a plethora of large and evolutionarily extremely flexible noncoding RNAs that exert a vastly diverse array of molecule functions. In this chapter we provide a-necessarily incomplete-overview of the current state of comparative analysis of noncoding RNAs, emphasizing computational approaches as a means to gain a global picture of the modern RNA world.

The emerging epitranscriptomics of long noncoding RNAs

The pervasive transcription of genomes into long noncoding RNAs has been amply demonstrated in recent years and garnered much attention. Similarly, emerging 'epitranscriptomics' research has shown that chemically modified nucleosides, thought to be largely the domain of tRNAs and other infrastructural RNAs, are far more widespread and can exert unexpected influence on RNA utilization. Both areas are characterized by the often-ephemeral nature of the subject matter in that few individual examples have been fully assessed for their molecular or cellular function, and effects might often be subtle and cumulative. Here we review available information at the intersection of these two exciting areas of biology, by focusing on four RNA modifications that have been mapped transcriptome-wide: 5-methylcytidine, N6-methyladenosine, pseudouridine as well as adenosine to inosine (A-to-I) editing, and their incidence and function in long noncoding RNAs. This article is part of a Special Issue entitled: Clues to long noncoding RNA taxonomy, edited by Dr. Tetsuro Hirose and Dr. Shinichi Nakagawa.