title: Robustly detecting differential expression in RNA sequencing data using observation weights.
abstract: A popular approach for comparing gene expression levels between (replicated) conditions of RNA sequencing data relies on counting reads that map to features of interest. Within such count-based methods, many flexible and advanced statistical approaches now exist and offer the ability to adjust for covariates (e.g. batch effects). Often, these methods include some sort of 'sharing of information' across features to improve inferences in small samples. It is important to achieve an appropriate tradeoff between statistical power and protection against outliers. Here, we study the robustness of existing approaches for count-based differential expression analysis and propose a new strategy based on observation weights that can be used within existing frameworks. The results suggest that outliers can have a global effect on differential analyses. We demonstrate the effectiveness of our new approach with real data and simulated data that reflects properties of real datasets (e.g. dispersion-mean trend) and develop an extensible framework for comprehensive testing of current and future methods. In addition, we explore the origin of such outliers, in some cases highlighting additional biological or technical factors within the experiment. Further details can be downloaded from the project website: http://imlspenticton.uzh.ch/robinson_lab/edgeR_robust/.
authors.affiliation: Institute of Molecular Life Sciences, University of Zurich, CH-8057 Zurich, Switzerland SIB Swiss Institute of Bioinformatics, University of Zurich, CH-8057 Zurich, Switzerland.
authors.forename: Xiaobei
authors.initials: X
authors.lastname: Zhou
authors.affiliation: Institute of Molecular Life Sciences, University of Zurich, CH-8057 Zurich, Switzerland SIB Swiss Institute of Bioinformatics, University of Zurich, CH-8057 Zurich, Switzerland.
authors.forename: Helen
authors.initials: H
authors.lastname: Lindsay
authors.affiliation: Institute of Molecular Life Sciences, University of Zurich, CH-8057 Zurich, Switzerland SIB Swiss Institute of Bioinformatics, University of Zurich, CH-8057 Zurich, Switzerland mark.robinson@imls.uzh.ch.
authors.forename: Mark D
authors.initials: MD
authors.lastname: Robinson
mesh_terms: D003198:Computer Simulation; D020869:Gene Expression Profiling / Q000379:methods*; D017423:Sequence Analysis, RNA / Q000379:methods*
references.citation: Biostatistics. 2008 Apr;9(2):321-32
references.pmid: 17728317
references.citation: Biostatistics. 2013 Jan;14(1):113-28
references.pmid: 22988280
references.citation: BMC Bioinformatics. 2011 Nov 16;12:449
references.pmid: 22087737
references.citation: Bioinformatics. 2013 Apr 15;29(8):1035-43
references.pmid: 23428641
references.citation: BMC Bioinformatics. 2013 Mar 09;14:91
references.pmid: 23497356
references.citation: BMC Bioinformatics. 2010 Aug 10;11:422
references.pmid: 20698981
references.citation: Gastroenterology. 2005 Dec;129(6):2064-75
references.pmid: 16344072
references.citation: BMC Biol. 2010 May 11;8:58
references.pmid: 20459774
references.citation: Stat Appl Genet Mol Biol. 2013 Mar 26;12(1):49-70
references.pmid: 23502340
references.citation: Nucleic Acids Res. 2012 May;40(10):4288-97
references.pmid: 22287627
references.citation: Nat Biotechnol. 2011 May 15;29(7):644-52
references.pmid: 21572440
references.citation: Nat Protoc. 2013 Sep;8(9):1765-86
references.pmid: 23975260
references.citation: Nature. 2010 Apr 1;464(7289):773-7
references.pmid: 20220756
references.citation: Genome Biol. 2010;11(3):R25
references.pmid: 20196867
references.citation: Bioinformatics. 2007 Nov 1;23(21):2881-7
references.pmid: 17881408
references.citation: PLoS Biol. 2010 Sep 14;8(9):
references.pmid: 20856902
references.citation: BMC Bioinformatics. 2011 Aug 04;12:323
references.pmid: 21816040
references.citation: Genome Biol. 2010;11(10):R106
references.pmid: 20979621
references.citation: Nat Rev Genet. 2009 Jan;10(1):57-63
references.pmid: 19015660
references.citation: Genome Biol. 2004;5(10):R80
references.pmid: 15461798
references.citation: Bioinformatics. 2010 Jan 1;26(1):139-40
references.pmid: 19910308
references.citation: Biostatistics. 2013 Apr;14(2):232-43
references.pmid: 23001152
references.citation: Genome Res. 2010 Feb;20(2):180-9
references.pmid: 20009012
references.citation: Genome Biol. 2014 Feb 03;15(2):R29
references.pmid: 24485249
references.citation: Genome Biol. 2013;14(9):R95
references.pmid: 24020486
references.citation: Biostatistics. 2012 Apr;13(2):204-16
references.pmid: 22285995
references.citation: Front Genet. 2013 Sep 17;4:178
references.pmid: 24062766
references.citation: Am J Pathol. 1997 May;150(5):1647-59
references.pmid: 9137090
references.citation: Stat Methods Med Res. 2013 Oct;22(5):519-36
references.pmid: 22127579