Comprehensive detection of single base changes in human genomic DNA using denaturing gradient gel electrophoresis and a GC clamp

We present a simple, efficient extension of denaturing gradient gel electrophoresis that allows the detection of nearly any sequence change in a defined fragment of DNA. The fragment can be obtained either by means of the polymerase chain reaction or by restriction digestion of genomic DNA. With restriction fragments of genomic DNA, sequence information is not required, and covalent modifications in genomic DNA that are lost in a PCR, such as methylation, are detectable. We describe how a GC clamp (an arbitrary, G+C-rich sequence of 30 to 60 bp) can be attached to a selected restriction fragment present in a digest of genomic DNA. The GC clamp alters the melting properties of the fragment; this change greatly increases the fraction of possible mutations that is detectable. In a 272-bp HaeIII fragment from the human beta-globin gene, we were able to detect 13 of 13 mutations tested in human genomic DNA. Four additional mutations in cloned plasmids were analyzed. The data agree with a simple theoretical model for DGGE, which predicts how two fragments, differing at a single (specified) base pair, are resolved in a gradient gel as a function of running time for the gel. The calculation assists in the design of probes and gel conditions that aid in the detection of sequence changes.     

Analysis of sequence alterations in a defined DNA region: comparison of temperature-modulated heteroduplex analysis and denaturing gradient gel electrophoresis

The ability to detect DNA sequence heterogeneity quickly and reliably is becoming increasingly important as more genes involved in disease processes are discovered. We have assessed the ability of a high pressure liquid chromatography technique (HPLC) termed temperature-modulated heteroduplex analysis (TMHA) to detect a collection of 20 point mutations distributed throughout a 279 base pair fragment spanning the exon 8 region of the human HPRT gene. All mutant/wild type heteroduplexes formed from mutations in the lowest temperature melting domain of the fragment were easily resolved from the corresponding mutant and wt homoduplexes, while those generated from mutants in the next higher melting domain barely resolved from their parental homoduplexes. For comparison, identical heteroduplex samples were subjected to denaturing gradient gel electrophoresis (DGGE). Heteroduplexes in the lowest temperature melting domain were easily resolved, while no resolution was achieved with those in the next higher melting domain. These results suggest that TMHA and DGGE are measuring similar melting characteristics in heteroduplex molecules. TMHA appears to be a robust approach for detecting and/or purifying a wide variety of mutations in a defined region of DNA, provided that the melting characteristics of the fragment under study are carefully considered.     

Modification of the melting properties of duplex DNA by attachment of a GC-rich DNA sequence as determined by denaturing gradient gel electrophoresis.

The melting behavior of a DNA fragment carrying the mouse beta maj-globin promoter was investigated as a means of establishing procedures for separating DNA fragments differing by any single base substitution using the denaturing gradient gel electrophoresis procedure of Fischer and Lerman (1,2). We find that attachment of a 300 base pair GC-rich DNA sequence, termed a GC-clamp, to a 135 bp DNA fragment carrying the mouse beta-globin promoter significantly alters the pattern of DNA melting within the promoter. When the promoter is attached to the clamp, the promoter sequences melt without undergoing strand dissociation. The calculated distribution of melting domains within the promoter differs markedly according to the relative orientation of the clamp and promoter sequences. We find that the behavior of DNA fragments containing the promoter and clamp sequences on denaturing gradient polyacrylamide gels is in close agreement with the theoretical melting calculations. These studies provide the basis for critical evaluation of the parameters for DNA melting calculations, and they establish conditions for determining whether all single base substitutions within the promoter can be separated on denaturing gradient gels.     

Alterations in DNA helix stability due to base modifications can be evaluated using denaturing gradient gel electrophoresis

DNA molecules that differ by a single base-pair can be separated by denaturing gradient gel electrophoresis due to the sequence-specific melting properties of DNA. Base modifications such as methylation are also known to affect the melting temperature of DNA. We examined the final position of DNA fragments containing either 5-methyl-cytosine or 6-methyl-adenine in denaturing gradient gels. The presence of a single methylated base within an early melting domain resulted in a well-resolved shift in fragment position relative to the unmethylated sequence. In addition, fragments containing hemimethylated and fully methylated sites could be distinguished, and a proportionally larger shift was observed with an increasing number of methylated bases. Denaturing gradient gel electrophoresis thus provides a sensitive method for analyzing the methylation state of DNA, which is not dependent on the presence of restriction enzyme cleavage sites. We also demonstrate that denaturing gradient gel electrophoresis can be used to obtain a quantitative estimate of the change in helix stability caused by modification of one or two bases in a complex DNA sequence. Such estimates should allow more accurate modeling of melting of natural DNA sequences.     

Nearly all single base substitutions in DNA fragments joined to a GC-clamp can be detected by denaturing gradient gel electrophoresis.

Duplex DNA fragments differing by single base substitutions can be separated by electrophoresis in denaturing gradient polyacrylamide gels, but only substitutions in a restricted part of the molecule lead to a separation (1). In an effort to circumvent this problem, we demonstrated that the melting properties and electrophoretic behavior of a 135 base pair DNA fragment containing a beta-globin promoter are changed by attaching a GC-rich sequence, called a 'GC-clamp' (2). We predicted that these changes should make it possible to resolve most, if not all, single base substitutions within fragments attached to the clamp. To test this possibility we examined the effect of several different single base substitutions on the electrophoretic behavior of the beta-globin promoter fragment in denaturing gradient gels. We find that the GC-clamp allows the separation of fragments containing substitutions throughout the promoter fragment. Many of these substitutions do not lead to a separation when the fragment is not attached to the clamp. Theoretical calculations and analysis of a large number of different mutations indicate that approximately 95% of all possible single base substitutions should be separable when attached to a GC-clamp.     

Improvement of fragment and primer selection for mutation detection by denaturing gradient gel electrophoresis.

Denaturing gradient gel electrophoresis (DGGE) is one of the most powerful methods for mutation detection currently available. For successful application the appropriate selection of PCR fragments and PCR primers is crucial. The sequence of interest should always be within the domain with the lowest melting temperature. When more than one melting domain is present the fragment is generally divided into several smaller ones. This, however, is not always necessary. We found that simple modifications of PCR fragments and primer sequences may substantially reduce the number of amplicons required. Furthermore, by plotting the (natural) melting curves of fragments without a GC-clamp, we could explain why fragments theoretically perfect for DGGE in practice failed to reveal mutations. Alternative fragment selection and the use of modified primers (addition of T/A or G/C tails) result in the detection of mutations that originally remained undetected. Our studies extend the utility of DGGE by using a minimum of PCR fragments and achieving a maximum of mutation detection.     

Detecting base pair substitutions in DNA fragments by temperature-gradient gel electrophoresis.

A vertical gel electrophoresis apparatus is described which can distinguish DNA fragments differing by single base pair substitutions. The system employs a homogenous polyacrylamide gel containing urea-formamide and a temperature gradient which runs either perpendicular or parallel to the direction of electrophoresis. The temperature-gradient system simplifies several features of the denaturant-gradient system (1) and is relatively inexpensive to construct. Eight homologous 373 bp DNAs differing by one, two, or nine base pair substitutions were examined. DNA electrophoretic mobility changed abruptly with the temperature induced unwinding of DNA domains. GC to AT substitutions at different locations within the first melting domain, as well as an AT to TA transversion were separated with temperature gradients parallel to the electrophoretic direction. The relative stabilities of the DNAs observed in the gels were compared to predictions of DNA melting theory. General agreement was observed however complete correspondence was not obtained.     

Site and sequence specificity of the daunomycin-DNA interaction

The site and sequence specificity of the daunomycin-DNA interaction was examined by equilibrium binding methods, by deoxyribonuclease I footprinting studies, and by examination of the effect of the antibiotic on the cleavage of linearized pBR322 DNA by restriction endonucleases PvuI and EcoRI. These three experimental approaches provide mutually consistent results showing that daunomycin indeed recognizes specific sites along the DNA lattice. The affinity of daunomycin toward natural DNA increases with increasing GC content. The quantitative results are most readily explained by binding models in which daunomycin interacts with sites containing two adjacent GC base pairs, possibly occurring as part of a triplet recognition sequence. Deoxyribonuclease I footprinting studies utilizing the 160 base pair (bp) tyrT DNA fragment and 61 and 53 bp restriction fragments isolated from pBR322 DNA further define the sequence specificity of daunomycin binding. Specific, reproducible protection patterns were obtained for each DNA fragment at 4 degrees C. Seven protected sequences, ranging in size from 4 to 14 bp, were identified within the tyrT fragment. Relative to the overall tyrT sequence, these protected sequences were GC rich and contained a more limited and distinct distribution of di- and trinucleotides. Within all of the protected sequences, a triplet containing adjacent GC base pairs flanked by an AT base pair could be found in one or more copies. Nowhere in the tyrT fragment did that triplet occur outside a protected sequence. The same triplet occurred within seven out of nine protected sequences observed in the fragments isolated from pBR322 DNA. In the two remaining cases, three contiguous GC base pairs were found. We conclude that the preferred daunomycin triplet binding site contains adjacent GC base pairs, of variable sequence, flanked by an AT base pair. This conclusion is consistent with the results of a recent theoretical study of daunomycin sequence specificity [Chen, K.-X., Gresh, N., & Pullman, B. (1985) J. Biomol. Struct. Dyn. 3, 445-466]. Adriamycin and the beta-anomer of adriamycin produce the same qualitative pattern of protection as daunomycin with the tyrT fragment. Daunomycin inhibits the rate of digestion of pBR322 DNA by PvuI (recognition sequence 5'-CGATCG-3') to a greater extent than it does EcoRI (recognition sequence 5'-GAATTC-3'), a finding consistent with the conclusions derived from our footprinting studies. Our results, as a whole, are the clearest indication to date that daunomycin recognizes a specific DNA sequence as a preferred binding site.     

Z-curve: a computer program calculating DNA helical axis coordinates for three-dimensional graphic presentation of curvature.

In order to predict curvature of DNA fragments, we previously developed a computer program for simply calculating a vectorial sum of all individual roll, tilt and twist wedge angles between the nearest base pairs for a given DNA fragment [Lee et al., (1991)]. Now, a new program, called Z-curve, was developed to calculate three-dimensional coordinates of the helical center of each base pair along the DNA, using helical axis deviations from B-form DNA by wedge angles. The output file of the new program was designed to become an input file for a graphics program, Insight II. Thus, we were able to obtain three-dimensional graphic presentations of DNA helical axis curvatures of any length. It visualized spatial details of the DNA curvature, where and how much it curves, and to which direction. It also allowed calculation of the three-dimensional distance between two ends of a DNA fragment, which could provide a measure of its curvature. Here, three DNA fragments, both curved and straight, were subjected to the Z-curve and Insight II programs. The results showed that their curvature details could be visualized to the level of the base pair, whether the DNA fragments contained an oligo(A) track or not. Their estimated curvatures were consistent with the experimental results of permutation gel mobility assay.     

Influence of nearest neighbor sequence on the stability of base pair mismatches in long DNA; determination by temperature-gradient gel electrophoresis.

Temperature-gradient gel electrophoresis (TGGE) was employed to determine the thermal stabilities of 48 DNA fragments that differ by single base pair mismatches. The approach provides a rapid way for studying how specific base mismatches effect the stability of a long DNA fragment. Homologous 373 bp DNA fragments differing by single base pair substitutions in their first melting domain were employed. Heteroduplexes were formed by melting and reannealing pairs of DNAs, one of which was 32P-labeled on its 5'-end. Product DNAs were separated based on their thermal stability by parallel and perpendicular temperature-gradient gel electrophoresis. The order of stability was determined for all common base pairs and mismatched bases in four different nearest neighbor environments; d(GXT).d(AYC), d(GXG).d(CYC), d(CXA).d(TYG), and d(TXT).d(AYA) with X,Y = A, T, C, or G. DNA fragments containing a single mismatch were destabilized by 1 to 5 degrees C with respect to homologous DNAs with complete Watson-Crick base pairing. Both the bases at the mismatch site and neighboring stacking interactions influence the destabilization caused by a mismatch. G.T, G.G and G.A mismatches were always among the most stable mismatches for all nearest neighbor environments examined. Purine.purine mismatches were generally more stable than pyrimidine.pyrimidine mispairs. Our results are in very good agreement with data where available from solution studies of short DNA oligomers.     

Restriction endonucleases and their applications

Restriction endonucleases are deoxyribonucleases which cleave double-stranded DNA into fragments. With only one exception, all restriction endonucleases recognize short, non-methylated DNA sequences. Restriction endonucleases can be divided into two groups based on the position of the cleavage site relative to the recognition sequence. Class I restriction endonucleases cleave double-stranded DNA at positions outside the recognition sequence and generate fragments of random size. The cleavage sites of Class II restriction endonucleases are located, in most cases, within the recognition sequence. Most of the Class II restriction endonucleases recognize 4, 5, or 6 base pair palindromes and generate fragments with either flush ends or staggered ends. DNA fragments with staggered ends contain 3, 4, or 5 nucleotide single-stranded tails called ‘sticky ends’. DNA fragments produced by Class II restriction endonuclease cleavage can be separated on gels according to their molecular weight. The fragments can be isolated from the gel and used for sequence analysis to elucidate genetic information stored in DNA. Further, an isolated fragment can be inserted into a small extrachromosomal DNA, e.g. plasmid, phage or viral DNA, and its replication and expression can be studied in clones of prokaryotic or eukaryotic cells. Restriction endonucleases and cloning technology are powerful modern tools for attacking genetic problems in medicine, agriculture and industrial microbiology.     

Sequence-dependent variability of DNA structure. Influence of flanking sequences and fragment length on digestion by conformationally sensitive nucleases

DNase I and 1,10-phenanthroline-copper are two nucleolytic activities which are sequence-dependent in their scission reaction yet are not nucleotide-specific at their site of cutting. When these two nucleases are used to digest identical sequences in 18-base pair oligonucleotides and in restriction fragments 10-fold longer, the digestion patterns are similar at sequence positions in the interior of the fragment. Changes in reactivity to 1,10-phenanthroline-copper associated with mutational changes in the lac promoter in biochemically functional restriction fragments are duplicated in 18-base pair oligonucleotides. The structural variability of a given DNA sequence detected by these conformationally sensitive nucleolytic activities is therefore encoded in local sequence and not sensitive to fragment length. Digestion patterns of a repeated 7-base pair sequence within a longer sequence have the same characteristic except for the two nucleotides at the 5' periphery of the direct repeat. This conclusion is based on the digestion pattern of a restriction fragment which contains the polyadenylation site of the mouse immunoglobulin mu heavy chain gene. Two pairs of different 7-base pair sequences repeated in this fragment retain their distinctive digestion patterns. DNA sequences which comprise the binding sites of regulatory proteins, retain a characteristic structure only influenced at their peripheries by two to three bases of the flanking sequence.     

Application of constant denaturant capillary electrophoresis (CDCE) to mutation detection in humans

Constant denaturant electrophoresis is a DNA separation technique based on the principle of cooperative melting equilibrium. DNA sequences with distinct high and low melting domains can be utilized to separate and identify molecules differing by only one base pair in the lower melting domain. Combined with capillary gel electrophoresis and when coupled with high fidelity DNA amplification, this approach can detect mutants at a fraction of 10⁻⁶. Modifications to the capillary elecctrophoretic system have also increased DNA loading capacity which allows for analysis of rare mutations in a large, heterogeneous population such as DNA samples derived from human tissues. Employment of this technology has determined the first mutational spectrum in human cells and tissues in a mitochondrial sequence without phenotypic selection of mutants.     

Effects of mutations within the herpes simplex virus type 1 DNA encapsidation signal on packaging efficiency

The cis-acting signals required for cleavage and encapsidation of the herpes simplex virus type 1 genome lie within the terminally redundant region or a sequence. The a sequence is flanked by short direct repeats (DR1) containing the site of cleavage, and quasi-unique regions, Uc and Ub, occupy positions adjacent to the genomic L and S termini, respectively, such that a novel fragment, Uc-DR1-Ub, is generated upon ligation of the genomic ends. The Uc-DR1-Ub fragment can function as a minimal packaging signal, and motifs have been identified within Uc and Ub that are conserved near the ends of other herpesvirus genomes (pac2 and pac1, respectively). We have introduced deletion and substitution mutations within the pac regions of the Uc-DR1-Ub fragment and assessed their effects on DNA packaging in an amplicon-based transient transfection assay. Within pac2, mutations affecting the T tract had the greatest inhibitory effect, but deletion of sequences on either side of this element also reduced packaging, suggesting that its position relative to other sequences within the Uc-DR1-Ub fragment is likely to be important. No single region essential for DNA packaging was detected within pac1. However, mutants lacking the G tracts on either side of the pac1 T-rich motif exhibited a reduced efficiency of serial propagation, and alteration of the sequences between DR1 and the pac1 T element also resulted in defective generation of Ub-containing terminal fragments. The data are consistent with a model in which initiation and termination of packaging are specified by sequences within Uc and Ub, respectively.     

Denaturing gradient gel method for mapping single base changes in human mitochondrial DNA.

A denaturing gradient gel electrophoresis (DGGE) method is described that detects even single base pair changes in mitochondrial DNA (mtDNA). In this method, restriction fragments of mtDNA are electrophoresed in a urea/formamide gradient gel at 60 degrees C. Migration distance of each mtDNA fragment in the gel depends on melting behavior which reflects base composition. Fragments are located by Southern blotting with specific mtDNA probes. With just four carefully chosen restriction enzymes and as little as 50-100 ng of mtDNA, the method covers almost the entire human mitochondrial genome. To demonstrate the method, human mtDNA was analyzed. In six normal individuals, DGGE revealed melting behavior polymorphisms (MBPs) in mtDNA fragments that were not detected by restriction fragment length polymorphism (RFLP) analysis in agarose gels. Another individual, shown to have a melting behavior polymorphism in the cytochrome b coding region, was studied in detail. By mapping, the mutation was deduced to lie between nt 14905 and 15370. The affected fragment was amplified by PCR and sequenced. Specific base changes were identified in the region predicted by the gel result. This method will be especially useful as a diagnostic tool in mitochondrial disease for rapid localization of mtDNA mutations to specific regions of the genome, but DGGE also could complement RFLP analysis as a more sensitive method to follow maternal lineage in human and animal populations in a variety of research fields.     

Specificity of mutagenesis by 4-aminobiphenyl. A possible role for N-(deoxyadenosin-8-yl)-4-aminobiphenyl as a premutational lesion

Mutagenesis by N-acetoxy-N-trifluoroacetyl-4-aminobiphenyl, a reactive form of the human bladder carcinogen 4-aminobiphenyl (ABP), was studied in Escherichia coli virus M13mp10. N-acetoxy-N-trifluoroacetyl-4-ABP-treated DNA containing 140 lesions/duplex genome, when introduced into excision repair-competent cells induced for SOS mutagenic processing, resulted in a 40-fold increase in mutation frequency over background in the lacZ alpha gene fragment. DNA sequence changes were determined for 20 independent mutants. G-C base pairs were the major targets for base pair substitution mutations, although significant mutagenic activity was also observed at certain A-T base pairs. Deletion and frameshift mutations also were found in this sample. The salient feature of this partial "mutational spectrum" was a hotspot that occurred at position 6357 (amino acid 30 of the beta-galactosidase fragment encoded by M13mp10); this A-T to T-A transversion appeared in 6 of the 20 mutants. The property of ABP to mutate A-T base pairs was consistent with the result that N-hydroxy-ABP reverted Salmonella typhimurium strain TA104, which is presumed to revert primarily due to mutations at these sites. The ability of the major carcinogen-DNA adduct formed by ABP in vivo and in vitro, N-(deoxyguanosin-8-yl)-4-aminobiphenyl, to cause base pair substitution mutations was also investigated. This adduct was positioned specifically in the minus strand at position 6270 in duplex M13mp10 DNA. In the presence of the mutagenesis-enhancing plasmid pGW16 and UV induction of SOS mutagenic processing, it was shown that fewer than 0.02% of the adducts resulted in transition or transversion mutations following transfection of DNA into excision-repair competent cells. Similar results were obtained in uvrA and uvrC backgrounds. Although the major adduct did not cause base substitution mutations under these experimental conditions, the contribution of this lesion to the entire spectrum of mutations in the lacZ alpha fragment seems likely.