|
|
||||||||
Original Research |
National Farm Medicine Center, Marshfield Clinic Research Foundation, Marshfield, Wisconsin
National Farm Medicine Center, Marshfield Clinic Research Foundation, Marshfield, Wisconsin
National Farm Medicine Center, Marshfield Clinic Research Foundation, Marshfield, Wisconsin
REPRINT REQUESTS: Anne R. Greenlee, PhD, Marshfield Clinic Research Foundation, 1000 North Oak Avenue Mailstop: 2R5, Marshfield, WI 54449, Telephone: 715-389-4012, Fax: 715-389-3808, Email: greenlee.anne{at}mcrf.mfldclin.edu
Received: January 21, 2004.
Accepted: April 20, 2004.
| Abstract |
|---|
|
|
|---|
The purpose of this study was to establish a system for rapidly detecting single nucleotide polymorphisms (SNPs) in mitochondrial DNA (mtDNA) using hybridization probes and melting temperature (Tm) analysis. This technology should prove useful for population-based studies on the interaction between genetic factors and environmental exposures and the risk of Parkinson disease (PD).
METHODS
Mitochondrial DNA (mtDNA) was extracted from whole blood. Rapid polymerase chain reaction (PCR) and melting curve analyses were performed with primers and fluorochrome-labeled probes on a LightCycler (Roche Molecular Biochemical, Mannheim, Germany). Genotyping of 10 SNPs in 15 subjects was based on the analysis of allele-specific Tm of detection probes. The results of melting curve analyses were verified by sequencing all 150 PCR products.
RESULTS
Real-time monitoring showed optimal PCR amplification of each mtDNA fragment. The nucleotide changes at positions 1719, 4580, 7028, 8251, 9055, 10398, 12308, 13368, 13708, and 16391 from wild-type to mutant genotype resulted in 6.51, 8.29, 3.26, 7.82, 4.79, 2.84, 2.73, 9.04, 8.53, and 9.52°C declines in Tm of the detection probes, respectively. Genotyping of all 150 samples was verified by 100% correspondence with the results of sequencing. Fourteen subjects were haplogrouped by combining results for all 10 SNPs.
CONCLUSION
A rapid and reliable detection system for identifying mitochondrial polymorphisms and haplotypes was developed based on hybridization probe technology. This method may be suitable for mitochondrial genotyping of samples from large-scale epidemiology studies, and may prove useful for exploring the molecular etiopathogenesis of PD, identifying markers of genetic susceptibility, and protecting susceptible individuals from PD.
Key Words: Parkinson disease Mitochondrial DNA DNA Polymorphism DNA Hybridization Genotype Haplogroup
| INTRODUCTION |
|---|
|
|
|---|
Growing evidence suggests that the interaction of environmental factors810 and genetic predisposition may initiate the neurodegenerative process in PD.1113 One area of investigation has been the role of mitochondrial dysfunction in the pathogenesis of PD.14,15 Studies suggest pesticides and other environmental toxins inhibiting complex I of the electron transport system may be involved in the pathogenesis of PD.16 Indirect evidence suggests that, in some cases, the mitochondrial dysfunction in idiopathic PD may be due to mutations in mitochondrial DNA (mtDNA).1720 A single nucleotide polymorphism (SNP) leading to alterations in complex I activity may play a role in the pathogenesis of PD.21 The activity of complex I is reduced within brains of patients with PD.22 Whether enzyme deficiency is due to the presence of neurotoxins or to genetics, and how the interaction between environmental exposure and genetic factors influence the development of PD have not been established.
Restriction fragment length polymorphism (RFLP) studies of mtDNA from European Caucasoid, Native American, and African populations have revealed a number of stable polymorphic sites in mtDNA coding regions.2325 Related groups of polymorphic sites in mtDNAs are called haplogroups. Most mutations observed in mtDNA coding and control regions occur within these pre-existing haplogroups and define the individual mtDNA types or haplotypes. In comparison to individuals carrying the most common haplogroup, haplogroup H, individuals classified as haplogroup J and K demonstrate less risk of developing PD.21 This protective effect is strongly associated with the SNP 10398G that defines haplogroups J and K. The absence of 10398G characterizes haplogroup U and may be associated with an increased risk of PD.26 SNP 9055A may be associated with a protective effect for women.21
Population-based, case-control studies may assist with determining the interaction between host genetic factors (e.g., mitochondrial polymorphisms) and environmental exposures (e.g., pesticides, solvents and heavy metals) on the risk of PD.27 Towards this goal, we developed a method to detect mtDNA polymorphisms that combines real-time polymerase chain reaction (PCR) with hybridization probes on the LightCycler (Roche Molecular Biochemical, Mannheim, Germany). The hybridization technology selected is based on the detection of two adjacent oligonucleotide probes, whose fluorescent labels communicate through fluorescence resonance energy transfer (FRET).28,29 The anchor and sensor hybridization probes are added to the PCR and allowed to hybridize to an internal sequence of the PCR template. One probe is labeled at the 5'-end with the fluorophore LC Red 640, the other probe is labeled at the 3'-end with fluorescein. When these two probes hybridize to the PCR template, they are in close proximity, resulting in FRET between the two fluorophores. During FRET, the anchor fluorophore (fluorescein) is excited by the light source of the LightCycler instrument, and part of the excitation energy is transferred to the sensor fluorophore (LC Red 640). The LightCycler instrument monitors emitted fluorescence in real-time.
The fluorescence signal enables genotyping by analysis of the allele-specific melting behaviors of the sensor and anchor hybridization probes. The probe that is shorter has a lower melting temperature (Tm), spans the predicted polymorphism site and functions as a detection probe. During melting curve analysis, the temperature is slowly increased. If the sensor probe with a lower Tm is perfectly matched to its target sequence, it will melt away at a relatively high temperature. If this probe is mismatched, melting will occur at a lower temperature.30,31 When the two dyes are no longer in close proximity, the fluorescence signal decreases.
Here, we report on a rapid detection system using hybridization probes to genotype 10 SNPs in the human mitochondrial genome to define mtDNA haplogroups. This methodology should assist with epidemiologic studies designed to gain a better understanding of populations at risk of PD.
| MATERIALS AND METHODS |
|---|
|
|
|---|
97% Caucasian.32 The mtDNA was isolated from 100 µL of blood (mtDNA Extractor WB Kit, Wako Pure Chemical Industries, Ltd., Japan) and suspended in 10 mmol/L Tris (pH 7.4) containing 1 mmol/L EDTA. Any mtDNA not used immediately was stored at 4°C for no longer than 7 days or frozen at 20°C.
Design of primers and probes
Ten SNPs that define European mtDNA haplogroups located in sequences coding for genes and transcriptional control elements were chosen for genotyping according to the previous work of van der Walt et al. and Torroni et al.21,23 Haplogroups were characterized by combining ten substitutions at nucleotide positions (np) 1719, 4580, 7028, 8251, 9055, 10398, 12308, 13368, 13708, and 16391. These nps were based on the revised Cambridge mitochondrial genomic sequence, Genbank accession number J01415.
Table 1
lists the sequences for PCR primers and hybridization probes used for SNP detection. Primers and probes were designed with LightCycler Probe Design software, version 1.0 (Roche Molecular Biochemical), except for the probes recognizing polymorphisms at np 10398 and 8251. The probes for these two nps were designed by TIB Molbiol LLC (Adelphia, NJ). Probes to detect SNPs at nps 1719, 4580, 8251, 9055, 10398, 13368, 13708, and 16391 were complementary to the antisense (reverse) strands of the mitochondrial genome, whereas probes to detect SNPs at nps 7028 and 12308 were complementary to the sense (forward) strands. A probe to detect the nucleotide at np 7028 corresponded to the mutation, while the other nine pairs of probes corresponded to the wild-type sequences.
|
|
The LightCycler was programmed for pre-incubation at 95°C for 10 minutes (to activate FastStart Polymerase), amplification of target mtDNA (2847 cycles of denaturation at 95°C for 10 seconds, annealing at 55°C for 10 seconds, and extension at 72°C for 10 seconds), melting curve analysis, and cooling of the sample carousel. LC Red 640 was monitored on channel 2 (F2). The number of amplification cycles was 28 for SNP 10398; 29 for SNP 16391; 30 for SNPs 1719, 7028, 12308, and 13368; 31 for SNP 9055; 32 for SNP 4580; 34 for SNP 13708; and 47 for SNP 8251. The temperature transition rates were programmed at 20°C/second. Fluorescence was measured in a single step at the end of the annealing period of each cycle to monitor amplification. The number of amplification cycles was varied to ensure that each amplification curve would begin a "hook effect," which is defined as an exponential rise of the fluorescence signal followed by a drop in later cycles. Total run time ranged between 60 and 90 minutes and varied with the optimal number of amplification and melting cycles.
Melting analysis was performed following template amplification. Fluorescence was measured continuously during the slow temperature rise to monitor the dissociation of the shorter probe from the amplicon. The melting curve was initiated by cooling the reaction mixture to 45°C at 20°C/second, holding it at 45°C for 45 seconds, and then slowly heating it to 80°C at 0.1°C/second. Multiple melting cycles were conducted with the same amplified product to optimize the melting curve shape and to better discriminate between allele-specific Tm. One melting cycle was performed for SNP 1719; 2 cycles for SNP 8251; 4 cycles for SNPs 10398, 12308, and 13368; 5 cycles for SNP 9055; 6 cycles for SNPs 7028 and 16391; and 7 cycles for SNPs 4580 and 13708. Samples were held at 4°C upon completion of cycling. Fluorescence (F2) was plotted against temperature (T) to produce melting curves for each sample. Melting curves were then converted to melting peaks by plotting the negative derivative of the fluorescence with respect to temperature against temperature (d(F2)/dT versus T) using version 3.5 of the LightCycler® software.
Sequencing
Results of genotyping with hybridization probes were compared to PCR amplicon sequences to verify melting curve analyses. PCR products were purified using the QIAquick PCR purification kit (QIAGEN Inc., Valencia, CA). Purified PCR products were then sequenced using BigDye Terminator V3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA). The sequencing reaction contained 1.5 µL (1.5 ng) of purified PCR product, 4 µL of BigDye Cycle Sequencing Mix, 1.28 µL of primer (1.6 pmol), and 3.22 µL H2O. Thermal cycling conditions were 25 cycles of 96°C for 10 seconds, 50°C for 5 seconds, and 60°C for 4 minutes. Unincorporated primers and dye terminators were removed using DyeEx 96 plates (QIAGEN GmbH, Hilden, Germany). Purified samples were dried at 70°C and rehydrated with 10 µL of Hi-DiTM Formamide (Applied Biosystems). Samples were denatured at 95°C for 2 minutes and immediately chilled on ice for 2 minutes. Sequencing was performed on an ABI PRISM 3100 Genetic Analyzer (Applied Biosystems) with 3100 Genetic Analyzer Data Collection software, version 1.1. Sequencing analysis was performed using Lasergene SeqMan II software (DNASTAR, Inc., Madison, WI). Amplicon sequences were aligned with the revised Cambridge mtDNA genome sequence using BLAST 2 SEQUENCES software (http://www.ncbi.nlm.nih.gov/blast/bl2seq/bl2.html) to locate and identify base changes.
| RESULTS |
|---|
|
|
|---|
|
|
99% homology with mitochondrial genome sequence (100% identity for wild-type samples,
99% identity for mutant samples). Single base changes at the targeted SNP sites corresponded with results of melting curve analyses. As shown in table 3
|
| DISCUSSION |
|---|
|
|
|---|
We were able to classify 14 of 15 mtDNA samples into 8 of 9 haplogroups based on SNP combinations encompassing Finnish, Swedish, Tuscan and European mtDNA populations.21,23 Included in our study population were mtDNA haplogroups H, J, K, T, U, V, W and X. Haplogroup I was not identified in our samples which differed from the two previously cited studies. Haplogroup H, found most commonly in Scandinavian and European populations, represented 5 of 15 (33%) of our mtDNA study samples. Haplogroups K and W each were identified in 13% of our study samples. The haplogroups J, T, X, V, and U were each identified at a frequency of 7% in our study samples. One mtDNA sample (6%) was beyond the definition of haplogroups associated with European ancestry. This is consistent with two previous studies in which 2.4% of Caucasian samples23 and 8.4% of Caucasian samples21 were beyond these defined haplogroups.
Genotyping of SNPs by using fluorescent probes offers several advantages in comparison to the more common methodologies of RFLP analysis and DNA sequencing.2325,33,21 RFLP analysis relies on amplification of template, digestion of amplicon DNA with restriction enzymes, and subsequent sizing of the fragments by agarose gel electrophoresis. This approach is time consuming, labor intensive, and dependent on the presence of an appropriate restriction site and complete digestion, which cannot always be achieved. Some researchers employ only sequencing to detect mitochondrial mutations.34,35 However, it is difficult to apply sequencing to large sample sizes because of cost and complexities of sample preparation.
Another platform for fluorescent probe chemistry and genotyping SNPs is the TaqMan Allelic Discrimination Assay. TaqMan has shown promise for large-scale genotyping of mtDNA polymorphisms when combined with DNA extraction systems and multi-well plates.21 The assay uses the 5' nuclease activity of Taq polymerase to generate a fluorescent reporter signal during PCR extension. For SNP genotyping, one pair of TaqMan probes and one pair of PCR primers are used. The assay uses two TaqMan probes that differ at the polymorphic site, with one probe complementary to the wild-type allele and the other to the variant allele. Reporter and quencher dyes are covalently linked to the 5' and 3' ends of the wild-type and variant probes. When the probes are intact, fluorescence is quenched because of the physical proximity of the reporter and quencher dyes. During the PCR annealing step, the TaqMan probes hybridize to either the wild-type or polymorphic site. During the PCR extension phase, the 5' reporter dye is cleaved by the 5' nuclease activity of the Taq polymerase, leading to an increase in fluorescence of the reporter dye. Specific genotyping is determined by measuring the signal intensity of the reporter dye after the PCR reaction.36 Limitations of this approach include difficulties designing specific probes without knowing the sequence of all possible allelic variants, background fluorescence due to nonspecific binding between probes and target sequences, and an inability to measure the size of amplicons. Negative reactions may require sequencing for verification.
LightCycler technology minimizes these difficulties by monitoring shifts in melting temperatures of hybrids comprised of probe and wild-type or variant alleles. Fluorescence due to FRET occurs only when both probes anneal to the target sequence, reducing nonspecific background signals. Amplicon size is predicted by its direct relationship to Tm. Individual or multiple base changes due to insertions, deletions, and substitutions can be detected simultaneously. The principal drawback is the size of the carousel, which limits the number of samples to 32 per run.
Single base mismatches may lower the Tm of hybridization probes by several degrees Celsius.37 The amount of destabilization depends on the specific mismatch, neighboring nucleotides, and the position of the mutation within the sequence of the probe.38 Minor differences in Tm (<1°C) can be distinguished because of the reproducibility of the nucleic acid Tm. We achieved differences in Tm ranging between 2.73°C and 9.52°C, suggesting our anchor and sensor probes were properly designed for screening single nucleotide variations. Our standard deviations were <0.2°C, consistent with previously reported values.39
Typically, total DNA fractions are used to analyze mtDNA, although mtDNA occupies only a small portion of nucleic acid isolates.21,40 Isolation of mtDNA using the mtDNA Extractor WB Kit is efficient and enables simultaneous preparation of mtDNA from multiple samples for screening assays. We extracted approximately 1 µg mtDNA from 100 µL blood samples. Robust amplification of mtDNA sequences was achieved using nanogram quantities of template and real-time PCR. Nonspecific amplification of genomic DNA sequences was minimized using template preparations enriched with mtDNA.
Hybridization probe technology and melting curve analysis may have clinical, as well as research applications.39 The procedure has been multiplexed for the simultaneous detection of five polymorphisms modifying drug sensitivity and for SNPs in human cytochrome P450 enzymes involved in drug metabolism.41,42 Optimizing conditions for 10 different SNPs required changing the number of amplification or melting cycles without modifying reagent source or concentration. Necessary modifications were identified in real-time without setting up de novo reactions.
In summary, we have developed an assay system for rapidly genotyping and haplogrouping mitochondrial polymorphisms based on hybridization probe technology and melting curve analysis. Genotyping of 10 SNPs in 32 samples was accomplished in approximately 1 to 1.5 hours, making this method well suited for typing a variety of mitochondrial polymorphisms for association studies.
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
GRANT SUPPORT: Funded by National Institute for Occupational Safety and Health (NIOSH grant 5U50 OH007543-02)
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |