|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Molecular Diagnostics and Genetics |
1 Department of Pediatrics, Yong Loo Lin School of Medicine, National University of Singapore, Singapore; 2 Childrens Medical Institute, National University Hospital, Singapore; 3 Preimplantation Genetic Diagnosis Center, National University Hospital, Singapore; 4 Department of Pediatric Medicine, KK Womens and Childrens Hospital, Singapore; 5 Department of Biochemistry, Yong Loo Lin School of Medicine, National University of Singapore, Singapore, Singapore.
aAddress correspondence to this author at: Department of Pediatrics, Yong Loo Lin School of Medicine, National University of Singapore, Level 4, NUH, 5 Lower Kent Ridge Rd., Singapore 119074, Singapore. Fax 65-6779-7486; e-mail paecs{at}nus.edu.sg.
| Abstract |
|---|
|
|
|---|
Methods: We evaluated a novel, rapid microplate-based screening method for HD that combines the amplification and detection procedures in a single-step, closed-tube format. We carried out both the PCR for the HTT CAG-repeat region and the subsequent automated melting-curve analysis of the amplicon in the same wells on the plate. To establish cutoff melting temperatures (Tms) for each allelic class, we used a panel of reference DNA samples of known CAG-repeat sizes that represent a range of HTT alleles [normal (
26 repeats), intermediate (27–35 repeats), reduced penetrance expanded (36–39 repeats), and fully penetrant expanded (
40 repeats)]. We also measured well-to-well variation in Tm across the thermal block and validated cutoff Tms with DNA samples from 5 different populations. We also conducted a blinded validation analysis of clinical samples from an additional 40 HD-affected and 30 unaffected individuals.
Results: We observed a strong correlation between CAG-repeat size and amplicon Tm among the reference DNA samples. Use of the Tm cutoffs we established revealed that 5 samples from unaffected individuals had been misclassified as affected (1.1% false-positive rate). All samples from HD-affected and unaffected individuals were correctly identified in the blinded analysis.
Conclusions: This simple and scalable homogeneous assay may serve as a convenient, rapid, and accurate screen to detect the presence of pathologic expanded HD alleles in symptomatic patients.
| Introduction |
|---|
|
|
|---|
26 CAG repeats are stable and considered normal. Intermediate alleles contain 27–35 triplets, which may expand into the disease range during meiosis, especially in paternal alleles. Alleles with 36–39 CAG repeats are meiotically unstable, and although such alleles have been associated with the HD phenotype, they have reduced penetrance. Individuals carrying alleles with full penetrance (
40 CAG repeats) will display the HD pathology. Once in the disease range, repeat length is unstable during vertical transmission, with the instability appearing more pronounced in alleles of paternal origin(2)(3)(4). The age of HD onset is inversely correlated with CAG-repeat length, although other factors also play a modulating role(5)(6). A few methods for the molecular diagnosis of HD have been proposed(7)(8)(9)(10)(11). The most common method for accurately sizing CAG repeats uses fluorescently tagged PCR primers that anneal to sequences flanking the CAG-repeat region and generate amplification products that can then be detected with automated genotyping systems(11)(12). The separate amplification and detection methods involved in these methods thus require the transfer of samples after the PCR. More recently, Gundry and Wittwer demonstrated that amplicons containing normal and expanded numbers of CAG repeats can be distinguished in the presence of SYBR Green I dye on the basis of their slightly different melting temperatures (Tms)(13). We used a 96-well plate format with a real-time PCR instrument to assess the utility of melting curve analysis (MCA) for distinguishing HTT allele classes. We used both automated capillary electrophoresis and MCA to analyze samples of known genotypes that represent a wide spectrum of CAG-repeat alleles and established Tm cutoffs for differentiating allelic classes. The assay was validated both with assumed wild-type DNA samples from 5 population groups and in a blinded analysis of clinical samples from HD-affected and unaffected individuals.
| Materials and Methods |
|---|
|
|
|---|
pcr and mca
All PCR-MCA assays were performed on a LightCycler® 480 instrument (Roche Applied Science). We designed a primer pair, HD-F (5'-CCTTCGAGTCCCTCAAGTCCTT-3') and HD-R (5'-GCGGTGGCGGCTGTTG-3'), to amplify only across the CAG-repeat stretch of the HTT gene, and not the adjacent downstream CCG repeats, and took into account previously reported single-nucleotide polymorphisms in the region(14)(15). Each well on a 96-well reaction plate contained a 20-µL PCR reaction mixture consisting of 0.4 µmol/L of each primer, 50 ng genomic DNA, 1x Q-Solution (Qiagen; provided with the Qiagen HotStarTaq polymerase at 5x concentration. Components are not disclosed in the product inserts.), and 1x LightCycler 480 SYBR Green I Master [Roche Applied Science; components of the LightCycler 480 SYBR Green I Master are FastStart Taq DNA polymerase, reaction buffer, dNTP mix (with dUTP instead of dTTP), SYBR Green I dye, and MgCl2; the concentrations of the components are not listed.]. An initial polymerase-activation step at 95 °C for 10 min was followed by 40 cycles of 98 °C for 45 s, 70 °C for 30 s, and 72 °C for 45 s. The PCR reaction was followed automatically with a melting program consisting of denaturation at 95 °C for 1 min, a temperature-hold step at 60 °C for 1 min, and then slow temperature ramping from 60 °C to 95 °C over approximately 27 min. We acquired 50 readings of fluorescence intensity per degree, a data-acquisition rate that we derived empirically. The same thermocycling and post-PCR melting conditions were used for amplification and MCA on a 384-well platform, with each well containing 5 µL of the PCR reaction mixture.
We used LightCycler 480 software to collect fluorescence intensities and convert the data into melting-curve and melting-peak plots. The Tm for each PCR product corresponds to the temperature at the center of the melting peak.
Selected DNA samples were used to assess run-to-run reproducibility and well-to-well uniformity in Tm. PCR-MCA runs were replicated for each DNA sample. For each run, identical aliquots of the PCR reaction mix (consisting of the components mentioned above: primer, Q-Solution, and LightCycler 480 SYBR Green I Master) containing the DNA sample were loaded into all 96 wells of the PCR reaction plate.
pcr and GENESCAN analysis
Genomic DNA was amplified in a final volume of 25 µL. Each reaction contained 0.4 µmol/L of each primer, 0.25 mmol/L of each deoxynucleoside triphosphate, 1.25 U of HotStarTaqTM DNA polymerase (Qiagen), 1x Q-Solution, 1x of the supplied PCR buffer (including 1.5 mmol/L MgCl2; Qiagen), and 50 ng genomic DNA. We used the primers described above, with the exception that we fluorescently labeled HD-R at the 5' end with HEX (6-carboxyhexachlorofluorescein). An initial denaturation step at 95 °C for 15 min was followed by 30 cycles of 98 °C for 45 s, 63 °C for 60 s, and 72 °C for 90 s.
One microliter of each PCR product was mixed with 0.3 µL of the GeneScan –500 ROXTM size calibrator and 9 µL of HiDiTM Formamide (Applied Biosystems). The mixtures were subjected to capillary electrophoresis on a 3130XL Genetic Analyzer (Applied Biosystems) for 50 min, and the electropherograms were analyzed with GeneMapper® software (version 4.0; Applied Biosystems).
statistical analysis
The Statistical Package for Social Sciences (version 15.0 for Windows; SPSS) and Microsoft Excel were used for ROC curve and correlation analyses, respectively. To determine cutoff Tms between allele classes, we plotted ROC curves by using the sizes of the CAG repeats (obtained by GeneScan analysis) and the Tms of the respective PCR products (obtained by PCR-MCA) as the test variables and using allele class (normal, intermediate, expanded with incomplete penetrance, or expanded with complete penetrance) as the classification variable.
| Results |
|---|
|
|
|---|
|
We noted that NA20206 displayed only a single Tm peak for one normal allele, although this individual is heterozygous for 2 normal alleles that differ by 1 CAG repeat (Fig. 1B
). Similarly, both HD-affected individuals (NA20208 and GM05542) showed only a single Tm peak for one expanded allele, although the alleles of these individuals differ by 11 repeats and 5 repeats, respectively (Fig. 1, G and H
). The rest of the HD-affected individuals each carried 1 normal and 1 expanded allele and showed 2 Tm peaks, one within the normal range of Tms and one within the expanded range (see Table 1 in the online Data Supplement; Fig. 1). The largest normal allele in the reference DNA series, a 23-repeat allele in NA20246, had a Tm of 87.00 °C, whereas the smallest expanded allele in the series, a 36-repeat allele in NA20248, had a Tm of 88.37 °C. An intermediate allele of 30 CAG repeats (in CD00022) had a Tm of 88.03 °C. Thus, a modest but distinct separation of >0.3 °C existed between the largest nonexpanded HD allele and the smallest expanded allele that we tested. These results suggest that PCR-MCA can distinguish most HD alleles with normal, intermediate, and expanded numbers of repeats and, more importantly, can clearly distinguish between pathogenic expanded HD alleles and nonpathogenic normal and intermediate HD alleles; however, the results also reveal that this system cannot distinguish normal alleles that differ by small numbers of repeats or distinguish large alleles that differ by as many as 11 repeats when both normal alleles or both expanded alleles are present in the same individual.
reproducibility of the pcr-mca TM
We used DNA from GM04282 (16 and 74 repeats) to assess run-to-run reproducibility and well-to-well uniformity of the Tm value and performed replicate PCR-MCA runs with the same DNA sample. For each run, we loaded identical aliquots of the PCR reaction mix containing GM04282 DNA into all 96 wells of the microplate. The PCR was followed immediately by MCA in the same well. The Tm peaks obtained remained generally similar across both runs, with 0.76 °C being the largest difference between runs (Fig. 2
, top row, left and center panels). We observed interwell variation in Tm peaks within each run, with the Tms of samples placed in the middle of the plate being approximately 0.6 °C lower than the Tms of samples at the plate edges. The variation in Tm was greatest across a row and less down a column (Fig. 2
, middle row, left and center panels). We observed similar trends in interwell variation across the plate for PCR-MCAs performed on a 384-well platform, but the Tm variance among wells increased slightly (SD = 0.2 °C; Fig. 2
, right panels).
|
To evaluate the risk of misdiagnosing expanded and intermediate alleles as normal alleles caused by interwell variation in Tm, we performed similar runs with NA20248 and NA20250 DNA samples as templates and determined the SD for Tms among wells to be approximately 0.2 °C for both the 36-repeat and 39-repeat alleles in NA20248 and NA20250, respectively (Fig. 3
). Similarly, for DNA from an anonymous cord blood sample with HD alleles of 16 and 27 repeats (based on GeneScan analysis), we measured the Tm SD for the 27-repeat allele across a plate. The limited amount of this DNA sample allowed us to test only 14 well positions distributed evenly across the 96-well plate, and we obtained a mean (SD) Tm of 87.8 °C (0.3 °C) for the 27-repeat allele (Fig. 3
).
|
determination of TM cutoffs separating normal, intermediate, and expanded allele classes
We conducted an ROC curve analysis of the results obtained with DNA samples supplied by Coriell Cell Repositories to identify the optimal Tm cutoff values separating the 4 allele classes and obtained Tm cutoffs of 87.3 °C for the boundary between normal and intermediate alleles, 88.2 °C for the boundary between intermediate alleles and expanded alleles with reduced penetrance, and 88.5 °C for the boundary between expanded alleles with reduced penetrance and expanded alleles with full penetrance (Fig. 1
).
None of the Tms for the 36-repeat and 39-repeat expanded alleles of NA20248 and NA20250 (Fig. 3
) fell below the calculated Tm cutoff of 88.2 °C separating intermediate alleles from expanded alleles, regardless of well position. Similarly, none of the Tms for the 27-repeat intermediate allele from the anonymous DNA sample (Fig. 3
) fell below the calculated Tm cutoff of 87.3 °C separating normal alleles from intermediate alleles. Thus, there was minimal risk of mistaking an intermediate or expanded allele for a nonexpanded allele due to between-well variation in Tm.
validation of pcr-mca screening for hd
To test the feasibility of PCR-MCA screening for HD with the 96-well plate format, we genotyped DNA samples from anonymous donors from 5 different populations (Caucasian, African American, Chinese, Malay, and Indian). In parallel, we precisely sized the CAG-repeat region of each sample via the PCR and GeneScan analysis (see Figs. 1–5 in the online Data Supplement). When we excluded samples in which the melting peaks of 2 alleles of different sizes were too close and merged to form a single peak (producing a single Tm that was the average of those of the 2 unresolved melting peaks), we found CAG-repeat size and amplicon Tm to be strongly correlated in all 5 populations (Fig. 4
). Only individuals with 2 normal alleles that differed by
6 repeats generated 2 distinct Tm peaks with this PCR-MCA assay. Including the Tms of merged peaks produced greater within-population variation and poorer correlation between CAG-repeat size and amplicon Tm (Fig. 4
). Of the 471 samples from unaffected individuals that we screened, 20 samples (4%) had the Tm of the larger allele exceed their allele class cutoffs. Of these 20 samples, only 5 of the samples from unaffected individuals were incorrectly called as being from HD-affected individuals (1% false-positive rate). This result was due to the presence of an intermediate allele with a Tm
88.2 °C (3 Caucasians and 2 African Americans with CAG-repeat lengths of 27, 29, 31, and 34) (Fig. 4
; see Figs. 1 and 2 in the online Data Supplement). The remaining 15 samples had a large normal allele (23–26 repeats) with a Tm
87.3 °C, the Tm cutoff separating the normal and intermediate allele classes. None of the normal alleles had Tms
88.2 °C, the Tm cutoff separating the intermediate and expanded allele classes. Of the 11 samples carrying intermediate alleles, only 1 sample with a 27-repeat allele had a Tm that fell below the 87.3 °C cutoff separating normal and intermediate alleles. There were no expanded alleles in the populations used in this validation study.
|
To test the accuracy of the PCR-MCA assay for distinguishing between samples from HD-affected and unaffected individuals, we performed an additional blinded analysis of 70 clinical samples, 40 from HD-affected individuals and 30 from unaffected individuals. The analysis was performed twice. In the first run, we placed samples randomly, starting from the center of the block and moving toward the periphery. In the second run, samples were placed randomly, beginning at the periphery and moving toward the center of the block (see Figs. 6 and 7 in the online Data Supplement). The results showed that this assay and the use of an 88.2 °C Tm cutoff allowed HD-affected individuals to be distinguished from unaffected individuals with 100% sensitivity and specificity in both of the runs (see Table 2 in the online Data Supplement). No expanded alleles were misdiagnosed as normal. Only 1 sample from an HD-unaffected individual with an intermediate allele (23 and 29 repeats) was classified as having normal alleles only. This result was due to merging of the melting peaks of the large normal allele and the small intermediate allele for this sample.
| Discussion |
|---|
|
|
|---|
The results of the validation study indicate that the PCR-MCA assay is able to distinguish between most normal and intermediate alleles and can differentiate all normal and expanded HD alleles with no dropout of any of the expanded alleles tested, up to 99 CAG repeats. The PCR-MCA assay appears to be more robust than capillary electrophoresis for detecting large expanded alleles in the presence of a normal allele, possibly because of the tendency of the SYBR Green I dye to migrate from the melted shorter amplicon and intercalate with the larger amplicon (Fig. 1
, D–F and I).
The assay was able to distinguish different normal alleles in the same individual that differ by
6 CAG repeats. Normal alleles that differ by 5 or fewer CAG repeats have similar Tms that cause the melting peaks to merge and form a single but slightly broader peak. Similarly, the melting peaks of 2 different expanded alleles in homozygous affected individuals could not be separated when such alleles differed by up to as many as 10 CAG repeats. However, this inability to distinguish clearly between 2 alleles in the fully penetrant range of expansion should not lead to misdiagnosis or a false-negative result.
CAG-repeat length and amplicon Tm were strongly correlated among the DNA samples from Coriell Cell Repositories (Fig. 4
). Normal alleles had Tms
87.00 °C, intermediate alleles had Tms between 87.73 °C and 88.03 °C, and expanded alleles had Tms
88.37 °C. We used these results to establish Tm cutoffs separating the different allele classes and validated these cutoffs both in an analysis of DNA samples from anonymous donors from 5 populations and in a blinded analysis of DNA from clinical samples.
None of the normal alleles were classified as expanded alleles among the 471 samples from the 5 populations. Fifteen large normal alleles (23–26 repeats) were misclassified as intermediate alleles, and 5 intermediate alleles were misclassified as expanded alleles with incomplete penetrance (see Fig. 4
and Figs. 1–5 in the online Data Supplement for details). The single misclassification of an intermediate allele (27 repeats) as normal was due to the co-presence of a normal allele that was 5 repeats shorter than the intermediate allele, which produced merged melting peaks with an apparent Tm that was the average of the Tms for the 2 alleles.
In the blinded validation analysis, samples were randomly distributed across the 96-well plate to determine if well position and DNA-extraction method affected the diagnostic accuracy of this assay (see Table 2 and Figs. 6 and 7 in the online Data Supplement). All samples from HD-affected individuals were detected with 100% sensitivity, regardless of DNA-extraction method and well position.
The observation of interwell variation in Tm for the same sample reflects an inherent property of all block-based regular and real-time thermal-cycling instruments—uneven heating across the reaction plate. Therefore, block-based instruments may not be ideal for precisely classifying alleles with the Tm cutoffs we have defined; however, the assay does represent an effective and accurate method for screening patients with clinical symptoms of HD, i.e., to differentiate those who carry a disease-causing expanded allele (and thus are truly affected with this disorder) from those who display HD-like symptoms but are unaffected by HD. Nonetheless, including control samples of known CAG-repeat sizes at strategic locations around the plate in the same run is necessary to detect appreciable instrument-related Tm shifts in the heating block over time. When available, we recommend that positive-control samples carry alleles with 27 and 36 repeats, which are the smallest repeat sizes for intermediate and expanded alleles, respectively.
Recent studies have shown that the reproducibility of results of single-nucleotide polymorphism genotyping across a reaction plate with high-resolution melting analysis, a recent derivation of standard amplicon melting analysis, can be greatly improved by including 2 temperature calibrators that melt above and below the expected Tm of the amplicon(17)(18). Well-to-well variation in Tm can be corrected by aligning the melting peaks of the temperature calibrators for all the samples in a reaction plate. Such temperature calibrators could potentially be applied to HD PCR-MCA to correct for well-to-well variation in Tm across a thermal block; however, unlike the single invariant amplicon size involved in genotyping of single-nucleotide polymorphisms with high-resolution melting analysis, appropriate temperature calibrators would have to be developed to bracket the entire Tm range of CAG-repeat amplicon sizes in HD. This goal may be difficult to achieve, given the 99 °C thermal limit of the LightCycler 480 block and the need for the upper limit of an amplicon Tm to be several degrees below the thermal limit because of the requirement that the fluorescence be at baseline before the final temperature cutoff. This consideration is further complicated by the fact that the Tm of a 99-repeat allele is already 90.88 °C, and the Tms of larger alleles are unknown. As with existing PCR-based diagnostic assays for HD, there is a risk of failure to detect the hyperexpanded allele with the PCR-MCA method. Hyperexpanded alleles are extremely rare and usually present atypical phenotypes(19)(20). If such hyperexpansion is suspected, because of a very early age of disease onset for example, we recommend that a negative test or screening result be followed up with Southern blotting and pedigree analysis.
Despite some limitations, the most distinctive advantage of this assay over existing molecular diagnostic assays for HD involving single-sample, capillary-based approaches is its technical simplicity and potential for scaled-up screening. This simple highly automated assay requires no probe or fluorescently tagged primer (except for CAG-sizing follow-up) and no sample transfer. Moreover, the single-step, closed-tube format reduces the potential for the sample mix-ups inherent in multistep analyses. PCR-MCA is thus attractive as a rapid first-line screening assay for HD because size confirmation will be required only for the small fraction of samples that test positive, as indicated by the Tm values. Any false-positive calls will be detected at this stage and will not lead to misdiagnosis. The results of preliminary trials with the new generation of DNA-saturating dyes (such as those used in high-resolution melting analysis) in the place of SYBR Green I suggest that better Tm (and allele) resolution can be achieved in some situations (data not shown), and these findings warrant further evaluations. The potential of adapting PCR-MCA to the molecular diagnosis of other trinucleotide-repeat disorders is also worth exploring.
| Acknowledgments |
|---|
Financial Disclosures: None declared.
Acknowledgments: We thank Shen Liang, Biostatistics Unit, Yong Loo Lin School of Medicine, for assistance with statistical analyses.
| Footnotes |
|---|
2 Human genes: HTT, huntingtin. ![]()
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |