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.     

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.     

A powerful global test for spliceQTL effects

Statistical methods to test for effects of single nucleotide polymorphisms (SNPs) on exon inclusion exist but often rely on testing of associations between multiple exon–SNP pairs, with sometimes subsequent summarization of results at the gene level. Such approaches require heavy multiple testing corrections and detect mostly events with large effect sizes. We propose here a test to find spliceQTL (splicing quantitative trait loci) effects that takes all exons and all SNPs into account simultaneously. For any chosen gene, this score-based test looks for an association between the set of exon expressions and the set of SNPs, via a random-effects model framework. It is efficient to compute and can be used if the number of SNPs is larger than the number of samples. In addition, the test is powerful in detecting effects that are relatively small for individual exon–SNP pairs but are observed for many pairs. Furthermore, test results are more often replicated across datasets than pairwise testing results. This makes our test more robust to exon–SNP pair-specific effects, which do not extend to multiple pairs within the same gene. We conclude that the test we propose here offers more power and better replicability in the search for spliceQTL effects.     

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.     

Hybridization thermodynamics of NimbleGen Microarrays

Background  

While microarrays are the predominant method for gene expression profiling, probe signal variation is still an area of active research. Probe signal is sequence dependent and affected by probe-target binding strength and the competing formation of probe-probe dimers and secondary structures in probes and targets.     

Analyzing Illumina Gene Expression Microarray Data from Different Tissues: Methodological Aspects of Data Analysis in the MetaXpress Consortium

Microarray profiling of gene expression is widely applied in molecular biology and functional genomics. Experimental and technical variations make meta-analysis of different studies challenging. In a total of 3358 samples, all from German population-based cohorts, we investigated the effect of data preprocessing and the variability due to sample processing in whole blood cell and blood monocyte gene expression data, measured on the Illumina HumanHT-12 v3 BeadChip array.Gene expression signal intensities were similar after applying the log₂ or the variance-stabilizing transformation. In all cohorts, the first principal component (PC) explained more than 95% of the total variation. Technical factors substantially influenced signal intensity values, especially the Illumina chip assignment (33–48% of the variance), the RNA amplification batch (12–24%), the RNA isolation batch (16%), and the sample storage time, in particular the time between blood donation and RNA isolation for the whole blood cell samples (2–3%), and the time between RNA isolation and amplification for the monocyte samples (2%). White blood cell composition parameters were the strongest biological factors influencing the expression signal intensities in the whole blood cell samples (3%), followed by sex (1–2%) in both sample types. Known single nucleotide polymorphisms (SNPs) were located in 38% of the analyzed probe sequences and 4% of them included common SNPs (minor allele frequency >5%). Out of the tested SNPs, 1.4% significantly modified the probe-specific expression signals (Bonferroni corrected p-value<0.05), but in almost half of these events the signal intensities were even increased despite the occurrence of the mismatch. Thus, the vast majority of SNPs within probes had no significant effect on hybridization efficiency.In summary, adjustment for a few selected technical factors greatly improved reliability of gene expression analyses. Such adjustments are particularly required for meta-analyses.     

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.     

Influence of Dangling Ends and Surface-Proximal Tails of Targets on Probe-Target Duplex Formation in 16S rRNA Gene-Based Diagnostic Arrays

Dangling ends and surface-proximal tails of gene targets influence probe-target duplex formation and affect the signal intensity of probes on diagnostic microarrays. This phenomenon was evaluated using an oligonucleotide microarray containing 18-mer probes corresponding to the 16S rRNA genes of 10 waterborne pathogens and a number of synthetic and PCR-amplified gene targets. Signal intensities for Klenow/random primer-labeled 16S rRNA gene targets were dissimilar from those for 45-mer synthetic targets for nearly 73% of the probes tested. Klenow/random primer-labeled targets resulted in an interaction with a complex mixture of 16S rRNA genes (used as the background) 3.7 times higher than the interaction of 45-mer targets with the same mixture. A 7-base-long dangling end sequence with perfect homology to another single-stranded background DNA sequence was sufficient to produce a cross-hybridization signal that was as strong as the signal obtained by the probe-target duplex itself. Gibbs free energy between the target and a well-defined background was found to be a better indicator of hybridization signal intensity than the sequence or length of the dangling end alone. The dangling end (Gibbs free energy of −7.6 kcal/mol) was found to be significantly more prone to target-background interaction than the surface-proximal tail (Gibbs free energy of −64.5 kcal/mol). This study underlines the need for careful target preparation and evaluation of signal intensities for diagnostic arrays using 16S rRNA and other gene targets due to the potential for target interaction with a complex background.     

Peptide nucleic acid (PNA) facilitates multistranded hybrid formation between linear double-stranded DNA targets and RecA protein-coated complementary single-stranded DNA probes

RecA protein-coated single-stranded DNA probes, known as RecA nucleoprotein filaments, bind specifically to homologous DNA sequences within double-stranded DNA targets, forming multistranded probe-target DNA hybrids. This DNA hybridization reaction can be used for sequence-specific gene capture, gene modification, and gene regulation. Thus, factors that enhance the efficiency of the hybridization reaction are of significant practical importance. We show here that the hybridization of a peptide nucleic acid (PNA) within or adjacent to the probe-target homology region significantly enhances the yield of hybrid DNA formed in the reaction between linear double-stranded DNA targets and RecA protein-coated complementary single-stranded (css)DNA probes. The possible mechanisms and the advantages of using RecA nucleoprotein filaments in combination with PNA for genomic DNA cloning and mutagenesis are presented.     

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.     

Influence of dangling ends and surface-proximal tails of targets on probe-target duplex formation in 16S rRNA gene-based diagnostic arrays

Dangling ends and surface-proximal tails of gene targets influence probe-target duplex formation and affect the signal intensity of probes on diagnostic microarrays. This phenomenon was evaluated using an oligonucleotide microarray containing 18-mer probes corresponding to the 16S rRNA genes of 10 waterborne pathogens and a number of synthetic and PCR-amplified gene targets. Signal intensities for Klenow/random primer-labeled 16S rRNA gene targets were dissimilar from those for 45-mer synthetic targets for nearly 73% of the probes tested. Klenow/random primer-labeled targets resulted in an interaction with a complex mixture of 16S rRNA genes (used as the background) 3.7 times higher than the interaction of 45-mer targets with the same mixture. A 7-base-long dangling end sequence with perfect homology to another single-stranded background DNA sequence was sufficient to produce a cross-hybridization signal that was as strong as the signal obtained by the probe-target duplex itself. Gibbs free energy between the target and a well-defined background was found to be a better indicator of hybridization signal intensity than the sequence or length of the dangling end alone. The dangling end (Gibbs free energy of -7.6 kcal/mol) was found to be significantly more prone to target-background interaction than the surface-proximal tail (Gibbs free energy of -64.5 kcal/mol). This study underlines the need for careful target preparation and evaluation of signal intensities for diagnostic arrays using 16S rRNA and other gene targets due to the potential for target interaction with a complex background.     

AffyMAPSDetector: a software tool to characterize Affymetrix GeneChip™ expression arrays with respect to SNPs

Background  

Affymetrix gene expression arrays incorporate paired perfect match (PM) and mismatch (MM) probes to distinguish true signals from those arising from cross-hybridization events. A MM signal often shows greater intensity than a PM signal; we propose that one underlying cause is the presence of allelic variants arising from single nucleotide polymorphisms (SNPs). To annotate and characterize SNP contributions to anomalous probe binding behavior we have developed a software tool called AffyMAPSDetector.     

Normalized Affymetrix expression data are biased by G-quadruplex formation

Probes with runs of four or more guanines (G-stacks) in their sequences can exhibit a level of hybridization that is unrelated to the expression levels of the mRNA that they are intended to measure. This is most likely caused by the formation of G-quadruplexes, where inter-probe guanines form Hoogsteen hydrogen bonds, which probes with G-stacks are capable of forming. We demonstrate that for a specific microarray data set using the Human HG_U133A Affymetrix GeneChip and RMA normalization there is significant bias in the expression levels, the fold change and the correlations between expression levels. These effects grow more pronounced as the number of G-stack probes in a probe set increases. Approximately 14% of the probe sets are directly affected. The analysis was repeated for a number of other normalization pipelines and two, FARMS and PLIER, minimized the bias to some extent. We estimate that ～15% of the data sets deposited in the GEO database are susceptible to the effect. The inclusion of G-stack probes in the affected data sets can bias key parameters used in the selection and clustering of genes. The elimination of these probes from any analysis in such affected data sets outweighs the increase of noise in the signal.     

Effect of secondary structure on single nucleotide polymorphism detection with a porous microarray matrix; implications for probe selection

Oligonucleotide arrays capable of detecting single nucleotide polymorphisms (SNPs) from amplified nucleic acid have many applications. The expected SNP is usually placed approximately in the center of the probe to ensure the maximum shift in Tm between complementary and SNP sequences. Unfortunately, different short probes (< 30 bases) selected using widely accepted criteria do not perform consistently in this type of assay. Here we present a systematic study on the effect of secondary structure on the ability of oligonucleotide probes to detect an SNP, using real-time array monitoring of a porous microarray substrate that incorporates a novel intra-array mixing system. These results demonstrate that, although positioning of an SNP in the middle of the probe is highly destabilizing, the effect of stable secondary structure on the signal obtained is so dramatic that such probes may be very insensitive. Therefore, if the SNP flanking sequence contains significant secondary structure, then more sensitive probes with good specificity may be obtained by positioning the mutation towards one end of the probe.     

Single nucleotide polymorphisms affect both cis- and trans-eQTLs

Single nucleotide polymorphisms (SNPs) between microarray probes and RNA targets can affect the performance of expression array by weakening the hybridization. In this paper, we examined the effect of the SNPs on Affymetrix GeneChip probe set summaries and the expression quantitative trait loci (eQTL) mapping results in two eQTL datasets, one from mouse and one from human. We showed that removing SNP-containing probes significantly changed the probe set summaries and the more SNP-containing probes we removed the greater the change. Comparison of the eQTL mapping results between with and without SNP-containing probes showed that less than 70% of the significant eQTL peaks were concordant regardless of the significance threshold. These results indicate that SNPs do affect both probe set summaries and eQTLs (both cis and trans), thus SNP-containing probes should be filtered out to improve the performance of eQTL mapping.     

DNA suspension arrays: silencing discrete artifacts for high-sensitivity applications

Detection of low frequency single nucleotide polymorphisms (SNPs) has important implications in early screening for tumorgenesis, genetic disorders and pathogen drug resistance. Nucleic acid arrays are a powerful tool for genome-scale SNP analysis, but detection of low-frequency SNPs in a mixed population on an array is problematic. We demonstrate a model assay for HIV-1 drug resistance mutations, wherein ligase discrimination products are collected on a suspension array. In developing this system, we discovered that signal from multiple polymorphisms was obscured by two discrete hybridization artifacts. Specifically: 1) tethering of unligated probes on the template DNA elicited false signal and 2) unpredictable probe secondary structures impaired probe capture and suppressed legitimate signal from the array. Two sets of oligonucleotides were used to disrupt these structures; one to displace unligated reporter labels from the bead-bound species and another to occupy sequences which interfered with array hybridization. This artifact silencing system resulted in a mean 21-fold increased sensitivity for 29 minority variants of 17 codons in our model assay for mutations most commonly associated with HIV-1 drug resistance. Furthermore, since the artifacts we characterized are not unique to our system, their specific inhibition might improve the quality of data from solid-state microarrays as well as from the growing number of multiple analyte suspension arrays relying on sequence-specific nucleic acid target capture.