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.

[1]  Susan R. Wilson,et al.  An adsorption model of hybridization behaviour on oligonucleotide microarrays , 2004 .

[2]  Raimond L Winslow,et al.  Gene expression profiles in end-stage human idiopathic dilated cardiomyopathy: altered expression of apoptotic and cytoskeletal genes. , 2004, Genomics.

[3]  G. Grinstein,et al.  Modeling of DNA microarray data by using physical properties of hybridization , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[4]  Felix Naef,et al.  Absolute mRNA concentrations from sequence-specific calibration of oligonucleotide arrays. , 2003, Nucleic acids research.

[5]  Rafael A Irizarry,et al.  Exploration, normalization, and summaries of high density oligonucleotide array probe level data. , 2003, Biostatistics.

[6]  Raymond J Carroll,et al.  DNA Microarray Experiments: Biological and Technological Aspects , 2002, Biometrics.

[7]  Eric R. Ziegel,et al.  Generalized Linear Models , 2002, Technometrics.

[8]  R. Georgiadis,et al.  The effect of surface probe density on DNA hybridization. , 2001, Nucleic acids research.

[9]  D. Stern,et al.  Thermodynamics of Duplex Formation and Mismatch Discrimination on Photolithographically Synthesized Oligonucleotide Arrays , 1997 .

[10]  Y. Chen,et al.  Ratio-based decisions and the quantitative analysis of cDNA microarray images. , 1997, Journal of biomedical optics.

[11]  Ganapati P. Patil,et al.  The gamma distribution and weighted multimodal gamma distributions as models of population abundance , 1984 .

[12]  R. Sips,et al.  On the Structure of a Catalyst Surface , 1948 .

[13]  Christina Kendziorski,et al.  On Differential Variability of Expression Ratios: Improving Statistical Inference about Gene Expression Changes from Microarray Data , 2001, J. Comput. Biol..

[14]  E. Mammen The Bootstrap and Edgeworth Expansion , 1997 .

[15]  M. Muir Physical Chemistry , 1888, Nature.

[16]  Stat Pairs,et al.  Statistical Algorithms Description Document , 2022 .