Rapid and enhanced detection of mitochondrial DNA variation using single-strand conformation analysis of superposed restriction enzyme fragments from polymerase chain reaction-amplified products

A reduced number of mtSNPs saturates mitochondrial DNA haplotype diversity of worldwide population groups

Background

The high levels of variation characterising the mitochondrial DNA (mtDNA) molecule are due ultimately to its high average mutation rate; moreover, mtDNA variation is deeply structured in different populations and ethnic groups. There is growing interest in selecting a reduced number of mtDNA single nucleotide polymorphisms (mtSNPs) that account for the maximum level of discrimination power in a given population. Applications of the selected mtSNP panel range from anthropologic and medical studies to forensic genetic casework.

Methodology/Principal Findings

This study proposes a new simulation-based method that explores the ability of different mtSNP panels to yield the maximum levels of discrimination power. The method explores subsets of mtSNPs of different sizes randomly chosen from a preselected panel of mtSNPs based on frequency. More than 2,000 complete genomes representing three main continental human population groups (Africa, Europe, and Asia) and two admixed populations (“African-Americans” and “Hispanics”) were collected from GenBank and the literature, and were used as training sets. Haplotype diversity was measured for each combination of mtSNP and compared with existing mtSNP panels available in the literature. The data indicates that only a reduced number of mtSNPs ranging from six to 22 are needed to account for 95% of the maximum haplotype diversity of a given population sample. However, only a small proportion of the best mtSNPs are shared between populations, indicating that there is not a perfect set of “universal” mtSNPs suitable for all population contexts. The discrimination power provided by these mtSNPs is much higher than the power of the mtSNP panels proposed in the literature to date. Some mtSNP combinations also yield high diversity values in admixed populations.

Conclusions/Significance

The proposed computational approach for exploring combinations of mtSNPs that optimise the discrimination power of a given set of mtSNPs is more efficient than previous empirical approaches. In contrast to precedent findings, the results seem to indicate that only few mtSNPs are needed to reach high levels of discrimination power in a population, independently of its ancestral background.

Testing the performance of mtSNP minisequencing in forensic samples

There is a growing interest among forensic geneticists in developing efficient protocols for genotyping coding region mitochondrial DNA (mtDNA) SNPs (mtSNPs). Minisequencing is becoming a popular method for SNP genotyping, but it is still used by few forensic laboratories. In part, this is due to the lack of studies testing its efficiency and reproducibility when applied to real and complex forensic samples. Here we tested a minisequencing design that consists of 71 mtSNPs (in three multiplexes) that are diagnostic of known branches of the R0 phylogeny, in real forensic samples, including degraded bones and teeth, hair shafts, and serial dilutions. The fact that amplicons are short coupled with the natural efficiency of the minisequencing technique allow these assays to perform well with all the samples tested either degraded and/or those containing low DNA amount. We did not observe phylogenetic inconsistencies in the 71 mtSNP haplotypes generated, indicating that the technique is robust against potential artefacts that could arise from unintended contamination and/or spurious amplification of nuclear mtDNA pseudogenes (NUMTs).

Comparative analysis of the HV1 and HV2 regions of human mitochondrial DNA by denaturing high-performance liquid chromatography

Denaturing high-performance liquid chromatography (DHPLC) was evaluated as a sequencing-independent means of detecting the presence of sequence differences in pair-wise mixtures of nonconcordant amplicons of human mitochondrial DNA (mtDNA). A total of 920 pair-wise combinations of HV1 and HV2 mtDNA amplicons from 95 individuals were assayed by DHPLC for sequence concordance/nonconcordance. For the 72 combinations of amplicons from different individuals who shared identical DNA sequences, DHPLC assays consistently indicated sequence concordance between the samples. This was in 100% agreement with sequencing data. For the 849 combinations of amplicons which differed in sequence, DHPLC detected the presence of sequence nonconcordance in all but 13 assays to yield 98.5% concordance with sequencing. Thus, DHPLC can be used to detect a diversity of sequence differences (transitions, transversions, insertions, and deletions) in the mtDNA D-loop. Accordingly, DHPLC may have utility as a presumptive indicator of mtDNA sequence concordance samples, as a screen for heteroplasmy/situational mixtures, and as a means for the physical fractionation of the individual contributors to an mtDNA mixture prior to sequencing.

Rapid identification of unknown heteroplasmic mitochondrial DNA mutations with mismatch-specific surveyor nuclease

Identification of mitochondrial DNA (mtDNA) mutations is essential for diagnosis and genetic counseling of mitochondrial diseases. In this chapter, we describe a strategy for the rapid identification of heteroplasmic mtDNA mutations that can be used routinely in molecular genetic laboratories. This protocol involves the following three steps: (i) PCR amplification of the entire human mitochondrial genome with 17 overlapping PCR products, (ii) localization of mtDNA mismatch(es) after digestion of the 17 amplicons by Surveyor Nuclease, a member of a family of plant DNA endonucleases that cleave double-strand DNA at any mismatch site, and (iii) identification of the mutation by sequencing the region containing the mismatch.

Minisequencing mitochondrial DNA pathogenic mutations

BACKGROUND

There are a number of well-known mutations responsible of common mitochondrial DNA (mtDNA) diseases. In order to overcome technical problems related to the analysis of complete mtDNA genomes, a variety of different techniques have been proposed that allow the screening of coding region pathogenic mutations.

METHODS

We here propose a minisequencing assay for the analysis of mtDNA mutations. In a single reaction, we interrogate a total of 25 pathogenic mutations distributed all around the whole mtDNA genome in a sample of patients suspected for mtDNA disease.

RESULTS

We have detected 11 causal homoplasmic mutations in patients suspected for Leber disease, which were further confirmed by standard automatic sequencing. Mutations m.11778G>A and m.14484T>C occur at higher frequency than expected by change in the Galician (northwest Spain) patients carrying haplogroup J lineages (Fisher's Exact test, P-value < 0.01). The assay performs well in mixture experiments of wild:mutant DNAs that emulate heteroplasmic conditions in mtDNA diseases.

CONCLUSION

We here developed a minisequencing genotyping method for the screening of the most common pathogenic mtDNA mutations which is simple, fast, and low-cost. The technique is robust and reproducible and can easily be implemented in standard clinical laboratories.

Rapid identification of mitochondrial DNA (mtDNA) mutations in neuromuscular disorders by using surveyor strategy

Mutations of mitochondrial genome are responsible for respiratory chain defects in numerous patients. We have used a strategy, based on the use of a mismatch-specific DNA endonuclease named " Surveyor Nuclease", for screening the entire mtDNA in a group of 50 patients with neuromuscular features, suggesting a respiratory chain dysfunction. We identified mtDNA mutations in 20% of patients (10/50). Among the identified mutations, four are not found in any mitochondrial database and have not been reported previously. We also confirm that mtDNA polymorphisms are frequently found in a heteroplasmic state (15 different polymorphisms were identified among which five were novel).

Rapid identification of unknown heteroplasmic mutations across the entire human mitochondrial genome with mismatch-specific Surveyor Nuclease

Mitochondrial DNA (mtDNA) mutations are responsible for mitochondrial diseases in numerous patients. But, until now, no rapid method was available for the identification of unknown deleterious point mutations. Here, we describe a new strategy for the rapid identification of heteroplasmic mtDNA mutations using mismatch-specific Surveyor Nuclease. This protocol involves the following three steps: (i) PCR amplification of the entire human mitochondrial genome in 17 overlapping fragments; (ii) localization of mtDNA mismatch(es) after digestion of the 17 amplicons by Surveyor Nuclease; and (iii) identification of the mutation by sequencing the region containing the mismatch. This Surveyor Nuclease-based strategy allows a systematic screening of the entire mtDNA to identify a mutation within 2 days. It represents an important diagnostic approach for mitochondrial diseases that can be routinely used in molecular diagnostic laboratories.

Coding region mitochondrial DNA SNPs: targeting East Asian and Native American haplogroups

We have developed a single PCR multiplex SNaPshot reaction that consists of 32 coding region SNPs that allows (i) increasing the discrimination power of the mitochondrial DNA (mtDNA) typing in forensic casework, and (ii) haplogroup assignments of mtDNA profiles in both human population studies (e.g. anthropological) and medical research. The selected SNPs target the East Asian phylogeny, including its Native American derived branches. We have validated this multiplex assay by genotyping a sample of East Asians (Taiwanese) and Native Americans (Argentineans). In addition to the coding SNP typing, we have sequenced the complete control region for the same samples. The genotyping results (control region plus SNaPshot profiles) are in good agreement with previous human population genetic studies (based on e.g. complete sequencing) and the known mtDNA phylogeny. We observe that the SNaPshot method is reliable, rapid, and cost effective in comparison with other techniques of multiplex SNP genotyping. We discuss the advantages of our SNP genotyping selection with respect to previous attempts, and we highlight the importance of using the known mtDNA phylogeny as a framework for SNP profile interpretation and as a tool to minimize genotyping errors.

Progressive accumulation of mitochondrial DNA mutations and decline in mitochondrial function lead to beta-cell failure

A key adaptation enabling the fetus to survive in a limited energy environment may be the reprogramming of mitochondrial function, which can have deleterious effects. Critical questions are whether mitochondrial dysfunction progressively declines after birth, and if so, what mechanism might underlie this process. To address this, we developed a model of intrauterine growth retardation (IUGR) in the rat that leads to diabetes in adulthood. Reactive oxygen species (ROS) production and oxidative stress gradually increased in IUGR islets. ATP production was impaired and continued to deteriorate with age. The activities of complex I and III of the electron transport chain progressively declined in IUGR islets. Mitochondrial DNA point mutations accumulated with age and were associated with decreased mitochondrial DNA content and reduced expression of mitochondria-encoded genes in IUGR islets. Mitochondrial dysfunction resulted in impaired insulin secretion. These results demonstrate that IUGR induces mitochondrial dysfunction in the fetal beta-cell, leading to increased production of ROS, which in turn damage mitochondrial DNA. A self-reinforcing cycle of progressive deterioration in mitochondrial function leads to a corresponding decline in beta-cell function. Finally, a threshold in mitochondrial dysfunction and ROS production is reached, and diabetes ensues.

Separating human DNA mixtures using denaturing high-performance liquid chromatography

DNA mixtures represent challenging samples that are rarely amenable to direct DNA sequence analysis and many of the strategies available to separate mixtures are both labor and time intensive. Denaturing high-performance liquid chromatography is an accurate and rapid approach for the detection and scoring of mutations. It can also be used to separate DNA mixtures. The technique relies on the chromatographic separation of crosshybridization products to isolate the individual components of a mixture. By eliminating secondary amplification and excessive manipulation prior to sequencing, denaturing high-performance liquid chromatography can streamline the analysis of conditions ranging from somatic mosaicism, microchimerism and mitochondrial heteroplasmy to evidentiary material containing mixtures of DNA encountered in forensic investigations.

Surveyor™ Nuclease: A new strategy for a rapid identification of heteroplasmic mitochondrial DNA mutations in patients with respiratory chain defects

Molecular analysis of mitochondrial DNA (mtDNA) is a critical step in diagnosis and genetic counseling of respiratory chain defects. No fast method is currently available for the identification of unknown mtDNA point mutations. We have developed a new strategy based on complete mtDNA PCR amplification followed by digestion with a mismatch-specific DNA endonuclease, Surveyor Nuclease. This enzyme, a member of the CEL nuclease family of plant DNA endonucleases, cleaves double-strand DNA at any mismatch site including base substitutions and small insertions/deletions. After digestion, cleavage products are separated and analyzed by agarose gel electrophoresis. The size of the digestion products indicates the location of the mutation, which is then confirmed and characterized by sequencing. Although this method allows the analysis of 2 kb mtDNA amplicons and the detection of multiple mutations within the same fragment, it does not lead to the identification of homoplasmic base substitutions. Homoplasmic pathogenic mutations have been described. Nevertheless, most homoplasmic base substitutions are neutral polymorphisms while deleterious mutations are typically heteroplasmic. Here, we report that this method can be used to detect mtDNA mutations such as m.3243A>G tRNA(Leu) and m.14709T>C tRNA(Glu) even when they are present at levels as low as 3% in DNA samples derived from patients with respiratory chain defects. Then, we tested five patients suffering from a mitochondrial respiratory chain defect and we identified a variant (m.16189T>C) in two of them, which was previously associated with susceptibility to diabetes and cardiomyopathy. In conclusion, this method can be effectively used to rapidly and completely screen the entire human mitochondrial genome for heteroplasmic mutations and in this context represents an important advance for the diagnosis of mitochondrial diseases.

Techniques and pitfalls in the detection of pathogenic mitochondrial DNA mutations

Mutations in the mitochondrial DNA (mtDNA) are now recognized as major contributors to human pathologies and possibly to normal aging. A large number of rearrangements and point mutations in protein coding and tRNA genes have been identified in patients with mitochondrial disorders. In this review, we discuss genotype-phenotype correlations in mitochondrial diseases and common techniques used to identify pathogenic mtDNA mutations in human tissues. Although most of these approaches employ standard molecular biology tools, the co-existence of wild-type and mutated mtDNA (mtDNA heteroplasmy) in diseased tissues complicates both the detection and accurate determination of the size of the mutated fractions. To address these problems, novel approaches were developed and are discussed in this review.

Single-strand conformation polymorphism of hyper-variable regions HV1 and HV2 of human mitochondrial DNA: detection by silver staining

Human mitochondrial DNA from 50 trios consisting of mother (M), child (C) and father (F) was PCR amplified with primers flanking the hyper-variable regions, HVR1 and HVR2. The amplified products were then fractionated under non-denaturing conditions, silver-stained and compared by single-stranded conformational polymorphism (SSCP). In all but one case, mother and child displayed identical patterns, which could be promptly distinguished from that of the father. For the remaining cases, either set of primers was sufficient to resolve the familial ties. In no instance, M displayed alleles different from those of C within each trio, demonstrating that no false exclusions occurred. The SSCP approach proved to be a robust technique suitable as a preliminary screening in cases requiring identification of multiple samples.

Fluorescent SSCP of overlapping fragments (FSSCP-OF): a highly sensitive method for the screening of mitochondrial DNA variation

The mtDNA analysis (mtDNA) is increasingly being demanded for forensic purposes due to the fact that many times the use of standard nuclear marker fails to analyze degraded samples (such as bones) and specially for the analysis of hair shafts (a common sample in the crime scene). However, analysis of mtDNA sequencing implies a great lab effort when a high number of samples must be analyzed. The present work introduces a novel and reliable method for the screening of mtDNA variation in the first and second hypervariables (HV1 and HV2) regions which we have denominated fluorescent single strand conformation polymorphism (SSCP) of overlapping fragments (FSSCP-OF). FSSCP-OF is based on the basic theory of SSCP analysis and combines two complementary strategies: the use of PCR amplified overlapping fragments and fluorescent detection technology. The overlap region contains a high percentage (50%) of the d-loop mtDNA variation and for this reason, the probability to detect a polymorphic position by SSCP analysis is clearly increased in comparison to conventional SSCP methods due to the fact that the same polymorphic position is usually placed in a different "relative" position in the two overlapped fragments. The use of multicolor fluorescent technology allows also the multiplex amplification of overlapping fragment and its subsequent analysis in an automatic sequencer. We have analyzed 50 samples of unrelated individuals through the FSSCP-OF technique and we have found that using this methodology the probability to distinguish two samples with different sequences is close to 100%. FSSCP-OF has other important advantages with respect to previous screening methods, such as the automation and standardization of the protocols, which is of special interest for the forensic routine.

Mutation analysis of the entire mitochondrial genome using denaturing high performance liquid chromatography

In patients with mitochondrial disease a continuously increasing number of mitochondrial DNA (mtDNA) mutations and polymorphisms have been identified. Most pathogenic mtDNA mutations are heteroplasmic, resulting in heteroduplexes after PCR amplification of mtDNA. To detect these heteroduplexes, we used the technique of denaturing high performance liquid chromatography (DHPLC). The complete mitochondrial genome was amplified in 13 fragments of 1-2 kb, digested in fragments of 90-600 bp and resolved at their optimal melting temperature. The sensitivity of the DHPLC system was high with a lowest detection of 0.5% for the A8344G mutation. The muscle mtDNA from six patients with mitochondrial disease was screened and three mutations were identified. The first patient with a limb-girdle-type myopathy carried an A3302G substitution in the tRNA(Leu(UUR)) gene (70% heteroplasmy), the second patient with mitochondrial myopathy and cardiomyopathy carried a T3271C mutation in the tRNA(Leu(UUR)) gene (80% heteroplasmy) and the third patient with Leigh syndrome carried a T9176C mutation in the ATPase6 gene (93% heteroplasmy). We conclude that DHPLC analysis is a sensitive and specific method to detect heteroplasmic mtDNA mutations. The entire automatic procedure can be completed within 2 days and can also be applied to exclude mtDNA involvement, providing a basis for subsequent investigation of nuclear genes.

Reproducibility of mtDNA analysis between laboratories: a report of the European DNA Profiling Group (EDNAP)

The aim of this collaborative exercise was to determine whether uniformity of mtDNA sequencing results could be achieved among different EDNAP laboratories. Laboratories were asked to sequence mtDNAHV1 region (16024-16365) from three bloodstains, proceeding in accordance with the protocol and strategies currently used in each individual laboratory. Cycle sequencing was used by 11 laboratories and solid phase single stranded sequencing was used by one laboratory. Different PCR strategies and PCR conditions were used by the different laboratories. Three laboratories used semi-nested PCR, two nested PCR, three direct amplification of HV1 and four amplification of overlapping fragments covering the HV1 region. Despite the diversity of methodologies used, all the laboratories reported the same results. The successful result of this exercise shows that PCR based mtDNA typing by automated sequencing is a valid, robust and reliable means of forensic identification despite the different strategies and methodologies used by the different laboratories.

Rapid mitochondrial DNA typing using restriction enzyme digestion of polymerase chain reaction amplicons followed by capillary electrophoresis separation with laser-induced fluorescence detection

The polymorphic control region of mitochondrial DNA (mtDNA) is becoming more commonly used in forensic applications to differentiate among individuals in a population. Two hypervariable regions (HV1 and HV2) are often sequenced following amplification of the mtDNA via the polymerase chain reaction (PCR). More rapid screening assays would reduce both the effort and the expense of comparing two samples. A methodology has been developed that first uses restriction endonuclease digestion of the PCR-amplified mtDNA using RsaI and MnlI and then capillary electrophoresis (CE) to separate and size the PCR-RFLP fragments. This rapid procedure offers an alternative method for screening of polymorphisms in amplified mtDNA samples. In addition, the presence of a T-->C transition at position 16189, which gives rise to the so-called "C-stretch" in HV1, may be predicted from the presence of nonspecific PCR products in the CE results.