Statistical design of reverse dye microarrays

MOTIVATION: In cDNA microarray experiments all samples are labelled with either Cy3 dye or Cy5 dye. Certain genes exhibit dye bias-a tendency to bind more efficiently to one of the dyes. The common reference design avoids the problem of dye bias by running all arrays 'forward', so that the samples being compared are always labelled with the same dye. But comparison of samples labelled with different dyes is sometimes of interest. In these situations, it is necessary to run some arrays 'reverse'-with the dye labelling reversed-in order to correct for the dye bias. The design of these experiments will impact one's ability to identify genes that are differentially expressed in different tissues or conditions. We address the design issue of how many specimens are needed, how many forward and reverse labelled arrays to perform, and how to optimally assign Cy3 and Cy5 labels to the specimens. RESULTS: We consider three types of experiments for which some reverse labelling is needed: paired samples, samples from two predefined groups, and reference design data when comparison with the reference is of interest. We present simple probability models for the data, derive optimal estimators for relative gene expression, and compare the efficiency of the estimators for a range of designs. In each case, we present the optimal design and sample size formulas. We show that reverse labelling of individual arrays is generally not required.     

A new approach to intensity-dependent normalization of two-channel microarrays

A two-channel microarray measures the relative expression levels of thousands of genes from a pair of biological samples. In order to reliably compare gene expression levels between and within arrays, it is necessary to remove systematic errors that distort the biological signal of interest. The standard for accomplishing this is smoothing "MA-plots" to remove intensity-dependent dye bias and array-specific effects. However, MA methods require strong assumptions, which limit their general applicability. We review these assumptions and derive several practical scenarios in which they fail. The "dye-swap" normalization method has been much less frequently used because it requires two arrays per pair of samples. We show that a dye-swap is accurate under general assumptions, even under intensity-dependent dye bias, and that a dye-swap removes dye bias from a single pair of samples in general. Based on a flexible model of the relationship between mRNA amount and single-channel fluorescence intensity, we demonstrate the general applicability of a dye-swap approach. We then propose a common array dye-swap (CADS) method for the normalization of two-channel microarrays. We show that CADS removes both dye bias and array-specific effects, and preserves the true differential expression signal for every gene under the assumptions of the model.     

Evaluation of the gene-specific dye bias in cDNA microarray experiments

MOTIVATION: In cDNA microarray experiments all samples are labeled with either Cy3 or Cy5. Systematic and gene-specific dye bias effects have been observed in dual-color experiments. In contrast to systematic effects which can be corrected by a normalization method, the gene-specific dye bias is not completely suppressed and may alter the conclusions about the differentially expressed genes. METHODS: The gene-specific dye bias is taken into account using an analysis of variance model. We propose an index, named label bias index, to measure the gene-specific dye bias. It requires at least two self-self hybridization cDNA microarrays. RESULTS: After lowess normalization we have found that the gene-specific dye bias is the major source of experimental variability between replicates. The ratio (R/G) may exceed 2. As a consequence false positive genes may be found in direct comparison without dye-swap. The stability of this artifact and its consequences on gene variance and on direct or indirect comparisons are addressed. AVAILABILITY: http://www.inapg.inra.fr/ens_rech/mathinfo/recherche/mathematique     

SynteBase/SynteView: a tool to visualize gene order conservation in prokaryotic genomes

Background

Most microarray studies are made using labelling with one or two dyes which allows the hybridization of one or two samples on the same slide. In such experiments, the most frequently used dyes areCy3 andCy5. Recent improvements in the technology (dye-labelling, scanner and, image analysis) allow hybridization up to four samples simultaneously. The two additional dyes areAlexa488 andAlexa494. The triple-target or four-target technology is very promising, since it allows more flexibility in the design of experiments, an increase in the statistical power when comparing gene expressions induced by different conditions and a scaled down number of slides. However, there have been few methods proposed for statistical analysis of such data. Moreover the lowess correction of the global dye effect is available for only two-color experiments, and even if its application can be derived, it does not allow simultaneous correction of the raw data.

Results

We propose a two-step normalization procedure for triple-target experiments. First the dye bleeding is evaluated and corrected if necessary. Then the signal in each channel is normalized using a generalized lowess procedure to correct a global dye bias. The normalization procedure is validated using triple-self experiments and by comparing the results of triple-target and two-color experiments. Although the focus is on triple-target microarrays, the proposed method can be used to normalizepdifferently labelled targets co-hybridized on a same array, for any value ofpgreater than 2.

Conclusion

The proposed normalization procedure is effective: the technical biases are reduced, the number of false positives is under control in the analysis of differentially expressed genes, and the triple-target experiments are more powerful than the corresponding two-color experiments. There is room for improving the microarray experiments by simultaneously hybridizing more than two samples.     

Functional analysis and comparative genomics of expressed sequence tags from the lycophyte Selaginella moellendorffii

Background

DNA microarrays are widely used in gene expression analyses. To increase throughput and minimize costs without reducing gene expression data obtained, we investigated whether four mRNA samples can be analyzed simultaneously by applying four different fluorescent dyes.

Results

Following tests for cross-talk of fluorescence signals, Alexa 488, Alexa 594, Cyanine 3 and Cyanine 5 were selected for hybridizations. For self-hybridizations, a single RNA sample was labelled with all dyes and hybridized on commercial cDNA arrays or on in-house spotted oligonucleotide arrays. Correlation coefficients for all combinations of dyes were above 0.9 on the cDNA array. On the oligonucleotide array they were above 0.8, except combinations with Alexa 488, which were approximately 0.5. Standard deviation of expression differences for replicate spots were similar on the cDNA array for all dye combinations, but on the oligonucleotide array combinations with Alexa 488 showed a higher variation.

Conclusion

In conclusion, the four dyes can be used simultaneously for gene expression experiments on the tested cDNA array, but only three dyes can be used on the tested oligonucleotide array. This was confirmed by hybridizations of control with test samples, as all combinations returned similar numbers of differentially expressed genes with comparable effects on gene expression.     

Comparison of algorithms for the analysis of Affymetrix microarray data as evaluated by co-expression of genes in known operons

Oligonucleotide microarrays are an informative tool to elucidate gene regulatory networks. In order for gene expression levels to be comparable across microarrays, normalization procedures have to be invoked. A large number of methods have been described to correct for systematic biases in microarray experiments. The performance of these methods has been tested only to a limited extend. Here, we evaluate two different types of microarray analyses: (i) the same gene in replicate samples and (ii) different, but co-expressed genes in the same sample. The reliability of the latter analysis needs to be determined for the analysis of regulatory networks and our report is the first attempt to evaluate for the accuracy of different microarray normalization methods in this respect. Consistent with previous results we observed a large effect of the normalization method on the outcome of the expression analyses. Our analyses indicate that different normalization methods should be performed depending on whether a study is aiming to detect differential gene expression between independent samples or whether co-expressed genes should be identified. We make recommendations about the most appropriate method to use.     

Comparison of Alexa Fluor<Superscript>®</Superscript> and CyDye<Superscript>™</Superscript> for practical DNA microarray use

Microarrays are a powerful tool for comparison and understanding of gene expression levels in healthy and diseased states. The method relies upon the assumption that signals from microarray features are a reflection of relative gene expression levels of the cell types under investigation. It has previously been reported that the classical fluorescent dyes used for microarray technology, Cy3 and Cy5, are not ideal due to the decreased stability and fluorescence intensity of the Cy5 dye relative to the Cy3, such that dye bias is an accepted phenomena necessitating dye swap experimental protocols and analysis of differential dye affects. The incentive to find new fluorophores is based on alleviating the problem of dye bias through synonymous performance between counterpart dyes. Alexa Fluor 555 and Alexa Fluor 647 are increasingly promoted as replacements for CyDye in microarray experiments. Performance relates to the molecular and steric similarities, which will vary for each new pair of dyes as well as the spectral integrity for the specific application required. Comparative analysis of the performance of these two competitive dye pairs in practical microarray applications is warranted towards this end. The findings of our study showed that both dye pairs were comparable but that conventional CyDye resulted in significantly higher signal intensities (P < 0.05) and signal minus background levels (P < 0.05) with no significant difference in background values (P > 0.05). This translated to greater levels of differential gene expression with CyDye than with the Alexa Fluor counterparts. However, CyDye fluorophores and in particular Cy5, were found to be less photostable over time and following repeated scans in microarray experiments. These results suggest that precautions against potential dye affects will continue to be necessary and that no one dye pair negates this need.     

Normalization of two-channel microarray experiments: a semiparametric approach

MOTIVATION: An important underlying assumption of any experiment is that the experimental subjects are similar across levels of the treatment variable, so that changes in the response variable can be attributed to exposure to the treatment under study. This assumption is often not valid in the analysis of a microarray experiment due to systematic biases in the measured expression levels related to experimental factors such as spot location (often referred to as a print-tip effect), arrays, dyes, and various interactions of these effects. Thus, normalization is a critical initial step in the analysis of a microarray experiment, where the objective is to balance the individual signal intensity levels across the experimental factors, while maintaining the effect due to the treatment under investigation. RESULTS: Various normalization strategies have been developed including log-median centering, analysis of variance modeling, and local regression smoothing methods for removing linear and/or intensity-dependent systematic effects in two-channel microarray experiments. We describe a method that incorporates many of these into a single strategy, referred to as two-channel fastlo, and is derived from a normalization procedure that was developed for single-channel arrays. The proposed normalization procedure is applied to a two-channel dose-response experiment.     

Correcting for gene-specific dye bias in DNA microarrays using the method of maximum likelihood

MOTIVATION: In two-color microarray experiments, well-known differences exist in the labeling and hybridization efficiency of Cy3 and Cy5 dyes. Previous reports have revealed that these differences can vary on a gene-by-gene basis, an effect termed gene-specific dye bias. If uncorrected, this bias can influence the determination of differentially expressed genes. RESULTS: We show that the magnitude of the bias scales multiplicatively with signal intensity and is dependent on which nucleotide has been conjugated to the fluorescent dye. A method is proposed to account for gene-specific dye bias within a maximum-likelihood error modeling framework. Using two different labeling schemes, we show that correcting for gene-specific dye bias results in the superior identification of differentially expressed genes within this framework. Improvement is also possible in related ANOVA approaches. AVAILABILITY: A software implementation of this procedure is freely available at http://cellcircuits.org/VERA     

An adaptive method for cDNA microarray normalization

Background  

Normalization is a critical step in analysis of gene expression profiles. For dual-labeled arrays, global normalization assumes that the majority of the genes on the array are non-differentially expressed between the two channels and that the number of over-expressed genes approximately equals the number of under-expressed genes. These assumptions can be inappropriate for custom arrays or arrays in which the reference RNA is very different from the experimental samples.     

寡聚核苷酸芯片表达谱系统偏移的校正算法

多芯片对比实验中,由于多方面的变异因素,使得芯片间存在明显的系统偏移.因此,芯片表达谱数据的校正处理是关键的数据预处理步骤.当前,已经提出了很多校正算法,比如:比例常数校正、非线性校正、分位数校正等.提出了一种新的校正算法.在选择的最小秩差异探针集上,进行非线性M-A校正.并采用迭代策略减弱基准芯片方法对基准芯片选择的敏感性.在标准测试集上,同几种已知的方法进行了对比分析.     

Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation

There are many sources of systematic variation in cDNA microarray experiments which affect the measured gene expression levels (e.g. differences in labeling efficiency between the two fluorescent dyes). The term normalization refers to the process of removing such variation. A constant adjustment is often used to force the distribution of the intensity log ratios to have a median of zero for each slide. However, such global normalization approaches are not adequate in situations where dye biases can depend on spot overall intensity and/or spatial location within the array. This article proposes normalization methods that are based on robust local regression and account for intensity and spatial dependence in dye biases for different types of cDNA microarray experiments. The selection of appropriate controls for normalization is discussed and a novel set of controls (microarray sample pool, MSP) is introduced to aid in intensity-dependent normalization. Lastly, to allow for comparisons of expression levels across slides, a robust method based on maximum likelihood estimation is proposed to adjust for scale differences among slides.     

Biosorption of a reactive dye (Rhodamine-B) from an aqueous solution using dried biomass of activated sludge

Low cost, locally available biomaterial was tested for its ability to remove reactive dyes from aqueous solution. Granules prepared from dried activated sludge (DAS) were utilized as a sorbent for the uptake of Rhodamine-B (Rh-B) dye. The effects of various experimental parameters (dye concentration, sludge concentrations, swelling, pretreatment and other factors) were investigated and optimal experimental conditions were ascertained. Nearly 15min was required for the equilibrium adsorption, and Rh-B dyes could be removed effectively. Dye removal performance of Rh-B and DAS increased with increasing concentrations. The acid pretreated biomass exhibited a slightly better biosorption capacity than alkali pretreated or non-pretreated biomass. The optimum swelling time for dye adsorption of the DAS within the swelling time range studied was 12h. Both the Freundlich and Langmuir isotherm models could describe the adsorption equilibrium of the reactive dye onto the activated sludge with the Langmuir isotherm showing the better agreement of the two. Second-order kinetic models confirmed the agreement.     

A practical question based on cross-platform microarray data normalization: are BOEC more like large vessel or microvascular endothelial cells or neither of them?

Since the available microarray data of BOEC (human blood outgrowth endothelial cells), large vessel, and microvascular endothelial cells were from two different platforms, a working cross-platform normalization method was needed to make these data comparable. With six HUVEC (human umbilical vein endothelial cells) samples hybridized on two-channel cDNA arrays and six HUVEC samples on Affymetrix arrays, 64 possible combinations of a three-step normalization procedure were investigated to search for the best normalization method, which was selected, based on two criteria measuring the extent to which expression profiles of biological samples of the same cell type arrayed on two platforms were indistinguishable. Next, three discriminative gene lists between the large vessel and the microvascular endothelial cells were achieved by SAM (significant analysis of microarrays), PAM (prediction analysis for microarrays), and a combination of SAM and PAM lists. The final discriminative gene list was selected by SVM (support vector machine). Based on this discriminative gene list, SVM classification analysis with best tuning parameters and 10,000 times of validations showed that BOEC were far from large vessel cells, they either formed their own class, or fell into the microvascular class. Based on all the common genes between the two platforms, SVM analysis further confirmed this conclusion.     

Van Gieson's picrofuchsin. The staining mechanisms for collagen and cytoplasm,and an examination of the dye diffusion rate model of differential staining

The staining mechanism of van Gieson's picrofuchsin was studied by use of simple protein model systems and tissue sections, and by spectrophotometry and dialysis experiments. At the endpoint of the staining reaction (equilibrium) cytoplasm is yellow. Dye dilution experiments demonstrated that the highest affinity in the tissue section — picrofuchsin system is between binding sites in cytoplasmic protein and acid fuchsin. Nevertheless sections that were first stained in acid fuchsin (AcF) and then in picrofuchsin ended up with cytoplasm stained yellow. It was concluded that differences in the dye diffusion rates and differences in the permeability of tissue components cannot be invoked to explain the differential staining result. Model experiments with dissolved proteins demonstrated a positive relationship between protein concentration and uptake of picric acid (PA) from picrofuchsin. From this and experiments with additives (sodium dodecylsulphate, urea etc.) and organic solvents, it is proposed that coagulant interchain cross-linking at the high protein concentration of the cytoplasm masks potential dye-binding sites. This affects high affinity dyes with multiple binding sites more than small dyes, and so puts AcF at a disadvantage compared to PA. Staining of non-collagen proteins is mainly by hydrophobic bonding, involving ionic attractions, apolar bonds, and release of water. This mode of binding is relatively strong, decreases swelling and leads to slow dye exchange. Dye binding to collagen is mostly by hydrogen bonds, but in aqueous dye solvent nonpolar residues and charged residues may also participate. This structure remains relatively open during and after dye-binding, and the bound dye ions are therefore easily exchanged for other dye ions. Address at which the main part of the investigation was carried out