MitoRCA-seq reveals unbalanced cytocine to thymine transition in Polg mutant mice

Comparison of mitochondrial DNA sequences from whole blood and lymphoblastoid cell lines

Lymphoblastoid cell lines (LCLs) provide an unlimited source of genomic DNA for genetic studies. Here, we compared mtDNA sequence variants, heteroplasmic or homplasmic, between LCL (sequenced by mitoRCA-seq method) and whole blood samples (sequenced through whole genome sequencing approach) of the same 130 participants in the Framingham Heart Study. We applied harmonization of sequence coverages and consistent quality control to mtDNA sequences. We identified 866 variation sites in the 130 LCL samples and 666 sites in the 130 blood samples. More than 94% of the identified homoplasmies were present in both LCL and blood samples while more than 70% of heteroplasmic sites were uniquely present either in LCL or in blood samples. The LCL and whole blood samples carried a similar number of homoplasmic variants (p = 0.45) per sample while the LCL carried a greater number of heteroplasmic variants than whole blood per sample (p < 2.2e−16). Furthermore, the LCL samples tended to accumulate low level heteroplasmies (heteroplasmy level in 3–25%) than their paired blood samples (p = 0.001). These results suggest that cautions should be taken in the interpretation and comparison of findings when different tissues/cell types or different sequencing technologies are applied to obtain mtDNA sequences.

Mitophagy in the basolateral amygdala mediates increased anxiety induced by aversive social experience

Psychosocial stress is a common risk factor for anxiety disorders. The cellular mechanism for the anxiogenic effect of psychosocial stress is largely unclear. Here, we show that chronic social defeat (CSD) stress in mice causes mitochondrial impairment, which triggers the PINK1-Parkin mitophagy pathway selectively in the amygdala. This mitophagy elevation causes excessive mitochondrial elimination and consequent mitochondrial deficiency. Mitochondrial deficiency in the basolateral amygdalae (BLA) causes weakening of synaptic transmission in the BLA-BNST (bed nucleus of the stria terminalis) anxiolytic pathway and increased anxiety. The CSD-induced increase in anxiety-like behaviors is abolished in Pink1-/- and Park2-/- mice and alleviated by optogenetic activation of the BLA-BNST synapse. This study identifies an unsuspected role of mitophagy in psychogenetic-stress-induced anxiety elevation and reveals that mitochondrial deficiency is sufficient to increase anxiety and underlies the psychosocial-stress-induced anxiety increase. Mitochondria and mitophagy, therefore, can be potentially targeted to ameliorate anxiety.

The mtDNA mutation spectrum in the PolG mutator mouse reveals germline and somatic selection

BACKGROUND

Mitochondrial DNA (mtDNA) codes for products necessary for electron transport and mitochondrial gene translation. mtDNA mutations can lead to human disease and influence organismal fitness. The PolG mutator mouse lacks mtDNA proofreading function and rapidly accumulates mtDNA mutations, making it a model for examining the causes and consequences of mitochondrial mutations. Premature aging in PolG mice and their physiology have been examined in depth, but the location, frequency, and diversity of their mtDNA mutations remain understudied. Identifying the locations and spectra of mtDNA mutations in PolG mice can shed light on how selection shapes mtDNA, both within and across organisms.

RESULTS

Here, we characterized somatic and germline mtDNA mutations in brain and liver tissue of PolG mice to quantify mutation count (number of unique mutations) and frequency (mutation prevalence). Overall, mtDNA mutation count and frequency were the lowest in the D-loop, where an mtDNA origin of replication is located, but otherwise uniform across the mitochondrial genome. Somatic mtDNA mutations have a higher mutation count than germline mutations. However, germline mutations maintain a higher frequency and were also more likely to be silent. Cytosine to thymine mutations characteristic of replication errors were the plurality of basepair changes, and missense C to T mutations primarily resulted in increased protein hydrophobicity. Unlike wild type mice, PolG mice do not appear to show strand asymmetry in mtDNA mutations. Indel mutations had a lower count and frequency than point mutations and tended to be short, frameshift deletions.

CONCLUSIONS

Our results provide strong evidence that purifying selection plays a major role in the mtDNA of PolG mice. Missense mutations were less likely to be passed down in the germline, and they were less likely to spread to high frequencies. The D-loop appears to have resistance to mutations, either through selection or as a by-product of replication processes. Missense mutations that decrease hydrophobicity also tend to be selected against, reflecting the membrane-bound nature of mtDNA-encoded proteins. The abundance of mutations from polymerase errors compared with reactive oxygen species (ROS) damage supports previous studies suggesting ROS plays a minimal role in exacerbating the PolG phenotype, but our findings on strand asymmetry provide discussion for the role of polymerase errors in wild type organisms. Our results provide further insight on how selection shapes mtDNA mutations and on the aging mechanisms in PolG mice.

Mitochondrial mutations in Caenorhabditis elegans show signatures of oxidative damage and an AT-bias

Rapid mutation rates are typical of mitochondrial genomes (mtDNAs) in animals, but it is not clear why. The difficulty of obtaining measurements of mtDNA mutation that are not biased by natural selection has stymied efforts to distinguish between competing hypotheses about the causes of high mtDNA mutation rates. Several studies which have measured mtDNA mutations in nematodes have yielded small datasets with conflicting conclusions about the relative abundance of different substitution classes (i.e., the mutation spectrum). We therefore leveraged Duplex Sequencing, a high-fidelity DNA sequencing technique, to characterize de novo mtDNA mutations in Caenorhabditis elegans. This approach detected nearly an order of magnitude more mtDNA mutations than documented in any previous nematode mutation study. Despite an existing extreme AT bias in the C. elegans mtDNA (75.6% AT), we found that a significant majority of mutations increase genomic AT content. Compared to some prior studies in nematodes and other animals, the mutation spectrum reported here contains an abundance of CG→AT transversions, supporting the hypothesis that oxidative damage may be a driver of mtDNA mutations in nematodes. Furthermore, we found an excess of G→T and C→T changes on the coding DNA strand relative to the template strand, consistent with increased exposure to oxidative damage. Analysis of the distribution of mutations across the mtDNA revealed significant variation among protein-coding genes and as well as among neighboring nucleotides. This high-resolution view of mitochondrial mutations in C. elegans highlights the value of this system for understanding relationships among oxidative damage, replication error, and mtDNA mutation.

Unbiased PCR-free spatio-temporal mapping of the mtDNA mutation spectrum reveals brain region-specific responses to replication instability

Background

The accumulation of mtDNA mutations in different tissues from various mouse models has been widely studied especially in the context of mtDNA mutation-driven ageing but has been confounded by the inherent limitations of the most widely used approaches. By implementing a method to sequence mtDNA without PCR amplification prior to library preparation, we map the full unbiased mtDNA mutation spectrum across six distinct brain regions from mice.

Results

We demonstrate that ageing-induced levels of mtDNA mutations (single nucleotide variants and deletions) reach stable levels at 50 weeks of age but can be further elevated specifically in the cortex, nucleus accumbens (NAc), and paraventricular thalamic nucleus (PVT) by expression of a proof-reading-deficient mitochondrial DNA polymerase, Polg^D181A. The increase in single nucleotide variants increases the fraction of shared SNVs as well as their frequency, while characteristics of deletions remain largely unaffected. In addition, Polg^D181A also induces an ageing-dependent accumulation of non-coding control-region multimers in NAc and PVT, a feature that appears almost non-existent in wild-type mice.

Conclusions

Our data provide a novel view of the spatio-temporal accumulation of mtDNA mutations using very limited tissue input. The differential response of brain regions to a state of replication instability provides insight into a possible heterogenic mitochondrial landscape across the brain that may be involved in the ageing phenotype and mitochondria-associated disorders.

Identification of unique and shared mitochondrial DNA mutations in neurodegeneration and cancer by single-cell mitochondrial DNA structural variation sequencing (MitoSV-seq)

Background

Point mutations and structural variations (SVs) in mitochondrial DNA (mtDNA) contribute to many neurodegenerative diseases. Technical limitations and heteroplasmy, however, have impeded their identification, preventing these changes from being examined in neurons in healthy and disease states.

Methods

We have developed a high-resolution technique—Mitochondrial DNA Structural Variation Sequencing (MitoSV-seq)—that identifies all types of mtDNA SVs and single-nucleotide variations (SNVs) in single neurons and novel variations that have been undetectable with conventional techniques.

Findings

Using MitoSV-seq, we discovered SVs/SNVs in dopaminergic neurons in the Ifnar1^−/− murine model of Parkinson disease. Further, MitoSV-seq was found to have broad applicability, delivering high-quality, full-length mtDNA sequences in a species-independent manner from human PBMCs, haematological cancers, and tumour cell lines, regardless of heteroplasmy. We characterised several common SVs in haematological cancers (AML and MDS) that were linked to the same mtDNA region, MT-ND5, using only 10 cells, indicating the power of MitoSV-seq in determining single-cancer-cell ontologies. Notably, the MT-ND5 hotspot, shared between all examined cancers and Ifnar1^−/− dopaminergic neurons, suggests that its mutations have clinical value as disease biomarkers.

Interpretation

MitoSV-seq identifies disease-relevant mtDNA mutations in single cells with high resolution, rendering it a potential drug screening platform in neurodegenerative diseases and cancers.

Funding

The Lundbeck Foundation, Danish Council for Independent Research-Medicine, and European Union Horizon 2020 Research and Innovation Programme.

Concentration of mitochondrial DNA mutations by cytoplasmic transfer from platelets to cultured mouse cells

Accumulation of mutations in mitochondrial DNA (mtDNA) is thought to be responsible for mitochondrial, and other, diseases and biological phenomena, such as diabetes, cancer, neurodegenerative diseases, and aging. Mouse models may elucidate the relationship between mutations in mtDNA and these abnormalities. However, because of the difficulty of mtDNA manipulation, generation of mouse models has not sufficiently progressed to enable such studies. To overcome this difficulty and to establish a source of diverse mtDNA mutations, we here generated cultured mouse cells containing mtDNA derived from an mtDNA mutator mouse that accumulates random mtDNA mutations with age. Mutation analysis of the obtained transmitochondrial cytoplasmic hybrid cells (cybrids) revealed that the cells harbored diverse mtDNA mutations occurring at a higher frequency than in mouse tissues, and exhibited severe respiration defects that would be lethal in tissues or organs. Abnormal respiratory complex formation and high stress on the mitochondrial protein quality control system appeared to be involved in these severe respiration defects. The mutation rates of the majority of highly accumulated mutations converged to either approximately 5%, 10%, or 40%, suggesting that these mutations are linked on the respective mtDNA molecules, and mtDNA in cybrid cells likely consisted of mtDNA molecules clonally expanded from the small population of introduced mtDNAs. Thus, the linked mutations in these cybrid cells cannot be evaluated individually. In addition, mtDNA mutations homologous to confirmed pathogenic mutations in human were rarely observed in our generated cybrids. However, the transmitochondrial cybrids constitute a useful tool for concentrating pathogenic mtDNA mutations and as a source of diverse mtDNA mutations to elucidate the relationship between mtDNA mutations and diseases.

Lineage Tracing in Humans Enabled by Mitochondrial Mutations and Single-Cell Genomics

Lineage tracing provides key insights into the fate of individual cells in complex organisms. Although effective genetic labeling approaches are available in model systems, in humans, most approaches require detection of nuclear somatic mutations, which have high error rates, limited scale, and do not capture cell state information. Here, we show that somatic mutations in mtDNA can be tracked by single-cell RNA or assay for transposase accessible chromatin (ATAC) sequencing. We leverage somatic mtDNA mutations as natural genetic barcodes and demonstrate their utility as highly accurate clonal markers to infer cellular relationships. We track native human cells both in vitro and in vivo and relate clonal dynamics to gene expression and chromatin accessibility. Our approach should allow clonal tracking at a 1,000-fold greater scale than with nuclear genome sequencing, with simultaneous information on cell state, opening the way to chart cellular dynamics in human health and disease.

The mitochondrial type IB topoisomerase drives mitochondrial translation and carcinogenesis

Mitochondrial topoisomerase IB (TOP1MT) is a nuclear-encoded topoisomerase, exclusively localized to mitochondria, which resolves topological stress generated during mtDNA replication and transcription. Here, we report that TOP1MT is overexpressed in cancer tissues and demonstrate that TOP1MT deficiency attenuates tumor growth in human and mouse models of colon and liver cancer. Due to their mitochondrial dysfunction, TOP1MT-KO cells become addicted to glycolysis, which limits synthetic building blocks and energy supply required for the proliferation of cancer cells in a nutrient-deprived tumor microenvironment. Mechanistically, we show that TOP1MT associates with mitoribosomal subunits, ensuring optimal mitochondrial translation and assembly of oxidative phosphorylation complexes that are critical for sustaining tumor growth. The TOP1MT genomic signature profile, based on Top1mt-KO liver cancers, is correlated with enhanced survival of hepatocellular carcinoma patients. Our results highlight the importance of TOP1MT for tumor development, providing a potential rationale to develop TOP1MT-targeted drugs as anticancer therapies.

Deleterious mitochondrial DNA point mutations are overrepresented in Drosophila expressing a proofreading-defective DNA polymerase γ

Mitochondrial DNA (mtDNA) mutations cause severe maternally inherited syndromes and the accumulation of somatic mtDNA mutations is implicated in aging and common diseases. However, the mechanisms that influence the frequency and pathogenicity of mtDNA mutations are poorly understood. To address this matter, we created a Drosophila mtDNA mutator strain expressing a proofreading-deficient form of the mitochondrial DNA polymerase. Mutator flies have a dramatically increased somatic mtDNA mutation frequency that correlates with the dosage of the proofreading-deficient polymerase. Mutator flies also exhibit mitochondrial dysfunction, shortened lifespan, a progressive locomotor deficit, and loss of dopaminergic neurons. Surprisingly, the frequency of nonsynonymous, pathogenic, and conserved-site mutations in mutator flies exceeded predictions of a neutral mutational model, indicating the existence of a positive selection mechanism that favors deleterious mtDNA variants. We propose from these findings that deleterious mtDNA mutations are overrepresented because they selectively evade quality control surveillance or because they are amplified through compensatory mitochondrial biogenesis.

The energy needs of an animal cell are supplied by tiny organelles known as mitochondria. Each of the many mitochondria in a cell has a set of blueprints for making more mitochondria, known as mitochondrial DNA (mtDNA). As animals age, their mtDNA acquires irreversible defects called mutations, which accumulate and may cause aging. Cells can selectively destroy malfunctioning mitochondria, so we hypothesized that mitochondria with harmful mutations would be selectively destroyed. To test our theory, we created a fruit fly strain with a high mtDNA mutation rate. Our hypothesis predicts that, because mitochondria bearing harmful mtDNA mutations would be destroyed, we should detect primarily less harmful mutations in our strain. However, the mtDNA mutations we detected were more harmful than expected by chance. We suggest two possible explanations: First, mitochondria with harmful mtDNA mutations may be degraded less often because they generate little energy and are not damaged by toxic byproducts of energy production. Second, cells may compensate for harmful mtDNA mutations by stimulating mitochondria to multiply, creating more healthy mitochondria but also more mitochondria with harmful mtDNA mutations. Future studies will distinguish between these models and further advance our understanding of aging and aging related disease.

Calorie restriction does not influence oocyte quality in oocytes from POLG mitochondrial mutator mice

It has recently been demonstrated that moderate adult onset caloric restriction (e.g. calorie restriction; CR) has a positive impact on female fertility in aged mice, due in large to preventing the age-associated decline in the quality of oocytes available for fertilization. The impact of CR on oocyte quality has been attributed, at least in part, to mitochondrial functions. In mitochondrial DNA (mtDNA) mutator mice (Polg^D257A/D257A), which harbor a mutation in the proofreading mtDNA polymerase-gamma (POLG), mitochondrial mutations rapidly accumulate, resulting in a premature aging phenotype and female infertility. As CR has been shown to extend both lifespan and ‘healthspan’ as well as improve oocyte quality in aged mice, we investigated whether adult onset CR could improve oocyte quality in the POLG mouse. Female Polg^D257A/D257A mice exhibited infertility based on an inability to produce litters through natural mating. Analysis of oocytes from 8–9-month-old Polg^D257A/D257A mice on CR following hormone stimulation revealed no improvement in the number of oocytes ovulated. Furthermore, CR did not result in a greater percentage of metaphase II oocytes, with the majority of the oocytes prematurely arrested at the germinal vesicle stage. Finally, CR did not improve the abnormal mitochondrial distribution or pronounced defects in meiotic spindle assembly and chromosomal distribution observed in the ad libitum fed Polg^D257A/D257A. Taken together, these data suggest that although CR benefits oocyte quality and fertility outcomes in naturally aged female mice, it does not sufficiently improve oocyte quality in Polg^D257A/D257A.

Deep sequencing of the mitochondrial genome reveals common heteroplasmic sites in NADH dehydrogenase genes

Increasing evidence implicates mitochondrial dysfunction in aging and age-related conditions. But little is known about the molecular basis for this connection. A possible cause may be mutations in the mitochondrial DNA (mtDNA), which are often heteroplasmic-the joint presence of different alleles at a single locus in the same individual. However, the involvement of mtDNA heteroplasmy in aging and age-related conditions has not been investigated thoroughly. We deep-sequenced the complete mtDNA genomes of 356 Framingham Heart Study participants (52% women, mean age 43, mean coverage 4570-fold), identified 2880 unique mutations and comprehensively annotated them by MITOMAP and PolyPhen-2. We discovered 11 heteroplasmic "hot" spots [NADH dehydrogenase (ND) subunit 1, 4, 5 and 6 genes, n = 7; cytochrome c oxidase I (COI), n = 2; 16S rRNA, n = 1; D-loop, n = 1] for which the alternative-to-reference allele ratios significantly increased with advancing age (Bonferroni correction p < 0.001). Four of these heteroplasmic mutations in ND and COI genes were predicted to be deleterious nonsynonymous mutations which may have direct impact on ATP production. We confirmed previous findings that healthy individuals carry many low-frequency heteroplasmy mutations with potentially deleterious effects. We hypothesize that the effect of a single deleterious heteroplasmy may be minimal due to a low mutant-to-wildtype allele ratio, whereas the aggregate effects of many deleterious mutations may cause changes in mitochondrial function and contribute to age-related diseases. The identification of age-related mtDNA mutations is an important step to understand the genetic architecture of age-related diseases and may uncover novel therapeutic targets for such diseases.

Proliferation Cycle Causes Age Dependent Mitochondrial Deficiencies and Contributes to the Aging of Stem Cells

In addition to chronological aging, stem cells are also subject to proliferative aging during the adult life span. However, the consequences of proliferative cycle and their contributions to stem cells aging have not been well investigated. Using Drosophila female germ line stem cells as a model, we found that the replication cycle leads to the age dependent decline of female fecundity, and is a major factor causing developmental abnormalities in the progeny of old females. The proliferative aging does not cause telomere shortening, but causes an accumulation of mitochondrial DNA (mtDNA) mutations or rearrangements at the control region. We propose that damaging mutations on mtDNA caused by accumulation of proliferation cycles in aged stem cells may disrupt mitochondrial respiration chain and impair mtDNA replication and represent a conserved mechanism underlying stem cell aging.