Simulated annealing of microarray data reduces noise and enables cross-experimental comparisons

Microarrays are a powerful tool for assessing the genome-wide induction of a transcriptional response to internal or external stimuli, but are not considered quantitatively rigorous (i.e., the signal intensity of hybridized probe is normally used to quantify relative transcript abundance). Thus, it is difficult, if not impossible, to accurately compare separate microarray experiments without a reference standard. However, even among replicated microarray experiments, each gene varies significantly in the amount of signal detected, suggesting no single gene would be appropriate as a standard. We propose and test a method to "align" experimental transcription profiles to a set of reference experiments using simulated annealing (SA), essentially using the relative positions of all genes as a reference standard. SA attempts to find a globally optimal adjustment factor for the relative expression level of each experimental gene expression signal, given a previously observed range of gene expression measurements. By defining a relative dynamic range of gene expression under control conditions for all genes, we can more accurately compare transcription profiles between separate experiments and, potentially, between species--enabling comparative transcriptomics. Testing SA on a published dataset, we find that it significantly reduces interexperimental variation, suggesting it holds promise to accomplish this goal.