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

Weiling Xu, Yun Soo Hong, Bo Hu, Suzy A. A. Comhair, Allison J. Janocha, Joe G. Zein, Ruoying Chen, Deborah A. Meyers, David T. Mauger, Victor E. Ortega, Eugene R. Bleecker, Mario Castro, Loren C. Denlinger, John V. Fahy, Elliot Israel, Bruce D. Levy, Nizar N. Jarjour, Wendy C. Moore, Sally E. Wenzel, Benjamin Gaston, Chunyu Liu, Dan E. Arking, Serpil C. Erzurum, the National Heart, Lung, and Blood Institute (NHLBI) Severe Asthma Research Program (SARP), TOPMed mtDNA Working Group in NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium
doi: https://doi.org/10.1101/2023.12.05.23299392

Abstract

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.

Introduction

One of the most common chronic diseases in the United States, asthma is characterized by airway inflammation and hyper-reactivity. Severe disease, which occurs in 10% of asthmatics, is difficult to treat, and those with frequent exacerbations, which occur in 5% of asthmatics, suffer a disabling burden with asthma. Thus, prevention of exacerbations is the major goal of global asthma care guidelines. Yet, we lack understanding of the molecular mechanisms that lead to the exacerbation-prone asthma phenotype and thus, are unable to apply specific prevention interventions. Pre-clinical models and human studies identify excess production of reactive oxygen species, with bursts of oxidants during exacerbations (1-6). Although mitochondria are pivotal to oxidant generation and inflammation, and changes in mitochondria function are linked to asthma (3, 7), the relationship of asthma risk, severity, or exacerbations to mitochondria is unknown.

Mitochondria produce energy from oxidative-phosphorylation and regulate redox status. Each human cell contains hundreds to thousands of mitochondria, dependent on tissue and cell type (8). Mitochondria have their own circular DNA (mtDNA) that contains 13 genes encoding oxidative phosphorylation proteins crucial for bioenergetics, but the majority of oxidative phosphorylation proteins are encoded by nuclear DNA. Mitochondria independently replicate, transcribe, and translate their mtDNA using their own ribosomal and transfer RNA, which are encoded by mtDNA. Thus, reductions in mtDNA-CN lead to low levels of mitochondrial-encoded proteins in oxidative phosphorylation, decreased respiratory capacity, and greater production of reactive oxygen species, all of which lead to immune dysfunction, inflammation, and altered cell signaling (9, 10). Thus, mitochondrial DNA copy number (mtDNA-CN) is a good proxy of mitochondrial function, with low levels indicating mitochondrial dysfunction. Variation of mtDNA copy number (mtDNA-CN) has been successfully used as an indirect assessment of mitochondrial function in population studies (9, 10), i.e., low mtDNA-CN has been associated with cardiovascular disease, diabetes, and cancers (9, 11).

We hypothesize that mitochondrial dysfunction is associated with asthma and the parameters of severity and exacerbation. To test if mitochondrial dysfunction is associated with asthma, we analyzed mtDNA-CN estimated from a combination of whole exome sequencing data and genotyping data in participants of the UK Biobank (UKB). To carefully determine whether lower mtDNA-CN is associated with asthma severity and whether there is a relationship to exacerbations, we analyzed mtDNA-CN, measures of inflammation and oxidative stress, and asthma exacerbations in severe and non-severe asthmatics of the large-scale consortium-based Severe Asthma Research Program (SARP) through NHLBI Trans-Omics for Precision Medicine (TOPMed).

Methods

Study Population

The curated UKB dataset includes 268,926 unrelated individuals who self-identify as White (European) ethnicity including 29,768 asthmatics and 239,158 non-asthmatics. Asthmatics were defined based on self-report (12). The UKB curated exclusion criteria are detailed in the supplement and in prior reports [see UKB criteria in supplement] (13). Exclusions specific to these analyses include current smokers, cancer, other lung diseases, gastro-esophageal reflux disease, chronic sinusitis, sleep apnea, and kidney, liver, and heart failure [supplement, Table E1]. All study participants provided written informed consent, and all centers obtained approval from their institutional review boards (13). Normative lung functions were predicted to calculate the % FEV1 in the UKB (supplement, UKB lung function methods, Table E2). Mild or moderate-severe asthma were defined by asthma medication use id codes in the UKB dataset [supplement] (12).

The Severe Asthma Research Program (SARP) includes 1283 asthmatics, including 703 non-severe and 580 severe asthmatics as previously defined [see supplement, SARP dataset] (14, 15). Exclusion criteria similar to those in the UKB are detailed in the supplement. Fractional exhaled nitric oxide (FENO) was measured by an online method at a constant flow rate of 50 ml/second according to the standards published by ATS. Asthma exacerbation was defined as an event of a worsening of asthma requiring treatment with systemic corticosteroids for three or more days to prevent a serious outcome. Asthma exacerbations in the prior year were determined by participant recall for the past twelve months. For exacerbations in the one-year longitudinal follow up, exacerbation events were collected using standardized questionnaires at in-person visits or during phone calls. Those with three or more exacerbations in the one-year follow up have been designated as an exacerbation-prone asthma phenotype (16). The study protocol and procedures were approved by the IRB at each participating center and by an independent Data Safety Monitoring Board. All participants provided written informed consent. The study is registered on ClinicalTrials.gov [NCT01750411].

Mitochondrial DNA Copy Number (mtDNA-CN)

DNA used for mtDNA-CN estimation in the UKB was derived from buffy coat. mtDNA-CN was estimated from a combination of whole exome sequencing data and probe intensities from the Affymetrix Axiom array (13). DNA used for mtDNA-CN estimation in SARP was derived from whole blood. Samples underwent whole genome sequencing (WGS) at the New York Genome Center of TOPMed (https://topmed.nhlbi.nih.gov/). mtDNA-CN was estimated from WGS by comparing the coverage of mitochondrial reads to autosomal reads using the program fastMitoCalc in mitoCaller (17). Raw unadjusted means of mtDNA-CN are provided for the UKB and SARP and termed mtDNA-CN. Because there is variation in mtDNA-CN methodology (Affymetrix Axiom array in the UKB as opposed to mtDNA/nDNA in TOPMed) and because depending on study, mtDNA-CN and age, sex, and white blood cell (WBC) measurement can vary, we calculated mtDNA-CN residuals after adjusting for covariates that may influence mtDNA-CN and standardized these residuals to a mean of 0 and SD of 1. These include age, sex, center where sample was collected, and WBC as well as neutrophils, eosinophils, and/or platelets when available (Table 1) (see supplement, Adjusted mtDNA-CN).

View this table:
Table 1. Covariates used to adjust mtDNA-CN

Serum Antioxidant Analyses

Superoxide dismutase (SOD) activity and glutathione peroxidase activity (GPx) were determined as previously described (18, 19).

Statistical Analysis

Analyses were conducted with JMP Pro 15 (SAS Institute, Cary, NC) or R (version 4.2.2, cran-project.org). Data were summarized as mean ± SD or frequency (percentage) where appropriate. For comparisons of two groups, Student’s t-test or Wilcoxon test was used as appropriate. For three or more groups, analysis of variance (ANOVA) was used. Non-normally distributed variables were log transformed prior to linear regression analysis. The level of significance for P was chosen at 0.05.

Association with Asthma Phenotypes

Models were adjusted as appropriate to address hypotheses (Table 1) (see supplement, Adjusted mtDNA-CN). For the UKB, the differences in mtDNA-CN (95% confidence intervals [CI]) by presence of asthma and by asthma severity were evaluated using linear regression models adjusted for age, age2, sex, center, log-transformed neutrophil and eosinophil counts, and platelet count (supplement). To examine the relationship between unadjusted mtDNA-CN levels and age by asthma, models were used to include an interaction term between age and asthma, adjusting the same set of covariates (Table 1 and supplement). To also allow for non-linear association between age and mtDNA-CN, we used age as a continuous variable with restricted cubic splines with 3 knots (at 10th, 50th, and 90th percentiles). We performed separate analyses for asthma (no / yes) and for asthma severity (non-asthma / mild asthma / moderate-severe asthma). For SARP, the difference (95% CI) in mtDNA-CN by non-severe and severe asthma, and the difference (95% CI) in mtDNA-CN by exacerbation occurrence were examined using linear regression models adjusted for age, age2, sex, ethnicity, WBC and neutrophil counts, SARP cohort, and cohort and age interaction (Table 1 and supplement). The relationship between mtDNA-CN and age by asthma severity was examined using linear regression models adjusted for sex, ethnicity, WBC and neutrophil counts, and SARP cohort (Table 1 and supplement).

Results

Asthma in UK Biobank

The clinical characteristics of participants in the UKB are displayed in Table 2. The asthma group has a greater percentage of females, lower lung functions, and higher neutrophil and eosinophil counts (P < 0.001 for all) (Table 2), consistent with known characteristics of asthma. Based on medication use, among the asthmatics, 25,322 (85%) have mild asthma and 4,446 (15%) have moderate-severe asthma. Compared with mild asthmatics, moderate-severe asthmatics are older, have higher body mass index (BMI), lower lung functions, and higher eosinophil counts (Table 2), consistent with known characteristics of severe asthma. Asthmatics have lower unadjusted mtDNA-CN (Table 2). Even after taking into account covariates of age, sex, and blood cell counts (see Table 2 and supplement), mtDNA-CN is significantly lower in asthmatics compared to non-asthmatics (beta, -0.026 [95% CI, -0.038 to -0.014], P = 2.46×10-5) (Table 2) (Figure 1). Both mild and moderate-severe asthmatics have lower adjusted mtDNA-CN compared to non-asthmatics (mild asthma, beta, -0.024 [95% CI, -0.037 to -0.012], P = 1.90×10-4; moderate-severe asthma, beta, -0.034 [95% CI, -0.063 to -0.004], P = 0.02) (Figure 1). Moderate-severe asthmatics have lower unadjusted mtDNA-CN than mild asthmatics (Table 2), but there is no significant difference in mtDNA-CN by severity when adjusted for age, age2, sex, center, log-transformed neutrophil and eosinophil counts, and platelet counts (Table 2) (Figure 1). mtDNA-CN is inversely related to age in asthmatic and non-asthmatic individuals (all P < 0.001). The decline in mtDNA-CN with increasing age is similar among groups (all P > 0.05) (Figure 2A – D). Overall, mtDNA-CN is lower in asthmatics as compared to non-asthmatics and is independent of asthma severity defined by medication use. To more carefully assess mitochondrial dysfunction among non-severe and severe asthma, we examined mtDNA-CN SARP network data.

Figure 1.
Figure 1.

Asthmatics have lower mitochondrial DNA copy number (mtDNA-CN) compared to non-asthmatics in the UK Biobank (UKB). Forest plot depicting beta value and associated 95% confidence interval (CI) of mtDNA-CN in all asthmatics, mild or moderate-severe asthmatics as compared to non-asthmatics.

Figure 2.
Figure 2.

mtDNA-CN inversely related with age in the UKB. (A-B) The predicted mtDNA-CN residuals using the unadjusted model (including only non-linear age (restricted cubic splines with 3 knots), asthma status, and an interaction between non-linear age and asthma status) decreases with age in non-asthmatics and all asthmatics (A), and in mild asthmatics, and moderate-severe asthmatics (B). (C-D) The predicted mtDNA-CN residuals marginally adjusted for sex, center, log-transformed neutrophil and eosinophil counts, and platelet counts were significantly related with age in non-asthmatics and all asthmatics (C), and in mild asthmatics, and moderate-severe asthmatics (D).

View this table:
Table 2. Features of asthmatics and non-asthmatics in the UK Biobank (UKB)

mtDNA-CN in Asthmatics in SARP

SARP was designed to characterize severe asthma at the molecular, cellular, and clinical levels and includes detailed clinical and mechanistic data. Thus, this cohort allows investigation of mtDNA-CN in relation to asthma characteristics and mechanisms. The clinical characteristics of the participants of SARP are displayed in Table 3. Compared to UKB, the participants in SARP are younger, have higher BMI, worse lung functions, and higher eosinophil counts. The SARP cohort includes all racial and ethnic backgrounds. In the UKB, severity was defined based on asthma medication use in self-reported asthma; in contrast, in SARP, asthma was validated based upon American Thoracic Society (ATS) guidelines, which include a positive methacholine challenge test and/or reversible airflow obstruction, and asthma severity classified based on ATS and European Respiratory Society (ERS) guidelines (14, 15). As compared to non-severe asthmatics, severe asthmatics in SARP have lower lung functions, higher white blood cell (WBC), neutrophil, and eosinophil counts, and worse asthma control as measured with the Asthma Control Test (Table 3).

View this table:
Table 3. Features of participants in the Severe Asthma Research Program (SARP)

Variation in mtDNA-CN is observed by sex and ethnicity (Figure E1). After adjusting for age, sex, ethnicity, WBC and other covariates, severe and non-severe asthmatics have similar mtDNA-CN (Table 3). BMI is related to mtDNA-CN among all asthma (r = 0.060, P = 0.03, n = 1237). Linear regression analyses suggests that mtDNA-CN in severe asthma may tend to decline more with age than in non-severe asthma (P = 0.09 severe asthmatics vs. non-severe asthmatics.) (Figure 3A). However, when mtDNA-CN is adjusted for covariates, age-related decline of mtDNA-CN is no longer observed in this cohort that is smaller in number and relatively younger than the UKB cohort (Figure 3B).

Figure 3.
Figure 3.

mtDNA-CN related with age in SARP. (A) mtDNA-CN residuals are inversely related to age in severe asthmatics and all asthmatics but not in non-severe asthmatics, (B) while adjusted mtDNA-CN is not related to age in severe or non-severe asthmatics.

Increased Oxidative Stress in Asthma of SARP

A small number of healthy individuals without asthma were included in SARP as previously described (Table E3) (20), but these do not have mtDNA-CN available. Here, asthmatics have greater nitrative stress than healthy controls as measured by FENO as previously reported (Table E3) (20, 21). In an early report from SARP with 114 asthmatics and 20 healthy controls (22), asthmatics were reported to have lower serum SOD as compared with controls (4, 5, 18, 22). Similarly here, in the larger SARP cohort, SOD activity is lower in asthma as compared to healthy controls (SOD U/ml, controls 22.3 ± 10.8, n = 114, asthmatics 19.1 ± 10.7, n = 960; P = 0.003) (Table E3), but severe and non-severe asthma have similar SOD (Table 3). GPx is similar among asthmatics and controls, as previously published in a much smaller cohort (22) (GPx U/ml, controls 0.087 ± 0.020, n = 44, asthmatics 0.087 ± 0.023, n = 297; P = 0.9) (Table E3). These results support the paradigm that asthma inflammation is characterized by changes in redox balance as compared to healthy controls (4, 5, 23), but the severe and non-severe groups are similar in the measures of SOD, GPx, FENO, and mtDNA-CN (Table 3).

mtDNA-CN Predicts Asthma Exacerbations

Predicting and preventing exacerbation frequency is critical for the treatment of asthma (16, 24). Here, participants with asthma in one-year longitudinal follow-up in SARP III (n = 468) are analyzed for associations of mtDNA-CN and SOD with exacerbations. Based on previous studies in SARP (16, 25), asthmatics are classified into three groups [0, 1 – 2, or ≥3] by the number of exacerbations in one-year longitudinal follow-up (Table 4). As previously reported (16), of 468 asthmatics, 281 (60%) were without exacerbation in one-year longitudinal follow-up, while 126 (27%) had one or two exacerbations, and 61 (13%) had three or more exacerbations (Table 4). Fifty-five of the 61 (90.2%) asthmatics of the group with three or more exacerbations have severe disease; 131 of the 281 (46.6%) asthmatics of the group with no exacerbations, and 91 of the 126 (72.2%) asthmatics of the group with one or two exacerbations have severe disease (Table 4). As previously reported (16), FENO levels are different among the three groups (ANOVA P = 0.01) (Table 4). Here, mtDNA-CN is lowest in the group with three or more exacerbations in longitudinal follow-up (Table 4). In contrast to association with future exacerbations, the mtDNA-CN is not associated with numbers of exacerbations in the prior year [1-2 vs. 0 exacerbation group, P = 0.2; 3 or more vs. 0 exacerbation group, P = 0.7), suggesting that high mtDNA-CN levels may be protective from future exacerbations.

View this table:
Table 4. Features of study participants categorized by frequency of exacerbations‡ in SARP

To study the association of mtDNA-CN to the likelihood that asthma exacerbation will occur, odds ratios (OR) are calculated. In evaluation of mtDNA quartiles, the top quartile of mtDNA-CN has significantly lower odds (OR 0.257 [95% CI, 0.099 to 1.670], P = 0.005) of three or more exacerbations as compared with the bottom quartile (Figure E2). The top-quartile mtDNA-CN individuals have features similar to the bottom-quartile, including the numbers that have severe asthma, but differ in having a higher BMI (P = 0.005) (Table E4).

To determine cutoff values of mtDNA, the dichotomized biomarker at the cutoffs e.g., -0.033 cutoff for mtDNA, are examined to analyze the odds of having three or more exacerbations vs. zero exacerbations (Table E5). Asthmatics with mtDNA-CN below the cutoff are older than those with mtDNA-CN above the cutoff (P = 0.04) (Table E6). The odd ratio of mtDNA-CN is 0.352 ([95% CI, 0.164 to 0.753], P = 0.007) (Figure 4), indicating that asthmatics with levels of mtDNA-CN above the cutoff have ∼35% of the risk of three or more exacerbations compared to those with mtDNA-CN levels below the cutoff. Similarly, using continuous analyses per SD increase, mtDNA-CN predicts those who will have few exacerbations (OR 0.442 [95% CI, 0.224 to 0.875], P = 0.01), while FENO predicts those who will have more exacerbations (OR 1.387 [95% CI, 1.090 to 1.764], P = 0.008) (Figure E3B). FENO levels (higher than cutoff of 27.5 ppb) increase exacerbation risk twofold (P = 0.009) (Figure 4) as previously reported (16, 21). SOD is not predictive of exacerbations.

Figure 4.
Figure 4.

Odds ratios for risk of three or more exacerbations in the subsequent year for biomarkers of asthma. Cutoff values were determined (Table E5) for mtDNA-CN, serum superoxide dismutase activity, and fractional exhaled nitric oxide. A higher mtDNA-CN predicted a lower risk of three or more exacerbations, while a higher fractional exhaled nitric oxide predicted a greater risk of three or more exacerbations.

Discussion

This study shows that asthma is associated with low mtDNA-CN, an indicator of mitochondrial dysfunction. Interestingly, although asthmatics have lower mtDNA-CN as compared to non-asthmatics, the asthma severity groups did not differ with respect to mtDNA-CN. However, higher mtDNA-CN predicts less frequent exacerbations in one-year follow-up, suggesting mitochondrial dysfunction contributes to the exacerbation-prone asthma phenotype.

Mitochondria are principal producers of ROS and regulate redox status in the cell. Due to multiple copies of mtDNA, lack of a protective histone protein, and limited mtDNA repair capacity, mtDNA is also more susceptible to damage by excess ROS, which may lead to dysregulation of the respiratory chain and further increase ROS production (26-30). The non-coding regulatory D-loop region (mtDNA control region) is known to be more sensitive to ROS-mediated damage than the coding portion of the mitochondrial genome and the nuclear genome (31, 32). Increased mtDNA D-loop methylation induced by excess ROS was associated with reduced mtDNA-CN (33). In the Study of Asthma Phenotypes and Pharmacogenomic Interactions by Race-Ethnicity (SAPPHIRE), mtDNA-CN in whole blood was higher in African American asthmatics as compared with their non-asthmatic counterparts, but mtDNA-CN was not related to asthma severity (34). Here, we also found mtDNA-CN is not significantly different between severe and non-severe asthma after adjusting for age, sex, ethnicity, and blood cell counts in both UKB and SARP cohorts; however, asthmatics, regardless of severity, have significantly lower mtDNA-CN as compared with non-asthmatics. The UKB and SARP cohorts in this study differ from the SAPPHIRE cohort in genetic ancestries, ages, smoking status, BMI, WBC counts, and sample size of studies, all of which might lead to the discrepant findings. Here, asthmatics with African ancestry have higher levels than those with European ancestry. Further studies are needed to understand the influence of ancestry and mitochondrial function.

Previous studies indicated that mtDNA-CN is affected by and should be adjusted for age (35-38). Our results show that mtDNA-CN is inversely related to age in all asthmatics but does not account for differences among asthmatics and non-asthmatics. Asthmatics have lower mtDNA-CN at every age in comparison to non-asthmatics. The mtDNA-CN per cell is estimated from genomic DNA derived from buffy coat, which greatly depends on peripheral blood cell type, platelet count, and WBC counts (35, 39, 40). Previously, others reported that mtDNA-CN was inversely correlated with WBC count in their study populations (35, 39, 40). Here, WBC and neutrophil counts are significantly higher in severe asthmatics than in non-severe asthmatics and inversely related with mtDNA-CN. In SARP, blood eosinophil counts are higher in severe asthmatics than in non-severe asthmatics but are not related with mtDNA-CN, indicating that lower levels of mtDNA-CN may be related to neutrophilic inflammation and not Th2 inflammation. Mitochondria numbers can vary in different cells and/or tissues. Prior work has shown changes in mitochondrial function in platelets and airway epithelial cells from asthmatics (41). Interestingly, although asthma is associated with obesity (42), resting energy expenditure is elevated in asthmatics as compared to non-asthmatics (43), supporting the concept of changes in mitochondrial function.

Asthma is a common problem. Severe asthma consumes a great deal of resources. However, there are very few biomarkers at the molecular level that discriminate phenotypes of asthma, relative risks for asthma severity, and/or risks of frequent exacerbations. Previous findings showed that FENO (44, 45), blood eosinophils (16, 24, 45, 46), sputum eosinophils (16, 44), body mass index (16, 47), age (16, 47), and urinary total conjugated 3-bromotyrosine (48) are associated with uncontrolled asthma and exacerbation frequency. Existing biomarkers can be grouped into categories, e.g., related to redox pathways, arginine/NO pathway, and/or Th2 eosinophilic pathway (16, 24, 44-47, 49, 50). Although frequent exacerbations happen in a small number of patients, prediction of asthma exacerbations is important for optimal care of individuals (16, 24). Prior reports in SARP cohort indicate that blood eosinophils, BMI, and bronchodilator responsiveness were associated with exacerbation frequency in the prior year (16); high plasma IL-6 and blood eosinophils were associated with increased incident rate of exacerbations over three years of follow-up (24), and urinary total conjugated 3-bromotyrosine predicted exacerbation risk over one year of follow-up (48). Here, asthmatics with high mtDNA-CN levels have a 65% reduced risk of having three or more exacerbations than those with lower mtDNA-CN. Unlike other measures, we did not find any association of mtDNA-CN to prior year exacerbations. This may suggest that changes in mtDNA-CN occur in advance of exacerbations. The loss of SOD activity in asthma is well described and related to oxidative inactivation due to cleavage of the C57-C146 disulfide bond and exposure of usually unavailable cysteines (6). However, loss of SOD activity is not associated with exacerbations. A prior study suggests a relationship of SOD to severity of airflow obstruction (22).

Overall, this study identifies that mitochondrial dysfunction is present in asthma and may play a role in the molecular basis of asthma and exacerbations. Higher mtDNA-CN, an indicator of better mitochondrial function, is predictive of fewer future exacerbations, suggesting that targeting mitochondria and/or metabolism may be useful in future approaches for treatment of exacerbation-prone asthma.

Data Availability

All data produced in the present study are available upon reasonable request to the authors.

Acknowledgments

We gratefully acknowledge the studies and participants who provided biological samples and data for SARP, UKB, and TOPMed. The authors would like to thank SARP clinical coordinators and laboratory personnel for their contribution to fulfilling SARP’s scientific mission. The research findings resulting from SARP would not have been possible without their dedication and assistance with study visits, data collection, and biological sample processing. The authors also thank Peter Bazeley, Thomas W. Blackwell, Xiaoqi Geng, Goncalo Abecasis, Jun Ding for mitochondrial DNA copy number estimation.

SARP3 acknowledgement (industry): The following companies provided financial support for study activities at the Coordinating and Clinical Centers beyond the third year of patient follow-up: AstraZeneca, Boehringer-Ingelheim, Genentech, GlaxoSmithKline, Sanofi–Genzyme–Regeneron, and TEVA.

Footnotes

  • Conflicts of interest: The authors have declared that no conflict of interest exists.

  • Our study identifies that mitochondrial function plays a role in the molecular basis of asthma. Higher mtDNA-CN, an indicator of mitochondrial function, is predictive of fewer future exacerbations, suggesting that targeting mitochondria and/or metabolism may be useful in future approaches for treatment of exacerbation-prone asthma.

References

  1. 1.↵
    Comhair SA, Bhathena PR, Dweik RA, Kavuru M, Erzurum SC. Rapid loss of superoxide dismutase activity during antigen-induced asthmatic response. Lancet 2000; 355: 624.
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.
    Wu W, Samoszuk MK, Comhair SA, Thomassen MJ, Farver CF, Dweik RA, et al. Eosinophils generate brominating oxidants in allergen-induced asthma. J Clin Invest 2000; 105: 1455–1463.
    OpenUrlCrossRefPubMedWeb of Science
  3. 3.↵
    Xu W, Ghosh S, Comhair SA, Asosingh K, Janocha AJ, Mavrakis DA, et al. Increased mitochondrial arginine metabolism supports bioenergetics in asthma. J Clin Invest 2016; 126: 2465–2481.
    OpenUrlCrossRefPubMed
  4. 4.↵
    Comhair SA, Erzurum SC. Redox control of asthma: molecular mechanisms and therapeutic opportunities. Antioxid Redox Signal 2010; 12: 93–124.
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.↵
    Erzurum SC. New Insights in Oxidant Biology in Asthma. Ann Am Thorac Soc 2016; 13 Suppl 1: S35–39.
    OpenUrl
  6. 6.↵
    Ghosh S, Willard B, Comhair SA, Dibello P, Xu W, Shiva S, et al. Disulfide bond as a switch for copper-zinc superoxide dismutase activity in asthma. Antioxid Redox Signal 2013; 18: 412–423.
    OpenUrlCrossRefPubMed
  7. 7.↵
    Reddy PH. Mitochondrial Dysfunction and Oxidative Stress in Asthma: Implications for Mitochondria-Targeted Antioxidant Therapeutics. Pharmaceuticals (Basel) 2011; 4: 429–456.
    OpenUrl
  8. 8.↵
    Robin ED, Wong R. Mitochondrial DNA molecules and virtual number of mitochondria per cell in mammalian cells. J Cell Physiol 1988; 136: 507–513.
    OpenUrlCrossRefPubMedWeb of Science
  9. 9.↵
    Wallace DC. A mitochondrial bioenergetic etiology of disease. J Clin Invest 2013; 123: 1405–1412.
    OpenUrlCrossRefPubMedWeb of Science
  10. 10.↵
    Jeng JY, Yeh TS, Lee JW, Lin SH, Fong TH, Hsieh RH. Maintenance of mitochondrial DNA copy number and expression are essential for preservation of mitochondrial function and cell growth. J Cell Biochem 2008; 103: 347–357.
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    Castellani CA, Longchamps RJ, Sun J, Guallar E, Arking DE. Thinking outside the nucleus: Mitochondrial DNA copy number in health and disease. Mitochondrion 2020; 53: 214–223.
    OpenUrlCrossRef
  12. 12.↵
    Shrine N, Portelli MA, John C, Soler Artigas M, Bennett N, Hall R, et al. Moderate-to-severe asthma in individuals of European ancestry: a genome-wide association study. Lancet Respir Med 2019; 7: 20–34.
    OpenUrlCrossRef
  13. 13.↵
    Longchamps RJ, Yang SY, Castellani CA, Shi W, Lane J, Grove ML, et al. Genome-wide analysis of mitochondrial DNA copy number reveals loci implicated in nucleotide metabolism, platelet activation, and megakaryocyte proliferation. Hum Genet 2022; 141: 127–146.
    OpenUrl
  14. 14.↵
    Proceedings of the ATS workshop on refractory asthma: current understanding, recommendations, and unanswered questions. American Thoracic Society. Am J Respir Crit Care Med 2000; 162: 2341-2351.
    OpenUrlCrossRefPubMed
  15. 15.↵
    Chung KF, Wenzel SE, Brozek JL, Bush A, Castro M, Sterk PJ, et al. International ERS/ATS guidelines on definition, evaluation and treatment of severe asthma. Eur Respir J 2014; 43: 343–373.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    Denlinger LC, Phillips BR, Ramratnam S, Ross K, Bhakta NR, Cardet JC, et al. Inflammatory and Comorbid Features of Patients with Severe Asthma and Frequent Exacerbations. Am J Respir Crit Care Med 2017; 195: 302–313.
    OpenUrl
  17. 17.↵
    Ding J, Sidore C, Butler TJ, Wing MK, Qian Y, Meirelles O, et al. Assessing Mitochondrial DNA Variation and Copy Number in Lymphocytes of ∼2,000 Sardinians Using Tailored Sequencing Analysis Tools. PLoS Genet 2015; 11: e1005306.
    OpenUrlCrossRefPubMed
  18. 18.↵
    Comhair SA, Xu W, Ghosh S, Thunnissen FB, Almasan A, Calhoun WJ, et al. Superoxide dismutase inactivation in pathophysiology of asthmatic airway remodeling and reactivity. Am J Pathol 2005; 166: 663–674.
    OpenUrlCrossRefPubMedWeb of Science
  19. 19.↵
    Comhair SA, Erzurum SC. Antioxidant responses to oxidant-mediated lung diseases. Am J Physiol Lung Cell Mol Physiol 2002; 283: L246–255.
    OpenUrlCrossRefPubMedWeb of Science
  20. 20.↵
    Dweik RA, Sorkness RL, Wenzel S, Hammel J, Curran-Everett D, Comhair SA, et al. Use of exhaled nitric oxide measurement to identify a reactive, at-risk phenotype among patients with asthma. Am J Respir Crit Care Med 2010; 181: 1033–1041.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    Dweik RA, Boggs PB, Erzurum SC, Irvin CG, Leigh MW, Lundberg JO, et al. An official ATS clinical practice guideline: interpretation of exhaled nitric oxide levels (FENO) for clinical applications. Am J Respir Crit Care Med 2011; 184: 602–615.
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.↵
    Comhair SA, Ricci KS, Arroliga M, Lara AR, Dweik RA, Song W, et al. Correlation of systemic superoxide dismutase deficiency to airflow obstruction in asthma. Am J Respir Crit Care Med 2005; 172: 306–313.
    OpenUrlCrossRefPubMedWeb of Science
  23. 23.↵
    Dweik RA, Laskowski D, Abu-Soud HM, Kaneko F, Hutte R, Stuehr DJ, et al. Nitric oxide synthesis in the lung. Regulation by oxygen through a kinetic mechanism. J Clin Invest 1998; 101: 660–666.
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.↵
    Peters MC, Mauger D, Ross KR, Phillips B, Gaston B, Cardet JC, et al. Evidence for Exacerbation-Prone Asthma and Predictive Biomarkers of Exacerbation Frequency. Am J Respir Crit Care Med 2020; 202: 973–982.
    OpenUrl
  25. 25.↵
    Moore WC, Bleecker ER, Curran-Everett D, Erzurum SC, Ameredes BT, Bacharier L, et al. Characterization of the severe asthma phenotype by the National Heart, Lung, and Blood Institute’s Severe Asthma Research Program. J Allergy Clin Immunol 2007; 119: 405–413.
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.↵
    Enright H, Miller WJ, Hays R, Floyd RA, Hebbel RP. Preferential targeting of oxidative base damage to internucleosomal DNA. Carcinogenesis 1996; 17: 1175–1177.
    OpenUrlCrossRefPubMedWeb of Science
  27. 27.
    Enright HU, Miller WJ, Hebbel RP. Nucleosomal histone protein protects DNA from iron-mediated damage. Nucleic Acids Res 1992; 20: 3341–3346.
    OpenUrlCrossRefPubMedWeb of Science
  28. 28.
    Shokolenko IN, Alexeyev MF. Mitochondrial DNA: A disposable genome? Biochim Biophys Acta 2015; 1852: 1805–1809.
    OpenUrlCrossRefPubMed
  29. 29.
    Shokolenko IN, Wilson GL, Alexeyev MF. Persistent damage induces mitochondrial DNA degradation. DNA Repair (Amst) 2013; 12: 488–499.
    OpenUrlCrossRefPubMed
  30. 30.↵
    Yakes FM, Van Houten B. Mitochondrial DNA damage is more extensive and persists longer than nuclear DNA damage in human cells following oxidative stress. Proc Natl Acad Sci U S A 1997; 94: 514–519.
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    Mambo E, Gao X, Cohen Y, Guo Z, Talalay P, Sidransky D. Electrophile and oxidant damage of mitochondrial DNA leading to rapid evolution of homoplasmic mutations. Proc Natl Acad Sci U S A 2003; 100: 1838–1843.
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    Pastukh VM, Gorodnya OM, Gillespie MN, Ruchko MV. Regulation of mitochondrial genome replication by hypoxia: The role of DNA oxidation in D-loop region. Free Radic Biol Med 2016; 96: 78–88.
    OpenUrl
  33. 33.↵
    Zheng LD, Linarelli LE, Liu L, Wall SS, Greenawald MH, Seidel RW, et al. Insulin resistance is associated with epigenetic and genetic regulation of mitochondrial DNA in obese humans. Clin Epigenetics 2015; 7: 60.
  34. 34.↵
    Cocco MP, White E, Xiao S, Hu D, Mak A, Sleiman P, et al. Asthma and its relationship to mitochondrial copy number: Results from the Asthma Translational Genomics Collaborative (ATGC) of the Trans-Omics for Precision Medicine (TOPMed) program. PLoS One 2020; 15: e0242364.
    OpenUrl
  35. 35.↵
    Castellani CA, Longchamps RJ, Sumpter JA, Newcomb CE, Lane JA, Grove ML, et al. Mitochondrial DNA copy number can influence mortality and cardiovascular disease via methylation of nuclear DNA CpGs. Genome Med 2020; 12: 84.
    OpenUrl
  36. 36.
    Longchamps RJ, Castellani CA, Yang SY, Newcomb CE, Sumpter JA, Lane J, et al. Evaluation of mitochondrial DNA copy number estimation techniques. PLoS One 2020; 15: e0228166.
    OpenUrlCrossRef
  37. 37.
    Ashar FN, Moes A, Moore AZ, Grove ML, Chaves PHM, Coresh J, et al. Association of mitochondrial DNA levels with frailty and all-cause mortality. J Mol Med (Berl) 2015; 93: 177–186.
    OpenUrlCrossRefPubMed
  38. 38.↵
    Maremanda KP, Sundar IK, Li D, Rahman I. Age-Dependent Assessment of Genes Involved in Cellular Senescence, Telomere, and Mitochondrial Pathways in Human Lung Tissue of Smokers, COPD, and IPF: Associations With SARS-CoV-2 COVID-19 ACE2-TMPRSS2-Furin-DPP4 Axis. Front Pharmacol 2020; 11: 584637.
    OpenUrlCrossRefPubMed
  39. 39.↵
    Hulgan T, Kallianpur AR, Guo Y, Barnholtz-Sloan JS, Gittleman H, Brown TT, et al. Peripheral Blood Mitochondrial DNA Copy Number Obtained From Genome-Wide Genotype Data Is Associated With Neurocognitive Impairment in Persons With Chronic HIV Infection. J Acquir Immune Defic Syndr 2019; 80: e95–e102.
    OpenUrl
  40. 40.↵
    Knez J, Winckelmans E, Plusquin M, Thijs L, Cauwenberghs N, Gu Y, et al. Correlates of Peripheral Blood Mitochondrial DNA Content in a General Population. Am J Epidemiol 2016; 183: 138–146.
    OpenUrlCrossRefPubMed
  41. 41.↵
    Xu W, Cardenes N, Corey C, Erzurum SC, Shiva S. Platelets from Asthmatic Individuals Show Less Reliance on Glycolysis. PLoS One 2015; 10: e0132007.
    OpenUrlCrossRef
  42. 42.↵
    Taylor B, Mannino D, Brown C, Crocker D, Twum-Baah N, Holguin F. Body mass index and asthma severity in the National Asthma Survey. Thorax 2008; 63: 14–20.
    OpenUrlAbstract/FREE Full Text
  43. 43.↵
    Mey JT, Matuska B, Peterson L, Wyszynski P, Koo M, Sharp J, et al. Resting Energy Expenditure Is Elevated in Asthma. Nutrients 2021; 13.
  44. 44.↵
    Wedes SH, Wu W, Comhair SA, McDowell KM, DiDonato JA, Erzurum SC, et al. Urinary bromotyrosine measures asthma control and predicts asthma exacerbations in children. J Pediatr 2011; 159: 248–255 e241.
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    Malinovschi A, Fonseca JA, Jacinto T, Alving K, Janson C. Exhaled nitric oxide levels and blood eosinophil counts independently associate with wheeze and asthma events in National Health and Nutrition Examination Survey subjects. J Allergy Clin Immunol 2013; 132: 821–827 e821-825.
    OpenUrlCrossRefWeb of Science
  46. 46.↵
    Wenzel S, Ford L, Pearlman D, Spector S, Sher L, Skobieranda F, et al. Dupilumab in persistent asthma with elevated eosinophil levels. N Engl J Med 2013; 368: 2455–2466.
    OpenUrlCrossRefPubMedWeb of Science
  47. 47.↵
    Zein JG, Wu CP, Attaway AH, Zhang P, Nazha A. Novel Machine Learning Can Predict Acute Asthma Exacerbation. Chest 2021; 159: 1747–1757.
    OpenUrl
  48. 48.↵
    Wang Z, Xu W, Comhair SAA, Fu X, Shao Z, Bearden R, et al. Urinary total conjugated 3-bromotyrosine, asthma severity, and exacerbation risk. Am J Physiol Lung Cell Mol Physiol 2022; 323: L548–L557.
    OpenUrl
  49. 49.↵
    Comhair SAA, Bochenek G, Baicker-McKee S, Wang Z, Stachura T, Sanak M, et al. The utility of biomarkers in diagnosis of aspirin exacerbated respiratory disease. Respir Res 2018; 19: 210.
    OpenUrl
  50. 50.↵
    Price DB, Rigazio A, Campbell JD, Bleecker ER, Corrigan CJ, Thomas M, et al. Blood eosinophil count and prospective annual asthma disease burden: a UK cohort study. Lancet Respir Med 2015; 3: 849–858.
    OpenUrl
Back to top
Posted December 05, 2023.
Download PDF

Supplementary Material

Data/Code
Email
Mitochondrial DNA Copy Number Variation in Asthma Risk, Severity, and Exacerbations
Weiling Xu, Yun Soo Hong, Bo Hu, Suzy A. A. Comhair, Allison J. Janocha, Joe G. Zein, Ruoying Chen, Deborah A. Meyers, David T. Mauger, Victor E. Ortega, Eugene R. Bleecker, Mario Castro, Loren C. Denlinger, John V. Fahy, Elliot Israel, Bruce D. Levy, Nizar N. Jarjour, Wendy C. Moore, Sally E. Wenzel, Benjamin Gaston, Chunyu Liu, Dan E. Arking, Serpil C. Erzurum, the National Heart, Lung, and Blood Institute (NHLBI) Severe Asthma Research Program (SARP), TOPMed mtDNA Working Group in NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium
medRxiv 2023.12.05.23299392; doi: https://doi.org/10.1101/2023.12.05.23299392
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area