Abstract
BACKGROUND: A recent paper reported that low muscle mass in the erector spinae muscles (ESM) was strongly associated with poor prognosis and declining muscle mass over time in subjects with COPD. However, effects of pulmonary rehabilitation (PR), if any, on ESM mass have not been reported. We hypothesized that PR reduces the annual decline in ESM mass.
METHODS: This was a retrospective cohort study. Thirty-nine subjects with COPD who received PR and underwent chest computed tomography before and after PR were evaluated (rehabilitation group). We also evaluated 39 age-matched subjects with COPD who did not receive PR (nonrehabilitation group). Data were collected from August 2010 until March 2020 in both groups. The ESM cross-sectional area (ESMCSA) was measured using axial computed tomography images, and annual changes were calculated. The 6-min walk distance (6MWD) was measured before and after PR; the minimum clinically important difference was defined as 30 m.
RESULTS: ESMCSA declined in the nonrehabilitation group over time (−116.0 ± 141.2 mm2/y) but increased in the PR group (51.0 ± 95.3 mm2/y; P < .001). The annual increase in ESMCSA was significantly higher among subjects with an increase in 6MWD that exceeded the minimum clinically important difference compared with nonresponders in the rehabilitation group. The annual change in ESMCSA was negatively correlated with comorbidity index, and triple therapy (long-acting β2-agonist/long-acting muscarinic antagonist/inhaled corticosteroid) had a favorable effect on annual change in ESMCSA. Multiple regression analysis revealed that only PR was an independent factor for annual change in ESMCSA.
CONCLUSIONS: ESM mass was shown to decline yearly in subjects with COPD. The annual decline in muscle mass was reduced by PR.
Introduction
Skeletal muscle dysfunction and atrophy are key systemic consequences of COPD and have important clinical implications, such as reduced exercise tolerance, quality of life, and even survival.1 Sarcopenia describes age-related loss of skeletal muscle. The prevalence of sarcopenia is 14.5−24% in patients with COPD and increases with age and disease severity.2,3 Although body mass index (BMI) is indicative of basic nutritional status and has been reported to predict mortality in subjects with COPD,4 BMI is not sensitive to changes in body composition as it can be normal or even increased despite the occurrence of muscle wasting.5,6 McDonald et al7 reported that the pectoralis muscle area was more significantly associated with COPD-related traits, including pulmonary function, dyspnea sensation, and exercise capacity, compared with BMI using existing computed tomography (CT) scans in a previous epidemiologic study. The pectoralis muscle is an accessory respiratory muscle that contributes to inspiration,8 and the cross-sectional area (CSA) of muscles is correlated with pulmonary function parameters, such as FEV1 and FVC.7 Tanimura et al9 reported that the CSA of the mass of the erector spinae muscles (ESMCSA) is the strongest risk factor for mortality compared with other clinical parameters including lung function, , BMI, modified Medical Research Council (mMRC) dyspnea scale, comorbidity indexes, St George Respiratory Questionnaire score, prior COPD exacerbations, and smoking history in subjects with COPD. These investigators also reported that ESMCSA decreased yearly in subjects with COPD.10 Because the trunk muscles, including ESM, have an important role in stabilizing the body, maintaining posture, and controlling spinal and pelvic movement, atrophy of these muscles may lead to an increase in the risk of falling and disability in activities of daily living in elderly people and patients with COPD.11,12
Both pectoralis muscles and ESMs can be assessed with chest CT. Because chest CT is usually performed to diagnose and evaluate emphysema and to exclude other critical pulmonary diseases, such as lung cancer and pulmonary fibrosis, these images may simultaneously be used to evaluate local skeletal muscle mass without a need for additional radiation exposure.
A previous study reported that pulmonary rehabilitation (PR) enhanced lower limb muscle function and increased muscle mass13; however, the effect of PR on ESMs or pectoralis muscles has not been evaluated. We hypothesized that PR increases muscle mass for both ESM and pectoralis muscles and reduces the annual decline in muscle mass.
Quick Look
Current Knowledge
Skeletal muscle dysfunction and atrophy are key systemic consequences of COPD and have important clinical implications, such as reduced exercise tolerance, quality of life, and survival. Erector spinae muscle (ESM) mass is the strongest risk factor for mortality and decreases year by year.
What This Paper Contributes to Our Knowledge
We used computed tomography to evaluate the effect of pulmonary rehabilitation (PR) on the annual change of ESM mass in subjects with COPD. ESM cross-sectional areas declined in subjects who did not receive PR and increased in subjects who received PR. Our results indicate that PR can effectively improve ESM in subjects with COPD.
Methods
Subjects and Study Design
This was a retrospective cohort study. The Ethics Committee of Kindai University approved this study (approval number, R02-012) and authorized a waiver of consent. Confidentiality of subject information was guaranteed. The study was conducted in accordance with the Declaration of Helsinki. The basic information in this study has been registered with the University Hospital Medical Information Network Clinical Trials Registry (UMIN 00042233).
Medical records were reviewed to evaluate the effect of PR on ESM mass in subjects with COPD. A total of 89 subjects (age range, 51–91 y) with COPD were referred to the out-patient PR program at Kindai University Hospital (Osaka, Japan) from August 2010 until March 2020. Seventy-one subjects completed the PR program (Fig. 1). Subjects were considered to have completed the program if they attended at least 12 sessions, given that most professional society guidelines recommend a minimum of 12 supervised sessions.14,15 Among the 71 subjects, 41 had 2 CT scans (before and after the PR program). Chest CT scans were performed at the physician’s discretion to evaluate emphysema or to exclude the possibility of lung cancer, irrespective of this study. Therefore, a CT scan was not performed routinely in PR program attendees, and the interval between CT scans varied because of the retrospective nature of this observational study. Two subjects were excluded from the analysis because the CT scans were performed at the time of admission during a COPD exacerbation. A previous study reported that ESMs are affected by such exacerbations.10 Therefore, 39 subjects who completed the PR program were evaluated in this study (rehabilitation group). We also evaluated 39 age-matched out-patients with COPD who did not participate in the PR program and had 2 CT scans (non-rehabilitation group). All subjects in both groups were diagnosed with COPD and treated in the Department of Respiratory Medicine and Allergology of Kindai University Hospital. The diagnosis and classification of COPD were made according to the Global Initiative for Chronic Obstructive Lung Disease (GOLD) guidelines.16
The out-patient PR program was conducted twice a week for 12 weeks (24 sessions), and the PR maintenance program consisted of supervised endurance training once or twice per month. Each session consisted of respiratory muscle stretch gymnastics,17 followed by endurance training involving the lower limbs through exercise on a treadmill or ergometric bicycling at 60–70% of peak work load for 20–40 min, and strength training of the upper and lower limbs for 10−20 min, as previously described in previous study.14 Respiratory muscle stretch gymnastics is designed to decrease chest wall stiffness, particularly of the respiratory muscles, because the majority of subjects with COPD had reduced chest wall mobility.17,18
Clinical Variables
Clinical data were retrieved from the medical files and included pulmonary function, current medications, smoking history, medical history, past COPD exacerbations, and comorbidities. Comorbidity indices (Charlson comorbidity index and COPD-specific comorbidity test) were calculated, as reported previously.14,19,20 We monitored severe exacerbations requiring hospitalization. All patients included in the study underwent spirometry (CHEST AC-55V; Chest, Tokyo, Japan), according to the method described in the American Thoracic Society guidelines.21 The single-breath diffusing capacity of the lung for carbon monoxide (DLCO) was also measured.22 The values for total lung capacity, FVC, residual volume, FEV1, and DLCO were related to the percentage of the predicted values.23,24 Exercise tolerance was assessed before and after the PR program by measuring the 6-min walk distance (6MWD) in accordance with the American Thoracic Society guidelines.25 The 6MWD was measured in the rehabilitation group, but not in the non-rehabilitation group. The minimum clinically important difference in the 6MWD was defined as an increase > 30 m.26 Responders to pulmonary rehabilitation for exercise capacity were defined as having a 30-m increase in the 6MWD.
CT Scan Acquisition and Analysis
ESMCSA was measured using axial chest CT images at the level of the lower margin of the 12th thoracic vertebra with a manually drawn line depicting the ESMCSA.9,10,27 ESMCSA is presented as the sum of the bilateral sectional areas calculated using the medical imaging and information management system (SYNAPSE 4.4.3, Fujifilm Medical, Tokyo, Japan; see the supplementary materials at http://www.rcjournal.com).9,10,27 The CSA of the pectoralis muscle (PMCSA) was also determined based on a single axial CT scan slice above the aortic arch, as described previously7,9 (see the supplementary materials at http://www.rcjournal.com). Emphysema was quantified by calculating the percentage of the low-attenuation area, determined using a cutoff value of –950 HU on whole-lung CT images using Synapse Vincent (Fujifilm Medical).9 Annual changes in the CSA of muscles were calculated from 2 CT scans using the following formula: (CSA in second CT scan – CSA in first CT scan)/(interval between two CT scans [in years]).
Statistical Analysis
Descriptive statistics were calculated for continuous variables, as represented by mean ± SD (assumed to follow a normal distribution), median (interquartile range [IQR]; assumed to follow a non-normal distribution), and categorical variables by n (%). To compare 2 groups for parameters, we used the Student t test (for continuous variables with normal distributions), the Wilcoxon rank-sum test (for non-normal distributions), and the chi-square test (for categorical variables). The Shapiro-Wilk test was used to determine whether the distribution was normal. For comparison of the annual change of muscle CSA between groups (with vs without PR; responders vs non-responders), we performed linear regression using study groups and baseline values as explanatory variables. Correlations between parameters were assessed using Spearman’s correlation analysis. Correlation coefficients < 0.4 were considered to be weak, and those from 0.4 to 0.69 were considered to be moderate.28 We also used standard least squares fitting with the restricted maximum likelihood method for multiple regression analyses. Comparisons between COPD stages or mMRC dyspnea scales were performed using analysis of variance. The minimum clinically important difference for 6MWD was defined as an increase > 30 m.26 Analyses were performed using JMP 11.2.0 software (SAS Institute, Tokyo, Japan).
Results
Subject Characteristics
There were no significant differences between the rehabilitation and nonrehabilitation groups with respect to age, gender, BMI, smoking history, percent of predicted FEV1, percent of predicted FVC, and drug treatment (Table 1). However, the rehabilitation group had more severe dyspnea (per the mMRC), more advanced GOLD stage, a greater proportion requiring long-term oxygen therapy, and more frequent exacerbations, while the nonrehabilitation group had a greater prevalence of asthma-COPD overlap and higher comorbidity indices (ie, Charlson comorbidity index). The proportion of prescriptions for inhaled drugs did not differ between the groups. With the exception of 1 subject in the nonrehabilitation group, the subjects were treated with a long-acting β2 agonist (LABA) or a long-acting muscarinic antagonist (LAMA). The interval between the first and second CT scans did not differ between the groups (rehabilitation group 2.3 ± 1.3 y vs nonrehabilitation group 2.1 ± 1.2 y; P = .24). In the PR group, the interval between the first CT and the start of PR was 1.4 ± 1.0 y, and between the end of PR and the second CT the interval was 0.9 ± 0.7 y.
Effect of PR on Annual Changes in Muscle CSA
There were no significant differences between the rehabilitation and nonrehabilitation groups with respect to the mean values of ESMCSA or PMCSA at the time of the first or second CT scan (Table 1). ESMCSA decreased in 13 subjects in the rehabilitation group and in 33 subjects in the nonrehabilitation group. The annual changes in ESMCSA and PMCSA of the rehabilitation group (51.0 ± 95.3 and 15.3 ± 171.9 mm2/y, respectively) were significantly higher than the changes in ESMCSA and PMCSA of the nonrehabilitation group (−116.0 ± 141.2 mm2/y, P < .001; and −138.4 ± 317.4 mm2/y, P = .02, compared with the rehabilitation group) (Fig. 2, Table 2). We could not measure PMCSA in 2 subjects during the second CT scan because it was an abdominal CT scan. The annual changes in ESMCSA were significantly greater in subjects for whom the 6MWD increased more than the minimum clinically important difference (responders) (n = 17) compared with the other subjects (n = 22) in the rehabilitation group (responders 94.0 ± 88.0 vs nonresponders 17.0 ± 88.4 mm2/y, P < .001) (Fig. 3; see the supplementary materials at http://www.rcjournal.com). The annual changes in PMCSA were not different between responders and nonresponders (responders 52.1 ± 206.4 vs nonresponders −15.9 ± 133.7 mm2/y, P = .18). There was no significant difference between responders and nonresponders in terms of background characteristics (see the supplementary materials at http://www.rcjournal.com).
Effect of Clinical Background on Annual Changes in Muscle CSA
The annual changes in ESMCSA or PMCSA were not affected by treatment with inhaled corticosteroids (ICS) (P = .26 and .28, respectively; see the supplementary materials at http://www.rcjournal.com). The annual changes in ESMCSA were significantly greater in the 29 subjects who received triple therapy (LABA/LAMA/ICS) compared with the 49 subjects who did not (23.6 ± 131.8 vs −66.4 ± 145.2 mm2/y, respectively; P < .001). The median (IQR) of past COPD exacerbations was greater in subjects treated with triple therapy compared with other therapies, but this difference was not statistically significant (0 [IQR 0–2] vs 0 [IQR 0–1], respectively; P = .22), and a higher percentage of subjects with asthma-COPD overlap were treated with triple therapy (62.1%) compared with other therapies (16.3%) (P < .001). The annual changes in PMCSA were unaffected by triple therapy (P = .26). When combined with asthma (ie, asthma-COPD overlap), the annual changes in ESMCSA or PMCSA were unaffected (P = .81 and .21, respectively). Likewise, annual changes in ESMCSA or PMCSA were unaffected by long-term oxygen therapy (P = .82 and .40, respectively), nor were they different as a function of the GOLD stage (P = .82 and .40, respectively) or mMRC dyspnea scale (P = .12 and .68, respectively) (see the supplementary materials at http://www.rcjournal.com). In summary, in terms of clinical background, other than PR, only triple therapy affected annual changes in ESMCSA.
Correlation Between Clinical Parameters and Muscle CSA
ESMCSA at the first CT scan (baseline) was moderately and positively correlated with BMI, 6MWD, and PMCSA; a weak positive correlation was observed with inspiratory capacity; and a weak negative correlation was evident with residual volume (Table 3). PMCSA at baseline was moderately and positively correlated with BMI; there was a weak positive correlation with percent of predicted FEV1; and weak negative correlations were observed with residual volume and low-attenuation area percentage.
The annual changes in ESMCSA were weakly and positively correlated with mMRC dyspnea scale and COPD exacerbations during the period between the 2 CT scans; a moderate positive correlation was observed with annual changes in PMCSA; and a weak negative correlation existed with the Charlson comorbidity index (Table 4).
We performed multiple regression analysis of the annual changes in ESMCSA and determined that PR alone was an independent variable (Table 5).
Discussion
Our results indicate that ESMCSA and PMCSA decreased year by year, and the annual changes of ESMCSA and PMCSA were different between the PR group and the nonrehabilitation group, in favor of the PR group. This is the first study to report the effect of PR on the annual changes in ESMs, which are major trunk muscles. Skeletal muscle mass is an independent predictor of mortality and frailty.1,29 Exercise training alone or when combined with other treatments have favorable effects on muscle mass and function.13 However, there are no reports on the effect of PR on trunk muscle mass in subjects with COPD. A recent study reported that ESMCSA is the strongest predictor of mortality in subjects with COPD and declines year by year.10 The CSA of ESMs has also been reported to predict mortality after liver transplants and vascular surgery.30,31 The health trajectory of subjects with COPD may be improved as a result of a PR program that maintains ESMCSA. A recent report showed that PR after hospitalization improved survival in subjects with COPD.32 The mechanism by which a usual PR program can improve ESMCSA remains unclear. However, reports indicate that PR for subjects with COPD improves balancing ability,33-35 which suggests that trunk muscle improvement may be implicated. Usual PR programs that include endurance training may be able to improve trunk muscles similarly to a more specific trunk flexor rehabilitation program.36
ESMCSA and PMCSA declined in some subjects even after completing the PR program. The magnitude of response to exercise training in COPD is highly variable, with some subjects achieving little or no benefit in terms of strengthening of the lower limb muscles.37 The lack of response to exercise training may relate to the inability to tolerate sufficient intensity or duration of training, or to poor adherence to the training intervention.37 We measured the CSA of muscles before and after the PR program; however, because this was a retrospective study, the measurements were performed several months after the PR program in some subjects, and the CT scans were performed regardless of the PR. Therefore, the effect of the PR program might have diminished by the time measurements were recorded for some subjects.
Triple therapy (ICS/LABA/LAMA) had a favorable effect on the annual change in ESMCSA, although the effect was not significant based on multiple regression analysis. Triple therapy was prescribed for subjects with COPD and a concomitant asthmatic component or frequent COPD exacerbations.16 In other words, triple therapy was prescribed for subjects who had more severe disease and symptoms. We reasoned that the annual decline in ESMCSA was greater in subjects with severe COPD complicated by frequent exacerbations,10 but the annual change in ESMCSA was not correlated with the presence of asthma-COPD overlap or the number of exacerbations. Recent randomized trials have reported the superior effects of triple therapy on pulmonary function, quality of life score, and all-cause mortality compared with LAMA alone, LAMA/LABA, or LABA/ICS.38,39 The superiority of triple therapy over single- or double-agent therapy may affect the annual change in ESMCSA; however, a randomized controlled trial will be necessary to more definitively prove or disprove this hypothesis.
ESMCSA at baseline was moderately correlated with 6MWD, while PMCSA was weakly correlated with pulmonary function parameters such as FEV1, residual volume, and percentage of the low-attenuation area. This finding is not surprising because ESMs are important for the maintenance of normal posture during walking, and the PMs are important accessory inspiratory muscles of respiration. Because the trunk muscles, including ESMs, have an important role in stabilizing the body, maintaining posture, and controlling spinal and pelvic movements, their atrophy can increase the risk of falling and limitations of activities of daily living in elderly people11 and in those with COPD.33 However, these correlations were mostly weak. Muscle CSA does not directly reflect the strength of a particular muscle because no muscle works in isolation, rather groups of muscles are generally recruited to work together. The annual change in ESMCSA was weakly and negatively correlated with the Charlson comorbidity index. In this study, we included subjects with a history of cancer or cardiovascular disease; therefore, the Charlson comorbidity index was higher than for previous reports in which such patients were excluded.9,10 Comorbidities may accelerate a decline in ESMCSA. Against all expectations, COPD exacerbations during the period between the 2 CT scans were weakly and positively correlated with the annual change in ESMCSA. This result may be caused by the more frequent COPD exacerbations in the PR group compared with the nonrehabilitation group (Table 1).
There are several methods to evaluate muscle mass. Specifically, these include bioelectrical impedance analysis, dual-energy radiography absorptiometry, magnetic resonance imaging (MRI), and B-mode ultrasound. These methods are widely used to quantify both total and local skeletal muscle mass. Most reports of ESMCSA are based on CT scans or MRIs.9,10,27,30,31,40-43 Several studies have reported a moderate correlation between ESMCSA and isometric muscle strength using CT40,43 or MRI.42 The downside of CT scans is exposure to radiation, while MRIs are costly in terms of time and money. Chest CT scans are often used in Japan to diagnose lung disease, assess changes in emphysematous COPD, diagnose interstitial pneumonia, and for the early detection of lung cancer in the elderly because chest MRI remains inaccessible for the majority of subjects with COPD. Therefore, we used existing CT images to evaluate muscle mass without adding to radiation exposure. Ultrasound imaging may be useful for measuring CSA whenever CT scans are performed infrequently.44 We could not measure muscle mass in healthy control subjects because of the nature of this retrospective study. Tanimura et al9 reported ESMCSA in subjects without respiratory disease who were nonsmoking controls. ESMCSA in control subjects (3,920 ± 698 mm2) was greater than that in the subjects with COPD in this study. Trunk stabilization muscles are categorized as deep or superficial.45 ESMs are categorized as superficial muscles because they do not directly attach to vertebrae and stabilization is limited to pressurized spinal segments. The deep trunk muscles consist of the transversus abdominis, multifidus, quadratus lumborum, and rotatores muscles, which attach to individual lumbar vertebrae and provide control and rigidity for spinal segments.45 There are no reports pertaining to the stabilizing role of deep trunk muscles on the health trajectory of patients with COPD, and more research is needed to elucidate this functionality.
There are some limitations of our study. It is possible that it was underpowered, given the relatively small number of subjects. The retrospective study design included subjects with clinical backgrounds that were quite different in the rehabilitation and nonrehabilitation groups. Severe dyspnea symptoms (per mMRC) and severe disease (GOLD III−IV) were more prevalent in the rehabilitation group than in the nonrehabilitation group. The Charlson comorbidity index was higher in the nonrehabilitation group than the rehabilitation group. Symptomatic subjects with severe disease might have been referred for PR, and comorbidities could have interfered with their attendance at sessions during the PR program. Although these background differences may have affected the results, associated factors were excluded in the multiple regression analysis. The interval between the 2 CT scans also varied, and these differences might have influenced the results. We did not assess muscle areas using multi-slice CT to measure muscle volume. CSA was analyzed via manual shading of a specific muscle area, as done in previous studies.9,10,27 A standard method for measuring ESMCSA has yet to be established. Automated programs to measure ESMCSA are necessary to allow more objective and precise analysis. Lastly, there were some missing data, and we were unable to analyze low-attenuation area percentage and DLCO in some subjects; these missing data may have affected the study results.
Conclusions
ESMCSA and PMCSA declined year by year in subjects with COPD, although it increased in subjects in the pulmonary rehabilitation group, suggesting that a pulmonary rehabilitation program may improve the mass of these muscles.
Footnotes
- Correspondence: Yuji Higashimoto MD PhD, 377-2 Onohigashi, 5898511, Osakasayama, Japan. E-mail: yhigashi{at}med.kindai.ac.jp
Supplementary material related to this paper is available at http://www.rcjournal.com.
Dr Tohda has disclosed relationships with Kyorin Pharmaceutical, Teijin Pharma, Meiji Seika Pharma, Boehringer Ingelheim, Daiichi Sankyo, Astellas Pharma, and Pearl Therapeutics. The other authors have disclosed no conflicts of interest.
- Copyright © 2021 by Daedalus Enterprises