Model-based probe set optimization for high-performance microarrays

A major challenge in microarray design is the selection of highly specific oligonucleotide probes for all targeted genes of interest, while maintaining thermodynamic uniformity at the hybridization temperature. We introduce a novel microarray design framework (Thermodynamic Model-based Oligo Design Optimizer, TherMODO) that for the first time incorporates a number of advanced modelling features: (i) A model of position-dependent labelling effects that is quantitatively derived from experiment. (ii) Multi-state thermodynamic hybridization models of probe binding behaviour, including potential cross-hybridization reactions. (iii) A fast calibrated sequence-similarity-based heuristic for cross-hybridization prediction supporting large-scale designs. (iv) A novel compound score formulation for the integrated assessment of multiple probe design objectives. In contrast to a greedy search for probes meeting parameter thresholds, this approach permits an optimization at the probe set level and facilitates the selection of highly specific probe candidates while maintaining probe set uniformity. (v) Lastly, a flexible target grouping structure allows easy adaptation of the pipeline to a variety of microarray application scenarios. The algorithm and features are discussed and demonstrated on actual design runs. Source code is available on request.     

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.     

植物病毒检测芯片的杂交条件优化

利用芯片点样仪将5种侵染马铃薯的病毒/类病毒(苜蓿花叶病毒、黄瓜花叶病毒、黄瓜花叶病毒-卫星病毒、马铃薯病毒Y、马铃薯块茎纺锤状类病毒)的保守区寡核苷酸(Oligonucleotide,oligo)探针和PCR探针点样于玻片,并以植物18S rRNA作为内参照制成基因芯片。研究探针浓度、杂交时间、杂交温度以及点样液对芯片杂交的影响,并验证优化后病毒检测芯片的特异性。结果表明,寡核苷酸探针浓度介于5-20 ?mol/L之间对杂交信号强度影响不大,PCR探针浓度与杂交信号强度间呈线性关系;在45℃杂交4 h时,芯片的杂交信号最强,且该条件下进行杂交对两种探针芯片的影响趋势一致;点样液中以DMSO的杂交效果最好。经过整体条件优化后的两种探针芯片在杂交检测上具有较高的特异性,适于检测植物病毒。     

Synthetic oligonucleotide probes deduced from amino acid sequence data. Theoretical and practical considerations

Synthetic probes deduced from amino acid sequence data are widely used to detect cognate coding sequences in libraries of cloned DNA segments. The redundancy of the genetic code dictates that a choice must be made between (1) a mixture of probes reflecting all codon combinations, and (2) a single longer "optimal" probe. The second strategy is examined in detail. The frequency of sequences matching a given probe by chance alone can be determined and also the frequency of sequences closely resembling the probe and contributing to the hybridization background. Gene banks cannot be treated as random associations of the four nucleotides, and probe sequences deduced from amino acid sequence data occur more often than predicted by chance alone. Probe lengths must be increased to confer the necessary specificity. Examination of hybrids formed between unique homologous probes and their cognate targets reveals that short stretches of perfect homology occurring by chance make a significant contribution to the hybridization background. Statistical methods for improving homology are examined, taking human coding sequences as an example, and considerations of codon utilization and dinucleotide frequencies yield an overall homology of greater than 82%. Recommendations for probe design and hybridization are presented, and the choice between using multiple probes reflecting all codon possibilities and a unique optimal probe is discussed.     

Selection of optimal oligonucleotide probes for microarrays using multiple criteria, global alignment and parameter estimation

The oligonucleotide specificity for microarray hybridization can be predicted by its sequence identity to non-targets, continuous stretch to non-targets, and/or binding free energy to non-targets. Most currently available programs only use one or two of these criteria, which may choose ‘false’ specific oligonucleotides or miss ‘true’ optimal probes in a considerable proportion. We have developed a software tool, called CommOligo using new algorithms and all three criteria for selection of optimal oligonucleotide probes. A series of filters, including sequence identity, free energy, continuous stretch, GC content, self-annealing, distance to the 3′-untranslated region (3′-UTR) and melting temperature (T_m), are used to check each possible oligonucleotide. A sequence identity is calculated based on gapped global alignments. A traversal algorithm is used to generate alignments for free energy calculation. The optimal T_m interval is determined based on probe candidates that have passed all other filters. Final probes are picked using a combination of user-configurable piece-wise linear functions and an iterative process. The thresholds for identity, stretch and free energy filters are automatically determined from experimental data by an accessory software tool, CommOligo_PE (CommOligo Parameter Estimator). The program was used to design probes for both whole-genome and highly homologous sequence data. CommOligo and CommOligo_PE are freely available to academic users upon request.     

The effects of mismatches on hybridization in DNA microarrays: determination of nearest neighbor parameters

Quantifying interactions in DNA microarrays is of central importance for a better understanding of their functioning. Hybridization thermodynamics for nucleic acid strands in aqueous solution can be described by the so-called nearest neighbor model, which estimates the hybridization free energy of a given sequence as a sum of dinucleotide terms. Compared with its solution counterparts, hybridization in DNA microarrays may be hindered due to the presence of a solid surface and of a high density of DNA strands. We present here a study aimed at the determination of hybridization free energies in DNA microarrays. Experiments are performed on custom Agilent slides. The solution contains a single oligonucleotide. The microarray contains spots with a perfect matching (PM) complementary sequence and other spots with one or two mismatches (MM) : in total 1006 different probe spots, each replicated 15 times per microarray. The free energy parameters are directly fitted from microarray data. The experiments demonstrate a clear correlation between hybridization free energies in the microarray and in solution. The experiments are fully consistent with the Langmuir model at low intensities, but show a clear deviation at intermediate (non-saturating) intensities. These results provide new interesting insights for the quantification of molecular interactions in DNA microarrays.     

Optimization of the BLASTN substitution matrix for prediction of non-specific DNA microarray hybridization

DNA microarray measurements are susceptible to error caused by non-specific hybridization between a probe and a target (cross-hybridization), or between two targets (bulk-hybridization). Search algorithms such as BLASTN can quickly identify potentially hybridizing sequences. We set out to improve BLASTN accuracy by modifying the substitution matrix and gap penalties. We generated gene expression microarray data for samples in which 1 or 10% of the target mass was an exogenous spike of known sequence. We found that the 10% spike induced 2-fold intensity changes in 3% of the probes, two-third of which were decreases in intensity likely caused by bulk-hybridization. These changes were correlated with similarity between the spike and probe sequences. Interestingly, even very weak similarities tended to induce a change in probe intensity with the 10% spike. Using this data, we optimized the BLASTN substitution matrix to more accurately identify probes susceptible to non-specific hybridization with the spike. Relative to the default substitution matrix, the optimized matrix features a decreased score for A–T base pairs relative to G–C base pairs, resulting in a 5–15% increase in area under the ROC curve for identifying affected probes. This optimized matrix may be useful in the design of microarray probes, and in other BLASTN-based searches for hybridization partners.     

Determination of interactions between structured nucleic acids by fluorescence resonance energy transfer (FRET): selection of target sites for functional nucleic acids.

We previously developed a method for monitoring the integrity of oligonucleotides in vitro and in vivo by quantitating fluorescence resonance energy transfer (FRET) between two different fluorochromes attached to a single oligonucleotide. As an extension of this analysis, we examined changes in the extent of FRET in the presence or absence of target nucleic acids with a specific sequence and a higher-ordered structure. In this system FRET was maximal when probes were free in solution and a decrease in FRET was evidence of successful hybridization. We used a single-stranded oligodeoxyribonucleotide labeled at its 5'-end and its 3'-end with 6-carboxyfluorescein and 6-carboxytetramethylrhodamine, respectively. Incubation of the probe with a single-stranded complementary oligonucleotide reduced the FRET. Moreover, a small change in FRET was also observed when the probe was incubated with an oligonucleotide in which the target site had been embedded in a stable hairpin structure. The decrease in the extent of FRET depended on the length of the stem region of the hairpin structure and also on the higher-ordered structure of the probe. These results indicate that this spectrofluorometric method and FRET probes can be used to estimate the efficacy of hybridization between a probe and its target site within highly ordered structures. This conclusion based on changes in FRET was confirmed by gel-shift assays.     

On-chip hybridization kinetics for optimization of gene expression experiments

DNA microarray technology is a powerful tool for getting an overview of gene expression in biological samples. Although the successful use of microarray-based expression analysis was demonstrated in a number of applications, the main problem with this approach is the fact that expression levels deduced from hybridization experiments do not necessarily correlate with RNA concentrations. Moreover oligonucleotide probes corresponding to the same gene can give different hybridization signals. Apart from cross-hybridizations and differential splicing, this could be due to secondary structures of probes or targets. In addition, for low-copy genes, hybridization equilibrium may be reached after hybridization times much longer than the one commonly used (overnight, i.e., 15 h). Thus, hybridization signals could depend on kinetic properties of the probe, which may vary between different oligonucleotide probes immobilized on the same microarray. To validate this hypothesis, on-chip hybridization kinetics and duplex thermostability analysis were performed using oligonucleotide microarrays containing 50-mer probes corresponding to 10 mouse genes. We demonstrate that differences in hybridization kinetics between the probes exist and can influence the interpretation of expression data. In addition, we show that using on-chip hybridization kinetics, quantification of targets is feasible using calibration curves.     

小鼠细胞因子相关基因表达检测寡核苷酸芯片的制备及分析

生物芯片技术用于基因表达谱研究是近年来发展起来的一项新技术 ,该方法本质上是基于对一玻璃片或膜表面上固定的ｃＤＮＡ或寡核苷酸的分子杂交 ,这一新技术可同时测定成千上万个基因的作用方式 ,几周获得的信息用其它方法可能要几年才能得到 ,是以定量方式同时监测大量基因相对表达的强有力的新方法[1 ,2 ] 。国内外目前主要采用ｃＤＮＡ芯片进行基因表达的检测 ,芯片制备所用的ＤＮＡ探针一般为已知基因ｃＤＮＡ克隆的ＰＣＲ扩增产物或ＥＳＴ的扩增产物[3～ 8] 。对基因的表达检测来说 ,ｃＤＮＡ芯片技术是一条非常适用的检测方法 ,但在有…     

Analysis of crucial factors resulting in microarray hybridization failure

The factors that affect the formation and stability of DNA/DNA duplexes are complicated and still mostly unknown. In this study attempts were made to look for the crucial factor affecting hybridization failure in DNA microarray assays. A comprehensive range of factors were investigated simultaneously using a 25-mer oligonucleotide Potyvirus microarray. These included steric hindrance, direct/indirect labelling types, distance of a probe to the fluorescent labelling end, target (the DNA fragment used to hybridize with microarray probes) strand types either single strand or double strand, probes without mismatch and with different numbers of mismatch nucleotides (up to 36%) and different mismatch locations (5' end, centre and 3' end), probe GC content and T(m), secondary structures of probes and targets, different target lengths (0.277 kb to ~1.3 kb) and concentrations (0.1-30 nM). The results showed that whilst most of these known factors were unlikely to be the main causes of failed hybridization, there was strong evidence suggesting that the viral amplicon target structure is the most crucial factor. However, computing predicted target secondary structures by Mfold showed no correlation with the hybridization results. One explanation is that the predicted target secondary structures are different from the real structures. Here we postulate that the real target structure might be a combination of secondary structures resulting in a three-dimensional structure from exposure to three types of sub-structures: (1) a completely exposed linear structure to allow probes access for the successful hybridization and showing strong fluorescent signals; (2) a partially exposed structure to allow unstable binding and showing weak fluorescent signals; (3) a closed structure resulting in failed hybridization. These results are very important for microarray based studies as they not only provide an explanation for some current controversial results, but also provide potential resolution for the future studies. Due to the lack of available software for predicting the true target structure, development of microarrays should conduct an initial oligonucleotide probe selection procedure and those probes with capacity to hybridize with the target should be considered for the microarray development.     

Cooperativity of paired oligonucleotide probes for microarray hybridization assays

Synthetic DNA probes attached to microarrays usually range in length from 25 to 70 nucleotides. There is a compromise between short probes with lower sensitivity, which can be accurately synthesized in higher yields, and long probes with greater sensitivity but lower synthesis yields. Described here are microarrays printed with spots containing a mixture of two short probes, each designed to hybridize at noncontiguous sites in the same targeted sequence. We have shown that, for a printed microarray, mixed probe spots containing a pair of 30mers show significantly greater hybridization than spots containing a single 30mer and can approach the amount of hybridization to spots containing a 60mer or a 70mer. These spots with mixed oligonucleotide probes display cooperative hybridization signals greater than those that can be achieved by either probe alone. Both the higher synthesis yields of short probes and the greater sensitivity of long oligonucleotides can be utilized. This strategy provides new design options for microarray hybridization assays to detect RNA abundance, RNA splice variants, or sequence polymorphisms.     

Towards quality control for DNA microarrays.

We present a framework for detecting probes in oligonucleotide microarrays that may add significant error to measurements in hybridization experiments. Four types of so-called degenerate probe behavior are considered: secondary structure formation, self-dimerization, cross-hybridization, and dimerization. The framework uses a well-established model for computing the free energy of nucleic acid sequence hybridization and a novel method for the detection of patterns in hybridization experiment data. Our primary result is the identification of unique patterns in hybridization experiment data that are shown to correlate with each type of degenerate probe behavior. A support function for identifying degenerate probes from a large set of hybridization experiments is given and some preliminary experimental results are given for the Affymetrix HuGeneFL GeneChip. Finally, we show a strong relationship between the Affymetrix discrimination measure for a probe and the free-energy estimate from theoretical models of hybridization. In particular, probes on the HuGeneFL GeneChip with high free-energy estimates (weak hybridization) have almost always approximately zero discrimination. The framework can be applied to any Affymetrix oligonucleotide array, and the software is made freely available to the community.     

Improvement of oligonucleotide probe design criteria for functional gene microarrays in environmental applications

To optimize oligonucleotide probe design criteria, PCR products with different similarities to probes were hybridized to a functional gene microarray designed to detect homologous genes from different organisms. In contrast to more restrictive probe designs based on a single criterion, simultaneous consideration of the percent similarity (< or =90%), the length of identical sequence stretches (< or =20 bases), and the binding free energy (> or =-35 kcal mol(-1)) was found to be predictive of probe specificity.     

Identification of high-stringency DNA hairpin probes by partial gene folding

Hairpin DNA sequences are widely used as probes for oligonucleotides in a broad range of assays, often as "molecular beacons". A potential disadvantage of the standard methodology for molecular beacon design is the need to add several self-complementary bases to each end of the probe, since these do not correspond to the target sequence. We describe a conceptually new method of hairpin DNA probe identification, in which a secondary structure prediction algorithm is employed to identify oligonucleotide sequences within an expressed gene having the requisite hairpin structure. Intuitively, such probes should have significantly improved performance over "traditional" hairpin probes, because they are fully complementary with the target. We present experimental evidence verifying this hypothesis for a series of hairpin probes targeting the pag gene of Bacillus anthracis.     

Probe selection algorithms with applications in the analysis of microbial communities

We propose two efficient heuristics for minimizing the number of oligonucleotide probes needed for analyzing populations of ribosomal RNA gene (rDNA) clones by hybridization experiments on DNA microarrays. Such analyses have applications in the study of microbial communities. Unlike in the classical SBH (sequencing by hybridization) procedure, where multiple probes are on a DNA chip, in our applications we perform a series of experiments, each one consisting of applying a single probe to a DNA microarray containing a large sample of rDNA sequences from the studied population. The overall cost of the analysis is thus roughly proportional to the number of experiments, underscoring the need for minimizing the number of probes. Our algorithms are based on two well-known optimization techniques, i.e. simulated annealing and Lagrangian relaxation, and our preliminary tests demonstrate that both algorithms are able to find satisfactory probe sets for real rDNA data.     

Empirical Evaluation of Oligonucleotide Probe Selection for DNA Microarrays

DNA-based microarrays are increasingly central to biomedical research. Selecting oligonucleotide sequences that will behave consistently across experiments is essential to the design, production and performance of DNA microarrays. Here our aim was to improve on probe design parameters by empirically and systematically evaluating probe performance in a multivariate context. We used experimental data from 19 array CGH hybridizations to assess the probe performance of 385,474 probes tiled in the Duchenne muscular dystrophy (DMD) region of the X chromosome. Our results demonstrate that probe melting temperature, single nucleotide polymorphisms (SNPs), and homocytosine motifs all have a strong effect on probe behavior. These findings, when incorporated into future microarray probe selection algorithms, may improve microarray performance for a wide variety of applications.     

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.     

Oligonucleotide conjugates with minor groove ligands as probes for hybridization microarray chips

A possibility of using oligonucleotide conjugates with minor groove ligands as probes for hybridization microarray chips was studied. The oligonucleotide conjugates contain a hairpin ligand (MGB) composed of two tripyrrolcarboxamide residues with an aminocaproic acid residue as a linker and bound to the oligonucleotide duplex AT tract in a site-specific manner. We used as (5'-3') probes GACAAGAp, GACAAAAp, GACAAGA-MGB, and GACAAAA-MGB. The oligonucleotides labeled with Cy3 cyanine dye, Cy3-ACTAATTTTGTC and Cy3-ACTAATCTTGTC, were used as targets. The maximal MGB effect on the fluorescence level of microarray chip spots, which caused its fourfold increase as compared with the initial unmodified duplex, was observed for the duplex containing only AT pairs in the ligand binding site. The presence of A-C and G-T mutations in the binding site (imperfect duplexes) or a C-G pair (perfect duplex) affects the change in fluorescence level to a considerably lesser degree.     

Improvement of Oligonucleotide Probe Design Criteria for Functional Gene Microarrays in Environmental Applications

To optimize oligonucleotide probe design criteria, PCR products with different similarities to probes were hybridized to a functional gene microarray designed to detect homologous genes from different organisms. In contrast to more restrictive probe designs based on a single criterion, simultaneous consideration of the percent similarity (≤90%), the length of identical sequence stretches (≤20 bases), and the binding free energy (≥−35 kcal mol⁻¹) was found to be predictive of probe specificity.