MOTIVATION: Pixel saturation occurs when the pixel intensity exceeds a threshold and the recorded pixel intensity is truncated. Microarray experiments are commonly afflicted with saturated pixels. As a result, estimators of gene expression are biased, with the amount of bias increasing as a function of the proportion of pixels saturated. Saturation is directly related to the photomultiplier tube (PMT) voltage settings and RNA abundance and is not necessarily associated with poor array or poor spot quality. When choosing PMT settings, higher PMT settings are desired because of improved signal-to-noise ratios of low-intensity spots. This improved signal is somewhat offset by saturation of high-intensity spots. In practice, spots with saturated pixels are discarded or the biased value is used. Neither of these approaches is appealing, particularly the former approach when a highly expressed gene is discarded because of saturation. RESULTS: We present a method to correct for saturation using pixel-level data. The method is based on a censored regression model. Evaluations on several arrays indicate that the method performs well. Simulation studies suggest that the method is robust under certain model violations.

Abstract Motivation: Recent technological advances such as cDNA microarray technology have made it possible to simultaneously interrogate thousands of genes in a biological specimen. A cDNA microarray experiment produces a gene expression ‘profile’. Often interest lies in discovering novel subgroupings, or ‘clusters’, of specimens based on their profiles, for example identification of new tumor taxonomies. Cluster analysis techniques such as hierarchical clustering and self-organizing maps have frequently been used for investigating structure in microarray data. However, clustering algorithms always detect clusters, even on random data, and it is easy to misinterpret the results without some objective measure of the reproducibility of the clusters. Results: We present statistical methods for testing for overall clustering of gene expression profiles, and we define easily interpretable measures of cluster-specific reproducibility that facilitate understanding of the clustering structure. We apply these methods to elucidate structure in cDNA microarray gene expression profiles obtained on melanoma tumors and on prostate specimens. Availability: Software to implement these methods is contained in BRB ArrayTools microarray analysis package available from http://linus.nci.nih.gov./BRB-ArrayTools.html Contact: [email protected] * To whom correspondence should be addressed. † Present address: Human Genome Sciences Inc., Rockville, MD 20850, USA