Characterization of mismatch and high-signal intensity probes associated with Affymetrix genechips

MOTIVATION: For Affymetrix microarray platforms, gene expression is determined by computing the difference in signal intensities between perfect match (PM) and mismatch (MM) probesets. Although the use of PM is not controversial, MM probesets have been associated with variance and ultimately inaccurate gene expression calls. A principal focus of this study was to investigate the nature of the MM signal intensities and demonstrate its contribution to the experimental results. RESULTS: While most MM intensities were likely associated with random noise, a subset of approximately 20% (99,485) of the MM probes displayed relatively high signal intensities to the corresponding PM probes (MM > PM) in a non-random fashion; 13,440 of these probes demonstrated exceptionally high 'outlier' intensities. About 15,938 PM probes also demonstrated exceptionally high outlier intensities consistently across all hybridizations. About 92% of the MM > PM probes had either a dThymidine (dT) or a dCytidine (dC) at the 13th position of the probe sequence. MM and PM probes displaying extremely high outlier intensities contained high dC rich nucleotides, and low dA contents at other nucleotides positions along the 25mer probe sequence. Differentially expressed genes generated using Genechip Operating System (GCOS) or modified PM-only methods were also examined. Of those candidate genes identified in the PM-only method, 157 of them were designated by GCOS as absent across all datasets and many others contained probes with MM > PM signal intensities. Our data suggests that MM intensity from PM signal can be a major source of error analysis, leading to fewer potentially biologically important candidate genes. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Evaluating single-base-pair discriminating capability of planar oligonucleotide microchips using a non-equilibrium dissociation approach

The capability of planar rRNA-based oligonucleotide microarrays for single-base-pair discrimination was evaluated using an approach that compares the non-equilibrium dissociation profiles and dissociation temperatures (Tds) of all probe-target duplexes simultaneously. Three sets of 16S rRNA gene specific probes at different levels of specificity were used along with their counter probes for individual sets having either one or two mismatches (MM) to their targets at specific external (next to terminus) and various internal positions. Criteria based on the Td approach and a discrimination index (DI) were proven to be competent in discriminating PM from internal MM duplexes, but not always for external MM duplexes. Maximal DI for separating PM duplexes from ones with two and one internal MM usually occurred at temperatures approximately 5-10 degrees C and 10-15 degrees C, respectively, higher than the Tds of the PM duplexes. Washing buffer type and salt concentration, and MM number and position were shown statistically to affect dissociation profiles, Td, and single-base-pair discriminating capability. The reusability potential of the planar microchip was further demonstrated.     

Position of the fluorescent label is a crucial factor determining signal intensity in microarray hybridizations

A key issue in applications of short oligonucleotide-based microarrays is how to design specific probes with high sensitivity. Some details of the factors affecting microarray hybridization remain unclear, hampering a reliable quantification of target nucleic acids. We have evaluated the effect of the position of the fluorescent label [position of label (POL)] relative to the probe-target duplex on the signal output of oligonucleotide microarrays. End-labelled single-stranded DNA targets of different lengths were used for hybridization with perfect-match oligonucleotide probe sets targeting different positions of the same molecule. Hybridization results illustrated that probes targeting the labelled terminus of the target showed significantly higher signals than probes targeting other regions. This effect was independent of the target gene, the fluorophore and the slide surface chemistry. Comparison of microarray signal patterns of fluorescently end-labelled, fluorescently internally random-labelled and radioactively end-labelled target-DNAs with the same set of oligonucleotide probes identified POL as a critical factor affecting signal intensity rather than binding efficiency. Our observations define a novel determinant for large differences of signal intensities. Application of the POL effect may contribute to better probe design and data interpretation in microarray applications.     

Revision of the nonequilibrium thermal dissociation and stringent washing approaches for identification of mixed nucleic acid targets by microarrays

Microarray experiments typically involve washing steps that remove hybridized nonspecific targets with the purpose of improving the signal-to-noise ratio. The quality of washing ultimately affects downstream analysis of the microarray and interpretation. The paucity of fundamental studies directed towards understanding the dissociation of mixed targets from microarrays makes the development of meaningful washing/dissociation protocols difficult. To fill the void, we examined activation energies and preexponential coefficients of 47 perfect match (PM) and double-mismatch (MM) duplex pairs to discover that there was no statistical difference between the kinetics of the PM and MM duplexes. Based on these findings, we evaluated the nonequilibrium thermal dissociation (NTD) approach, which has been used to identify specific microbial targets in mixed target samples. We found that the major premises for various washing protocols and the NTD approach might be seriously compromised because: (i) nonspecific duplexes do not always dissociate before specific ones, and (ii) the relationship between dissociation rates of the PM and MM duplexes depends on temperature and duplex sequence. Specifically for the NTD, we show that previously suggested use of reference curves, indices of curves and temperature ramps lead to erroneous conclusions.     

Simultaneous quantification of multiple nucleic acid targets in complex rRNA mixtures using high density microarrays and nonspecific hybridization as a source of information

To date, it has been problematic to accurately quantify multiple nucleic acid sequences, representing microbial targets, in multi-target mixtures using oligonucleotide microarrays, primarily due to nonspecific target binding (i.e., cross-hybridization). While some studies ignore the effects of nonspecific binding, other studies have developed approaches to minimize nonspecific binding, such as physical modeling to design highly specific probes, subtracting nonspecific signal using mismatch probes, and/or removing nonspecific duplexes by scanning through a range of wash stringencies. We have developed an alternative approach that, in contrast to previous approaches, uses nonspecific target binding as a source of information. Specifically, the new approach uses hybridization patterns (fingerprints) to quantify specific nucleic acid targets in complex target mixtures. We evaluated the approach by mixing together in vitro transcribed 28S rRNA targets at varying concentrations (up to 1.0 nM), and hybridizing the 24 mixtures to microarrays (n=3160 probes, in duplicate). Three independent Latin-square-designed experiments revealed accurate quantification of the targets. The regression between actual concentration of targets and those determined by the approach were highly positively correlated with high R(2) values (e.g., R(2)=0.90, n=6 targets; R(2)=0.84, n=8 targets; R(2)=0.82, n=10 targets).     

A model of molecular interactions on short oligonucleotide microarrays

High-density short oligonucleotide microarrays have become a widely used tool for measuring gene expression on a large scale. However, details of the mechanism of binding on microarrays remain unclear. Short oligonucleotide probes currently synthesized on microarrays are often ineffective as a result of limited sequence specificity or low sensitivity. Here, we describe a model of binding interactions on microarrays that reveals how probe signals depend on probe sequences and why certain probes are ineffective. The model indicates that the amount of nonspecific binding can be estimated from a simple rule. Using this model, we have developed an improved measure of gene expression for use in data analysis.     

Empirical establishment of oligonucleotide probe design criteria

Criteria for the design of gene-specific and group-specific oligonucleotide probes were established experimentally via an oligonucleotide array that contained perfect match (PM) and mismatch probes (50-mers and 70-mers) based upon four genes. The effects of probe-target identity, continuous stretch, mismatch position, and hybridization free energy on specificity were tested. Little hybridization was observed at a probe-target identity of < or =85% for both 50-mer and 70-mer probes. PM signal intensities (33 to 48%) were detected at a probe-target identity of 94% for 50-mer oligonucleotides and 43 to 55% for 70-mer probes at a probe-target identity of 96%. When the effects of sequence identity and continuous stretch were considered independently, a stretch probe (>15 bases) contributed an additional 9% of the PM signal intensity compared to a nonstretch probe (< or =15 bases) at the same identity level. Cross-hybridization increased as the length of continuous stretch increased. A 35-base stretch for 50-mer probes or a 50-base stretch for 70-mer probes had approximately 55% of the PM signal. Little cross-hybridization was observed for probes with a minimal binding free energy greater than -30 kcal/mol for 50-mer probes or -40 kcal/mol for 70-mer probes. Based on the experimental results, a set of criteria are suggested for the design of gene-specific and group-specific oligonucleotide probes, and the experimentally established criteria should provide valuable information for new software and algorithms for microarray-based studies.     

Characterization of the expression ratio noise structure in high-density oligonucleotide arrays

Background

High-density oligonucleotide microarrays provide a powerful tool for assessing differential mRNA expression levels. Characterizing the noise resulting from the enzymatic and hybridization steps, called type I noise, is essential for attributing significance measures to the differential expression scores. We introduce scoring functions for expression ratios, and associated quality measures. Both the PM (Perfect Match) probes and PM-MM differentials (MM is the single MisMatch) are considered as raw intensities. We then characterize the log-ratio noise structure using robust estimates of their intensity dependent variance.

Results

We show the relationships between the obtained ratios and their quality measures. The complementarity of PM and PM-MM methods is emphasized by the probe sets signal to noise measures. Using a large set of replicate experiments, we demonstrate that the noise structure in the log-ratios very closely follows a local log-normal distribution for both the PM and PM-MM cases. Therefore, significance relative to the type I noise can be quantified reliably using the local STD. We discuss the intensity dependence of the STD and show that ratio scores >1.25 are significant in the mid- to high-intensity range.

Conclusions

The ratio noise structure inherent to high-density oligonucleotide arrays can be well described in terms of local log-normal ratio distributions with characteristic intensity dependence. Therefore, robust estimates of the local STD of these distributions provide a simple and powerful way for assessing significance (relative to type I noise) in differential gene expression. This approach will be helpful for improving the reliability of predictions from hybridization experiments in general.     

Catalysis of RNA cleavage by the Tetrahymena thermophila ribozyme. 1. Kinetic description of the reaction of an RNA substrate complementary to the active site

A ribozyme derived from the intervening sequence (IVS) of the Tetrahymena preribosomal RNA catalyzes a site-specific endonuclease reaction: G2CCCUCUA5 + G in equilibrium with G2CCCUCU + GA5 (G = guanosine). This reaction is analogous to the first step in self-splicing of the pre-rRNA, with the product G2CCCUCU analogous to the 5'-exon. The following mechanistic conclusions have been derived from pre-steady-state and steady-state kinetic measurements at 50 degrees C and neutral pH in the presence of 10 mM Mg2+. The value of kcat/Km = 9 x 10(7) M-1 min-1 for the oligonucleotide substrate with saturating G represents rate-limiting binding. This rate constant for binding is of the order expected for formation of a RNA.RNA duplex between oligonucleotides. (Phylogenetic and mutational analyses have shown that this substrate is recognized by base pairing to a complementary sequence within the IVS). The value of kcat = 0.1 min-1 represents rate-limiting dissociation of the 5'-exon analogue, G2CCCUCU. The product GA5 dissociates first from the ribozyme because of this slow off-rate for G2CCCUCU. The similar binding of the product, G2CCCUCU, and the substrate, G2CCCUCUA5, to the 5'-exon binding site of the ribozyme, with Kd = 1-2 nM, shows that the pA5 portion of the substrate makes no net contribution to binding. Both the substrate and product bind approximately 10(4)-fold (6 kcal/mol) stronger than expected from base pairing with the 5'-exon binding site. Thus, tertiary interactions are involved in binding. Binding of G2CCCUCU and binding of G are independent. These and other data suggest that binding of the oligonucleotide substrate, G2CCCUCUA5, and binding of G are essentially random and independent. The rate constant for reaction of the ternary complex is calculated to be kc approximately equal to 350 min-1, a rate constant that is not reflected in the steady-state rate parameters with saturating G. The simplest interpretation is adopted, in which kc represents the rate of the chemical step. A site-specific endonuclease reaction catalyzed by the Tetrahymena ribozyme in the absence of G was observed; the rate of the chemical step with solvent replacing guanosine, kc(-G) = 0.7 min-1, is approximately 500-fold slower than that with saturating guanosine. The value of kcat/Km = 6 x 10(7) M-1 min-1 for this hydrolysis reaction is only slightly smaller than that with saturating guanosine, because the binding of the oligonucleotide substrate is predominantly rate-limiting in both cases.(ABSTRACT TRUNCATED AT 400 WORDS)     

Empirical Establishment of Oligonucleotide Probe Design Criteria

Criteria for the design of gene-specific and group-specific oligonucleotide probes were established experimentally via an oligonucleotide array that contained perfect match (PM) and mismatch probes (50-mers and 70-mers) based upon four genes. The effects of probe-target identity, continuous stretch, mismatch position, and hybridization free energy on specificity were tested. Little hybridization was observed at a probe-target identity of ≤85% for both 50-mer and 70-mer probes. PM signal intensities (33 to 48%) were detected at a probe-target identity of 94% for 50-mer oligonucleotides and 43 to 55% for 70-mer probes at a probe-target identity of 96%. When the effects of sequence identity and continuous stretch were considered independently, a stretch probe (>15 bases) contributed an additional 9% of the PM signal intensity compared to a nonstretch probe (≤15 bases) at the same identity level. Cross-hybridization increased as the length of continuous stretch increased. A 35-base stretch for 50-mer probes or a 50-base stretch for 70-mer probes had approximately 55% of the PM signal. Little cross-hybridization was observed for probes with a minimal binding free energy greater than −30 kcal/mol for 50-mer probes or −40 kcal/mol for 70-mer probes. Based on the experimental results, a set of criteria are suggested for the design of gene-specific and group-specific oligonucleotide probes, and the experimentally established criteria should provide valuable information for new software and algorithms for microarray-based studies.     

DNA exhibits multi-stranded binding recognition on glass microarrays

In the course of exploring the hybridization properties of glass DNA microarrays, multi-stranded DNA structures were observed that could not be accounted for by classical Watson-Crick base pairing. Non-denatured double-stranded DNA array elements were shown to hybridize to single-stranded (ss)DNA probes. Similarly, ssDNA array elements were shown to bind duplex DNA probes. This led to a series of experiments demonstrating the formation of multi-stranded DNA structures on the surface of microarrays. These structures were observed with a number of heterogeneous sequences, including both purine and pyrimidine bases, with shared sequence identity between the ssDNA and one of the duplex strands. Furthermore, we observed a strong binding preference near the ends of duplexes containing a 3'-homologous strand. We suggest that such binding interactions on cationic solid surfaces could serve as a model for a number of biological processes mediated through multi-stranded DNA.     

Unrestricted hybridization of oligonucleotides to structure-free DNA

The existence of secondary structure in long single-stranded DNA and RNA is a serious obstacle to the practical use of short oligonucleotide probes (<20-mers). Here, we show that replication of a highly structured DNA in the presence of a unique set of dNTP analogues leads to synthesis of daughter DNA with a significantly reduced level of secondary structure. This replicated DNA, composed of 2-aminoadenine, 2-thiothymine, 7-deazaguanine, and cytosine bases, was readily accessible to tiled 8-mer LNA and 15-mer DNA probes, whereas an unmodified version of the same DNA was inaccessible. Importantly, while the base analogues enhanced probe-target stability, they did not significantly reduce the specificity of base pairing. The availability of structure-free DNA targets should facilitate the use of short oligonucleotide probes and promote development of generic oligonucleotide microarrays.     

Immunochemical detection of multiple conformations within a 36 base pair oligonucleotide.

A 36 base pair chimeric oligonucleotide containing a central core of DNA duplex flanked by RNA/DNA hybrid at each end was synthesized. These distinct regions of the oligonucleotide adopt different conformations which were detected with antibody probes. Enzyme linked immunosorbent assays (ELISA) and a gel electrophoresis retardation assay were used to demonstrate the binding of antibodies which recognize B-DNA, Z-DNA and RNA/DNA hybrid. The DNA duplex core of this oligonucleotide adopts the B-conformation in 0.14 M NaCl. In high salt solution (4 M NaCl) the DNA core adopts the Z-conformation. The RNA/DNA hybrid at the ends of the oligomer adopt a conformation which is distinct from both B-DNA and A-RNA.     

Hybridization of 2'-O-methyl and 2'-deoxy molecular beacons to RNA and DNA targets

Molecular beacons are stem-loop hairpin oligonucleotide probes labeled with a fluorescent dye at one end and a fluorescence quencher at the other end; they can differentiate between bound and unbound probes in homogeneous hybridization assays with a high signal-to-background ratio and enhanced specificity compared with linear oligonucleotide probes. However, in performing cellular imaging and quantification of gene expression, degradation of unmodified molecular beacons by endogenous nucleases can significantly limit the detection sensitivity, and results in fluorescence signals unrelated to probe/target hybridization. To substantially reduce nuclease degradation of molecular beacons, it is possible to protect the probe by substituting 2'-O-methyl RNA for DNA. Here we report the analysis of the thermodynamic and kinetic properties of 2'-O-methyl and 2'-deoxy molecular beacons in the presence of RNA and DNA targets. We found that in terms of molecular beacon/target duplex stability, 2'-O-methyl/RNA > 2'-deoxy/RNA > 2'-deoxy/DNA > 2'-O-methyl/DNA. The improved stability of the 2'-O-methyl/RNA duplex was accompanied by a slightly reduced specificity compared with the duplex of 2'-deoxy molecular beacons and RNA targets. However, the 2'-O-methyl molecular beacons hybridized to RNA more quickly than 2'-deoxy molecular beacons. For the pairs tested, the 2'-deoxy-beacon/DNA-target duplex showed the fastest hybridization kinetics. These findings have significant implications for the design and application of molecular beacons.     

Tests of rRNA hybridization to microarrays suggest that hybridization characteristics of oligonucleotide probes for species discrimination cannot be predicted

Hybridization of rRNAs to microarrays is a promising approach for prokaryotic and eukaryotic species identification. Typically, the amount of bound target is measured by fluorescent intensity and it is assumed that the signal intensity is directly related to the target concentration. Using thirteen different eukaryotic LSU rRNA target sequences and 7693 short perfect match oligonucleotide probes, we have assessed current approaches for predicting signal intensities by comparing Gibbs free energy (ΔG°) calculations to experimental results. Our evaluation revealed a poor statistical relationship between predicted and actual intensities. Although signal intensities for a given target varied up to 70-fold, none of the predictors were able to fully explain this variation. Also, no combination of different free energy terms, as assessed by principal component and neural network analyses, provided a reliable predictor of hybridization efficiency. We also examined the effects of single-base pair mismatch (MM) (all possible types and positions) on signal intensities of duplexes. We found that the MM effects differ from those that were predicted from solution-based hybridizations. These results recommend against the application of probe design software tools that use thermodynamic parameters to assess probe quality for species identification. Our results imply that the thermodynamic properties of oligonucleotide hybridization are by far not yet understood.