Differential RNA methylation analysis for MeRIP-seq data under general experimental design

ALKBH5 regulates chicken adipogenesis by mediating LCAT mRNA stability depending on m6A modification

BACKGROUND

Previous studies have demonstrated the role of N6-methyladenosine (m6A) RNA methylation in various biological processes, our research is the first to elucidate its specific impact on LCAT mRNA stability and adipogenesis in poultry.

RESULTS

The 6 100-day-old female chickens were categorized into high (n = 3) and low-fat chickens (n = 3) based on their abdominal fat ratios, and their abdominal fat tissues were processed for MeRIP-seq and RNA-seq. An integrated analysis of MeRIP-seq and RNA-seq omics data revealed 16 differentially expressed genes associated with to differential m6A modifications. Among them, ELOVL fatty acid elongase 2 (ELOVL2), pyruvate dehydrogenase kinase 4 (PDK4), fatty acid binding protein 9 (PMP2), fatty acid binding protein 1 (FABP1), lysosomal associated membrane protein 3 (LAMP3), lecithin-cholesterol acyltransferase (LCAT) and solute carrier family 2 member 1 (SLC2A1) have ever been reported to be associated with adipogenesis. Interestingly, LCAT was down-regulated and expressed along with decreased levels of mRNA methylation methylation in the low-fat group. Mechanistically, the highly expressed ALKBH5 gene regulates LCAT RNA demethylation and affects LCAT mRNA stability. In addition, LCAT inhibits preadipocyte proliferation and promotes preadipocyte differentiation, and plays a key role in adipogenesis.

CONCLUSIONS

In conclusion, ALKBH5 mediates RNA stability of LCAT through demethylation and affects chicken adipogenesis. This study provides a theoretical basis for further understanding of RNA methylation regulation in chicken adipogenesis.

m6Aexpress-enet: Predicting the regulatory expression m6A sites by an enet-regularization negative binomial regression model

As the most abundant mRNA modification, m6A controls and influences many aspects of mRNA metabolism including the mRNA stability and degradation. However, the role of specific m6A sites in regulating gene expression still remains unclear. In additional, the multicollinearity problem caused by the correlation of methylation level of multiple m6A sites in each gene could influence the prediction performance. To address the above challenges, we propose an elastic-net regularized negative binomial regression model (called m6Aexpress-enet) to predict which m6A site could potentially regulate its gene expression. Comprehensive evaluations on simulated datasets demonstrate that m6Aexpress-enet could achieve the top prediction performance. Applying m6Aexpress-enet on real MeRIP-seq data from human lymphoblastoid cell lines, we have uncovered the complex regulatory pattern of predicted m6A sites and their unique enrichment pathway of the constructed co-methylation modules. m6Aexpress-enet proves itself as a powerful tool to enable biologists to discover the mechanism of m6A regulatory gene expression. Furthermore, the source code and the step-by-step implementation of m6Aexpress-enet is freely accessed at https://github.com/tengzhangs/m6Aexpress-enet.

magpie: A power evaluation method for differential RNA methylation analysis in N6-methyladenosine sequencing

Recently, novel biotechnologies to quantify RNA modifications became an increasingly popular choice for researchers who study epitranscriptome. When studying RNA methylations such as N6-methyladenosine (m6A), researchers need to make several decisions in its experimental design, especially the sample size and a proper statistical power. Due to the complexity and high-throughput nature of m6A sequencing measurements, methods for power calculation and study design are still currently unavailable. In this work, we propose a statistical power assessment tool, magpie, for power calculation and experimental design for epitranscriptome studies using m6A sequencing data. Our simulation-based power assessment tool will borrow information from real pilot data, and inspect various influential factors including sample size, sequencing depth, effect size, and basal expression ranges. We integrate two modules in magpie: (i) a flexible and realistic simulator module to synthesize m6A sequencing data based on real data; and (ii) a power assessment module to examine a set of comprehensive evaluation metrics.

m6A-Atlas v2.0: updated resources for unraveling the N6-methyladenosine (m6A) epitranscriptome among multiple species

N 6-Methyladenosine (m6A) is one of the most abundant internal chemical modifications on eukaryote mRNA and is involved in numerous essential molecular functions and biological processes. To facilitate the study of this important post-transcriptional modification, we present here m6A-Atlas v2.0, an updated version of m6A-Atlas. It was expanded to include a total of 797 091 reliable m6A sites from 13 high-resolution technologies and two single-cell m6A profiles. Additionally, three methods (exomePeaks2, MACS2 and TRESS) were used to identify >16 million m6A enrichment peaks from 2712 MeRIP-seq experiments covering 651 conditions in 42 species. Quality control results of MeRIP-seq samples were also provided to help users to select reliable peaks. We also estimated the condition-specific quantitative m6A profiles (i.e. differential methylation) under 172 experimental conditions for 19 species. Further, to provide insights into potential functional circuitry, the m6A epitranscriptomics were annotated with various genomic features, interactions with RNA-binding proteins and microRNA, potentially linked splicing events and single nucleotide polymorphisms. The collected m6A sites and their functional annotations can be freely queried and downloaded via a user-friendly graphical interface at: http://rnamd.org/m6a.