Fidelity of capture-enrichment for mtDNA genome sequencing: influence of NUMTs

Enriching target sequences in sequencing libraries via capture hybridization to bait/probes is an efficient means of leveraging the capabilities of next-generation sequencing for obtaining sequence data from target regions of interest. However, homologous sequences from non-target regions may also be enriched by such methods. Here we investigate the fidelity of capture enrichment for complete mitochondrial DNA (mtDNA) genome sequencing by analyzing sequence data for nuclear copies of mtDNA (NUMTs). Using capture-enriched sequencing data from a mitochondria-free cell line and the parental cell line, and from samples previously sequenced from long-range PCR products, we demonstrate that NUMT alleles are indeed present in capture-enriched sequence data, but at low enough levels to not influence calling the authentic mtDNA genome sequence. However, distinguishing NUMT alleles from true low-level mutations (e.g. heteroplasmy) is more challenging. We develop here a computational method to distinguish NUMT alleles from heteroplasmies, using sequence data from artificial mixtures to optimize the method.     

Mutations in mitochondrial DNA accumulate differentially in three different human tissues during ageing.

In 60 human tissue samples (encompassing skeletal muscle, heart and kidney) obtained from subjects aged from under 1 to 90 years, we used quantitative PCR procedures to quantify mitochondrial DNA (mtDNA) molecules carrying the 4977 bp deletion (mtDNA4977) and 3243 A-->G base substitution. In addition, the prevalence of multiple mtDNA deletions was assessed in a semi-quantitative manner. For all three tissues, the correlations between the accumulation of the particular mtDNA mutations and age of the subject are highly significant. However, differential extents of accumulation of the two specific mutations in the various tissues were observed. Thus, the mean abundance (percentage of mutant mtDNA out of total mtDNA) of mtDNA4977in a subset of age-matched adults is substantially higher in skeletal muscle than in heart and kidney. However, the mean abundance of the 3243 A-->G mutation in skeletal muscle was found to be lower than that in heart and kidney. Visualisation of arrays of PCR products arising from multiple mtDNA deletions in DNA extracted from adult skeletal muscle, was readily made after 30 cycles of PCR. By contrast, in DNA extracted from adult heart or kidney, amplification for 35 cycles of PCR was required to detect multiple mtDNA deletions. Although such multiple deletions are less abundant in heart and kidney than in skeletal muscle, in all tissue extracts there are unique patterns of bands, even from different tissues of the same subject. The differential accumulation of mtDNA4977, other mtDNA deletions and the 3243 A-->G mutation in the three tissues analysed presumably reflects different metabolic and senescence characteristics of these various tissues.     

Detection of mutations in mitochondrial DNA by droplet digital PCR

Mutations in mitochondrial DNA (mtDNA) may result in various pathological processes. Detection of mutant mtDNAs is a problem for diagnostic practice that is complicated by heteroplasmy – a phenomenon of the inferring presence of at least two allelic variants of the mitochondrial genome. Also, the level of heteroplasmy largely determines the profile and severity of clinical manifestations. Here we discuss detection of mutations in heteroplasmic mtDNA using up-todate methods that have not yet been introduced as routine clinical assays. These methods can be used for detecting mutations in mtDNA to verify diagnosis of “mitochondrial disease”, studying dynamics of mutant mtDNA in body tissues of patients, as well as investigating structural features of mtDNAs. Original data on allele-specific discrimination of m.11778G>A mutation by droplet digital PCR are presented, which demonstrate an opportunity for simultaneous detection and quantitative assessment of mutations in mtDNAs.     

Quantification of mitochondrial DNA mutation load

Mitochondrial DNA (mtDNA) mutations are an important cause of genetic disease and have been proposed to play a role in the ageing process. Quantification of total mtDNA mutation load in ageing tissues is difficult as mutational events are rare in a background of wild-type molecules, and detection of individual mutated molecules is beyond the sensitivity of most sequencing based techniques. The methods currently most commonly used to document the incidence of mtDNA point mutations in ageing include post-PCR cloning, single-molecule PCR and the random mutation capture assay. The mtDNA mutation load obtained by these different techniques varies by orders of magnitude, but direct comparison of the three techniques on the same ageing human tissue has not been performed. We assess the procedures and practicalities involved in each of these three assays and discuss the results obtained by investigation of mutation loads in colonic mucosal biopsies from ten human subjects.     

Molecular investigations of mitochondrial deletions: Evaluating the usefulness of different genetic tests

Deletions in mitochondrial DNA are a common cause of mitochondrial disorders. The molecular diagnosis of mtDNA deletions for years was based on Southern hybridization later replaced by PCR methods such as PCR with primers specific for a particular deletion (mainly the so-called common deletion of 4977bp) and long PCR. In order to evaluate the usefulness of MLPA (Multiplex Ligation-dependent Probe Amplification) in molecular diagnosis of large scale mtDNA deletions we compare four diagnostic methods: Southern hybridization, PCR, long-PCR and MLPA in a group of 16 patients with suspected deletions. Analysis was performed on blood, muscle and in one case hepatic tissue DNA. The MLPA was not able to confirm all the deletions detected by PCR methods, but due to its relative ease of processing, minimal equipment, low costs and the additional possibility to detect frequent point mtDNA mutations in one assay it is worth considering as a screening method. We recommend to always confirm MLPA results by PCR methods.     

No recombination of mtDNA after heteroplasmy for 50 generations in the mouse maternal germline

Variants of mitochondrial DNA (mtDNA) are commonly used as markers to track human evolution because of the high sequence divergence and exclusive maternal inheritance. It is assumed that the inheritance is clonal, i.e. that mtDNA is transmitted between generations without germline recombination. In contrast to this assumption, a number of studies have reported the presence of recombinant mtDNA molecules in cell lines and animal tissues, including humans. If germline recombination of mtDNA is frequent, it would strongly impact phylogenetic and population studies by altering estimates of coalescent time and branch lengths in phylogenetic trees. Unfortunately, this whole area is controversial and the experimental approaches have been widely criticized as they often depend on polymerase chain reaction (PCR) amplification of mtDNA and/or involve studies of transformed cell lines. In this study, we used an in vivo mouse model that has had germline heteroplasmy for a defined set of mtDNA mutations for more than 50 generations. To assess recombination, we adapted and validated a method based on cloning of single mtDNA molecules in the λ phage, without prior PCR amplification, followed by subsequent mutation analysis. We screened 2922 mtDNA molecules and found no germline recombination after transmission of mtDNA under genetically and evolutionary relevant conditions in mammals.     

Preferential amplification is minimised in long-PCR systems

The advent of long PCR (XL-PCR) has proven to be a major advance in PCR technology and is currently being utilised to investigate numerous biological systems. The analysis of mixed DNA populations is a particularly useful application for XL-PCR. For example, XL-PCR has been used to investigate the occurrence of heterogeneous mitochondrial DNA (mtDNA) rearrangement mutations. With XL-PCR it became possible to amplify the entire length of the mtDNA chromosome and detect any mtDNA deletion or insertion mutations based on a measurable change in overall sequence length. In the present communication, XL-PCR and conventional short-length PCR were used to amplify mitochondrial DNA sequences from several human vastus lateralis skeletal muscle samples. The experiments demonstrated that there was minimal preferential amplification of shorter DNA sequences with XL-PCR and was significantly less than the preferential amplification of shorter sequences observed with conventional PCR. Also, XL-PCR amplification of the complete mtDNA sequence from control DNA containing a single mtDNA template (leucocyte extracts) showed that the generation of PCR artefacts was not a predisposed failing of the system but was dependant on the standard rules that govern the set up and optimisation of any PCR reaction. In optimised systems, XL-PCR artefacts were not generated and a single PCR product was always recovered.     

A Standard Reference Material to determine the sensitivity of techniques for detecting low-frequency mutations, SNPs, and heteroplasmies in mitochondrial DNA

Human mitochondrial DNA (mtDNA) mutations are important for forensic identifications and mitochondrial disease diagnostics. Low-frequency mutations, heteroplasmies, or SNPs scattered throughout the DNA in the presence of a majority of mtDNA with the Cambridge Reference Sequence (CRS) are almost impossible to detect. Therefore, the National Institute of Science and Technology has developed heteroplasmic human mtDNA Standard Reference Material (SRM) 2394 to allow scientists to determine their sensitivity in detecting such differences. SRM 2394 is composed of mixtures ranging from 1/99 to 50/50 of two 285-bp PCR products from two cell lines that differ at one nucleotide position. Twelve laboratories using various mutation detection methods participated in a blind interlaboratory evaluation of a prototype of SRM 2394. Most of these procedures were unable to detect the mutation when present below 20%, an indication that, in many real-life cases, low-frequency mutations remain undetected and that more sensitive mutation detection techniques are urgently needed.     

Mitochondrial DNA deletion mutations in adult mouse cardiac side population cells

We investigated the presence and potential role of mitochondrial DNA (mtDNA) deletion mutations in adult cardiac stem cells. Cardiac side population (SP) cells were isolated from 12-week-old mice. Standard polymerase chain reaction (PCR) was used to screen for the presence of mtDNA deletion mutations in (a) freshly isolated SP cells and (b) SP cells cultured to passage 10. When present, the abundance of mtDNA deletion mutation was analyzed in single cell colonies. The effect of different levels of deletion mutations on SP cell growth and differentiation was determined. MtDNA deletion mutations were found in both freshly isolated and cultured cells from 12-week-old mice. While there was no significant difference in the number of single cell colonies with mtDNA deletion mutations from any of the groups mentioned above, the abundance of mtDNA deletion mutations was significantly higher in the cultured cells, as determined by quantitative PCR. Within a single clonal cell population, the detectable mtDNA deletion mutations were the same in all cells and unique when compared to deletions of other colonies. We also found that cells harboring high levels of mtDNA deletion mutations (i.e. where deleted mtDNA comprised more than 60% of total mtDNA) had slower proliferation rates and decreased differentiation capacities. Screening cultured adult stem cells for mtDNA deletion mutations as a routine assessment will benefit the biomedical application of adult stem cells.     

Rapid screening of the entire mitochondrial DNA for low-level heteroplasmic mutations

Alterations of the mitochondrial DNA (mtDNA) are implicated in various pathological conditions. In this study, we used denaturing high performance liquid chromatography (DHPLC) as a method to rapidly screen the entire mtDNA for mutations. Overlapping DNA fragments, amplified by one single cycling protocol from frozen pre-formulated PCR mixes, were subjected to DHPLC analysis. Single DHPLC injections of fragments yielded straightforward interpretation of results with a detection limit down to 1% mtDNA heteroplasmy. Furthermore, collection and re-amplification of low degree heteroduplex peak-fractions allowed sequence analysis of mtDNA mutations down to the detection limit of the DHPLC method. In order to demonstrate that the method has diagnostic value, we analyzed and confirmed known mtDNA mutations in patient samples.     

DNA extraction procedures meaningfully influence qPCR-based mtDNA copy number determination

Quantitative real time PCR (qPCR) is commonly used to determine cell mitochondrial DNA (mtDNA) copy number. This technique involves obtaining the ratio of an unknown variable (number of copies of an mtDNA gene) to a known parameter (number of copies of a nuclear DNA gene) within a genomic DNA sample. We considered the possibility that mtDNA:nuclear DNA (nDNA) ratio determinations could vary depending on the method of genomic DNA extraction used, and that these differences could substantively impact mtDNA copy number determination via qPCR. To test this we measured mtDNA:nDNA ratios in genomic DNA samples prepared using organic solvent (phenol–chloroform–isoamyl alcohol) extraction and two different silica-based column methods, and found mtDNA:nDNA ratio estimates were not uniform. We further evaluated whether different genomic DNA preparation methods could influence outcomes of experiments that use mtDNA:nDNA ratios as endpoints, and found the method of genomic DNA extraction can indeed alter experimental outcomes. We conclude genomic DNA sample preparation can meaningfully influence mtDNA copy number determination by qPCR.     

Defeating numts: semi-pure mitochondrial DNA from eggs and simple purification methods for field-collected wildlife tissues.

Mitochondrial DNA (mtDNA) continues to play a pivotal role in phylogeographic, phylogenetic, and population genetic studies. PCR amplification with mitochondrial primers often yields ambiguous sequences, in part because of the co-amplification of nuclear copies of mitochondrial genes (numts) and true mitochondrial heteroplasmy arising from mutations, hybridization with paternal leakage, gene duplications, and recombination. Failing to detect numts or to distinguish the origin of such homologous sequences results in the incorrect interpretation of data. However, few studies obtain purified mtDNA to confirm the mitochondrial origin of the first reference sequences for a species. Here, we demonstrate the importance and ease of obtaining semi-pure mtDNA from wildlife tissues, preserved under various typical field conditions, and investigate the success of 3 commercial extraction kits, cesium-chloride gradient mtDNA purification, long-template PCR amplification, cloning, and more species-specific degenerate primers. Using more detailed avian examples, we illustrate that unfertilized or undeveloped eggs provide the purest sources of mtDNA; that kits provide an alternative to cesium-chloride gradient methods; and that long-template PCR, cloning, and degenerate primers cannot be used to produce reliable mitochondrial reference sequences, but can be powerful tools when used in conjunction with purified mtDNA stocks to distinguish numts from true heteroplasmy.     

ND4 mutation in transitional cell carcinoma: Does mitochondrial mutation occur before tumorigenesis?

To investigate how mitochondrial mutation occurs in cancers, we analyzed ND4 mutation in 53 transitional cell carcinomas (TCCs) of the upper urinary tract and the normal counterpart (perirenal soft tissue). Three methods, i.e., DNA sequencing, restriction fragment length polymorphism (RFLP), and denaturing high-performance liquid chromatography (DHPLC), were employed because of their different sensitive of detecting mutation. The results of sequencing and RFLP showed that ND4 mutations were only found in 24.5% (13/53) of tumor. However, 11 of these mutations could also be identified in the normal tissue by DHPLC, indicating that most mitochondrial mutations identified in tumors preexist as minor components, which are too low in quantity to be detected by less sensitive methods such as DNA sequencing. The result suggests that mtDNA mutation occurs before tumorigenesis and become apparent in cancer cells.     

Origin of intra-individual variation in PCR-amplified mitochondrial cytochrome oxidase I of Thrips tabaci (Thysanoptera: Thripidae): mitochondrial heteroplasmy or nuclear integration?

The mitochondrial genome is increasingly being used as a species diagnostic marker in insects. Typically, genomic DNA is PCR amplified and then analysed by restriction analyses or sequencing. This analysis system may cause some serious problems for molecular diagnosis. Besides the errors introduced by the PCR process, mtDNA sequence variation of amplified fragments may originate from mtDNA heteroplasmy or from nuclear integrations of mtDNA fragments, both of which have been shown to occur in insects. Here we document abundant variation in PCR-amplified sequences of the mitochondrial cytochrome oxidase I gene of Thrips tabaci. We confirm that the most common haplotype is of mitochondrial origin. Some of the observed mutations were introduced by the amplification process. However, the occurrence of some haplotypes at elevated frequencies indicates that within-individual variation of the respective fragment exists at low levels in T. tabaci. The frequencies of these sequences are too low to negatively affect mtDNA-based molecular diagnosis of T. tabaci. The possible origin of these variant haplotypes is discussed.     

Accumulation of homoplasmic mtDNA point mutations in erythroblasts isolated from the bone marrow of patients with refractory anemia with ring sideroblasts (RARS)

It was hypothesised that mitochondrial iron overload in patients with refractory anemia with ring sideroblasts (RARS) results from mitochondrial DNA (mtDNA) mutations. To analyse the mtDNA sequence of iron storing mitochondria sensitively, we developed new protocols for selective erythroblasts isolation, mtDNA PCR amplification and sequencing. Using this approach, we found in each of the three RARS patients examined a unique spectrum of homoplasmic mtDNA point mutations affecting several mtDNA genes. Prediction analyses suggest that identified mutations do not result in major perturbations of mitochondrial functions and are tolerated. We discuss a mechanism explaining how the mutations identified may contribute to RARS pathogenesis.     

DNA supercoiling suppresses real-time PCR: a new approach to the quantification of mitochondrial DNA damage and repair

As a gold standard for quantification of starting amounts of nucleic acids, real-time PCR is increasingly used in quantitative analysis of mtDNA copy number in medical research. Using supercoiled plasmid DNA and mtDNA modified both in vitro and in cancer cells, we demonstrated that conformational changes in supercoiled DNA have profound influence on real-time PCR quantification. We showed that real-time PCR signal is a positive function of the relaxed forms (open circular and/or linear) rather than the supercoiled form of DNA, and that the conformation transitions mediated by DNA strand breaks are the main basis for sensitive detection of the relaxed DNA. This new finding was then used for sensitive detection of structure-mediated mtDNA damage and repair in stressed cancer cells, and for accurate quantification of total mtDNA copy number when all supercoiled DNA is converted into the relaxed forms using a prior heat-denaturation step. The new approach revealed a dynamic mtDNA response to oxidative stress in prostate cancer cells, which involves not only early structural damage and repair but also sustained copy number reduction induced by hydrogen peroxide. Finally, the supercoiling effect should raise caution in any DNA quantification using real-time PCR.     

Detection of mosaic pattern of mitochondrial DNA alterations in different populations of cells from the same endometrial tumor

Somatic mitochondrial DNA (mtDNA) alterations including point mutations and microsatellite instability (MSI) have been frequently detected in human cancers. To further explore the extensiveness of mtDNA alterations, we have analyzed the occurrence of somatic mtDNA mutations in different populations of endometrial cancer cells from the same tumor tissues as compared with adjacent non-tumor cells. Laser-captured micro-dissection was used to harvest endometrial cancer cells from separated areas of the same tumor and adjacent normal cells. Total DNA isolated from micro-dissected cells was PCR amplified and analyzed for mtDNA alterations by polyacrylamide gel electrophoresis and DNA sequencing. Multiple mtDNA alterations were detected in different portions of the same tumor. Different populations of endometrial cancer cells carried different patterns of mtDNA mutations. Interestingly, unlike previous reports, most mutations were found to be heteroplasmic. We have demonstrated the occurrence of hyper-variability of mtDNA alterations in a single piece of tumor tissue. Our observations support the hypothesis that the accumulation of mtDNA alterations is random and expands independently. The data presented here showed the heterogeneity of cancer cells in terms of mtDNA alterations in endometrial cancer.     

Detection and quantification of mitochondrial DNA deletions in individual cells by real-time PCR

Defects of mitochondrial DNA (mtDNA) are an important cause of disease and play a role in the ageing process. There are multiple copies of the mitochondrial genome in a single cell. In many patients with acquired or inherited mtDNA mutations, there exists a mixture of mutated and wild type genomes (termed heteroplasmy) within individual cells. As a biochemical and clinical defect is only observed when there are high levels of mutated mtDNA, a crucial investigation is to determine the level of heteroplasmic mutations within tissues and individual cells. We have developed an assay to determine the relative amount of deleted mtDNA using real-time fluorescence PCR. This assay detects the vast majority of deleted molecules, thus eliminating the need to develop specific probes. We have demonstrated an excellent correlation with other techniques (Southern blotting and three- primer competitive PCR), and have shown this technique to be sensitive to quantify the level of deleted mtDNA molecules in individual cells. Finally, we have used this assay to investigate patients with mitochondrial disease and shown in individual skeletal muscle fibres that there exist different patterns of abnormalities between patients with single or multiple mtDNA deletions. We believe that this technique has significant advantages over other methods to quantify deleted mtDNA and, employed alongside our method to sequence the mitochondrial genome from single cells, will further our understanding of the role of mtDNA mutations in human disease and ageing.     

The mutation rate of the human mtDNA deletion mtDNA4977.

The human mitochondrial mutation mtDNA4977 is a 4,977-bp deletion that originates between two 13-bp direct repeats. We grew 220 colonies of cells, each from a single human cell. For each colony, we counted the number of cells and amplified the DNA by PCR to test for the presence of a deletion. To estimate the mutation fate, we used a model that describes the relationship between the mutation rate and the probability that a colony of a given size will contain no mutants, taking into account such factors as possible mitochondrial turnover and mistyping due to PCR error. We estimate that the mutation rate for mtDNA4977 in cultured human cells is 5.95 x 10(-8) per mitochondrial genome replication. This method can be applied to specific chromosomal, as well as mitochondrial, mutations.     

Reduced Adolescent-Age Spatial Learning Ability Associated with Elevated Juvenile-Age Superoxide Levels in Complex I Mouse Mutants

Large-scale, heteroplasmic and generally pathogenic mtDNA defects (as induced by defective mitochondrial DNA polymerase, clonal mutations or DNA deletions) are known to negatively impact on life span and can result in apoptosis and tissue loss in, e.g., skeletal muscle or reduce learning abilities. The functional impact of homoplasmic specific mtDNA point mutations, e.g., in genes coding for the electron transport chain, however, remains a matter of debate. The present study contributes to this discussion and provides evidence that a single point mutation in complex I of the respiratory chain is associated with impairment of spatial navigation in adolescent (6-month-old) mice, i.e., reduced performance in the Morris Water Maze, which goes along with increased production of reactive oxygen species (ROS) in juvenile mice (3 months) but not at the age of phenotype expression. A point mutation in complex III goes along with only a mild and non-significant negative effect on cognitive performance and no significant changes in ROS production. These findings suggest to also consider the ontogenetic development of phenotypes when studying mtDNA mutations and highlights a possible impact of complex I dysfunction on the emergence of neurological deficits.