Biases induced by pooling samples in microarray experiments

MOTIVATION: If there is insufficient RNA from the tissues under investigation from one organism, then it is common practice to pool RNA. An important question is to determine whether pooling introduces biases, which can lead to inaccurate results. In this article, we describe two biases related to pooling, from a theoretical as well as a practical point of view. RESULTS: We model and quantify the respective parts of the pooling bias due to the log transform as well as the bias due to biological averaging of the samples. We also evaluate the impact of the bias on the statistical differential analysis of Affymetrix data.     

poolMC: Smart pooling of mRNA samples in microarray experiments

Background  

Typically, pooling of mRNA samples in microarray experiments implies mixing mRNA from several biological-replicate samples before hybridization onto a microarray chip. Here we describe an alternative smart pooling strategy in which different samples, not necessarily biological replicates, are pooled in an information theoretic efficient way. Further, each sample is tested on multiple chips, but always in pools made up of different samples. The end goal is to exploit the compressibility of microarray data to reduce the number of chips used and increase the robustness to noise in measurements.     

Questioning the utility of pooling samples in microarray experiments with cell lines

We describe a microarray experiment using the MCF-7 breast cancer cell line in two different experimental conditions for which the same number of independent pools as the number of individual samples was hybridized on Affymetrix GeneChips. Unexpectedly, when using individual samples, the number of probe sets found to be differentially expressed between treated and untreated cells was about three times greater than that found using pools. These findings indicate that pooling samples in microarray experiments where the biological variability is expected to be small might not be helpful and could even decrease one's ability to identify differentially expressed genes.     

The efficiency of pooling mRNA in microarray experiments

In a microarray experiment, messenger RNA samples are oftentimes pooled across subjects out of necessity, or in an effort to reduce the effect of biological variation. A basic problem in such experiments is to estimate the nominal expression levels of a large number of genes. Pooling samples will affect expression estimation, but the exact effects are not yet known as the approach has not been systematically studied in this context. We consider how mRNA pooling affects expression estimates by assessing the finite-sample performance of different estimators for designs with and without pooling. Conditions under which it is advantageous to pool mRNA are defined; and general properties of estimates from both pooled and non-pooled designs are derived under these conditions. A formula is given for the total number of subjects and arrays required in a pooled experiment to obtain gene expression estimates and confidence intervals comparable to those obtained from the no-pooling case. The formula demonstrates that by pooling a perhaps increased number of subjects, one can decrease the number of arrays required in an experiment without a loss of precision. The assumptions that facilitate derivation of this formula are considered using data from a quantitative real-time PCR experiment. The calculations are not specific to one particular method of quantifying gene expression as they assume only that a single, normalized, estimate of expression is obtained for each gene. As such, the results should be generally applicable to a number of technologies provided sufficient pre-processing and normalization methods are available and applied.     

Assessment of repeated microarray experiments using mixed tissue RNA reference samples

Genome-scale gene expression technologies are increasingly being applied for biological research as a whole and toxicological screening in particular. In order to monitor data quality and process drift, we adopted the use of two rat-tissue mixtures (brain, liver, kidney, and testis) previously introduced as RNA reference samples. These samples were processed every time a microarray experiment was hybridized, thereby verifying the comparability of the resulting expression data for cross-study comparison. This study presents the analysis of 21 technical replicates of these two mixed-tissue samples using Affymetrix RAE230_2 GeneChip over a period of 12 months. The results show that detection sensitivity, measured by the number of present and absent sequences, is robust, and data correlation, indicated by scatter plots, varies little over time. Receiver operating characteristic (ROC) curves show the sensitivity and specificity of the current measurements are consistent with arrays previously classified as well performing. Overall, this paper shows that the inclusion of standard samples during microarray labeling and hybridization experiments is useful to benchmark the performance of microarray experiments over time and allows discovery of any process drift that, if it occurs, may confound the comparison of these datasets.     

Statistical designs for two-color spotted microarray experiments

Two-color cDNA or oligonucleotide-based spotted microarrays have been commonly used in measuring the expression levels of thousands of genes simultaneously. To realize the immense potential of this powerful new technology, budgeted within limited resources or other constraints, practical designs with high efficiencies are in demand. In this study, we address the design issue concerning the arrangement of the mRNA samples labeled with fluorescent dyes and hybridized on the slides. A normalization model is proposed to characterize major sources of systematic variation in a two-color microarray experiment. This normalization model establishes a connection between designs for two-color microarray experiments with a particular class of classical row-column designs. A heuristic algorithm for constructing A-optimal or highly efficient designs is provided. Statistical optimality results are found for some of the designs generated from the algorithm. It is believed that the constructed designs are the best or very close to the best possible for estimating the relative gene expression levels among the mRNA samples of interest.     

Statistical process control for large scale microarray experiments

MOTIVATION: Maintaining and controlling data quality is a key problem in large scale microarray studies. In particular systematic changes in experimental conditions across multiple chips can seriously affect quality and even lead to false biological conclusions. Traditionally the influence of these effects can be minimized only by expensive repeated measurements, because a detailed understanding of all process relevant parameters seems impossible. RESULTS: We introduce a novel method for microarray process control that estimates quality based solely on the distribution of the actual measurements without requiring repeated experiments. A robust version of principle component analysis detects single outlier microarrays and thereby enables the use of techniques from multivariate statistical process control. In particular, the T(2) control chart reliably tracks undesired changes in process relevant parameters. This can be used to improve the microarray process itself, limits necessary repetitions to only affected samples and therefore maintains quality in a cost effective way. We prove the power of the approach on 3 large sets of DNA methylation microarray data.     

Statistical methodology for the analysis of dye-switch microarray experiments

Background  

In individually dye-balanced microarray designs, each biological sample is hybridized on two different slides, once with Cy3 and once with Cy5. While this strategy ensures an automatic correction of the gene-specific labelling bias, it also induces dependencies between log-ratio measurements that must be taken into account in the statistical analysis.     

Statistical tests for differential expression in cDNA microarray experiments

Extracting biological information from microarray data requires appropriate statistical methods. The simplest statistical method for detecting differential expression is the t test, which can be used to compare two conditions when there is replication of samples. With more than two conditions, analysis of variance (ANOVA) can be used, and the mixed ANOVA model is a general and powerful approach for microarray experiments with multiple factors and/or several sources of variation.     

Statistical tests for identifying differentially expressed genes in time-course microarray experiments

MOTIVATION: Microarray technology allows the monitoring of expression levels for thousands of genes simultaneously. In time-course experiments in which gene expression is monitored over time, we are interested in testing gene expression profiles for different experimental groups. However, no sophisticated analytic methods have yet been proposed to handle time-course experiment data. RESULTS: We propose a statistical test procedure based on the ANOVA model to identify genes that have different gene expression profiles among experimental groups in time-course experiments. Especially, we propose a permutation test which does not require the normality assumption. For this test, we use residuals from the ANOVA model only with time-effects. Using this test, we detect genes that have different gene expression profiles among experimental groups. The proposed model is illustrated using cDNA microarrays of 3840 genes obtained in an experiment to search for changes in gene expression profiles during neuronal differentiation of cortical stem cells.     

Systematic analysis of T7 RNA polymerase basedin vitrolinear RNA amplification for use in microarray experiments

Background  

The requirement of a large amount of high-quality RNA is a major limiting factor for microarray experiments using biopsies. An average microarray experiment requires 10–100 μg of RNA. However, due to their small size, most biopsies do not yield this amount. Several different approaches for RNA amplificationin vitrohave been described and applied for microarray studies. In most of these, systematic analyses of the potential bias introduced by the enzymatic modifications are lacking.     

Impact of RNA extraction from limited samples on microarray results

To move microarray technology into the diagnostic realm, the impact of technical parameters, such as sample preparation and RNA extraction, needs to be understood and minimized. We evaluated the impact of two RNA extraction methods, DNase treatment and the amount of hybridized cDNA probe, on the outcome of microarray results. The results for both RNA extraction methods were comparable, although one method resulted in residual DNA that slightly affected the microarray results. As little as one microgram of total RNA could be used to synthesize a cDNA probe and resulted in a gene expression profile that was similar to one produced using 5 micrograms total RNA, even though the overall signal intensity was lower. These experiments illustrate that microarray technology holds great promise for the use of limited clinical samples in the diagnostic setting.     

Statistical methods for microarray assays

The paper shortly reviews statistical methods used in the area of DNA microarray studies. All stages of the experiment are taken into account: planning, data collection, data preprocessing, analysis and validation. Among the methods of data analysis, the algorithms for estimating differential expression, multivariate approaches, clustering methods, as well as classification and discrimination are reviewed. The need is stressed for routine statistical data processing protocols and for the search of links of microarray data analysis with quantitative genetic models.     

Bayesian models for pooling microarray studies with multiple sources of replications

Background  

Biologists often conduct multiple but different cDNA microarray studies that all target the same biological system or pathway. Within each study, replicate slides within repeated identical experiments are often produced. Pooling information across studies can help more accurately identify true target genes. Here, we introduce a method to integrate multiple independent studies efficiently.