Statistical analysis of multiplex brain gene expression images

Analysis of variance (ANOVA) was employed to investigate 9,000 gene expression patterns from brains of both normal mice and mice with a pharmacological model of Parkinson's disease (PD). The data set was obtained using voxelation, a method that allows high-throughput acquisition of 3D gene expression patterns through analysis of spatially registered voxels (cubes). This method produces multiple volumetric maps of gene expression analogous to the images reconstructed in biomedical imaging systems. The ANOVA model was compared to the results from singular value decomposition (SVD) by using the first 42 singular vectors of the data matrix, a number equal to the rank of the ANOVA model. The ANOVA was also compared to the results from non-parametric statistics. Lastly, images were obtained for a subset of genes that emerged from the ANOVA as significant. The results suggest that ANOVA will be a valuable framework for insights into the large number of gene expression patterns obtained from voxelation.