RADAR: differential analysis of MeRIP-seq data with a random effect model.
Ontology highlight
ABSTRACT:
SUBMITTER: Zhang Z
PROVIDER: S-EPMC6927177 | biostudies-literature | 2019 Dec
REPOSITORIES: biostudies-literature
Publications
RADAR: differential analysis of MeRIP-seq data with a random effect model.
Zhang Zijie Z Zhan Qi Q Eckert Mark M Zhu Allen A Chryplewicz Agnieszka A De Jesus Dario F DF Ren Decheng D Kulkarni Rohit N RN Lengyel Ernst E He Chuan C Chen Mengjie M
Genome biology 20191223 1
Epitranscriptome profiling using MeRIP-seq is a powerful technique for in vivo functional studies of reversible RNA modifications. We develop RADAR, a comprehensive analytical tool for detecting differentially methylated loci in MeRIP-seq data. RADAR enables accurate identification of altered methylation sites by accommodating variability of pre-immunoprecipitation expression level and post-immunoprecipitation count using different strategies. In addition, it is compatible with complex study des ...[more]
Similar Datasets
Project description:Despite the prevalent studies of DNA/Chromatin related epigenetics, such as, histone modifications and DNA methylation, RNA epigenetics has not drawn deserved attention until a new affinity-based sequencing approach MeRIP-Seq was developed and applied to survey the global mRNA N6-methyladenosine (m(6)A) in mammalian cells. As a marriage of ChIP-Seq and RNA-Seq, MeRIP-Seq has the potential to study the transcriptome-wide distribution of various post-transcriptional RNA modifications. We have previously developed an R/Bioconductor package 'exomePeak' for detecting RNA methylation sites under a specific experimental condition or the identifying the differential RNA methylation sites in a case control study from MeRIP-Seq data. Compared with other relatively well studied data types such as ChIP-Seq and RNA-Seq, the study of MeRIP-Seq data is still at very early stage, and existing protocols are not optimized for dealing with the intrinsic characteristic of MeRIP-Seq data. We therein provide here a detailed and easy-to-use protocol of using exomePeak R/Bioconductor package along with other software programs for analysis of MeRIP-Seq data, which covers raw reads alignment, RNA methylation site detection, motif discovery, differential RNA methylation analysis, and functional analysis. Particularly, the rationales behind each processing step as well as the specific method used, the best practice, and possible alternative strategies are briefly discussed. The exomePeak R/Bioconductor package is freely available from Bioconductor: http://www.bioconductor.org/packages/release/bioc/html/exomePeak.html.
| S-EPMC4194139 | biostudies-literature
Project description:The recent advent of the state-of-art high throughput sequencing technology, known as Methylated RNA Immunoprecipitation combined with RNA sequencing (MeRIP-seq) revolutionizes the area of mRNA epigenetics and enables the biologists and biomedical researchers to have a global view of N (6)-Methyladenosine (m(6)A) on transcriptome. Yet there is a significant need for new computation tools for processing and analysing MeRIP-Seq data to gain a further insight into the function and m(6)A mRNA methylation.We developed a novel algorithm and an open source R package ( http://compgenomics.utsa.edu/metcluster ) for uncovering the potential types of m(6)A methylation by clustering the degree of m(6)A methylation peaks in MeRIP-Seq data. This algorithm utilizes a hierarchical graphical model to model the reads account variance and the underlying clusters of the methylation peaks. Rigorous statistical inference is performed to estimate the model parameter and detect the number of clusters. MeTCluster is evaluated on both simulated and real MeRIP-seq datasets and the results demonstrate its high accuracy in characterizing the clusters of methylation peaks. Our algorithm was applied to two different sets of real MeRIP-seq datasets and reveals a novel pattern that methylation peaks with less peak enrichment tend to clustered in the 5' end of both in both mRNAs and lncRNAs, whereas those with higher peak enrichment are more likely to be distributed in CDS and towards the 3'end of mRNAs and lncRNAs. This result might suggest that m(6)A's functions could be location specific.In this paper, a novel hierarchical graphical model based algorithm was developed for clustering the enrichment of methylation peaks in MeRIP-seq data. MeTCluster is written in R and is publicly available.
| S-EPMC5001242 | biostudies-literature
Project description:RNA modifications, especially methylation of the N(6) position of adenosine (A)-m(6)A, represent an emerging research frontier in RNA biology. With the rapid development of high-throughput sequencing technology, in-depth study of m(6)A distribution and function relevance becomes feasible. However, a robust method to effectively identify m(6)A-modified regions has not been available yet. Here, we present a novel high-efficiency and user-friendly analysis pipeline called MeRIP-PF for the signal identification of MeRIP-Seq data in reference to controls. MeRIP-PF provides a statistical P-value for each identified m(6)A region based on the difference of read distribution when compared to the controls and also calculates false discovery rate (FDR) as a cut off to differentiate reliable m(6)A regions from the background. Furthermore, MeRIP-PF also achieves gene annotation of m(6)A signals or peaks and produce outputs in both XLS and graphical format, which are useful for further study. MeRIP-PF is implemented in Perl and is freely available at http://software.big.ac.cn/MeRIP-PF.html.
| S-EPMC4357668 | biostudies-literature
Project description:RNA-Seq is increasingly being used for gene expression profiling. In this approach, next-generation sequencing (NGS) platforms are used for sequencing. Due to highly parallel nature, millions of reads are generated in a short time and at low cost. Therefore analysis of the data is a major challenge and development of statistical and computational methods is essential for drawing meaningful conclusions from this huge data. In here, we assessed three different types of normalization (transcript parts per million, trimmed mean of M values, quantile normalization) and evaluated if normalized data reduces technical variability across replicates. In addition, we also proposed two novel methods for detecting differentially expressed genes between two biological conditions: (i) likelihood ratio method, and (ii) Bayesian method. Our proposed methods for finding differentially expressed genes were tested on three real datasets. Our methods performed at least as well as, and often better than, the existing methods for analysis of differential expression.
| S-EPMC4393055 | biostudies-literature
Project description:BackgroundRNA-seq is a tool for measuring gene expression and is commonly used to identify differentially expressed genes (DEGs). Gene clustering is used to classify DEGs with similar expression patterns for the subsequent analyses of data from experiments such as time-courses or multi-group comparisons. However, gene clustering has rarely been used for analyzing simple two-group data or differential expression (DE). In this study, we report that a model-based clustering algorithm implemented in an R package, MBCluster.Seq, can also be used for DE analysis.ResultsThe input data originally used by MBCluster.Seq is DEGs, and the proposed method (called MBCdeg) uses all genes for the analysis. The method uses posterior probabilities of genes assigned to a cluster displaying non-DEG pattern for overall gene ranking. We compared the performance of MBCdeg with conventional R packages such as edgeR, DESeq2, and TCC that are specialized for DE analysis using simulated and real data. Our results showed that MBCdeg outperformed other methods when the proportion of DEG (PDEG) was less than 50%. However, the DEG identification using MBCdeg was less consistent than with conventional methods. We compared the effects of different normalization algorithms using MBCdeg, and performed an analysis using MBCdeg in combination with a robust normalization algorithm (called DEGES) that was not implemented in MBCluster.Seq. The new analysis method showed greater stability than using the original MBCdeg with the default normalization algorithm.ConclusionsMBCdeg with DEGES normalization can be used in the identification of DEGs when the PDEG is relatively low. As the method is based on gene clustering, the DE result includes information on which expression pattern the gene belongs to. The new method may be useful for the analysis of time-course and multi-group data, where the classification of expression patterns is often required.
| S-EPMC8527798 | biostudies-literature
Project description:RNA-Seq technology measures the transcript abundance by generating sequence reads and counting their frequencies across different biological conditions. To identify differentially expressed genes between two conditions, it is important to consider the experimental design as well as the distributional property of the data. In many RNA-Seq studies, the expression data are obtained as multiple pairs, e.g., pre- vs. post-treatment samples from the same individual. We seek to incorporate paired structure into analysis.We present a Bayesian hierarchical mixture model for RNA-Seq data to separately account for the variability within and between individuals from a paired data structure. The method assumes a Poisson distribution for the data mixed with a gamma distribution to account variability between pairs. The effect of differential expression is modeled by two-component mixture model. The performance of this approach is examined by simulated and real data.In this setting, our proposed model provides higher sensitivity than existing methods to detect differential expression. Application to real RNA-Seq data demonstrates the usefulness of this method for detecting expression alteration for genes with low average expression levels or shorter transcript length.
| S-EPMC3663822 | biostudies-literature
Project description:Uplands represent unique landscapes that provide a range of vital benefits to society, but are under increasing pressure from the management needs of a diverse number of stakeholders (e.g. farmers, conservationists, foresters, government agencies and recreational users). Mapping the spatial distribution of upland vegetation could benefit management and conservation programmes and allow for the impacts of environmental change (natural and anthropogenic) in these areas to be reliably estimated. The aim of this study was to evaluate the use of medium spatial resolution optical and radar satellite data, together with ancillary soil and topographic data, for identifying and mapping upland vegetation using the Random Forests (RF) algorithm. Intensive field survey data collected at three study sites in Ireland as part of the National Parks and Wildlife Service (NPWS) funded survey of upland habitats was used in the calibration and validation of different RF models. Eight different datasets were analysed for each site to compare the change in classification accuracy depending on the input variables. The overall accuracy values varied from 59.8% to 94.3% across the three study locations and the inclusion of ancillary datasets containing information on the soil and elevation further improved the classification accuracies (between 5 and 27%, depending on the input classification dataset). The classification results were consistent across the three different study areas, confirming the applicability of the approach under different environmental contexts.
| S-EPMC6686255 | biostudies-literature
Project description:RNA m5C methylation profile of MCF10A and MDA486 by using MeRIP-Seq protocol Immunoprecipitation of Methylated mRNA at Cytosine (m5C) residues: Affinity purified of anti-methyl cytosine (m5C) polyclonal antibody 7ug (Zymo Research, Catalog#A3001-50) was conjugated with protein-A magnetic beads for 2 h at 4°C in end to end rotator. After that, conjugated beads were extensively washed with RNA immunoprecipitation (RIP) wash buffer to remove unbound antibody. Fragmented 25 ug polyA RNA (mRNA) was incubated with m5C conjugated beads for overnight at 4°C in in the rotating platform in RIP buffer. RIP was done using Megna RNA Immunoprecipitation kit (Millipore, Catalog#17-700). m5C mRNA-immune bead complex was treated with proteinase K buffer to release m5C mRNA from the conjugated antibody. To isolate m5C, mRNA was treated with phenol:chloroform:isoamyl and mixed with 400 ul of chloroform, which was centrifuged at 14000 rpm for 10 minutes to separate aqueous phase. The aqueous phase was ethanol precipitated at -80°C for overnight, to get m5C mRNA. This precipitated m5C mRNA pellet was washed twice with 70% ethanol and air dried. Finally, m5C mRNA pellet was dissolved in nuclease free Water. The m5C mRNA integrity and conentration was quantified by bioanalyzer (Agilent) and Qubit 2.0 flurometer (Invitrogen). The fragmented mRNA was used by following TruSeq RNA Sample Preparation Guide to develop RNA-Seq library for sequencing.
| S-ECPF-GEOD-53370 | biostudies-other
Project description:We propose a functional random effect time-varying coefficient model to establish the dynamic relationship between the response and predictor variables in longitudinal data. This model allows us not only to interpret time-varying covariate effects, but also to depict random effects via time-varying profiles that are characterized by functional principal components. We develop the functional profiling-backfitting method to estimate model components, which includes the profiling and backfitting procedures via a set of least squares type estimating equations. Asymptotic properties of the resulting estimator are obtained. Furthermore, we investigate the finite sample performance of the proposed method through simulation studies and present an application to primary biliary cirrhosis data.
| S-EPMC3640349 | biostudies-literature
Project description:We propose eight data transformations (r, r2, rv, rv2, l, l2, lv, and lv2) for RNA-seq data analysis aiming to make the transformed sample mean to be representative of the distribution center since it is not always possible to transform count data to satisfy the normality assumption. Simulation studies showed that for data sets with small (e.g., nCases = nControls = 3) or large sample size (e.g., nCases = nControls = 100) limma based on data from the l, l2, and r2 transformations performed better than limma based on data from the voom transformation in term of accuracy, FDR, and FNR. For datasets with moderate sample size (e.g., nCases = nControls = 30 or 50), limma with the rv and rv2 transformations performed similarly to limma with the voom transformation. Real data analysis results are consistent with simulation analysis results: limma with the r, l, r2, and l2 transformation performed better than limma with the voom transformation when sample sizes are small or large; limma with the rv and rv2 transformations performed similarly to limma with the voom transformation when sample sizes are moderate. We also observed from our data analyses that for datasets with large sample size, the gene-selection via the Wilcoxon rank sum test (a non-parametric two sample test method) based on the raw data outperformed limma based on the transformed data.
| S-EPMC6423143 | biostudies-literature