A benchmark for RNA-seq deconvolution analysis under dynamic testing environments.
Ontology highlight
ABSTRACT: Background
Deconvolution analyses have been widely used to track compositional alterations of cell types in gene expression data. Although a large number of novel methods have been developed, due to a lack of understanding of the effects of modeling assumptions and tuning parameters, it is challenging for researchers to select an optimal deconvolution method suitable for the targeted biological conditions.Results
To systematically reveal the pitfalls and challenges of deconvolution analyses, we investigate the impact of several technical and biological factors including simulation model, quantification unit, component number, weight matrix, and unknown content by constructing three benchmarking frameworks. These frameworks cover comparative analysis of 11 popular deconvolution methods under 1766 conditions.Conclusions
We provide new insights to researchers for future application, standardization, and development of deconvolution tools on RNA-seq data.
SUBMITTER: Jin H
PROVIDER: S-EPMC8042713 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Obtaining RNA-seq measurements involves a complex data analytical process with a large number of competing algorithms as options. There is much debate about which of these methods provides the best approach. Unfortunately, it is currently difficult to evaluate their performance due in part to a lack of sensitive assessment metrics. We present a series of statistical summaries and plots to evaluate the performance in terms of specificity and sensitivity, available as a R/Bioconductor package ( http://bioconductor.org/packages/rnaseqcomp ). Using two independent datasets, we assessed seven competing pipelines. Performance was generally poor, with two methods clearly underperforming and RSEM slightly outperforming the rest.
| S-EPMC4842274 | biostudies-literature
Project description:By measuring messenger RNA levels for all genes in a sample, RNA-seq provides an attractive option to characterize the global changes in transcription. RNA-seq is becoming the widely used platform for gene expression profiling. However, real transcription signals in the RNA-seq data are confounded with measurement and sequencing errors and other random biological/technical variation. To extract biologically useful transcription process from the RNA-seq data, we propose to use the second ODE for modeling the RNA-seq data. We use differential principal analysis to develop statistical methods for estimation of location-varying coefficients of the ODE. We validate the accuracy of the ODE model to fit the RNA-seq data by prediction analysis and 5-fold cross validation. To further evaluate the performance of the ODE model for RNA-seq data analysis, we used the location-varying coefficients of the second ODE as features to classify the normal and tumor cells. We demonstrate that even using the ODE model for single gene we can achieve high classification accuracy. We also conduct response analysis to investigate how the transcription process responds to the perturbation of the external signals and identify dozens of genes that are related to cancer.
| S-EPMC4539434 | biostudies-literature
Project description:Thousands of RNA-binding proteins orchestrate RNA processing and altered protein-RNA interactions frequently lead to disease. Here, we present experimental and computational analysis pipelines of fractionated eCLIP-seq (freCLIP-seq), a modification of enhanced UV-crosslinking and RNA immunoprecipitation followed by sequencing. FreCLIP-seq allows transcriptome-wide analysis of protein-RNA interactions at single-nucleotide level and provides an additional level of resolution by isolating binding signals of individual RNA-binding proteins within a multicomponent complex. Binding occupancy can be inferred from read counts and crosslinking events. For complete details on the use and execution of this protocol, please refer to Biancon et al. (2022).
| S-EPMC9676202 | biostudies-literature
Project description:Hyperbaric oxygen (HBO) is widely applied to treat several hypoxia-related diseases. Previous studies have focused on the immediate effect of HBO-exposure induced oxidative stress on the lungs, but knowledge regarding the chronic effects from repetitive HBO exposure is limited, especially at the gene expression level. We found that repetitive HBO exposure did not alter the morphology of murine lungs. However, by deconvolution of RNA-seq from those mice lungs using CIBERSORTx and the expression profile matrices of 8 mesenchymal cell subtypes obtained from bleomycin-treated mouse lungs, we identify several mesenchymal cell subtype changes. These include increases in Col13a1 matrix fibroblasts, mesenchymal progenitors and mesothelial cell populations and decreases in lipofibroblasts, endothelial and Pdgfrb high cell populations. Our data suggest that repetitive HBO exposure may affect biological processes in the lungs such as response to wounding, extracellular matrix, vasculature development and immune response.
| S-EPMC7039706 | biostudies-literature
Project description:RNA-seq, a next-generation sequencing based method for transcriptome analysis, is rapidly emerging as the method of choice for comprehensive transcript abundance estimation. The accuracy of RNA-seq can be highly impacted by the purity of samples. A prominent, outstanding problem in RNA-seq is how to estimate transcript abundances in heterogeneous tissues, where a sample is composed of more than one cell type and the inhomogeneity can substantially confound the transcript abundance estimation of each individual cell type. Although experimental methods have been proposed to dissect multiple distinct cell types, computationally "deconvoluting" heterogeneous tissues provides an attractive alternative, since it keeps the tissue sample as well as the subsequent molecular content yield intact.Here we propose a probabilistic model-based approach, Transcript Estimation from Mixed Tissue samples (TEMT), to estimate the transcript abundances of each cell type of interest from RNA-seq data of heterogeneous tissue samples. TEMT incorporates positional and sequence-specific biases, and its online EM algorithm only requires a runtime proportional to the data size and a small constant memory. We test the proposed method on both simulation data and recently released ENCODE data, and show that TEMT significantly outperforms current state-of-the-art methods that do not take tissue heterogeneity into account. Currently, TEMT only resolves the tissue heterogeneity resulting from two cell types, but it can be extended to handle tissue heterogeneity resulting from multi cell types. TEMT is written in python, and is freely available at https://github.com/uci-cbcl/TEMT.The probabilistic model-based approach proposed here provides a new method for analyzing RNA-seq data from heterogeneous tissue samples. By applying the method to both simulation data and ENCODE data, we show that explicitly accounting for tissue heterogeneity can significantly improve the accuracy of transcript abundance estimation.
| S-EPMC3622628 | biostudies-literature
Project description:Because of the advantages of RNA sequencing (RNA-Seq) over microarrays, it is gaining widespread popularity for highly parallel gene expression analysis. For example, RNA-Seq is expected to be able to provide accurate identification and quantification of full-length splice forms. A number of informatics packages have been developed for this purpose, but short reads make it a difficult problem in principle. Sequencing error and polymorphisms add further complications. It has become necessary to perform studies to determine which algorithms perform best and which if any algorithms perform adequately. However, there is a dearth of independent and unbiased benchmarking studies. Here we take an approach using both simulated and experimental benchmark data to evaluate their accuracy.We conclude that most methods are inaccurate even using idealized data, and that no method is highly accurate once multiple splice forms, polymorphisms, intron signal, sequencing errors, alignment errors, annotation errors and other complicating factors are present. These results point to the pressing need for further algorithm development.Simulated datasets and other supporting information can be found at http://bioinf.itmat.upenn.edu/BEERS/bp2.Supplementary data are available at Bioinformatics online.
| S-EPMC4673975 | biostudies-literature
Project description:RNA-Seq analysis carried out in three different backcrosses based on Iberian breed with Landrace, Pietrain and Duroc breeds
2017-10-31 | PRJEB23068 | EVA
Project description:The differentiation and maturation trajectories of fetal liver stem/progenitor cells (LSPCs) are not fully understood at single-cell resolution, and a priori knowledge of limited biomarkers could restrict trajectory tracking.We employed marker-free single-cell RNA-Seq to characterize comprehensive transcriptional profiles of 507 cells randomly selected from seven stages between embryonic day 11.5 and postnatal day 2.5 during mouse liver development, and also 52 Epcam-positive cholangiocytes from postnatal day 3.25 mouse livers. LSPCs in developing mouse livers were identified via marker-free transcriptomic profiling. Single-cell resolution dynamic developmental trajectories of LSPCs exhibited contiguous but discrete genetic control through transcription factors and signaling pathways. The gene expression profiles of cholangiocytes were more close to that of embryonic day 11.5 rather than other later staged LSPCs, cuing the fate decision stage of LSPCs. Our marker-free approach also allows systematic assessment and prediction of isolation biomarkers for LSPCs.Our data provide not only a valuable resource but also novel insights into the fate decision and transcriptional control of self-renewal, differentiation and maturation of LSPCs.
| S-EPMC5715535 | biostudies-literature
Project description:Structural RNA domains are widely involved in the regulation of biological functions, such as gene expression, gene modification, and gene repair. Activity of these dynamic regions depends sensitively on the global fold of the RNA, in particular, on the binding affinity of individual conformations to effector molecules in solution. Consequently, both the 1) structure and 2) conformational dynamics of noncoding RNAs prove to be essential in understanding the coupling that results in biological function. Toward this end, we recently reported observation of three conformational states in the metal-induced folding pathway of the tRNA-like structure domain of Brome Mosaic Virus, via single-molecule fluorescence resonance energy transfer studies. We report herein selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE)-directed structure predictions as a function of metal ion concentrations ([Mn+]) to confirm the three-state folding model, as well as test 2° structure models from the literature. Specifically, SHAPE reactivity data mapped onto literature models agrees well with the secondary structures observed at 0-10 mM [Mg2+], with only minor discrepancies in the E hairpin domain at low [Mg2+]. SHAPE probing and SHAPE-directed structure predictions further confirm the stepwise unfolding pathway previously observed in our single-molecule studies. Of special relevance, this means that reduction in the metal-ion concentration unfolds the 3' pseudoknot interaction before unfolding the long-range stem interaction. This work highlights the synergistic power of combining 1) single-molecule Förster resonance energy transfer and 2) SHAPE-directed structure-probing studies for detailed analysis of multiple RNA conformational states. In particular, single-molecule guided deconvolution of the SHAPE reactivities permits 2° structure predictions of isolated RNA conformations, thereby substantially improving on traditional limitations associated with current structure prediction algorithms.
| S-EPMC5937115 | biostudies-literature
Project description:Studies on the crosstalk between muscle and adipose tissue can provide beneficial help in elucidating the pathogenesis and treatment of obesity-related diseases [1]. In this data article, we performed RNA sequence analysis of mRNA isolated from epididymal adipose tissue and gastrocnemius muscle tissue in obese rats. Twenty-two samples were selected for gene expression analysis. Raw data from the Illumina Hiseq™ platform sequencer was used for differential gene expression analysis using DESeq and deposited in the GEO public repository under accession number GSE237950. With the economic development and the change of people's lifestyle, obesity has become a major public health problem that endangers global health. Obesity is a metabolic disorder caused by excessive accumulation of white adipose tissue, which can further induce metabolic syndrome such as insulin resistance, type 2 diabetes, and cardiovascular and cerebrovascular diseases. Studies have shown that altitude hypoxic exercise can not only improve muscle buffering capacity and body performance, but also reduce body weight and body fat more significantly. In many countries, it has been used as a treatment program for obesity diseases [2]. Hypoxic exercise can improve lipid metabolism, reduce blood lipid levels, inhibit fatty acid synthesis, and promote fatty acid decomposition and oxidation, which is the mechanism of hypoxic exercise to significantly reduce weight and fat. However, the mechanism of the cross-talk between muscle and fat tissue is not well understood under hypoxia exercise and normoxia exercise conditions. The data contained rat's four different states: normoxia quiet, normoxia exercise, hypoxic quiet, and hypoxic exercise. RNA-seq data will provide insights into the cross-talk between muscle and fat, and the mechanisms of fat metabolism. The data of this study have not been published and are hereby published on this platform to study the cross talk between muscle tissue and adipose tissue in rats under different oxygen content and exercise environment.
| S-EPMC10859296 | biostudies-literature