Variability of DNA microarray gene expression profiles in cultured rat primary hepatocytes.

ABSTRACT

DNA microarray is a powerful tool in biomedical research. However, transcriptomic profiling using DNA microarray is subject to many variations including biological variability. To evaluate the different sources of variation in mRNA gene expression profiles, gene expression profiles were monitored using the Affymetrix RatTox U34 arrays in cultured primary hepatocytes derived from six rats over a 26 hour period at 6 time points (0 h, 2h, 5h, 8h, 14 h and 26 h) with two replicate arrays at each time point for each animal. In addition, the impact of sample size on the variability of differentially expressed gene lists and the consistency of biological responses were also investigated. Excellent intra-animal reproducibility was obtained at all time points with 0 out of 370 present probe sets across all time points showing significant difference between the 2 replicate arrays (3-way ANOVA, p animal biological variation in mRNA expression profiles was observed with 337 out of 370 present probe sets showing significant differences among 6 animals (3-way ANOVA, p Principal Component Analysis (PCA) revealed that time effect (PC1) in this data set accounted for 47.4% of total variance indicating the dynamics of transcriptomics. The second and third largest effects came from animal difference, which accounted for 16.9% (PC2 and PC3) of the total variance. The reproducibility of gene lists and their functional classification was declined considerably when the sample size was decreased. Overall, our results strongly support that there is significant inter-animal variability in the time-course gene expression profiles, which is a confounding factor that must be carefully evaluated to correctly interpret microarray gene expression studies. The consistency of the gene lists and their biological functional classification are also sensitive to sample size with the reproducibility decreasing considerably under small sample size.