Length bias correction for RNA-seq data in gene set analyses.
BACKGROUND
Next-generation sequencing technologies are being rapidly applied to quantifying transcripts (RNA-seq). However, due to the unique properties of the RNA-seq data, the differential expression of longer transcripts is more likely to be identified than that of shorter transcripts with the same effect size. This bias complicates the downstream gene set analysis (GSA) because the methods for GSA previously developed for microarray data are based on the assumption that genes with same effect size have equal probability (power) to be identified as significantly differentially expressed. Since transcript length is not related to gene expression, adjusting for such length dependency in GSA becomes necessary.


RESULTS
In this article, we proposed two approaches for transcript-length adjustment for analyses based on Poisson models: (i) At individual gene level, we adjusted each gene's test statistic using the square root of transcript length followed by testing for gene set using the Wilcoxon rank-sum test. (ii) At gene set level, we adjusted the null distribution for the Fisher's exact test by weighting the identification probability of each gene using the square root of its transcript length. We evaluated these two approaches using simulations and a real dataset, and showed that these methods can effectively reduce the transcript-length biases. The top-ranked GO terms obtained from the proposed adjustments show more overlaps with the microarray results.


BACKGROUND
R scripts are at http://www.soph.uab.edu/Statgenetics/People/XCui/r-codes/.
Department of Biostatistics, University of Alabama at Birmingham, Birmingham, AL 35294, USA.
Liyan
L
Gao
Zhide
Z
Fang
Kui
K
Zhang
Degui
D
Zhi
Xiangqin
X
Cui
D015982:Bias; D003198:Computer Simulation; D020869:Gene Expression Profiling / Q000379:methods*; D006801:Humans; D020411:Oligonucleotide Array Sequence Analysis; D012313:RNA / Q000032:analysis*; D017423:Sequence Analysis, RNA / Q000379:methods*
D012313:RNA
Biostatistics. 2008 Apr;9(2):321-32
17728317
Genome Biol. 2010;11(2):R14
20132535
Cell Mol Life Sci. 2010 Feb;67(4):569-79
19859660
Nat Rev Genet. 2006 Jan;7(1):55-65
16369572
Nat Rev Genet. 2009 Jan;10(1):57-63
19015660
Nat Genet. 2006 Jun;38(6):636-43
16715097
Nat Methods. 2008 Jul;5(7):621-8
18516045
Bioinformatics. 2003 Jan 22;19(2):185-93
12538238
Biol Direct. 2009 Apr 16;4:14
19371405
BMC Bioinformatics. 2010 Feb 18;11:94
20167110
Bioinformatics. 2009 Apr 15;25(8):1026-32
19244387
Bioinformatics. 2009 Mar 1;25(5):578-84
19168912
Biostatistics. 2003 Apr;4(2):249-64
12925520
Nat Methods. 2008 Jul;5(7):613-9
18516046
Bioinformatics. 2007 Nov 1;23(21):2881-7
17881408
Genome Res. 2008 Sep;18(9):1509-17
18550803
Genome Biol. 2010;11(10):R106
20979621
Nature. 2006 Jun 15;441(7095):840-6
16699522
Stat Appl Genet Mol Biol. 2004;3:Article3
16646809