Technical and biological variance structure in mRNA-Seq data: life in the real world
1 Division of Biomedical Statistics and Informatics, Department of Health Sciences Research, Mayo Clinic, 200 1st St SW, Rochester, MN, 55905, USA
2 Mayo Vaccine Research Group, Mayo Clinic, 200 1st St SW, Rochester, MN, 55905, USA
3 Statistical Genetics, Sage Bionetworks, 1100 Fairview Ave N, M1-C108, Seattle, WA, 98109, USA
4 Program in Translational Immunovirology and Biodefense, Mayo Clinic, 200 1st St SW, Rochester, MN, 55905, USA
5 Department of Medicine, Mayo Clinic, 200 1st St SW, Rochester, MN, 55905, USA
BMC Genomics 2012, 13:304 doi:10.1186/1471-2164-13-304Published: 7 July 2012
mRNA expression data from next generation sequencing platforms is obtained in the form of counts per gene or exon. Counts have classically been assumed to follow a Poisson distribution in which the variance is equal to the mean. The Negative Binomial distribution which allows for over-dispersion, i.e., for the variance to be greater than the mean, is commonly used to model count data as well.
In mRNA-Seq data from 25 subjects, we found technical variation to generally follow a Poisson distribution as has been reported previously and biological variability was over-dispersed relative to the Poisson model. The mean-variance relationship across all genes was quadratic, in keeping with a Negative Binomial (NB) distribution. Over-dispersed Poisson and NB distributional assumptions demonstrated marked improvements in goodness-of-fit (GOF) over the standard Poisson model assumptions, but with evidence of over-fitting in some genes. Modeling of experimental effects improved GOF for high variance genes but increased the over-fitting problem.
These conclusions will guide development of analytical strategies for accurate modeling of variance structure in these data and sample size determination which in turn will aid in the identification of true biological signals that inform our understanding of biological systems.