|
|
|
|
Genome Res. 13:2316-2324, 2003 ©2003 by Cold Spring Harbor Laboratory Press; ISSN 1088-9051/03 $5.00 Methods DNA Enrichment by Allele-Specific Hybridization (DEASH): A Novel Method for Haplotyping and for Detecting Low-Frequency Base Substitutional Variants and Recombinant DNA MoleculesDepartment of Genetics, University of Leicester, Leicester LE1 7RH, UK
Detecting rare sequence variants in genomic DNA is central to the analysis of de novo mutation and recombination events and the detection of rare pathological mutations in mixed cell populations. Current PCR techniques suffer from noise that limits detection to variants present at a frequency of at least 10-4-10-5 per cell. We nowdescribe an alternative approach that recovers genomic DNA molecules containing a known single-nucleotide variant by hybridization selection using a biotinylated allele-specific oligonucleotide, followed by hybrid capture on streptavidin-coated paramagnetic beads and subsequent analysis by PCR. This technique of DNA enrichment by allele-specific hybridization (DEASH) is fast, effective for all tested single-nucleotide polymorphisms (SNPs), and can recover large (>10 kb) single-stranded molecules. A single round of DEASH is effective in separating haplotypes from genomic DNA and can not only readily detect and validate DNA molecules containing a single base change at a frequency of 10-5 per cell, but can also place these changes within the context of an extended haplotype. This technique offers a newapproach to the analysis of mutation and recombination, and has the potential to detect very rare de novo base substitutions.
The detection of sequence variants present at very low frequency in genomic DNA remains a major challenge in molecular genetics, of direct relevance to the detection of germline and somatic de novo mutation events and other forms of spontaneous genomic rearrangement (Tiemann-Boege et al. 2002
Many methods for detecting low-frequency variants in genomic DNA use the polymerase chain reaction (PCR) to amplify mutant plus wild-type target, followed by interrogating these products by sequencing, oligonucleotide ligation, restriction digestion, mass spectrometry, or hybridization with allele-specific oligonucleotides to identify the variant against the background of wild-type DNA (see, e.g., Saiki et al. 1986
Most mutation detection systems yield an assay signal that is difficult to validate in terms of the number of mutant molecules detected. This can be overcome in part by analyzing multiple samples, each containing limited DNA (typically 50 genome equivalents), yielding direct information on numbers of mutant molecules (digital PCR; Vogelstein and Kinzler 1999
The unifying problem behind all of these PCR approaches for detecting rare variants is replication infidelity during amplification. The solution therefore is to collect mutant DNA molecules from genomic DNA prior to analyzing them by PCR. We now describe a rapid method that achieves this goal and which should be applicable to any base-substitutional variant. This method is a modification of traditional techniques for enriching nucleic acids by hybridization with biotinylated DNA probes (see, e.g., Abe 1992
The Principles of DEASH The DEASH procedure is outlined in Figure 1. The input is double-stranded DNA containing a mixture of molecules that differ in the chosen single-base substitution; this mixture could be due to SNP heterozygosity in a single individual, or could represent a mixed DNA sample, or be due to de novo mutation. This DNA is mixed with a biotinylated allele-specific oligonucleotide (bio-ASO) directed to the chosen variant, plus a competitor ASO complementary to the other allele to increase the specificity of hybridization. The mixture is thermally denatured and hybridized, and bio-ASO/target hybrids are captured on streptavidin-coated paramagnetic Dynabeads. The captured hybrid is then thermally eluted from the beads in low-salt buffer, yielding single-stranded target enriched for the chosen variant and free from bio-ASO. Additional target can be recovered by adding bio-ASO to the unbound DNA and subjecting it to further cycles of re-extraction using DEASH. In addition, the enriched target can be purified further simply by adding buffer, bio-ASO, and competitor ASO and proceeding with a second round of DEASH (reenrichment; Fig. 1).
To minimize damage to target DNA, the DEASH protocol uses a modified high-pH buffer based on a standard PCR buffer (Jeffreys et al. 1988
Separation of Haplotypes From PCR Products
Separating Haplotypes From Genomic DNA
DEASH recovers only one of the two complementary strands of a target, with the other strand remaining in the unbound genomic DNA fraction (Fig. 1). Because both strands serve as templates for the PCR assay, it was not possible to assay target depletion from the unbound DNA during purification. To circumvent this problem and directly test levels of depletion, we recovered single-stranded DNA by selection with bio-BbG and used this as target in an additional round of enrichment with re-extractions using bio-BbG. Assaying both unbound and bound DNA showed that most of the target was recovered by a single cycle of enrichment; after two further cycles of re-extraction, about 50% of starting target was recovered in a PCR-amplifiable form, with less than 5% of target DNA remaining in the unbound DNA fraction (Fig. 3C). Proportionately more target was recovered in the first cycle of extraction compared to genomic DNA (Fig. 3, cf. B and C). This suggests that low yields seen particularly with unrestricted genomic DNA may be due to inter-molecular hybridization between repeat DNA that could sterically hinder binding of the bio-ASO/target hybrid to the streptavidin beads. Such hybridization would also create high-molecular-weight aggregates that might promote shearing of the bound target, again reducing yields of intact target DNA.
Bulk Fractionation of Genomic DNA
This MHC target region contained multiple additional SNP heterozygosities in the individual tested (Fig. 4A). Linkage phase of these SNPs relative to the JJK8 selector site could be easily determined by ASO hybridization. Analysis of each SNP on dot blots of PCR products from the enriched target allowed complete haplotypes to be rapidly established (data not shown).
Properties of Bio-ASOs Used in DEASH
Use of DEASH to Detect Low-Frequency Changes in Human DNA
DEASH is a powerful, flexible method for fractionating not only PCR products but also genomic DNA for known base substitutional variants. It can be applied equally efficiently to small amounts (100 ng or less) of genomic DNA and to bulk fractionation at the 10 µg level (data not shown), and could be readily scaled to process even larger amounts of DNA. DEASH provides a simple, rapid method for haplotype separation at the 10-kb scale, requiring only the presence of a known SNP heterozygosity and providing an approach to haplotype determination that complements traditional methods using pedigrees, allele-specific PCR, or chromosome separation in somatic cell hybrids (Papadopoulos et al. 1995 Analysis of sperm recombination events has shown that just one round of DEASH enrichment giving 100-fold purification is sufficient to detect single base changes (arising by gene conversion) at a frequency of 10-5 per cell (Fig. 5; additional data not shown). Without fractionation, the conversion events detected in Figure 5 would have required the analysis of 2000 PCR reactions, each containing 50 molecules (using ASO hybridization, this is the pool size limit for reliably detecting a mutant molecule against the background of wild-type molecules). Furthermore, such large numbers of reactions would generate the occasional false "mutant" signal due to PCR misincorporation during the earliest rounds of PCR amplification, generating a background noise that would mask any low-frequency events. DEASH circumvents these problems by greatly reducing the numbers of molecules screened for mutation, and further validates the mutants by showing that a given mutant is only detected if DNA has been enriched using the relevant bio-ASO.
The sensitivity of DEASH in detecting rare variants reaches or surpasses the detection limits of all other mutation detection systems based on direct analysis of DNA. It is already sufficiently great to enable in principle the direct detection of the de novo germline mutations in sperm DNA that cause achondroplasia and Apert syndrome, both of which show remarkably high spontaneous mutation rates, of 3x10-5 and 1x10-5-5x10-5, respectively (Cohen et al. 1992 Quantitation of the rate of gene conversion events detected in sperm DNA was straightforward, given that the enrichment separated the two progenitor haplotypes and that the number of progenitor haplotype molecules in the purified DNA assayed for recombination was known (Fig. 5). In contrast, the enrichment and quantitation of de novo mutant molecules, or low-frequency variant molecules present in a mixed DNA sample, will involve the removal of almost all wild-type DNA and will not give information on the efficiency of recovery of mutant molecules. It will therefore be necessary to develop reporter constructs carrying the mutation that can be added to the starting genomic DNA, providing an internal control from which yields and thus mutation rates can be determined. DEASH enrichment is flexible and could be readily extended to sequential enrichment using two or more different bio-ASOs. For instance, it should be possible to recover recombinant molecules without using PCR. Thus, if recombinants are required between haplotypes A-B and a-b, then bio-ASO fractionation with marker A followed by the repulsion-phase marker b would enrich for A-b recombinants. Another feature of DEASH is that it should be readily amenable to automation, allowing the parallel processing of multiple samples in, for example, the detection of low-frequency mutations in patient samples. Finally, DEASH is currently the only technology with the potential to detect directly very-low-frequency de novo base substitutional mutations in the human germline, an area of major importance in understanding the dynamics of heri-table mutation and the origin of pathological mutations.
Oligonucleotides All bio-ASOs and competitor ASOs were 18-nt-long, with the allele-specific base 11 nt from the 5' end. 5' biotinylated ASOs were HPLC-purified before use. Single- and triple-biotinylated ASOs (Oswel) gave very similar yields in DEASH, and most fractionations were performed with single-biotinylated ASOs. MS32 PCR primers used were: A 5'-GTC AAG GCA GAT TTG AAG GGA G-3'; B 5'-GGT AGC CAC ATC TCT GCA TGA TCT C-3'; C 5'-CCA TGC CCA GCC TTC ATT CTT TTC-3'; D 5'-CTT CCT CGT TCT CCT CAG CCC TAG-3'. HLA-DMB primers were: E 5'-AGG AGC CCC CCA TCT GGG AGT AGT-3'; F 5'-GTC TGC CTC AGG TGC TGA CC-3'; G 5'-ATG GCT CCC ATG CTC CCT GG-3'; JJ6C 5'-CTG CTC TGG TGG TGT GGC-3'; JJ6T 5'-GCT CTG GTG GTG TGG T-3'; JJ7C 5'-CCC CCT TGC TTT GAA ATG AGG C-3'.
The DEASH Procedure
Analysis of Fractionated DNA
This work is the subject of a UK Patent Application. We thank J. Blower for providing semen samples, S. Mistry for oligonucleotide synthesis, and colleagues for helpful discussions. This work was supported by grants to A.J.J. from the Medical Research Council and Royal Society. The publication costs of this article were defrayed in part by payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 USC section 1734 solely to indicate this fact.
1 Corresponding author. E-MAIL ajj{at}le.ac.uk; FAX 44-116-252-3378. Article and publication are at http://www.genome.org/cgi/doi/10.1101/gr.1214603.
Abe, K. 1992. Rapid isolation of desired sequences from lone linker PCR amplified cDNA mixtures: Application to identification and recovery of expressed sequences in cloned genomic DNA. Mamm. Genome 2: 252-259.[CrossRef][Medline] Adinolfi, M., Camporese, C., and Carr, T. 1989. Gene amplification to detect fetal nucleated cells in pregnant women. Lancet 2: 328-329. Bellus, G.A., Hefferon, T.W., Ortiz de Luna, R.I., Hecht, J.T., Horton, W.A., Machado, M., Kaitila, I., McIntosh, I., and Francomano, C.A. 1995. Achondroplasia is defined by recurrent G380R mutations of FGFR3. Am. J. Hum. Genet. 56: 368-373.[Medline] Berdanier, C.D. and Everts, H.B. 2001. Mitochondrial DNA in aging and degenerative disease. Mutat. Res. 475: 169-183.[Medline] Boyle, J.O., Mao, L., Brennan, J.A., Koch, W.M., Eisele, D.W., Saunders, J.R., and Sidransky, D. 1994. Gene mutations in saliva as molecular markers for head and neck squamous cell carcinomas. Am. J. Surg. 168: 429-432.[CrossRef][Medline] Cha, R.S., Zarbl, H., Keohavong, P., and Thilly, W.G. 1992. Mismatch amplification mutation assay (MAMA): Application to the c-H-ras gene. PCR Methods Appl. 2: 14-20. Cohen, M.M., Kreiborg, S., Lammer, E.J., Cordero, J.F., Mastroiacovo, P., Erickson, J.D., Roeper, P., and Martinez-Frias, M.L. 1992. Birth prevalence study of the Apert syndrome. Am. J. Med. Genet. 42: 655-659.[CrossRef][Medline] Crow, J.F. 1999. Spontaneous mutation in man. Mutat. Res. 437: 5-9.[CrossRef][Medline] Czeizel, A.E., Elek, C., and Susanszky, E. 1993. Birth prevalence study of Apert syndrome. Am. J. Med. Genet. 45: 392-393.[CrossRef][Medline] Fox, J.C., England, J., White, P., Ellison, G., Callaghan, K., Charlesworth, N.R., Hehir, J., McCarthy, T.L., Smith-Ravin, J., Talbot, I.C., et al. 1998. The detection of K-ras mutations in colorectal cancer using the amplification-refractory mutation system. Br. J. Cancer 77: 1267-1274.[Medline] Gerry, N.P., Witowski, N.E., Day, J., Hammer, R.P., Barany, G., and Barany, F. 1999. Universal DNA microarray method for multiplex detection of low abundance point mutations. J. Mol. Biol. 292: 251-262.[CrossRef][Medline]
Gibbs, R.A., Nguyen, P.N., and Caskey, C.T. 1989. Detection of single DNA base differences by competitive oligonucleotide priming. Nucleic Acids Res. 17: 2437-2448.
Haliassos, A., Chomel, J.C., Grandjouan, S., Kruh, J., Kaplan, J.C., and Kitzis, A. 1989. Detection of minority point mutations by modified PCR technique: A new approach for a sensitive diagnosis of tumor-progression markers. Nucleic Acids Res. 17: 8093-8099. Haque, K., Hehir, J., Fox, J.C., Newton, C.R., and Little, S. 1998. Amplification refractory mutation system linear extension: A novel, gel-free, enzyme-linked immunoassay method for DNA genotyping. Diagn. Mol. Pathol. 7: 248-252.[CrossRef][Medline] Hengstschlager, M., Holzl, G., Ulm, B., and Bernaschek, G. 1997. Raising the sensitivity of fetal RhD typing and sex determination from maternal blood. J. Med. Genet. 34: 350-351.
Jeffreys, A.J., Wilson, V., Neumann, R., and Keyte, J. 1988. Amplification of human minisatellites by the polymerase chain reaction: Towards DNA fingerprinting of single cells. Nucleic Acids Res. 16: 10953-10971. Jeffreys, A.J., Murray, J., and Neumann, R. 1998a. High-resolution mapping of crossovers in human sperm defines a minisatellite-associated recombination hotspot. Mol. Cell 2: 267-273.[CrossRef][Medline] Jeffreys, A.J., Neil, D.L., and Neumann, R. 1998b. Repeat instability at human minisatellites arising from meiotic recombination. EMBO J. 17: 4147-4157.[CrossRef][Medline]
Jeffreys, A.J., Ritchie, A., and Neumann, R. 2000. High-resolution analysis of haplotype diversity and meiotic crossover in the human TAP2 recombination hotspot. Hum. Mol. Genet. 9: 725-733. Jeffreys, A.J., Kauppi, L., and Neumann, R. 2001. Intensely punctate meiotic recombination in the class II region of the major histocompatibility complex. Nat. Genet. 29: 217-222.[CrossRef][Medline] Knight, J.A. 2000. The biochemistry of aging. Adv. Clin. Chem. 35: 1-62.[Medline] Laken, S.J., Jackson, P.E., Kinzler, K.W., Vogelstein, B., Strickland, P.T., Groopman, J.D., and Friesen, M.D. 1998. Genotyping by mass spectrometric analysis of short DNA fragments. Nat. Biotechnol. 16: 1352-1356.[CrossRef][Medline] Lo, Y.M., Bowell, P.J., Selinger, M., Mackenzie, I.Z., Chamberlain, P., Gillmer, M.D., Elliott, P., Pratt, G., Littlewood, T.J., Fleming, K.A., et al. 1994. Prenatal determination of fetal rhesus D status by DNA amplification of peripheral blood of rhesus-negative mothers. Ann. NY Acad. Sci. 731: 229-236.[Abstract]
Mao, L., Hruban, R.H., Boyle, J.O., Tockman, M., and Sidransky, D. 1994. Detection of oncogene mutations in sputum precedes diagnosis of lung cancer. Cancer Res. 54: 1634-1637. May, C.A., Shone, A.C., Kalaydjieva, L., Sajantila, A., and Jeffreys, A.J. 2002. Crossover clustering and rapid decay of linkage disequilibrium in the Xp/Yp pseudoautosomal SHOX gene. Nat. Genet. 31: 272-275.[CrossRef][Medline]
Mills, N.E., Fishman, C.L., Scholes, J., Anderson, S.E., Rom, W.N., and Jacobson, D.R. 1995. Detection of K-ras oncogene mutations in bronchoalveolar lavage fluid for lung cancer diagnosis. J. Natl. Cancer Inst. 87: 1056-1060.
Nachman, M.W., and Crowell, S.L. 2000. Estimate of the mutation rate per nucleotide in humans. Genetics 156: 297-304.
Newton, C.R., Graham, A., Heptinstall, L.E., Powell, S.J., Summers, C., Kalsheker, N., Smith, J.C., and Markham, A.F. 1989. Analysis of any point mutation in DNA. The amplification refractory mutation system (ARMS). Nucleic Acids Res. 17: 2503-2516. Papadopoulos, N., Leach, F.S., Kinzler, K.W., and Vogelstein, B. 1995. Monoallelic mutation analysis (MAMA) for identifying germline mutations. Nat. Genet. 11: 99-102.[CrossRef][Medline] Reed, W., Kong, D.Z., Lee, T.H., Cowan, M.J., Busch, M.P., and Baxter-Lowe, L.A. 2002. Noninvasive determination of the paternal HLA haplotype of a fetus using kinetic PCR to detect fetal microchimerism in maternal plasma. Bone Marrow Transplant 29: 527-529.[CrossRef][Medline]
Saiki, R.K., Bugawan, T.L., Horn, G.T., Mullis, K.B., and Erlich, H.A. 1986. Analysis of enzymatically amplified
Sidransky, D., Von Eschenbach, A., Tsai, Y.C., Jones, P., Summerhayes, I., Marshall, F., Paul, M., Green, P., Hamilton, S.R., Frost, P., et al. 1991. Identification of p53 gene mutations in bladder cancers and urine samples. Science 252: 706-709.
Sidransky, D., Tokino, T., Hamilton, S.R., Kinzler, K.W., Levin, B., Frost, P., and Vogelstein, B. 1992. Identification of ras oncogene mutations in the stool of patients with curable colorectal tumors. Science 256: 102-105.
Smith-Ravin, J., England, J., Talbot, I.C., and Bodmer, W. 1995. Detection of c-Ki-ras mutations in faecal samples from sporadic colorectal cancer patients. Gut 36: 81-86.
Tada, M., Omata, M., Kawai, S., Saisho, H., Ohto, M., Saiki, R.K., and Sninsky, J.J. 1993. Detection of ras gene mutations in pancreatic juice and peripheral blood of patients with pancreatic adenocarcinoma. Cancer Res. 53: 2472-2474.
Tiemann-Boege, I., Navidi, W., Grewal, R., Cohn, D., Eskenazi, B., Wyrobek, A.J., and Arnheim, N. 2002. The observed human sperm mutation frequency cannot explain the achondroplasia paternal age effect. Proc. Natl. Acad. Sci. 99: 14952-14957. Vijg, J. 2000. Somatic mutations and aging: A re-evaluation. Mutat. Res. 447: 117-135.[Medline]
Vogelstein, B. and Kinzler, K.W. 1999. Digital PCR. Proc. Natl. Acad. Sci. 96: 9236-9241.
Zhang, J., Fidler, C., Murphy, M.F., Chamberlain, P.F., Sargent, I.L., Redman, C.W., Hjelm, N.M., Wainscoat, J.S., and Lo, Y.M. 2000. Determination of fetal RhD status by maternal plasma DNA analysis. Ann. NY Acad. Sci. 906: 153-155.
http://www.le.ac.uk/genetics/ajj/HLA/index.html; MHC SNP data at the University of Leicester. http://www.le.ac.uk/ge/ajj/SHOX; SHOX gene SNP data at the University of Leicester.
Received January 23, 2003;
accepted in revised format August 4, 2003.
This article has been cited by other articles:
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||