Analyzing mRNA Epigenetic Sequencing Data with TRESS.

RNA epigenetics has emerged as an active topic to study gene regulation mechanisms. In this regard, the MeRIP-seq technology allows profiling transcriptome-wide mRNA modifications, in particular m6A. The primary goals for the analysis of MeRIP-seq data are the identification of m6A-methylated regions under each condition and across different biological conditions. Here we describe detailed procedures to guide researchers in MeRIP-seq data analyses by providing step-by-step instructions of the dedicated bioconductor package TRESS.

Recruitment of CTCF to an Fto enhancer is responsible for transgenerational inheritance of BPA-induced obesity.

The mechanisms by which environmentally-induced epiphenotypes are transmitted transgenerationally in mammals are poorly understood. Here we show that exposure of pregnant mouse females to bisphenol A (BPA) results in obesity in the F2 progeny due to increased food intake. This epiphenotype can be transmitted up to the F6 generation. Analysis of chromatin accessibility in sperm of the F1-F6 generations reveals alterations at sites containing binding motifs for CCCTC-binding factor (CTCF) at two cis-regulatory elements (CREs) of the Fto gene that correlate with transmission of obesity. These CREs show increased interactions in sperm of obese mice with the Irx3 and Irx5 genes, which are involved in the differentiation of appetite-controlling neurons. Deletion of the CTCF site in Fto results in mice that have normal food intake and fail to become obese when ancestrally exposed to BPA. The results suggest that epigenetic alterations of Fto can lead to the same phenotypes as genetic variants.

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

MOTIVATION: RNA epigenetics is an emerging field to study the post-transcriptional gene regulation. The dynamics of RNA epigenetic modification have been reported to associate with many human diseases. Recently developed high-throughput technology named Methylated RNA Immunoprecipitation Sequencing (MeRIP-seq) enables the transcriptome-wide profiling of N6-methyladenosine (m6A) modification and comparison of RNA epigenetic modifications. There are a few computational methods for the comparison of mRNA modifications under different conditions but they all suffer from serious limitations. RESULTS: In this work, we develop a novel statistical method to detect differentially methylated mRNA regions from MeRIP-seq data. We model the sequence count data by a hierarchical negative binomial model that accounts for various sources of variations and derive parameter estimation and statistical testing procedures for flexible statistical inferences under general experimental designs. Extensive benchmark evaluations in simulation and real data analyses demonstrate that our method is more accurate, robust and flexible compared to existing methods. AVAILABILITY AND IMPLEMENTATION: Our method TRESS is implemented as an R/Bioconductor package and is available at https://bioconductor.org/packages/devel/TRESS. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Epitranscriptomic dynamics in brain development and disease.

Distinct cell types are generated at specific times during brain development and are regulated by epigenetic, transcriptional, and newly emerging epitranscriptomic mechanisms. RNA modifications are known to affect many aspects of RNA metabolism and have been implicated in the regulation of various biological processes and in disease. Recent studies imply that dysregulation of the epitranscriptome may be significantly associated with neuropsychiatric, neurodevelopmental, and neurodegenerative disorders. Here we review the current knowledge surrounding the role of the RNA modifications N6-methyladenosine, 5-methylcytidine, pseudouridine, A-to-I RNA editing, 2'O-methylation, and their associated machinery, in brain development and human diseases. We also highlight the need for the development of new technologies in the pursuit of directly mapping RNA modifications in both genome- and single-molecule-level approach.

Dysregulated mitochondrial and cytosolic tRNA m1A methylation in Alzheimer's disease.

RNA modifications affect many aspects of RNA metabolism and are involved in the regulation of many different biological processes. Mono-methylation of adenosine in the N1 position, N1-methyladensoine (m1A), is a reversible modification that is known to target rRNAs and tRNAs. m1A has been shown to increase tRNA structural stability and induce correct tRNA folding. Recent studies have begun to associate the dysregulation of epitranscriptomic control with age-related disorders such as Alzheimer's disease. Here, we applied the newly developed m1A-quant-seq approach to map the brain abundant m1A RNA modification in the cortex of an Alzheimer's disease mouse model, 5XFAD. We observed hypomethylation in both mitochondrial and cytosolic tRNAs in 5XFAD mice compared with wild type. Furthermore, the main enzymes responsible for the addition of m1A in mitochondrial (TRMT10C, HSD17B10) and cytosolic tRNAs (TRMT61A) displayed decreased expression in 5XFAD compared with wild-type mice. Knockdown of these enzymes results in a more severe phenotype in a Drosophila tau model, and differential m1A methylation is correlated with differences in mature mitochondrial tRNA expression. Collectively, this work suggests that hypo m1A modification in tRNAs may play a role in Alzheimer's disease pathogenesis.

Detecting m6A methylation regions from Methylated RNA Immunoprecipitation Sequencing.

MOTIVATION: The post-transcriptional epigenetic modification on mRNA is an emerging field to study the gene regulatory mechanism and their association with diseases. Recently developed high-throughput sequencing technology named Methylated RNA Immunoprecipitation Sequencing (MeRIP-seq) enables one to profile mRNA epigenetic modification transcriptome wide. A few computational methods are available to identify transcriptome-wide mRNA modification, but they are either limited by over-simplified model ignoring the biological variance across replicates or suffer from low accuracy and efficiency. RESULTS: In this work, we develop a novel statistical method, based on an empirical Bayesian hierarchical model, to identify mRNA epigenetic modification regions from MeRIP-seq data. Our method accounts for various sources of variations in the data through rigorous modeling and applies shrinkage estimation by borrowing information from transcriptome-wide data to stabilize the parameter estimation. Simulation and real data analyses demonstrate that our method is more accurate, robust and efficient than the existing peak calling methods. AVAILABILITY AND IMPLEMENTATION: Our method TRES is implemented as an R package and is freely available on Github at https://github.com/ZhenxingGuo0015/TRES. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

N6-methyladenosine dynamics in neurodevelopment and aging, and its potential role in Alzheimer's disease.

BACKGROUND: N6-methyladenosine (m6A) modification is known to impact many aspects of RNA metabolism, including mRNA stability and translation, and is highly prevalent in the brain. RESULTS: We show that m6A modification displays temporal and spatial dynamics during neurodevelopment and aging. Genes that are temporally differentially methylated are more prone to have mRNA expression changes and affect many pathways associated with nervous system development. Furthermore, m6A shows a distinct tissue-specific methylation profile, which is most pronounced in the hypothalamus. Tissue-specific methylation is associated with an increase in mRNA expression and is associated with tissue-specific developmental processes. During the aging process, we observe significantly more m6A sites as age increases, in both mouse and human. We show a high level of overlap between mouse and human; however, humans at both young and old ages consistently show more m6A sites compared to mice. Differential m6A sites are found to be enriched in alternative untranslated regions of genes that affect aging-related pathways. These m6A sites are associated with a strong negative effect on mRNA expression. We also show that many Alzheimer-related transcripts exhibit decreased m6A methylation in a mouse model of Alzheimer's disease, which is correlated with reduced protein levels. CONCLUSIONS: Our results suggest that m6A exerts a critical function in both early and late brain development in a spatio-temporal fashion. Furthermore, m6A controls protein levels of key genes involved in Alzheimer's disease-associated pathways, suggesting that m6A plays an important role in aging and neurodegenerative disease.

Ago2-Dependent Processing Allows miR-451 to Evade the Global MicroRNA Turnover Elicited during Erythropoiesis.

MicroRNAs (miRNAs) are sequentially processed by two RNase III enzymes, Drosha and Dicer. miR-451 is the only known miRNA whose processing bypasses Dicer and instead relies on the slicer activity of Argonaute-2 (Ago2). miR-451 is highly conserved in vertebrates and regulates erythrocyte maturation, where it becomes the most abundant miRNA. However, the basis for the non-canonical biogenesis of miR-451 is unclear. Here, we show that Ago2 is less efficient than Dicer in processing pre-miRNAs, but this deficit is overcome when miR-144 represses Dicer in a negative-feedback loop during erythropoiesis. Loss of miR-144-mediated Dicer repression in zebrafish embryos and human cells leads to increased canonical miRNA production and impaired miR-451 maturation. Overexpression of Ago2 rescues some of the defects of miR-451 processing. Thus, the evolution of Ago2-dependent processing allows miR-451 to circumvent the global repression of canonical miRNAs elicited, in part, by the miR-144 targeting of Dicer during erythropoiesis.

Dynamic N6-methyladenosine RNA methylation in brain and diseases.

N6-methyladenosine (m6A) is a dynamic RNA modification that regulates various aspects of RNA metabolism and has been implicated in many biological processes and transitions. m6A is highly abundant in the brain; however, only recently has the role of m6A in brain development been a focus. The machinery that controls m6A is critically important for proper neurodevelopment, and the precise mechanisms by which m6A regulates these processes are starting to emerge. However, the role of m6A in neurodegenerative and neuropsychiatric diseases still requires much elucidation. This review discusses and summarizes the current body of knowledge surrounding the function of the m6A modification in regulating normal brain development, neurodegenerative diseases and outlines possible future directions.

Assessing miR-451 Activity and Its Role in Erythropoiesis.

The ability to microinject small RNAs and mRNAs into zebrafish embryos, of different genetic backgrounds, allows for the precise dissection of microRNA processing pathways at the molecular level, while simultaneously provides insight into their physiologic role. Here, we apply such an approach to determine the impact of Argonaute 2 in the processing of miR-451, a vertebrate-specific microRNA required for terminal erythrocyte differentiation. This was achieved using fluorescent microRNA reporter sensor assays and phenotype rescue experiments.