Statistical Analysis of Adsorption Models for Oligonucleotide Microarrays

Recent analyses have shown that the relationship between intensity measurements from high density oligonucleotide microarrays and known concentration is non linear. Thus many measurements of so-called gene expression are neither measures of transcript nor mRNA concentration as might be expected.Intensity as measured in such microarrays is a measurement of fluorescent dye attached to probe-target duplexes formed during hybridization of a sample to the probes on the microarray. We develop several dynamic adsorption models relating fluorescent dye intensity to target RNA concentration, the simplest of which is the equilibrium Langmuir isotherm, or hyperbolic response function. Using data from the Affymerix HG-U95A Latin Square experiment, we evaluate various physical models, including equilibrium and non-equilibrium models, by applying maximum likelihood methods. We show that for these data, equilibrium Langmuir isotherms with probe dependent parameters are appropriate. We describe how probe sequence information may then be used to estimate the parameters of the Langmuir isotherm in order to provide an improved measure of absolute target concentration.