Mitochondrial DNA Content Is Linked to Cardiovascular Disease Patient Phenotypes.

Background We sought to determine whether mitochondrial DNA (mtDNA) content can be used as markers for 12 key phenotypes among cardiovascular disease patients, and whether these markers are valid across patients with diverse ancestries. Methods and Results DNA was collected from the peripheral blood of 996 cardiovascular disease patients at the Cleveland Clinic. The mtDNA copy number and DNA-level variation were assessed from whole-genome sequence. Patients were also ascertained retrospectively for histories of 10 clinical events, as well as for maximum stenosis and extent of disease at baseline. Self-reported race and maternal ancestry inferred from mtDNA sequence were recorded. MtDNA copy number and overall mtDNA rare variant load were significantly lower in patients with histories of various adverse clinical events, and mtDNA copy number was inversely correlated with extent of disease. Strong associations were also found between absence of rare variants in the genes MT-ATP6 and MT-COII and patient histories of hyperlipidemia and myocardial infarction, respectively. Importantly, associations were not ancestry dependent. Conclusions This study provides evidence that mtDNA copy number in circulation is associated with a variety of cardiovascular disease patient phenotypes. Results also suggest a protective role for some rare inherited mtDNA variants. Overall, the study supports the potential of mtDNA content and abundance as biomarkers in heart disease, in a manner that is valid across diverse ancestries.

Detection and quantification of mitochondrial DNA deletions from next-generation sequence data.

BACKGROUND: Chromosomal deletions represent an important class of human genetic variation. Various methods have been developed to mine "next-generation" sequencing (NGS) data to detect deletions and quantify their clonal abundances. These methods have focused almost exclusively on the nuclear genome, ignoring the mitochondrial chromosome (mtDNA). Detecting mtDNA deletions requires special care. First, the chromosome's relatively small size (16,569 bp) necessitates the ability to detect extremely focal events. Second, the chromosome can be present at thousands of copies in a single cell (in contrast to two copies of nuclear chromosomes), and mtDNA deletions may be present on only a very small percentage of chromosomes. Here we present a method, termed MitoDel, to detect mtDNA deletions from NGS data. RESULTS: We validate the method on simulated and real data, and show that MitoDel can detect novel and previously-reported mtDNA deletions. We establish that MitoDel can find deletions such as the "common deletion" at heteroplasmy levels well below 1%. CONCLUSIONS: MitoDel is a tool for detecting large mitochondrial deletions at low heteroplasmy levels. The tool can be downloaded at http://mendel.gene.cwru.edu/laframboiselab/ .

Phylogenetic sequence analysis and functional studies reveal compensatory amino acid substitutions in loop 2 of human ribonucleotide reductase.

Eukaryotic class I ribonucleotide reductases (RRs) generate deoxyribonucleotides for DNA synthesis. Binding of dNTP effectors is coupled to the formation of active dimers and induces conformational changes in a short loop (loop 2) to regulate RR specificity among its nucleoside diphosphate substrates. Moreover, ATP and dATP bind at an additional allosteric site 40 Å away from loop 2 and thereby drive formation of activated or inactive hexamers, respectively. To better understand how dNTP binding influences specificity, activity, and oligomerization of human RR, we aligned >300 eukaryotic RR sequences to examine natural sequence variation in loop 2. We found that most amino acids in eukaryotic loop 2 were nearly invariant in this sample; however, two positions co-varied as nonconservative substitutions (N291G and P294K; human numbering). We also found that the individual N291G and P294K substitutions in human RR additively affect substrate specificity. The P294K substitution significantly impaired effector-induced oligomerization required for enzyme activity, and oligomerization was rescued in the N291G/P294K enzyme. None of the other mutants exhibited altered ATP-mediated hexamerization; however, certain combinations of loop 2 mutations and dNTP effectors perturbed ATP's role as an allosteric activator. Our results demonstrate that the observed compensatory covariation of amino acids in eukaryotic loop 2 is essential for its role in dNTP-induced dimerization. In contrast, defects in substrate specificity are not rescued in the double mutant, implying that functional sequence variation elsewhere in the protein is necessary. These findings yield insight into loop 2's roles in regulating RR specificity, allostery, and oligomerization.

Heteroplasmic shifts in tumor mitochondrial genomes reveal tissue-specific signals of relaxed and positive selection.

Although mitochondrial genomes (mtDNA) accumulate elevated levels of mutations in cancer cells, the origin and functional impact of these mutations remain controversial. Here, we queried whole-genome sequence data from 1,916 patients across 24 cancer types to characterize patterns of mtDNA mutations and elucidate the selective constraints driving their fate. Given that mitochondrial genomes are polyploid, cells with advantageous levels of mtDNA mutations can be selected for depending on their cellular environment. Therefore, we tracked changes in per-cell abundances of mtDNA mutations from normal to tumor cells in the same patient. Tumor mitochondrial genomes show distinct mutational patterns and are disproportionately enriched for protein-altering changes. Moreover, protein-altering mtDNA variants that are initially present at low frequencies in normal cells preferentially expand in the altered tumor environment, suggesting selective advantage. We also perform these analyses with attention to the cancer's tissue of origin, which revealed tissue-specific differences in selective signals. The mitochondrial genomes in renal chromophobe and thyroid cancers show particularly strong signals of positive selection, indicated by higher proportions and per-cell abundances of truncating variants. Dramatic tumor- and tissue-specific variations in selective pressures suggest that cancer cells with advantageous levels of damaged mitochondrial genomes will selectively proliferate to facilitate the tumorigenic process.

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

[This corrects the article DOI: 10.1371/journal.pone.0139253.].

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.