Optimization of Single-Base-Pair Mismatch Discrimination in Oligonucleotide Microarrays

The discrimination between perfect-match and single-base-pair-mismatched nucleic acid duplexes was investigated by using oligonucleotide DNA microarrays and nonequilibrium dissociation rates (melting profiles). DNA and RNA versions of two synthetic targets corresponding to the 16S rRNA sequences of Staphylococcus epidermidis (38 nucleotides) and Nitrosomonas eutropha (39 nucleotides) were hybridized to perfect-match probes (18-mer and 19-mer) and to a set of probes having all possible single-base-pair mismatches. The melting profiles of all probe-target duplexes were determined in parallel by using an imposed temperature step gradient. We derived an optimum wash temperature for each probe and target by using a simple formula to calculate a discrimination index for each temperature of the step gradient. This optimum corresponded to the output of an independent analysis using a customized neural network program. These results together provide an experimental and analytical framework for optimizing mismatch discrimination among all probes on a DNA microarray.     

Optimization of single-base-pair mismatch discrimination in oligonucleotide microarrays

The discrimination between perfect-match and single-base-pair-mismatched nucleic acid duplexes was investigated by using oligonucleotide DNA microarrays and nonequilibrium dissociation rates (melting profiles). DNA and RNA versions of two synthetic targets corresponding to the 16S rRNA sequences of Staphylococcus epidermidis (38 nucleotides) and Nitrosomonas eutropha (39 nucleotides) were hybridized to perfect-match probes (18-mer and 19-mer) and to a set of probes having all possible single-base-pair mismatches. The melting profiles of all probe-target duplexes were determined in parallel by using an imposed temperature step gradient. We derived an optimum wash temperature for each probe and target by using a simple formula to calculate a discrimination index for each temperature of the step gradient. This optimum corresponded to the output of an independent analysis using a customized neural network program. These results together provide an experimental and analytical framework for optimizing mismatch discrimination among all probes on a DNA microarray.     

Optimization of an oligonucleotide microchip for microbial identification studies: a non-equilibrium dissociation approach

The utility of a high-density oligonucleotide microarray (microchip) for identifying strains of five closely related bacilli ( Bacillus anthracis , Bacillus cereus , Bacillus mycoides , Bacillus medusa and Bacillus subtilis ) was demonstrated using an approach that compares the non-equilibrium dissociation rates ('melting curves') of all probe–target duplexes simultaneously. For this study, a hierarchical set of 30 oligonucleotide probes targeting the 16S ribosomal RNA of these bacilli at multiple levels of specificity (approximate taxonomic ranks of domain, kingdom, order, genus and species) was designed and immobilized in a high-density matrix of gel pads on a glass slide. Reproducible melting curves for probes with different levels of specificity were obtained using an optimized salt concentration. Clear discrimination between perfect match (PM) and mismatch (MM) duplexes was achieved. By normalizing the signals to an internal standard (a universal probe), a more than twofold discrimination (> 2.4×) was achieved between PM and 1-MM duplexes at the dissociation temperature at which 50% of the probe–target duplexes remained intact. This provided excellent differentiation among representatives of different Bacillus species, both individually and in mixtures of two or three. The overall pattern of hybridization derived from this hierarchical probe set also provided a clear 'chip fingerprint' for each of these closely related Bacillus species.     

Single-base-pair discrimination of terminal mismatches by using oligonucleotide microarrays and neural network analyses.

The effects of single-base-pair near-terminal and terminal mismatches on the dissociation temperature (T(d)) and signal intensity of short DNA duplexes were determined by using oligonucleotide microarrays and neural network (NN) analyses. Two perfect-match probes and 29 probes having a single-base-pair mismatch at positions 1 to 5 from the 5' terminus of the probe were designed to target one of two short sequences representing 16S rRNA. Nonequilibrium dissociation rates (i.e., melting profiles) of all probe-target duplexes were determined simultaneously. Analysis of variance revealed that position of the mismatch, type of mismatch, and formamide concentration significantly affected the T(d) and signal intensity. Increasing the concentration of formamide in the washing buffer decreased the T(d) and signal intensity, and it decreased the variability of the signal. Although T(d)s of probe-target duplexes with mismatches in the first or second position were not significantly different from one another, duplexes with mismatches in the third to fifth positions had significantly lower T(d)s than those with mismatches in the first or second position. The trained NNs predicted the T(d) with high accuracies (R(2) = 0.93). However, the NNs predicted the signal intensity only moderately accurately (R(2) = 0.67), presumably due to increased noise in the signal intensity at low formamide concentrations. Sensitivity analysis revealed that the concentration of formamide explained most (75%) of the variability in T(d)s, followed by position of the mismatch (19%) and type of mismatch (6%). The results suggest that position of the mismatch at or near the 5' terminus plays a greater role in determining the T(d) and signal intensity of duplexes than the type of mismatch.     

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.     

Evaluation of gel-pad oligonucleotide microarray technology by using artificial neural networks

Past studies have suggested that thermal dissociation analysis of nucleic acids hybridized to DNA microarrays would improve discrimination among duplex types by scanning through a broad range of stringency conditions. To more fully constrain the utility of this approach using a previously described gel-pad microarray format, artificial neural networks (NNs) were trained to recognize noisy or low-quality data, as might derive from nonspecific fluorescence, poor hybridization, or compromised data collection. The NNs were trained to classify dissociation profiles (melts) into groups based on selected characteristics (e.g., initial signal intensity, area under the curve) using a data set of 21,044 profiles derived from 186 probes hybridized to a study set of RNA extracted from 32 microbes common to the human oral cavity. Three melt profile groups were identified: one group consisted mostly of ideal melt profiles; another group consisted mostly of poor melt profiles; and, the remainder were difficult to classify. Screening of melting profiles of perfect-match hybrids revealed inconsistencies in the form of melting profiles even for identical probes on the same microarray hybridized to same target rRNA. Approximately 18% of perfect-match duplex types were correctly classified as poor. Experimental variability and deviation from ideal melt behavior were shown to be attributable primarily to a method of local background subtraction that was very sensitive to displacement of the grid frames used for image capture (both determined by the image analysis system) and duplexes with low binding constants. Additional results showed that long RNA fragments limit the discriminating power among duplex types.     

DNA fragment conformations IV - Helix-coil transition and conformation of d-CCATGG in aqueous solution by 1H-NMR spectroscopy.

The helix-coil transition and conformation of d-CCATGG were investigated using 1H-NMR spectroscopy at various frequencies (90, 276, 400 MHz). The changes in the chemical shifts and linewidths of imino protons between 5 degrees and 35 degrees C show that the d-CCATGG fraying process consists of two stages: the external dC.dG base pairs open at first, th internal dC.dG and central dA.dT base pairs then open simultaneously at higher temperatures similar to the case of d-ACATGT. The midpoint temperatures, the helix and coil proportions and the dissociation constant were determined from the sigma = f(t degree) curves of the base and sugar protons. The results indicate that the midpoint temperature increases with the number of the dG.dC base pair in a given size sequence, while the dissociation enthalpy appears to be independent. The difference between the T1 value of a base proton of the external and internal residues of the same nature is found to be a good criterion for base proton assignment. The high predominance of the S conformation for all residues shows that d-CCATGG duplexes adopt the B-helical conformation.     

The use of oligodeoxynucleotide probes in chaotrope-based hybridization solutions.

Hybridization solutions containing chaotropes may be used to modulate the thermal stability (Tm or Td) of oligodeoxynucleotide (ODN) duplexes or hybrids over a 90 degrees C range. Modulation of Td allows formulation of hybridization solutions that permit ambient temperature hybridization using most combinations of probe length, probe composition, target type, and facilitates development of convenient and rapid assay formats. The conditions required to achieve ODN duplex fidelity, and optimal yields of hybridized product, are described for trichloroacetate, thiocyanate, guanidinium salts and other chaotropic salts. The effects of different solid supports on Td are described. Also, a method is presented that uses chaotropic compounds to reduce background arising from signal ODN probes in a sandwich assay hybridization format.     

Specific and nonspecific hybridization of oligonucleotide probes on microarrays

Gene expression analysis by means of microarrays is based on the sequence-specific binding of RNA to DNA oligonucleotide probes and its measurement using fluorescent labels. The binding of RNA fragments involving sequences other than the intended target is problematic because it adds a chemical background to the signal, which is not related to the expression degree of the target gene. The article presents a molecular signature of specific and nonspecific hybridization with potential consequences for gene expression analysis. We analyzed the signal intensities of perfect match (PM) and mismatch (MM) probes of GeneChip microarrays to specify the effect of specific and nonspecific hybridization. We found that these events give rise to different relations between the PM and MM intensities as function of the middle base of the PM, namely a triplet-like (C > G approximately T > A > 0) and a duplet-like (C approximately T > 0 > G approximately A) pattern of the PM-MM log-intensity difference upon binding of specific and nonspecific RNA fragments, respectively. The systematic behavior of the intensity difference can be rationalized on the level of basepairings of DNA/RNA oligonucleotide duplexes in the middle of the probe sequence. Nonspecific binding is characterized by the reversal of the central Watson-Crick (WC) pairing for each PM/MM probe pair, whereas specific binding refers to the combination of a WC and a self-complementary (SC) pairing in PM and MM probes, respectively. The Gibbs free energy contribution of WC pairs to duplex stability is asymmetric for purines and pyrimidines of the PM and decreases according to C > G approximately T > A. SC pairings on the average only weakly contribute to duplex stability. The intensity of complementary MM introduces a systematic source of variation which decreases the precision of expression measures based on the MM intensities.     

Single-Base-Pair Discrimination of Terminal Mismatches by Using Oligonucleotide Microarrays and Neural Network Analyses

The effects of single-base-pair near-terminal and terminal mismatches on the dissociation temperature (T_d) and signal intensity of short DNA duplexes were determined by using oligonucleotide microarrays and neural network (NN) analyses. Two perfect-match probes and 29 probes having a single-base-pair mismatch at positions 1 to 5 from the 5′ terminus of the probe were designed to target one of two short sequences representing 16S rRNA. Nonequilibrium dissociation rates (i.e., melting profiles) of all probe-target duplexes were determined simultaneously. Analysis of variance revealed that position of the mismatch, type of mismatch, and formamide concentration significantly affected the T_d and signal intensity. Increasing the concentration of formamide in the washing buffer decreased the T_d and signal intensity, and it decreased the variability of the signal. Although T_ds of probe-target duplexes with mismatches in the first or second position were not significantly different from one another, duplexes with mismatches in the third to fifth positions had significantly lower T_ds than those with mismatches in the first or second position. The trained NNs predicted the T_d with high accuracies (R² = 0.93). However, the NNs predicted the signal intensity only moderately accurately (R² = 0.67), presumably due to increased noise in the signal intensity at low formamide concentrations. Sensitivity analysis revealed that the concentration of formamide explained most (75%) of the variability in T_ds, followed by position of the mismatch (19%) and type of mismatch (6%). The results suggest that position of the mismatch at or near the 5′ terminus plays a greater role in determining the T_d and signal intensity of duplexes than the type of mismatch.     

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.     

Discrimination between perfect and mismatched duplexes with oligonucleotide gel microchips: role of thermodynamic and kinetic effects during hybridization

The efficiency of discrimination between perfect and mismatched duplexes during hybridization on microchips depends on the concentrations of target DNA in solution and immobilized probes, buffer composition, and temperature of hybridization and is determined by both thermodynamic relationships and hybridization kinetics. In this work, optimal conditions of discrimination were studied using hybridization of fluorescently labeled target DNA with custom-made gel-based oligonucleotide microchips. The higher the concentration of immobilized probes and the higher the association constant, the higher the concentration of the formed duplexes and the stronger the corresponding fluorescence signal, but, simultaneously, the longer the time needed to reach equilibrium. Since mismatched duplexes hybridize faster than their perfect counterparts, perfect-to-mismatch signal ratio is lower in transient regime, and short hybridization times may hamper the detection of mutations. The saturation time can be shortened by decreasing the probe concentration or augmenting the gel porosity. This improves the detection of mutations in transient regime. It is shown that the decrease in the initial concentration of oligonucleotide probes by an order of magnitude causes only 1.5-2.5-fold decrease of fluorescence signals after hybridization of perfect duplexes for 3-12 h. At the same time, these conditions improve the discrimination between perfect and mismatched duplexes more than two-fold. A similar improvement may be obtained using an optimized dissociation procedure.     

Enhanced mismatch selectivity of T4 DNA ligase far above the probe: Target duplex dissociation temperature

T4 DNA ligase is a widely used ligase in many applications; yet in single nucleotide polymorphism analysis, it has been found generally lacking owing to its tendency to ligate mismatches quite efficiently. To address this lack of selectivity, we explored the effect of temperature on the selectivity of the ligase in discriminating single base pair mismatches at the 3′‐terminus of the ligating strand using short ligation probes (9‐mers). Remarkably, we observe outstanding selectivities when the assay temperature is increased to 7 °C to 13 °C above the dissociation temperature of the matched probe:target duplexes using commercially available enzyme at low concentration. Higher enzyme concentration shifts the temperature range to 13 °C to 19 °C above the probe:target dissociation temperatures. Finally, substituting the 5′‐phosphate terminus with an abasic nucleotide decreases the optimal temperature range to 7 °C to 10 °C above the matched probe:target duplex. We compare the temperature dependence of the T4 DNA ligase catalyzed ligation and a nonenzymatic ligation system to contrast the origin of their modes of selectivity. For the latter, temperatures above the probe:target duplex dissociation lead to lower ligation conversions even for the perfect matched system. This difference between the two ligation systems reveals the uniqueness of the T4 DNA ligase's ability to maintain excellent ligation yields for the matched system at elevated temperatures. Although our observations are consistent with previous mechanistic work on T4 DNA ligase, by mapping out the temperature dependence for different ligase concentrations and probe modifications, we identify simple strategies for introducing greater selectivity into SNP discrimination based on ligation yields.     

A comparison of three molecular typing methods for the discrimination of Salmonella enterica serovar Infantis

Seventy-six epidemiologically unrelated Salmonella enterica serovar Infantis (S. Infantis) isolates were typed by pulsed-field gel electrophoresis (PFGE), multiple amplification of phage loci typing (MAPLT) and multiple-locus variable-number tandem-repeat analysis (MLVA). PFGE, using the restriction endonuclease XbaI, generated 23 different profiles for the 76 isolates (DI=0.848). MAPLT was undertaken using a combination of 11 primer sets based on bacteriophage sequences and generated 28 different profiles (DI=0.938). By contrast, MLVA only produced nine profiles (DI=0.668) with 13 different primer sets, including the five primer sets routinely used for S. Typhimurium typing. Reducing the number of MAPLT primer sets to four still provided a diversity index of 0.838. All three typing methods revealed two distinct lineages of S. Infantis, with most isolates demonstrating genetic traits of either lineage but not both. The results demonstrate that MAPLT can potentially provide greater discrimination and separation of S. Infantis isolates than both PFGE and MLVA. Furthermore, MAPLT data can be generated much more rapidly and with reduced labour input than PFGE and without the need for expensive PFGE electrophoresis equipment, nor does it require capillary sequencing of PCR fragments to accurately determine PCR fragment lengths as is the case with MLVA.     

A comparison of solution and membrane-bound DNA × DNA hybridization,as used to infer phylogeny

A new method of membrane-bound DNA × DNA hybridization was devised to accommodate the study of small quantities of DNA obtained from museum specimens for phylogeny reconstruction. Membranebound, single-stranded target genomic DNAs were competitively hybridized with a total genomic DNA probe to form hybrid duplexes required for the DNA dissociation experiments. We compared the thermal elution profiles derived from dissociating duplexes made with probes of whole genomic, single-copy, and repetitive DNA, as well as solution DNA × DNA hybridization using sc tracer. Quantitatively, pairwise indices of genetic distance derived from duplexes made with genomic probes depended entirely on hybridization of repetitive sequences, but a parallel set of experiments using repetitive and sc probes produced qualitatively similar results. The indices of genetic distance generated by the membrane-bound hybrids form an internally consistent, resolved tree which is in agreement with the solution DNA × DNA hybridization trials and traditional views of the phylogeny of the taxa under study.Correspondence to: P. Houde     

Enhanced thermal stability and mismatch discrimination of mutation-carrying DNA duplexes and their kinetic and thermodynamic properties in microchannel laminar flow

This article reports the enhancement of thermal stability involving normal duplex and mutation-carrying DNA duplexes in microchannel laminar flow. The application of an in-house temperature-controllable microchannel-type flow cell is demonstrated for improved discrimination of mismatch base pairs such as A-G and T-G that are difficult to distinguish due to the rather small thermal destabilizations. Enhancement in thermal stability is reflected by an increased thermal melting temperature achieved in microchannel laminar flow as compared with batch reactions. To examine the kinetics and thermodynamics of duplex-coil equilibrium of DNA oligomers, denaturation-renaturation hysteresis curves were measured. The influence of microchannel laminar flow on DNA base mismatch analysis was described from the kinetic and thermodynamic perspectives. An increasing trend was observed for association rate constant as flow rate increased. In contrast, an apparent decrease in dissociation rate constant was observed with increasing flow rate. The magnitudes of the activation energies of dissociation were nearly constant for both the batch and microchannel laminar flow systems at all flow rates. In contrast, the magnitudes of activation energies of association decreased as flow rate increased. These results clearly show how microchannel laminar flow induces change in reaction rate by effecting change in activation energy. We anticipate, therefore, that this approach based on microchannel laminar flow system holds great promise for improved mismatch discrimination in DNA analyses, particularly on single-base-pair mismatch, by pronouncedly enhancing thermal stability.     

Dissociation analysis in polymerase chain reaction and 1X SSC buffer as a prerequisite for selection of 13mer microarray probe sets with uniform hybridization behavior

Homogeneous probe hybridization is a prerequisite for the robust design of microarrays. Elaborate algorithms were developed to select for probe sets with uniform melting temperatures (Tm). However, at least short oligonucleotides (<20 bp) show large variation in the on-chip hybridization efficiency even if designed with state-of-the-art algorithms. This variation can be explained by steric effects and interferences on the solid surface as well as by chemical conditions that may deviate from conditions used to develop the algorithms. We designed 412 random 13mer duplexes to study the differences between the Tm of nearest-neighbor algorithms and the Tm values measured by dissociation analysis in polymerase chain reaction (PCR) and 1X SSC buffer. We tested the effects of theoretical vs empirical Tm values on the hybridization variation of 40 duplexes on-chip. Although the empirical approach resulted in a slightly better prediction of hybridization efficiency, less than one-fifth (17%) of the observed variation could be explained by the factor Tm alone. We conclude that state-of-the-art algorithms can be used for a first selection of short oligonucleotide probes, but that it is then necessary to perform an on-chip selection to obtain a probe set with a uniform hybridization behavior.     

Titration of variant DNA sequences differing by a single point-mutation by selective dot-blot hybridization with synthetic oligonucleotides

A dot-blot hybridization procedure with synthetic oligonucleotide probes is reported, which allows the quantitative titration in genomic DNA of variant forms of repeated genes differing by a single nucleotide change. It involves the utilization of a pair of 22-base long oligonucleotides matching the two variant sequences and the choice of an hybridization temperature very close to the Td of the oligonucleotide/DNA duplexes. The selectivity is achieved through a competition between the cognate labeled and the non-cognate unlabeled probes in the hybridization mixture.     

On-chip non-equilibrium dissociation curves and dissociation rate constants as methods to assess specificity of oligonucleotide probes

Nucleic acid hybridization serves as backbone for many high-throughput systems for detection, expression analysis, comparative genomics and re-sequencing. Specificity of hybridization between probes and intended targets is always critical. Approaches to ensure and evaluate specificity include use of mismatch probes, obtaining dissociation curves rather than single temperature hybridizations, and comparative hybridizations. In this study, we quantify effects of mismatch type and position on intensity of hybridization signals and provide a new approach based on dissociation rate constants to evaluate specificity of hybridized signals in complex target mixtures. Using an extensive set of 18mer oligonucleotide probes on an in situ synthesized biochip platform, we demonstrate that mismatches in the center of the probe are more discriminating than mismatches toward the extremities of the probe and mismatches toward the attached end are less discriminating than those toward the loose end. The observed destabilizing effect of a mismatch type agreed in general with predictions using the nearest neighbor model. Use of a new parameter, specific dissociation temperature (T_d-w, temperature of maximum specific dissociation rate constant), obtained from probe–target duplex dissociation profiles considerably improved the evaluation of specificity. These results have broad implications for hybridization data obtained from complex mixtures of nucleic acids.