scholarly journals zingeR: unlocking RNA-seq tools for zero-inflation and single cell applications

2017 ◽  
Author(s):  
Koen Van den Berge ◽  
Charlotte Soneson ◽  
Michael I. Love ◽  
Mark D. Robinson ◽  
Lieven Clement
Keyword(s):  
Single Cell ◽  
Differential Expression ◽  
Expression Analysis ◽  
Negative Binomial ◽  
Differential Expression Analysis ◽  
Negative Binomial Model ◽  
Binomial Model ◽  
Rna Seq ◽  
Zero Inflation ◽  
Zero Counts

AbstractDropout in single cell RNA-seq (scRNA-seq) applications causes many transcripts to go undetected. It induces excess zero counts, which leads to power issues in differential expression (DE) analysis and has triggered the development of bespoke scRNA-seq DE tools that cope with zero-inflation. Recent evaluations, however, have shown that dedicated scRNA-seq tools provide no advantage compared to traditional bulk RNA-seq tools. We introduce zingeR, a zero-inflated negative binomial model that identifies excess zero counts and generates observation weights to unlock bulk RNA-seq pipelines for zero-inflation, boosting performance in scRNA-seq differential expression analysis.

scholarly journals Observation weights to unlock bulk RNA-seq tools for zero inflation and single-cell applications

2018 ◽  
Author(s):  
Koen Van den Berge ◽  
Fanny Perraudeau ◽  
Charlotte Soneson ◽  
Michael I. Love ◽  
Davide Risso ◽  
...  
Keyword(s):  
Single Cell ◽  
Differential Expression ◽  
Transcriptome Sequencing ◽  
Negative Binomial ◽  
Rna Seq ◽  
Zero Inflation ◽  
Weighting Strategy ◽  
Cell Transcriptome ◽  
Single Cell Transcriptome ◽  
Zero Counts

AbstractDropout events in single-cell transcriptome sequencing (scRNA-seq) cause many transcripts to go undetected and induce an excess of zero read counts, leading to power issues in differential expression (DE) analysis. This has triggered the development of bespoke scRNA-seq DE methods to cope with zero inflation. Recent evaluations, however, have shown that dedicated scRNA-seq tools provide no advantage compared to traditional bulk RNA-seq tools. We introduce a weighting strategy, based on a zero-inflated negative binomial (ZINB) model, that identifies excess zero counts and generates gene and cell-specific weights to unlock bulk RNA-seq DE pipelines for zero-inflated data, boosting performance for scRNA-seq.

scholarly journals DEsingle for detecting three types of differential expression in single-cell RNA-seq data

2017 ◽  
Author(s):  
Zhun Miao ◽  
Ke Deng ◽  
Xiaowo Wang ◽  
Xuegong Zhang
Keyword(s):  
Single Cell ◽  
Differential Expression ◽  
Negative Binomial ◽  
Single Cells ◽  
R Package ◽  
Supplementary Information ◽  
Binomial Model ◽  
Supplementary Data ◽  
Rna Seq ◽  
Real Zeros

AbstractSummaryThe excessive amount of zeros in single-cell RNA-seq data include “real” zeros due to the on-off nature of gene transcription in single cells and “dropout” zeros due to technical reasons. Existing differential expression (DE) analysis methods cannot distinguish these two types of zeros. We developed an R package DEsingle which employed Zero-Inflated Negative Binomial model to estimate the proportion of real and dropout zeros and to define and detect 3 types of DE genes in single-cell RNA-seq data with higher accuracy.Availability and ImplementationThe R package DEsingle is freely available at https://github.com/miaozhun/DEsingle and is under Bioconductor’s consideration [email protected] informationSupplementary data are available at bioRxiv online.

scholarly journals Two-phase differential expression analysis for single cell RNA-seq

2018 ◽  
Vol 34 (19) ◽  
pp. 3340-3348 ◽  
Author(s):  
Zhijin Wu ◽  
Yi Zhang ◽  
Michael L Stitzel ◽  
Hao Wu
Keyword(s):  
Single Cell ◽  
Differential Expression ◽  
Expression Analysis ◽  
Differential Expression Analysis ◽  
Rna Seq ◽  
Two Phase

scholarly journals Bias, robustness and scalability in differential expression analysis of single-cell RNA-seq data

2017 ◽  
Author(s):  
Charlotte Soneson ◽  
Mark D. Robinson
Keyword(s):  
Single Cell ◽  
Differential Expression ◽  
Statistical Methods ◽  
Expression Analysis ◽  
Method Development ◽  
Differential Expression Analysis ◽  
Data Sets ◽  
Rna Seq ◽  
Data Set ◽  
Extensive Evaluation

AbstractBackgroundAs single-cell RNA-seq (scRNA-seq) is becoming increasingly common, the amount of publicly available data grows rapidly, generating a useful resource for computational method development and extension of published results. Although processed data matrices are typically made available in public repositories, the procedure to obtain these varies widely between data sets, which may complicate reuse and cross-data set comparison. Moreover, while many statistical methods for performing differential expression analysis of scRNA-seq data are becoming available, their relative merits and the performance compared to methods developed for bulk RNA-seq data are not sufficiently well understood.ResultsWe present conquer, a collection of consistently processed, analysis-ready public single-cell RNA-seq data sets. Each data set has count and transcripts per million (TPM) estimates for genes and transcripts, as well as quality control and exploratory analysis reports. We use a subset of the data sets available in conquer to perform an extensive evaluation of the performance and characteristics of statistical methods for differential gene expression analysis, evaluating a total of 30 statistical approaches on both experimental and simulated scRNA-seq data.ConclusionsConsiderable differences are found between the methods in terms of the number and characteristics of the genes that are called differentially expressed. Pre-filtering of lowly expressed genes can have important effects on the results, particularly for some of the methods originally developed for analysis of bulk RNA-seq data. Generally, however, methods developed for bulk RNA-seq analysis do not perform notably worse than those developed specifically for scRNA-seq.

scholarly journals Adjusted Sample Size Calculation for RNA-seq Data in the Presence of Confounding Covariates

2021 ◽  
Vol 1 (2) ◽  
pp. 47-63
Author(s):  
Xiaohong Li ◽  
Shesh N. Rai ◽  
Eric C. Rouchka ◽  
Timothy E. O’Toole ◽  
Nigel G. F. Cooper
Keyword(s):  
Sample Size ◽  
Differential Expression ◽  
Expression Analysis ◽  
Negative Binomial ◽  
Differential Expression Analysis ◽  
Sample Size Calculation ◽  
Size Estimation ◽  
Rna Seq ◽  
Simulation Based ◽  
The Relationship

Sample size calculation for adequate power analysis is critical in optimizing RNA-seq experimental design. However, the complexity increases for directly estimating sample size when taking into consideration confounding covariates. Although a number of approaches for sample size calculation have been proposed for RNA-seq data, most ignore any potential heterogeneity. In this study, we implemented a simulation-based and confounder-adjusted method to provide sample size recommendations for RNA-seq differential expression analysis. The data was generated using Monte Carlo simulation, given an underlined distribution of confounding covariates and parameters for a negative binomial distribution. The relationship between the sample size with the power and parameters, such as dispersion, fold change and mean read counts, can be visualized. We demonstrate that the adjusted sample size for a desired power and type one error rate of α is usually larger when taking confounding covariates into account. More importantly, our simulation study reveals that sample size may be underestimated by existing methods if a confounding covariate exists in RNA-seq data. Consequently, this underestimate could affect the detection power for the differential expression analysis. Therefore, we introduce confounding covariates for sample size estimation for heterogeneous RNA-seq data.

scholarly journals Valid post-clustering differential analysis for single-cell RNA-Seq

2018 ◽  
Author(s):  
Jesse M. Zhang ◽  
Govinda M. Kamath ◽  
David N. Tse
Keyword(s):  
Single Cell ◽  
Differential Expression ◽  
Expression Analysis ◽  
State Of The Art ◽  
Differential Expression Analysis ◽  
Differential Analysis ◽  
Rna Seq ◽  
Analysis Framework ◽  
Link Type ◽  
False Discoveries

SummarySingle-cell computational pipelines involve two critical steps: organizing cells (clustering) and identifying the markers driving this organization (differential expression analysis). State-of-the-art pipelines perform differential analysis after clustering on the same dataset. We observe that because clustering forces separation, reusing the same dataset generates artificially low p-values and hence false discoveries. We introduce a valid post-clustering differential analysis framework which corrects for this problem. We provide software at https://github.com/jessemzhang/tn_test.

scholarly journals Analysis of Single-Cell RNA-Sequencing Data: A Step-by-Step Guide

2021 ◽  
Vol 2 (1) ◽  
pp. 43-61
Author(s):  
Aanchal Malhotra ◽  
Samarendra Das ◽  
Shesh N. Rai
Keyword(s):  
Single Cell ◽  
Differential Expression ◽  
Rna Sequencing ◽  
Count Data ◽  
Negative Binomial ◽  
Expression Profiles ◽  
Differential Expression Analysis ◽  
Sequencing Data ◽  
Zero Inflation ◽  
Single Cell Rna Sequencing

Single-cell RNA-sequencing (scRNA-seq) technology provides an excellent platform for measuring the expression profiles of genes in heterogeneous cell populations. Multiple tools for the analysis of scRNA-seq data have been developed over the years. The tools require complicated commands and steps to analyze the underlying data, which are not easy to follow by genome researchers and experimental biologists. Therefore, we describe a step-by-step workflow for processing and analyzing the scRNA-seq unique molecular identifier (UMI) data from Human Lung Adenocarcinoma cell lines. We demonstrate the basic analyses including quality check, mapping and quantification of transcript abundance through suitable real data example to obtain UMI count data. Further, we performed basic statistical analyses, such as zero-inflation, differential expression and clustering analyses on the obtained count data. We studied the effects of excess zero-inflation present in scRNA-seq data on the downstream analyses. Our findings indicate that the zero-inflation associated with UMI data had no or minimal role in clustering, while it had significant effect on identifying differentially expressed genes. We also provide an insight into the comparative analysis for differential expression analysis tools based on zero-inflated negative binomial and negative binomial models on scRNA-seq data. The sensitivity analysis enhanced our findings in that the negative binomial model-based tool did not provide an accurate and efficient way to analyze the scRNA-seq data. This study provides a set of guidelines for the users to handle and analyze real scRNA-seq data more easily.

scholarly journals Differential expression analysis of RNA sequencing data by incorporating non-exonic mapped reads

2015 ◽  
Author(s):  
Hung-I Harry Chen ◽  
Yuanhang Liu ◽  
Yi Zou ◽  
Zhao Lai ◽  
Devanand Sarkar ◽  
...  
Keyword(s):  
Differential Expression ◽  
Rna Sequencing ◽  
Expression Analysis ◽  
Background Noise ◽  
Negative Binomial ◽  
Differential Expression Analysis ◽  
Rna Seq ◽  
Sequencing Data ◽  
Expression Levels ◽  
Differential Expressed Genes

Background RNA sequencing (RNA-seq) is a powerful tool for genome-wide expression profiling of biological samples with the advantage of high-throughput and high resolution. There are many existing algorithms nowadays for quantifying expression levels and detecting differential gene expression, but none of them takes the misaligned reads that are mapped to non-exonic regions into account. We developed a novel algorithm, XBSeq, where a statistical model was established based on the assumption that observed signals are the convolution of true expression signals and sequencing noises. The mapped reads in non-exonic regions are considered as sequencing noises, which follows a Poisson distribution. Given measureable observed and noise signals from RNA-seq data, true expression signals, assuming governed by the negative binomial distribution, can be delineated and thus the accurate detection of differential expressed genes. Results We implemented our novel XBSeq algorithm and evaluated it by using a set of simulated expression datasets under different conditions, using a combination of negative binomial and Poisson distributions with parameters derived from real RNA-seq data. We compared the performance of our method with other commonly used differential expression analysis algorithms. We also evaluated the changes in true and false positive rates with variations in biological replicates, differential fold changes, and expression levels in non-exonic regions. We also tested the algorithm on a set of real RNA-seq data where the common and different detection results from different algorithms were reported. Conclusions In this paper, we proposed a novel XBSeq, a differential expression analysis algorithm for RNA-seq data that takes non-exonic mapped reads into consideration. When background noise is at baseline level, the performance of XBSeq and DESeq are mostly equivalent. However, our method surpasses DESeq and other algorithms with the increase of non-exonic mapped reads. Only in very low read count condition XBSeq had a slightly higher false discovery rate, which may be improved by adjusting the background noise effect in this situation. Taken together, by considering non-exonic mapped reads, XBSeq can provide accurate expression measurement and thus detect differential expressed genes even in noisy conditions.

Sign in / Sign up

Export Citation Format

Share Document