Multiplex genotyping and allele frequency estimation in pooled DNAs using non-gel capillary electrophoresis

We present a novel approach of single-nucleotide polymorphism (SNP) analysis in which allele-specific oligonucleotide hybridization is followed by non-gel capillary electrophoresis (ASOH-NGCE) in conjunction with laser-induced fluorescence (LIF). This allows rapid multiplex allelotyping and allele frequency estimation. This method, based on site separation of the hybridization duplexes, retains the simplicity and specificity of ASOH and the homogeneous feature of NGCE with poly(N,N-dimethylacrylamide) (PDMA) as a sieving medium. ASOH-NGCE can be applied to multiplex SNP loci genotyping with excellent separation of hybridization mixtures. Average relative standard deviations (RSDs) were low for within-day (1.10%) and between-day (2.41%) reproducibility. Moreover, the allele frequencies in pooled DNAs were accurately determined from peak areas and equilibrium dissociation constants. Our method was highly sensitive in detecting alleles with frequency as low as 1% and in distinguishing allele frequencies differing by 1% between pools. The average value of differences between real and estimated frequencies (accuracy) was only 0.004.     

DNA multiplex hybridization on microarrays and thermodynamic stability in solution: a direct comparison

Hybridization intensities of 30 distinct short duplex DNAs measured on spotted microarrays, were directly compared with thermodynamic stabilities measured in solution. DNA sequences were designed to promote formation of perfect match, or hybrid duplexes containing tandem mismatches. Thermodynamic parameters ΔH°, ΔS° and ΔG° of melting transitions in solution were evaluated directly using differential scanning calorimetry. Quantitative comparison with results from 63 multiplex microarray hybridization experiments provided a linear relationship for perfect match and most mismatch duplexes. Examination of outliers suggests that both duplex length and relative position of tandem mismatches could be important factors contributing to observed deviations from linearity. A detailed comparison of measured thermodynamic parameters with those calculated using the nearest-neighbor model was performed. Analysis revealed the nearest-neighbor model generally predicts mismatch duplexes to be less stable than experimentally observed. Results also show the relative stability of a tandem mismatch is highly dependent on the identity of the flanking Watson–Crick (w/c) base pairs. Thus, specifying the stability contribution of a tandem mismatch requires consideration of the sequence identity of at least four base pair units (tandem mismatch and flanking w/c base pairs). These observations underscore the need for rigorous evaluation of thermodynamic parameters describing tandem mismatch stability.     

Use of hybridization melting kinetics for detecting Phytophthora species using three-dimensional microarrays: demonstration of a novel concept for the differentiation of detection targets

Microarray-based detection is limited by variable and inconsistent hybridization intensities across the diversity of probes used in each array. In this paper, we introduce a novel concept for the differentiation of detection targets using duplex melting kinetics. A microarray assay was developed on a PamChip microarray enabling the differentiation of target Phytophthora species using the melting kinetics of probe-target duplexes. In the majority of cases the hybridization kinetics of target and non-target duplexes differed significantly. Analysis of the melting kinetics of duplexes formed by probes with target and non-target DNA was found to be an effective method for determining specific hybridization and was independent of fluctuations in hybridization signal intensity. This form of analysis was more robust than the traditional approach based on hybridization intensity, and enabled the detection of individual Phytophthora species and mixtures thereof.     

Combined SSCP/duplex analysis by capillary electrophoresis for more efficient mutation detection

SSCP and heteroduplex analysis (HA) continue to be the most popular methods of mutation detection due to their simplicity, high sensitivity and low cost. The advantages of these methods are most clearly visible when large genes, such as BRCA1 and BRCA2, are scanned for scattered unknown mutations and/or when a large number of DNA samples is screened for specific mutations. Here we describe a novel combined SSCP/duplex analysis adapted to the modern capillary electrophoresis (CE) system, which takes advantage of multicolor labeling of DNA fragments and laser-induced fluorescence detection. In developing this method, we first established the optimum conditions for homoduplex and heteroduplex analysis by CE. These were determined based on comprehensive analysis of representative Tamra-500 markers and BRCA1 fragments at different concentrations of sieving polymer and temperatures in the presence or absence of glycerol. The intrinsic features of DNA duplex structures are discussed in detail to explain differences in the migration rates between various types of duplexes. When combined SSCP/duplex analysis was carried out in single conditions, those found to be optimal for analysis of duplexes, all 31 BRCA1 and BRCA2 mutations, polymorphisms and variants tested were detected. It is worth noting that the panel of analyzed sequence variants was enriched in base substitutions, which are usually more difficult to detect. The sensitivity of mutation detection in the SSCP portion alone was 90%, and that in the duplex portion was 81% in the single conditions of electrophoresis. As is also shown here, the proposed combined SSCP/duplex analysis by CE has the potential of being applied to the analysis of pooled genomic DNA samples, and to multiplex analysis of amplicons from different gene fragments. These modifications may further reduce the costs of analysis, making the method attractive for large scale application in SNP scanning and screening.     

Smoothing of the thermal stability of DNA duplexes by using modified nucleosides and chaotropic agents

The effect of alkyltrimethylammonium ions on the thermostability of natural and modified DNA duplexes has been investigated. We have shown that the use of tetramethylammonium ions TMA+along with the chemical modification of duplexes allow the fine adjustment of T m and the possibility of obtaining several duplex systems with varied isostabilizedtemperatures, some of which show greater stability than those of natural DNA. This approach could be very useful for DNA sequencing by hybridization.     

Crystal structure of d(gcGXGAgc) with X=G: a mutation at X is possible to occur in a base-intercalated duplex for multiplex formation

DNA fragments with the sequences d(gcGX[Y]n Agc) (n=1, X=A, and Y=A, T, or G)form base-intercalated duplexes, which is a basic unit for formation of multiplexes such as octaplex and hexaplex. To examine the stability of multiplexes, a DNA with X=Y=G and n=1 was crystallized under conditions different from those of the previously determined sequences, and its crystal structure has been determined. The two strands are coupled in an anti-parallel fashion to form a base-intercalated duplex, in which the first and second residues form Watson-Crick type G:C pairs and the third and sixth residues form a sheared G:A pairs at both ends of the duplex. The G4 and G5 bases are stacked alternatively on those of the counter strand to form a long G column of G3-G4-G5*-G5-G4*-G3*, the central four Gs being protruded. In addition, the three duplexes are associated to form a hexaplex around a mixture of calcium and sodium cations on the crystallographic threefold axis. These structural features are similar to those of the previous crystals, though slightly different in detail. The present study indicates that mutation at the 4th position is possible to occur in a base-intercalated duplex for multiplex formations, suggesting that DNA fragments with any sequence sandwiched between the two triplets gcG and Agc can form a multiplex.     

Parallel thermodynamic analysis of duplexes on oligodeoxyribonucleotide microchips.

A microchip method has been developed for massive and parallel thermodynamic analyses of DNA duplexes. Fluorescently labeled oligonucleotides were hybridized with oligonucleotides immobilized in the 100 x 100 x 20 mum gel pads of the microchips. The equilibrium melting curves for all microchip duplexes were measured in real time in parallel for all microchip duplexes. Thermodynamic data for perfect and mismatched duplexes that were obtained using the microchip method directly correlated with data obtained in solution. Fluorescent labels or longer linkers between the gel and the oligonucleotides appeared to have no significant effect on duplex stability. Extending the immobilized oligonucleotides with a four-base mixture from the 3'-end or one or two universal bases (5-nitroindole) from the 3'- and/or 5'-end increased the stabilities of their duplexes. These extensions were applied to increase the stabilities of the duplexes formed with short oligonucleotides in microchips, to significantly lessen the differences in melting curves of the AT- and GC-rich duplexes, and to improve discrimination of perfect duplexes from those containing poorly recognized terminal mismatches. This study explored a way to increase the efficiency of sequencing by hybridization on oligonucleotide microchips.     

Analysis of Perfect and Mismatched DNA Duplexes with a Generic Hexanucleotide Microchip

Thermodynamic analysis was performed for the duplexes formed by fluorescently labeled oligonucleotide targets on a generic hexanucleotide microchip. All 4096 different hexanucleotide chains were immobilized as probes in individual gel pads of the microchip. To strengthen the hybridization, each hexamer was extended at both ends by one nucleotide from the equimolar mixture of all four nucleotides to serve as nonselective linkers. It has been shown that the melting curves for oligonucleotide duplexes formed on the microchip and in a solution are quite similar. The influence of ionic surrounding has been studied in terms of the hybridization efficiency and discrimination between the mismatched and perfect duplexes. Different approaches have been tested to compensate the dependence of duplex stability on the GC content. It has been demonstrated that the use of chaotropic agents, addition of nonlabeled GC-rich competitor oligonucleotides, as well as creation of a temperature gradient along the microchip reproducing the distribution of melting temperatures, efficiently level out the AT/GC differences. The use of tetramethylammonium chloride for the same purpose was accompanied by weakening to some extent the discrimination between the mismatched duplexes and the perfect ones.     

Analysis of perfect and mismatched DNA duplexes by a generic hexanucleotide microchip

The thermodynamic analysis was done for the duplexes formed by fluorescently labeled oligonucleotide targets on a genetic hexanucleotide microchip. All 4096 different hexanucleotide chains were immobilized as probes in individual gel pads of the microchip. To strengthen the hybridization, each hexamer was extended at both ends by one nucleotide from the equimolar mixture of all four nucleotides to serve as nonselective linkers. It has been shown that the melting curves for oligonucleotide duplexes formed on the microchip and in a solution are quite similar. The influence of ionic surrounding has been studied in terms of the hybridization efficiency and discrimination between the mismatched and perfect duplexes. Different approaches have been tested to compensate the dependence of duplex stability on the GC content. It has been demonstrated that the use of chaotropic agents, addition of nonlabeled GC-rich competitor oligonucleotides, as well as creation of a temperature gradient along the microchip reproducing the distribution of melting temperatures, efficiently level out the AT/GC differences. The use of tetramethylammonium chloride for the same purpose was accompanied by weakening to some extent the discrimination between the mismatched duplexes and the perfect ones.     

A method for HLA genotyping using the specific cleavage of DNA-rN1-DNA/DNA with RNase HII from Chlamydia pneumoniae

Single nucleotide polymorphisms (SNPs) provide a great opportunity for the study of human disease and bacterial drug resistance. However, many SNP typing techniques require dedicated instruments and high cost. Here, we develop a novel method for SNP genotyping based on specific cleavage properties of RNase HII from Chlamydia pneumoniae (CpRNase HII), termed the "CpRNase HII-based method." CpRNase HII cleaves the DNA-rN(1)-DNA/DNA duplex at the 5'-side of the ribonucleotide (rN(1) = one ribonucleotide). Moreover, the cleavage efficiencies of the perfectly matched DNA-rN(1)-DNA/DNA duplexes are higher than those carrying a mismatched ribonucleotide. DNA-rN(1)-DNA fragments are modified with a fluorophore at the 5'-end and a quencher at the 3'-end to generate molecular beacons (MBs), which hybridize with single-stranded DNA (analyte) to be cleaved by CpRNase HII. As perfectly matched duplexes can be cleaved efficiently and mismatched duplexes cannot, CpRNase HII-catalyzed reactions can differentiate between one-nucleotide variations on the DNA-rN(1)-DNA/DNA duplexes. We have validated this method with nine SNPs of the HLA gene, which were successfully determined by endpoint measurements of fluorescence intensity. The new method is simple and effective, because the design of MBs is easy, and all steps of the genotyping consist of simple additions of solutions and incubation. This method will be suitable for large-scale genotyping.     

Effect of thiazole orange doubly labeled thymidine on DNA duplex formation

Nucleic acid oligonucleotides are widely used in hybridization experiments for specific detection of complementary nucleic acid sequences. For design and application of oligonucleotides, an understanding of their thermodynamic properties is essential. Recently, exciton-controlled hybridization-sensitive fluorescent oligonucleotides (ECHOs) were developed as uniquely labeled DNA oligomers containing commonly one thymidine having two covalently linked thiazole orange dye moieties. The fluorescent signal of an ECHO is strictly hybridization-controlled, where the dye moieties have to intercalate into double-stranded DNA for signal generation. Here we analyzed the hybridization thermodynamics of ECHO/DNA duplexes, and thermodynamic parameters were obtained from melting curves of 64 ECHO/DNA duplexes measured by ultraviolet absorbance and fluorescence. Both methods demonstrated a substantial increase in duplex stability (ΔΔG°(37) ～ -2.6 ± 0.7 kcal mol(-1)) compared to that of DNA/DNA duplexes of the same sequence. With the exception of T·G mismatches, this increased stability was mostly unaffected by other mismatches in the position opposite the labeled nucleotide. A nearest neighbor model was constructed for predicting thermodynamic parameters for duplex stability. Evaluation of the nearest neighbor parameters by cross validation tests showed higher predictive reliability for the fluorescence-based than the absorbance-based parameters. Using our experimental data, a tool for predicting the thermodynamics of formation of ECHO/DNA duplexes was developed that is freely available at http://genome.gsc.riken.jp/echo/thermodynamics/ . It provides reliable thermodynamic data for using the unique features of ECHOs in fluorescence-based experiments.     

Free energy of DNA duplex formation on short oligonucleotide microarrays

DNA/DNA duplex formation is the basic mechanism that is used in genome tiling arrays and SNP arrays manufactured by Affymetrix. However, detailed knowledge of the physical process is still lacking. In this study, we show a free energy analysis of DNA/DNA duplex formation these arrays based on the positional-dependent nearest-neighbor (PDNN) model, which was developed previously for describing DNA/RNA duplex formation on expression microarrays. Our results showed that the two ends of a probe contribute less to the stability of the duplexes and that there is a microarray surface effect on binding affinities. We also showed that free energy cost of a single mismatch depends on the bases adjacent to the mismatch site and obtained a comprehensive table of the cost of a single mismatch under all possible combination of adjacent bases. The mismatch costs were found to be correlated with those determined in aqueous solution. We further demonstrate that the DNA copy number estimated from the SNP array correlates negatively with the target length; this is presumably caused by inefficient PCR amplification for long fragments. These results provide important insights into the molecular mechanisms of microarray technology and have implications for microarray design and the interpretation of observed data.     

DNA sequence analysis by hybridization with oligonucleotide microchips: MALDI mass spectrometry identification of 5mers contiguously stacked to microchip oligonucleotides

Matrix-assisted laser desorption ionization mass spectrometry (MALDI MS) has been applied to increase the informational output from DNA sequence analysis. It has been used to analyze DNA by hybridization with microarrays of gel-immobilized oligonucleotides extended with stacked 5mers. In model experiments, a 28 nt long DNA fragment was hybridized with 10 immobilized, overlapping 8mers. Then, in a second round of hybridization DNA–8mer duplexes were hybridized with a mixture of 10 5mers. The stability of the 5mer complex with DNA was increased to raise the melting temperature of the duplex by 10–15°C as a result of stacking interaction with 8mers. Contiguous 13 bp duplexes containing an internal break were formed. MALDI MS identified one or, in some cases, two 5mers contiguously stacked to each DNA–8mer duplex formed on the microchip. Incorporating a mass label into 5mers optimized MALDI MS monitoring. This procedure enabled us to reconstitute the sequence of a model DNA fragment and identify polymorphic nucleotides. The application of MALDI MS identification of contiguously stacked 5mers to increase the length of DNA for sequence analysis is discussed.     

Thermodynamic study of hybridization properties of heterochiral nucleic acids

Heterochiral DNA and RNA heptamers, which contained an unnatural L-nucleotide, were synthesized, and thermodynamic analyses of their hybridization properties with complementary DNA and RNA strands were systematically conducted by UV melting experiments. The results clearly demonstrated that the incorporation of an L-ribonucleotide into the RNA strand leads to more significant destabilization of the duplexes than that of an L-deoxyribonucleotide into the DNA strand, regardless of whether the complementary strand is DNA or RNA. The destabilization of the duplexes by the substitution of D-thymidine with L-thymidine in the DNA strand is entropically driven, whereas that by the substitution of D-uridine with L-uridine in the RNA strand is enthalpically driven. The thermodynamic characteristic that the stability of homochiral duplex is far superior to that of heterochiral duplex is much more remarkable in RNA than in DNA. Thus, RNA might have been a self-replicating system superior to DNA to exclude the chiral antipode.     

Thermodynamic, counterion, and hydration effects for the incorporation of locked nucleic acid nucleotides into DNA duplexes

A locked nucleic acid (LNA) monomer is a conformationally restricted nucleotide analogue with an extra 2'-O, 4'-C-methylene bridge added to the ribose ring. LNA-modified oligonucleotides are known to exhibit enhanced hybridization affinity toward complementary DNA and RNA. In this work, we have evaluated the hybridization thermodynamics of a series of LNA-substituted DNA octamers, modified to various extents by one to three LNA substitutions, introduced at either adenine (5'-AGCACCAG) or thymine (5'-TGCTCCTG) nucleotides. To understand the energetics, counterion effects, and the hydration contribution of the incorporation of LNA modification, a combination of spectroscopic and calorimetric techniques was used. The CD spectra of the corresponding duplexes showed that the modified duplexes adopt an A-type conformation. UV and DSC melting studies revealed that each type of duplex unfolds in a two-state transition. A complete thermodynamic profile at 5 degrees C indicated that the net effect of modification on thermodynamic parameters might be positional and that the neighboring bases flanking the modification might influence the favorable formation of the modified duplexes. Furthermore, relative to the formation of the unmodified reference duplexes, the formation of modified duplexes is accompanied by a higher uptake of counterions and a lower uptake of water molecules.     

Simultaneous detection and identification of trichothecene- and moniliformin-producing Fusarium species based on multiplex SNP analysis

AIM: To develop a multiplex identification method for trichothecene- and moniliformin-producing Fusarium species. METHOD AND RESULTS: In this article, we present a single nucleotide polymorphism (SNP) assay to simultaneously detect and identify 16 trichothecene- and moniliformin-producing Fusarium species. A number of SNP primers are designed to detect clades of species with particular mycotoxigenic synthetic abilities. The assay is based on minisequencing using SNaPshot reactions and the SNP primers are designed based on motifs derived from phylogenetic analyses of translation elongation factor-1alpha sequences. The present version of the Fusarium SNP assay can distinguish major groups of trichothecene producers; the strict-type-A, the strict-type-B, the type-A and type-B trichothecene producers and the putative moniliformin producers. The SNP assay was validated against five naturally infected cereal samples that previously had been analysed morphologically, chemically and by a multiplex DNA array hybridization. CONCLUSIONS: The Fusarium SNP assay reveals the advantages of using SNPs for multiplex species identification. Significance AND IMPACT OF THE STUDY: The current assay may qualify as a high-throughput screening method for small-grain cereals in the feed and food chain, and may facilitate detection of new or introduced Fusarium species.     

The thermodynamic advantage of DNA oligonucleotide 'stacking hybridization' reactions: energetics of a DNA nick.

'Stacking hybridization reactions' wherein two or more short DNA oligomers hybridize in a contiguous tandem orientation onto a longer complementary DNA single strand have been employed to enhance a variety of analytical oligonucleotide hybridization schemes. If the short oligomers anneal in perfect head-to-tail register the resulting duplex contains a nick at every boundary between hybridized oligomers. Alternatively, if the short oligomers do not hybridize precisely in register, i.e. single strand regions on the longer strand are left unbound, gaps are formed between regions where short oligomers bind. The resulting gapped DNA duplexes are considerably less stable than their nicked duplex analogs. Formation of base pair stacking interactions between neighboring oligomers at the nicks that do not occur in gapped duplexes has been proposed as the source of the observed added stability. However, quantitative evidence supporting this hypothesis for DNA has not been reported. Until now, a direct comparison of the thermodynamics of DNA nicks versus DNA gaps has not been performed. In this communication we report such a comparison. Analysis of optical melting experiments in a well defined molecular context enabled quantitative evaluations of the relative thermodynamic difference between nicked and gapped DNA duplexes. Results of the analysis reveal that a nick may be energetically favored over a gap by at least 1.4 kcal/mol and perhaps as much as 2.4 kcal/mol. The presence of a 5'phosphate at a nick or gap fails to significantly affect their stabilities.     

A measurement of the fraction of chloroplast DNA transcribed in Euglena

The fraction of the chloroplast DNA transcribed in the single celled alga $\text{Euglena}$ has been determined by RNA-DNA hybridization. A vast excess of total cell RNA from cells which were rapidly dividing in the light was hybridized in liquid to [¹²⁵I] — chloroplast DNA, and the resulting duplexes separated on hydroxyapatite columns. The contribution of DNA-DNA duplex formation was determined separately and was used to calculate that portion of the duplex which was actually a RNA-DNA hybrid. Sixteen percent of the single stranded chloroplast DNA forms a duplex with this RNA suggesting that 32 percent of the double stranded DNA molecule is being transcribed into RNA under these conditions of cell growth.     

Cooperative lectin recognition of periodical glycoclusters along DNA duplexes: alternate hybridization and full hybridization

We describe herein the construction of periodically, spatially controlled glycoclusters along DNA duplexes and their cooperative lectin recognition. Site-specifically alpha-mannosylated oligodeoxynucleotide 20-mer (Man-ODN20) was synthesized via the phosphoramidite solid-phase synthesis. Alternate hybridization of the Man-ODN20 with the half-sliding complementary ODN 20-mer (hscODN20) gave an alternately prolonged Man-cluster Man-ODN20/hscODN20. The binding of the Man-cluster to FITC-labeled ConA lectin showed sigmoidal fluorescence dependency on the concentration of Man-ODN, indicating that some mannose residues along the repeating DNA duplex were cooperatively bound to ConA (apparent affinity constant: K(af)=2.4 x 10(4)M(-1) and Hill coefficient: n=3.5). The duplex of Man-ODN20 with full complementary ODN 20-mer (fcODN20) was little bound to ConA. The binding behavior of Man-ODN20/hscODN20 is compared with that of the alternately prolonged Gal-cluster Gal-ODN20/hscODN20 previously reported. Duplexes 20-mer, 40-mer, and 60-mer presenting one, two, and three periodic galactoses were also prepared by full hybridization of 20-mer beta-galactosylated oligodeoxynucleotide (Gal-ODN20) with the periodically repeating full complementary 20-mer, 40-mer, and 60-mer ODNs. RCA(120) lectin was found to little bind the 20-mer and 40-mer duplexes and to bind weakly and non-cooperatively the 60-mer duplex (K(af)=1.1 x 10(4)M(-1)). The cooperative lectin recognition of these glycoclusters in relation with the degree of association (DA) of ODN and the numbers of glycosides along the DNA duplex is discussed.