Shrinkage estimation of dispersion in Negative Binomial models for RNA-seq experiments with small sample size

Motivation: RNA-seq experiments produce digital counts of reads that are affected by both biological and technical variation. To distinguish the systematic changes in expression between conditions from noise, the counts are frequently modeled by the Negative Binomial distribution. However, in experiments with small sample size, the per-gene estimates of the dispersion parameter are unreliable. Method: We propose a simple and effective approach for estimating the dispersions. First, we obtain the initial estimates for each gene using the method of moments. Second, the estimates are regularized, i.e. shrunk towards a common value that minimizes the average squared difference between the initial estimates and the shrinkage estimates. The approach does not require extra modeling assumptions, is easy to compute and is compatible with the exact test of differential expression. Results: We evaluated the proposed approach using 10 simulated and experimental datasets and compared its performance with that of currently popular packages edgeR, DESeq, baySeq, BBSeq and SAMseq. For these datasets, sSeq performed favorably for experiments with small sample size in sensitivity, specificity and computational time. Availability: http://www.stat.purdue.edu/∼ovitek/Software.html and Bioconductor. Contact: ovitek@purdue.edu Supplementary information: Supplementary data are available at Bioinformatics online.

• Publication year

  2013

• Venue

  Bioinform.

• Publication date

  2013-04-14

• Fields of study

  Biology, Mathematics, Computer Science, Medicine

• Identifiers
  DOI 10.1093/bioinformatics/btt143PMID 23589650PMCID 3654711
• External record

  Open on Semantic Scholar

• Source metadata

  Semantic Scholar, PubMed

• No claims are published for this paper.

• No concepts are published for this paper.

• A comparison of methods for differential expression analysis of RNA-seq data

  2013cited by this paper
• Finding consistent patterns: A nonparametric approach for identifying differential expression in RNA-Seq data

  2013influential reference
• Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation

  2012cited by this paper
• Normalization, testing, and false discovery rate estimation for RNA-sequencing data.

  2012cited by this paper
• Computational methods for transcriptome annotation and quantification using RNA-seq

  2011cited by this paper
• Differential expression analysis for sequence count data

  2011influential reference
• ReCount: A multi-experiment resource of analysis-ready RNA-seq gene count datasets

  2011cited by this paper
• Conservation of an RNA regulatory map between Drosophila and mammals.

  2011cited by this paper
• Statistical Applications in Genetics and Molecular Biology A Two-Stage Poisson Model for Testing RNA-Seq Data

  2011cited by this paper
• Statistical Applications in Genetics and Molecular Biology High-Dimensional Regression and Variable Selection Using CAR Scores

  2011cited by this paper
• Evaluating Gene Expression in C57BL/6J and DBA/2J Mouse Striatum Using RNA-Seq and Microarrays

  2011cited by this paper
• baySeq: Empirical Bayesian methods for identifying differential expression in sequence count data

  2010influential reference
• From RNA-seq reads to differential expression results

  2010influential reference
• Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments

  2010cited by this paper
• A scaling normalization method for differential expression analysis of RNA-seq data

  2010cited by this paper
• Evaluation of gene expression data generated from expired Affymetrix GeneChip® microarrays using MAQC reference RNA samples

  2010cited by this paper
• Sequencing technologies — the next generation

  2010cited by this paper
• Tumor Transcriptome Sequencing Reveals Allelic Expression Imbalances Associated with Copy Number Alterations

  2010cited by this paper
• Differential expression analysis for sequence count data

  2010cited by this paper
• mRNA-seq with agnostic splice site discovery for nervous system transcriptomics tested in chronic pain.

  2010cited by this paper
• Alternative expression analysis by RNA sequencing

  2010cited by this paper
• DEGseq: an R package for identifying differentially expressed genes from RNA-seq data

  2010cited by this paper
• Ultrafast and memory-efficient alignment of short DNA sequences to the human genome

  2009cited by this paper
• edgeR: a Bioconductor package for differential expression analysis of digital gene expression data

  2009cited by this paper
• Computation for ChIP-seq and RNA-seq studies

  2009cited by this paper
• How to Understand the Cell by Breaking It: Network Analysis of Gene Perturbation Screens

  2009influential reference
• RNA-Seq: a revolutionary tool for transcriptomics

  2009cited by this paper
• Cross-platform comparison of SYBR® Green real-time PCR with TaqMan PCR, microarrays and other gene expression measurement technologies evaluated in the MicroArray Quality Control (MAQC) study

  2008cited by this paper
• RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays.

  2008influential reference
• A Global View of Gene Activity and Alternative Splicing by Deep Sequencing of the Human Transcriptome

  2008cited by this paper
• Next-generation DNA sequencing methods.

  2008cited by this paper
• Maximum Likelihood Estimation of the Negative Binomial Dispersion Parameter for Highly Overdispersed Data, with Applications to Infectious Diseases

  2007cited by this paper
• Moderated statistical tests for assessing differences in tag abundance

  2007influential reference
• Performance comparison of one-color and two-color platforms within the Microarray Quality Control (MAQC) project

  2006cited by this paper
• The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements

  2006cited by this paper
• The Gamma-Poisson model as a statistical method to determine if micro-organisms are randomly distributed in a food matrix.

  2006cited by this paper
• Statistical practice in high-throughput screening data analysis

  2006influential reference
• INADMISSIBILITY OF THE USUAL ESTIMATOR FOR THE MEAN OF A MULTIVARIATE NORMAL DISTRIBUTION

  2005cited by this paper
• limma: Linear Models for Microarray Data

  2005cited by this paper
• A systematic comparison of methods for combining p

  2004influential reference
• Linear Models and Empirical Bayes Methods for Assessing Differential Expression in Microarray Experiments

  2004cited by this paper
• NIST/SEMATECH e-Handbook of Statistical Methods; Chapter 1: Exploratory Data Analysis

  2003cited by this paper
• Generalized Linear Models

  2002cited by this paper
• Regression Analysis of Count Data: Author Index

  1998cited by this paper
• Regression Analysis of Count Data

  1998influential reference
• Controlling the false discovery rate: a practical and powerful approach to multiple testing

  1995cited by this paper
• Generalized Linear Models

  1992cited by this paper
• Estimation with Quadratic Loss

  1992cited by this paper
• Maximum likelihood estimation for the negative binomial dispersion parameter.

  1990cited by this paper
• Estimation of the Negative Binomial Parameter κ by Maximum Quasi -Likelihood

  1989cited by this paper
• Extended moment series and the parameters of the negative binomial distribution

  1984cited by this paper
• Multistage estimation compared with fixed-sample-size estimation of the negative binomial parameter k

  1984cited by this paper
• Methods in medical research

  1967cited by this paper
• Proceedings of the Third Berkeley Symposium on Mathematical Statistics and Probability : held at the Statistical Laboratory, University of California, December, 1954, July and August, 1955

  1958cited by this paper
• Inadmissibility of the Usual Estimator for the Mean of a Multivariate Normal Distribution

  1956cited by this paper
• Theory of point estimation

  1950cited by this paper
• BIOINFORMATICS ORIGINAL PAPER

  year unknowncited by this paper

• EGFR-STAT1 pathway drives fibrosis initiation in fibroinflammatory skin diseases

  2025cites this paper
• Atherosclerosis enhances the efficacy of liposome-encapsulated bromocriptine in reducing the incidence of prolactinemia in pituitary tumors

  2025cites this paper
• Chromosomal instability shapes the tumor microenvironment of esophageal adenocarcinoma via a cGAS–chemokine–myeloid axis

  2025cites this paper
• Comprehensive single-cell RNA-sequencing study of Tollip deficiency effect in IL-13-stimulated human airway epithelial cells

  2025influential citation
• HAND1 level controls the specification of multipotent cardiac and extraembryonic progenitors from human pluripotent stem cells

  2025influential citation
• Enhancing Gradient-based Discrete Sampling via Parallel Tempering

  2025cites this paper
• Empirical partially Bayes multiple testing and compound χ2 decisions

  2025cites this paper
• Kölliker's Organ Functions as a Developmental Hub in Mouse Cochlea Regulating Spiral Limbus and Tectorial Membrane Development

  2025cites this paper
• Polymeric nano-formulation of spectrum selective RTK inhibitor strengthens anti-cancer effects via immune remodeling by endoplasmic reticulum stress-modulating mitochondrial metabolism

  2024cites this paper
• Comparative Methods for Demystifying Spatial Transcriptomics.

  2024cites this paper
• Spatial Transcriptomic Analysis Reveals HDAC Inhibition Modulates Microglial Dynamics to Protect Against Ischemic Stroke in Mice

  2024influential citation
• Identifying immune signatures of sepsis to increase diagnostic accuracy in very preterm babies

  2024cites this paper
• Changes revealed by histology, physiology, enteric microbiota, and transcriptome after ovarian ripening of shrimp, Litopenaeus vannamei, through using fresh baits support the innovation of artificial ripening of shrimps

  2024cites this paper
• Comparing Some Negative Binomial Regression with Simulation

  2024cites this paper
• The level of HAND1 controls the specification of multipotent cardiac and extraembryonic progenitors

  2024cites this paper
• An ex vivo model of Toxoplasma recrudescence reveals developmental plasticity of the bradyzoite stage

  2023cites this paper
• High resolution mapping of the tumor microenvironment using integrated single-cell, spatial and in situ analysis

  2023cites this paper
• Identification of new markers of angiogenic sprouting using transcriptomics: New role for RND3

  2023influential citation
• Single-Cell RNA Sequencing (scRNA-seq) Identifies L1CAM as a Key Mediator between Epithelial Tuft Cell and Innate Lymphoid Cell in the Colon of Hnrnp I Knockout Mice

  2023influential citation
• Dysregulation of murine long noncoding single-cell transcriptome in nonalcoholic steatohepatitis and liver fibrosis

  2023cites this paper
• Direct androgen receptor control of sexually dimorphic gene expression in the mammalian kidney.

  2023cites this paper
• IL-6 trans-signaling in a humanized mouse model of scleroderma

  2023cites this paper
• Spatial transcriptomics of a giant pilomatricoma

  2023cites this paper
• Single cell RNA-sequencing of human precision-cut lung slices: A novel approach to study the effect of vaping and viral infection on lung health

  2023cites this paper
• TCDD dysregulation of lncRNA expression, liver zonation and intercellular communication across the liver lobule.

  2023cites this paper
• Sexual dimorphism in obesity is governed by RELMα regulation of adipose macrophages and eosinophils

  2023cites this paper
• Direct androgen receptor regulation of sexually dimorphic gene expression in the mammalian kidney

  2023cites this paper
• TCDD dysregulation of lncRNA expression, liver zonation and intercellular communication across the liver lobule

  2023cites this paper
• Dysregulation of murine long non-coding single cell transcriptome in non-alcoholic steatohepatitis and liver fibrosis

  2022cites this paper
• Epiregulin is a dendritic cell-derived EGFR ligand that maintains skin and lung fibrosis

  2022cites this paper
• A missense variant in the nuclear localization signal of DKC1 causes Hoyeraal-Hreidarsson syndrome

  2022cites this paper
• AhR Activation Leads to Attenuation of Murine Autoimmune Hepatitis: Single-Cell RNA-Seq Analysis Reveals Unique Immune Cell Phenotypes and Gene Expression Changes in the Liver

  2022cites this paper
• Interplay between GH-regulated, sex-biased liver transcriptome and hepatic zonation revealed by single nucleus RNAseq.

  2022cites this paper
• Interplay between murine sex-biased gene expression and hepatic zonation revealed by single nucleus RNA sequencing

  2022cites this paper
• Single-cell transcriptomics unveils xylem cell development and evolution

  2021cites this paper
• Kidney pathology alters cell-clustering in single cell RNA sequencing

  2021cites this paper
• Reciprocal Transcription Factor Networks Govern Tissue-Resident ILC3 Subset Function and Identity

  2021cites this paper
• The Warburg effect is necessary to promote glycosylation in the blastema during zebrafish tail regeneration

  2021cites this paper
• A resampling strategy for studying robustness in virus detection pipelines

  2021cites this paper
• Title : RNA-seq analyses : Benchmarking differential expression analyses tools reveals the effect 2 of higher number of replicates on performance

  2021cites this paper
• Spatial transcriptomics analysis of uterine gene expression in enhancer of zeste homolog 2 conditional knockout mice

  2021cites this paper
• Effect of high variation in transcript expression on identifying differentially expressed genes in RNA‐seq analysis

  2021cites this paper
• Pre-capture multiplexing provides additional power to detect copy number variation in exome sequencing

  2021cites this paper
• Early satellite cell communication creates a permissive environment for long-term muscle growth

  2021cites this paper
• Single Cell Gene Expression Analysis in a 3D Microtissue Liver Model Reveals Cell Type-Specific Responses to Pro-Fibrotic TGF-β1 Stimulation

  2021cites this paper
• Where Do We Stand in Regularization for Life Science Studies?

  2021cites this paper
• Inference of 3D genome architecture by modeling overdispersion of Hi-C data

  2021cites this paper
• Identification of differentially distributed gene expression and distinct sets of cancer-related genes identified by changes in mean and variability

  2021cites this paper
• Bulk and Single-Cell Profiling of Breast Tumors Identifies TREM-1 as a Dominant Immune Suppressive Marker Associated With Poor Outcomes

  2021cites this paper
• Spatial transcriptomics reveals gene expression characteristics in invasive micropapillary carcinoma of the breast

  2021cites this paper
• Precisely modeling zero‐inflated count phenotype for rare variants

  2021cites this paper
• Single-Cell RNA Sequencing of Bone Marrow Mesenchymal Stem Cells from the Elderly People

  2021cites this paper
• Single-cell RNA Sequencing Reveals Heterogeneity of Cultured Bovine Satellite Cells

  2021cites this paper
• Conversion of human cardiac progenitor cells into cardiac pacemaker-like cells.

  2020cites this paper
• Choice of library size normalization and statistical methods for differential gene expression analysis in balanced two-group comparisons for RNA-seq studies

  2020cites this paper
• Transcription factor AP2 controls cnidarian germ cell induction

  2020cites this paper
• Single-Cell Transcriptomic Analysis.

  2020cites this paper
• A metabolic shift to glycolysis promotes zebrafish tail regeneration through TGF–β dependent dedifferentiation of notochord cells to form the blastema

  2020cites this paper
• A parsimonious family of multivariate Poisson-lognormal distributions for clustering multivariate count data

  2020cites this paper
• Chronic myelomonocytic leukaemia stem cell transcriptomes anticipate disease morphology and outcome

  2020cites this paper
• Single-Cell RNA-Sequencing From Mouse Incisor Reveals Dental Epithelial Cell-Type Specific Genes

  2020cites this paper
• A comprehensive overview of computational tools for RNA-seq analysis

  2020cites this paper
• Overexpression of PKMYT1 Facilitates Tumor Development and Is Correlated with Poor Prognosis in Clear Cell Renal Cell Carcinoma

  2020cites this paper
• A family of parsimonious mixtures of multivariate Poisson‐lognormal distributions for clustering multivariate count data

  2020cites this paper
• Model-Based Approach for Diffuse Glioma Classification, Grading, and Patient Survival Prediction

  2020cites this paper
• Developmental genetics of color pattern establishment in cats

  2020cites this paper
• A lineage-resolved molecular atlas of C. elegans embryogenesis at single-cell resolution

  2019cites this paper
• Identification of the lymphangioleiomyomatosis cell and its uterine origin

  2019influential citation
• R-spondin-3 induces secretory, antimicrobial Lgr5+ cells in the stomach

  2019cites this paper
• Stability of methods for differential expression analysis of RNA-seq data

  2019cites this paper
• Single-Cell RNA-Sequencing Identifies Activation of TP53 and STAT1 Pathways in Human T Lymphocyte Subpopulations in Response to Ex Vivo Radiation Exposure

  2019cites this paper
• Variance de l'expression des microARN et des ARN messagers dans le cancer

  2019cites this paper
• Mechanism of Enhanced MerTK-Dependent Macrophage Efferocytosis by Extracellular Vesicles.

  2019cites this paper
• Dissecting the heterogeneity of DENV vaccine-elicited cellular immunity using single-cell RNA sequencing and metabolic profiling

  2019cites this paper
• Analysis of high-throughput biological data using their rank values

  2019cites this paper
• A Review on Recent Statistical Models for RNA-Seq Data

  2019cites this paper
• Adaptive plasticity of IL-10+ and IL-35+ Treg cells cooperatively promotes tumor T cell exhaustion

  2019cites this paper
• Power and sample size calculations for high-throughput sequencing-based experiments

  2018cites this paper
• Finding the Association of mRNA and miRNA Using Next Generation Sequencing Data of Kidney Renal Cell Carcinoma

  2018cites this paper
• Single-cell transcriptomics reveal that PD-1 mediates immune tolerance by regulating proliferation of regulatory T cells

  2018cites this paper
• Deconstructing Adipogenesis Induced by β3-Adrenergic Receptor Activation with Single-Cell Expression Profiling.

  2018cites this paper
• Modeling Over-Dispersed Microbiome Data

  2018cites this paper
• Single-cell stabilization method identifies gonadotrope transcriptional dynamics and pituitary cell type heterogeneity

  2018cites this paper
• Non-contact competition between soft and hard corals: a transcriptomic perspective

  2018cites this paper
• Zero-Inflated Models for RNA-Seq Count Data

  2018cites this paper
• Regulatory Architecture of the LβT2 Gonadotrope Cell Underlying the Response to Gonadotropin-Releasing Hormone

  2018cites this paper
• Modeling and analysis of RNA-seq data: a review from a statistical perspective

  2018cites this paper
• Optimization of an RNA-Seq Differential Gene Expression Analysis Depending on Biological Replicate Number and Library Size

  2018cites this paper
• Single-Cell Deconvolution of Fibroblast Heterogeneity in Mouse Pulmonary Fibrosis

  2018cites this paper
• Classification of RNA-Seq Data via Gaussian Copulas

  2017cites this paper
• Non-equivalence of Wnt and R-spondin ligands during Lgr5+ intestinal stem cell self-renewal

  2017cites this paper
• Intestinal Enteroendocrine Lineage Cells Possess Homeostatic and Injury-Inducible Stem Cell Activity.

  2017cites this paper
• Normalization and microbial differential abundance strategies depend upon data characteristics

  2017cites this paper
• Evolution and characterization of carboplatin resistance at single cell resolution

  2017cites this paper
• RNA-seq 2G: online analysis of differential gene expression with comprehensive options of statistical methods

  2017cites this paper
• New Statistical Methods of Single-subject Transcriptome Analysis for Precision Medicine

  2017cites this paper
• CIRH1A augments the proliferation of RKO colorectal cancer cells.

  2017cites this paper
• Feasibility of sample size calculation for RNA‐seq studies

  2017cites this paper
• voomDDA: discovery of diagnostic biomarkers and classification of RNA-seq data

  2017cites this paper
• Proteobacteria explain significant functional variability in the human gut microbiome

  2017cites this paper