Global rank-invariant set normalization (GRSN) to reduce systematic distortions in microarray data
1 Department of Molecular and Medical Genetics, Oregon Health and Sciences University, Portland, OR 97239-3098, USA
2 Department of Dermatology, Oregon Health and Sciences University, Portland, OR 97239-3098, USA
3 Department of Cell and Developmental Biology, Oregon Health and Sciences University, Portland, OR 97239-3098, USA
4 Department of Medicine, Oregon Health and Sciences University, Portland, OR 97239-3098, USA
5 OHSU Knight Cancer Institute, Oregon Health and Sciences University, Portland, OR 97239-3098, USA
BMC Bioinformatics 2008, 9:520 doi:10.1186/1471-2105-9-520Published: 4 December 2008
Microarray technology has become very popular for globally evaluating gene expression in biological samples. However, non-linear variation associated with the technology can make data interpretation unreliable. Therefore, methods to correct this kind of technical variation are critical. Here we consider a method to reduce this type of variation applied after three common procedures for processing microarray data: MAS 5.0, RMA, and dChip®.
We commonly observe intensity-dependent technical variation between samples in a single microarray experiment. This is most common when MAS 5.0 is used to process probe level data, but we also see this type of technical variation with RMA and dChip® processed data. Datasets with unbalanced numbers of up and down regulated genes seem to be particularly susceptible to this type of intensity-dependent technical variation. Unbalanced gene regulation is common when studying cancer samples or genetically manipulated animal models and preservation of this biologically relevant information, while removing technical variation has not been well addressed in the literature. We propose a method based on using rank-invariant, endogenous transcripts as reference points for normalization (GRSN). While the use of rank-invariant transcripts has been described previously, we have added to this concept by the creation of a global rank-invariant set of transcripts used to generate a robust average reference that is used to normalize all samples within a dataset. The global rank-invariant set is selected in an iterative manner so as to preserve unbalanced gene expression. Moreover, our method works well as an overlay that can be applied to data already processed with other probe set summary methods. We demonstrate that this additional normalization step at the "probe set level" effectively corrects a specific type of technical variation that often distorts samples in datasets.
We have developed a simple post-processing tool to help detect and correct non-linear technical variation in microarray data and demonstrate how it can reduce technical variation and improve the results of downstream statistical gene selection and pathway identification methods.