Clinical Chemistry 45:7 982–986 (1999)
Molecular Diagnostics and Genetics
TaqMan PCR-based Gene Dosage Assay for Predictive Testing in Individuals from a Cancer Family with INK4 Locus Haploinsufficiency Ingrid Laurendeau, Michel Bahuau, Nicolas Vodovar, Claire Larramendy, Martine Olivi, Ivan Bieche, Michel Vidaud, and Dominique Vidaud*
Background: A genetic syndrome of cutaneous malignant melanoma and nervous system tumors recently has been characterized and shown to be linked to the INK4 locus in the 9p21 region. Hemizygosity at adjacent physically mapped microsatellite markers indicated deletion of p16, p19, and p15 clustered tumor suppressors. Because individuals from this family could benefit from predictive testing in terms of cancer prevention, we developed a direct test without need to analyze parental DNAs to comply with the rules of individual consent and secrecy. Methods: We developed an assay using TaqManTM real-time quantitative PCR, with p15 as the test sequence and albumin (ALB) as the reference gene. The normalized ratio of p15/ALB is expected to yield a value of ;1 in individuals without the deletion, whereas a ratio of ;0.5, indicating p15 haploinsufficiency, is expected in predisposed individuals. Results: All patients harboring the previously defined at-risk haplotype were correctly identified using this approach. In six individuals with deletions, the p15/ALB ratios were 0.472– 0.556 (SD, 0.013– 0.078). In the five individuals without deletions, the ratios were 0.919 – 1.019 (SD, 0.006 – 0.075). Conclusions: This is the first report of a high-throughput, automatable gene dosage assay successfully applied to the identification of a germ-line deletion. This approach, not limited by marker informativeness or the need for harvesting live cells, can be applied to any condition with haploinsufficiency and extended to the characterization of most abnormalities of the ploidy.
A genetic syndrome consisting of cutaneous malignant melanoma (CMM)1 and nervous system tumors (NSTs) recently has been defined on the basis of a five-generation family with a predisposing allele to CMM/dysplastic nevi, astrocytoma/glioblastoma, neurofibroma, schwannoma, and meningioma (1 ). This CMM-NST syndrome recently has been shown to be linked to the INK4 locus, with all patients sharing a common haplotype as defined by genetic markers of the 9p21 region, and has been demonstrated to result from a large germ-line deletion that ablates the whole cluster containing the p16, p19, and p15 tumor suppressor genes (2 ). Such a large DNA lesion was characterized by hemizygosity mapping based on the segregation of adjacent physically mapped microsatellite markers. We are now addressing the question of direct genetic testing aimed at identification of subjects with individual risk among the 100 initially surveyed. Indeed, a direct test that did not require the analysis of parental DNAs, as opposed to indirect microsatellite-based diagnosis, was a prerequisite for complying with the rules of individual consent and secrecy. The recent development of real-time quantitative PCR based on the 59-39 exonuclease activity of the Taq polymerase (3 ), referred as TaqManTM, has offered the opportunity to set up an original gene dosage assay using the ABI PRISM 7700 Sequence Detection System (PerkinElmer Biosystems), in which p15 was selected as the test sequence and albumin (ALB) as the reference disomic gene used to normalize the amounts of input genomic DNAs. Briefly, during the PCR process, a dual-labeled TaqMan probe annealed to the target sequence is cleaved by the 59-39 exonuclease activity of the Taq polymerase, releasing the reporter dye (FAM) from the quencher dye
© 1999 American Association for Clinical Chemistry
Laboratoire de Ge´ne´tique Mole´culaire, Faculte´ des Sciences Pharmaceutiques et Biologiques, Universite´ Paris V, 4 Avenue de l’Observatoire, 75006 Paris, France. *Author for correspondence. Fax 33 1 44 07 17 54; e-mail
[email protected]. Received February 26, 1999; accepted May 3, 1999.
1 Nonstandard abbreviations: CMM, cutaneous malignant melanoma; NST, nervous system tumor; and Ct, threshold cycle.
982
983
Clinical Chemistry 45, No. 7, 1999
(TAMRA). Upon excitation of by an argon laser, the release of the FAM produces an increase in the fluorescent emission, which is captured at 518 nm by a chargedcouple device camera and analyzed through the algorithms of the ABI PRISM 7700 Sequence Detection System computer software (4, 5 ). To develop a direct automatable predictive test in the family with INK4 locus haploinsufficiency, this p15 gene dosage assay was applied to a panel of 11 test, haplotypepositive or -negative, individuals from this large kindred.
Patients and Methods patients, nucleic acids, and reference dna High-molecular weight DNAs were prepared using standard proteinase K digestions followed by phenol-chloroform extractions (6 ) from whole-blood leukocytes or lymphoblastoid cell lines from 11 individuals, including 3 clinically affected and 3 healthy, all haplotype-sharing, relatives, as well as in 5 healthy haplotype-negative relatives, including 1 spouse (Fig. 1). The haplotypes have been defined elsewhere (1 ). High-quality human genomic DNA (Boehringer Mannheim) served as reference DNA.
real-time pcr In the TaqMan approach, one of the main variables, referred as Ct, is defined as the fractional cycle number at which the fluorescence generated by cleavage of the probe (DRn) crosses a fixed threshold. The target gene copy
Table 1. Primers and probes. Primer and probe
p15 gene p15-F p15-R p15-probe ALB gene Alb-F Alb-R Alb-probe
Sequence
59AGG CCA TAG CAG GCA GAT CAC CT 39 59AGC TCA CTG CAA CCT CTT CCA GCA T 39 59AGG AGG CTG AGG CAG GAG AAT CGC TT 39 59GCT GTC ATC TCT TGT GGG CTG T 39 59ACT CAT GGG AGC TGC TGG TTC 39 59GGA GAG ATT TGT GTG GGC ATG ACA GG 39
number in unknown samples is inferred by plotting the Ct value against a calibration curve. To correctly determine the starting copy number regardless of the precise amounts and qualities of input genomic DNAs, we also quantified an internal control gene (ALB) in each single reaction. The normalized gene dose, N, is given by the following ratio: N 5
Starting copy number of target gene (p15) Starting copy number of reference gene (ALB)
The 7700 Sequence Detection System software automatically determines the Ct value and infers the starting copy number in each sample. The real-time PCR and the calibration curve are presented in Fig. 2 for ALB. Normalized gene doses N were then determined for each sample in three independent assays, and the corresponding means were calculated. In this method, a nondeleted test
Fig. 1. Test individuals from the CMM-NST family. Squares refer to males, and circles refer to females. Hatched symbols denote the development of CMM or NST. The haplotype status of each test relative is indicated beneath his or her individual symbol, INK4 refers to the at-risk haplotype, and wt refers to any other wild-type haplotype, as previously defined and according to an individual numbering system based on a pedigree published elsewhere (2 ).
984
Laurendeau et al.: TaqMan Gene Dosage for Direct Predictive Testing
Fig. 2. Albumin (ALB) gene dosage by real-time PCR. Top, amplification plots for reactions with starting ALB gene copy number of 33 000 (A1, 100 ng), 8250 (A4, 25 ng), 2062 (A7, 6.25 ng), or 515 (A10, 1.56 ng). The cycle number is plotted vs the change in normalized reporter signal (DRn). For each reaction tube, the fluorescence signal of the reporter dye (FAM) is divided by the fluorescence signal of the passive reference dye (ROX) to obtain a ratio defined as the normalized reporter signal (Rn). DRn represents the normalized reporter signal (Rn) minus the baseline signal established in the first 15 PCR cycles. DRn increases during PCR as ALB PCR product copy number increases until the reaction reaches a plateau. Ct represents the fractional cycle number at which a significant increase in Rn above a baseline signal (horizontal black line) can first be detected. Three replicates were performed for each reference DNA sample, but the data for only one are shown here. Bottom, calibration curve plotting log starting copy number vs Ct. The black symbols represent the triplicate PCR amplification of the reference DNA samples and red symbols the triplicate PCR amplification of unknown genomic DNA, all included inside the calibration curve. The copy number of ALB (x) can be calculated as follows: y 5 23.374x 1 40.593, where the Ct value is substituted as y.
sequence is expected to yield a ratio of N 5 1, as opposed to N 5 0.5 when it is heterozygously deleted.
primers and probes Primers and probes were chosen with assistance of the computer programs OligoTM, Ver. 4.0 (National Bio-
sciences) and Primer ExpressTM, Ver. 1.0 (PE Biosystems). Their nucleotide sequences are shown in Table 1. Primers were purchased from Life Technologies, and probes were purchased from PE Biosystems. The TaqMan PCR Core Reagent Kit, MicroAmp Optical TubesTM, and MicroAmp Optical CapsTM were from PE Biosystems.
985
Clinical Chemistry 45, No. 7, 1999
Table 2. Predictive testing based on a p15 TaqMan assay.a Assay 1
Assay 2
b
ALB 3103
b
p15 3103
p15/ ALB
b
ALB 3103
b
p15 3103
Assay 3
p15/ ALB
b
ALB 3103
p15 b 3103
p15/ ALB
CMM/ NST
Haplotype
p15/ALB 6 SD
Deletion
NAc
NA
1.92
1.90
0.987
1.26
1.49
1.187
2.99
2.93
0.980
1.051 6 0.112
NA
II-9
2
2
1.76
1.81
1.024
5.06
5.17
1.021
2.12
2.15
1.012
1.019 6 0.006
2
III-9
2
1
3.97
2.08
0.525
7.14
3.65
0.510
4.49
2.16
0.480
0.505 6 0.045
1
III-10
2
1
2.66
1.28
0.479
3.98
2.03
0.509
3.34
1.45
0.435
0.474 6 0.078
1
III-13
2
2
5.55
5.25
0.947
6.11
6.46
1.056
5.56
5.43
0.977
0.993 6 0.057
2
III-15
1
1
3.35
1.73
0.516
3.05
1.74
0.570
2.58
1.50
0.581
0.556 6 0.063
1
III-21
2
2
4.96
4.97
1.002
2.89
2.71
0.934
4.88
4.87
0.998
0.978 6 0.039
2
III-22
1
1
2.15
1.11
0.515
2.86
1.64
0.573
2.69
1.33
0.495
0.528 6 0.077
1
III-24
2
1
4.98
2.37
0.476
6.64
3.16
0.476
5.64
2.62
0.465
0.472 6 0.013
1
IV-15
2
2
4.17
4.21
1.009
7.42
7.06
0.952
4.93
4.77
0.868
0.943 6 0.075
2
IV-32
1
1
2.19
1.15
0.525
3.17
1.74
0.549
2.80
1.42
0.508
0.527 6 0.039
1
IV-36
2
2
17.78
16.49
0.927
10.46
9.71
0.928
17.83
16.05
0.901
0.919 6 0.017
2
Reference DNA
The clinical status (1, presence of CMM and/or NST; 2, healthy), the presence (1) or absence (2) of the predisposing haplotype, the ALB and p15 copy numbers, and the p15/ALB mean ratios are indicated for three independent assays of samples from 11 individuals of the CMM-NST family [pedigree positions as in Ref. (2)], compared with a reference genomic DNA (Boehringer Mannheim). The inferred deletional status (1, deleted; 2, not deleted) is indicated based on the mean of the three p15/ALB ratios (6SD) of ;0.5 or ;1, respectively. All individuals harboring the predisposing haplotype were correctly found to have the deletion, using the p15 gene dosage assay, irrespective of their clinical status. b The SDs for ALB and p15 were ,1% to 12% of the means (not shown, available from authors). c NA, not applicable. a
pcr amplification TaqMan amplification reactions were carried out in reaction volumes of 50 mL, using the components as supplied in the TaqMan PCR Core Reagent Kit. Each reaction contained 13 TaqMan buffer; 200 nmol/L each primer; 100 nmol/L each corresponding fluorogenic probe; 5 mmol/L MgCl2, 200 mmol/L each dATP, dCTP, and dGTP; 400 mmol/L dUTP; 1.25 U of AmpliTaqTM Gold, and 0.5 U of AmpEraseTM uracil N-glycosylase. Each sample was analyzed in triplicate, using 20 ng of DNA in each reaction. The reference DNA was serially diluted and run in parallel to establish the calibration curve and to infer copy numbers from the cycle thresholds (Cts), assuming a conversion factor of 6.6 pg of DNA per diploid genome. Thermal cycling was initiated with a 2-min incubation at 50 °C, followed by a first denaturation step of 10 min at 95 °C, and then 40 cycles of 95 °C for 15 s and 65 °C for 1 min.
Results The results are presented in Table 2 as p15 gene dosage performed in triplicate and in three independent assays. Although the ALB and p15 genes were quantified on the same plate, the low interassay variability, which the use of a calibration curve helped reduce, allowed us to ascertain these genes in different plates. In addition, accurate results can be obtained with duplicate amplification. The six haplotype-sharing affected or healthy individuals showed INK4 haploinsufficiency, i.e., were heterozygotes for the INK4 locus deletion [N ' 0.5 (0.472– 0.556; SD, 0.013– 0.078)]. This is in contrast to healthy, haplotype-negative individuals who were shown to contain
two copies of the p15 gene per cell [N ' 1 (0.919 –1.019; SD, 0.006 – 0.075)]. Therefore, all results were clear-cut with relatively small SDs, rendering the deletional status unambiguous.
Discussion Since 1985, the identification of point mutations has been greatly facilitated by the development of PCR-based assays, but diagnosis of large gene rearrangements still awaited rapid and efficient technology. Recently, analysis of the loss of heterozygosity at a given locus by the use of highly polymorphic microsatellite markers allowed the characterization of a large number of gene deletions, especially in DNAs derived from tumor tissues. Nevertheless, this technology is limited by the presence of polymorphic markers and by their informativeness. In 1998, Cairns et al. (7 ) published a comparative study between microsatellite and quantitative PCR analyses to detect p16 copy number in primary bladder tumors and concluded that quantitative PCR, which is rapid, accurate, and sensitive, was certainly the most efficient method to identify genetic changes leading to tumorigenesis. In the aim to identify subjects with individual risk in a large French family diagnosed with a CMM-NST syndrome associated with a whole INK4 locus germ-line deletion, a combination of microsatellite markers of the 9p21 region, flanking the deletion and defining the at-risk haplotype was initially studied. However, this indirect test, which was limited by marker informativeness and required the analysis of nuclear families, possibly including both deceased parents and minor children, was not compatible with the present rules of individual consent
986
Laurendeau et al.: TaqMan Gene Dosage for Direct Predictive Testing
and secrecy. Thus, a TaqMan PCR-based gene dosage assay based on the quantification of the p15 gene was developed. This direct and accurate test allowed us to rapidly analyze all the family members willing to benefit from predictive testing and to identify those at risk. After an initial optimization step for adequate primer pairs and probes, the TaqMan gene dosage method can be used efficiently to assay any chromosome, gene, or exon and thus has a wide variety of applications in both clinical and research settings. This is especially the case with autosomal dominant conditions with haploinsufficiency, such as familial breast and ovarian cancer (8 ), neurofibromatosis type 1 (9 ), or in carrier detection in X-linked disorders such as Duchenne muscular dystrophy (10 ). We believe this high-throughput, fast-turnaround time (up to 96 samples/2 h, which corresponds to 12 individuals tested using our approach), and simple method is clearly a cost-effective alternative to fluorescence in situ hybridization in clinical applications. In this respect, TaqMan lends itself particularly well to routine diagnosis of common microdeletion syndromes. Because of its accuracy, this novel approach, more reliable than end-point PCR gene dosage (11 ) and not limited by short tandem repeat marker informativeness (12 ), also strengthens the potentialities of DNA diagnosis in the field of human aneuploidies.
This work was supported by the Association pour la Recherche sur le Cancer and the French Ministe`re de L’Enseignement Supe´rieur et de la Recherche.
References 1. Bahuau M, Vidaud D, Kujas M, Palangie´ A, Assouline B, Chaignaud-Lebreton M, et al. Familial aggregation of malignant melanoma/dysplastic naevi and tumours of the nervous system: an original syndrome of tumour proneness. Ann Genet 1997;40:78–91.
2. Bahuau M, Vidaud D, Jenkins RB, Bie`che I, Kimmel DW, Assouline B, et al. Germ-line deletion involving the INK4 locus in familial proneness to melanoma and nervous system tumors. Cancer Res 1998;58:2298 –303. 3. Holland PM, Abramson RD, Watson R, Gelfand DH. Detection of specific polymerase chain reaction product by utilizing the 59 to 39 exonuclease activity of Thermus aquaticus DNA polymerase. Proc Natl Acad Sci U S A 1991;88:7276 – 80. 4. Lie YS, Petropoulos CJ. Advances in quantitative PCR technology: 59 nuclease assays. Curr Opin Biotechnol 1998;9:43– 8. 5. Lo YMD, Tein MSC, Lau TK, Haines CJ, Leung TN, Poon PMK, et al. Quantitative analysis of fetal DNA in maternal plasma and serum: implications for noninvasive prenatal diagnosis. Am J Hum Genet 1998;62:768 –75. 6. Sambrook J, Fritsch EF, Maniatis T, eds. Molecular cloning: a laboratory manual. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory, 1989. 7. Cairns JP, Chiang PW, Ramamoorthy S, Kurnit DM, Sidransky D. A comparison between microsatellite and quantitative PCR analyses to detect frequent p16 copy number changes in primary bladder tumors. Clin Cancer Res 1998;4:441– 4. 8. Petrij-Bosch A, Peelen T, van Vliet M, van Eijk R, Olmer R, Drusedeau M, et al. BRCA1 genomic deletions are major founder mutations in Dutch breast cancer patients. Nat Genet 1997;17: 341–5. 9. Shen MH, Harper PS, Upadhyaya M. Molecular genetics of neurofibromatosis type 1 (NF1). J Med Genet 1996;33:2–17. 10. van Essen AJ, Kneppers AL, van der Hout AH, Scheffer H, Ginjaar IB, ten Kate LP, et al. The clinical and molecular genetic approach to Duchenne and Becker muscular dystrophy: an updated protocol. J Med Genet 1997;34:805–12. 11. Celi FS, Cohen MM, Antonarakis SE, Wertheimer E, Roth J, Shuldiner AR. Determination of gene dosage by a quantitative adaptation of the polymerase chain reaction (gd-PCR): rapid detection of deletions and duplications of gene sequences. Genomics 1994;21:304 –10. 12. Verma L, Macdonald F, Leedham P, McConachie M, Dhanjal S, Hulte´n M. Rapid and simple prenatal DNA diagnosis of Down’s syndrome. Lancet 1998;352:9 –12.