title: A flexible count data model to fit the wide diversity of expression profiles arising from extensively replicated RNA-seq experiments.
abstract: BACKGROUND
  High-throughput RNA sequencing (RNA-seq) offers unprecedented power to capture the real dynamics of gene expression. Experimental designs with extensive biological replication present a unique opportunity to exploit this feature and distinguish expression profiles with higher resolution. RNA-seq data analysis methods so far have been mostly applied to data sets with few replicates and their default settings try to provide the best performance under this constraint. These methods are based on two well-known count data distributions: the Poisson and the negative binomial. The way to properly calibrate them with large RNA-seq data sets is not trivial for the non-expert bioinformatics user.
  
  
  RESULTS
  Here we show that expression profiles produced by extensively-replicated RNA-seq experiments lead to a rich diversity of count data distributions beyond the Poisson and the negative binomial, such as Poisson-Inverse Gaussian or Pólya-Aeppli, which can be captured by a more general family of count data distributions called the Poisson-Tweedie. The flexibility of the Poisson-Tweedie family enables a direct fitting of emerging features of large expression profiles, such as heavy-tails or zero-inflation, without the need to alter a single configuration parameter. We provide a software package for R called tweeDEseq implementing a new test for differential expression based on the Poisson-Tweedie family. Using simulations on synthetic and real RNA-seq data we show that tweeDEseq yields P-values that are equally or more accurate than competing methods under different configuration parameters. By surveying the tiny fraction of sex-specific gene expression changes in human lymphoblastoid cell lines, we also show that tweeDEseq accurately detects differentially expressed genes in a real large RNA-seq data set with improved performance and reproducibility over the previously compared methodologies. Finally, we compared the results with those obtained from microarrays in order to check for reproducibility.
  
  
  CONCLUSIONS
  RNA-seq data with many replicates leads to a handful of count data distributions which can be accurately estimated with the statistical model illustrated in this paper. This method provides a better fit to the underlying biological variability; this may be critical when comparing groups of RNA-seq samples with markedly different count data distributions. The tweeDEseq package forms part of the Bioconductor project and it is available for download at http://www.bioconductor.org.
authors: Center for Research in Environmental Epidemiology (CREAL), Barcelona, Spain. robert.castelo@upf.edu.
authors: Mikel
authors: M
authors: Esnaola
authors: Pedro
authors: P
authors: Puig
authors: David
authors: D
authors: Gonzalez
authors: Robert
authors: R
authors: Castelo
authors: Juan R
authors: JR
authors: Gonzalez
mesh_terms: D019295:Computational Biology / Q000379:methods*; D020869:Gene Expression Profiling / Q000379:methods*; D059014:High-Throughput Nucleotide Sequencing; D006801:Humans; D015233:Models, Statistical; D015203:Reproducibility of Results; D017423:Sequence Analysis, RNA / Q000379:methods*; D012984:Software
references: Nature. 2010 Apr 1;464(7289):768-72
references: 20220758
references: Trends Genet. 2003 Jul;19(7):362-5
references: 12850439
references: BMC Bioinformatics. 2011 Dec 17;12:480
references: 22177264
references: Nat Biotechnol. 2011 Jul 11;29(7):572-3
references: 21747377
references: Biostatistics. 2008 Apr;9(2):321-32
references: 17728317
references: Biostatistics. 2013 Jan;14(1):113-28
references: 22988280
references: Biostatistics. 2013 Apr;14(2):232-43
references: 23001152
references: Proc Natl Acad Sci U S A. 2007 Jun 5;104(23):9758-63
references: 17537913
references: Stat Appl Genet Mol Biol. 2012 Oct 22;11(5):
references: 23104842
references: Biostatistics. 2012 Apr;13(2):204-16
references: 22285995
references: Stat Appl Genet Mol Biol. 2004;3:Article3
references: 16646809
references: Nucleic Acids Res. 2011 Jan;39(Database issue):D1011-5
references: 21177656
references: Nature. 2005 Mar 17;434(7031):400-4
references: 15772666
references: Genome Biol. 2006;7 Suppl 1:S4.1-9
references: 16925838
references: Nucleic Acids Res. 2012 May;40(10):4288-97
references: 22287627
references: Bioinformatics. 2013 Mar 1;29(5):614-21
references: 23329413
references: Nature. 2003 Jun 19;423(6942):825-37
references: 12815422
references: Genome Biol. 2010;11(3):R25
references: 20196867
references: Biometrics. 1997 Dec;53(4):1225-38
references: 9423246
references: Genome Res. 2008 Sep;18(9):1509-17
references: 18550803
references: Genome Biol. 2010;11(10):R106
references: 20979621
references: Genome Biol. 2004;5(10):R80
references: 15461798
references: Bioinformatics. 2010 Jan 1;26(1):139-40
references: 19910308
references: Nat Methods. 2012 Dec;9(12):1185-8
references: 23103880
references: Nat Genet. 2006 Jan;38(1):47-53
references: 16341221
references: Nat Methods. 2008 Jul;5(7):621-8
references: 18516045
references: BMC Bioinformatics. 2010 Aug 10;11:422
references: 20698981
references: Proc Natl Acad Sci U S A. 2003 Aug 5;100(16):9440-5
references: 12883005
references: PLoS Genet. 2007 Sep;3(9):1724-35
references: 17907809