Transcriptome-wide Mapping of Internal N7-Methylguanosine Methylome in Mammalian mRNA

Nanopore-Based Detection of Viral RNA Modifications

The chemical modification of ribonucleotides plays an integral role in the biology of diverse viruses and their eukaryotic host cells. Mapping the precise identity, location, and abundance of modified ribonucleotides remains a key goal of many studies aimed at characterizing the function and importance of a given modification. While mapping of specific RNA modifications through short-read sequencing approaches has powered a wealth of new discoveries in the past decade, this approach is limited by inherent biases and an absence of linkage information. Moreover, in viral contexts, the challenge is increased due to the compact nature of viral genomes giving rise to many overlapping transcript isoforms that cannot be adequately resolved using short-read sequencing approaches. The recent emergence of nanopore sequencing, specifically the ability to directly sequence native RNAs from virus-infected host cells, provides not just a new methodology for mapping modified ribonucleotides but also a new conceptual framework for what can be derived from the resulting sequencing data. In this minireview, we provide a detailed overview of how nanopore direct RNA sequencing works, the computational approaches applied to identify modified ribonucleotides, and the core concepts underlying both. We further highlight recent studies that have applied this approach to interrogating viral biology and finish by discussing key experimental considerations and how we predict that these methodologies will continue to evolve.

An Effective Deep Learning-Based Architecture for Prediction of N7-Methylguanosine Sites in Health Systems

N7-methylguanosine (m7G) is one of the most important epigenetic modifications found in rRNA, mRNA, and tRNA, and performs a promising role in gene expression regulation. Owing to its significance, well-equipped traditional laboratory-based techniques have been performed for the identification of N7-methylguanosine (m7G). Consequently, these approaches were found to be time-consuming and cost-ineffective. To move on from these traditional approaches to predict N7-methylguanosine sites with high precision, the concept of artificial intelligence has been adopted. In this study, an intelligent computational model called N7-methylguanosine-Long short-term memory (m7G-LSTM) is introduced for the prediction of N7-methylguanosine sites. One-hot encoding and word2vec feature schemes are used to express the biological sequences while the LSTM and CNN algorithms have been employed for classification. The proposed “m7G-LSTM” model obtained an accuracy value of 95.95%, a specificity value of 95.94%, a sensitivity value of 95.97%, and Matthew’s correlation coefficient (MCC) value of 0.919. The proposed predictive m7G-LSTM model has significantly achieved better outcomes than previous models in terms of all evaluation parameters. The proposed m7G-LSTM computational system aims to support the drug industry and help researchers in the fields of bioinformatics to enhance innovation for the prediction of the behavior of N7-methylguanosine sites.

N7-Methylguanosine Regulatory Genes Profoundly Affect the Prognosis, Progression, and Antitumor Immune Response of Hepatocellular Carcinoma

Background

Hepatocellular carcinoma (HCC) is a common abdominal cancer with poor survival outcomes. Although there is growing evidence that N7-methylguanosine (m7G) is closely associated with tumor prognosis, development, and immune response, few studies focus on this topic.

Methods

The novel m7G risk signature was constructed through the Lasso regression analysis. Its prognostic value was evaluated through a series of survival analyses and was tested in ICGC-LIRI, GSE14520, and GSE116174 cohorts. CIBERSORT, ssGSEA, and ESTIMATE methods were applied to explore the effects of the m7G risk score on tumor immune microenvironment (TIM). The GSEA method was used to evaluate the impacts of the m7G risk score on glycolysis, ferroptosis, and pyroptosis. The human protein atlas (HPA) database was used to clarify the histological expression levels of five m7G signature genes. The biofunctions of NCBP2 in hepatocellular cancer (HC) cells were confirmed through qPCR, CCK8, and transwell assays.

Results

Five m7G regulatory genes comprised the novel risk signature. The m7G risk score was identified as an independent prognostic factor of HCC and could increase the decision-making benefit of traditional prognostic models. Besides, we established a nomogram containing the clinical stage and m7G risk score to predict the survival rates of HCC patients. The prognostic value of the m7G model was successfully validated in ICGC and GSE116174 cohorts. Moreover, high m7G risk led to a decreased infiltration level of CD8+ T cells, whereas it increased the infiltration levels of Tregs and macrophages. The glycolysis and pyroptosis processes were found to be enriched in the HCC patients with high m7G risk. Finally, overexpression of NCBP2 could promote the proliferation, migration, and invasion of HC cells.

Conclusions

The m7G risk score was closely related to the prognosis, antitumor immune process, glycolysis, and malignant progression of HCC. NCBP2 has pro-oncogenic abilities, showing promise as a novel treatment target.

Construction of a Novel Prognostic Model in Lung Adenocarcinoma Based on 7-Methylguanosine-Related Gene Signatures

Increasing evidence has implicated the modification of 7-methylguanosine (m7G), a type of RNA modification, in tumor progression. However, no comprehensive analysis to date has summarized the predicted role of m7G-related gene signatures in lung adenocarcinoma (LUAD). Herein, we aimed to develop a novel prognostic model in LUAD based on m7G-related gene signatures. The LUAD transcriptome profiling data and corresponding clinical data were acquired from the Cancer Genome Atlas (TCGA) and two Gene Expression Omnibus datasets. After screening, we first obtained 29 m7G-related genes, most of which were upregulated in tumor tissues and negatively associated with overall survival (OS). According to the expression similarity of m7G-related genes, the combined samples from the TCGA-LUAD and GSE68465 datasets were further classified as two clusters that exhibit distinct OS rates and genetic heterogeneity. Then, we constructed a novel prognostic model involving four genes by using 130 differentially expressed genes among the two clusters. The combined samples were randomly divided into a training cohort and an internal validation cohort in a 1:1 ratio, and the GSE72094 dataset was used as an external validation cohort. The samples were divided into high- and low-risk groups. We demonstrated that a higher risk score was an independent negative prognostic factor and predicted poor OS. A nomogram was further constructed to better predict the survival of LUAD patients. Functional enrichment analyses indicated that cell cycle and DNA replication-related biological processes and pathways were enriched in the high-risk group. More importantly, the low-risk group had greater infiltration and enrichment of most immune cells, as well as higher ESTIMATE, immune, and stromal scores. In addition, the high-risk group had a lower TIDE score and higher expressions of most immune checkpoint-related genes. We finally noticed that patients in the high-risk group were more sensitive to chemotherapeutic agents commonly used in LUAD. In conclusion, we herein summarized for the first time the alterations and prognostic role of m7G-related genes in LUAD and then constructed a prognostic model based on m7G-related gene signatures that could accurately and stably predict survival and guide individualized treatment decision-making in LUAD patients.

Rapid Determination of RNA Modifications in Consensus Motifs by Nuclease Protection with Ion-Tagged Oligonucleotide Probes and Matrix-Assisted Laser Desorption Ionization Mass Spectrometry

The reversible and substoichiometric modification of RNA has recently emerged as an additional layer of translational regulation in normal biological function and disease. Modifications are often enzymatically deposited in and removed from short (~5 nt) consensus motif sequences to carefully control the translational output of the cell. Although characterization of modification occupancy at consensus motifs can be accomplished using RNA sequencing methods, these approaches are generally time-consuming and do not directly detect post-transcriptional modifications. Here, we present a nuclease protection assay coupled with matrix-assisted laser desorption ionization mass spectrometry (MALDI-MS) to rapidly characterize modifications in consensus motifs, such as GGACU, which frequently harbor N6-methyladenosine (m6A). While conventional nuclease protection methods rely on long (~30 nt) oligonucleotide probes that preclude the global assessment of consensus motif modification stoichiometry, we investigated a series of ion-tagged oligonucleotide (ITO) probes and found that a benzylimidazolium-functionalized ITO (ABzIM-ITO) conferred significantly improved nuclease resistance for GGACU targets. After optimizing the conditions of the nuclease protection assay, we applied the ITO and MALDI-MS-based method for determining the stoichiometry of GG(m6A)CU and GGACU in RNA mixtures. Overall, the ITO-based nuclease protection and MALDI-MS method constitutes a rapid and promising approach for determining modification stoichiometries of consensus motifs.

Chemical modifications to mRNA nucleobases impact translation elongation and termination

Messenger RNAs (mRNAs) serve as blueprints for protein synthesis by the molecular machine the ribosome. The ribosome relies on hydrogen bonding interactions between adaptor aminoacyl-transfer RNA molecules and mRNAs to ensure the rapid and faithful translation of the genetic code into protein. There is a growing body of evidence suggesting that chemical modifications to mRNA nucleosides impact the speed and accuracy of protein synthesis by the ribosome. Modulations in translation rates have downstream effects beyond protein production, influencing protein folding and mRNA stability. Given the prevalence of such modifications in mRNA coding regions, it is imperative to understand the consequences of individual modifications on translation. In this review we present the current state of our knowledge regarding how individual mRNA modifications influence ribosome function. Our comprehensive comparison of the impacts of 16 different mRNA modifications on translation reveals that most modifications can alter the elongation step in the protein synthesis pathway. Additionally, we discuss the context dependence of these effects, highlighting the necessity of further study to uncover the rules that govern how any given chemical modification in an mRNA codon is read by the ribosome.

Characterization of RNA modifications in gastric cancer to identify prognosis-relevant gene signatures

BACKGROUND

Most human genes have diverse transcript isoforms, which mainly arise from alternative cleavage and polyadenylation (APA) at 3' ends. N7-methylguanosine (m⁷ G) is also an essential epigenetic modification at the 5' end. However, the contribution of these two RNA modifications to the development, prognosis, regulation mechanisms, and drug sensitivity of gastric cancer (GC) is unclear.

METHODS

The expression data of 2412 patients were extracted from 12 cohorts and the RNA modification patterns of 20 marker genes were systematically identified into phenotypic clusters using the unsupervised clustering approach. Following that, we developed an RNA modification model (RMscore) to quantify each GC patient's RNA modification index. Finally, we examined the correlation between RMscore and clinical features such as survival outcomes, molecular subtypes identified by the Asian Cancer Research Group (ACRG), posttranscriptional regulation, and chemotherapeutic sensitivity in GC.

RESULTS

The samples were categorized into two groups on the basis of their RMscore: high and low. The group with a low RMscore had a bad prognosis. Moreover, the low RMscore was associated with KRAS, Hedgehog, EMT, and TGF-β signaling, whereas a high RMscore was related to abnormal cell cycle signaling pathway activation. The findings also revealed that the RMscore contributes to the regulation of the miRNA-mRNA network. Drug sensitivity analysis revealed that RMscore is associated with the response to some anticancer drugs.

CONCLUSIONS

The RMscore model has the potential to be a useful tool for prognosis prediction in patients with GC. A comprehensive investigation of APA-RNA and m⁷ G-RNA modifications may reveal novel insights into the epigenetics of GC and aid in the development of more effective treatment strategies.

Emerging Role of Epitranscriptomics in Diabetes Mellitus and Its Complications

Diabetes mellitus (DM) and its related complications are among the leading causes of disability and mortality worldwide. Substantial studies have explored epigenetic regulation that is involved in the modifications of DNA and proteins, but RNA modifications in diabetes are still poorly investigated. In recent years, posttranscriptional epigenetic modification of RNA (the so-called 'epitranscriptome') has emerged as an interesting field of research. Numerous modifications, mainly N6 -methyladenosine (m6A), have been identified in nearly all types of RNAs and have been demonstrated to have an indispensable effect in a variety of human diseases, such as cancer, obesity, and diabetes. Therefore, it is particularly important to understand the molecular basis of RNA modifications, which might provide a new perspective for the pathogenesis of diabetes mellitus and the discovery of new therapeutic targets. In this review, we aim to summarize the recent progress in the epitranscriptomics involved in diabetes and diabetes-related complications. We hope to provide some insights for enriching the understanding of the epitranscriptomic regulatory mechanisms of this disease as well as the development of novel therapeutic targets for future clinical benefit.

The Progression of N6-methyladenosine Study and Its Role in Neuropsychiatric Disorders

Epitranscriptomic modifications can affect every aspect of RNA biology, including stability, transport, splicing, and translation, participate in global intracellular mRNA metabolism, and regulate gene expression and a variety of biological processes. N6-methyladenosine (m6A) as the most prevalent modification contributes to normal embryonic brain development and memory formation. However, changes in the level of m6A modification and the expression of its related proteins cause abnormal nervous system functions, including brain tissue development retardation, axon regeneration disorders, memory changes, and neural stem cell renewal and differentiation disorders. Recent studies have revealed that m6A modification and its related proteins play key roles in the development of various neuropsychiatric disorders, such as depression, Alzheimer’s disease, and Parkinson’s disease. In this review, we summarize the research progresses of the m6A modification regulation mechanism in the central nervous system and discuss the effects of gene expression regulation mediated by m6A modification on the biological functions of the neuropsychiatric disorders, thereby providing some insight into new research targets and treatment directions for human diseases.

Identification and Verification of m7G Modification Patterns and Characterization of Tumor Microenvironment Infiltration via Multi-Omics Analysis in Clear Cell Renal Cell Carcinoma

The epigenetic modification of tumorigenesis and progression in neoplasm has been demonstrated in recent studies. Nevertheless, the underlying association of N7-methylguanosine (m⁷G) regulation with molecular heterogeneity and tumor microenvironment (TME) in clear cell renal cell carcinoma (ccRCC) remains unknown. We explored the expression profiles and genetic variation features of m⁷G regulators and identified their correlations with patient outcomes in pan-cancer. Three distinct m⁷G modification patterns, including MGCS1, MGCS2, and MGCS3, were further determined and systematically characterized via multi-omics data in ccRCC. Compared with the other two subtypes, patients in MGCS3 exhibited a lower clinical stage/grade and better prognosis. MGCS1 showed the lowest enrichment of metabolic activities. MGCS2 was characterized by the suppression of immunity. We then established and validated a scoring tool named m7Sig, which could predict the prognosis of ccRCC patients. This study revealed that m⁷G modification played a vital role in the formation of the tumor microenvironment in ccRCC. Evaluating the m⁷G modification landscape helps us to raise awareness and strengthen the understanding of ccRCC’s characterization and, furthermore, to guide future clinical decision making.

The potential role of N7-methylguanosine (m7G) in cancer

N⁷-methylguanosine (m7G), one of the most prevalent RNA modifications, has recently attracted significant attention. The m7G modification actively participates in biological and pathological functions by affecting the metabolism of various RNA molecules, including messenger RNA, ribosomal RNA, microRNA, and transfer RNA. Increasing evidence indicates a critical role for m7G in human disease development, especially cancer, and aberrant m7G levels are closely associated with tumorigenesis and progression via regulation of the expression of multiple oncogenes and tumor suppressor genes. Currently, the underlying molecular mechanisms of m7G modification in cancer are not comprehensively understood. Here, we review the current knowledge regarding the potential function of m7G modifications in cancer and discuss future m7G-related diagnostic and therapeutic strategies.

i5hmCVec: Identifying 5-Hydroxymethylcytosine Sites of Drosophila RNA Using Sequence Feature Embeddings

5-Hydroxymethylcytosine (5hmC), one of the most important RNA modifications, plays an important role in many biological processes. Accurately identifying RNA modification sites helps understand the function of RNA modification. In this work, we propose a computational method for identifying 5hmC-modified regions using machine learning algorithms. We applied a sequence feature embedding method based on the dna2vec algorithm to represent the RNA sequence. The results showed that the performance of our model is better that of than state-of-art methods. All dataset and source codes used in this study are available at: https://github.com/liu-h-y/5hmC_model.

The functional roles of m6A modification in T lymphocyte responses and autoimmune diseases

RNA N6-methyladenosine (m⁶A) modification is abundant in eukaryotes, bacteria and archaea. It is an RNA modification mainly existing in messenger RNA (mRNAs) and has a significant effect on the metabolism and function of mRNAs. m⁶A modification is controlled by three types of proteins, namely methyltransferase as the "writers", demethylase as the "erasers", and specific m⁶A recognized protein (YTHDF1-3) as the "readers". Recent studies have shown that m⁶A modification plays an important role in cancer, viral infection and autoimmune diseases. In this review, we will elaborate on the m⁶A modifications in the homeostasis and differentiation of T cells. Then we will further summarize the effects of m⁶A modification on the T cell responses and T cell-mediated autoimmune diseases. This will advance T cell epigenetics research and provide potential biomarkers and therapeutic targets for autoimmune diseases.

A Comprehensive Analysis of METTL1 to Immunity and Stemness in Pan-Cancer

Background

Previous studies have reported the effect of N⁷-methylguanosine (m⁷G) regulator methyltransferase like-1 protein (METTL1) in tumor initiation, metastasis, and chemosensitivity. However, the relationship between METTL1 and cancer immune infiltration is not validated and the prognostic significance of METTL1 in pan-cancer remains unclear.

Methods

Clinical parameters, including gender, age, lifetime, stage, and treatment response were analyzed to evaluate the prognostic significance of METTL1. To evaluate protein level of METTL1, the METTL1 activity was generated by single sample gene set enrichment analysis. The one-class logistic regression algorithm was used to calculate the stemness indices based on transcriptomics and methylation data of pan-cancer and pluripotent stem cells. The relationship between METTL1 expression or activity and tumor immune infiltration were analyzed to explore the significance of METTL1 in tumor immunotherapy. Meanwhile, the correlation between three immunotherapeutic biomarkers and METTL1 was investigated. Finally, to calculate the association between drug sensitivity and METTL1 expression, spearman correlation analysis was performed.

Results

METTL1 was not intimately related to gender, age, tumor stage, or treatment outcome of the various cancers, but it displayed potential prognostic significance for evaluating patient survival. High METTL1 expression was related to tumor progression-relevant pathways. Moreover, METTL1 exhibited a distinct correlation with tumor immune microenvironment infiltration and stemness indices. In the anti-PD-L1 cohort, patients in treatment response group exhibited significantly higher METTL1 expression than those in the no/limited response group. Further analysis showed that tumor cell lines with higher METTL1 expression were more sensitive to drugs targeting chromatin histone methylation, ERK-MAPK and WNT signaling pathways.

Conclusion

This study provides insight into the correlation of METTL1 with tumor immune infiltration and stemness in pan-cancer, revealing the significance of METTL1 for cancer progression and guiding more effective and generalized therapy strategies.

Advances in COVID-19 mRNA vaccine development

To date, the coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has determined 399,600,607 cases and 5,757,562 deaths worldwide. COVID-19 is a serious threat to human health globally. The World Health Organization (WHO) has declared COVID-19 pandemic a major public health emergency. Vaccination is the most effective and economical intervention for controlling the spread of epidemics, and consequently saving lives and protecting the health of the population. Various techniques have been employed in the development of COVID-19 vaccines. Among these, the COVID-19 messenger RNA (mRNA) vaccine has been drawing increasing attention owing to its great application prospects and advantages, which include short development cycle, easy industrialization, simple production process, flexibility to respond to new variants, and the capacity to induce better immune response. This review summarizes current knowledge on the structural characteristics, antigen design strategies, delivery systems, industrialization potential, quality control, latest clinical trials and real-world data of COVID-19 mRNA vaccines as well as mRNA technology. Current challenges and future directions in the development of preventive mRNA vaccines for major infectious diseases are also discussed.

Challenges with Simulating Modified RNA: Insights into Role and Reciprocity of Experimental and Computational Approaches

RNA is critical to a broad spectrum of biological and viral processes. This functional diversity is a result of their dynamic nature; the variety of three-dimensional structures that they can fold into; and a host of post-transcriptional chemical modifications. While there are many experimental techniques to study the structural dynamics of biomolecules, molecular dynamics simulations (MDS) play a significant role in complementing experimental data and providing mechanistic insights. The accuracy of the results obtained from MDS is determined by the underlying physical models i.e., the force-fields, that steer the simulations. Though RNA force-fields have received a lot of attention in the last decade, they still lag compared to their protein counterparts. The chemical diversity imparted by the RNA modifications adds another layer of complexity to an already challenging problem. Insight into the effect of RNA modifications upon RNA folding and dynamics is lacking due to the insufficiency or absence of relevant experimental data. This review provides an overview of the state of MDS of modified RNA, focusing on the challenges in parameterization of RNA modifications as well as insights into relevant reference experiments necessary for their calibration.

N7-methylguanosine tRNA modification promotes esophageal squamous cell carcinoma tumorigenesis via the RPTOR/ULK1/autophagy axis

Mis-regulated RNA modifications promote the processing and translation of oncogenic mRNAs to facilitate cancer progression, while the molecular mechanisms remain unclear. Here we reveal that tRNA m⁷G methyltransferase complex proteins METTL1 and WDR4 are significantly up-regulated in esophageal squamous cell carcinoma (ESCC) tissues and associated with poor ESCC prognosis. In addition, METTL1 and WDR4 promote ESCC progression via the tRNA m⁷G methyltransferase activity in vitro and in vivo. Mechanistically, METTL1 or WDR4 knockdown leads to decreased expression of m⁷G-modified tRNAs and reduces the translation of a subset of oncogenic transcripts enriched in RPTOR/ULK1/autophagy pathway. Furthermore, ESCC models using Mettl1 conditional knockout and knockin mice uncover the essential function of METTL1 in promoting ESCC tumorigenesis in vivo. Our study demonstrates the important oncogenic function of mis-regulated tRNA m⁷G modification in ESCC, and suggest that targeting METTL1 and its downstream signaling axis could be a promising therapeutic target for ESCC treatment.

Deregulation of METTL1-mediated N7- methylguanosine tRNA modification can promote oncogenesis. Here, the authors report that this modification regulates the translation of proteins in both the mTOR and negative regulators of autophagy pathways, resulting in the progression of esophageal squamous cell carcinoma.

The Non-N6-Methyladenosine Epitranscriptome Patterns and Characteristics of Tumor Microenvironment Infiltration and Mesenchymal Transition in Glioblastoma

Background

An increasing number of RNA modification types other than N⁶-methyladenosine (m⁶A) modification have been detected. Nonetheless, the probable functions of RNA modifications beyond m⁶A in the tumor microenvironment (TME), mesenchymal (MES) transition, immunotherapy, and drug sensitivity remain unclear.

Methods

We analyzed the characteristics of 32 non-m⁶A RNA modification regulators in 539 glioblastoma (GBM) patients and the TME cell infiltration and MES transition patterns. Using principal component analysis, a non-m⁶A epitranscriptome regulator score (RM score) model was established. We estimated the association between RM score and clinical characteristics, TME status, GBM subtypes, and drug and immunotherapy response.

Results

Three definite non-m⁶A RNA modification patterns associated with diverse biological pathways and clinical characteristics were identified. The high RM score group was characterized by a poor prognosis, enhanced immune infiltration, and MES subtype. Further analysis indicated that the high RM score group had a lower tumor mutation burden as well as a weaker response to immunotherapy. The higher RM score group may benefit more from drugs targeting the EGFR and WNT signaling pathways.

Conclusion

Our study exposed the potential relationship of non-m⁶A RNA modification regulators with clinical features, TME status, and GBM subtype and clarified its therapeutic value.

The epitranscriptome toolbox

In the last decade, the notion that mRNA modifications are involved in regulation of gene expression was demonstrated in thousands of studies. To date, new technologies and methods allow accurate identification, transcriptome-wide mapping, and functional characterization of a growing number of RNA modifications, providing important insights into the biology of these marks. Most of the methods and approaches were developed for studying m⁶A, the most prevalent internal mRNA modification. However, unique properties of other RNA modifications stimulated the development of additional approaches. In this technical primer, we will discuss the available tools and approaches for detecting and studying different RNA modifications.

Mapping the m1A, m5C, m6A and m7G methylation atlas in zebrafish brain under hypoxic conditions by MeRIP-seq

BACKGROUND

The epigenetic modifications play important regulatory roles in tissue development, maintenance of physiological functions and pathological process. RNA methylations, including newly identified m1A, m5C, m6A and m7G, are important epigenetic modifications. However, how these modifications are distributed in the transcriptome of vertebrate brains and whether their abundance is altered under pathological conditions are still poorly understood. In this study, we chose the model animal of zebrafish to conduct a systematic study to investigate the mRNA methylation atlas in the brain.

RESULTS

By performing unbiased analyses of the m1A, m5C, m6A and m7G methylation of mRNA, we found that within the whole brain transcriptome, with the increase of the gene expression levels, the overall level of each of these four modifications on the related genes was also progressively increased. Further bioinformatics analysis indicated that the zebrafish brain has an abundance of m1A modifications. In the hypoxia-treated zebrafish brains, the proportion of m1A is decreased, affecting the RNA splicing and zebrafish endogenous retroviruses.

CONCLUSIONS

Our study presents the first comprehensive atlas of m1A, m5C, m6A and m7G in the epitranscriptome of the zebrafish brain and reveals the distribution of these modifications in mRNA under hypoxic conditions. These data provide an invaluable resource for further research on the involvement of m1A, m5C, m6A and m7G in the regulation of miRNA and repeat elements in vertebrates, and provide new thoughts to study the brain hypoxic injury on the aspect of epitranscriptome.

Single-Base Resolution Mapping Reveals Distinct 5-Formylcytidine in Saccharomyces cerevisiae mRNAs

5-Formylcytidine (f⁵C) is one type of post-transcriptional RNA modification, which is known at the wobble position of tRNA in mitochondria and essential for mitochondrial protein synthesis. Here, we show a method to detect f⁵C modifications in RNA and a transcriptome-wide f⁵C mapping technique, named f⁵C-seq. It is developed based on the treatment of pyridine borane, which can reduce f⁵C to 5,6-dihydrouracil, thus inducing C-to-T transition in f⁵C sites during PCR to achieve single-base resolution detection. More than 1000 f⁵C sites were identified after mapping in Saccharomyces cerevisiae by f⁵C-seq. Moreover, codon composition demonstrated a preference for f⁵C within wobble sites in mRNA, suggesting the potential role in regulation of translation. These findings expand the scope of the understanding of cytosine modifications in mRNA.

The METTL5-TRMT112 N6-methyladenosine methyltransferase complex regulates mRNA translation via 18S rRNA methylation

Ribosomal RNAs (rRNAs) have long been known to carry chemical modifications, including 2'O-methylation, pseudouridylation, N⁶-methyladenosine (m⁶A), and N^6,6-dimethyladenosine. While the functions of many of these modifications are unclear, some are highly conserved and occur in regions of the ribosome critical for mRNA decoding. Both 28S rRNA and 18S rRNA carry single m⁶A sites, and while the methyltransferase ZCCHC4 has been identified as the enzyme responsible for the 28S rRNA m⁶A modification, the methyltransferase responsible for the 18S rRNA m⁶A modification has remained unclear. Here, we show that the METTL5-TRMT112 methyltransferase complex installs the m⁶A modification at position 1832 of human 18S rRNA. Our work supports findings that TRMT112 is required for METTL5 stability and reveals that human METTL5 mutations associated with microcephaly and intellectual disability disrupt this interaction. We show that loss of METTL5 in human cancer cell lines and in mice regulates gene expression at the translational level; additionally, Mettl5 knockout mice display reduced body size and evidence of metabolic defects. While recent work has focused heavily on m⁶A modifications in mRNA and their roles in mRNA processing and translation, we demonstrate here that deorphanizing putative methyltransferase enzymes can reveal previously unappreciated regulatory roles for m⁶A in noncoding RNAs.

N7-methylguanosine modification of lncRNAs in a rat model of hypoxic pulmonary hypertension: a comprehensive analysis

Background

Long non-coding RNAs (lncRNAs) play a critical role in the pathogenesis of hypoxic pulmonary hypertension (HPH). The role of N7-methylguanosine (m7G) modification in lncRNAs has received increased attentions in recent years. However, the m7G-methylation of lncRNA in HPH has yet to be determined. We have therefore performed a transcriptome-wide analysis of m7G lncRNAs in HPH.

Results

Differentially-expressed m7Gs were detected in HPH, and m7G lncRNAs were significantly upregulated compared with non-m7G lncRNAs in HPH. Importantly, this was the first time that the upregulated m7G lncXR_591973 and m7G lncXR_592398 were identified in HPH.

Conclusion

This study provides the first m7G transcriptome-wide analysis of HPH. Importantly, two HPH-associated m7G lncRNAs were identified, although their clinical significance requires further validation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08188-8.

Modulation of DNA/RNA Methylation Signaling Mediating Metabolic Homeostasis in Cancer

Nucleic acid methylation is a fundamental epigenetic mechanism that impinges upon several cellular attributes, including metabolism and energy production. The dysregulation of deoxyribonucleic acid (DNA)/ribonucleic acid (RNA) methylation can lead to metabolic rewiring in the cell, which in turn facilitates tumor development. Here, we review the current knowledge on the interplay between DNA/RNA methylation and metabolic programs in cancer cells. We also discuss the mechanistic role of these pathways in tumor development and progression.

RNA nucleotide methylation: 2021 update

Among RNA modifications, transfer of methylgroups from the typical cofactor S-adenosyl-l-methionine by methyltransferases (MTases) to RNA is by far the most common reaction. Since our last review about a decade ago, the field has witnessed the re-emergence of mRNA methylation as an important mechanism in gene regulation. Attention has then spread to many other RNA species; all being included into the newly coined concept of the "epitranscriptome." The focus moved from prokaryotes and single cell eukaryotes as model organisms to higher eukaryotes, in particular to mammals. The perception of the field has dramatically changed over the past decade. A previous lack of phenotypes in knockouts in single cell organisms has been replaced by the apparition of MTases in numerous disease models and clinical investigations. Major driving forces of the field include methylation mapping techniques, as well as the characterization of the various MTases, termed "writers." The latter term has spilled over from DNA modification in the neighboring epigenetics field, along with the designations "readers," applied to mediators of biological effects upon specific binding to a methylated RNA. Furthermore "eraser" enzymes effect the newly discovered oxidative removal of methylgroups. A sense of reversibility and dynamics has replaced the older perception of RNA modification as a concrete-cast, irreversible part of RNA maturation. A related concept concerns incompletely methylated residues, which, through permutation of each site, lead to inhomogeneous populations of numerous modivariants. This review recapitulates the major developments of the past decade outlined above, and attempts a prediction of upcoming trends. This article is categorized under: RNA Processing > RNA Editing and Modification.

Chemical methods and advanced sequencing technologies for deciphering mRNA modifications

RNA modification, like other epigenetic modifications such as DNA modification and histone modification, is an emerging player in the field of the posttranscriptional regulation of gene expression. More than 160 kinds of RNA modifications have been identified, and they are widely distributed in different types of RNA. Recently, researchers have increasingly used advanced technologies to study modified nucleic acids in order to elucidate their biological functions and expand the understanding of the central laws of epigenetics. In this tutorial review, we comprehensively outline current advanced techniques for decoding RNA modifications, highlighting some of the bottlenecks in existing approaches as well as new opportunities that may lead to innovations. With this review, we expect to provide chemistry and biology students and researchers with ideas for solving some challenging problems, such as how to simultaneously detect multiple types of modifications within the same system. Moreover, some low-coverage modifications that may act as 'candidates' in important transcriptional processes need to be further explored. These novel approaches have the potential to lay a foundation for understanding the nuanced complexities of the biological functions of RNA modification.

Balancing of mitochondrial translation through METTL8-mediated m3C modification of mitochondrial tRNAs

Mitochondria contain a specific translation machinery for the synthesis of mitochondria-encoded respiratory chain components. Mitochondrial tRNAs (mt-tRNAs) are also generated from the mitochondrial DNA and, similar to their cytoplasmic counterparts, are post-transcriptionally modified. Here, we find that the RNA methyltransferase METTL8 is a mitochondrial protein that facilitates 3-methyl-cytidine (m3C) methylation at position C32 of the mt-tRNASer(UCN) and mt-tRNAThr. METTL8 knockout cells show a reduction in respiratory chain activity, whereas overexpression increases activity. In pancreatic cancer, METTL8 levels are high, which correlates with lower patient survival and an enhanced respiratory chain activity. Mitochondrial ribosome profiling uncovered mitoribosome stalling on mt-tRNASer(UCN)- and mt-tRNAThr-dependent codons. Further analysis of the respiratory chain complexes using mass spectrometry revealed reduced incorporation of the mitochondrially encoded proteins ND6 and ND1 into complex I. The well-balanced translation of mt-tRNASer(UCN)- and mt-tRNAThr-dependent codons through METTL8-mediated m3C32 methylation might, therefore, facilitate the optimal composition and function of the mitochondrial respiratory chain.

METTL1 promotes hepatocarcinogenesis via m7 G tRNA modification-dependent translation control

BACKGROUND

N7 -methylguanosine (m7 G) modification is one of the most common transfer RNA (tRNA) modifications in humans. The precise function and molecular mechanism of m7 G tRNA modification in hepatocellular carcinoma (HCC) remain poorly understood.

METHODS

The prognostic value and expression level of m7 G tRNA methyltransferase complex components methyltransferase-like protein-1 (METTL1) and WD repeat domain 4 (WDR4) in HCC were evaluated using clinical samples and TCGA data. The biological functions and mechanisms of m7 G tRNA modification in HCC progression were studied in vitro and in vivo using cell culture, xenograft model, knockin and knockout mouse models. The m7 G reduction and cleavage sequencing (TRAC-seq), polysome profiling and polyribosome-associated mRNA sequencing methods were used to study the levels of m7 G tRNA modification, tRNA expression and mRNA translation efficiency.

RESULTS

The levels of METTL1 and WDR4 are elevated in HCC and associated with advanced tumour stages and poor patient survival. Functionally, silencing METTL1 or WDR4 inhibits HCC cell proliferation, migration and invasion, while forced expression of wild-type METTL1 but not its catalytic dead mutant promotes HCC progression. Knockdown of METTL1 reduces m7 G tRNA modification and decreases m7 G-modified tRNA expression in HCC cells. Mechanistically, METTL1-mediated tRNA m7 G modification promotes the translation of target mRNAs with higher frequencies of m7 G-related codons. Furthermore, in vivo studies with Mettl1 knockin and conditional knockout mice reveal the essential physiological function of Mettl1 in hepatocarcinogenesis using hydrodynamics transfection HCC model.

CONCLUSIONS

Our work reveals new insights into the role of the misregulated tRNA modifications in liver cancer and provides molecular basis for HCC diagnosis and treatment.

METTL1-m7 G-EGFR/EFEMP1 axis promotes the bladder cancer development

BACKGROUND

The posttranscriptional modifications of transfer RNA (tRNA) are critical for all aspects of the tRNA function and have been implicated in the tumourigenesis and progression of many human cancers. By contrast, the biological functions of methyltransferase-like 1 (METTL1)-regulated m7 G tRNA modification in bladder cancer (BC) remain obscure.

RESULTS

In this research, we show that METTL1 was highly expressed in BC, and its level was correlated with poor patient prognosis. Silencing METTL1 suppresses the proliferation, migration and invasion of BC cells in vitro and in vivo. Multi-omics analysis reveals that METTL1-mediated m7 G tRNA modification altered expression of certain target genes, including EGFR/EFEMP1. Mechanistically, METTL1 regulates the translation of EGFR/EFEMP1 via modifying certain tRNAs. Furthermore, forced expression of EGFR/EFEMP1 partially rescues the effect of METTL1 deletion on BC cells.

CONCLUSIONS

Our findings demonstrate the oncogenic role of METTL1 and the pathological significance of the METTL1-m7 G-EGFR/EFEMP1 axis in the BC development, thus providing potential therapeutic targets for the BC treatment.

The epitranscriptome of small non-coding RNAs

Small non-coding RNAs are short RNA molecules and involved in many biological processes, including cell proliferation and differentiation, immune response, cell death, epigenetic regulation, metabolic control. A diversity of RNA modifications have been identified in these small non-coding RNAs, including transfer RNAs (tRNAs), microRNAs (miRNAs), PIWI-interacting RNAs (piRNAs), small nuclear RNA (snRNA), small nucleolar RNAs (snoRNAs), and tRNA-derived small RNAs (tsRNAs). These post-transcriptional modifications are involved in the biogenesis and function of these small non-coding RNAs. In this review, we will summarize the existence of RNA modifications in the small non-coding RNAs and the emerging roles of these epitranscriptomic marks.

Technical challenges in defining RNA modifications

It is well established that DNA base modifications play a key role in gene regulation during development and in response to environmental stress. This type of epigenetic control of development and environmental responses has been intensively studied over the past few decades. Similar to DNA, various RNA species also undergo modifications that play important roles in, for example, RNA splicing, protein translation, and the avoidance of immune surveillance by host. More than 160 different types of RNA modifications have been identified. In addition to base modifications, RNA modification also involves splicing of pre-mRNAs, leading to as many as tens of transcript isoforms from a single pre-RNA, especially in higher organisms. However, the function, prevalence and distribution of RNA modifications are poorly understood. The lack of a suitable method for the reliable identification of RNA modifications constitutes a significant challenge to studying their functions. This review focuses on the technologies that enable de novo identification of RNA base modifications and the alternatively spliced mRNA transcripts.

RNA modifications as emerging therapeutic targets

The field of epitranscriptome, posttranscriptional modifications to RNAs, is still growing up and has presented substantial evidences for the role of RNA modifications in diseases. In terms of new drug development, RNA modifications have a great promise for therapy. For example, more than 170 type of modifications exist in various types of RNAs. Regulatory genes and their roles in critical biological process have been identified and they are associated with several diseases. Current data, for example, identification of inhibitors targeting RNA modifications regulatory genes, strongly support the idea that RNA modifications have potential as emerging therapeutic targets. Therefore, in this review, RNA modifications and regulatory genes were comprehensively documented in terms of drug development by summarizing the findings from previous studies. It was discussed how RNA modifications or regulatory genes can be targeted by altering molecular mechanisms. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Processing > RNA Editing and Modification.

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.

Inhibition of AlkB Nucleic Acid Demethylases: Promising New Epigenetic Targets

The AlkB family of nucleic acid demethylases is currently of intense chemical, biological, and medical interest because of its critical roles in several key cellular processes, including epigenetic gene regulation, RNA metabolism, and DNA repair. Emerging evidence suggests that dysregulation of AlkB demethylases may underlie the pathogenesis of several human diseases, particularly obesity, diabetes, and cancer. Hence there is strong interest in developing selective inhibitors for these enzymes to facilitate their mechanistic and functional studies and to validate their therapeutic potential. Herein we review the remarkable advances made over the past 20 years in AlkB demethylase inhibition research. We discuss the rational design of reported inhibitors, their mode-of-binding, selectivity, cellular activity, and therapeutic opportunities. We further discuss unexplored structural elements of the AlkB subfamilies and propose potential strategies to enable subfamily selectivity. It is hoped that this perspective will inspire novel inhibitor design and advance drug discovery research in this field.

The ribosome epitranscriptome: inert-or a platform for functional plasticity?

A universal property of all rRNAs explored to date is the prevalence of post-transcriptional ("epitranscriptional") modifications, which expand the chemical and topological properties of the four standard nucleosides. Are these modifications an inert, constitutive part of the ribosome? Or could they, in part, also regulate the structure or function of the ribosome? In this review, we summarize emerging evidence that rRNA modifications are more heterogeneous than previously thought, and that they can also vary from one condition to another, such as in the context of a cellular response or a developmental trajectory. We discuss the implications of these results and key open questions on the path toward connecting such heterogeneity with function.

Advances in RNA Epigenetic Modifications in Hepatocellular Carcinoma and Potential Targeted Intervention Strategies

The current interventions for hepatocellular carcinoma (HCC) are not satisfactory, and more precise targets and promising strategies need to be explored. Recent research has demonstrated the non-negligible roles of RNA epigenetic modifications such as N6-methyladenosine (m6A) and 5-methylcytosine (m5C) in various cancers, including HCC. However, the specific targeting mechanisms are not well elucidated. In this review, we focus on the occurrence and detailed physiopathological roles of multiple RNA modifications on diverse RNAs closely related to the HCC process. In particular, we highlight fresh insights into the impact mechanisms of these posttranscriptional modifications on the whole progression of HCC. Furthermore, we analyzed the possibilities and significance of these modifications and regulators as potential therapeutic targets in HCC treatment, which provides the foundation for exploring targeted intervention strategies. This review will propel the identification of promising therapeutic targets and novel strategies that can be translated into clinical applications for HCC treatment.

Identifying N7-methylguanosine sites by integrating multiple features

Recent studies reported that N7-methylguanosine (m7G) plays a vital role in gene expression regulation. As a consequence, determining the distribution of m7G is a crucial step towards further understanding its biological functions. Although biological experimental approaches are capable of accurately locating m7G sites, they are labor-intensive, costly, and time-consuming. Therefore, it is necessary to develop more effective and robust computational methods to replace, or at least complement current experimental methods. In this study, we developed a novel sequence-based computational tool to identify RNA m7G sites. In this model, 22 kinds of dinucleotide physicochemical (PC) properties were employed to encode the RNA sequence. Three types of descriptors, including auto-covariance, cross-covariance, and discrete wavelet transform were adopted to extract effective features from the PC matrix. The least absolute shrinkage and selection operator (LASSO) algorithm was utilized to reduce the influence of irrelevant or redundant features. Finally, these selected features were fed into a support vector machine (SVM) for distinguishing m7G from non-m7G sites. The proposed method significantly outperforms existing predictors across all evaluation metrics. It indicates that the approach is effective in identifying RNA m7G sites.

Strategies to Avoid Artifacts in Mass Spectrometry-Based Epitranscriptome Analyses

In this report, we perform structure validation of recently reported RNA phosphorothioate (PT) modifications, a new set of epitranscriptome marks found in bacteria and eukaryotes including humans. By comparing synthetic PT-containing diribonucleotides with native species in RNA hydrolysates by high-resolution mass spectrometry (MS), metabolic stable isotope labeling, and PT-specific iodine-desulfurization, we disprove the existence of PTs in RNA from E. coli, S. cerevisiae, human cell lines, and mouse brain. Furthermore, we discuss how an MS artifact led to the initial misidentification of 2'-O-methylated diribonucleotides as RNA phosphorothioates. To aid structure validation of new nucleic acid modifications, we present a detailed guideline for MS analysis of RNA hydrolysates, emphasizing how the chosen RNA hydrolysis protocol can be a decisive factor in discovering and quantifying RNA modifications in biological samples.

Aberrant RNA methylation triggers recruitment of an alkylation repair complex

Central to genotoxic responses is their ability to sense highly specific signals to activate the appropriate repair response. We previously reported that the activation of the ASCC-ALKBH3 repair pathway is exquisitely specific to alkylation damage in human cells. Yet the mechanistic basis for the selectivity of this pathway was not immediately obvious. Here, we demonstrate that RNA but not DNA alkylation is the initiating signal for this process. Aberrantly methylated RNA is sufficient to recruit ASCC, while an RNA dealkylase suppresses ASCC recruitment during chemical alkylation. In turn, recruitment of ASCC during alkylation damage, which is mediated by the E3 ubiquitin ligase RNF113A, suppresses transcription and R-loop formation. We further show that alkylated pre-mRNA is sufficient to activate RNF113A E3 ligase in vitro in a manner dependent on its RNA binding Zn-finger domain. Together, our work identifies an unexpected role for RNA damage in eliciting a specific response to genotoxins.

Metabolic turnover and dynamics of modified ribonucleosides by 13C labeling

Tandem mass spectrometry (MS/MS) is an accurate tool to assess modified ribonucleosides and their dynamics in mammalian cells. However, MS/MS quantification of lowly abundant modifications in non-ribosomal RNAs is unreliable, and the dynamic features of various modifications are poorly understood. Here, we developed a ¹³C labeling approach, called ¹³C-dynamods, to quantify the turnover of base modifications in newly transcribed RNA. This turnover-based approach helped to resolve mRNA from ncRNA modifications in purified RNA or free ribonucleoside samples and showed the distinct kinetics of the N6-methyladenosine (m⁶A) versus 7-methylguanosine (m⁷G) modification in polyA+-purified RNA. We uncovered that N6,N6-dimethyladenosine (m⁶₂A) exhibits distinct turnover in small RNAs and free ribonucleosides when compared to known m⁶₂A-modified large rRNAs. Finally, combined measurements of turnover and abundance of these modifications informed on the transcriptional versus posttranscriptional sensitivity of modified ncRNAs and mRNAs, respectively, to stress conditions. Thus, ¹³C-dynamods enables studies of the origin of modified RNAs at steady-state and subsequent dynamics under nonstationary conditions. These results open new directions to probe the presence and biological regulation of modifications in particular RNAs.

METTL3 Inhibitors for Epitranscriptomic Modulation of Cellular Processes

The methylase METTL3 is the writer enzyme of the N6 -methyladenosine (m6 A) modification of RNA. Using a structure-based drug discovery approach, we identified a METTL3 inhibitor with potency in a biochemical assay of 280 nM, while its enantiomer is 100 times less active. We observed a dose-dependent reduction in the m6 A methylation level of mRNA in several cell lines treated with the inhibitor already after 16 h of treatment, which lasted for at least 6 days. Importantly, the prolonged incubation (up to 6 days) with the METTL3 inhibitor did not alter levels of other RNA modifications (i. e., m1 A, m6 Am , m7 G), suggesting selectivity of the developed compound towards other RNA methyltransferases.

Machine learning applications in RNA modification sites prediction

Ribonucleic acid (RNA) modifications are post-transcriptional chemical composition changes that have a fundamental role in regulating the main aspect of RNA function. Recently, large datasets have become available thanks to the recent development in deep sequencing and large-scale profiling. This availability of transcriptomic datasets has led to increased use of machine learning based approaches in epitranscriptomics, particularly in identifying RNA modifications. In this review, we comprehensively explore machine learning based approaches used for the prediction of 11 RNA modification types, namely,m 1 A ,m 6 A ,m 5 C , 5 hmC , ψ , 2 ' - O - Me , ac 4 C ,m 7 G , A - to - I ,m 2 G , and D . This review covers the life cycle of machine learning methods to predict RNA modification sites including available benchmark datasets, feature extraction, and classification algorithms. We compare available methods in terms of datasets, target species, approach, and accuracy for each RNA modification type. Finally, we discuss the advantages and limitations of the reviewed approaches and suggest future perspectives.

General Principles for the Detection of Modified Nucleotides in RNA by Specific Reagents

Epitranscriptomics heavily rely on chemical reagents for the detection, quantification, and localization of modified nucleotides in transcriptomes. Recent years have seen a surge in mapping methods that use innovative and rediscovered organic chemistry in high throughput approaches. While this has brought about a leap of progress in this young field, it has also become clear that the different chemistries feature variegated specificity and selectivity. The associated error rates, e.g., in terms of false positives and false negatives, are in large part inherent to the chemistry employed. This means that even assuming technically perfect execution, the interpretation of mapping results issuing from the application of such chemistries are limited by intrinsic features of chemical reactivity. An important but often ignored fact is that the huge stochiometric excess of unmodified over-modified nucleotides is not inert to any of the reagents employed. Consequently, any reaction aimed at chemical discrimination of modified versus unmodified nucleotides has optimal conditions for selectivity that are ultimately anchored in relative reaction rates, whose ratio imposes intrinsic limits to selectivity. Here chemical reactivities of canonical and modified ribonucleosides are revisited as a basis for an understanding of the limits of selectivity achievable with chemical methods.

1,4,9-Triazaspiro[5.5]undecan-2-one Derivatives as Potent and Selective METTL3 Inhibitors

N⁶-methyladenosine (m⁶A) is the most frequent of the 160 RNA modifications reported so far. Accumulating evidence suggests that the METTL3/METTL14 protein complex, part of the m⁶A regulation machinery, is a key player in a variety of diseases including several types of cancer, type 2 diabetes, and viral infections. Here we report on a protein crystallography-based medicinal chemistry optimization of a METTL3 hit compound that has resulted in a 1400-fold potency improvement (IC₅₀ of 5 nM for the lead compound 22 (UZH2) in a time-resolved Förster resonance energy transfer (TR-FRET) assay). The series has favorable ADME properties as physicochemical characteristics were taken into account during hit optimization. UZH2 shows target engagement in cells and is able to reduce the m⁶A/A level of polyadenylated RNA in MOLM-13 (acute myeloid leukemia) and PC-3 (prostate cancer) cell lines.

BERT-m7G: A Transformer Architecture Based on BERT and Stacking Ensemble to Identify RNA N7-Methylguanosine Sites from Sequence Information

As one of the most prevalent posttranscriptional modifications of RNA, N7-methylguanosine (m7G) plays an essential role in the regulation of gene expression. Accurate identification of m7G sites in the transcriptome is invaluable for better revealing their potential functional mechanisms. Although high-throughput experimental methods can locate m7G sites precisely, they are overpriced and time-consuming. Hence, it is imperative to design an efficient computational method that can accurately identify the m7G sites. In this study, we propose a novel method via incorporating BERT-based multilingual model in bioinformatics to represent the information of RNA sequences. Firstly, we treat RNA sequences as natural sentences and then employ bidirectional encoder representations from transformers (BERT) model to transform them into fixed-length numerical matrices. Secondly, a feature selection scheme based on the elastic net method is constructed to eliminate redundant features and retain important features. Finally, the selected feature subset is input into a stacking ensemble classifier to predict m7G sites, and the hyperparameters of the classifier are tuned with tree-structured Parzen estimator (TPE) approach. By 10-fold cross-validation, the performance of BERT-m7G is measured with an ACC of 95.48% and an MCC of 0.9100. The experimental results indicate that the proposed method significantly outperforms state-of-the-art prediction methods in the identification of m7G modifications.

Evolution of methyltransferase like (METTL) proteins in Metazoa: A complex gene family involved in epitranscriptomic regulation and other epigenetic processes

The methyltransferase-like (METTL) proteins constitute a family of seven-beta-strand methyltransferases with S-adenosyl methionine-binding domains that modify DNA, RNA, and proteins. Methylation by METTL proteins contributes to the epigenetic, and in the case of RNA modifications, epitranscriptomic regulation of a variety of biological processes. Despite their functional importance, most investigations of the substrates and functions of METTLs within metazoans have been restricted to model vertebrate taxa. In the present work, we explore the evolutionary mechanisms driving the diversification and functional differentiation of 33 individual METTL proteins across Metazoa. Our results show that METTLs are nearly ubiquitous across the animal kingdom, with most having arisen early in metazoan evolution (i.e., occur in basal metazoan phyla). Individual METTL lineages each originated from single independent ancestors, constituting monophyletic clades, which suggests that each METTL was subject to strong selective constraints driving its structural and/or functional specialization. Interestingly, a similar process did not extend to the differentiation of nucleoside-modifying and protein-modifying METTLs (i.e., each METTL type did not form a unique monophyletic clade). The members of these two types of METTLs also exhibited differences in their rates of evolution. Overall, we provide evidence that the long-term evolution of METTL family members was driven by strong purifying selection, which in combination with adaptive selection episodes, led to the functional specialization of individual METTL lineages. This work contributes useful information regarding the evolution of a gene family that fulfills a variety of epigenetic functions, and can have profound influences on molecular processes and phenotypic traits.