Home LiteratureArticle Details
PMID: 25516281 Published · ppublish English Journal Article Research Support, N.I.H., Extramural Research Support, Non-U.S. Gov't

Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2.

Genome biology ·Vol. 15 ·No. 12 ·2014-00-00 ·Pages 550

Love MI, Huber W, Anders S

Abstract

In comparative high-throughput sequencing assays, a fundamental task is the analysis of count data, such as read counts per gene in RNA-seq, for evidence of systematic changes across experimental conditions. Small replicate numbers, discreteness, large dynamic range and the presence of outliers require a suitable statistical approach. We present DESeq2, a method for differential analysis of count data, using shrinkage estimation for dispersions and fold changes to improve stability and interpretability of estimates. This enables a more quantitative analysis focused on the strength rather than the mere presence of differential expression. The DESeq2 package is available at http://www.bioconductor.org/packages/release/bioc/html/DESeq2.html webcite.

MeSH Terms
Algorithms Computational Biology/methods High-Throughput Nucleotide Sequencing Models, Genetic RNA/analysis Sequence Analysis, RNA Software
Chemicals
RNA
Authors & Affiliations
3 authors, click to expand affiliations / ORCID
Love Michael I
Huber Wolfgang
Anders Simon
References (47)
47 references, click to expand
  1. HTSeq--a Python framework to work with high-throughput sequencing data.
    Bioinformatics. 2015 Jan 15;31(2):166-9 PMID: 25260700
  2. Design and analysis of Bar-seq experiments.
    G3 (Bethesda). 2014 Jan;4(1):11-8 PMID: 24192834
  3. Identification and classification of differentially expressed genes in renal cell carcinoma by expression profiling on a global human 31,500-element cDNA array.
    Genome Res. 2001 Nov;11(11):1861-70 PMID: 11691851
  4. Therapeutic targeting of BET bromodomain proteins in castration-resistant prostate cancer.
    Nature. 2014 Jun 12;510(7504):278-82 PMID: 24759320
  5. A powerful and flexible approach to the analysis of RNA sequence count data.
    Bioinformatics. 2011 Oct 1;27(19):2672-8 PMID: 21810900
  6. ReCount: a multi-experiment resource of analysis-ready RNA-seq gene count datasets.
    BMC Bioinformatics. 2011;12:449 PMID: 22087737
  7. High-throughput screening of a CRISPR/Cas9 library for functional genomics in human cells.
    Nature. 2014 May 22;509(7501):487-91 PMID: 24717434
  8. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data.
    Bioinformatics. 2010 Jan 1;26(1):139-40 PMID: 19910308
  9. EBSeq: an empirical Bayes hierarchical model for inference in RNA-seq experiments.
    Bioinformatics. 2013 Apr 15;29(8):1035-43 PMID: 23428641
  10. Finding consistent patterns: a nonparametric approach for identifying differential expression in RNA-Seq data.
    Stat Methods Med Res. 2013 Oct;22(5):519-36 PMID: 22127579
  11. Regularization Paths for Generalized Linear Models via Coordinate Descent.
    J Stat Softw. 2010;33(1):1-22 PMID: 20808728
  12. easyRNASeq: a bioconductor package for processing RNA-Seq data.
    Bioinformatics. 2012 Oct 1;28(19):2532-3 PMID: 22847932
  13. Comparative ribosome profiling reveals extensive translational complexity in different Trypanosoma brucei life cycle stages.
    Nucleic Acids Res. 2014 Apr;42(6):3623-37 PMID: 24442674
  14. Removing technical variability in RNA-seq data using conditional quantile normalization.
    Biostatistics. 2012 Apr;13(2):204-16 PMID: 22285995
  15. Differential expression analysis for sequence count data.
    Genome Biol. 2010;11(10):R106 PMID: 20979621
  16. Robustly detecting differential expression in RNA sequencing data using observation weights.
    Nucleic Acids Res. 2014 Jun;42(11):e91 PMID: 24753412
  17. Complete alternative splicing events are bubbles in splicing graphs.
    J Comput Biol. 2009 Aug;16(8):1117-40 PMID: 19689216
  18. Differential analysis of gene regulation at transcript resolution with RNA-seq.
    Nat Biotechnol. 2013 Jan;31(1):46-53 PMID: 23222703
  19. Detecting differential usage of exons from RNA-seq data.
    Genome Res. 2012 Oct;22(10):2008-17 PMID: 22722343
  20. Evaluating gene expression in C57BL/6J and DBA/2J mouse striatum using RNA-Seq and microarrays.
    PLoS One. 2011;6(3):e17820 PMID: 21455293
  21. Bioconductor: open software development for computational biology and bioinformatics.
    Genome Biol. 2004;5(10):R80 PMID: 15461798
  22. A variance-stabilizing transformation for gene-expression microarray data.
    Bioinformatics. 2002;18 Suppl 1:S105-10 PMID: 12169537
  23. Testing significance relative to a fold-change threshold is a TREAT.
    Bioinformatics. 2009 Mar 15;25(6):765-71 PMID: 19176553
  24. voom: Precision weights unlock linear model analysis tools for RNA-seq read counts.
    Genome Biol. 2014;15(2):R29 PMID: 24485249
  25. Waste not, want not: why rarefying microbiome data is inadmissible.
    PLoS Comput Biol. 2014 Apr;10(4):e1003531 PMID: 24699258
  26. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions.
    Genome Biol. 2013;14(4):R36 PMID: 23618408
  27. GFOLD: a generalized fold change for ranking differentially expressed genes from RNA-seq data.
    Bioinformatics. 2012 Nov 1;28(21):2782-8 PMID: 22923299
  28. Linear models and empirical bayes methods for assessing differential expression in microarray experiments.
    Stat Appl Genet Mol Biol. 2004;3:Article3 PMID: 16646809
  29. NPEBseq: nonparametric empirical bayesian-based procedure for differential expression analysis of RNA-seq data.
    BMC Bioinformatics. 2013;14:262 PMID: 23981227
  30. Bayesian analysis of RNA sequencing data by estimating multiple shrinkage priors.
    Biostatistics. 2013 Jan;14(1):113-28 PMID: 22988280
  31. A new shrinkage estimator for dispersion improves differential expression detection in RNA-seq data.
    Biostatistics. 2013 Apr;14(2):232-43 PMID: 23001152
  32. Comparison of Affymetrix GeneChip expression measures.
    Bioinformatics. 2006 Apr 1;22(7):789-94 PMID: 16410320
  33. Identifying differentially expressed transcripts from RNA-seq data with biological variation.
    Bioinformatics. 2012 Jul 1;28(13):1721-8 PMID: 22563066
  34. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation.
    Nucleic Acids Res. 2012 May;40(10):4288-97 PMID: 22287627
  35. Small-sample estimation of negative binomial dispersion, with applications to SAGE data.
    Biostatistics. 2008 Apr;9(2):321-32 PMID: 17728317
  36. Differential oestrogen receptor binding is associated with clinical outcome in breast cancer.
    Nature. 2012 Jan 19;481(7381):389-93 PMID: 22217937
  37. Understanding mechanisms underlying human gene expression variation with RNA sequencing.
    Nature. 2010 Apr 1;464(7289):768-72 PMID: 20220758
  38. Moderated statistical tests for assessing differences in tag abundance.
    Bioinformatics. 2007 Nov 1;23(21):2881-7 PMID: 17881408
  39. Significance testing in ridge regression for genetic data.
    BMC Bioinformatics. 2011;12:372 PMID: 21929786
  40. baySeq: empirical Bayesian methods for identifying differential expression in sequence count data.
    BMC Bioinformatics. 2010;11:422 PMID: 20698981
  41. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features.
    Bioinformatics. 2014 Apr 1;30(7):923-30 PMID: 24227677
  42. GC-content normalization for RNA-Seq data.
    BMC Bioinformatics. 2011;12:480 PMID: 22177264
  43. mRNA-seq with agnostic splice site discovery for nervous system transcriptomics tested in chronic pain.
    Genome Res. 2010 Jun;20(6):847-60 PMID: 20452967
  44. Variance stabilization applied to microarray data calibration and to the quantification of differential expression.
    Bioinformatics. 2002;18 Suppl 1:S96-104 PMID: 12169536
  45. On differential variability of expression ratios: improving statistical inference about gene expression changes from microarray data.
    J Comput Biol. 2001;8(1):37-52 PMID: 11339905
  46. Software for computing and annotating genomic ranges.
    PLoS Comput Biol. 2013;9(8):e1003118 PMID: 23950696
  47. Independent filtering increases detection power for high-throughput experiments.
    Proc Natl Acad Sci U S A. 2010 May 25;107(21):9546-51 PMID: 20460310
Article Info
Journal
Genome biology
Abbr.
Genome Biol
ISSN
1474-760X
Published
2014-00-00
Pages
550
Language
English
Region
England
NLM ID
100960660
PMCID
PMC4302049
Subset
IM
Grants
NCI NIH HHS · T32 CA009337 · United States
NCI NIH HHS · 5T32CA009337-33 · United States
Analysis Services
Analysis Services

Contact

No. 2 Wenbo Road, Zhangqiu District, Jinan, Shandong

Qilu Normal University · Genelibs Bioinformatics Lab

750 Shunhua Rd, Jinan

2F, Bldg F, University Science Park

Tel: 0531-88819269

WeChat Official Account

Follow our WeChat subscription account for real-time updates and the latest in medical and biological research.


Business Email

E-mail: product@genelibs.com