A new efficient statistical test for detecting variability in the gene expression data

Abstract

DNA microarray technology allows researchers to monitor the expressions of thousands of genes under different conditions. The detection of differential gene expression under two different conditions is very important in microarray studies. Microarray experiments are multi-step procedures and each step is a potential source of variance. This makes the measurement of variability difficult because approach based on gene-by-gene estimation of variance will have few degrees of freedom. It is highly possible that the assumption of equal variance for all the expression levels may not hold. Also, the assumption of normality of gene expressions may not hold. Thus it is essential to have a statistical procedure which is not based on the normality assumption and also it can detect genes with differential variance efficiently. The detection of differential gene expression variance will allow us to identify experimental variables that affect different biological processes and accuracy of DNA microarray measurements.

In this article, a new nonparametric test for scale is developed based on the arctangent of the ratio of two expression levels. Most of the tests available in literature require the assumption of normal distribution, which makes them inapplicable in many situations, and it is also hard to verify the suitability of the normal distribution assumption for the given data set. The proposed test does not require the assumption of the distribution for the underlying population and hence makes it more practical and widely applicable. The asymptotic relative efficiency is calculated under different distributions, which show that the proposed test is very powerful when the assumption of normality breaks down. Monte Carlo simulation studies are performed to compare the power of the proposed test with some of the existing procedures. It is found that the proposed test is more powerful than commonly used tests under almost all the distributions considered in the study. A microarray data is used to illustrate the working of the proposed test. Results indicate that the proposed test is very powerful in detecting the smallest change in differential expression variance with high degree of confidence than some of its competitors.

Get full access to this article

View all access options for this article.

References

Novak JP , Sladek R. , Hudson TJ Characterization of variability in large-scale gene expression data: implications for study design. Genomics 2002 ; 79: 104-13.

Churchill GA Fundamentals of experimental design for cDNA microarray. Nature Genetics 2002; 32(Suppl): 490-95.

Yang YH , Speed TP Design issues for cDNA microarray experiments. Nature Review Genetics 2002; 3: 579-83.

Mood AM On the asymptotic efficiency of certain nonparametric two-sample tests . Annals of Mathematical Statistics 1954; 25: 514-33.

Freund JE , Ansari AR Two-way rank sum test for variances. Technical Report No. 34, Virginia Polytechnic Institute, Blacksburg, Virginia, 1957.

Klotz J. Nonparametric tests for scale. Annals of Mathematical Statistics 1962; 33: 498-512.

Siegel S. , Tukey JW A nonparametric sum of ranks procedure for relative spread in unpaired samples . Journal of American Statistical Association 1960; 55: 429-44.

Levene H. Robust tests for equality of variances. In Olkin I ed. Contributions to Probability and Statistics. Palo Alto, CA: Stanford University Press , 1960: 278-92.

Brown MB , Forsythe AB Robust tests for equality of variances. Journal of the American Statistical Association 1974; 69: 364-67.

10.

Conover WJ , Johnson ME , Johnson MM A comparative study of tests for homogeneity variances, with applications to the outer continental shelf bidding data. Technometrics , 1981; 23: 351-61.

11.

Lim TS , Loh W. A comparison of test of equality of variance. Computational Statistics and Data Analysis 1996; 22: 287-301.

12.

Kerr M. , Afshari C. , Bennet L. , Bushel P. , Martinej J. , Walker N. , Churchil G. Statistical analysis of a gene expression microarray experiment with replication . Statistica Sinica 2002; 12: 203-18.

13.

Dragichi S. , Kulaeva O. , Hoff B. , Petrov A. , Shams S. , Tainsky MA Noise sampling method: an ANOVA approach allowing robust selection of differentially regulated genes measured by DNA microarrays . Bioinformatics 2003; 19: 1348-59.

14.

Wolfinger RD , Gibson G. , Wolfinger ED , Bennet L. , Hamadeh H. , Bushel P. , Afshasi C. , Paules RS Assessing gene significance from cDNA microarray expression data via mixed models. Journal of Computational Biology 2001; 8: 625-37.

15.

Tusher V. , Tibshirani R. , Chu G. Significance analysis of microarrays applied to the ionizing radiation response . Proceeding of the National Academy of Sciences 2001; 98: 5116-24.

16.

Baldi P. , Long AD A Bayesian framework for the analysis of microarray expression data: regularized t-test and statistical inferences of gene changes. Bioinformatics 2001; 17: 509-19.

17.

Efron B. , Tibshirani R. , Storey JD , Tusher V. Empirical Bayes analysis of a microarray experiment. Journal of American Statistical Association 2001; 96: 1151-60.

18.

Lonnstedt I. , Speed T. Replicated microarray data. Statistica Sinica 2002; 12: 31-46.

19.

Delmar P. , Robin S. , Daudin JJ VarMixt: efficient variance modeling for the differential analysis of replicated gene expression data. Bioinformatics 2005; 21(4): 502-8.

20.

Chen JJ , Delongchamp RR , Tsai C. , Hsueh H. , Sistare F. , Thompson KL , Desai VG , Fuscoe JC Analysis of variance components in gene expression data. Bioinformatics 2004; 20(9): 1436-46.

21.

Wright GW , Simon RM A random variance model for detection of differential gene expression in small microarray experiments. Bioinformatics 2003 ; 19(18): 2448-55.

22.

Mathur SK A robust statistical method for identification of differentially expressed genes. Applied Bioinformatics 2005; 4(4): 247-52.

23.

Dudoit S. , Laan MJ Vander , Pollard KS Multiple testing. Part I. Single-Step procedures for control of general Type I error rates. Statistical Applications in Genetics and Molecular Biology 2004; 3: 1, Article 13.

24.

Arfin SM , Long AD , Ito ET , Tollesi L. , Richle MM , Paegle ES , Hatfield GW Global gene expression profiling in Escherichia coli K12 the effect of integration host factor. Journal of Biological Chemistry 2000; 275: 29672-864.