Quantitative Heteroduplex Analysis

General DNA analysis includes detection of targets, as well as identification and quantification of specific variants. Genomic DNA exists in the form of homoduplexes, with all corresponding base pairs being complementary, A:T and C:G. We call a double-stranded DNA (dsDNA) molecule a heteroduplex when it contains any noncomplementary base pairs. Conformational and thermodynamic changes produced by mismatched bases facilitate the detection of heteroduplexes. Denaturation of samples followed by hybridization to promote heteroduplex formation has been used to screen diploid DNA for heterozygous variations. This form of heteroduplex analysis has focused on detection, and although it may involve a quantitative threshold, the result is essentially binary. Recent work in this area (1 ) and elsewhere (2 ) has explored new implementations and applications for a more truly quantitative form of heteroduplex analysis. One setting in which heteroduplexes commonly arise is during PCR. Although the extension process of PCR replicates genomic homoduplexes with high accuracy, if the original sample is heterozygous, heteroduplexes arise during the plateau stage when dsDNA is more likely to be formed by hybridization than extension. Denaturation and hybridization associates complementary and nearcomplementary strands almost randomly. Any of various methods (listed below) that can detect the presence of heteroduplexes can then be used to detect heterozygous sequence variations. If heteroduplexes are not detected, variations, if any, are known to be homozygous. Homozygous variants may sometimes be distinguished by differences in melting temperature, but these differences may be near the limit of detection even with high-resolution techniques. Indeed, 4% of human single nucleotide polymorphisms (SNPs) exhibit a symmetry that makes homozygous variants thermodynamically identical up to the nearest-neighbor approximation (3 ). In this case, a 2-stage mixing strategy allows heteroduplex analysis to be used instead of melting temperature analysis for the detection of homozygous changes. Samples that are identified as not heterozygous in the first analysis are mixed with known wild-type DNA, denatured, and reannealed. Detection of the heteroduplexes in the mixed sample then indicates the presence of a distinct homozygous variant of the wild-type DNA. The 2 analyses can be performed in parallel, but then all samples must be prepared mixed and unmixed, then amplified, denatured, hybridized, and analyzed. Alternatively, the mixing can be performed after determining which samples are not heterozygous. The specific methods that have been employed for detection of heteroduplexes include postPCR separation techniques such as conventional gel electrophoresis (4 –9 ), electron microscopy and electrophoresis (10 ), denaturing HPLC (dHPLC) (11–13 ), and temperaturegradient capillary electrophoresis (TGCE) (14 ). Closedtube methods without separation have also been used to detect heteroduplexes by changes in the shape of high-resolution melting curves with fluorescently labeled primers (15 ) or a saturating DNA dye (16 ). Direct sequencing has also been used to identify and quantify heteroduplexes in gene dosage studies (17 ). Other studies (18 –20 ) have analyzed the effects of heteroduplex formation on molecular diagnostic tests including quantitative PCR. In quantitative heteroduplex analysis, the heteroduplex content of a product is indicated by the magnitude of the changes detected by these methods after PCR. Then the proportions of dsDNA genotypes in an original sample before PCR can be determined. This analysis involves 2 primary factors, isolating the effect of the heteroduplex content on the measured quantity (such as the melting curve or elution trace), and establishing the relationship of the measured heteroduplex content to the duplex proportions in the original sample. The methods for systematically isolating and quantifying the portion of the signal that depends on heteroduplex content are somewhat platform specific. These methods typically involve separation of a background component, followed by identification of measurable quantities such as heights or areas of particular peaks in the normalized data that are directly proportional to heteroduplex content. For dHPLC, it has been shown (1 ) that temperature conditions can be optimized by use of predictive software to improve the separation of heteroduplex and homoduplex peaks and achieve high-resolution results. For high-resolution melting, the procedure has been carried out (2 ) using a novel exponential fit to 2 slopes of the raw curve, at temperatures below and above its transition, and reported elsewhere (21 ). As an independent confirmation, quantitative heteroduplex analysis was performed on TGCE traces taken from the same PCR products in (2). The reports by Lim et al. (1 ) and Palais et al. (2 ) both confirm that the mathematical model relating the proportions of 2 different types of homozygous DNA in a sample before PCR to the heteroduplex content after PCR is simple yet accurate. This model describes the heteroduplex proportion after PCR is performed on a mixture that before amplification consists of 2 types of dsDNA, type A:A in proportion x and type B:B in proportion 1 x. These reports support the assumption that different strands are amplified in proportion to their initial concentrations, a theory that has been confirmed in multiple studies (5–7 ). Also, thermodynamic analyses predict that for small mutations, the strands A and B are sufficiently similar that they hybridize randomly with their own and each other’s complements, again in proportion to their concentrations. Based on these considerations, after PCR, the total heteroduplex proportion will be the sum of the Editorial