Enhancing Epitranscriptome Module Detection from m6A-Seq Data Using Threshold-Based Measurement Weighting Strategy

Recent advances in functional annotation and prediction of the epitranscriptome

RNA modifications, in particular N⁶-methyladenosine (m⁶A), participate in every stages of RNA metabolism and play diverse roles in essential biological processes and disease pathogenesis. Thanks to the advances in sequencing technology, tens of thousands of RNA modification sites can be identified in a typical high-throughput experiment; however, it remains a major challenge to decipher the functional relevance of these sites, such as, affecting alternative splicing, regulation circuit in essential biological processes or association to diseases. As the focus of RNA epigenetics gradually shifts from site discovery to functional studies, we review here recent progress in functional annotation and prediction of RNA modification sites from a bioinformatics perspective. The review covers naïve annotation with associated biological events, e.g., single nucleotide polymorphism (SNP), RNA binding protein (RBP) and alternative splicing, prediction of key sites and their regulatory functions, inference of disease association, and mining the diagnosis and prognosis value of RNA modification regulators. We further discussed the limitations of existing approaches and some future perspectives.

BDBB: A Novel Beta-distribution-based Biclustering Algorithm for Revealing Local Co-methylation Patterns in Epi-transcriptome Profiling Data

N6-methyladenosine (m⁶A) has been shown to play crucial roles in RNA metabolism, physiology, and pathological processes. However, the specific regulatory mechanisms of most methylation sites remain uncharted due to the complexity of life processes. Biological experimental methods are costly to solve this problem, and computational methods are relatively lacking. The discovery of local co-methylation patterns (LCPs) of m6A epi-transcriptome data can benefit to solve the above problems. Based on this, we propose a novel biclustering algorithm based on the beta distribution (BDBB), which realizes the mining of LCPs of m⁶A epi-transcriptome data. BDBB employs the Gibbs sampling method to complete parameter estimation. In the process of modeling, LCPs are recognized as sharp beta distributions compared to the background distribution. Simulation study showed BDBB can extract all the three actual LCPs implanted in the background data and the overlap conditions between them with considerable accuracy (almost close to 100%). On MeRIP-Seq data of 69,446 methylation sites under 32 experimental conditions from 10 human cell lines, BDBB unveiled two LCPs, and Gene Ontology (GO) enrichment analysis showed that they were enriched in histone modification and embryo development, etc. important biological processes respectively. The GOE_Score scoring indicated that the biclustering results of BDBB in the m⁶A epi-transcriptome data are more biologically meaningful than the results of other biclustering algorithms.

m6A Reader: Epitranscriptome Target Prediction and Functional Characterization of N 6-Methyladenosine (m6A) Readers

N ⁶-methyladenosine (m⁶A) is the most abundant post-transcriptional modification in mRNA, and regulates critical biological functions via m⁶A reader proteins that bind to m⁶A-containing transcripts. There exist multiple m⁶A reader proteins in the human genome, but their respective binding specificity and functional relevance under different biological contexts are not yet fully understood due to the limitation of experimental approaches. An in silico study was devised to unveil the target specificity and regulatory functions of different m⁶A readers. We established a support vector machine-based computational framework to predict the epitranscriptome-wide targets of six m⁶A reader proteins (YTHDF1-3, YTHDC1-2, and EIF3A) based on 58 genomic features as well as the conventional sequence-derived features. Our model achieved an average AUC of 0.981 and 0.893 under the full-transcript and mature mRNA model, respectively, marking a substantial improvement in accuracy compared to the sequence encoding schemes tested. Additionally, the distinct biological characteristics of each individual m⁶A reader were explored via the distribution, conservation, Gene Ontology enrichment, cellular components and molecular functions of their target m⁶A sites. A web server was constructed for predicting the putative binding readers of m⁶A sites to serve the research community, and is freely accessible at: http://m6areader.rnamd.com.

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

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

m6Acomet: large-scale functional prediction of individual m6A RNA methylation sites from an RNA co-methylation network

Background

Over one hundred different types of post-transcriptional RNA modifications have been identified in human. Researchers discovered that RNA modifications can regulate various biological processes, and RNA methylation, especially N6-methyladenosine, has become one of the most researched topics in epigenetics.

Results

To date, the study of epitranscriptome layer gene regulation is mostly focused on the function of mediator proteins of RNA methylation, i.e., the readers, writers and erasers. There is limited investigation of the functional relevance of individual m⁶A RNA methylation site. To address this, we annotated human m⁶A sites in large-scale based on the guilt-by-association principle from an RNA co-methylation network. It is constructed based on public human MeRIP-Seq datasets profiling the m⁶A epitranscriptome under 32 independent experimental conditions. By systematically examining the network characteristics obtained from the RNA methylation profiles, a total of 339,158 putative gene ontology functions associated with 1446 human m⁶A sites were identified. These are biological functions that may be regulated at epitranscriptome layer via reversible m⁶A RNA methylation. The results were further validated on a soft benchmark by comparing to a random predictor.

Conclusions

An online web server m6Acomet was constructed to support direct query for the predicted biological functions of m⁶A sites as well as the sites exhibiting co-methylated patterns at the epitranscriptome layer. The m6Acomet web server is freely available at: www.xjtlu.edu.cn/biologicalsciences/m6acomet.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2840-3) contains supplementary material, which is available to authorized users.

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.