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 31 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. Published 2021. This article is a U.S. Government work and is in the public domain in the USA. Aging Cell published by Anatomical Society and John Wiley & Sons Ltd.

INTRODUCTION

Mitochondria produce energy for metabolic and functional activity through aerobic metabolism and are also linked to a broad range of cellular processes including apoptosis, and immune signaling, iron and calcium homeostasis, and reactive oxygen species signaling (Gonzalez‐Freire et al., 2015). Dysfunction of these activities has been implicated in the development of chronic disease and is also considered a hallmark of aging (Lopez‐Otin et al., 2013). Muscle is an energetically demanding tissue that is central to decline in physical function in later life. However, commonly used methods to assess mitochondrial function in skeletal muscle, such as respirometry in muscle biopsies and 31P magnetic resonance spectroscopy (MRS), are resource‐intensive and impractical for population‐based studies (Coen et al., 2013; Conley et al., 2000; Short et al., 2005). Recent studies have shown that measurements of mitochondrial oxidative capacity in human skeletal muscle via 31P‐MRS are associated with multiple chronic diseases and morbidity (AlGhatrif et al., 2017; Brown et al., 2019; Zampino et al., 2020). Techniques have also been developed to estimate mitochondrial DNA copy number (mtDNA‐CN) and heteroplasmy load from whole‐genome sequencing (WGS), most often in blood samples, such as buffy coat specimens (Ding et al., 2015; Qian et al., 2017). However, the relationship between mtDNA‐CN and heteroplasmy load estimated from WGS and mitochondrial function in skeletal muscle is unknown, as mitochondrial characteristics, number, and volume vary by organ, tissue, and cell types. If such a relationship was established, mtDNA‐CN and heteroplasmy load would be invaluable for clinical research, including studies of aging, as proxies for muscle bioenergetic status. In addition, they would form therapeutic targets for drugs and interventions designed to address chronic disease and improve function in aging through improved mitochondrial function (Andreux et al., 2019).

Using data collected from 230 participants of the Baltimore Longitudinal Study of Aging (BLSA), a study of community‐dwelling individuals, we tested the hypothesis that mtDNA‐CN and heteroplasmy load estimated from WGS would be associated with muscle bioenergetic status as assessed with the phosphocreatine (PCr) exponential recovery time constant (τPCr) determined by 31P‐MRS (Coen et al., 2013). The mean age of the study sample was 74.7 years, and 53% of participants were women (Table 1). Study participants were free of dismobility, mean usual gait speed of 1.20 m/s with a minimal value of 0.63 m/s (Table 1) (Cummings et al., 2014). In this sample, higher values of τPCr were associated with older age, indicating age‐related decline in mitochondrial function (Figure S1). BLSA protocols are approved by the National Institutes of Health Institutional Review Board, and all participants provided written informed consent.

TABLE 1

Participant characteristics

			Heteroplasmy load ≤3		Heteroplasmy load >3
	Overall (n = 230)		High mtDNA‐CN (n = 71)	Low mtDNA‐CN (n = 75)	High mtDNA‐CN (n = 44)	Low mtDNA‐CN (n = 40)
	Mean ± standard deviation or N (%)	range	Mean ± standard deviation or N (%)	Mean ± standard deviation or N (%)	Mean ± standard deviation or N (%)	Mean ± standard deviation or N (%)
Age, years	74.7 ± 9.8	50 – 93	71.3 ± 9.1	75.2 ± 9.9	75.5 ± 9.8	78.7 ± 9.6
Women	121 (53)	‐	53 (75)	27 (36)	27 (61)	14 (35)
Body mass index, kg/m2	26.2 ± 4.3	17.1–42.9	26.6 ± 4.1	26.9 ± 4.9	24.8 ± 3.8	25.9 ± 3.5
6‐meter usual gait speed, m/s	1.20 ± 0.20	0.63–1.81	1.23 ± 0.19	1.20 ± 0.19	1.21 ± 0.23	1.15 ± 0.21
Any cancer	97 (42)	‐	27 (39)	37 (49)	13 (30)	20 (53)
Non‐skin cancer	26 (11)	‐	4 (6)	9 (12)	7 (16)	6 (15)
Diabetes	11 (5)		2 (3)	5 (7)	2 (5)	2 (5)
Current or former smokers	6 (3)	‐	3 (4)	1 (1)	1 (2)	1 (3)
High physical activity category (≥150 min/week)	63 (27)	‐	21 (30)	18 (24)	14 (32)	10 (25)
τPCr	49.4 ± 10.4	22.6 – 3.8	49.4 ± 10.0	48.8 ± 10.1	46.2 ± 10.6	54.1 ± 10.2
PCr depletion, %	59.3 ± 10.2	33.1–91.0	59 ± 10	59 ± 11	59 ± 10	59 ± 10
mtDNA‐CN	219 ± 42	127–366	250 ± 29	188 ± 18	257 ± 34	181 ± 23
Heteroplasmy load	Median (IQR)	range	Median (IQR)	Median (IQR)	Median (IQR)	Median (IQR)
	3 (2)	1–96	2 (1)	2 (1)	5 (2)	5 (3)

mtDNA‐CN was estimated using the mitoCalc algorithm ((Ding et al., 2015; Qian et al., 2017) in WGS data from buffy coat samples. Heteroplasmic variants (mtDNA variants with more than one allele at a DNA site) were identified in the same sequencing data using the mitoCaller algorithm. Heteroplasmy load is represented by the total number of heteroplasmic variants in each individual. In vivo 31P‐MRS measurements of the concentrations of phosphorus‐containing metabolites including phosphocreatine (PCr) were obtained from the vastus lateralis muscle using 31P MRS at 3T, following a standardized protocol (Choi et al., 2016). τPCr, the PCr exponential recovery time constant measured in seconds, was calculated by fitting time‐dependent changes in PCr peak area to the monoexponential recovery function: where PCr(0) is the end‐of‐exercise PCr signal area and ΔPCr is the decrease in signal area from its pre‐exercise value (Choi et al., 2016). Higher values of τPCr indicate longer recovery and lower oxidative capacity. Associations of mtDNA‐CN and heteroplasmy load with τPCr were tested using multivariable linear regression (SAS v9.4; SAS Institute, Inc). After adjustment for age, sex, extent of 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 (Table S1, Model 1, Model 2).

Since recent data suggest that mtDNA‐CN and heteroplasmy load provide complementary information on mitochondrial function in patients with peripheral artery disease (Gonzalez‐Freire et al., 2020), we also tested the hypothesis that the association between mtDNA‐CN and τPCr would be different according to levels of heteroplasmy load through the evaluation of an interaction term in multivariable models. There was a significant interaction between mtDNA‐CN and heteroplasmy load (Table S1, Model 3); after stratifying by a median split of heteroplasmy load of 3, higher mtDNA‐CN was significantly associated with lower τPCr in participants with high heteroplasmy load (Table S2; Figure 1), while there was not a significant association between mtDNA‐CN and τPCr in those with low heteroplasmy load (Table S2; Figure 1). At high heteroplasmy load, there was also a relationship between mtDNA‐CN and/or τPCr and gait speed (p = 0.06 and <0.001, respectively). Results were not substantially changed after excluding participants with diabetes and non‐skin cancer after full adjustment (β (SE), p‐value for those with high (n = 68) and low (n = 126) heteroplasmy load: −0.315 (0.128), 0.017; and 0.117 (0.113), 0.302, respectively). Results remained similar in sex‐stratified analysis (Table S3).

FIGURE 1

Scatter plots (A) (red: women; blue: men) and boxplots (B) for the unadjusted relationship between mitochondrial DNA copy number and τPCr in participants with low and high heteroplasmy load. Caption: τPCr, the PCr exponential recovery time constant measured in seconds, higher values of which indicate longer recovery time and lower oxidative capacity. In scatter plots, the specific amount of mitochondrial DNA copy number attributed to low and high heteroplasmy load ranged from 136 to 356, and from 127 to 366, respectively. In boxplots, mitochondrial DNA copy number was categorized by a median split at 216

In this sample of community‐dwelling older adults, we demonstrate that mtDNA‐CN and heteroplasmy load provide complementary information on mitochondrial oxidative capacity measured in skeletal muscle. The significant interaction between mtDNA‐CN and heteroplasmy load indicated that the relationship of mtDNA‐CN with muscle bioenergetics was different according to levels of heteroplasmy load. In particular, we found that mtDNA‐CN was associated with skeletal muscle oxidative capacity only in individuals with high heteroplasmy load. While the mechanism for this interaction remains unclear, we have previously demonstrated that mtDNA‐CN can be associated with both positive and negative health outcomes, for example in participants with and without diabetes (Moore et al., 2018). Our results support the notion that mtDNA‐CN can be assessed at the population level and it can be correlated with other parameters and used to compare subgroups. It may also serve as a measure sensitive to mitochondrial mass in skeletal muscle. However, if higher heteroplasmy load leads to altered mitochondrial damage, increased mtDNA‐CN may indicate a homeostatic response toward increasing mitochondrial biogenesis (Filograna et al., 2020; Qian & Van Houten, 2010). It is also possible that increased mtDNA‐CN may be observed in conjunction with a high level of heteroplasmy load within the context of damage induced by oxidative stress, a strong correlate of impaired mitochondrial function. Thus, increased mtDNA‐CN may reflect the spilling of mtDNA from damaged mitochondria, and in this case, mitochondrial mass may not correlate with function. It can be hypothesized that reduced ATP production may be sensed by AMP kinase which in turn stimulates mitochondrial biogenesis through PGC‐1 alpha leading to increased mtDNA synthesis through transcription factor A, which is required for maintenance of normal levels of mtDNA (Kasashima et al., 2012). The additional mtDNA that is not incorporated in functioning mitochondria may enter the circulation. These data suggest that mtDNA‐CN and heteroplasmy load should be used in conjunction to obtain insight into mitochondrial function.

Our study has certain limitations, including a modest sample size; the BLSA population from which the study sample is drawn is healthier, less diverse, and more well‐educated than the general population. Future studies with larger samples encompassing a greater range of functional status would provide validation of our results. In addition, our study has several strengths. The study sample is comprised of well‐characterized community‐dwelling older men and women, allowing us to control for multiple known covariates: Our observations were robust to adjustment for differential white blood cell and platelet count as well as to exclusion of participants with diabetes and non‐skin cancer, both of which are known to affect mtDNA‐CN.

In sum, when taken together mtDNA parameters, mitochondrial DNA copy number and heteroplasmy load provide an indication of muscle bioenergetics. Assessing mitochondrial DNA copy number and heteroplasmy load may provide cost‐effective and accessible insight into muscle mitochondrial function in large epidemiological studies.

CONFLICT OF INTEREST

None declared.

AUTHOR CONTRIBUTIONS

QT developed the study concept and statistical analysis strategy, analyzed, interpreted, and drafted the manuscript. AZM developed the study concept and statistical analysis strategy, interpreted, and reviewed the manuscript. JD, RO, and MP involved in data collection and reviewed the manuscript. KWF and RGS involved in P31 MR Spectroscopy study protocol and reviewed the manuscript. LF designed the study and reviewed the manuscript. All authors edited and approved the manuscript.

Supplementary Material

Click here for additional data file.