|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Molecular Diagnostics and Genetics |
1 Institute for Molecular and Human Genetics, Georgetown University Medical Center, Washington, DC.
aAddress correspondence to this author at: Institute for Molecular and Human Genetics, Georgetown University Medical Center, M4000, 3800 Reservoir Rd. NW, Washington, DC 20007. Fax 202-444-1770; e-mail wonglj{at}georgetown.edu.
| Abstract |
|---|
|
|
|---|
Methods: To reliably quantify low proportions of the mutant mtDNA, we developed a real-time amplification refractory mutation system quantitative PCR (ARMS-qPCR) assay. We validated the method with experimental samples containing known proportions of mutant A3243G mtDNA generated by mixing known amounts of cloned plasmid DNA containing either the wild-type or the mutant sequences.
Results: A correlation coefficient of 0.9995 between the expected and observed values for the proportions of mutant A3243G in the experimental samples was found. Evaluation of a total of 36 patient DNA samples demonstrated consistent results between PCRrestriction fragment length polymorphism (RFLP) analysis and real-time ARMS-qPCR. However, the latter method was much more sensitive for detecting low percentages of mutant heteroplasmy. Three samples contained allele-specific oligonucleotide-detectable but RFLP-undetectable mutations.
Conclusions: The real-time ARMS-qPCR method provides rapid, reliable, one-step quantitative detection of heteroplasmic mutant mtDNA.
| Introduction |
|---|
|
|
|---|
Numerous nucleotide substitutions responsible for respiratory chain defects have been identified (2)(5)(6)(7). The most common point mutation is the A3243G mutation in tRNA leu(UUR), which is responsible for mitochondrial encephalopathy lactic acidosis, and stroke-like episodes (MELAS) syndrome and maternally inherited diabetes and deafness (MIDD). The proportion of mutant mitochondria must reach the threshold to cause clinical phenotype. Thus, to predict a patients clinical outcome, it is important to determine the extent of mutant load in the affected or relevant tissue. In the case of MIDD, patients usually carry low percentages of mutant mtDNA in the blood that may not be detected by conventional restriction fragment length polymorphism (RFLP) analysis. This becomes a problem, especially when the determination of carrier status is important to confirm the diagnosis of a maternally inherited mitochondrial disorder and to provide appropriate genetic counseling. The detection limit for RFLPethidium bromide (EtBr) gel analysis is
510%, and that for radioactive allele-specific oligonucleotide (ASO) dot-blot analysis is
2% (8)(9). However, ASO analysis alone cannot be used to determine the percentage of mutant heteroplasmy.
We recently developed a real-time quantitative PCR (RT-qPCR) that uses allele-specific TaqMan probes to measure the percentage of heteroplasmic A3243G mutation (10), Although the results were consistent with the results from RFLP analysis, because of nonspecific binding of the mutant TaqMan probe to wild-type target DNA sequences and vice versa, the heteroplasmy showed a significant background (unpublished observation). Thus, it was not of value in the measurement of low proportions of mutant heteroplasmy.
To measure low proportions of mutant heteroplasmy, we designed mismatched primers for an amplification refractory mutation system (ARMS) assay. The introduction of one or two mismatched nucleotides immediately 5' to the mutation site greatly increased the binding specificity of the allele-specific modified primers toward either the wild-type or the mutant sequence targets (11). Here we report the validation of the RT ARMS-qPCR assay for the detection and quantification of known DNA samples with the A3243G mutation. The results demonstrate that the RT ARMS-qPCR assay is a rapid, sensitive, reliable, and cost-effective method with potential for use in the detection and quantification of heteroplasmic mtDNA mutants.
| Materials and Methods |
|---|
|
|
|---|
primers for arms
In ARMS PCR, the reverse primer (5'-TGGCCATGGGTATGTTGTTA-3') was at nucleotide positions 33193300 for the amplification of both wild-type mtDNA and the A3243G mutation. The forward primers for RT ARMS-PCR analysis of the A3243G mutation are listed in Table 1
. The numbers correspond to the nucleotide positions in GenBank accession no. NC001807 and the MITOMAP database (2). Originally, a mismatch at the penultimate nucleotide position of the mutation site was introduced in the forward primers (Table 1
, set A) to increase the specificity of the ARMS reaction based on the principles described by Newton et al. (11). Nevertheless, the specificity was not sufficient, and there were background signals at 0% mutant. Two mismatches at the two nucleotides immediately 5' to the mutation site were then introduced (Table 1
, set B), which greatly improved the specificity. Thus, modified forward primers from set B were used for this study. At the mutation site, the wild-type primer containing an "A" would perfectly match the wild-type target sequence but would be a weak "AC" mismatch with the mutant target sequence. Similarly, the mutant primer containing a "G" would be a perfect match with the mutant target sequence but would be a weak "GT" mismatch with the wild-type target sequence. Therefore, introduction of a strong CC mismatch at the penultimate nucleotide is expected to increase the primer specificity (11).
|
RT-QPCR
The RT ARMS-qPCR assay was performed in triplicate for each wild-type and mutant target sequence. The 20-µL PCR reaction mixture contained 1x Platinum SYBR Green qPCR SuperMix UDG (Invitrogen), 500 nM each primer, Rox dye, and
4 ng of total genomic DNA extract, according to the manufactures instructions. Real-time PCR conditions were 2 min at 50 °C and 10 min at 95 °C, followed by 45 cycles of denaturation for 15 s at 95 °C and annealing/extension for 60 s at 63 °C. SYBR Green dye binds to the minor groove of double-stranded DNA and increases the intensity of the fluorescent emissions while the amplicons are produced in each amplification cycle. The fluorescent signal intensities were recorded and analyzed during PCR in an ABI Prism 7700 sequence detector system (Applied Biosystems) using the SDS (Ver. 1.91) software. Dissociation curves for the amplicons were generated after each run to confirm that the increased fluorescence intensities were not attributable to nonspecific signals (primer-dimers). The increase in fluorescent signal is associated with an exponential growth of PCR product during the linear-log phase. The threshold cycle (CT) is the cycle at which a significant increase in the reaction product is first detected. The higher the initial amount of DNA, the sooner accumulated product is detected in the PCR process and the lower the CT value. Thus, the CT values within the linear exponential increase phase are used to measure the original DNA template copy numbers and to construct the calibration curve. If a sample contained >100 000 or <100 copies based on the CT, the assay was repeated at a higher or lower dilution of the DNA extract so that the measurement would fall within a linear DNA copy number range.
preparation of dna for calibration curves
DNA for the wild-type and mutant target sequences was generated from cloned plasmid DNA containing pCR2.1-TOPO vector (Invitrogen) and PCR products of primers mtF3212 and mtR3471. The copy numbers of the wild-type and mutant DNA sequences were calculated based on the size and molecular weight of the plasmid DNA. Serial dilutions were made, and RT ARMS-qPCR reactions were performed to construct the calibration curves for the wild-type and mutant A3243G DNA sequences (Fig. 1
).
|
measurement of mutant heteroplasmy
Calibration curves for both wild-type and mutant mtDNA were always included in each run. The copy number of the target sequence in the sample was calculated from the CT number and the calibration curve. The proportion of mutant A3243G sequence was calculated from the copy number of the wild-type and mutant sequences. Alternatively, the proportion of the mutant mtDNA could be calculated from
CT (CTwild-type CTmutant), using the formula: Proportion of mutant = 1/(1 + 1/2
CT).
rflp analysis with agarose gel electrophoresisETBRstaining
For RFLP analysis, the primers used were mtF3085 and mtR3406, which produced a PCR product of 322 bp. In the presence of A3243G mutation, the PCR product was cleaved by ApaI to two fragments of 159 and 163 bp (14). The ApaI-digested PCR products were separated on a 3% agarose gel. Free software, Scion Image for Windows (Beta 4.02;http://www.scioncorp.com), was used to analyze the ratio of the signal intensity for the 159/163-bp band to the sum of the intensities for the 322-bp and 159/163-bp bands. RFLP analysis was also performed by electrophoresis of HaeIII digest of the PCR products on 5% polyacrylamide gels and calculation of ratios by the same software (8).
| Results |
|---|
|
|
|---|
For RT ARMS-qPCR analysis, calibration curves for the wild-type A3243 and mutant 3243G sequences were constructed. There was excellent correlation between the cycle number and mtDNA copy number with correlation coefficients of 1.000 and 0.999 for the wild-type and mutant target sequences, respectively. To validate the method, we mixed known amounts of the wild-type and mutant mtDNA to obtain 13 samples containing various proportions of mutant mtDNA. Results of the RT ARMS-qPCR are shown in Fig. 2
. The correlation coefficient (R2) for the observed vs expected proportion of mutant was 0.9995. A proportion of mutant mtDNA in samples as low as 0.1% was detected by this method.
|
To assess the potential usefulness of this method for the detection of unknown mutant A3243G heteroplasmy, we reevaluated the samples that had been analyzed by both the ASO and RFLP methods. We obtained mutant heteroplasmy proportions of 0.5%, 3%, and 4% for samples 29, 17, and 20, respectively, which had not been detected by RFLP with EtBr staining. For samples 3032, which were negative by both radioactive ASO and RFLPEtBr analysis, the qPCR analysis gave a value of 0.03%, close to the mean (SD) value of 0.03 (0.003)% for 12 wild-type controls. When we compared the results obtained by qPCR analysis with the results obtained by two independent RFLP analyses (Fig. 3
), the correlation coefficient was 0.935.
|
| Discussion |
|---|
|
|
|---|
Use of RT ARMS-qPCR for the measurement of mutant mtDNA heteroplasmy has several advantages. Generally, target sequences are amplified, and sequence-specific TaqMan (16) or Molecular Beacon (17) probes for the mutant or wild type are used for detection and quantification. These methods require the synthesis of each specific labeled probe, which can be costly. Our approach is to design target-sequence-specific ARMS primers for specific amplification of the target sequence followed by SYBR Green binding and reporting of the PCR products. It therefore eliminates the need to synthesize expensive fluorescent probes and can be easily expanded for the analysis of other common mtDNA mutations, such as A8344G for the MERRF (myoclonic epilepsy and ragged red fiber) syndrome, T8993G and T8993C for NARP (neuropathy, ataxia, and retinitis pigmentosa) syndrome, and the detection of homoplasmic LHON (Leber hereditary optic neuropathy) mutations. It not only detects and quantifies the heteroplasmic mtDNA mutation in a single step, but also makes it possible to determine the total amount of mtDNA genome if a standard nuclear gene target sequence is included in the analysis for normalization (17). This is particularly useful if mitochondrial proliferation is being considered, as in the case of the A8344G mutation (17). The method could allow simultaneous quantification of total mtDNA genome content in addition to detection and quantification of point mutations. Using this method, we have recently demonstrated the compensatory amplification of mtDNA in a mildly affected patient who harbored 92% deletion mutant mtDNA (18). RT-qPCR has been widely used to monitor mtDNA depletion in HIV-infected patients receiving prolonged nucleoside analog reverse transcriptase inhibitor therapy (10)(19). It has also been used to quantify mtDNA deletions (18)(20). We believe that RT ARMS-qPCR could ultimately provide a one-step technique for both qualitative and quantitative molecular analysis of mtDNA defects.
In conclusion, we present a rapid, sensitive, and reliable one-step RT-qPCR method for quantification of mutant mtDNA heteroplasmies. This method could be used for the quantification of heteroplasmy for any given mtDNA point mutation.
| Acknowledgments |
|---|
| Footnotes |
|---|
| References |
|---|
|
|
|---|
The following articles in journals at HighWire Press have cited this article:
![]() |
S M Ware, N El-Hassan, S G Kahler, Q Zhang, Y-W, E Miller, B Wong, R L Spicer, W J Craigen, B A Kozel, et al. Infantile cardiomyopathy caused by a mutation in the overlapping region of mitochondrial ATPase 6 and 8 genes J. Med. Genet., May 1, 2009; 46(5): 308 - 314. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. K. Yoon, A. Roorda, Y. Zhang, C. Nakanishi, L.-J. C. Wong, Q. Zhang, L. Gillum, A. Green, and J. L. Duncan Adaptive Optics Scanning Laser Ophthalmoscopy Images in a Family with the Mitochondrial DNA T8993C Mutation Invest. Ophthalmol. Vis. Sci., April 1, 2009; 50(4): 1838 - 1847. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Mochizuki, C. Sugimori, Z. Qi, X. Lu, A. Takami, K. Ishiyama, Y. Kondo, H. Yamazaki, H. Okumura, and S. Nakao Expansion of donor-derived hematopoietic stem cells with PIGA mutation associated with late graft failure after allogeneic stem cell transplantation Blood, September 1, 2008; 112(5): 2160 - 2162. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. S. Jahangir Tafrechi, F. M. van de Rijke, A. Allallou, C. Larsson, W. C.R. Sloos, M. van de Sande, C. Wahlby, G. M.C. Janssen, and A. K. Raap Single-cell A3243G Mitochondrial DNA Mutation Load Assays for Segregation Analysis J. Histochem. Cytochem., November 1, 2007; 55(11): 1159 - 1166. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Genasetti, M. L. Valentino, V. Carelli, D. Vigetti, M. Viola, E. G. Karousou, G. V. Melzi d'Eril, G. De Luca, A. Passi, and F. Pallotti Assessing Heteroplasmic Load in Leber's Hereditary Optic Neuropathy Mutation 3460G->A/MT-ND1 with A Real-Time PCR Quantitative Approach J. Mol. Diagn., September 1, 2007; 9(4): 538 - 545. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. S. Lim, R. K. Naviaux, and R. H. Haas Quantitative Mitochondrial DNA Mutation Analysis by Denaturing HPLC Clin. Chem., June 1, 2007; 53(6): 1046 - 1052. [Abstract] [Full Text] [PDF] |
||||
![]() |
L-J C Wong, D Yim, R-K Bai, H Kwon, M M Vacek, J Zane, C L Hoppel, and D S Kerr A novel mutation in the mitochondrial tRNASer(AGY) gene associated with mitochondrial myopathy, encephalopathy, and complex I deficiency. J. Med. Genet., September 1, 2006; 43(9): e46 - e46. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Fan, C. Civalier, J. K. Booker, M. L. Gulley, T. W. Prior, and R. A. Farber Detection of Common Disease-Causing Mutations in Mitochondrial DNA (Mitochondrial Encephalomyopathy, Lactic Acidosis with Stroke-Like Episodes MTTL1 3243 A>G and Myoclonic Epilepsy Associated with Ragged-Red Fibers MTTK 8344A>G) by Real-Time Polymerase Chain Reaction J. Mol. Diagn., May 1, 2006; 8(2): 277 - 281. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Waku-Kouomou, A. Alla, B. Blanquier, D. Jeantet, H. Caidi, A. Rguig, F. Freymuth, and F. T. Wild Genotyping Measles Virus by Real-Time Amplification Refractory Mutation System PCR Represents a Rapid Approach for Measles Outbreak Investigations J. Clin. Microbiol., February 1, 2006; 44(2): 487 - 494. [Abstract] [Full Text] [PDF] |
||||
![]() |
R.-K. Bai and L.-J. C. Wong Simultaneous Detection and Quantification of Mitochondrial DNA Deletion(s), Depletion, and Over-Replication in Patients with Mitochondrial Disease J. Mol. Diagn., November 1, 2005; 7(5): 613 - 622. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |