Abstract
Purpose. This study aimed to investigate the relationship of partial pressure of carbon dioxide (PaCO2) with BODE and GOLD in stable COPD subjects and to explore the predictive value of PaCO2 for severe COPD (BODE index score ≥5 or GOLD index score ≥3). Patients and Methods. In total, 80 participants with COPD and free from other conditions affecting PaCO2 were recruited. Arterial blood gases, BODE, GOLD, SGRQ, lung function, and other data were collected. The BODE index was calculated, and patients were divided into two groups according to the BODE index and PaCO2 median, respectively. We used Pearson’s correlation test and the receiver operating characteristic curves to evaluate the utility of PaCO2. Besides, the univariate and multivariate logistic regression analyses were conducted to verify whether PaCO2 was an independent factor associated with BODE grades. Results. COPD subjects with BODE ≥5 and GOLD ≥3 had significantly higher levels of PaCO2 ( = 0.004, = 0.001, respectively). In the high PaCO2 group, patients underwent poorer outcomes than the low PaCO2 group. PaCO2 was negatively correlated with forced expiratory volume in 1 second in percent of the predicted value (FEV1%) (r = −0.612, < 0.001). The performance of PaCO2 levels in predicting BODE ≥5 and GOLD ≥3 was 0.748 and 0.755, respectively. The logistic regression analyses proved that PaCO2 was associated with BODE ≥5 in COPD patients (odds ratio = 1.160, 95% CI: 1.025–1.313, = 0.019). Conclusions. A higher level of PaCO2 was associated with a higher index for BODE or GOLD in COPD and had the predictive value for severe COPD.
1. Introduction
Chronic obstructive pulmonary disease (COPD), characterized by persistent respiratory symptoms, airflow limitation, and recurrent exacerbations [1, 2], is currently one of the leading causes of morbidity and mortality worldwide in the adult population [3]. Murphy et al. reported that persistent hypercapnia after an exacerbation was associated with excess mortality and early rehospitalization [2]. In the study of Dreher et al., a quarter of patients with COPD GOLD III and IV had chronic hypercapnic respiratory failure (CHRF) [4]. The previous study proved the benefits of domiciliary noninvasive ventilation (NIV) in COPD patients with CHRF in terms of COPD Assessment Test (CAT), arterial blood gases (including PaCO2), the number of acute exacerbations, and the BODE index of the first 6 months [5]. Home noninvasive ventilation combined with home oxygen therapy prolonged the time to readmission or death among individuals with persistent hypercapnia in acute exacerbation of COPD (AECOPD) [2]. Moreover, a clinical trial reported that the technique veno-venous extracorporeal CO2 removal (vv-ECCO2R) had the potential to rapidly correct respiratory acidosis and could decrease elevated mean pulmonary artery pressure values in severe COPD significantly by reducing high PaCO2 [6]. Almagro et al. stated that partial arterial carbon dioxide tension (PaCO2) at discharge was one of the significant predictors of mortality [7]. Additionally, Quintana et al. showed that elevated PaCO2 and decreased pH upon emergency department were predictors of ICU or IRCU admission [8]. Multidimensional indexes, such as the body mass index, airflow obstruction, dyspnea, and exercise capacity (BODE index), characterize the severity of COPD in a complex way (based on body mass index (BMI), forced expiratory volume in 1 second (FEV1), modified Medical Research Council (mMRC) dyspnea scores, and 6-minute walking distance (6MWD)) [9, 10]. The Global Initiative for Chronic Obstructive Lung Disease (GOLD) stages define the varying severity of COPD [11]. Worse disease severity of COPD can be defined by BODE ≥5 or GOLD ≥3 [9, 11–13]. COPD with hypercapnia was related to severely impaired lung function, and the negative correlation between PaCO2 and discharge ambulation distance may reflect the severity of airway obstruction [14]. However, no research has explored the association between PaCO2 and severity of patients in stable COPD directly.
Taking these facts into account, the aims of this study were to investigate the relationship of PaCO2 with BODE and GOLD in stable COPD subjects and to explore the predictive value of PaCO2 for severe COPD (BODE index score ≥5 or GOLD index score ≥3) [9, 11–13].
2. Materials and Methods
2.1. Study Population
In this cross-sectional research, 80 stable COPD patients were enrolled from the Department of Respiratory Medicine, the Third Affiliated Hospital of Wenzhou Medical University. The recruitment started in February 2018 and ended in February 2019. Inclusion criteria were as follows: (1) age more than 40 years and (2) diagnosis of COPD as defined in the Global Initiative for Chronic Obstructive Lung Disease guidelines; exclusion criteria were as follows: (1) bronchiectasis; (2) asthma; (3) post-tuberculosis sequelae; (4) malignant tumor; (5) chronic heart failure; (6) hepatic and renal insufficiency; (7) clinically significant chest wall deformity; and (8) neuromuscular weakness.
The study protocol was approved by the Ethics Committee of the Third Affiliated Hospital of Wenzhou Medical University, and the registration number of the Ethics Committee was YJ20170015. The study was conducted in accordance with the Declaration of Helsinki. Written informed consent was obtained from all study subjects.
2.2. Data Collection
Data such as age, gender, duration of disease, smoking status, and body mass index (BMI) were recorded upon hospital admission. We collected blood samples on 24 hours of admission to analyze blood routine parameters and blood biochemistry. Freshly drawn arterial blood was used for the measurement of pH, partial pressure of carbon dioxide, and oxygen in arterial blood (PaCO2 and PaO2). Spirometry tests were taken to measure FEV1, forced vital capacity (FVC), FEV1/FVC, and FEV1 in percent of the predicted value (FEV1%), which were significant indicators of baseline severity of COPD. Spirometry was conducted to determine the GOLD grades of airflow limitation, which was performed after the administration of an adequate dose of a short-acting inhaled bronchodilator to minimize variability [15]. GOLD 1-2 (FEV1 ≥ 50% predicted) indicate low risk, while GOLD 3-4 (FEV1 < 50% predicted) indicate high risk. Moreover, composite markers of disease including GOLD (Global Initiative for Chronic Obstructive Lung Disease), BODE (BMI, airway obstruction, dyspnea, and exercise capacity ), mMRC (modified Medical Research Council ), and SGRQ (St. George’s Respiratory Questionnaire) were evaluated for every participant. Six-minute walk distance (6MWD) tests were conducted according to the ATS guidelines [16]. It measures the distance that a patient can quickly walk on a 100 ft hallway in a period of 6 minutes, which reflects the exercise capacity for daily physical activities [16].
2.3. Statistical Analysis
All statistical analysis was performed using SPSS 25.0 (IBM Analytics). Normally distributed variables were presented as mean ± SD, while variables with skewed distribution were expressed as median (25th–75th percentile). We conducted the independent-samples t-test and the Mann–Whitney U-test to compare the differences in clinical characteristics between the low BODE group and the high BODE group, and the low PaCO2 group and the high PaCO2 group. Categorical variables were expressed as counts and percentages, and the chi-squared (χ2) test or Fisher’s exact test was performed for the intergroup comparison. The correlations between PaCO2 and numerical variables were estimated using Pearson’s or Spearman’s correlation coefficients. Pearson’s correlation analysis was used to evaluate the correlation between PaCO2 and FEV1%. The predictive accuracy of PaCO2 for higher GOLD and BODE index was determined by calculating the area under ROC curves (AUC). To identify the risk factors associated with COPD severity (higher GOLD and BODE index), we performed univariate analyses. Additionally, the multivariate logistic regression models were carried out to explain the contribution of PaCO2 in severe COPD with related confounders adjusted. All effects were considered significant at < 0.05.
3. Results
3.1. Baseline Characteristics of the Study Participants with Different COPD Severities
We categorized 80 subjects into two groups according to the BODE index (BODE <5, n = 50; BODE ≥5, n = 23, 7 patients without BODE index). As shown in Table 1, there were no statistically significant differences in age, sex, duration of disease, proportions of current smokers, pH on day 7 of admission, and PaO2 between two groups. Statistical differences were found when it comes to BMI, pH, and FEV1%. In addition, the indicator PaCO2 of the high BODE group was remarkably higher than the low BODE group, which is described in Figure 1(a) (41.77 ± 5.22 vs. 49.75 ± 11.55, = 0.004). Moreover, the high GOLD group had a drastically higher level of PaCO2 compared with the low GOLD group (39.69 ± 3.89 vs. 46.47 ± 9.08, = 0.001, Figure 1(b)).

(a)

(b)

(c)

(d)

(e)

(f)
For the purpose of investigating the relationship between PaCO2 and disease severity in COPD patients further, subjects were divided into two groups according to the median of PaCO2 (PaCO2 < 42.7, n = 39; PaCO2 ≥ 42.7, n = 41). Demographic variables such as age, gender, duration of disease, and proportions of current smokers were insignificantly different between the two groups, as well as the blood routine parameter WBC. There were also no statistical significances in terms of BMI, albumin, mMRC, and PaO2. Patients with higher levels of PaCO2 had lower levels of pH, pH on day 7 of admission, and all indicators of pulmonary function. Figure 1(c) shows that subjects with a higher level of PaCO2 had a lower level of FEV1%, which implies poorer pulmonary function ( < 0.001). Besides, GOLD, BODE, and SGRQ were significantly higher in the high PaCO2 group, which are described in Table 2.
3.2. Correlation of PaCO2 with COPD Outcomes
To explore the correlation between PaCO2 and scales of COPD outcomes, we used Pearson’s or Spearman’s correlation analysis. Table 3 manifested the results that age, BMI, pH, and pH on day 7 of admission were negatively related to PaCO2. As the PaCO2 increased, the indicators of respiratory function such as FEV1, FVC, FEV1/FVC, and FEV1% dramatically descended. Figure 2 vividly illustrates the decreasing tendency of FEV1% with the PaCO2 elevated (r = −0.612, < 0.001). Conversely, the scores of GOLD, BODE, and SGRQ grew up markedly with the increased PaCO2. No other statistical correlations were found between the rest of the parameters and PaCO2.

On the basis of the study results above, we observed the intimate link of PaCO2 and GOLD or BODE, so we plotted the ROC curve analysis to evaluate the utility of PaCO2 for predicting poorer pulmonary function (GOLD ≥3) and higher disease severity (BODE ≥5) of COPD. The results are shown in Figure 3 that at a cut-off value of 0.399,the sensitivity and specificity of PaCO2 in predicting BODE ≥5 were 73.9% and 66.0% , respectively, with an AUC of 0.748 (95% CI: 0.630–0.866, = 0.001). At a cut-off value of 0.399, the sensitivity and specificity of PaCO2 for GOLD ≥3 were 44.2% sensitivity and 95.7% specificity, respectively, with an AUC of 0.755 (95% CI: 0.643–0.867, < 0.001). It elucidated that PaCO2 had high predictive value for severe COPD (GOLD ≥3 and BODE ≥5).

(a)

(b)
3.3. Increased PaCO2 Level Is Related to Higher COPD Severity
In terms of the low PaCO2 group, 17.14% suffered from severe COPD (BODE ≥5), while the proportion rose up to 44.74% in the high PaCO2 group ( = 0.030, Figure 1(d)). Additionally, participants with GOLD ≥3 took up a proportion of 80.00% in the high PaCO2 group in contrast with 57.14% in the low PaCO2 group ( = 0.002, Figure 1(e)). Moreover, in comparison with the low PaCO2 group, the high PaCO2 group comprised a larger proportion of patients with SGRQ ≥ 25, which represents poorer life quality ( = 0.030, Figure 1(f)).
Afterwards, we used univariate logistic regression analyses to determine the independent risk factors for COPD severity according to BODE ≥5. Consequently, BMI ( = 0.001), pH ( = 0.002), FEV1/FVC ( = 0.001), FEV1% ( < 0.001), GOLD ≥3 ( = 0.015), mMRC ≥2 ( < 0.001), and PaCO2 ( = 0.003) were observed to have significant relevance with BODE ≥5 in Table 4. Furthermore, the multivariate logistic regression analyses were performed with potential confounders controlled. In Model 1, nothing was adjusted (odds ratio (OR) = 1.146, 95% CI: 1.047–1.253, = 0.003). In Model 2, adjustment for age, sex, smoking status, and duration of disease slightly improved the magnitude of the OR (odds ratio (OR) = 1.160, 95% CI: 1.049–1.284, = 0.004). Even when adjusted for age, sex, smoking status, duration of disease, and FEV1/FVC (Model 3), the relationship between PaCO2 and BODE ≥5 remained still significant (odds ratio (OR) = 1.160, 95% CI: 1.025–1.313, = 0.019). As can been seen from Table 5, we could draw a conclusion that PaCO2 is an independent risk factor for BODE ≥5.
4. Discussion
To our knowledge, this study is the first to demonstrate the association of PaCO2 with BODE and GOLD index and indicators of pulmonary function in stable COPD subjects. In this study, we found a higher level of PaCO2 in COPD patients with BODE ≥5 than BODE <5. Besides, lower pH, higher level of BODE index, and higher level of GOLD index could be observed in patients with higher level of PaCO2. Furthermore, our data suggested that elevated PaCO2 was associated with other outcome indicators of COPD such as pulmonary indexes, especially lower levels of FEV1 and FEV1% and higher SGRQ index (SGRQ ≥25), which imply poorer pulmonary function and health-related quality of life in stable COPD, respectively [13, 17]. Univariate and multivariate logistic regression analyses illustrated that increased PaCO2 was an independent marker of BODE ≥5 and associated with reduced outcomes in stable COPD. We conducted the ROC curves to estimate the predictive value of PaCO2 for severe COPD (defined by BODE ≥5 and GOLD ≥3). The results showed the high accuracy of PaCO2 for predicting severe COPD, which may provide some clues to clinic.
Severe COPD patients with GOLD III and IV will develop CHRF [18], which provokes more often exacerbations, shown to be a poor prognostic factor for survival [5, 14, 19]. Ventilation-perfusion inequality always causes an increased arterial PaCO2 or hypercapnia in COPD patients [20]. Murphy et al. observed that hypercapnia decreases the secretion of IL-6 and tumor necrosis factor in the lungs and impairs lung neutrophil function in an animal model of lung infection [2]. Hypercapnic patients had lessened ambulation distance and worse function capacity, while the weaker respiratory muscles may be to blame [14].
A previous study showed that prescription of long-term oxygen therapy (LTOT) in patients with COPD and reversible hypercapnia slowed down the natural decline in exercise performance and improved exertional dyspnea [21]. Short-term NIV has become an accepted management approach for patients with CHRF [2]. Moreover, noninvasive positive pressure ventilation in patients with stable hypercapnic and high inspiratory pressures targeted to reduce partial pressure of CO2 in arterial blood to some extent can improve survival [22]. Clinical trials have reported the effects of NIV in terms of blood gases and pulmonary function status including FVC, FEV1, and quality of life in COPD [2]. The improvement of FVC reflects an amelioration of lung hyperinflation [2]. Our findings were consistent with the literature, which implicated the association of PaCO2 and poor outcomes of COPD. Zikyri et al. [5] reported that domiciliary NIV reduced the levels of PaCO2 of stable COPD patients with CHRF with a significant improvement in the BODE index. NIV can not only alleviate the airway wall edema to improve the lung function but also relieve fatiguing respiratory muscles by reducing degree of hyperinflation to provide patients a better exercise capacity [2]. Multiple mechanisms contain the improvement of sleep hypoventilation and better ventilation-perfusion (V/Q) matching could attribute to the amelioration of exacerbation number in COPD subjects, especially with hypoxemia during sleep [5]. Moreover, COPD patients with hypercapnia, despite severe ventilatory impairment and weak respiratory muscles, tolerate exercise well and benefit significantly from intensive inpatient pulmonary rehabilitation [14]. Furthermore, vv-ECCO2R has been shown to correct severe respiratory acidosis rapidly and decrease elevated systolic pulmonary artery pressure [6]. Besides, it was found previously that there was a strong correlation between pulmonary artery dilatation and severe COPD exacerbation [6]. This clinical trial Karagiannidis et al. [6] conducted could expound the close connection of PaCO2 and severity of COPD more directly.
However, there are still several limitations in our study. First, this study is a cross-sectional survey limited to retrospective cohort; thus, further research of PaCO2 associated with longitudinal outcomes needs to be conducted to verify the causality. Second, the sample size is relatively small. Third, the inherent mechanisms of the impact of PaCO2 on COPD remain to be clarified.
In conclusion, as a simple, reproducible, widely available, and inexpensive predictive tool, PaCO2 has a homogeneous and strong effect on severity of stable COPD. Patients with higher levels of PaCO2 were associated with poorer outcomes in terms of BODE or GOLD index. Moreover, PaCO2 has predictive value for severe COPD (BODE ≥5 or GOLD ≥3). Our preliminary findings seem worthy of attention because PaCO2 as a biomarker of COPD severity can be obtained from arterial blood gas analysis fast and conveniently in clinic.
Abbreviations
COPD: | Chronic obstructive pulmonary disease |
PaCO2: | Partial pressure of carbon dioxide in arterial blood |
PaO2: | Partial pressure of oxygen in arterial blood |
FEV1%: | Forced expiratory volume in 1 second in percent of the predicted value |
FEV1: | Forced expiratory volume in 1 second |
FVC: | Forced vital capacity |
BODE: | Body mass index, airway obstruction, dyspnea, and exercise capacity |
GOLD: | Global Initiative for Chronic Obstructive Lung Disease |
SGRQ: | St. George’s Respiratory Questionnaire |
mMRC: | Modified Medical Research Council dyspnea score |
CHRF: | Chronic hypercapnic respiratory failure |
NIV: | Noninvasive ventilation |
AECOPD: | Acute exacerbation of chronic obstructive pulmonary disease |
LTOT: | Long-term oxygen therapy |
vv-ECCO2R: | Veno-venous extracorporeal CO2 removal |
BMI: | Body mass index |
6MWD: | 6-Minute walking distance |
ROC: | Receiver operating characteristics |
AUC: | The area under the curve. |
Data Availability
The datasets generated and/or analyzed during this study are not publicly available due to individual privacy but are available from the corresponding author on reasonable request.
Disclosure
Xiaodiao Zhang, Xiaqi Miao, and Keke Ding are co-first authors.
Conflicts of Interest
The authors report no conflicts of interest.
Authors’ Contributions
All authors made substantial contributions to conception and design, acquisition of data, or analysis and interpretation of data; took part in drafting the article or revising it critically for important intellectual content; agreed to submit to the current journal; gave final approval of the version to be published; and agreed to be accountable for all aspects of the work. Xiaodiao Zhang, Xiaqi Miao and Keke Ding contributed equally to this work.
Acknowledgments
This study was supported by the Wenzhou Municipal Sci-Tech Bureau Program (no. Y20210842). The authors thank all the participants in the study and the staff at the Third Affiliated Hospital of Wenzhou Medical University for their contribution in obtaining the data and assisting in the successful completion of the COPD study.