PubMed Health⌕ Search

PubMed · 15127888

Statistical analysis of oligonucleotide microarray data.

Abstract

Microchip arrays have become one of the most rapidly growing techniques for monitoring gene expression at the genomic level and thereby gaining valuable insight about various important biological mechanisms. Examples of such mechanisms are: identifying disease-causing genes, genes involved in the regulation of some aspect of the cell cycle, etc. In this article, we discuss the problem of estimating gene expression based on a proper statistical model. More precisely, we show how the model introduced by Li and Wong can be used in its full bivariate generality to provide a new measure of gene expression from high-density oligonucleotide arrays. We also present a second gene expression index based on a new way of reducing the model into a simpler univariate model. In both cases, the gene expression indices are shown to be unbiased and to have lower variance than the established ones. Moreover, we present a bootstrap method aiming at providing non-parametric confidence intervals for the expression index.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Ziad Taib. 2004. Statistical analysis of oligonucleotide microarray data.. https://doi.org/10.1016/j.crvi.2003.05.003

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

Protocol to improve isoform-level quantification of low-abundance transcripts via STALARD pre-amplification.

STALARD (selective target amplification for low-abundance RNA detection) enables isoform-level quantification of low-abundance RNAs using conventional laboratory equipment. Here, we describe steps for RNA isolation, primer design, reverse transcription, selective target amplification, and downstream analysis. The protocol couples selective pre-amplification with a quantitative reverse-transcription PCR (RT-qPCR) readout and optional nanopore sequencing. Using 1 μg input RNA and 12 pre-amplification cycles, STALARD reduces Cq values by approximately 10-12 cycles, bringing the target into a reliably quantifiable range. For complete details on the use and execution of this protocol, please refer to Jeong et al.1.

Gene Expression↗

Gene expression profiles in microdissected neurons from human hippocampal subregions.

Pyramidal neurons in hippocampal subregions are selectively vulnerable in certain disease states. To investigate, we tested the hypothesis that selective vulnerability in human hippocampus is related to regional differences in neuronal cell death and cell receptor gene expression in CA1 vs. CA3 subregions. We used laser capture microdissection to remove approximately 600 CA1 and 600 CA3 pyramidal neurons each from five fresh-frozen normal post-mortem brains, extracted total RNA and double-amplified mRNA. This was reverse transcribed and labeled for hybridization onto human cDNA array chips containing probes to 10,174 genes and unknown ESTs. RNA from additional microdissections was pooled for replicate hybridizations and quantitative RT-PCR validation. Gene expression differences were few (< 1%). We found 43 enriched genes in CA1 neuronal samples that included peripheral benzodiazipine receptor-associated protein, nicotinic cholinergic receptor, two chemokine receptors (CCR1 and CCR5) and several transcriptional factors. We found 17 enriched genes in the CA3 neuronal samples that included fibroblast growth factor receptor and prostaglandin-endoperoxide synthase 1. We found no differential gene expression for 23 calcium channel proteins; nine transporter proteins; 55 cell death and apoptotic regulator proteins; and an additional 497 cell receptors, including 24 glutamate receptors. Quantitative RT-PCR of four differentially expressed genes confirmed the microarray data. The results confirm the ability to examine gene expression profiles in microdissected neurons from human autopsy brain. They show only minor gene expression differences between two distinct neuronal populations in the hippocampus and suggest that selective hippocampal vulnerability is due to factors other than intrinsic differential expression in glutamate receptors and cell death genes.

Gene Expression↗

A random-periods model for expression of cell-cycle genes.

We propose a nonlinear regression model for quantitatively analyzing periodic gene expression in studies of experimentally synchronized cells. Our model accounts for the observed attenuation in cycle amplitude by a simple and biologically plausible mechanism. We represent the expression level for each gene as an average across a large number of cells. For a given cell-cycle gene, we model its expression in each cell in the culture as following the same sinusoidal function except that the period, which in any individual cell must be the same for all cell-cycle genes, varies randomly across cells. We model these random periods by using a lognormal distribution. The variability in period causes the measured amplitude of the cyclic expression trajectory to attenuate over time as cells fall increasingly out of synchrony. Gene-specific parameters include initial amplitude and phase angle. Applying the model to data from Whitfield et al. [Whitfield, M. L., Sherlock, G., Saldanha, A. J., Murray, J, I., Ball, C. A., et al. (2002) Mol. Biol. Cell 13, 1977-2000], we fit the trajectories of 18 well characterized phase-marker genes and find that the fit does not suffer when a common lognormal distribution is assumed for all 18 genes compared with a separate distribution for each. We then use the model to identify 337 periodically expressed transcripts, including the 18 phase-marker genes. The model permits estimation of and hypothesis testing about biologically meaningful parameters that characterize cycling genes.

Gene Expression↗