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

Regulatory roles of RNA methylation in vascular lesions in ocular and cardiopulmonary diseases

RNA methylation is a widespread regulatory mechanism that controls gene expression in physiological processes. In recent years, the mechanisms and functions of RNA methylation under diseased conditions have been increasingly unveiled by RNA sequencing technologies with large scale and high resolution. In this review, the fundamental concept of RNA methylation is introduced, and the common types of transcript methylation and their machineries are described. Then, the regulatory roles of RNA methylation, particularly N6-methyladenosine and 5-methylcytosine, in the vascular lesions of ocular and cardiopulmonary diseases are discussed and compared. The ocular diseases include corneal neovascularization, retinopathy of prematurity, diabetic retinopathy, and pathologic myopia; whereas the cardiopulmonary ailments involve atherosclerosis and pulmonary hypertension. This review hopes to shed light on the common regulatory mechanisms underlying the vascular lesions in these ocular and cardiopulmonary diseases, which may be conducive to developing therapeutic strategies in clinical practice.

Epigenetic addition of m5C to HBV transcripts promotes viral replication and evasion of innate antiviral responses

Eukaryotic five-methylcytosine (m5C) is an important regulator of viral RNA splicing, stability, and translation. However, its role in HBV replication remains largely unknown. In this study, functional m5C sites are identified in hepatitis B virus (HBV) mRNA. The m5C modification at nt 1291 is not only indispensable for Aly/REF export factor (ALYREF) recognition to promote viral mRNA export and HBx translation but also for the inhibition of RIG-I binding to suppress interferon-β (IFN-β) production. Moreover, NOP2/Sun RNA methyltransferase 2 (NSUN2) catalyzes the addition of m5C to HBV mRNA and is transcriptionally downregulated by the viral protein HBx, which suppresses the binding of EGR1 to the NSUN2 promoter. Additionally, NSUN2 expression correlates with m5C modification of type I IFN mRNA in host cells, thus, positively regulating IFN expression. Hence, the delicate regulation of NSUN2 expression induces m5C modification of HBV mRNA while decreasing the levels of m5C in host IFN mRNA, making it a vital component of the HBV life cycle. These findings provide new molecular insights into the mechanism of HBV-mediated IFN inhibition and may inform the development of new IFN-α based therapies.

Underexplored reciprocity between genome-wide methylation status and long non-coding RNA expression reflected in breast cancer research: potential impacts for the disease management in the framework of 3P medicine

Breast cancer (BC) is the most common female malignancy reaching a pandemic scale worldwide. A comprehensive interplay between genetic alterations and shifted epigenetic regions synergistically leads to disease development and progression into metastatic BC. DNA and histones methylations, as the most studied epigenetic modifications, represent frequent and early events in the process of carcinogenesis. To this end, long non-coding RNAs (lncRNAs) are recognized as potent epigenetic modulators in pathomechanisms of BC by contributing to the regulation of DNA, RNA, and histones' methylation. In turn, the methylation status of DNA, RNA, and histones can affect the level of lncRNAs expression demonstrating the reciprocity of mechanisms involved. Furthermore, lncRNAs might undergo methylation in response to actual medical conditions such as tumor development and treated malignancies. The reciprocity between genome-wide methylation status and long non-coding RNA expression levels in BC remains largely unexplored. Since the bio/medical research in the area is, per evidence, strongly fragmented, the relevance of this reciprocity for BC development and progression has not yet been systematically analyzed. Contextually, the article aims at:consolidating the accumulated knowledge on both-the genome-wide methylation status and corresponding lncRNA expression patterns in BC andhighlighting the potential benefits of this consolidated multi-professional approach for advanced BC management. Based on a big data analysis and machine learning for individualized data interpretation, the proposed approach demonstrates a great potential to promote predictive diagnostics and targeted prevention in the cost-effective primary healthcare (sub-optimal health conditions and protection against the health-to-disease transition) as well as advanced treatment algorithms tailored to the individualized patient profiles in secondary BC care (effective protection against metastatic disease). Clinically relevant examples are provided, including mitochondrial health control and epigenetic regulatory mechanisms involved.

The RNA m5C Methylase NSUN2 Modulates Corneal Epithelial Wound Healing

Purpose

The emerging epitranscriptomics offers insights into the physiopathological roles of various RNA modifications. The RNA methylase NOP2/Sun domain family member 2 (NSUN2) catalyzes 5-methylcytosine (m5C) modification of mRNAs. However, the role of NSUN2 in corneal epithelial wound healing (CEWH) remains unknown. Here we describe the functional mechanisms of NSUN2 in mediating CEWH.

Methods

RT-qPCR, Western blot, dot blot, and ELISA were used to determine the NSUN2 expression and overall RNA m5C level during CEWH. NSUN2 silencing or overexpression was performed to explore its involvement in CEWH both in vivo and in vitro. Multi-omics was integrated to reveal the downstream target of NSUN2. MeRIP-qPCR, RIP-qPCR, and luciferase assay, as well as in vivo and in vitro functional assays, clarified the molecular mechanism of NSUN2 in CEWH.

Results

The NSUN2 expression and RNA m5C level increased significantly during CEWH. NSUN2 knockdown significantly delayed CEWH in vivo and inhibited human corneal epithelial cells (HCEC) proliferation and migration in vitro, whereas NSUN2 overexpression prominently enhanced HCEC proliferation and migration. Mechanistically, we found that NSUN2 increased ubiquitin-like containing PHD and RING finger domains 1 (UHRF1) translation through the binding of RNA m5C reader Aly/REF export factor. Accordingly, UHRF1 knockdown significantly delayed CEWH in vivo and inhibited HCEC proliferation and migration in vitro. Furthermore, UHRF1 overexpression effectively rescued the inhibitory effect of NSUN2 silencing on HCEC proliferation and migration.

Conclusions

NSUN2-mediated m5C modification of UHRF1 mRNA modulates CEWH. This finding highlights the critical importance of this novel epitranscriptomic mechanism in control of CEWH.

The emergent role of mitochondrial RNA modifications in metabolic alterations

Mitochondrial epitranscriptomics refers to the modifications occurring in all the different RNA types of mitochondria. Although the number of mitochondrial RNA modifications is less than those in cytoplasm, substantial evidence indicates that they play a critical role in accurate protein synthesis. Recent evidence supported those modifications in mitochondrial RNAs also have crucial implications in mitochondrial-related diseases. In the light of current knowledge about the involvement, the association between mitochondrial RNA modifications and diseases arises from studies focusing on mutations in both mitochondrial and nuclear DNA genes encoding enzymes involved in such modifications. Here, we review the current evidence available for mitochondrial RNA modifications and their role in metabolic disorders, and we also explore the possibility of using them as promising targets for prevention and early detection. Finally, we discuss future directions of mitochondrial epitranscriptomics in these metabolic alterations, and how these RNA modifications may offer a new diagnostic and theragnostic avenue for preventive purposes. This article is categorized under: RNA Processing > RNA Editing and Modification.

Advances and Trends in Omics Technology Development

The human history has witnessed the rapid development of technologies such as high-throughput sequencing and mass spectrometry that led to the concept of “omics” and methodological advancement in systematically interrogating a cellular system. Yet, the ever-growing types of molecules and regulatory mechanisms being discovered have been persistently transforming our understandings on the cellular machinery. This renders cell omics seemingly, like the universe, expand with no limit and our goal toward the complete harness of the cellular system merely impossible. Therefore, it is imperative to review what has been done and is being done to predict what can be done toward the translation of omics information to disease control with minimal cell perturbation. With a focus on the “four big omics,” i.e., genomics, transcriptomics, proteomics, metabolomics, we delineate hierarchies of these omics together with their epiomics and interactomics, and review technologies developed for interrogation. We predict, among others, redoxomics as an emerging omics layer that views cell decision toward the physiological or pathological state as a fine-tuned redox balance.

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.

m5CRegpred: Epitranscriptome Target Prediction of 5-Methylcytosine (m5C) Regulators Based on Sequencing Features

5-methylcytosine (m5C) is a common post-transcriptional modification observed in a variety of RNAs. m5C has been demonstrated to be important in a variety of biological processes, including RNA structural stability and metabolism. Driven by the importance of m5C modification, many projects focused on the m5C sites prediction were reported before. To better understand the upstream and downstream regulation of m5C, we present a bioinformatics framework, m5CRegpred, to predict the substrate of m5C writer NSUN2 and m5C readers YBX1 and ALYREF for the first time. After features comparison, window lengths selection and algorism comparison on the mature mRNA model, our model achieved AUROC scores 0.869, 0.724 and 0.889 for NSUN2, YBX1 and ALYREF, respectively in an independent test. Our work suggests the substrate of m5C regulators can be distinguished and may help the research of m5C regulators in a special condition, such as substrates prediction of hyper- or hypo-expressed m5C regulators in human disease.

Biological roles of RNA m5C modification and its implications in Cancer immunotherapy

Epigenetics including DNA and RNA modifications have always been the hotspot field of life sciences in the post-genome era. Since the first mapping of N6-methyladenosine (m⁶A) and the discovery of its widespread presence in mRNA, there are at least 160-170 RNA modifications have been discovered. These methylations occur in different RNA types, and their distribution is species-specific. 5-methylcytosine (m⁵C) has been found in mRNA, rRNA and tRNA of representative organisms from all kinds of species. As reversible epigenetic modifications, m⁵C modifications of RNA affect the fate of the modified RNA molecules and play important roles in various biological processes including RNA stability control, protein synthesis, and transcriptional regulation. Furthermore, accumulative evidence also implicates the role of RNA m⁵C in tumorigenesis. Here, we review the latest progresses in the biological roles of m⁵C modifications and how it is regulated by corresponding “writers”, “readers” and “erasers” proteins, as well as the potential molecular mechanism in tumorigenesis and cancer immunotherapy.

A mark of disease: how mRNA modifications shape genetic and acquired pathologies

RNA modifications have recently emerged as a widespread and complex facet of gene expression regulation. Counting more than 170 distinct chemical modifications with far-reaching implications for RNA fate, they are collectively referred to as the epitranscriptome. These modifications can occur in all RNA species, including messenger RNAs (mRNAs) and noncoding RNAs (ncRNAs). In mRNAs the deposition, removal, and recognition of chemical marks by writers, erasers and readers influence their structure, localization, stability, and translation. In turn, this modulates key molecular and cellular processes such as RNA metabolism, cell cycle, apoptosis, and others. Unsurprisingly, given their relevance for cellular and organismal functions, alterations of epitranscriptomic marks have been observed in a broad range of human diseases, including cancer, neurological and metabolic disorders. Here, we will review the major types of mRNA modifications and editing processes in conjunction with the enzymes involved in their metabolism and describe their impact on human diseases. We present the current knowledge in an updated catalog. We will also discuss the emerging evidence on the crosstalk of epitranscriptomic marks and what this interplay could imply for the dynamics of mRNA modifications. Understanding how this complex regulatory layer can affect the course of human pathologies will ultimately lead to its exploitation toward novel epitranscriptomic therapeutic 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.

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.

Autophagy Induced by Simian Retrovirus Infection Controls Viral Replication and Apoptosis of Jurkat T Lymphocytes

Autophagy and apoptosis are two important evolutionarily conserved host defense mechanisms against viral invasion and pathogenesis. However, the association between the two pathways during the viral infection of T lymphocytes remains to be elucidated. Simian type D retrovirus (SRV) is an etiological agent of fatal simian acquired immunodeficiency syndrome (SAIDS), which can display disease features that are similar to acquired immunodeficiency syndrome in humans. In this study, we demonstrate that infection with SRV-8, a newly isolated subtype of SRV, triggered both autophagic and apoptotic pathways in Jurkat T lymphocytes. Following infection with SRV-8, the autophagic proteins LC3 and p62/SQSTM1 interacted with procaspase-8, which might be responsible for the activation of the caspase-8/-3 cascade and apoptosis in SRV-8-infected Jurkat cells. Our findings indicate that autophagic responses to SRV infection of T lymphocytes promote the apoptosis of T lymphocytes, which, in turn, might be a potential pathogenetic mechanism for the loss of T lymphocytes during SRV infection.

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.

Messenger RNA Modifications in Plants

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

Methods for RNA Modification Mapping Using Deep Sequencing: Established and New Emerging Technologies

New analytics of post-transcriptional RNA modifications have paved the way for a tremendous upswing of the biological and biomedical research in this field. This especially applies to methods that included RNA-Seq techniques, and which typically result in what is termed global scale modification mapping. In this process, positions inside a cell’s transcriptome are receiving a status of potential modification sites (so called modification calling), typically based on a score of some kind that issues from the particular method applied. The resulting data are thought to represent information that goes beyond what is contained in typical transcriptome data, and hence the field has taken to use the term “epitranscriptome”. Due to the high rate of newly published mapping techniques, a significant number of chemically distinct RNA modifications have become amenable to mapping, albeit with variegated accuracy and precision, depending on the nature of the technique. This review gives a brief overview of known techniques, and how they were applied to modification calling.

trumpet: transcriptome-guided quality assessment of m6A-seq data

BACKGROUND

Methylated RNA immunoprecipitation sequencing (MeRIP-seq or m6A-seq) has been extensively used for profiling transcriptome-wide distribution of RNA N6-Methyl-Adnosine methylation. However, due to the intrinsic properties of RNA molecules and the intricate procedures of this technique, m6A-seq data often suffer from various flaws. A convenient and comprehensive tool is needed to assess the quality of m6A-seq data to ensure that they are suitable for subsequent analysis.

RESULTS

From a technical perspective, m6A-seq can be considered as a combination of ChIP-seq and RNA-seq; hence, by effectively combing the data quality assessment metrics of the two techniques, we developed the trumpet R package for evaluation of m6A-seq data quality. The trumpet package takes the aligned BAM files from m6A-seq data together with the transcriptome information as the inputs to generate a quality assessment report in the HTML format.

CONCLUSIONS

The trumpet R package makes a valuable tool for assessing the data quality of m6A-seq, and it is also applicable to other fragmented RNA immunoprecipitation sequencing techniques, including m1A-seq, CeU-Seq, Ψ-seq, etc.