Your new experience awaits. Try the new design now and help us make it even better

SYSTEMATIC REVIEW article

Front. Med., 22 September 2025

Sec. Pulmonary Medicine

Volume 12 - 2025 | https://doi.org/10.3389/fmed.2025.1640084

A meta-analysis of the effects of long-term oxygen therapy combined with exercise rehabilitation on exercise capacity, cardiopulmonary function, and quality of life in patients with COPD


Qi Muge,Qi Muge1,2SurigugaSuriguga1YuqingYuqing1AronggaowaAronggaowa3Taojin
Taojin1*Lanying Chen
Lanying Chen4*
  • 1Department of Geriatrics and Respiratory Medicine, Affiliated Hospital of Inner Mongolia Minzu University, Tongliao, China
  • 2School of Pharmacy, Jiangxi University of Traditional Chinese Medicine, Nanchang, Jiangxi, China
  • 3Clinical College, Inner Mongolia Minzu University, Tongliao, China
  • 4National Pharmaceutical Engineering Center for Solid Preparation of Chinese Herbal Medicine, Jiangxi University of Traditional Chinese Medicine, Nanchang, Jiangxi, China

Objective: To assess the effectiveness of long-term oxygen therapy (LTOT) combined with exercise rehabilitation vs. exercise rehabilitation alone in improving exercise capacity, cardiopulmonary function, and quality of life in chronic obstructive pulmonary disease (COPD) patients.

Methods: A comprehensive literature search was conducted in the Chinese Biomedical Literature Database (CBM), Wanfang, China National Knowledge Infrastructure (CNKI), Cochrane Library, EMBASE, ScienceDirect, and PubMed for studies published from January 2010 to the present. Controlled clinical trials comparing oxygen therapy and/or exercise rehabilitation in COPD patients were included. Two independent reviewers extracted data and assessed risk of bias using the Cochrane Handbook (version 5.3). Meta-analysis was performed using RevMan 5.3.

Results: Nine studies (N = 703) met inclusion criteria. Compared with the control group (CG), the combined LTOT and exercise group showed significant improvements in 6-min walk distance (6MWD), forced expiratory volume in 1 s (FEV1), and FEV1/FVC ratio (P < 0.05). PaO2 levels tended to be higher but showed substantial heterogeneity. No significant differences were observed in blood oxygen saturation, heart rate, or PaCO2. Quality of life significantly improved in the combined therapy group.

Conclusion: LTOT combined with exercise rehabilitation is more effective than exercise alone in improving exercise capacity, pulmonary function, and quality of life in COPD patients. However, cardiac benefits remain unclear, warranting further studies.

1 Introduction

Chronic obstructive pulmonary disease (COPD) is a primarily characterized by symptoms such as excessive sputum production, chronic cough, and dyspnea (1). As the disease progresses, patients experience a gradual decline in physical activity tolerance and cardiopulmonary function, leading to significant challenges such as difficulty breathing during daily activities. These deteriorations can severely impact individuals' quality of life and pose life-threatening risks (2).

Long-term oxygen therapy (LTOT) is recommended for severe resting hypoxemia and improves survival; it does not routinely benefit moderate resting or isolated exertional desaturation (3). LTOT is particularly indicated for patients with severe COPD, especially those who exhibit hypoxemia even at rest (4). Specifically, according to the Global Initiative for Chronic Obstructive Lung Disease (GOLD) guidelines, LTOT is recommended for patients with resting arterial oxygen saturation (SpO2) ≤ 88% or partial pressure of arterial oxygen (PaO2) ≤ 55 mm Hg (5). Additionally, LTOT may be considered for those with PaO2 ≤ 60 mm Hg if they have signs of cor pulmonale, polycythemia (hematocrit ≥55%), or pulmonary hypertension. Although LTOT has shown significant effects in improving blood oxygen levels, its direct effects on exercise tolerance and daily function are limited. Therefore, oxygen therapy alone may be insufficient to address the multifaceted rehabilitation needs of COPD patients. At the same time, exercise-based rehabilitation is increasingly vital in COPD management, combining aerobic exercise, strength training, and breathing exercises to enhance muscle strength, endurance, cardiopulmonary function, and overall quality of life (6). However, the heterogeneous nature of COPD patient conditions means that a single rehabilitation model often fails to address the unique needs of all individuals. Consequently, personalized or multifaceted rehabilitation approaches may be necessary to effectively meet the diverse requirements of COPD patients.

The combined use of LTOT and exercise rehabilitation has attracted increasing attention as a synergistic strategy to enhance physical capacity, slow disease progression, and reduce the healthcare burden in COPD. Yet, variability in exercise protocols, differing in type, intensity, frequency, and duration, has led to inconsistent outcomes, posing challenges for clinical practice. This study therefore conducts a systematic meta-analysis to evaluate the effects of LTOT combined with exercise rehabilitation on exercise capacity, cardiopulmonary function, and quality of life in stable COPD, with the goal of informing more standardized and effective rehabilitation strategies.

2 Research content and methods

2.1 Sources and search methods for literature resources

The following electronic databases were used to do a thorough literature search: the Chinese Biomedical Literature Database (CBM), China National Knowledge Infrastructure (CNKI), ScienceDirect, EMBASE, Cochrane Library, VIP Full-Text Database, Wanfang Database, and PubMed. The search encompassed relevant Chinese and international journals, dissertations, conference papers, news articles, and manual retrieval sources, and was further supplemented by literature tracing methods. Studies eligible for inclusion were clinical control trials that evaluated the effects of oxygen therapy and exercise rehabilitation in patients with COPD.

The search strategy employed both free-text keywords and subject-specific terms, including but not limited to “oxygen therapy,” “exercise rehabilitation,” “COPD,” “exercise capacity,” “cardiopulmonary function,” and “quality of life.” The search period spanned from January 2010 to the present. The quality and relevance of the data gathered for the meta-analysis were guaranteed by eliminating duplicate research and choosing pertinent publications according to predetermined inclusion and exclusion criteria.

2.2 Inclusion criteria for literature

(1) Study type: randomized controlled trials (RCTs) and controlled cohort studies that investigated oxygen therapy and/or exercise rehabilitation as interventions for COPD patients.

(2) Study subjects: patients diagnosed with COPD; diagnostic criteria as described in relevant literature (7).

(3) Intervention measures: the study group (SG) received LTOT combined with exercise rehabilitation, while the control group (CG) received LTOT alone or exercise rehabilitation alone. Eligible rehabilitation programs included aerobic exercise, strength training, breathing training, or combinations thereof, with clearly defined duration.

(4) Reported one or more of the following outcome indicators: (1) 6-Min Walk Distance (6MWD); (2) blood gas indicators; (3) heart rate; (4) pulmonary function indicators; (5) quality of life.

2.3 Literature exclusion criteria

These criteria were used to determine which studies were not included in the meta-analysis:

(1) Non-controlled studies: studies without a parallel comparator group (e.g., single-arm designs).

(2) Incomplete data reporting: the reported data were incomplete or insufficient for analysis.

(3) Duplicate publications: studies that duplicated research content; in such cases, only the most recent study was included.

(4) Insignificant therapeutic effect evaluation: studies that did not include evaluable therapeutic outcomes relevant to this analysis.

(5) Review articles: publications that were review articles or related literature rather than original research.

(6) Case reports/case series: case reports and small case series without a parallel comparator.

2.4 Quality evaluation and data extraction

(1) Bias risk assessment of included studies: the evaluation was conducted using the “bias risk Assessment” technique suggested in the Cochrane Manual of Systematic Reviews 5.3. (2) Literature screening and data extraction: to maintain impartiality and reduce bias, two researchers independently carried out data extraction, literature screening, and quality rating. A third researcher was contacted to make the final judgment if the researchers were unable to achieve a consensus after discussing and resolving any disagreements. Microsoft Excel and Note Express literature management software made it easier to organize and extract research resources. For studies with incomplete data, the corresponding authors were contacted in order to acquire the necessary information. (3) The data extraction process encompassed the following elements: (1) basic information: author(s), year of publication, and number of participants; (2) intervention measures: details of the intervention plan and treatment course; (3) outcome indicators: 6MWD, blood gas indicators, heart rate, pulmonary function indicators, and quality of life, etc.

2.5 Statistical analysis

Using RevMan 5.3 software, the meta-analysis was completed in accordance with the Cochrane Collaboration's guidelines. For each outcome measure, including exercise endurance, blood gas indices, heart rate, lung function indices, and quality of life, the mean values and standard deviations for both the SG and the CG were extracted and input into RevMan 5.3. For outcomes measured using the same scale across all studies (e.g., 6MWD, PaO2, PaCO2, FEV1, FEV1/FVC), Mean difference (MD) with 95% confidence intervals (CIs) was calculated. For outcomes assessed using different measurement scales across studies (e.g., quality of life), Standardized Mean Difference (SMD) with 95% CIs was used. First, to ascertain if there was statistical heterogeneity, the Chi2 test was employed. Studies were considered homogenous if the P-value for the Chi2 test was more than 0.05 and the I2 statistic was less than 50%. In such cases, for the meta-analysis, a fixed-effect model was used. If the P-value from the Chi2 test was less than 0.05 and the I2 statistic was equal to or greater than 50%, demonstrating considerable heterogeneity, a random effects model (REM) was selected for the analysis. In instances where significant heterogeneity was detected (P < 0.05) and the sources of heterogeneity couldn't be identified; a meta-analysis was deemed inappropriate. Instead, a descriptive analysis was performed to summarize the findings qualitatively. Subgroup differences were examined using the Chi2 test for interaction. To account for clinical/intervention heterogeneity, subgroup analyses were prespecified. For exercise capacity (6MWD), studies were categorized by exercise regimen (mixed-modality vs. aerobic-only). For other outcomes (e.g., blood gases, pulmonary function, QoL), subgroup analyses were planned by exercise type, LTOT duration, or baseline hypoxemia severity (severe vs. non-severe/unclear). However, the analyses were not feasible due to the small number of studies and inconsistent reporting of intervention characteristics. To further examine the publication bias of the included literature, an inverted funnel plot was created. The asymmetry of the funnel plot was examined using Eggers' test. In cases where the P-value of this test was less than 0.1, the Trim and Fill method may be utilized to adjust the impact size of the potential publication bias and correct the funnel plot.

3 Results and analysis

3.1 The findings of literature search and the basic information of the included literature

The Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) standards were used for conducting the literature search. Initially, a total of 1,309 records were retrieved from computer databases. After removing duplicate studies, 934 unique records remained. Of these, 715 unnecessary studies were eliminated after these data were filtered based on their abstracts and titles. Subsequently, 423 studies were preliminarily included after excluding case reports, reviews, and uncontrolled studies.

The full texts of these 423 studies were meticulously reviewed to identify those with complete data and relevant primary outcome indicators. During this stage, 414 studies were excluded due to the absence of primary outcome measures. Nine clinical controlled trials, totaling 703 participants, were ultimately included in the meta-analysis after meeting all inclusion criteria. Figure 1 shows the literature screening flow chart. Table 1 displays the fundamental attributes of the included literature.

Figure 1
Flowchart illustrating the identification of studies via databases. Initially, 1309 records are identified, with 375 removed before screening. After screening 934 records, 219 are excluded. From 715 reports sought for retrieval, 292 are not retrieved. Out of 423 assessed for eligibility, 414 are excluded due to incomplete data, failure to meet criteria, or unclear indicators. Ultimately, 9 studies are included.

Figure 1. Literature screening flow chart.

Table 1
www.frontiersin.org

Table 1. Basic characteristics of literature.

3.2 Evaluation of methodological quality of included literature

The methodological quality of the nine included studies was assessed using the Cochrane Risk of Bias Tool (version 5.3) across seven domains. Random sequence generation was low risk in 8/9 (88.9%) studies and unclear in 1/9 (11.1%). Allocation concealment was explicitly reported (low risk) in 2/9 (22.2%), unclear in 3/9 (33.3%), and high risk in 4/9 (44.4%) due to insufficient methodological detail. Blinding of participants/personnel (performance bias) was rated high in 3/9 (33.3%) and unclear in 6/9 (66.7%), reflecting the practical difficulty of blinding exercise/oxygen interventions and incomplete reporting. Blinding of outcome assessment (detection bias) was low risk in 4/9 (44.4%), unclear in 3/9 (33.3%), and high in 2/9 (22.2%). Attrition bias was low in 8/9 (88.9%) and unclear in 1/9 (11.1%). Selective reporting was low risk in 7/9 (77.8%), unclear in 1/9 (11.1%), and high in 1/9 (11.1%). Other bias was rated low risk in 9/9 (100%). These distributions correspond to the domain-level plots in Figures 2, 3.

Figure 2
Bar chart depicting risk of bias across seven categories: random sequence generation, allocation concealment, blinding of participants and personnel, blinding of outcome assessment, incomplete outcome data, selective reporting, and other bias. Bars show proportions of low (green), unclear (yellow), and high (red) risk of bias. Most categories have a majority low risk of bias.

Figure 2. Risk bias diagram.

Figure 3
Risk of bias summary table for various studies from 2014 to 2024, displaying evaluation across seven categories: random sequence generation, allocation concealment, blinding of participants and personnel, blinding of outcome assessment, incomplete outcome data, selective reporting, and other bias. Green plus signs indicate low risk, yellow question marks indicate unclear risk, and red minus signs indicate high risk.

Figure 3. Summary of risk bias.

3.3 Exercise endurance

The heterogeneity test showed χ2 = 6.56, df = 4, P = 0.16, I2 = 39%, indicating no significant heterogeneity. The fixed-effect model analysis (Figure 4) indicated that 6MWD in the SG was significantly longer than that in the CG after intervention (P < 0.05). When stratified by rehabilitation protocol, studies incorporating mixed exercise components (e.g., aerobic exercise, breathing training, or combined modalities) delivered alongside nasal catheter oxygen inhalation showed a mean improvement of 55.84 m (95% CI: 45.32–66.35) in 6MWD compared with the CG, while studies combining aerobic exercise with oxygen therapy showed no significant difference [−7.35 m (95% CI: −70.58 to 55.88)]. Across both approaches, the SG demonstrated longer 6MWD than the CG after intervention (P < 0.05; Figure 5).

Figure 4
Forest plot comparing experimental and control groups across five studies. Each study shows mean, standard deviation, total participants, and weight. Mean differences with confidence intervals are plotted, indicating a combined mean difference of 54.14 favoring the experimental group. Heterogeneity chi-squared is 6.56 with a P-value of 0.16, and overall effect Z-score is 10.23 with a P-value less than 0.00001.

Figure 4. Forest analysis diagram of 6MWD comparison between the two groups after intervention. “Experimental” refers to the Study Group (SG) receiving LTOT + exercise rehabilitation; “Control” refers to the Control Group (CG) receiving either intervention alone.

Figure 5
Forest plot showing the comparison of experimental and control groups in studies of nasal catheter oxygen inhalation and aerobic exercise. The plot includes a summary of mean differences with 95% confidence intervals. The studies under nasal catheter oxygen indicate a positive mean difference favoring the experimental group. Aerobic exercise shows varying results, with a subgroup not significantly favoring either group. Overall, the total mean difference indicates a slight favor towards the experimental interventions with heterogeneity details included.

Figure 5. Forest plot analyzing the impact of different oxygen therapy modalities on the 6-Min Walk Distance (6MWD) in two groups of COPD patients. Experimental = SG; Control = CG.

3.4 Blood oxygen saturation

A total of 102 participants from two studies reported SpO2. Heterogeneity was moderate (χ2 = 2.29, df = 1, P = 0.13; I2 = 56%). The REM analysis (Figure 6) showed no significant difference between groups (MD = 1.25%, 95% CI: −0.35 to 2.85; P = 0.13).

Figure 6
Forest plot comparing experimental and control groups from two studies, Marrara KT 2018 and Xu YL 2014, examining mean differences with values and 95% confidence intervals. The overall mean difference is 0.70 with a confidence interval of [-1.78, 3.18]. The plot reflects heterogeneity with I² = 56% and overall effect test, Z = 0.55 with P = 0.58.

Figure 6. Comparison of blood oxygen saturation between the two groups after intervention. Experimental = SG; Control = CG.

3.5 Arterial partial pressure of oxygen (PaO2)

Five studies including 461 participants evaluated PaO2. Heterogeneity was substantial (χ2 = 42.58, df = 4, P < 0.00001; I2 = 91%). The REM analysis (Figure 7) indicated that PaO2 was significantly higher in the SG (MD = 5.42 mmHg, 95% CI: 1.85–9.00; P = 0.003).

Figure 7
Forest plot showing meta-analysis results of five studies comparing experimental and control groups. Mean differences are displayed with 95% confidence intervals. Overall effect favors the experimental group with a mean difference of 8.16 [6.29, 10.03]. Significant heterogeneity is present (I² = 91%).

Figure 7. Forest analysis diagram of PaO2 comparison involving the two groupings after intervention. Experimental = SG; Control = CG.

3.6 Arterial partial pressure of carbon dioxide (PaCO2)

A meta-analysis was conducted on PaCO2 in the two groups after intervention, and a total of 461 samples were included into five literatures. The heterogeneity test results showed χ2 = 1,818.03, df =4, P < 0.00001, I2 = 100%, indicating that there was significant heterogeneity among the included study data. The REM analysis (Figure 8) showed no significant difference (MD = −0.86 mmHg, 95% CI: −3.20 to 1.48; P = 0.47).

Figure 8
Forest plot showing the mean differences between experimental and control groups from five studies. The mean differences with 95% confidence intervals are plotted: negative values favor control and positive values favor experimental. Studies include He YM 2017, Tian RY 2019, Wang XZ 2024, Wu SG 2020, Yan WP 2022. The overall effect size is 4.80 with a confidence interval of negative 7.42 to 17.02, indicating no significant overall effect (p-value = 0.44). Heterogeneity is high (Tau-squared = 193.74).

Figure 8. Forest analysis diagram of PaCO2 comparison involving the two groupings after intervention. Experimental = SG; Control = CG.

3.7 Heart rate

Two studies including 102 participants assessed heart rate. Heterogeneity was negligible (χ2 = 0.32, df = 1, P = 0.57; I2 = 0%), showing that the included study data showed no discernible heterogeneity. The fixed-effect model (Figure 9) showed no significant difference (MD = −1.15 beats/min, 95% CI: −4.32 to 2.02; P = 0.48).

Figure 9
Forest plot comparing experimental and control groups from two studies: Marrara KT 2018 and Xu YL 2014. The plot shows the mean differences with 95% confidence intervals for each study and combined. The overall mean difference is -0.51, with a confidence interval from -4.58 to 3.56. The heterogeneity is low with Chi² = 0.32 and I² = 0%. The test for overall effect shows Z = 0.24.

Figure 9. Forest analysis diagram of heart rate comparison between the two groups after intervention. Experimental = SG; Control = CG.

3.8 FEV1 forced expiratory volume in 1 s (FEV1)

A meta-analysis of FVE1 after intervention was conducted in the two groups, including seven articles with a total of 656 samples. The heterogeneity test results showed χ2 = 116.87, df = 6, P ≤ 0.00001, I2 = 95%, indicating considerable heterogeneity among the included research data. The REM analysis (Figure 10) showed that FEV1 was significantly higher in the SG (MD = 0.23 L, 95% CI: 0.08 to 0.38; P = 0.002).

Figure 10
Forest plot displaying mean differences between experimental and control groups in seven studies. Each study shows mean, standard deviation, total, and weight. Confidence intervals are visualized with lines and squares, indicating variability and effect size. Overall effect shows a mean difference of 0.27 with a confidence interval of 0.05 to 0.48, favoring the experimental group. Heterogeneity statistics are included.

Figure 10. Forest analysis diagram for comparison of FEV1 values between the two groups after intervention. Experimental = SG; Control = CG.

3.9 FEV1/FVC [FEV1/forced vital capacity ratio (FVC)]

Seven studies with 656 participants assessed FEV1/FVC. Heterogeneity was substantial (χ2 = 44.36, df = 6, P < 0.00001; I2 = 86%). The REM analysis (Figure 11) showed that the SG had significantly higher FEV1/FVC values (MD = 3.86%, 95% CI: 1.31–6.40; P = 0.003).

Figure 11
Forest plot displaying the mean differences and 95% confidence intervals for seven studies comparing experimental and control groups. Each study shows mean, standard deviation, and total for both groups. The overall mean difference is 3.86 with a confidence interval of 1.31 to 6.40, favoring the experimental group. Heterogeneity tests: Tau² = 9.52; Chi² = 44.36, degrees of freedom = 6, with a p-value < 0.00001; I² = 86%. The test for overall effect shows Z = 2.97, p = 0.003.

Figure 11. Forest analysis diagram comparing FEV1/FVC values between the two groups after intervention. Experimental = SG; Control = CG.

3.10 Quality of life

A total of 165 patients from two trials [Jianmei et al. (8); Yanli et al. (9)] reported quality-of-life outcomes across four subdomains (breathing, tiredness, emotion, and disease control). The heterogeneity test showed χ2 = 9.80, df = 7, P = 0.20, I2 = 29%, indicating no significant heterogeneity. Pooled analysis using a fixed-effect model (Figure 12) demonstrated a SMD of 0.90 (95% CI: 0.74–1.06; P < 0.001), favoring the SG over the CG. This reflects a large effect size and suggests a clinically meaningful improvement in quality of life.

Figure 12
Forest plot displaying standardized mean differences with 95% confidence intervals across four subgroups: Breather, Tired, Emotion, and Disease Control. The summary effect sizes and heterogeneity tests are provided for each subgroup. All subgroups favor the experimental group over the control, as indicated by the plot's alignment and effect size values. Overall, the total effect also favors the experimental group, with a notable statistical significance.

Figure 12. Forest analysis of the comparison of quality-of-life scores involving the two groupings after intervention. Experimental = SG; Control = CG.

3.11 Publication bias analysis

Funnel plots were prepared for each outcome (Figures 1315; Supplementary Figures S1, S2). Across outcomes, the number of independent studies (k_studies) ranged from 2 to 7, and no outcome reached k_studies ≥10. For QoL, two trials reported CRQ domain scores; the funnel shows 7 domain-level estimates from these two studies, so k_studies = 2 and estimates = 7. In accordance with methodological guidance, Egger's regression was not performed for any outcome. Funnel plots are therefore descriptive only, and any apparent asymmetry should be interpreted cautiously given small k_studies and between-study heterogeneity.

Figure 13
Funnel plot displaying standard error (SE) of mean difference (MD) against MD values ranging from negative five hundred to five hundred. The plot has a symmetrical triangular shape formed by dashed lines, with circles indicating data points distributed around the vertical line at zero MD.

Figure 13. Funnel plot based on 6WMD.

Figure 14
Three scatter plots labeled a, b, and c, each showing different distributions of data points along x and y axes marked as MD and CE, respectively. Vertical lines are present in each graph.

Figure 14. Funnel plot based on blood gas index. (a) Blood oxygen saturation; (b) PaO2; (c) PaCO2.

Figure 15
Funnel plot displaying the standard error (SE) of mean difference (MD) on the vertical axis against the MD on the horizontal axis. Two open circles are plotted within a symmetrical inverted funnel shape formed by dashed lines centered around zero on the horizontal axis.

Figure 15. Funnel plot based on heart rate.

4 Analysis and discussion

COPD is a common respiratory disease in clinical practice, with the characteristics of prolonged and difficult to cure, easy to repeat. According to relevant epidemiological study (10), the prevalence rate of COPD in adults over 20 years old in China is 8.6%, and that in adults over 40 years old can reach 13.7%, which has become one of the main disease factors threatening the life and health of middle-aged and elderly people in China. The course of COPD has two stages: acute exacerbation and stable stage. Compared with the acute stage, the severity of the condition of patients in stable stage is significantly reduced, but its potential risks cannot be ignored (11). Patients with stable COPD are susceptible to several physiological impairments, including reduced airway compliance, increased airway resistance, compromised ventilatory function, and diminished exercise endurance. These problems may result in hypoxemia and, in more extreme situations, respiratory failure, thereby posing significant risks to patients' lives and overall health (12, 13). Consequently, it is essential to implement scientifically grounded treatments and appropriate interventions at all stages of COPD. Such strategies aim to control the progression of the disease, prevent exacerbations, and improve the quality of life of patients (14, 15).

The improvement in exercise endurance is a key marker of rehabilitation efficacy in COPD (16). The patients often have problems such as limited respiratory function, decreased muscle strength, and chronic fatigue, which result in their exercise endurance being significantly lower than that of healthy people (17, 18). In our analysis, the 6MWD in the SG was consistently longer than in the CG by a pooled mean difference of 54.14 m (95% CI 43.77–64.51; P < 0.001; I2 = 39.1%; k = 5), a magnitude consistent with clinically meaningful improvement (P < 0.05). The gain was driven primarily by programs that delivered nasal-cannula oxygen during mixed-modality training (aerobic with breathing and/or strength work), which improved 6MWD by 55.84 m (95% CI 45.32–66.35); in contrast, aerobic-only exercise combined with oxygen did not show a clear advantage (MD −7.35 m, 95% CI −70.58 to 55.88). Both approaches were associated with improvement in walking distance, suggesting the beneficial impact of LTOT-supported rehabilitation across different structured exercise regimens. This finding may be attributed to the fact that exercise rehabilitation exercises enhance the muscle strength and peripheral blood flow in COPD patients, and improve the oxygen utilization efficiency of the whole body. In addition, the continuous oxygen supply of oxygen therapy supports the oxygen demand during exercise, thereby reducing the sense of dyspnea generated during exercise (19, 20).

COPD is primarily characterized by restricted airflow and impaired pulmonary ventilation. Consequently, improvements in pulmonary function indicators directly reflect the recovery of a patient's respiratory system (21). Commonly utilized pulmonary function metrics include Forced Vital Capacity (FVC), Forced Expiratory Volume in One Second (FEV1), and the FEV1/FVC ratio. These indicators serve not only as diagnostic and grading criteria for COPD but also play a vital part in evaluating the efficacy of rehabilitation interventions. Although the meta-analysis demonstrated a statistically significant improvement in the FEV1/FVC ratio following LTOT combined with exercise rehabilitation (WMD = 3.86%, 95% CI: 1.31–6.40; P = 0.003), this change lies near the lower threshold of clinical relevance. The FEV1/FVC ratio is widely used for diagnosis but is less sensitive for detecting short-term therapeutic effects in COPD. Current guidelines and expert consensus do not define a minimal clinically important difference (MCID) for this parameter. Prior studies suggest that changes smaller than 3% are unlikely to be clinically meaningful, whereas increases of 4%−5% or more may reflect relevant physiological improvement, particularly when accompanied by gains in FEV1 or quality of life (2224). Our observed effect approaches this threshold, suggesting possible clinical importance, but should be interpreted alongside other outcome measures. Moreover, significant improvements in pulmonary function were observed following the combined intervention of LTOT and exercise rehabilitation. This enhancement is likely attributable to the synergistic effects of these treatments in promoting lung function recovery. Specifically, LTOT improves tissue oxygenation by providing sufficient oxygen, thereby alleviating breathing difficulties and reducing the cardiorespiratory burden caused by hypoxia (25, 26). Concurrently, exercise rehabilitation enhances the strength and endurance of respiratory muscles, which improves lung ventilation efficiency and further facilitates the recovery of pulmonary function (27).

For blood gases, PaO2 was higher in the study group, but with substantial between-study heterogeneity (χ2 = 42.58, df = 4; I2 = 91%; P < 0.00001; k = 5), whereas PaCO2 showed no statistically significant between-group difference with very high heterogeneity (I2 = 100%; k = 5). SpO2 likewise did not differ significantly (k = 2). These very high I2 values indicate that true effects vary across studies, so a single average effect should be interpreted cautiously. Plausible contributors to variability include differences in oxygen modality, oxygen titration and flow targets, program duration, and baseline hypoxemia severity across trials. In addition, structural lung damage in COPD (e.g., airway stenosis and alveolar destruction) limits gas exchange capacity, making large and uniform improvements in PaO2 or PaCO2 less likely in the short term (2830). PaCO2 stability is also influenced by systemic factors such as metabolic status and respiratory drive (31, 32). These factors likely underlie the null effects and high heterogeneity, emphasizing the need for standardized protocols and larger trials to define the impact of LTOT plus rehabilitation on blood gases. Similarly, quality of life showed a pooled effect size of SMD = 0.90, suggesting a potentially large and clinically relevant improvement. However, this finding is derived from only two studies (n = 165), and should therefore be interpreted with caution. Compared with lung rehabilitation training alone, LTOT combined with exercise rehabilitation was associated with effective outcomes. This combination helps reduce fear and anxiety related to hypoxia, alleviates breathing difficulties, enhances exercise compliance, and ultimately improves the quality of life in COPD patients. Although there are some differences in exercise rehabilitation exercise standards and exercise time among the literatures included in this study, the overall results suggest that LTOT combined with exercise rehabilitation may be associated with improvements in exercise capacity, cardiopulmonary function, and quality of life of COPD findings. These benefits appeared across different program types, including mixed-modality and aerobic-focused protocols, indicating a consistent trend despite variability. However, given the small number of trials, heterogeneity in protocols, and overall moderate methodological quality, these findings should be regarded as tentative, and further high-quality, standardized studies are needed to confirm the robustness and generalizability of the effects.

This review is limited by the moderate methodological quality (grade B) of the included trials and by incomplete reporting in key domains. In particular, allocation concealment was adequately reported in 2/9 studies (22%), and blinded outcome assessment in 4/9 (44%), while performance bias was high or unclear in all studies given the practical challenges of blinding exercise/oxygen interventions; these features introduce potential selection and performance/detection bias that may influence effect estimates. Between-study heterogeneity was substantial for several outcomes, most notably PaCO2 (I2 ≈ 100%), and likely reflects variation in oxygen-delivery approach and titration, baseline hypoxemia, and rehabilitation dose/content (type, intensity, and duration). Exploratory subgroup analysis was possible for 6MWD, which suggested benefits in mixed-modality programs with LTOT, but subgroup analyses for other outcomes (e.g., PaO2, FEV1, QoL) were not feasible due to the small number of studies and inconsistent reporting. Future trials should incorporate pre-specified subgroup analyses to clarify these differences. Quality of life outcomes were available from only two small regional trials [Jianmei et al. (8); Yanli et al. (9)]. These were non-indexed studies with limited methodological reporting, which substantially reduces certainty and generalizability of the pooled QoL estimate. Assessment of small-study effects was also constrained: no outcome met k ≥ 10, so formal asymmetry tests (e.g., Egger) were not performed and funnel plots are descriptive only. Moreover, two trials [Jianmei et al. (8); Yanli et al. (9)] were published in non-indexed regional journals, which may limit external generalizability. Finally, this review was not prospectively registered in PROSPERO or another registry. While this limits transparency, we adhered closely to PRISMA guidelines and predefined our eligibility criteria, search strategy, and analytic approach to minimize bias. Taken together with generally small sample sizes, moderate methodological quality, and reliance on two non-indexed regional trials, these factors reduce certainty around pooled effects. Future research should prioritize multicenter randomized controlled trials with rigorous allocation concealment and blinded outcome assessment, prospective protocol registration (e.g., PROSPERO), and standardized oxygen delivery targets and rehabilitation protocols. In addition, adequately powered studies should specifically evaluate outcomes with high heterogeneity (e.g., PaCO2, FEV1) and confirm the tentative quality-of-life findings, which currently rest on only two small studies. Such work is needed to refine both the magnitude and the generalizability of benefit.

In summary, LTOT combined with exercise rehabilitation was associated with improvements in pulmonary function, exercise capacity, and quality of life in patients with COPD. While these findings are encouraging, they should be interpreted with caution due to the limited number of studies, variability in intervention protocols, and the moderate methodological quality of the included trials. Future large-scale, rigorously designed RCTs are needed to confirm these benefits and define optimal intervention.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

QM: Conceptualization, Data curation, Investigation, Methodology, Project administration, Resources, Software, Supervision, Visualization, Writing – original draft. S: Data curation, Methodology, Supervision, Validation, Writing – original draft. Y: Formal analysis, Investigation, Supervision, Validation, Writing – review & editing. A: Formal analysis, Supervision, Visualization, Writing – review & editing. T: Conceptualization, Data curation, Funding acquisition, Investigation, Project administration, Resources, Software, Visualization, Writing – original draft. LC: Conceptualization, Formal analysis, Funding acquisition, Methodology, Resources, Software, Validation, Visualization, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was sponsored in part by the sub project “Clinical Research on the Treatment of Chronic Obstructive Pulmonary Disease with Mongolian Medicine” under the National R&D Program “Modernization Research of Traditional Chinese Medicine” (Project No. 2018YFC1708002) focuses on the prevention and treatment of major diseases with ethnic minority medicine such as Tibetan medicine, Mongolian medicine, and Uyghur medicine; Medical science and Technology project of Inner Mongolia Health Commission (Grant number 20220164 2) and the Youth Innovation Fund project of Affiliated Hospital of Inner Mongolia Minzu University (Grant number 2020QNJJ01).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2025.1640084/full#supplementary-material

References

1. Christenson SA. Chronic obstructive pulmonary disease. Lancet. (2022) 399:2227–42. doi: 10.1016/S0140-6736(22)00470-6

PubMed Abstract | Crossref Full Text | Google Scholar

2. Sandelowsky H, Weinreich UM, Aarli BB, Sundh J, Hoines K, Stratelis G, et al. COPD – do the right thing. BMC Fam Pract. (2021) 22:244. doi: 10.1186/s12875-021-01583-w

PubMed Abstract | Crossref Full Text | Google Scholar

3. Hess MW. Oxygen therapy in COPD. Respir Care. (2023) 68:998–1012. doi: 10.4187/respcare.10876

PubMed Abstract | Crossref Full Text | Google Scholar

4. Taichman DB, Drazen JM. Making long-term oxygen therapy less burdensome. N Engl J Med. (2024) 391:1050–1. doi: 10.1056/NEJMe2410129

PubMed Abstract | Crossref Full Text | Google Scholar

5. Long term domiciliary oxygen therapy in chronic hypoxic cor pulmonale complicating chronic bronchitis and emphysema. Report of the Medical Research Council Working Party. Lancet. (1981) 1:681–6. doi: 10.1016/S0140-6736(81)91970-X

Crossref Full Text | Google Scholar

6. Kahnert K, Jorres RA, Behr J, Welte T. The diagnosis and treatment of COPD and its comorbidities. Dtsch Arztebl Int. (2023) 120:434–44. doi: 10.3238/arztebl.m2023.0027

Crossref Full Text | Google Scholar

7. Rabe KF, Hurd S, Anzueto A, Barnes PJ, Buist SA, Calverley P, et al. Global Initiative for chronic obstructive lung, global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease: GOLD executive summary. Am J Respir Crit Care Med. (2007) 176:532–55. doi: 10.1164/rccm.200703-456SO

PubMed Abstract | Crossref Full Text | Google Scholar

8. Jianmei F, Jie B, Chuang D, Xinyun K. Effect of joint family oxygen therapy on exercise ability and quality of life of elderly COPD patients in stable stage (in chinese). Health Home. (2021) 11:11–2.

Google Scholar

9. Yanli X. Effect of joint family oxygen therapy on exercise ability and quality of life of elderly COPD patients in stable stage. Southwest Mil Doc. (2014) 16:554–6.

Google Scholar

10. Upadhyay P, Wu CW, Pham A, Zeki AA, Royer CM, Kodavanti UP, et al. Animal models and mechanisms of tobacco smoke-induced chronic obstructive pulmonary disease (COPD). J Toxicol Environ Health B Crit Rev. (2023) 26:275–305. doi: 10.1080/10937404.2023.2208886

PubMed Abstract | Crossref Full Text | Google Scholar

11. Agusti A, Bohm M, Celli B, Criner GJ, Garcia-Alvarez A, Martinez F, et al. 2023: a brief update for practicing cardiologists. Clin Res Cardiol. (2024) 113:195–204. doi: 10.1007/s00392-023-02217-0

PubMed Abstract | Crossref Full Text | Google Scholar

12. Windisch W, Magnet FS. Long-term oxygen therapy in COPD: what is the evidence? Lancet Respir Med. (2022) 10:1010–1. doi: 10.1016/S2213-2600(22)00225-9

PubMed Abstract | Crossref Full Text | Google Scholar

13. Agusti A, Vogelmeier C, Faner R, COPD. 2020: changes and challenges. Am J Physiol Lung Cell Mol Physiol. (2020) 319:L879–83. doi: 10.1152/ajplung.00429.2020

PubMed Abstract | Crossref Full Text | Google Scholar

14. Barnes PJ. COPD 2020: new directions needed. Am J Physiol Lung Cell Mol Physiol. (2020) 319:L884–6. doi: 10.1152/ajplung.00473.2020

PubMed Abstract | Crossref Full Text | Google Scholar

15. Murgia N, Gambelunghe A. Occupational COPD-The most under-recognized occupational lung disease? Respirology. (2022) 27:399–410. doi: 10.1111/resp.14272

PubMed Abstract | Crossref Full Text | Google Scholar

16. O'Donnell DE, Milne KM, James MD, de Torres JP, Neder JA. Dyspnea in COPD: new mechanistic insights and management implications. Adv Ther. (2020) 37:41–60. doi: 10.1007/s12325-019-01128-9

PubMed Abstract | Crossref Full Text | Google Scholar

17. Ferri S, Paoletti G, Pelaia C, Heffler E, Canonica GW, Puggioni F, et al. COPD and biologic treatment: state of the art. Curr Opin Allergy Clin Immunol. (2023) 23:309–18. doi: 10.1097/ACI.0000000000000920

PubMed Abstract | Crossref Full Text | Google Scholar

18. Terry PD, Dhand R. Inhalation therapy for stable COPD: 20 years of GOLD reports. Adv Ther. (2020) 37:1812–28. doi: 10.1007/s12325-020-01289-y

PubMed Abstract | Crossref Full Text | Google Scholar

19. Calverley PMA, Walker PP. Contemporary concise review 2022: chronic obstructive pulmonary disease. Respirology. (2023) 28:428–36. doi: 10.1111/resp.14489

PubMed Abstract | Crossref Full Text | Google Scholar

20. Mummler C, Milger-Kneidinger K. [COPD-exacerbation: therapeutic consequences]. MMW Fortschr Med. (2020) 162:46–9. doi: 10.1007/s15006-020-4430-8

PubMed Abstract | Crossref Full Text | Google Scholar

21. Egenod T, Deslee G, Degano B. [Bronchoscopic COPD treatments]. Rev Mal Respir. (2023) 40:820–33. doi: 10.1016/j.rmr.2023.08.003

PubMed Abstract | Crossref Full Text | Google Scholar

22. Donohue JF. Minimal clinically important differences in COPD lung function. COPD. (2005) 2:111–24. doi: 10.1081/COPD-200053377

PubMed Abstract | Crossref Full Text | Google Scholar

23. Jones PW, Beeh KM, Chapman KR, Decramer M, Mahler DA, Wedzicha JA. Minimal clinically important differences in pharmacological trials. Am J Respir Crit Care Med. (2014) 189:250–5. doi: 10.1164/rccm.201310-1863PP

PubMed Abstract | Crossref Full Text | Google Scholar

24. Celli BR, MacNee W, Agusti A, Anzueto A, Berg B, Buist AS, et al. Standards for the diagnosis and treatment of patients with COPD: a summary of the ATS/ERS position paper. Eur Res J. 23:932–46. doi: 10.1183/09031936.04.00014304

PubMed Abstract | Crossref Full Text | Google Scholar

25. Papaporfyriou A, Bakakos P, Hillas G, Papaioannou AI, Loukides S. Blood eosinophils in COPD: friend or foe? Expert Rev Respir Med. (2022) 16:35–41. doi: 10.1080/17476348.2021.2011219

PubMed Abstract | Crossref Full Text | Google Scholar

26. Patel N. An update on COPD prevention, diagnosis, and management: the 2024 GOLD Report. Nurse Pract. (2024) 49:29–36. doi: 10.1097/01.NPR.0000000000000180

PubMed Abstract | Crossref Full Text | Google Scholar

27. Mycroft K, Krenke R, Gorska K. Eosinophils in COPD-current concepts and clinical implications. J Allergy Clin Immunol Pract. (2020) 8:2565–74. doi: 10.1016/j.jaip.2020.03.017

PubMed Abstract | Crossref Full Text | Google Scholar

28. Hogg JC, Chu F, Utokaparch S, Woods R, Elliott WM, Buzatu L, et al. The nature of small-airway obstruction in chronic obstructive pulmonary disease. N Engl J Med. (2004) 350:2645–53. doi: 10.1056/NEJMoa032158

PubMed Abstract | Crossref Full Text | Google Scholar

29. McDonough JE, Yuan R, Suzuki M, Seyednejad N, Elliott WM, Sanchez PG, et al. Small-airway obstruction and emphysema in chronic obstructive pulmonary disease. N Engl J Med. (2011) 365:1567–75. doi: 10.1056/NEJMoa1106955

PubMed Abstract | Crossref Full Text | Google Scholar

30. Yu X, Xiong Z, Cao J, Wang F, Bao X, Shi J. Pulmonary arterial morphological markers on non-contrast CT predicted acute exacerbations and disease progression in chronic obstructive pulmonary disease: a longitudinal cohort study. Jpn J Radiol. (2025). doi: 10.1007/s11604-025-01841-2. [Epub ahead of print].

PubMed Abstract | Crossref Full Text | Google Scholar

31. Mathews AM, Wysham NG, Xie J, Qin X, Giovacchini CX, Ekström M, et al. Hypercapnia in advanced chronic obstructive pulmonary disease: a secondary analysis of the national emphysema treatment trial. Chronic Obstr Pulm Dis. (2020) 7:336–45. doi: 10.15326/jcopdf.7.4.2020.0176

PubMed Abstract | Crossref Full Text | Google Scholar

32. Pilcher J, Perrin K, Beasley R. The effect of high concentration oxygen therapy on PaCO2 in acute and chronic respiratory disorders. Transl Respir Med. (2013) 1:8. doi: 10.1186/2213-0802-1-8

PubMed Abstract | Crossref Full Text | Google Scholar

33. Xuezhen W. Analysis of curative effect of oxygen therapy combined with pulmonary rehabilitation exercise on patients with chronic obstructive pulmonary disease (in chinese). Gansu Medical Journal. (2024) 43:151–2+156.

Google Scholar

34. Renyi T. Clinical analysis of ltot combined with pulmonary rehabilitation exercise on chronic obstructive pulmonary disease complicated with cor pulmonale (in chinese). World Latest Med Inform. (2019) 19:73+76.

Google Scholar

35. He Yiman YX. Observation on the therapeutic effect of ltot combined with pulmonary rehabilitation exercise on chronic obstructive pulmonary disease complicated with cor pulmonale (in chinese). Chin J Geriatric Care. (2017) 15:21–2.

Google Scholar

36. Shangguang W. Evaluation of ltot combined with pulmonary rehabilitation exercise on chronic obstructive pulmonary disease complicated with cor pulmonale (in chinese). Chin J Mod Drug App. (2020) 14:226–8.

Google Scholar

37. Weiping Y. Effect of oxygen therapy combined with pulmonary rehabilitation exercise on patients with chronic obstructive pulmonary disease complicated with cor pulmonale (in chinese). Chin Commun Doc. (2022) 38:155–7.

Google Scholar

38. Marrara KT, Di Lorenzo VAP, Jaenisch RB, Cabiddu R, de Oliveira Sato T, Mendes RG, et al. Noninvasive ventilation as an important adjunct to an exercise training program in subjects with moderate to severe COPD. Respir Care. (2018) 63:1388–98. doi: 10.4187/respcare.05763

PubMed Abstract | Crossref Full Text | Google Scholar

39. Vitacca M, Kaymaz D, Lanini B, Vagheggini G, Ergun P, Gigliotti F, et al. Non-invasive ventilation during cycle exercise training in patients with chronic respiratory failure on long-term ventilatory support: a randomized controlled trial. Respirology. (2018) 23:182–9. doi: 10.1111/resp.13181

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: oxygen therapy, exercise rehabilitation, chronic obstructive pulmonary disease, exercise capacity, cardiorespiratory function

Citation: Muge Q, Suriguga, Yuqing, Aronggaowa, Taojin and Chen L (2025) A meta-analysis of the effects of long-term oxygen therapy combined with exercise rehabilitation on exercise capacity, cardiopulmonary function, and quality of life in patients with COPD. Front. Med. 12:1640084. doi: 10.3389/fmed.2025.1640084

Received: 06 June 2025; Accepted: 04 September 2025;
Published: 22 September 2025.

Edited by:

Timothy Dean Bryson, Henry Ford Health System, United States

Reviewed by:

Longlong Wang, ETH Zürich, Switzerland
Cláudia Da Costa, Rio de Janeiro State University, Brazil

Copyright © 2025 Muge, Suriguga, Yuqing, Aronggaowa, Taojin and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Taojin, MTUwNDg1MDAzODhAMTYzLmNvbQ==; Lanying Chen, Y2x5NTgzMUAxNjMuY29t

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.