Gender modifies the effect of body mass index on lung function decline in mild-to-moderate COPD patients: a pooled analysis
Respiratory Research volume 22, Article number: 59 (2021)
Low body weight is associated with poor prognosis in patients with chronic obstructive pulmonary disease (COPD). However, it is not known whether gender modifies this relationship.
We pooled data of 8686 COPD patients from 7 studies with a median length of 36-months of follow up. Using a longitudinal natural cubic spline regression model, we examined the dose–response relationship between body mass index (BMI) and the rate of decline in forced expiratory volume in one second (FEV1) in patients with GOLD 1 and 2 disease, stratified by gender and adjusted for age, smoking status, and cohort effects.
There was an inverse linear relationship between BMI and the rate of FEV1 decline in GOLD Grades 1 and 2, which was modified by gender (p < 0.001). In male patients, an increase of BMI by 1 kg/m2 reduced FEV1 decline by 1.05 mL/year (95% CI 0.96, 1.14). However, in female patients, BMI status did not have a clinically meaningful impact on FEV1 decline: an increase of baseline BMI by 1 kg/m2 reduced FEV1 decline by 0.16 ml/year (95% CI 0.11, 0.21). These gender-modified relationships were similar between GOLD 1 and 2 patients, and between current and former smokers.
In mild to moderate COPD, higher BMI was associated with a less rapid decline of FEV1 in male patients whereas this association was minimal in females patients. This gender-specific BMI effect was independent of COPD severity and smoking status.
Chronic obstructive pulmonary disease (COPD) is one of the leading causes of morbidity and mortality in the world . According to the World Health Organization, over the next 20 years COPD will rise from the fifth to third leading cause of death worldwide . Among newly diagnosed COPD patients, over half had mild disease, one third had moderate, and 10% of patients had severe or very severe disease . The progression of COPD is characterized by an accelerated decline in lung function as indicated by forced expiratory volume in one second (FEV1). On average, patients with mild to moderate COPD experience a faster decline in FEV1 over time than those with the more severe form of COPD . Modification of lung function trajectory early in its course provides an opportunity to ameliorate patient’s quality of life and extend their life expectancy. It is now well-recognized that cachexia is a significant risk factor for poor outcomes including mortality in COPD patients [6,7,8]. However, while multiple studies report that those with low body mass index (BMI) may be at risk of COPD progression and those who are obese may be protected [5, 6], others have shown no significant association between BMI and FEV1 decline [11, 12]. One reason for the controversy is the fact that a decline in lung function varies according to certain factors, most notably smoking, gender, and disease stage [5,6,7]. Knowing whether BMI is a significant risk factor for COPD progression in patients with mild and moderate COPD is important, because it represents a potentially modifiable trait. Thus, we pooled individual-level data from 7 large international studies into a single, combined dataset, and examined the dose–response relationship between BMI and the rate of FEV1 decline in patients with mild and moderate COPD, according to gender and smoking status.
Study design and settings
This was a pooled analysis of patient-level data from 6 randomized controlled trials (RCTs) and 1 non-interventional prospective study. The RCTs were part of the Inhaled Steroids Effect Evaluation in COPD study (ISEEC)  and included the Lung Health Study (LHS, n = 5 594 patients, 11 years of follow-up) , the European Respiratory Society study on COPD (EUROSCOP, n = 1039 patients, 3 years of follow-up) , Inhaled Steroids in Obstructive Lung Disease in Europe (ISOLDE, n = 591 patients, 3 years of follow-up) , Copenhagen City Lung Study (CCLS, n = 225 patients, 3 years) , studies by Calverley et al. (n = 336 patients, 12 months of follow-up)  and Szafranski et al. (n = 292 patients, 12 months of follow-up) . We additionally included data from a non-RCT, the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE) study , which was a 3-year non-interventional study of 2652 patients with stable COPD from 46 centers across 12 countries . By pooling data from these 7 cohorts, the final study sample included patients with mostly mild and moderate COPD according to the Global Initiative for Chronic Obstructive Lung Disease (GOLD) severity grades at baseline (GOLD Grade 1, FEV1 ≥ 80 percent of predicted at 25 years of age, mild COPD; and Grade 2, FEV1 50–79 percent predicted, moderate COPD). We excluded patients in GOLD Grades 3 and 4 because they demonstrated different FEV1 trajectories over time compared with GOLD 1 and 2 patients (Table S1). For parsimony, we combined GOLD Grades 1 and 2 into one group because their FEV1 decline rates were similar (see Additional file 1: Table S1 and Figures S1, S2 for the observed rate of FEV1 decline between GOLD Grades 1 and 2). Further, patients were included if they had a valid measurement of BMI (≥ 10 kg/m2) at baseline, and 3 or more measurements of FEV1 across 3 or more different time points, which enabled a stable estimate of the FEV1 decline slope. Of note, spirometry standards and quality assurance were similar between the ISEEC trials and ECLIPSE. The Hankinson’s prediction equation was used to calculate percent of predicted FEV1 across all studies, using reference values for Caucasians, African–Americans, and Mexican–Americans derived from 7429 asymptomatic, non-smoking participants in the Third National Health and Nutrition Examination Survey (NHANES III) . Because the provision of inhaled corticosteroids (ICS) did not affect FEV1 decline over time in ISEEC , this treatment intervention was not considered a confounder in the present study. Thus, we included patients from both the treatment and placebo arms of the RCTs.
The primary outcome was the change in the absolute value of post-bronchodilator FEV1 over time, which was assessed in all included studies through standardized spirometry measurements as described previously . We did not evaluate changes in the percent predicted FEV1 over time because this parameter is adjusted for height, which is a key component of the BMI measurement.
The primary exposure, BMI at the baseline visit, was expressed as a continuous variable, and was obtained by dividing patient’s weight (Kg) by height squared (m2). We also examined the effects of gender (male, female) and their 2-way interaction effects with BMI, controlling for the confounding effect of cigarette smoking. As there could be significant between-study heterogeneity in the distribution of risk factors, study design, laboratory protocols, and enrollment period across the individual studies, we included relevant patient-level covariates in the model. This included baseline age as a continuous variable, follow-up years, and a categorical variable indicating cohort membership to account for other unobserved between-study heterogeneity.
All analyses were performed using SAS 9.3 (SAS Institute Inc, Cary, NC, United States). Detailed descriptions of the statistical methods are provided in Additional file 1: Sect. 1. The criterion for statistical significance was a two-tailed p-value of 0.05 or less. Descriptive statistics were calculated, and comparisons were made using Pearson Chi-square tests for categorical variables and Kruskal–Wallis tests for continuous variables.
First, to examine the unadjusted relationship between BMI and the rate of FEV1 decline, we fitted a linear mixed-effects model with two predictors: a random intercept (corresponding to baseline FEV1 value) and a random slope of follow-up time (corresponding to the rate of FEV1 decline). From this analysis, we obtained the individual rate of FEV1 decline (i.e., the random slope), and plotted it against the individual’s baseline BMI value in a scatter plot.
Next, to evaluate the covariate-adjusted, averaged, potentially non-linear relationship between BMI and the rate of FEV1 decline, we applied natural cubic spline models. This longitudinal regression model had FEV1 values as the dependent variable, and contained the following independent variables: a spline function of BMI (primary exposure in which knots were placed across every 5th percentile), age, gender, smoking status, follow-up years, cohort status, and 2nd-order and 3rd-oder interactions terms between BMI, gender, smoking status, and follow-up years. We did not include an interaction term between BMI and age because based on a preliminary variable selection process, we found that its inclusion reduced model fit. These natural cubic splines produced smooth curves, which took into account the nonlinear components of the relationship between the exposure variables and the outcome. Next, outcomes were derived as the covariate-adjusted dose–response curves of BMI-rate of FEV1 decline using a robust causal inference technique named the G-computation . In specific, based on the regression results, we obtained the individual FEV1 values at the 12-month follow-up, and calculated the individual rate of FEV1 decline as the instantaneous rate of change in FEV1 over a minimal duration of time (0.5E−5). We then evaluated the averaged changes in the rate of FEV1 decline per unit increase in baseline BMI value in dose–response curves according to gender. To enable the construction of 95% confidence bands of the dose–response curves, which accounted for the longitudinal, correlated data structure, we applied 1000 rounds of bootstrapping, which has been previously shown to efficiently handle correlated time-series data . In a secondary analysis, the results were further stratified by smoking status to determine the impact of smoking on the relationship between BMI and FEV1 decline.
Characteristics of the study population
This pooled analysis comprised of 8686 COPD patients (Fig. 1, cohort selection). Table 1 presents the baseline characteristics of the study population. The mean baseline age was 51.9 years (SD = 9.1); 37% were women; and 56% were current smokers. The average BMI at baseline was 25.7 kg/m2 (SD = 4.3). The median follow-up time was 36 months. A total of 3674 (42%) patients were in GOLD Grade 1; 5012 (58%) were in GOLD Grade 2.
Table 2 presents the gender-specific baseline characteristics and the observed rate of FEV1 decline of the combined GOLD Grades 1 and 2 samples. Male patients experienced a significantly faster decline in absolute FEV1 values compared to female patients (mL/year, − 36.6 vs − 29.2, p-value < 0.001). Male and female patients had similar mean ages (51.8 vs. 51.9 years) and proportionality of current smokers (54% vs. 57%). However, the majority of males had normal body weight (58% of BMI between 18.5 and 25.0 kg/m2), whereas the majority of females were more likely to be overweight or obese (60% of had a BMI above 25.1 kg/m2). Additional file 1: Table S2 shows the observed baseline FEV1 and the rate of FEV1 decline by smoking status and BMI level. Underweight individuals (BMI < 18.5 kg/m2) had a lower baseline FEV1 compared to other BMI groups (L, 2.17 vs 2.56–2.73), while obese individuals (BMI ≥ 30.1 kg/m2) had a slower rate of decline than those who were overweight (BMI 25.1–30.0 kg/m2), normal (BMI 18.5–25.0 kg/m2) or underweight (mL/year, − 27.5 vs − 34.4, − 35.3, − 34.6, respectively). Current smokers had a much more rapid decline than ex-smokers (mL/year, − 40.9 vs − 24.9).
Gender-modified effects of BMI on the rate of FEV1 decline
Figure 2 illustrates the fitted dose–response relationship between BMI (x-axis) and the rate of decline in absolute FEV1 values (y-axis, with confidence intervals shown in error bars) in combined GOLD Grades 1 and 2. The confidence interval of the curve was wider at both ends of the BMI scale, reflecting a greater variance of the association at extreme BMI values, probably due to the smaller number of patients at the extremes.
Gender significantly modified the relationship between BMI and FEV1 decline (p < 0.001). In male patients, the dose–response curve depicted a reverse association between BMI and rate of decline in absolute FEV1 values. This relationship was mostly linear, except for a slight fluctuation in the line at a BMI of 25 kg/m2 (Fig. 2). The slope of this curve showed that an increase of BMI by 1 kg/m2 reduced FEV1 decline by approximately 1.05 mL/year (95% CI: 0.96, 1.14). In female patients, the slope of the curve showed that BMI had a very small (and clinically insignificant) effect: an increase of BMI of 1 kg/m2 reduced FEV1 decline by only 0.16 ml/year (95% CI: 0.11, 0.21).
In a secondary analysis, the effects of BMI were stratified by GOLD Grades and smoking status (Fig. 3, upper panel, male patients [left, GOLD 1; right, GOLD 2], lower panel, female patients [left, GOLD 1, right, GOLD 2]). As for the effects of smoking, within GOLD Grade 1, the gender-specific dose–response curves appeared in parallel between current smokers and ex-smokers, though the curves were statistically different from each other (p < 0.001). Within GOLD Grade 2, there was a notable narrowing of the gap at higher BMI levels between current smokers and ex-smokers. This suggests that, conditional on COPD severity and gender, the additional impact of smoking on the relationship between BMI and FEV1 decline was small in GOLD Grade 1, whereas higher BMI level appeared to have a slight protective effect on FEV1 decline in GOLD Grade 2 current smokers, in particular among male smokers. Of note, the highest risk of decline was observed in underweight male smokers with GOLD 1 disease, who experienced, on average > 70 ml/year decline in FEV1 (Fig. 3).
In a large pooled analysis of seven multinational prospective studies, we showed that the relationship between BMI and FEV1 decline in mild to moderate COPD was significantly modified by gender. In females, BMI had no material impact on FEV1 decline; whereas in males, a 1 kg/m2 increase in BMI was associated with a reduction of approximately 1 mL/year in the rate of decline in absolute FEV1 values. The protective effect of high BMI appeared to be more prominent in GOLD 2 male smokers. Of note, underweight male current smokers in GOLD 1 were at the highest risk of disease progression. We did not include percent of predicted FEV1 as an outcome given that it was derived from height.
Previous studies have largely focused on effects of smoking as a risk factor for COPD disease progression. Consistent with previous findings, we found that smokers experienced over 10 mL/year faster decline in FEV1 compared with sustained quitters . We extend these previous findings in this large pooled analysis by showing that BMI significantly impacts FEV1 decline in males but has only a minimal effect in females with mild to moderate COPD. For example, the annual decline of FEV1 was at least 6 mL/year faster in underweight (BMI < 19 kg/m2) male patients than those who were overweight (BMI ≥ 25 kg/m2).
The mechanisms by which BMI modifies FEV1 decline have not been fully elucidated. High BMI may represent better nutrition status , increased body fat, muscle mass, and/or bone mineral density , while lower BMI may indicate poor nutrition and skeletal muscle loss that leads to accelerated lung function loss . The lower baseline FEV1 of underweight patients may also indicate childhood conditions such as restricted early-life growth that leads to poor lung development . In addition, emphysema might have also played a role, because emphysema is significantly associated with reduced BMI , and is also more common in male patients . Furthermore, sex-modified effects of BMI on FEV1 decline may be related to anatomical differences in the large airways that have been observed between men and women, such as smaller central airways and lower thoracic volume in females [26, 27]. Due to concerns regarding reverse causality, we were unable to study the effects of dynamic changes in BMI (i.e., increased lung burden caused weight loss). Notwithstanding these important mechanistic issues, BMI is easy to measure, accurate and reproducible. As such, BMI may be used clinically to identify COPD patients (particularly males) at risk for rapid disease progression.
A major strength of this analysis was that it pooled individual-level data from seven high-quality long-term studies, which reduced heterogeneity and yielded more reliable results than the previously reported meta-analysis . The inclusion of ECLIPSE, a real-world prospective cohort, added to the external validity and generalisability of results. Another strength was the use of a robust and powerful statistical approach, which adapted longitudinal analyses to principles of restricted cubic splines. This enabled us to extend measures of static, cross-sectional dose–response relationship to more pragmatic metrics representing the impact of BMI on the progression of COPD (ie., the rate of FEV1 decline). Importantly, our analyses were stratified by gender and adjusted for various confounders including GOLD grades of severity, cohort and calendar effects, smoking status and its interactions with BMI. This enabled us to tease out the differential effects of BMI between genders, which were not well known previously.
Our findings need to be interpreted within the context of certain limitations. First, BMI is an approximate measure for nutritional status, because it is unable to distinguish between fat and fat-free mass or its distribution. Future studies should consider other anthropometric measurements in female patients with mild to moderate COPD and investigate the potential roles of muscle and fat mass in the gender-specific progression of COPD. Second, we assessed BMI at baseline and weight may change dynamically over time. This was not necessarily a limitation, because this “intention-to-treat” approach protects against reverse causality and provides more valid inference to our research question. Third, this analysis could not adjust for unrecorded potential confounders such as comorbidities and exacerbations. However, while the comorbidities of low BMI COPD patients are different from those of patients with a high BMI , the direct impact of comorbidities on lung function decline is largely unknown. Fourth, non-smokers were not included in this study as data were unavailable, which may limit the generalizability of results. Finally, we limited the sample to patients with 3 or more FEV1 measurements to evaluate a possible dose–response curve using a cubic spline analysis. Only 3% of the entire cohort comprised of patients with 2 (or fewer) FEV1 measurements and these patients had shorter follow-up period compared to those with 3 or more FEV1 measurements (median follow up, 13 vs. 36 months), but the former has similar baseline FEV1 and rate of FEV1 decline compared to the latter whose follow-up period was up to 13 months.
We conclude that reduced BMI is a significant risk factor for accelerated decline in lung function but is modified by gender. Underweight male smokers with GOLD 1 disease are at the highest risk of rapid COPD progression and thus should be followed closely and be strongly counseled for smoking cessation.
Availability of data and materials
The data that support the findings of this pooled study are available from the individual trials, including CCLS, ECLIPSE, EUROSCOP, ISOLDE, LHS, and studies conducted by Calverley et al., Szafranski et al., but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of these trials.
Vogelmeier CF, Criner GJ, Martinez FJ, Anzueto A, Barnes PJ, Bourbeau J, et al. Global strategy for the diagnosis, management, and prevention of chronic obstructive lung disease 2017 report. GOLD Executive Summary. Am J Respir Crit Care Med. 2017;195(5):557–82.
WHO | Burden of COPD. WHO. https://www.who.int/respiratory/copd/burden/en/. Accessed 12 Mar 2019.
van den Boom G, Rutten-van Mölken MP, Tirimanna PR, van Schayck CP, Folgering H, van Weel C. Association between health-related quality of life and consultation for respiratory symptoms: results from the DIMCA programme. Eur Respir J. 1998;11(1):67–72.
Tantucci C, Modina D. Lung function decline in COPD. Int J Chron Obstruct Pulmon Dis. 2012;7:95–9.
Pérez-Padilla R, Fernandez-Plata R, Montes de Oca M, Lopez-Varela MV, Jardim JR, Muiño A, et al. Lung function decline in subjects with and without COPD in a population-based cohort in Latin-America. PLoS ONE. 2017;12(5):e0177032.
Casanova C, de Torres JP, Aguirre-Jaíme A, Pinto-Plata V, Marin JM, Cordoba E, et al. The progression of chronic obstructive pulmonary disease is heterogeneous: the experience of the BODE cohort. Am J Respir Crit Care Med. 2011;184(9):1015–21.
Sin DD, Wu L, Anderson JA, Anthonisen NR, Buist AS, Burge PS, et al. Inhaled corticosteroids and mortality in chronic obstructive pulmonary disease. Thorax. 2005;60(12):992–7.
Connett JE, Kusek JW, Bailey WC, O’Hara P, Wu M. Design of the Lung Health Study: a randomized clinical trial of early intervention for chronic obstructive pulmonary disease. Control Clin Trials. 1993;14(2 Suppl):3S-19S.
Pauwels RA, Löfdahl C-G, Laitinen LA, Schouten JP, Postma DS, Pride NB, et al. Long-term treatment with inhaled budesonide in persons with mild chronic obstructive pulmonary disease who continue smoking. N Engl J Med. 1999;340(25):1948–53.
Burge PS, Calverley PMA, Jones PW, Spencer S, Anderson JA, Maslen TK. Randomised, double blind, placebo controlled study of fluticasone propionate in patients with moderate to severe chronic obstructive pulmonary disease: the ISOLDE trial. BMJ. 2000;320(7245):1297–303.
Vestbo J, Søorensen T, Lange P, Brix A, Torre P, Viskum K. Long-term effect of inhaled budesonide in mild and moderate chronic obstructive pulmonary disease: a randomised controlled trial. Lancet. 1999;353(9167):1819–23.
Calverley PM, Boonsawat W, Cseke Z, Zhong N, Peterson S, Olsson H. Maintenance therapy with budesonide and formoterol in chronicobstructive pulmonary disease. Eur Respir J. 2003;22(6):912–9.
Szafranski W, Cukier A, Ramirez A, Menga G, Sansores R, Nahabedian S, et al. Efficacy and safety of budesonide/formoterol in the management of chronic obstructive pulmonary disease. Eur Respir J. 2003;21(1):74–81.
Vestbo J, Anderson W, Coxson HO, Crim C, Dawber F, Edwards L, et al. Evaluation of COPD longitudinally to identify predictive surrogate end-points (ECLIPSE). Eur Respir J. 2008;31(4):869–73.
Papi A, Magnoni MS, Muzzio CC, Benso G, Rizzi A. Phenomenology of COPD: interpreting phenotypes with the ECLIPSE study. Monaldi Arch Chest Dis. 2016;83(1–2):721.
Hankinson JL, Odencrantz JR, Fedan KB. Spirometric reference values from a sample of the general US population. Am J Respir Crit Care Med. 1999;159(1):179–87.
Austin PC, Urbach DR. Using G-computation to estimate the effect of regionalization of surgical services on the absolute reduction in the occurrence of adverse patient outcomes. Med Care. 2013;51(9):797–805.
Feng Z, McLerran D, Grizzle J. A comparison of statistical methods for clustered data analysis with Gaussian error. Stat Med. 1996;15(16):1793–806.
Lee PN, Fry JS. Systematic review of the evidence relating FEV1 decline to giving up smoking. BMC Med. 2010;8:84.
Agustí AGN. Systemic effects of chronic obstructive pulmonary disease. Proc Am Thorac Soc. 2005;2(4):367–70.
Schols AM, Ferreira IM, Franssen FM, Gosker HR, Janssens W, Muscaritoli M, et al. Nutritional assessment and therapy in COPD: a European Respiratory Society statement. Eur Respir J. 2014;44(6):1504–20.
Schols AM, Broekhuizen R, Weling-Scheepers CA, Wouters EF. Body composition and mortality in chronic obstructive pulmonary disease. Am J Clin Nutr. 2005;82(1):53–9.
den Dekker HT, Jaddoe VWV, Reiss IK, de Jongste JC, Duijts L. Fetal and infant growth patterns and risk of lower lung function and asthma. The Generation R Study. Am J Respir Crit Care Med. 2017;197(2):183–92.
Ogawa E, Nakano Y, Ohara T, Muro S, Hirai T, Sato S, et al. Body mass index in male patients with COPD: correlation with low attenuation areas on CT. Thorax. 2009;64(1):20–5.
Hardin M, Foreman M, Dransfield MT, Hansel N, Han MK, Cho MH, et al. Sex-specific features of emphysema among current and former smokers with COPD. Eur Respir J. 2016;47(1):104–12.
Christou S, Chatziathanasiou T, Angeli S, Koullapis P, Stylianou F, Sznitman J, et al. Anatomical variability in the upper tracheobronchial tree: sex-based differences and implications for personalized inhalation therapies. J Appl Physiol. 2020. https://doi.org/10.1152/japplphysiol.00144.2020.
Dominelli PB, Ripoll JG, Cross TJ, Baker SE, Wiggins CC, Welch BT, et al. Sex differences in large conducting airway anatomy. J Appl Physiol. 2018;125(3):960–5.
Blettner M, Sauerbrei W, Schlehofer B, Scheuchenpflug T, Friedenreich C. Traditional reviews, meta-analyses and pooled analyses in epidemiology. Int J Epidemiol. 1999;28(1):1–9.
Divo MJ, Cabrera C, Casanova C, Pinto-Plata VM, Marin JM, de-Torres JP, et al. Comorbidity distribution, clinical expression and survival in COPD patients with different body mass index. Chronic Obstr Pulm Dis COPD Found. 2014;1(2):229–38.
The study was sponsored by the Canadian Institutes of Health Research (FDN-143226). The study sponsor played no role in the study design, collection, analysis and interpretation of data, nor in the writing of the report or the decision to submit the paper for publication.
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Observed rate of FEV1 decline across GOLD grades of severity, gender and smoking status. Table S2. Observed rate of FEV1 decline across gender, smoking status and BMI category in combined GOLD Grades 1 and 2. Figure S1. Scatter plot of BMI on individual rate of FEV1 decline in GOLD Grades 1 and 2.
About this article
Cite this article
Chen, W., Sadatsafavi, M., FitzGerald, J.M. et al. Gender modifies the effect of body mass index on lung function decline in mild-to-moderate COPD patients: a pooled analysis. Respir Res 22, 59 (2021). https://doi.org/10.1186/s12931-021-01656-5