Assessing Mitochondrial DNA Variation and Copy Number in Lymphocytes of ~2,000 Sardinians Using Tailored Sequencing Analysis Tools

Approaches and challenges in identifying, quantifying, and manipulating dynamic mitochondrial genome variations

Mitochondria, the cellular powerhouses, possess their own unique genetic system, including replication, transcription, and translation. Studying these processes is crucial for comprehending mitochondrial disorders, energy production, and their related diseases. Over the past decades, various approaches have been applied in detecting and quantifying mitochondrial genome variations with also the purpose of manipulation of mitochondria or mitochondrial genome for therapeutics. Understanding the scope and limitations of above strategies is not only fundamental to the understanding of basic biology but also critical for exploring disease-related novel target(s), as well to develop innovative therapies. Here, this review provides an overview of different tools and techniques for accurate mitochondrial genome variations identification, quantification, and discuss novel strategies for the manipulation of mitochondria to develop innovative therapeutic interventions, through combining the insights gained from the study of mitochondrial genetics with ongoing single cell omics combined with advanced single molecular tools.

Subtype-Specific Association of Mitochondrial DNA Copy Number With Poststroke/TIA Outcomes in 10 241 Patients in China

BACKGROUND

Mitochondrial DNA copy number (mtDNA-CN) is associated with the severity and mortality in patients with stroke, but the associations in different stroke subtypes remain unexplored.

METHODS

We conducted an observational prospective cohort analysis on patients with ischemic stroke or transient ischemic attack enrolled in the Third China National Stroke Registry. We applied logistic models to assess the association of mtDNA-CN with functional outcome (modified Rankin Scale score, 3-6 versus 0-2) and Cox proportional hazard models to assess the association with stroke recurrence (treating mortality as a competing risk) and mortality during a 12-month follow-up, adjusting for sex, age, physical activity, National Institutes of Health Stroke Scale at admission, history of stroke and peripheral artery disease, small artery occlusion, and interleukin-6. Subgroup analyses stratified by age and stroke subtypes were conducted.

RESULTS

The Third China National Stroke Registry enrolled 15 166 patients, of which 10 241 with whole-genome sequencing data were retained (mean age, 62.2 [SD, 11.2] years; 68.8% men). The associations between mtDNA-CN and poststroke/transient ischemic attack outcomes were specific to patients aged ≤65 years, with lower mtDNA-CN significantly associated with stroke recurrence in 12 months (subdistribution hazard ratio, 1.15 per SD lower mtDNA-CN [95% CI, 1.04-1.27]; P=5.2×10-3) and higher all-cause mortality in 3 months (hazard ratio, 2.19 [95% CI, 1.41-3.39]; P=5.0×10-4). Across subtypes, the associations of mtDNA-CN with stroke recurrence were specific to stroke of undetermined cause (subdistribution hazard ratio, 1.28 [95% CI, 1.11-1.48]; P=6.6×10-4). In particular, lower mtDNA-CN was associated with poorer functional outcomes in stroke of undetermined cause patients diagnosed with embolic stroke of undetermined source (odds ratio, 1.53 [95% CI, 1.20-1.94]; P=5.4×10-4), which remained significant after excluding patients with recurrent stroke (odds ratio, 1.49 [95% CI, 1.14-1.94]; P=3.0×10-3).

CONCLUSIONS

Lower mtDNA-CN is associated with higher stroke recurrence rate and all-cause mortality, as well as poorer functional outcome at follow-up, among stroke of undetermined cause, embolic stroke of undetermined source, and younger patients.

Associations between mitochondrial copy number, exercise capacity, physiologic cost of walking, and cardiac strain in young adult survivors of childhood cancer

PURPOSE

Childhood cancer survivors are at risk for cardiac dysfunction and impaired physical performance, though underlying cellular mechanisms are not well studied. In this cross-sectional study, we examined the association between peripheral blood mitochondrial DNA copy number (mtDNA-CN, a proxy for mitochondrial function) and markers of performance impairment and cardiac dysfunction.

METHODS

Whole-genome sequencing, validated by quantitative polymerase chain reaction, was used to estimate mtDNA-CN in 1720 adult survivors of childhood cancer (48.5% female; mean age = 30.7 years, standard deviation (SD) = 9.0). Multivariable logistic regression was performed to evaluate the associations between mtDNA-CN and exercise intolerance, walking inefficiency, and abnormal global longitudinal strain (GLS), adjusting for treatment exposures, age, sex, and race and ethnicity.

RESULTS

The prevalence of exercise intolerance, walking inefficiency, and abnormal GLS among survivors was 25.7%, 10.7%, and 31.7%, respectively. Each SD increase of mtDNA-CN was associated with decreased odds of abnormal GLS (adjusted odds ratio (OR) = 0.88, p = 0.04) but was not associated with exercise intolerance (OR = 1.02, p = 0.76) or walking inefficiency (OR = 1.06, p = 0.46). Alkylating agent exposure was associated with increased odds of exercise intolerance (OR = 2.25, p < 0.0001), walking inefficiency (OR = 2.37, p < 0.0001), and abnormal GLS (OR = 1.78, p = 0.0002).

CONCLUSIONS

Increased mtDNA-CN is associated with decreased odds of abnormal cardiac function in childhood cancer survivors.

IMPLICATIONS FOR CANCER SURVIVORS

These findings demonstrate a potential role for mtDNA-CN as a biomarker of early cardiac dysfunction in this population.

Mitochondrial stress and aging: Lessons from C. elegans

Aging is accompanied by a progressive decline in mitochondrial function, which in turn contributes to a variety of age-related diseases. Counterintuitively, a growing number of studies have found that disruption of mitochondrial function often leads to increased lifespan. This seemingly contradictory observation has inspired extensive research into genetic pathways underlying the mitochondrial basis of aging, particularly within the model organism Caenorhabditis elegans. The complex and antagonistic roles of mitochondria in the aging process have altered the view of mitochondria, which not only serve as simple bioenergetic factories but also as signaling platforms for the maintenance of cellular homeostasis and organismal health. Here, we review the contributions of C. elegans to our understanding of mitochondrial function in the aging process over the past decades. In addition, we explore how these insights may promote future research of mitochondrial-targeted strategies in higher organisms to potentially slow aging and delay age-related disease progression.

A naturally occurring variant of SHLP2 is a protective factor in Parkinson's disease

Mitochondrial DNA single nucleotide polymorphisms (mtSNPs) have been associated with a reduced risk of developing Parkinson's disease (PD), yet the underlying mechanisms remain elusive. In this study, we investigate the functional role of a PD-associated mtSNP that impacts the mitochondrial-derived peptide (MDP) Small Humanin-like Peptide 2 (SHLP2). We identify m.2158 T > C, a mtSNP associated with reduced PD risk, within the small open reading frame encoding SHLP2. This mtSNP results in an alternative form of SHLP2 (lysine 4 replaced with arginine; K4R). Using targeted mass spectrometry, we detect specific tryptic fragments of SHLP2 in neuronal cells and demonstrate its binding to mitochondrial complex 1. Notably, we observe that the K4R variant, associated with reduced PD risk, exhibits increased stability compared to WT SHLP2. Additionally, both WT and K4R SHLP2 show enhanced protection against mitochondrial dysfunction in in vitro experiments and confer protection against a PD-inducing toxin, a mitochondrial complex 1 inhibitor, in a mouse model. This study sheds light on the functional consequences of the m.2158 T > C mtSNP on SHLP2 and provides insights into the potential mechanisms by which this mtSNP may reduce the risk of PD.

Association analysis of mitochondrial DNA heteroplasmic variants: methods and application

We rigorously assessed a comprehensive association testing framework for heteroplasmy, employing both simulated and real-world data. This framework employed a variant allele fraction (VAF) threshold and harnessed multiple gene-based tests for robust identification and association testing of heteroplasmy. Our simulation studies demonstrated that gene-based tests maintained an appropriate type I error rate at α=0.001. Notably, when 5% or more heteroplasmic variants within a target region were linked to an outcome, burden-extension tests (including the adaptive burden test, variable threshold burden test, and z-score weighting burden test) outperformed the sequence kernel association test (SKAT) and the original burden test. Applying this framework, we conducted association analyses on whole-blood derived heteroplasmy in 17,507 individuals of African and European ancestries (31% of African Ancestry, mean age of 62, with 58% women) with whole genome sequencing data. We performed both cohort- and ancestry-specific association analyses, followed by meta-analysis on both pooled samples and within each ancestry group. Our results suggest that mtDNA-encoded genes/regions are likely to exhibit varying rates in somatic aging, with the notably strong associations observed between heteroplasmy in the RNR1 and RNR2 genes (p<0.001) and advance aging by the Original Burden test. In contrast, SKAT identified significant associations (p<0.001) between diabetes and the aggregated effects of heteroplasmy in several protein-coding genes. Further research is warranted to validate these findings. In summary, our proposed statistical framework represents a valuable tool for facilitating association testing of heteroplasmy with disease traits in large human populations.

Mitochondrial DNA Copy Number Variation in Asthma Risk, Severity, and Exacerbations

Rationale

Although airway oxidative stress and inflammation are central to asthma pathogenesis, there is limited knowledge of the relationship of asthma risk, severity, or exacerbations to mitochondrial dysfunction, which is pivotal to oxidant generation and inflammation.

Objectives

We investigated whether mitochondrial DNA copy number (mtDNA-CN) as a measure of mitochondrial function is associated with asthma diagnosis, severity, oxidative stress, and exacerbations.

Methods

We measured mtDNA-CN in blood in two cohorts. In the UK Biobank (UKB), we compared mtDNA-CN in mild and moderate-severe asthmatics to non-asthmatics. In the Severe Asthma Research Program (SARP), we evaluated mtDNA-CN in relation to asthma severity, biomarkers of oxidative stress and inflammation, and exacerbations.

Measures and Main Results

In UK Biobank, asthmatics (n = 29,768) have lower mtDNA-CN compared to non-asthmatics (n = 239,158) (beta, -0.026 [95% CI, -0.038 to -0.014], P = 2.46×10-5). While lower mtDNA-CN is associated with asthma, mtDNA-CN did not differ by asthma severity in either UKB or SARP. Biomarkers of inflammation show that asthmatics have higher white blood cells (WBC), neutrophils, eosinophils, fraction exhaled nitric oxide (FENO), and lower superoxide dismutase (SOD) than non-asthmatics, confirming greater oxidative stress in asthma. In one year follow-up in SARP, higher mtDNA-CN is associated with reduced risk of three or more exacerbations in the subsequent year (OR 0.352 [95% CI, 0.164 to 0.753], P = 0.007).

Conclusions

Asthma is characterized by mitochondrial dysfunction. Higher mtDNA-CN identifies an exacerbation-resistant asthma phenotype, suggesting mitochondrial function is important in exacerbation risk.

Plausible Role of Mitochondrial DNA Copy Number in Neurodegeneration-a Need for Therapeutic Approach in Parkinson's Disease (PD)

Parkinson's disease (PD) is an advancing age-associated progressive brain disorder which has various diverse factors, among them mitochondrial dysfunction involves in dopaminergic (DA) degeneration. Aging causes a rise in mitochondrial abnormalities which leads to structural and functional modifications in neuronal activity and cell death in PD. This ends in deterioration of mitochondrial function, mitochondrial alterations, mitochondrial DNA copy number (mtDNA CN) and oxidative phosphorylation (OXPHOS) capacity. mtDNA levels or mtDNA CN in PD have reported that mtDNA depletion would be a predisposing factor in PD pathogenesis. To maintain the mtDNA levels, therapeutic approaches have been focused on mitochondrial biogenesis in PD. The depletion of mtDNA levels in PD can be influenced by autophagic dysregulation, apoptosis, neuroinflammation, oxidative stress, sirtuins, and calcium homeostasis. The current review describes the regulation of mtDNA levels and discusses the plausible molecular pathways in mtDNA CN depletion in PD pathogenesis. We conclude by suggesting further research on mtDNA depletion which might show a promising effect in predicting and diagnosing PD.

Inclusion body myositis with early onset: a population-based study

INTRODUCTION

Inclusion body myositis (IBM), an inflammatory myopathy with progressive weakness without efficient treatment, typically presents after 45 years of age and younger patients are sparsely studied.

METHODS

In a population-based study during a 33-year period, 142 patients with IBM were identified in western Sweden. Six patients fell outside the European Neuromuscular Centre 2011 criteria for IBM due to young age at symptom onset, verified by a muscle biopsy < 50 years of age. These were defined as early-onset IBM and included in this study. Medical records, muscle strength, comorbidities, muscle biopsies, and nuclear- and mitochondrial DNA were examined and compared with patients with IBM and age matched controls from the same population.

RESULTS

The median age at symptom onset was 36 (range 34-45) years and at diagnosis 43 (range 38-58) years. Four patients were deceased at a median age of 59 (range 50-75) years. The median survival from diagnosis was 14 (range 10-18) years. The prevalence December 31 2017 was 1.2 per million inhabitants and the mean incidence 0.12 patients per million inhabitants and year. The mean decline in quadriceps strength ± 1 standard deviation was 1.21 ± 0.2 Newton or 0.91 ± 0.2% per month and correlated to time from diagnosis (p < 0.001). Five patients had swallowing difficulties. All patients displayed mitochondrial changes in muscle including cytochrome c oxidase deficiency and the mitochondrial DNA mutation load was high.

CONCLUSIONS

Early-onset IBM is a severe disease, causing progressive muscle weakness, high muscle mitochondrial DNA mutation load and a reduced cumulative survival in young and middle-aged individuals.

Association Between Whole Blood-Derived Mitochondrial DNA Copy Number, Low-Density Lipoprotein Cholesterol, and Cardiovascular Disease Risk

Background The relationship between mitochondrial DNA copy number (mtDNA CN) and cardiovascular disease remains elusive. Methods and Results We performed cross-sectional and prospective association analyses of blood-derived mtDNA CN and cardiovascular disease outcomes in 27 316 participants in 8 cohorts of multiple racial and ethnic groups with whole-genome sequencing. We also performed Mendelian randomization to explore causal relationships of mtDNA CN with coronary heart disease (CHD) and cardiometabolic risk factors (obesity, diabetes, hypertension, and hyperlipidemia). P<0.01 was used for significance. We validated most of the previously reported associations between mtDNA CN and cardiovascular disease outcomes. For example, 1-SD unit lower level of mtDNA CN was associated with 1.08 (95% CI, 1.04-1.12; P<0.001) times the hazard for developing incident CHD, adjusting for covariates. Mendelian randomization analyses showed no causal effect from a lower level of mtDNA CN to a higher CHD risk (β=0.091; P=0.11) or in the reverse direction (β=-0.012; P=0.076). Additional bidirectional Mendelian randomization analyses revealed that low-density lipoprotein cholesterol had a causal effect on mtDNA CN (β=-0.084; P<0.001), but the reverse direction was not significant (P=0.059). No causal associations were observed between mtDNA CN and obesity, diabetes, and hypertension, in either direction. Multivariable Mendelian randomization analyses showed no causal effect of CHD on mtDNA CN, controlling for low-density lipoprotein cholesterol level (P=0.52), whereas there was a strong direct causal effect of higher low-density lipoprotein cholesterol on lower mtDNA CN, adjusting for CHD status (β=-0.092; P<0.001). Conclusions Our findings indicate that high low-density lipoprotein cholesterol may underlie the complex relationships between mtDNA CN and vascular atherosclerosis.

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans

Mitochondrial DNA (mtDNA) is a maternally inherited, high-copy-number genome required for oxidative phosphorylation1. Heteroplasmy refers to the presence of a mixture of mtDNA alleles in an individual and has been associated with disease and ageing. Mechanisms underlying common variation in human heteroplasmy, and the influence of the nuclear genome on this variation, remain insufficiently explored. Here we quantify mtDNA copy number (mtCN) and heteroplasmy using blood-derived whole-genome sequences from 274,832 individuals and perform genome-wide association studies to identify associated nuclear loci. Following blood cell composition correction, we find that mtCN declines linearly with age and is associated with variants at 92 nuclear loci. We observe that nearly everyone harbours heteroplasmic mtDNA variants obeying two principles: (1) heteroplasmic single nucleotide variants tend to arise somatically and accumulate sharply after the age of 70 years, whereas (2) heteroplasmic indels are maternally inherited as mixtures with relative levels associated with 42 nuclear loci involved in mtDNA replication, maintenance and novel pathways. These loci may act by conferring a replicative advantage to certain mtDNA alleles. As an illustrative example, we identify a length variant carried by more than 50% of humans at position chrM:302 within a G-quadruplex previously proposed to mediate mtDNA transcription/replication switching2,3. We find that this variant exerts cis-acting genetic control over mtDNA abundance and is itself associated in-trans with nuclear loci encoding machinery for this regulatory switch. Our study suggests that common variation in the nuclear genome can shape variation in mtCN and heteroplasmy dynamics across the human population.

Energizing Mitochondria to Prevent Mobility Loss in Aging: Rationale and Hypotheses

Based on recent studies from our group and others, we hypothesize that mitochondrial dysfunction during aging may be the root cause of mobility decline through deficits in the musculoskeletal and central nervous systems. Mitochondrial dysfunction could be a therapeutic target to prevent mobility decline in aging.

Association of Mitochondrial DNA Copy Number With Brain MRI Markers and Cognitive Function: A Meta-analysis of Community-Based Cohorts

BACKGROUND AND OBJECTIVES

Previous studies suggest that lower mitochondrial DNA (mtDNA) copy number (CN) is associated with neurodegenerative diseases. However, whether mtDNA CN in whole blood is related to endophenotypes of Alzheimer disease (AD) and AD-related dementia (AD/ADRD) needs further investigation. We assessed the association of mtDNA CN with cognitive function and MRI measures in community-based samples of middle-aged to older adults.

METHODS

We included dementia-free participants from 9 diverse community-based cohorts with whole-genome sequencing in the Trans-Omics for Precision Medicine (TOPMed) program. Circulating mtDNA CN was estimated as twice the ratio of the average coverage of mtDNA to nuclear DNA. Brain MRI markers included total brain, hippocampal, and white matter hyperintensity volumes. General cognitive function was derived from distinct cognitive domains. We performed cohort-specific association analyses of mtDNA CN with AD/ADRD endophenotypes assessed within ±5 years (i.e., cross-sectional analyses) or 5-20 years after blood draw (i.e., prospective analyses) adjusting for potential confounders. We further explored associations stratified by sex and age (<60 vs ≥60 years). Fixed-effects or sample size-weighted meta-analyses were performed to combine results. Finally, we performed mendelian randomization (MR) analyses to assess causality.

RESULTS

We included up to 19,152 participants (mean age 59 years, 57% women). Higher mtDNA CN was cross-sectionally associated with better general cognitive function (β = 0.04; 95% CI 0.02-0.06) independent of age, sex, batch effects, race/ethnicity, time between blood draw and cognitive evaluation, cohort-specific variables, and education. Additional adjustment for blood cell counts or cardiometabolic traits led to slightly attenuated results. We observed similar significant associations with cognition in prospective analyses, although of reduced magnitude. We found no significant associations between mtDNA CN and brain MRI measures in meta-analyses. MR analyses did not reveal a causal relation between mtDNA CN in blood and cognition.

DISCUSSION

Higher mtDNA CN in blood is associated with better current and future general cognitive function in large and diverse communities across the United States. Although MR analyses did not support a causal role, additional research is needed to assess causality. Circulating mtDNA CN could serve nevertheless as a biomarker of current and future cognitive function in the community.

Biomarkers of aging

Aging biomarkers are a combination of biological parameters to (i) assess age-related changes, (ii) track the physiological aging process, and (iii) predict the transition into a pathological status. Although a broad spectrum of aging biomarkers has been developed, their potential uses and limitations remain poorly characterized. An immediate goal of biomarkers is to help us answer the following three fundamental questions in aging research: How old are we? Why do we get old? And how can we age slower? This review aims to address this need. Here, we summarize our current knowledge of biomarkers developed for cellular, organ, and organismal levels of aging, comprising six pillars: physiological characteristics, medical imaging, histological features, cellular alterations, molecular changes, and secretory factors. To fulfill all these requisites, we propose that aging biomarkers should qualify for being specific, systemic, and clinically relevant.

Association of telomere length and mitochondrial DNA copy number, two biomarkers of biological aging, with the risk of venous thromboembolism

BACKGROUND

Venous thromboembolism (VTE) is the third most common cardiovascular disease and occurs in all age groups, albeit the risk increases considerably with age. Previous research indicates mitochondrial dysfunction and telomere shortening in cardiovascular aging. However, in the context of VTE this has not been investigated in detail.

AIM

We aimed to explore biomarkers reflecting biological aging (i.e. human mitochondrial DNA copy number (mtDNA) and telomere length) and their association with VTE.

METHODS

mtDNA and telomere length were measured in a case-control study of 116 patients with a history of VTE and 128 age- and sex-matched healthy individuals from isolated blood using a qPCR-based assay kit. Cases had at least one unprovoked VTE event and were enrolled no earlier than 3 months after the last VTE event.

RESULTS

The mtDNA copy number was significantly lower in VTE cases compared to controls (median [IQR]: 663 per diploid cells [78.75-2204.5] vs. 2832 per diploid cells [724-4350]; p < 0.001). After adjustment for age, sex, BMI, and smoking, mtDNA copy number was independently associated with VTE risk (odds ratio per increase in 400 mtDNA per diploid cell: 0.889, 95%CI 0.834-0.947). mtDNA copy numbers were significantly different between women and men (2375 [455-3737] women vs. 893 [152-3154] men; p < 0.001). The analysis of telomere length showed no significant difference between patients and healthy controls.

CONCLUSION

Lower mtDNA levels were found in patients with VTE compared to controls, indicating an association of biological aging with risk of VTE.

Mitochondrial DNA Copy Number Is a Potential Biomarker for Treatment Choice Between Metformin and Acarbose

Metformin is the first-line drug for type 2 diabetes (T2D) while acarbose is suggested as a viable alternative in Chinese patients with newly diagnosed T2D. However, few biomarkers have been established to guide the choice between these two agents. Mitochondrial DNA (mtDNA) copy number (mtDNA-CN) is a biomarker of mitochondrial function, which is associated with various metabolic outcomes. Using data from the trial of Metformin and Acarbose in Chinese as the Initial Hypoglycaemic Treatment (MARCH) (metformin n = 214; acarbose n = 198), we examined whether mtDNA-CN was associated with response to the drugs in terms of glycemic response and β-cell function protection response. The glycemic response is defined as the maximum glucose reduction of glycated hemoglobin A_1c , fasting plasma glucose, or postprandial blood glucose during 48 weeks. β-cell function protection response is defined as the maximum increment of insulinogenic index (IGI) or disposition index (DI). For all three glycemic responses, mtDNA-CN was not significantly associated with either metformin or acarbose. Importantly, for β-cell function protection response, we found the increased mtDNA-CN was significantly associated with more IGI increment (beta: 0.84; 95% confidence interval (CI), 0.02 to 1.66) in the metformin group, but less IGI increment (beta: -1.38; 95% CI, -2.52 to -0.23) in the acarbose group. A significant interaction (P = 0.008) between mtDNA-CN and the treatment group was observed. Consistent results were also obtained when DI increment was used as a measure of β-cell function response. This study demonstrated the potential application of mtDNA-CN in guiding the treatment choice between metformin and acarbose based on β-cell protection.

Higher mitochondrial DNA copy number is associated with metformin-induced weight loss

BACKGROUND

Considerable variability exists in response to metformin with few effective biomarkers to guide the treatment. Here we evaluated whether whole blood derived mitochondrial DNA copy number (mtDNA-CN) is a biomarker of metformin response as measured by glucose reduction or weight loss.

METHODS

Using data from the trial of Metformin (n = 304) and AcaRbose (n = 300) in Chinese as the initial Hypoglycaemic treatment (MARCH), we examined the association between mtDNA-CN and two metformin response outcomes of HbA1c reduction and weight loss. The acarbose arm was used as a comparator group. Whole blood mtDNA-CN was estimated by deep whole genome sequencing with adjustments for confounders. Multiple linear regression and repeated measurement analyses were used to evaluate the association between mtDNA-CN and drug response outcomes.

RESULTS

Here we show that glucose reduction is not significantly associated with mtDNA-CN and in either treatment arm. In the metformin arm, each increase of 1 SD in mtDNA-CN is significantly (P = 0.006) associated with a 0.43 kg more weight loss. Repeated measurement analysis shows that after 16 weeks of metformin monotherapy, patients in the top tertile of mtDNA-CN consistently lost 1.21 kg more weight than those in the bottom tertile (P < 0.001). In comparison, mtDNA-CN is not significantly associated with acarbose-induced weight loss.

CONCLUSIONS

Patients with higher mtDNA-CN are likely to lose more weight upon metformin treatment, suggesting mtDNA-CN as a potential novel biomarker for more effective weight management in type 2 diabetes.

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans

Human mitochondria contain a high copy number, maternally transmitted genome (mtDNA) that encodes 13 proteins required for oxidative phosphorylation. Heteroplasmy arises when multiple mtDNA variants co-exist in an individual and can exhibit complex dynamics in disease and in aging. As all proteins involved in mtDNA replication and maintenance are nuclear-encoded, heteroplasmy levels can, in principle, be under nuclear genetic control, however this has never been shown in humans. Here, we develop algorithms to quantify mtDNA copy number (mtCN) and heteroplasmy levels using blood-derived whole genome sequences from 274,832 individuals of diverse ancestry and perform GWAS to identify nuclear loci controlling these traits. After careful correction for blood cell composition, we observe that mtCN declines linearly with age and is associated with 92 independent nuclear genetic loci. We find that nearly every individual carries heteroplasmic variants that obey two key patterns: (1) heteroplasmic single nucleotide variants are somatic mutations that accumulate sharply after age 70, while (2) heteroplasmic indels are maternally transmitted as mtDNA mixtures with resulting levels influenced by 42 independent nuclear loci involved in mtDNA replication, maintenance, and novel pathways. These nuclear loci do not appear to act by mtDNA mutagenesis, but rather, likely act by conferring a replicative advantage to specific mtDNA molecules. As an illustrative example, the most common heteroplasmy we identify is a length variant carried by >50% of humans at position m.302 within a G-quadruplex known to serve as a replication switch. We find that this heteroplasmic variant exerts cis -acting genetic control over mtDNA abundance and is itself under trans -acting genetic control of nuclear loci encoding protein components of this regulatory switch. Our study showcases how nuclear haplotype can privilege the replication of specific mtDNA molecules to shape mtCN and heteroplasmy dynamics in the human population.

Whole-exome sequencing in 415,422 individuals identifies rare variants associated with mitochondrial DNA copy number

Inter-individual variation in the number of copies of the mitochondrial genome, called mitochondrial DNA copy number (mtDNA-CN), reflects mitochondrial function and has been associated with various aging-related diseases. We examined 415,422 exomes of self-reported White ancestry individuals from the UK Biobank and tested the impact of rare variants, at the level of single variants and through aggregate variant-set tests, on mtDNA-CN. A survey across nine variant sets tested enrichment of putatively causal variants and identified 14 genes at experiment-wide significance and three genes at marginal significance. These included associations at known mtDNA depletion syndrome genes (mtDNA helicase TWNK, p = 1.1 × 10-30; mitochondrial transcription factor TFAM, p = 4.3 × 10-15; mtDNA maintenance exonuclease MGME1, p = 2.0 × 10-6) and the V617F dominant gain-of-function mutation in the tyrosine kinase JAK2 (p = 2.7 × 10-17), associated with myeloproliferative disease. Novel genes included the ATP-dependent protease CLPX (p = 8.4 × 10-9), involved in mitochondrial proteome quality, and the mitochondrial adenylate kinase AK2 (p = 4.7 × 10-8), involved in hematopoiesis. The most significant association was a missense variant in SAMHD1 (p = 4.2 × 10-28), found on a rare, 1.2-Mb shared ancestral haplotype on chromosome 20. SAMHD1 encodes a cytoplasmic host restriction factor involved in viral defense response and the mitochondrial nucleotide salvage pathway, and is associated with Aicardi-Goutières syndrome 5, a childhood encephalopathy and chronic inflammatory response disorder. Rare variants were enriched in Mendelian mtDNA depletion syndrome loci, and these variants implicated core processes in mtDNA replication, nucleoid structure formation, and maintenance. These data indicate that strong-effect mutations from the nuclear genome contribute to the genetic architecture of mtDNA-CN.

Respiratory chain dysfunction in perifascicular muscle fibres in patients with dermatomyositis is associated with mitochondrial DNA depletion

AIMS

Patients with dermatomyositis (DM) suffer from reduced aerobic metabolism contributing to impaired muscle function, which has been linked to cytochrome c oxidase (COX) deficiency in muscle tissue. This mitochondrial respiratory chain dysfunction is typically seen in perifascicular regions, which also show the most intense inflammatory reaction along with capillary loss and muscle fibre atrophy. The objective of this study was to investigate the pathobiology of the oxidative phosphorylation deficiency in DM.

METHODS

Muscle biopsy specimens with perifascicular COX deficiency from five juveniles and seven adults with DM were investigated. We combined immunohistochemical analyses of subunits in the respiratory chain including complex I (subunit NDUFB8), complex II (succinate dehydrogenase, subunit SDHB) and complex IV (COX, subunit MTCO1) with in situ hybridisation, next generation deep sequencing and quantitative polymerase chain reaction (PCR).

RESULTS

There was a profound deficiency of complexes I and IV in the perifascicular regions with enzyme histochemical COX deficiency, whereas succinate dehydrogenase activity and complex II were preserved. In situ hybridisation of mitochondrial RNA showed depletion of mitochondrial DNA (mtDNA) transcripts in the perifascicular regions. Analysis of mtDNA by next generation deep sequencing and quantitative PCR in affected muscle regions showed an overall reduction of mtDNA copy number particularly in the perifascicular regions.

CONCLUSION

The respiratory chain dysfunction in DM muscle is associated with mtDNA depletion causing deficiency of complexes I and IV, which are partially encoded by mtDNA, whereas complex II, which is entirely encoded by nuclear DNA, is preserved. The depletion of mtDNA indicates a perturbed replication of mtDNA explaining the muscle pathology and the disturbed aerobic metabolism.

Two type I topoisomerases maintain DNA topology in human mitochondria

Genetic processes require the activity of multiple topoisomerases, essential enzymes that remove topological tension and intermolecular linkages in DNA. We have investigated the subcellular localisation and activity of the six human topoisomerases with a view to understanding the topological maintenance of human mitochondrial DNA. Our results indicate that mitochondria contain two topoisomerases, TOP1MT and TOP3A. Using molecular, genomic and biochemical methods we find that both proteins contribute to mtDNA replication, in addition to the decatenation role of TOP3A, and that TOP1MT is stimulated by mtSSB. Loss of TOP3A or TOP1MT also dysregulates mitochondrial gene expression, and both proteins promote transcription elongation in vitro. We find no evidence for TOP2 localisation to mitochondria, and TOP2B knockout does not affect mtDNA maintenance or expression. Our results suggest a division of labour between TOP3A and TOP1MT in mtDNA topology control that is required for the proper maintenance and expression of human mtDNA.

Telomere length and mitochondrial DNA copy number in bipolar disorder: A sibling study

OBJECTIVES

An accelerated cellular ageing has been observed in bipolar disorder (BD) using biomarkers such as telomere length (TL) and mitochondrial DNA copy number (mtDNAcn). Several risk factors might drive premature ageing in individuals with BD, including a familial predisposition. This study compared TL and mtDNAcn between individuals with BD and their (un)-affected siblings, and explored factors that may explain proband-sibling differences.

METHODS

Sixty individuals with BD and seventy-four siblings (34 affected with BD or mood disorders and 40 unaffected) were included. Quantitative polymerase chain reaction (qPCR) was used to measure TL and mtDNAcn from peripheral blood genomic DNA.

RESULTS

TL and mtDNAcn did not significantly differ between probands and their siblings, whatever these latter were affected or not with mood disorders. However, the correlation plots of TL or mtDNAcn in proband-sibling pairs suggested that some pairs were discordant. The within proband-sibling pairs differences for TL and mtDNAcn were not explained by differences in all tested factors.

CONCLUSIONS

This study shows that probands with BD and their siblings are concordant for TL and mtDNAcn suggesting that they may share some environmental or genetic determinants of these two biomarkers of cellular ageing.

Mitochondrial DNA copy number, metabolic syndrome, and insulin sensitivity: Insights from the Sugar, Hypertension, and Physical Exercise studies

Mitochondrial DNA copy number (mtDNA-CN) measured in blood has been associated with many aging-related diseases, with higher mtDNA-CN typically associated with lower disease risk. Exercise training is an excellent preventative tool against aging-related disorders and has been shown to increase mitochondrial function in muscle. Using the Sugar, Hypertension, and Physical Exercise cohorts (N = 105), we evaluated the effect of 6-months of exercise intervention on mtDNA-CN measured in blood. Although there was no significant relationship between exercise intervention and mtDNA-CN change (P = 0.29), there was a nominally significant association between mtDNA-CN and metabolic syndrome (P = 0.04), which has been seen in previous literature. We also identified a nominally significant association between higher mtDNA-CN and higher insulin sensitivity (P = 0.02).

A bioinformatics pipeline for estimating mitochondrial DNA copy number and heteroplasmy levels from whole genome sequencing data

Mitochondrial diseases are a heterogeneous group of disorders that can be caused by mutations in the nuclear or mitochondrial genome. Mitochondrial DNA (mtDNA) variants may exist in a state of heteroplasmy, where a percentage of DNA molecules harbor a variant, or homoplasmy, where all DNA molecules have the same variant. The relative quantity of mtDNA in a cell, or copy number (mtDNA-CN), is associated with mitochondrial function, human disease, and mortality. To facilitate accurate identification of heteroplasmy and quantify mtDNA-CN, we built a bioinformatics pipeline that takes whole genome sequencing data and outputs mitochondrial variants, and mtDNA-CN. We incorporate variant annotations to facilitate determination of variant significance. Our pipeline yields uniform coverage by remapping to a circularized chrM and by recovering reads falsely mapped to nuclear-encoded mitochondrial sequences. Notably, we construct a consensus chrM sequence for each sample and recall heteroplasmy against the sample's unique mitochondrial genome. We observe an approximately 3-fold increased association with age for heteroplasmic variants in non-homopolymer regions and, are better able to capture genetic variation in the D-loop of chrM compared to existing software. Our bioinformatics pipeline more accurately captures features of mitochondrial genetics than existing pipelines that are important in understanding how mitochondrial dysfunction contributes to disease.

Neither cardiac mitochondrial DNA variation nor copy number contribute to congenital heart disease risk

The well-established manifestation of mitochondrial mutations in functional cardiac disease (e.g., mitochondrial cardiomyopathy) prompted the hypothesis that mitochondrial DNA (mtDNA) sequence and/or copy number (mtDNAcn) variation contribute to cardiac defects in congenital heart disease (CHD). MtDNAcns were calculated and rare, non-synonymous mtDNA mutations were identified in 1,837 CHD-affected proband-parent trios, 116 CHD-affected singletons, and 114 paired cardiovascular tissue/blood samples. The variant allele fraction (VAF) of heteroplasmic variants in mitochondrial RNA from 257 CHD cardiovascular tissue samples was also calculated. On average, mtDNA from blood had 0.14 rare variants and 52.9 mtDNA copies per nuclear genome per proband. No variation with parental age at proband birth or CHD-affected proband age was seen. mtDNAcns in valve/vessel tissue (320 ± 70) were lower than in atrial tissue (1,080 ± 320, p = 6.8E-21), which were lower than in ventricle tissue (1,340 ± 280, p = 1.4E-4). The frequency of rare variants in CHD-affected individual DNA was indistinguishable from the frequency in an unaffected cohort, and proband mtDNAcns did not vary from those of CHD cohort parents. In both the CHD and the comparison cohorts, mtDNAcns were significantly correlated between mother-child, father-child, and mother-father. mtDNAcns among people with European (mean = 52.0), African (53.0), and Asian haplogroups (53.5) were calculated and were significantly different for European and Asian haplogroups (p = 2.6E-3). Variant heteroplasmic fraction (HF) in blood correlated well with paired cardiovascular tissue HF (r = 0.975) and RNA VAF (r = 0.953), which suggests blood HF is a reasonable proxy for HF in heart tissue. We conclude that mtDNA mutations and mtDNAcns are unlikely to contribute significantly to CHD risk.

Benchmarking the Effectiveness and Accuracy of Multiple Mitochondrial DNA Variant Callers: Practical Implications for Clinical Application

Mitochondrial DNA (mtDNA) mutations contribute to human disease across a range of severity, from rare, highly penetrant mutations causal for monogenic disorders to mutations with milder contributions to phenotypes. mtDNA variation can exist in all copies of mtDNA or in a percentage of mtDNA copies and can be detected with levels as low as 1%. The large number of copies of mtDNA and the possibility of multiple alternative alleles at the same DNA nucleotide position make the task of identifying allelic variation in mtDNA very challenging. In recent years, specialized variant calling algorithms have been developed that are tailored to identify mtDNA variation from whole-genome sequencing (WGS) data. However, very few studies have systematically evaluated and compared these methods for the detection of both homoplasmy and heteroplasmy. A publicly available synthetic gold standard dataset was used to assess four mtDNA variant callers (Mutserve, mitoCaller, MitoSeek, and MToolBox), and the commonly used Genome Analysis Toolkit “best practices” pipeline, which is included in most current WGS pipelines. We also used WGS data from 126 trios and calculated the percentage of maternally inherited variants as a metric of calling accuracy, especially for homoplasmic variants. We additionally compared multiple pathogenicity prediction resources for mtDNA variants. Although the accuracy of homoplasmic variant detection was high for the majority of the callers with high concordance across callers, we found a very low concordance rate between mtDNA variant callers for heteroplasmic variants ranging from 2.8% to 3.6%, for heteroplasmy thresholds of 5% and 1%. Overall, Mutserve showed the best performance using the synthetic benchmark dataset. The analysis of mtDNA pathogenicity resources also showed low concordance in prediction results. We have shown that while homoplasmic variant calling is consistent between callers, there remains a significant discrepancy in heteroplasmic variant calling. We found that resources like population frequency databases and pathogenicity predictors are now available for variant annotation but still need refinement and improvement. With its peculiarities, the mitochondria require special considerations, and we advocate that caution needs to be taken when analyzing mtDNA data from WGS data.

Ultra-deep whole genome bisulfite sequencing reveals a single methylation hotspot in human brain mitochondrial DNA

While DNA methylation is established as a major regulator of gene expression in the nucleus, the existence of mitochondrial DNA (mtDNA) methylation remains controversial. Here, we characterized the mtDNA methylation landscape in the prefrontal cortex of neurological healthy individuals (n=26) and patients with Parkinson’s disease (n=27), using a combination of whole-genome bisulphite sequencing (WGBS) and bisulphite-independent methods. Accurate mtDNA mapping from WGBS data required alignment to an mtDNA reference only, to avoid misalignment to nuclear mitochondrial pseudogenes. Once correctly aligned, WGBS data provided ultra-deep mtDNA coverage (16,723 ± 7,711) and revealed overall very low levels of cytosine methylation. The highest methylation levels (5.49 ± 0.97%) were found on CpG position m.545, located in the heavy-strand promoter 1 region. The m.545 methylation was validated using a combination of methylation-sensitive DNA digestion and quantitative PCR analysis. We detected no association between mtDNA methylation profile and Parkinson’s disease. Interestingly, m.545 methylation correlated with the levels of mtDNA transcripts, suggesting a putative role in regulating mtDNA gene expression. In addition, we propose a robust framework for methylation analysis of mtDNA from WGBS data, which is less prone to false-positive findings due to misalignment of nuclear mitochondrial pseudogene sequences.

Circulating cell-free mitochondrial DNA in pregnancy

Circulating cell-free mitochondrial DNA (ccf-mtDNA) released upon cell injury or death stimulates diverse pattern recognition receptors to activate innate immune responses and initiate systemic inflammation. In this review, we discuss the temporal changes of ccf-mtDNA during pregnancy and its potential contribution to adverse pregnancy outcomes in pregnancy complications.

Identifying Somatic Mitochondrial DNA Mutations

Mitochondria are cellular organelles that play an essential role in eukaryotes, producing the energy needed for a cell to survive. Beyond the ~3.2 Gb of nuclear genomic DNA, each human cell has hundreds of mitochondria which carry one or a few copies of the 16.5 kb circular mitochondria DNA (mtDNA). Despite its small size, the circular genome encodes 37 genes, including 13 proteins that generate respiratory chain complexes together with other proteins of nuclear origin. Similar to nuclear genome, mtDNA in cancer cells frequently harbor somatically acquired alterations. Whole-genome or whole-exome sequencing of the tumor and its matched normal tissues (frequently blood or adjacent non-tumor tissues) enables sensitive and efficient detection of somatic mtDNA mutations. Because each cancer cell commonly carries hundreds to thousands of mtDNA copies, detection of mtDNA mutations is dependent on the heteroplasmic level of each mutation. Here, we describe strategies to accurately identify somatic mtDNA mutations in cancer genome studies.

Personality traits are consistently associated with blood mitochondrial DNA copy number estimated from genome sequences in two genetic cohort studies

Background

Mitochondrial DNA copy number (mtDNAcn) in tissues and blood can be altered in conditions like diabetes and major depression and may play a role in aging and longevity. However, little is known about the association between mtDNAcn and personality traits linked to emotional states, metabolic health, and longevity. This study tests the hypothesis that blood mtDNAcn is related to personality traits and mediates the association between personality and mortality.

Methods

We assessed the big five personality domains and facets using the Revised NEO Personality Inventory (NEO-PI-R), assessed depressive symptoms with the Center for Epidemiologic Studies Depression Scale (CES-D), estimated mtDNAcn levels from whole-genome sequencing, and tracked mortality in participants from the Baltimore Longitudinal Study of Aging. Results were replicated in the SardiNIA Project.

Results

We found that mtDNAcn was negatively associated with the Neuroticism domain and its facets and positively associated with facets from the other four domains. The direction and size of the effects were replicated in the SardiNIA cohort and were robust to adjustment for potential confounders in both samples. Consistent with the Neuroticism finding, higher depressive symptoms were associated with lower mtDNAcn. Finally, mtDNAcn mediated the association between personality and mortality risk.

Conclusions

To our knowledge, this is the first study to show a replicable association between mtDNAcn and personality. Furthermore, the results support our hypothesis that mtDNAcn is a biomarker of the biological process that explains part of the association between personality and mortality.

Funding

Support for this work was provided by the Intramural Research Program of the National Institute on Aging (Z01-AG000693, Z01-AG000970, and Z01-AG000949) and the National Institute of Neurological Disorders and Stroke of the National Institutes of Health. AT was also supported by the National Institute on Aging of the National Institutes of Health Grant R01AG068093.

Diet Quality Scores Are Positively Associated with Whole Blood-Derived Mitochondrial DNA Copy Number in the Framingham Heart Study

BACKGROUND

The association between diet quality and mitochondrial DNA copy number (mtDNA-CN) remains to be examined.

OBJECTIVES

We aimed to study the relation between diet quality and mtDNA-CN.

METHODS

We analyzed data from 2931 Framingham Heart Study (FHS) participants (mean age of 57 y, 55% females). Whole-genome sequencing was used to calculate mtDNA-CN from whole-blood samples. We examined the cross-sectional associations between 3 diet quality scores, the Dietary Approaches to Stop Hypertension (DASH) score, the Alternative Healthy Eating Index (AHEI), and the Mediterranean diet score (MDS), and mtDNA-CN. Linear mixed models were used to account for maternal lineage.

RESULTS

We observed that a higher DASH score was positively associated with mtDNA-CN after adjusting for sex, age, energy intake, smoking status, alcohol intake, and physical activity level. A 1-SD increase in the DASH score was associated with a 0.042-SD greater mtDNA-CN (95% CI: 0.007, 0.077; P = 0.02). Similarly, for each SD increase in AHEI and MDS, the mtDNA-CN SD increased by 0.056 (95% CI: 0.019, 0.092; P = 0.003) and 0.047 (95% CI: 0.01, 0.083; P = 0.01), respectively. Diet quality scores were associated with neutrophil and lymphocyte counts but not platelet counts, e.g., for a 1-SD increase in the DASH, neutrophils decreased by 0.8% (95% CI: 0.5%, 1.1%; P = 4.1 × 10-6), lymphocytes increased by 0.7% (95% CI: 0.4%, 1%, P = 1.2 × 10-5), and there was no significant change in platelet number (0.1 × 1000/μL; 95% CI: -1.6, 1.9; P = 0.89). Further adjustment for neutrophil, lymphocyte, and platelet counts and the associations between diet quality scores and mtDNA-CN were completely attenuated to nonsignificant (P = 0.95, 0.54, and 0.91, respectively).

CONCLUSIONS

We observed that higher diet quality is associated with a greater whole-blood derived mtDNA-CN in middle-aged to older adult FHS participants, and that blood cell composition, particularly neutrophil counts, attenuated the association between diet quality and mtDNA-CN.

Characterization of mitochondrial DNA quantity and quality in the human aged and Alzheimer's disease brain

BACKGROUND

Mitochondrial dysfunction is a feature of neurodegenerative diseases, including Alzheimer's disease (AD). Changes in the mitochondrial DNA copy number (mtDNAcn) and increased mitochondrial DNA mutation burden have both been associated with neurodegenerative diseases and cognitive decline. This study aims to systematically identify which common brain pathologies in the aged human brain are associated with mitochondrial recalibrations and to disentangle the relationship between these pathologies, mtDNAcn, mtDNA heteroplasmy, aging, neuronal loss, and cognitive function.

METHODS

Whole-genome sequencing data from n = 1361 human brain samples from 5 different regions were used to quantify mtDNAcn as well as heteroplasmic mtDNA point mutations and small indels. Brain samples were assessed for 10 common pathologies. Annual cognitive test results were used to assess cognitive function proximal to death. For a subset of samples, neuronal proportions were estimated from RNA-seq profiles, and mass spectrometry was used to quantify the mitochondrial protein content of the tissue.

RESULTS

mtDNAcn was 7-14% lower in AD relative to control participants. When accounting for all 10 common neuropathologies, only tau was significantly associated with lower mtDNAcn in the dorsolateral prefrontal cortex. In the posterior cingulate cortex, TDP-43 pathology demonstrated a distinct association with mtDNAcn. No changes were observed in the cerebellum, which is affected late by pathologies. Neither age nor gender was associated with mtDNAcn in the studied brain regions when adjusting for pathologies. Mitochondrial content and mtDNAcn independently explained variance in cognitive function unaccounted by pathologies, implicating complex mitochondrial recalibrations in cognitive decline. In contrast, mtDNA heteroplasmy levels increased by 1.5% per year of life in the cortical regions, but displayed no association with any of the pathologies or cognitive function.

CONCLUSIONS

We studied mtDNA quantity and quality in relation to mixed pathologies of aging and showed that tau and not amyloid-β is primarily associated with reduced mtDNAcn. In the posterior cingulate cortex, the association of TDP-43 with low mtDNAcn points to a vulnerability of this region in limbic-predominant age-related TDP-43 encephalopathy. While we found low mtDNAcn in brain regions affected by pathologies, the absence of associations with mtDNA heteroplasmy burden indicates that mtDNA point mutations and small indels are unlikely to be involved in the pathogenesis of late-onset neurodegenerative diseases.

The Association of Mitochondrial Copy Number With Sarcopenia in Adult Survivors of Childhood Cancer

BACKGROUND

Adult childhood cancer survivors are at risk for frailty, including low muscle mass and weakness (sarcopenia). Using peripheral blood mitochondrial DNA copy number (mtDNAcn) as a proxy for functional mitochondria, this study describes cross-sectional associations between mtDNAcn and sarcopenia among survivors.

METHODS

Among 1762 adult childhood cancer survivors (51.6% male; median age = 29.4 years, interquartile range [IQR] = 23.3-36.8), with a median of 20.6 years from diagnosis (IQR = 15.2-28.2), mtDNAcn estimates were derived from whole-genome sequencing. A subset was validated by quantitative polymerase chain reaction and evaluated cross-sectionally using multivariable logistic regression for their association with sarcopenia, defined by race-, age-, and sex-specific low lean muscle mass or weak grip strength. All statistical tests were 2-sided.

RESULTS

The prevalence of sarcopenia was 27.0%, higher among female than male survivors (31.5% vs 22.9%; P < .001) and associated with age at diagnosis; 51.7% of survivors with sarcopenia were diagnosed ages 4-13 years (P = .01). Sarcopenia was most prevalent (39.0%) among central nervous system tumor survivors. Cranial radiation (odds ratio [OR] = 1.84, 95% confidence interval [CI] = 1.32 to 2.59) and alkylating agents (OR = 1.34, 95% CI = 1.04 to 1.72) increased, whereas glucocorticoids decreased odds (OR = 0.72, 95% CI = 0.56 to 0.93) of sarcopenia. mtDNAcn decreased with age (β = -0.81, P = .002) and was higher among female survivors (β = 9.23, P = .01) and among survivors with a C allele at mt.204 (β = -17.9, P = .02). In adjusted models, every standard deviation decrease in mtDNAcn increased the odds of sarcopenia 20% (OR = 1.20, 95% CI = 1.07 to 1.34).

CONCLUSIONS

A growing body of evidence supports peripheral blood mtDNAcn as a biomarker for adverse health outcomes; however, this study is the first to report an association between mtDNAcn and sarcopenia among childhood cancer survivors.

Mitochondrial DNA copy number and heteroplasmy load correlate with skeletal muscle oxidative capacity by P31 MR spectroscopy

The association between blood‐based estimates of mitochondrial DNA parameters, mitochondrial DNA copy number (mtDNA‐CN) and heteroplasmy load, with skeletal muscle bioenergetic capacity was evaluated in 230 participants of the Baltimore Longitudinal Study of Aging (mean age:74.7 years, 53% women). Participants in the study sample had concurrent data on muscle oxidative capacity (τ_PCr) assessed by ³¹P magnetic resonance spectroscopy, and mitochondrial DNA parameters estimated from whole‐genome sequencing data. In multivariable linear regression models, adjusted for age, sex, extent of phosphocreatine (PCr) depletion, autosomal sequencing coverage, white blood cell total, and differential count, as well as platelet count, mtDNA‐CN and heteroplasmy load were not significantly associated with τ_PCr (both p > 0.05). However, in models evaluating whether the association between mtDNA‐CN and τ_PCr varied by heteroplasmy load, there was a significant interaction between mtDNA‐CN and heteroplasmy load (p = 0.037). In stratified analysis, higher mtDNA‐CN was significantly associated with lower τ_PCr among participants with high heteroplasmy load (n = 84, β (SE) = −0.236 (0.115), p‐value = 0.044), but not in those with low heteroplasmy load (n = 146, β (SE) = 0.046 (0.119), p‐value = 0.702). Taken together, mtDNA‐CN and heteroplasmy load provide information on muscle bioenergetics. Thus, mitochondrial DNA parameters may be considered proxy measures of mitochondrial function that can be used in large epidemiological studies, especially when comparing subgroups.

Association of mitochondrial DNA copy number with cardiometabolic diseases

Mitochondrial DNA (mtDNA) is present in multiple copies in human cells. We evaluated cross-sectional associations of whole blood mtDNA copy number (CN) with several cardiometabolic disease traits in 408,361 participants of multiple ancestries in TOPMed and UK Biobank. Age showed a threshold association with mtDNA CN: among younger participants (<65 years of age), each additional 10 years of age was associated with 0.03 standard deviation (s.d.) higher level of mtDNA CN (P = 0.0014) versus a 0.14 s.d. lower level of mtDNA CN (P = 1.82 × 10^-13) among older participants (≥65 years). At lower mtDNA CN levels, we found age-independent associations with increased odds of obesity (P = 5.6 × 10^-238), hypertension (P = 2.8 × 10^-50), diabetes (P = 3.6 × 10^-7), and hyperlipidemia (P = 6.3 × 10^-5). The observed decline in mtDNA CN after 65 years of age may be a key to understanding age-related diseases.

Emerging methods for and novel insights gained by absolute quantification of mitochondrial DNA copy number and its clinical applications

The past thirty years have seen a surge in interest in pathophysiological roles of mitochondria, and the accurate quantification of mitochondrial DNA copy number (mCN) in cells and tissue samples is a fundamental aspect of assessing changes in mitochondrial health and biogenesis. Quantification of mCN between studies is surprisingly variable due to a combination of physiological variability and diverse protocols being used to measure this endpoint. The advent of novel methods to quantify nucleic acids like digital polymerase chain reaction (dPCR) and high throughput sequencing offer the ability to measure absolute values of mCN. We conducted an in-depth survey of articles published between 1969 -- 2020 to create an overview of mCN values, to assess consensus values of tissue-specific mCN, and to evaluate consistency between methods of assessing mCN. We identify best practices for methods used to assess mCN, and we address the impact of using specific loci on the mitochondrial genome to determine mCN. Current data suggest that clinical measurement of mCN can provide diagnostic and prognostic value in a range of diseases and health conditions, with emphasis on cancer and cardiovascular disease, and the advent of means to measure absolute mCN should improve future clinical applications of mCN measurements.

Mitochondrial genetic variation is enriched in G-quadruplex regions that stall DNA synthesis in vitro

As the powerhouses of the eukaryotic cell, mitochondria must maintain their genomes which encode proteins essential for energy production. Mitochondria are characterized by guanine-rich DNA sequences that spontaneously form unusual three-dimensional structures known as G-quadruplexes (G4). G4 structures can be problematic for the essential processes of DNA replication and transcription because they deter normal progression of the enzymatic-driven processes. In this study, we addressed the hypothesis that mitochondrial G4 is a source of mutagenesis leading to base-pair substitutions. Our computational analysis of 2757 individual genomes from two Italian population cohorts (SardiNIA and InCHIANTI) revealed a statistically significant enrichment of mitochondrial mutations within sequences corresponding to stable G4 DNA structures. Guided by the computational analysis results, we designed biochemical reconstitution experiments and demonstrated that DNA synthesis by two known mitochondrial DNA polymerases (Pol γ, PrimPol) in vitro was strongly blocked by representative stable G4 mitochondrial DNA structures, which could be overcome in a specific manner by the ATP-dependent G4-resolving helicase Pif1. However, error-prone DNA synthesis by PrimPol using the G4 template sequence persisted even in the presence of Pif1. Altogether, our results suggest that genetic variation is enriched in G-quadruplex regions that impede mitochondrial DNA replication.

The memory of neuronal mitochondrial stress is inherited transgenerationally via elevated mitochondrial DNA levels

The memory of stresses experienced by parents can be passed on to descendants as a forecast of the challenges to come. Here, we discovered that the neuronal mitochondrial perturbation-induced systemic mitochondrial unfolded protein response (UPR^mt) in Caenorhabditis elegans can be transmitted to offspring over multiple generations. The transgenerational activation of UPR^mt is mediated by maternal inheritance of elevated levels of mitochondrial DNA (mtDNA), which causes the proteostasis stress within mitochondria. Furthermore, results from intercrossing studies using wild C. elegans strains further support that maternal inheritance of higher levels of mtDNA can induce the UPR^mt in descendants. The mitokine Wnt signalling pathway is required for the transmission of elevated mtDNA levels across generations, thereby conferring lifespan extension and stress resistance to offspring. Collectively, our results reveal that the nervous system can transmit stress signals across generations by increasing mtDNA in the germline, enabling descendants to better cope with anticipated challenges.

Association of Mitochondrial DNA Copy Number and Telomere Length with Prevalent and Incident Cancer and Cancer Mortality in Women: A Prospective Swedish Population-Based Study

Changes in mitochondrial DNA copy number (mtDNA-CN) and telomere length have, separately, been proposed as risk factors for various cancer types. However, those results are conflicting. Here, mtDNA-CN and relative telomere length were measured in 3225 middle-aged women included in a large population-based prospective cohort. The baseline mtDNA-CN in patients with prevalent breast cancer was significantly higher (12.39 copies/µL) than cancer-free individuals. During an average of 15.2 years of follow-up, 520 patients were diagnosed with cancer. Lower mtDNA-CN was associated with decreased risk of genital organ cancer (hazard ratio (HR), 0.84), and shorter telomere length was associated with increased risk of urinary system cancer (HR, 1.79). Furthermore, mtDNA-CN was inversely associated with all-cause (HR, 1.20) and cancer-specific mortality (HR, 1.21) when considering all cancer types. Surprisingly, shorter telomere length was associated with decreased risk of cancer-specific mortality when considering all cancer types (HR, 0.85). Finally, lower mtDNA-CN and shorter telomere length were associated with increased risk of both all-cause and cancer-specific mortality in genital organ cancer patients. In this study population, we found that mtDNA-CN and telomere length were significantly associated with prevalent and incident cancer and cancer mortality. However, these associations were cancer type specific and need further investigation.

Clinical Bioinformatics in Precise Diagnosis of Mitochondrial Disease

Clinical bioinformatics system is well-established for diagnosing genetic disease based on next-generation sequencing, but requires special considerations when being adapted for the next-generation sequencing-based genetic diagnosis of mitochondrial diseases. Challenges are caused by the involvement of mitochondrial DNA genome in disease etiology. Heteroplasmy and haplogroup are key factors in interpreting mitochondrial DNA variant effects. Data resources and tools for analyzing variant and sequencing data are available at MSeqDR, MitoMap, and HmtDB. Revised specifications of the American College of Medical Genetics/Association of Molecular Pathology standards and guidelines for mitochondrial DNA variant interpretation are proposed by the MSeqDr Consortium and community experts.

Mitochondrial function in development and disease

Mitochondria are organelles with vital functions in almost all eukaryotic cells. Often described as the cellular ‘powerhouses’ due to their essential role in aerobic oxidative phosphorylation, mitochondria perform many other essential functions beyond energy production. As signaling organelles, mitochondria communicate with the nucleus and other organelles to help maintain cellular homeostasis, allow cellular adaptation to diverse stresses, and help steer cell fate decisions during development. Mitochondria have taken center stage in the research of normal and pathological processes, including normal tissue homeostasis and metabolism, neurodegeneration, immunity and infectious diseases. The central role that mitochondria assume within cells is evidenced by the broad impact of mitochondrial diseases, caused by defects in either mitochondrial or nuclear genes encoding for mitochondrial proteins, on different organ systems. In this Review, we will provide the reader with a foundation of the mitochondrial ‘hardware’, the mitochondrion itself, with its specific dynamics, quality control mechanisms and cross-organelle communication, including its roles as a driver of an innate immune response, all with a focus on development, disease and aging. We will further discuss how mitochondrial DNA is inherited, how its mutation affects cell and organismal fitness, and current therapeutic approaches for mitochondrial diseases in both model organisms and humans.

Summary: Mitochondria have a plethora of functions beyond metabolism. This Review discusses the emerging and multifaceted roles of mitochondria in different model organisms and human disease biology.

Mitochondrial genome copy number measured by DNA sequencing in human blood is strongly associated with metabolic traits via cell-type composition differences

BACKGROUND

Mitochondrial genome copy number (MT-CN) varies among humans and across tissues and is highly heritable, but its causes and consequences are not well understood. When measured by bulk DNA sequencing in blood, MT-CN may reflect a combination of the number of mitochondria per cell and cell-type composition. Here, we studied MT-CN variation in blood-derived DNA from 19184 Finnish individuals using a combination of genome (N = 4163) and exome sequencing (N = 19034) data as well as imputed genotypes (N = 17718).

RESULTS

We identified two loci significantly associated with MT-CN variation: a common variant at the MYB-HBS1L locus (P = 1.6 × 10-8), which has previously been associated with numerous hematological parameters; and a burden of rare variants in the TMBIM1 gene (P = 3.0 × 10-8), which has been reported to protect against non-alcoholic fatty liver disease. We also found that MT-CN is strongly associated with insulin levels (P = 2.0 × 10-21) and other metabolic syndrome (metS)-related traits. Using a Mendelian randomization framework, we show evidence that MT-CN measured in blood is causally related to insulin levels. We then applied an MT-CN polygenic risk score (PRS) derived from Finnish data to the UK Biobank, where the association between the PRS and metS traits was replicated. Adjusting for cell counts largely eliminated these signals, suggesting that MT-CN affects metS via cell-type composition.

CONCLUSION

These results suggest that measurements of MT-CN in blood-derived DNA partially reflect differences in cell-type composition and that these differences are causally linked to insulin and related traits.

Mitochondrial DNA variants in inclusion body myositis characterized by deep sequencing

Muscle pathology in inclusion body myositis (IBM) typically includes inflammatory cell infiltration, muscle fibers with rimmed vacuoles and cytochrome c oxidase (COX)‐deficient fibers. Previous studies have revealed clonal expansion of large mitochondrial DNA (mtDNA) deletions in the COX‐deficient muscle fibers. Technical limitations have prevented complete investigations of the mtDNA deletions and other mtDNA variants. Detailed characterization by deep sequencing of mtDNA in muscle samples from 21 IBM patients and 10 age‐matched controls was performed after whole genome sequencing with a mean depth of mtDNA coverage of 46,000x. Multiple large mtDNA deletions and duplications were identified in all IBM and control muscle samples. In general, the IBM muscles demonstrated a larger number of deletions and duplications with a mean heteroplasmy level of 10% (range 1%‐35%) compared to controls (1%, range 0.2%‐3%). There was also a small increase in the number of somatic single nucleotide variants in IBM muscle. More than 200 rearrangements were recurrent in at least two or more IBM muscles while 26 were found in both IBM and control muscles. The deletions and duplications, with a high recurrence rate, were mainly observed in three mtDNA regions, m.534‐4429, m.6330‐13993, and m.8636‐16072, where some were flanked by repetitive sequences. The mtDNA copy number in IBM muscle was reduced to 42% of controls. Immunohistochemical and western blot analyses of IBM muscle revealed combined complex I and complex IV deficiency affecting the COX‐deficient fibers. In conclusion, deep sequencing and quantitation of mtDNA variants revealed that IBM muscles had markedly increased levels of large deletions and duplications, and there were also indications of increased somatic single nucleotide variants and reduced mtDNA copy numbers compared to age‐matched controls. The distribution and type of variants were similar in IBM muscle and controls indicating an accelerated aging process in IBM muscle, possibly associated with chronic inflammation.

Comparison of whole genome sequencing and targeted sequencing for mitochondrial DNA

Mitochondrial dysfunction has emerged to be associated with a broad spectrum of diseases, and there is an increasing demand for accurate detection of mitochondrial DNA (mtDNA) variants. Whole genome sequencing (WGS) has been the dominant sequencing approach to identify genetic variants in recent decades, but most studies focus on variants on the nuclear genome. Whole genome sequencing is also costly and time consuming. Sequencing specifically targeted for mtDNA is commonly used in the diagnostic settings and has lower costs. However, there is a lack of pairwise comparisons between these two sequencing approaches for calling mtDNA variants on a population basis. In this study, we compared WGS and mtDNA-targeted sequencing (targeted-seq) in analyzing mitochondrial DNA from 1499 participants recruited into the Severe Asthma Research Program (SARP). Our study reveals that targeted-sequencing and WGS have comparable capacity to determine genotypes and to call haplogroups and homoplasmies on mtDNA. However, there exists a large variability in calling heteroplasmies, especially for low-frequency heteroplasmies, which indicates that investigators should be cautious about heteroplasmies acquired from different sequencing methods. Further research is highly desired to improve variant detection methods for mitochondrial DNA.

Validation of whole genome sequencing from dried blood spots

BACKGROUND

Dried blood spots (DBS) are a relatively inexpensive source of nucleic acids and are easy to collect, transport, and store in large-scale field surveys, especially in resource-limited settings. However, their performance in whole-genome sequencing (WGS) relative to that of venous blood DNA has not been analyzed for various downstream applications.

METHODS

This study compares the WGS performance of DBS paired with venous blood samples collected from 12 subjects.

RESULTS

Results of standard quality checks of coverage, base quality, and mapping quality were found to be near identical between DBS and venous blood. Concordance for single-nucleotide variants, insertions and deletions, and copy number variants was high between these two sample types. Additionally, downstream analyses typical of population-based studies were performed, such as mitochondrial heteroplasmy detection, haplotype analysis, mitochondrial copy number changes, and determination of telomere lengths. The absolute mitochondrial copy number values were higher for DBS than for venous blood, though the trend in sample-to-sample variation was similar between DBS and blood. Telomere length estimates in most DBS samples were on par with those from venous blood.

CONCLUSION

DBS samples can serve as a robust and feasible alternative to venous blood for studies requiring WGS analysis.

Nuclear genome-wide associations with mitochondrial heteroplasmy

The role of the nuclear genome in maintaining the stability of the mitochondrial genome (mtDNA) is incompletely known. mtDNA sequence variants can exist in a state of heteroplasmy, which denotes the coexistence of organellar genomes with different sequences. Heteroplasmic variants that impair mitochondrial capacity cause disease, and the state of heteroplasmy itself is deleterious. However, mitochondrial heteroplasmy may provide an intermediate state in the emergence of novel mitochondrial haplogroups. We used genome-wide genotyping data from 982,072 European ancestry individuals to evaluate variation in mitochondrial heteroplasmy and to identify the regions of the nuclear genome that affect it. Age, sex, and mitochondrial haplogroup were associated with the extent of heteroplasmy. GWAS identified 20 loci for heteroplasmy that exceeded genome-wide significance. This included a region overlapping mitochondrial transcription factor A (TFAM), which has multiple roles in mtDNA packaging, replication, and transcription. These results show that mitochondrial heteroplasmy has a heritable nuclear component.