PCR-Free Enrichment of Mitochondrial DNA from Human Blood and Cell Lines for High Quality Next-Generation DNA Sequencing

Recent advances in sequencing technology allow for accurate detection of mitochondrial sequence variants, even those in low abundance at heteroplasmic sites. Considerable sequencing cost savings can be achieved by enriching samples for mitochondrial (relative to nuclear) DNA. Reduction in nuclear DNA (nDNA) content can also help to avoid false positive variants resulting from nuclear mitochondrial sequences (numts). We isolate intact mitochondrial organelles from both human cell lines and blood components using two separate methods: a magnetic bead binding protocol and differential centrifugation. DNA is extracted and further enriched for mitochondrial DNA (mtDNA) by an enzyme digest. Only 1 ng of the purified DNA is necessary for library preparation and next generation sequence (NGS) analysis. Enrichment methods are assessed and compared using mtDNA (versus nDNA) content as a metric, measured by using real-time quantitative PCR and NGS read analysis. Among the various strategies examined, the optimal is differential centrifugation isolation followed by exonuclease digest. This strategy yields >35% mtDNA reads in blood and cell lines, which corresponds to hundreds-fold enrichment over baseline. The strategy also avoids false variant calls that, as we show, can be induced by the long-range PCR approaches that are the current standard in enrichment procedures. This optimization procedure allows mtDNA enrichment for efficient and accurate massively parallel sequencing, enabling NGS from samples with small amounts of starting material. This will decrease costs by increasing the number of samples that may be multiplexed, ultimately facilitating efforts to better understand mitochondria-related diseases.     

Targeted multiplex next‐generation sequencing: advances in techniques of mitochondrial and nuclear DNA sequencing for population genomics

Next‐generation sequencing (NGS) is emerging as an efficient and cost‐effective tool in population genomic analyses of nonmodel organisms, allowing simultaneous resequencing of many regions of multi‐genomic DNA from multiplexed samples. Here, we detail our synthesis of protocols for targeted resequencing of mitochondrial and nuclear loci by generating indexed genomic libraries for multiplexing up to 100 individuals in a single sequencing pool, and then enriching the pooled library using custom DNA capture arrays. Our use of DNA sequence from one species to capture and enrich the sequencing libraries of another species (i.e. cross‐species DNA capture) indicates that efficient enrichment occurs when sequences are up to about 12% divergent, allowing us to take advantage of genomic information in one species to sequence orthologous regions in related species. In addition to a complete mitochondrial genome on each array, we have included between 43 and 118 nuclear loci for low‐coverage sequencing of between 18 kb and 87 kb of DNA sequence per individual for single nucleotide polymorphisms discovery from 50 to 100 individuals in a single sequencing lane. Using this method, we have generated a total of over 500 whole mitochondrial genomes from seven cetacean species and green sea turtles. The greater variation detected in mitogenomes relative to short mtDNA sequences is helping to resolve genetic structure ranging from geographic to species‐level differences. These NGS and analysis techniques have allowed for simultaneous population genomic studies of mtDNA and nDNA with greater genomic coverage and phylogeographic resolution than has previously been possible in marine mammals and turtles.     

Isolation and introgression in the Intermountain West: contrasting gene genealogies reveal the complex biogeographic history of the American pika (Ochotona princeps)

Aim We studied the history of colonization, diversification and introgression among major phylogroups in the American pika, Ochotona princeps (Lagomorpha), using comparative and statistical phylogeographic methods. Our goal was to understand how Pleistocene climatic fluctuations have shaped the distribution of diversity at mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) loci in this alpine specialist. Location North America’s Intermountain West. Methods We accumulated mtDNA sequence data (c. 560–1700 bp) from 232 pikas representing 64 localities, and sequenced two nuclear introns (mast cell growth factor, c. 550 bp, n = 148; protein kinase C iota, c. 660 bp, n = 139) from a subset of individuals. To determine the distribution of major mtDNA lineages, we conducted a phylogenetic analysis on the mtDNA sequence data, and we calculated divergence times among the lineages using a Bayesian Markov chain Monte Carlo approach. Relationships among nuclear alleles were explored with minimum spanning networks. Finally, we conducted coalescent simulations of alternative models of population history to test for congruence between nDNA and mtDNA responses to Pleistocene glacial cycles. Results We found that: (1) all individuals could be assigned to one of five allopatric mtDNA lineages; (2) lineages are associated with separate mountain provinces; (3) lineages originated from at least two rounds of differentiation; (4) nDNA and mtDNA markers exhibited overall phylogeographic congruence; and (5) introgression among phylogroups has occurred at nuclear loci since their initial isolation. Main conclusions Pika populations associated with different mountain systems have followed separate but not completely independent evolutionary trajectories through multiple glacial cycles. Range expansion associated with climate cooling (i.e. glaciations) promoted genetic admixture among populations within mountain ranges. It also permitted periodic contact and introgression between phylogroups associated with different mountain systems, the record of which is retained at nDNA but not mtDNA loci. Evidence for different histories at nuclear and mtDNA loci (i.e. periodic introgression versus deep isolation, respectively) emphasizes the importance of multilocus perspectives for reconstructing complete population histories.     

A decreased mitochondrial DNA content is related to insulin resistance in adolescents

The aim of this study was to investigate whether mitochondrial DNA (mtDNA) content is associated with insulin resistance (IR) in a sample of adolescents with features of metabolic syndrome. We further studied the link between polymorphisms in three genes involved in mitochondrial biogenesis and the presence of deleted mtDNA and mtDNA content. Data and blood samples were collected from 175 adolescents out of a cross-sectional, population-based study of 934 high school students. On the basis of the median value of homeostasis model assessment of IR (HOMA-IR) of the whole sample (2.2), the population was divided into two groups: noninsulin resistance (NIR) and IR. mtDNA quantification using nuclear DNA (nDNA) as a reference was carried out using a real-time quantitative PCR method. Genotyping for peroxisome proliferator-activated receptor-gamma (PPAR-gamma) (pro12Ala), PPAR- gamma coactivator-1alpha (PGC-1alpha) (Gly482Ser), and Tfam (rs1937 and rs12247015) polymorphisms was performed by PCR-based restriction fragment length polymorphism. Long-extension PCR was performed to amplify the whole mitochondrial genome. The mtDNA/nDNA ratio was significantly lower in the IR group (median: 9.08, range: 68.94) in comparison with the NIR group (12.24, 71.92) (P<0.03). Besides, the mtDNA/nDNA ratio was inversely correlated with HOMA (R: -0.18, P<0.02), glucose (R: -0.21, P<0.008), and uric acid (R: -0.18, P<0.03). Genotypes for the PPAR- gamma, PGC-1alpha, and Tfam variants were not associated with the mtDNA/nDNA ratio. Long-extension PCR did not show significant levels of mtDNA deletions. In conclusion, our findings indicate that reduced mtDNA content in peripheral leukocytes is associated with IR. This result seems not to be related with the previously mentioned variants in genes involved in the regulation of mitochondrial biogenesis.     

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.     

The effect of photodynamic treatment with thiopyronine on yeast nDNA and mtDNA

Summary In order to ascertain the main target of the photodynamic effect with the sensitizer thiopyronine (TP) and its interference with cellular DNA, the uptake of TP into the yeastSaccharomyces cerevisiae and the distribution of the dye within the cells were studied. After fractionation of the cellular components, about 1% of the TP was found to be bound to the nucleus in anaerobically grown yeast cells; in aerobically grown cells, about 7% could be detected at the mitochondria.After careful isolation of the DNA from the organelles, only 0.024% of the dye once taken up by the cells was detectable at nuclear DNA, whereas 0.652% was bound to mitochondrial DNA. A calculation of the number of TP-molecules bound per nucleotides revealed a ration of one molecule TP per about 6,000 nucleotides of nuclear DNA and one molecule TP per 53 nucleotides of mitochondrial DNA. In vitro, the maximal binding capacity was estimated to be about one molecule TP per nucleotide.The induction of single strand breaks in the DNA after photodynamic treatment in vivo and in vitro was investigated by comparing the sedimentation of nDNA and mtDNA through alkaline sucrose-gradients. No differences in the sedimentation profiles of nDNA after photodynamic treatment of the cells in vivo as compared to the untreated control could be observed. In contrast, the sedimentation coefficient of mtDNA was significally decreased after photodynamic treatment, indicating the induction of single strand breaks in vivo only in mtDNA and not in nDNA.     

PERCHED AT THE MITO‐NUCLEAR CROSSROADS: DIVERGENT MITOCHONDRIAL LINEAGES CORRELATE WITH ENVIRONMENT IN THE FACE OF ONGOING NUCLEAR GENE FLOW IN AN AUSTRALIAN BIRD

Relationships among multilocus genetic variation, geography, and environment can reveal how evolutionary processes affect genomes. We examined the evolution of an Australian bird, the eastern yellow robin Eopsaltria australis, using mitochondrial (mtDNA) and nuclear (nDNA) genetic markers, and bioclimatic variables. In southeastern Australia, two divergent mtDNA lineages occur east and west of the Great Dividing Range, perpendicular to latitudinal nDNA structure. We evaluated alternative scenarios to explain this striking discordance in landscape genetic patterning. Stochastic mtDNA lineage sorting can be rejected because the mtDNA lineages are essentially distinct geographically for > 1500 km. Vicariance is unlikely: the Great Dividing Range is neither a current barrier nor was it at the Last Glacial Maximum according to species distribution modeling; nuclear gene flow inferred from coalescent analysis affirms this. Female philopatry contradicts known female‐biased dispersal. Contrasting mtDNA and nDNA demographies indicate their evolutionary histories are decoupled. Distance‐based redundancy analysis, in which environmental temperatures explain mtDNA variance above that explained by geographic position and isolation‐by‐distance, favors a nonneutral explanation for mitochondrial phylogeographic patterning. Thus, observed mito‐nuclear discordance accords with environmental selection on a female‐linked trait, such as mtDNA, mtDNA–nDNA interactions or genes on W‐chromosome, driving mitochondrial divergence in the presence of nuclear gene flow.     

Nuclear and mitochondrial DNA repair: similar pathways?

Mitochondrial DNA (mtDNA) alterations are implicated in a broad range of human diseases and alterations of the mitochondrial genome are assumed to be a result of its high susceptibility to oxidative damage and its limited DNA repair compared to nuclear DNA (nDNA). Characterization of DNA repair mechanisms has generally focused on these processes in nDNA but increasing interest and research effort have contributed to our knowledge of the mechanisms underlying DNA repair in mitochondria. In this review, we make comparisons between nDNA and mtDNA repair pathways and propose a model for how these pathways interact in mitochondria.     

Introgression and selection shaped the evolutionary history of sympatric sister‐species of coral reef fishes (genus: Haemulon)

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister‐species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site‐associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister‐species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174–0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.     

Evaluating DNA degradation rates in faecal pellets of the endangered pygmy rabbit

Noninvasive genetic sampling of faecal pellets can be a valuable method for monitoring rare and cryptic wildlife populations, like the pygmy rabbit (Brachylagus idahoensis). To investigate this method's efficiency for pygmy rabbit monitoring, we evaluated the effect of sample age on DNA degradation in faecal pellets under summer field conditions. We placed 275 samples from known individuals in natural field conditions for 1–60 days and assessed DNA quality by amplifying a 294‐base‐pair (bp) mitochondrial DNA (mtDNA) locus and five nuclear DNA (nDNA) microsatellite loci (111–221 bp). DNA degradation was influenced by sample age, DNA type, locus length and rabbit sex. Both mtDNA and nDNA exhibited high PCR success rates (94.4%) in samples <1 day old. Success rates for microsatellite loci declined rapidly from 80.0% to 42.7% between days 5 and 7, likely due to increased environmental temperature. Success rates for mtDNA amplification remained higher than nDNA over time, with moderate success (66.7%) at 21 days. Allelic dropout rates were relatively high (17.6% at <1 day) and increased to 100% at 60 days. False allele rates ranged from 0 to 30.0% and increased gradually over time. We recommend collecting samples as fresh as possible for individual identification during summer field conditions. Our study suggests that this method can be useful for future monitoring efforts, including occupancy surveys, individual identification, population estimation, parentage analysis and monitoring of genetic diversity both of a re‐introduced population in central Washington and across their range.     

Hidden support from unpromising data sets strongly unites snakes with anguimorph ‘lizards’

A combined analysis of nuclear, mitochondrial and morphological data robustly resolves snakes as the sister taxon to anguimorph ‘lizards’. Analysed in isolation, nuclear DNA (nDNA) produces a trichotomy between snakes, iguanians and anguimorphs, mitochondrial DNA (mtDNA) is largely uninformative at deeper levels, and morphology tends to nest snakes deep within anguimorphs or with various legless squamate groups. When analysed simultaneously, the nuclear signal is sufficiently strong that mtDNA and morphology are constrained to choose between alternative resolutions of the iguanian–anguimorph–snake trichotomy (generated by the nDNA) – and both support the snake–anguimorph solution. Combined analyses of fast‐evolving or idiosyncratically evolving markers (mtNDA, morphology) with conservative traits (e.g. nuclear genes) might be the best way to resolve ancient, closely spaced divergences. Fast or idiosyncratic markers potentially provide the most information about short, ancient internodes, but can converge on spurious trees if analysed in isolation. However, if constrained to only choosing between plausible trees, such data can contribute unique and valuable phylogenetic signal that resolves such problematic divergences.     

Diversity and demography in Beringia: multilocus tests of paleodistribution models reveal the complex history of arctic ground squirrels

To assess effects of historical climate change on northern species, we quantified the population history of the arctic ground squirrel (Spermophilus parryii), an arctic-adapted rodent that evolved in Beringia and was strongly influenced by climatic oscillations of the Quaternary. Competing hypotheses for the species' population history were derived from patterns of mitochondrial (mtDNA) structure and a bioclimatic envelope model (BEM). Hypotheses invoked (1) sequential isolation of regional populations beginning with the Arctic, (2) deep isolation only across central Alaska, and (3) widespread panmixia, and were tested using coalescent methods applied to eight nuclear (nDNA) loci. The data rejected strict interpretations of all three hypotheses, but perspectives underlying each encompassed aspects of the species' history. Concordance between mtDNA and nDNA geographic structure revealed three semi-independently evolving phylogroups, whereas signatures of gene flow at nDNA loci were consistent with a historical contact between certain populations as inferred by the BEM. Demographic growth was inferred for all regions despite expectations of postglacial habitat contraction for parts of Beringia. Our results highlight the complementary perspectives on species' histories that multiple lines of evidence provide, and underscore the utility of multilocus data for resolving complex population histories relevant to understanding effects of climate change.     

A duplex real-time PCR assay for detection of drug-induced mitochondrial DNA depletion in HepG2 cells

Screening of drug-induced mitochondrial DNA (mtDNA) depletion during early preclinical drug development is of major interest. Here we describe the establishment of a novel duplex calibrator-normalized real-time polymerase chain reaction (PCR) assay for rapid and reliable quantification of mtDNA in HepG2 cells. This assay involves quantification of an mtDNA target gene (cytochrome b) relative to a nuclear DNA (nDNA) reference gene (β-actin) in one tube. The assay was evaluated for its precision, linearity, and reproducibility, and reliable detection of mtDNA depletion was demonstrated. Using this novel real-time PCR assay, drug-induced mtDNA depletion could be accurately detected.     

An efficient method for screening faecal DNA genotypes and detecting new individuals and hybrids in the red wolf (<Emphasis Type="Italic">Canis rufus</Emphasis>) experimental population area

Previously, sequencing of mitochondrial DNA (mtDNA) from non-invasively collected faecal material (scat) has been used to help manage hybridization in the wild red wolf (Canis rufus) population. This method is limited by the maternal inheritance of mtDNA and the inability to obtain individual identification. Here, we optimize the use of nuclear DNA microsatellite markers on red wolf scat DNA to distinguish between individuals and detect hybrids. We develop a data filtering method in which scat genotypes are compared to known blood genotypes to reduce the number of PCR amplifications needed. We apply our data filtering method and the more conservative maximum likelihood ratio method (MLR) of Miller et al. (2002 Genetics 160:357–366) to a scat dataset previously screened for hybrids by sequencing of mtDNA. Using seven microsatellite loci, we obtained genotypes for 105 scats, which were matched to 17 individuals. The PCR amplification success rate was 50% and genotyping error rates ranged from 6.6% to 52.1% per locus. Our data filtering method produced comparable results to the MLR method, and decreased the time and cost of analysis by 25%. Analysis of this dataset using our data filtering method verified that no hybrid individuals were present in the Alligator River National Wildlife Refuge, North Carolina in 2000. Our results demonstrate that nuclear DNA microsatellite analysis of red wolf scats provides an efficient and accurate approach to screen for new individuals and hybrids.     

Simultaneous quantitative assessment of circulating cell-free mitochondrial and nuclear DNA by multiplex real-time PCR

Quantification of circulating nucleic acids in plasma and serum could be used as a non-invasive diagnostic tool for monitoring a wide variety of diseases and conditions. We describe here a rapid, simple and accurate multiplex real-time PCR method for direct synchronized analysis of circulating cell-free (ccf) mitochondrial (mtDNA) and nuclear (nDNA) DNA in plasma and serum samples. The method is based on one-step multiplex real-time PCR using a FAM-labeled MGB probe and primers to amplify the mtDNA sequence of the ATP 8 gene, and a VIC-labeled MGB probe and primers to amplify the nDNA sequence of the glycerinaldehyde-3-phosphate-dehydrogenase (GAPDH) gene, in plasma and serum samples simultaneously. The efficiencies of the multiplex assays were measured in serial dilutions. Based on the simulation of the PCR reaction kinetics, the relative quantities of ccf mtDNA were calculated using a very simple equation. Using our optimised real-time PCR conditions, close to 100% efficiency was obtained from the two assays. The two assays performed in the dilution series showed very good and reproducible correlation to each other. This optimised multiplex real-time PCR protocol can be widely used for synchronized quantification of mtDNA and nDNA in different samples, with a very high rate of efficiency.     

Disturbed population genetics: suspected introgressive hybridization between two Mnais damselfly species (Odonata)

Mnais costalis and M. pruinosa are damselflies (Odonata: Calopterygidae) with low dispersal abilities, both during their aquatic stream-living immature stage and their flying adult stage. A previous nuclear DNA (nDNA) sequencing and morphology study showed that these two species are very closely related, and cohabit widely in western Japan. The two species, however, segregate microhabitats along a stream: M. costalis lives in the lower reaches, and M. pruinosa in the upper reaches. In this study, our analyses were based on mitochondrial DNA (mtDNA), which usually mutates faster and is more variable among individuals than nDNA, and which is inherited maternally. We found that most COI haplotypes were shared between the two species, and that for most study sites interspecific riverine genetic structures were not clarified by mtDNA analysis. Incongruent population genetic structures based on nDNA and mtDNA suggested hybridization and introgression of mtDNA between the two species.     

The long amplicon quantitative PCR for DNA damage assay as a sensitive method of assessing DNA damage in the environmental model, Atlantic killifish (Fundulus heteroclitus)

DNA damage is an important mechanism of toxicity for a variety of pollutants, and therefore, is often used as an indicator of pollutant effects in ecotoxicological studies. Here, we adapted a PCR-based assay for nuclear and mitochondrial DNA damage for use in an important environmental model, the Atlantic killifish (Fundulus heteroclitus). We refer to this assay as the long amplicon quantitative PCR (LA-QPCR) assay. To validate this method in killifish, DNA damage was measured in liver, brain, and muscle of fish dosed with 10 mg/kg benzo[a]pyrene. This exposure caused 0.4-0.8 lesions/10 kb. We also measured DNA damage in liver and muscle tissues from killifish inhabiting a Superfund site, confirming the utility of this method for biomonitoring. In both cases, damage levels were comparable in nuclear DNA (nDNA) and mitochondrial DNA (mtDNA). Since extensive nDNA sequence data are not readily available for many environmentally relevant species, but mitochondrial genomes are frequently fully sequenced, this assay can be adapted to examine mtDNA damage in virtually any species with little development. Therefore, we argue that this assay will be a valuable tool in assessing DNA damage in ecotoxicological studies.     

Cardioselective and cumulative oxidation of mitochondrial DNA following subchronic doxorubicin administration

We recently reported the preferential accumulation of 8-hydroxydeoxyguanosine (8OHdG) adducts in cardiac mitochondrial DNA (mtDNA) following acute intoxication of rats with doxorubicin (C.M. Palmeira et al., Biochim. Biophys. Acta, 1321 (1997) 101-106). The concentration of 8OHdG adducts decreased to control values within 2 weeks. Since conventional antineoplastic therapy entails repeated administration of small doses of doxorubicin, it was of interest to characterize the kinetics for the accumulation and repair of 8OHdG adducts in the various DNA fractions. Weekly injections of doxorubicin (2 mg/kg, i.p.) to adult male Sprague-Dawley rats caused a cumulative dose-dependent increase in the concentration of 8OHdG adducts in both mtDNA and nuclear DNA (nDNA) from heart and liver. Following six weekly injections, the concentration of 8OHdG in cardiac mtDNA was 50% higher than liver mtDNA and twice that of cardiac nDNA. In contrast to the rapid repair of 8OHdG observed during the first days following an acute intoxicating dose of doxorubicin, the concentration of 8OHdG adducts remained constant between 1 and 5 weeks following the last injection. This was true for all DNA fractions examined. The cardioselective accumulation and persistence of 8OHdG adducts to mtDNA is consistent with the implication of mitochondrial dysfunction in the cumulative and irreversible cardiotoxicity observed clinically in patients receiving doxorubicin cancer chemotherapy.