Chronic obstructive pulmonary disease

Vitamin D status and chronic obstructive pulmonary disease risk: a prospective UK Biobank study

Abstract

Background Low vitamin D status has been linked to an increased risk for various inflammatory diseases. Conflicting results have been reported regarding chronic obstructive pulmonary disease (COPD). This study aims to investigate the associations of serum 25-hydroxyvitamin D (25(OH)D) concentrations with COPD risk and survival.

Methods We included 403 648 participants with serum 25(OH)D measurements and free of COPD at baseline from UK Biobank. Follow-up was until 30 September 2021. Multivariable-adjusted cox regression models were applied to estimate HRs and 95% CIs for the associations of season-standardised 25(OH)D concentrations with COPD risk and survival. The restricted cubic splines were used to assess dose–response relationship. Kaplan-Meier estimation was used to create graphs of the survival curves.

Results During a median follow-up of 12.3 (IQR: 11.4–13.2) years, 11 008 cases of COPD were recorded. We observed a non-linear inverse association between 25(OH)D concentrations and COPD risk. Compared with participants in the fourth quintile of 25(OH)D, those in the lowest quintile were associated with a 23% higher risk (HR, 1.23; 95% CI, 1.16 to 1.31). Stronger associations were observed for the risk in men and current smokers (Both p for interaction <0.05). In survival analyses, compared with the fourth quintile, cases in the lowest quintile had a 38% higher risk for overall death (HR, 1.38; 95% CI, 1.22 to 1.56).

Conclusion Our findings indicate that serum 25(OH)D concentrations are non-linearly negatively associated with incidence and mortality of COPD, suggesting a potential protective role of vitamin D in the pathogenesis of COPD.

What is already known on this topic

  • Low vitamin D status has been linked to an increased risk for various inflammatory diseases. Conflicting results have been reported regarding chronic obstructive pulmonary disease (COPD).

What this study adds

  • The decreasing risk of COPD appeared to be lowest at 55 nmol/L of 25-hydroxyvitamin D (25(OH)D). Lower pre-diagnostic 25(OH)D concentrations were associated with a significant decrease in overall and COPD-specific survival.

How this study might affect research, practice or policy

  • Whether lower concentrations of 25(OH)D are causal or contributory to COPD risk may spur future long-duration and large-scale randomised controlled trials.

Introduction

Chronic obstructive pulmonary disease (COPD) is a significant cause of disease burden worldwide.1 The pathological changes observed in COPD cases include chronic inflammation and structural changes due to repeated injury and repair.2 Abundant laboratory research has shown that vitamin D is essential in the acquired and innate immune systems.3 Meanwhile, observational studies have related vitamin D deficiency to a higher risk of infectious, autoimmune and inflammatory diseases.4 5 Although vitamin D status is traditionally known for its effect on skeletal health, existing studies suggest vitamin D status assessed by serum 25-hydroxyvitamin D (25(OH)D) has been related to the regulation of pathogenic mechanisms in COPD.5 6 However, epidemiological evidence regarding the associations of 25(OH)D concentrations with COPD incidence and survival remains inconclusive.

Cross-sectional studies have observed that vitamin D deficiency is highly prevalent among patients with COPD.7 Previous evidence has also associated lower 25(OH)D concentrations with worse lung function and a higher prevalence of COPD.8–10 As patients with COPD mostly exhibit increased skin ageing related to smoking behaviour and reduced sunlight exposure due to less outdoor activity, reductions in 25(OH)D concentrations are considered to be correlated with COPD development.11 Prospective cohort studies associating baseline 25(OH)D with declined lung function or incident COPD have yielded inconsistent results.12–14 Furthermore, results from the general population in the Copenhagen City Heart Study indicated an association between lower 25(OH)D concentrations and higher COPD risk,15 but another Danish study using national register data found no such association.16

Although vitamin D deficiency has been associated with higher mortality in the general population,17 18 the effect on patients with COPD remains unclear. Some studies based on patients with COPD found no relationship between 25(OH)D and mortality.19 20 However, post-diagnostic blood vitamin D status in these studies could be primarily modified by COPD induced systemic inflammatory responses.11 Thus, it seems more stable to investigate the association using pre-diagnostic 25(OH)D in patients with COPD.

Therefore, leveraging the serum vitamin D measurements in the UK Biobank, we aimed to assess the associations of 25(OH)D concentrations with COPD incidence in 403 648 participants and survival in 11 008 COPD cases.

Methods

Patient and public involvement

This study was conducted on UK Biobank database and patients or public were not involved in the development of study design or recruitment of participants.

Study population

UK Biobank is a large prospective cohort study that recruited over 500 000 participants aged 40–69 years in 22 assessment centres across the UK from 2006 to 2010. At the baseline recruitment visit, individuals who signed consent completed a questionnaire assessment and underwent physical measurements. Blood samples were collected from all participants at recruitment and from a subset of participants with a repeat assessment. All participants were registered with the UK National Health Service (NHS). The present study was conducted under the UK Biobank application number 84525.

We excluded participants with airway obstruction (n=49 733), including prevalent COPD at recruitment (n=3516), self-reported emphysema/chronic bronchitis (n=8349), forced expiratory volume at 1 s/forced vital capacity (FEV1/FVC) <lower limit of normal (calculated according to the Global Lung Function Initiative 2012 criteria21) (n=41 779) and participants with unavailable measurements on 25(OH)D (n=49 032). Totally, 403 648 participants were included in the study (online supplemental figure S1).

Assessment of serum 25(OH)D

As part of the UKB Biomarker Project, serum 25(OH)D concentrations were measured by a direct competitive chemiluminescent immunoassay (DiaSorin Liaison XL), with a detection range of 10–375 nmol/L. The mean coefficients of variation were 5.04%, 5.39% and 6.14% for low, medium and high internal quality control levels of 25(OH)D. To further verify the accuracy, the assay of serum 25(OH)D was registered with an external quality assurance (EQA) scheme (RIQAS Immunoassay Specialty 1). According to EQA results, 100% of participated distributions (n=108) were considered good or acceptable.

Ascertainment of COPD

Incident COPD data were available through linkage to the NHS register until 30 September 2021 for England, 31 July 2021 for Scotland and 28 February 2018 for Wales. Cause and date of death after COPD diagnosis were obtained from death certificates until 31 October 2021 for Scotland and 30 September 2021 for England and Wales. The registers are coded according to the 10th revision of the International Classification of Diseases (ICD-10). We used ICD-10 codes to define COPD as J40–44.

Ascertainment of covariates

Covariates were collected from baseline survey and selected according to possible effects on serum 25(OH)D and COPD, including sociodemographic characteristics (age, sex, ethnicity, education and Townsend Deprivation Index), lifestyles (smoking status, alcohol consumption, physical activity and dietary intake), environmental pollution (passive smoking and fine particulate matter, PM2.5), body mass index (BMI) (defined as the weight (kg) divided by height (m) squared), medical history (prevalent asthma and family history of respiratory system diseases) and female-specific factors (menopausal status and hormone replacement therapy). Physical activity was measured as metabolic equivalent task hours per week for moderate activity. PM2.5 estimated for the year 2010 was modelled for each address using a land use regression model. Asthma at baseline was defined using self-report and hospital admission data (ICD-10 codes J45–46). Further details on covariate measurements were recorded in the UK Biobank online protocol.22

Statistical analysis

We calculated season-standardised 25(OH)D concentrations by fitting a cosinor model.23 Baseline characteristics were described among participants grouped by quintiles of 25(OH)D concentrations. The intraclass correlation coefficient (ICC) in two measurements of 25(OH)D between the baseline and the repeated assessments was computed.

In the analysis of COPD risk, participants contributed person-time from the date of recruitment until the date of the initial diagnosis of COPD, lost to follow-up, death or the end of the study, whichever occurred first. Cox proportional hazard models were used to estimate HRs and 95% CIs for the associations of 25(OH)D with COPD risk among total participants, never-smokers and smokers. The proportional hazards assumption was tested using Schoenfeld residuals and no evidence of non-proportionality was detected. We used multivariate cubic regression splines with four knots to visually explore non-linear associations between serum 25(OH)D concentrations and the risk of COPD.

We classified serum 25(OH)D into quintiles and estimated the HRs using the fourth quintile reflected the typical vitamin D level in UK population as the reference in the primary analysis. We considered three analytic models. Model 1 was adjusted for age, sex, ethnicity and assessment centre. Model 2 was additionally adjusted for education degree, Townsend deprivation index, BMI, smoking status, alcohol drinking, physical activity, oily fish intake, regular use of vitamin D supplements, prevalent asthma, family history of respiratory diseases and female-specific factors. Model 3 was further adjusted for passive smoking, PM2.5 and occupations at risk of COPD. For covariates with missing data, we used the missing indicator method for categorical variables and used the median imputation by gender for continuous variables.

In addition, stratified analyses were conducted by age (≤60, >60 years), sex (male, female), BMI category (<25, 25–30, ≥30 kg/m2), smoking status (never, former and current smokers), passive smoking (never, ever), physical activity (≤median, >median) and use of vitamin D supplements (No, Yes) in the model 3. A likelihood ratio test was performed to examine possible effect modification between 25(OH)D concentrations and the stratification variables, comparing the models with and without interaction terms. Sensitivity analyses were conducted by excluding COPD cases diagnosed within 2 years after baseline assessment, excluding the participants with overall poor self-rated health in baseline questionnaires, and excluding the participants with use of vitamin D supplements.

In survival analyses, we calculated overall and COPD-specific survival time from the date of first diagnosis of COPD to the date of death or end of follow-up. Kaplan-Meier estimation was used to create the graphs of survival curves and the log-rank test was used to compare curves from different quintiles of 25(OH)D. We also use multivariate cox regression to estimate HRs and 95% CIs for overall and COPD-specific survival. The survival analyses were adjusted for the same covariates as the previous analyses.

All analyses were performed using SAS V.9.4 (SAS Institute). Statistical tests were all two-sided and p<0.05 was considered statistically significant.

Results

Among 403 648 participants (185 840 men, 217 808 women) followed for a median of 12.3 years (IQR, 11.4–13.2 years), we documented 11 008 COPD cases (5956 men and 5052 women). Of these COPD cases, 2773 deaths (225 from COPD) occurred within a median survival time of 3.8 years (IQR, 1.6–6.6 years). The ICC in serum 25(OH)D between the initial and repeated assessment was 0.59 (95% CI, 0.58 to 0.60).

The baseline characteristics of participants by quintiles of serum 25(OH)D concentrations are shown in table 1. Participants with higher 25(OH)D had higher economic standards, proportion of oily fish intake, and tended to use vitamin D supplements. The proportion of current smokers and prevalent asthma cases was lower among these participants.

Table 1
|
Baseline characteristics of participants according to quintiles of serum 25(OH)D concentrations

Figure 1 shows a non-linear inverse relationship between 25(OH)D concentrations and COPD incidence among total participants, never-smokers and smokers (all p values for non-linearity <0.0001). We visually observed that the risk of COPD decreased with increasing 25(OH)D concentrations and was lowest at 55 nmol/L in total participants. Table 2 shows the associations between 25(OH)D and COPD risk among total participants, never-smokers and smokers. In the fully adjusted model, compared with the fourth quintile, the lowest quintile groups increase a 23% risk (HR, 1.23; 95% CI, 1.16 to 1.31) in total participants, a 25% risk (HR, 1.25; 95% CI, 1.10 to 1.42) in never-smokers and a 23% risk (HR, 1.23; 95% CI, 1.15 to 1.32) in smokers.

Figure 1
Figure 1

Non-linear inverse associations of serum 25(OH)D concentrations with risk of COPD among total participants (A), never-smokers (B) and smokers (C). The associations were examined by multivariate cox regression models based on restricted cubic splines. Participants with 25(OH)D concentrations above the 99.9th percentile were excluded in figure 2(A). Never-smokers and smokers with 25(OH)D concentrations above the 99.9th percentile were excluded in figure 2(B) and (C). Solid line represents estimates of HRs and dashed lines represent 95% CIs. 25(OH)D, 25-hydroxyvitamin D; COPD, chronic obstructive pulmonary disease.

Figure 2
Figure 2

Survival curves for overall survival (A) and COPD-specific survival (B) according to quintiles of serum 25(OH)D concentrations. 25(OH)D, 25-hydroxyvitamin D; COPD, chronic obstructive pulmonary disease.

Table 2
|
Association between serum 25(OH)D concentrations and COPD risk among total participants, never-smokers and smokers

The stratified analysis (online supplemental figure S2) showed that the above associations for COPD generally remained consistent across subgroups. Effect modification by sex and smoking status was observed for COPD risk (both p values for interaction <0.05), and the HRs were stronger in men and current smokers.

Sensitivity analyses showed the results were unchanged after excluding 822 COPD cases diagnosed within 2 years prior to follow-up (online supplemental table S1), excluding 15 695 participants who self-rated overall poor health at baseline assessments (online supplemental table S2) or excluding 98 637 participants with use of vitamin D supplements (online supplemental table S3).

In survival analyses, we observed an association of higher pre-diagnostic 25(OH)D concentrations with a significant increase in overall and COPD-specific survival as analysed by a log-rank test (both p values <0.05) (figure 2). In the fully adjusted models, compared with the fourth quintile of 25(OH)D concentrations, cases in the lowest quintile had a 38% higher risk for overall death (HR, 1.38; 95% CI, 1.22 to 1.56) and 57% higher risk for COPD-specific death (HR, 1.57; 95% CI, 1.03 to 2.40) (table 3).

Table 3
|
Associations of serum 25(OH)D concentrations with all-cause or COPD-specific death in COPD cases

Discussion

In this large prospective study based on European population, we observed non-linear inverse associations between serum 25(OH)D concentrations and the risk of COPD. The decreasing risk of COPD appeared to be lowest at 55 nmol/L of 25(OH)D. The robustness of the associations was supported by comprehensive stratified and sensitivity analyses. In addition, we also found that lower pre-diagnostic 25(OH)D concentrations were associated with a significant decrease in overall and COPD-specific survival. Our findings imply that vitamin D might play a role in progression of COPD.

Although a growing body of evidence reported that vitamin D deficiency is associated with multiple diseases, the association between vitamin D status and COPD risk remained inconsistent.24 Cross-sectional and case–control studies indicated that vitamin D deficiency was highly prevalent in patients with COPD.25 26 Of note, the prior studies were conducted in patients with COPD but the results were inconclusive in generally healthy individuals. An initial analysis of the Third National Health and Nutrition Examination Survey (NHANES III) spirometry data found a strong association between 25(OH)D concentrations and lung function assessed by FEV1 and FVC, but no association could be found with airway obstruction.9 Several prospective cohort studies reported that vitamin D deficiency or insufficiency was associated with decreased lung function or increased risk of COPD.12 13 27 However, some other studies found no associations.14 28 Such inconsistency may be attributed to smaller sample size, short follow-up period and residual confounding caused by the incomprehensive adjustment. In addition, multiple randomised controlled trials (RCTs) had demonstrated that vitamin D supplementation reduced the rate of worsening moderate or severe COPD exacerbations in patients with low baseline 25(OH)D concentrations.29 30 Furthermore, a meta-analysis based on vitamin D binding protein gene polymorphisms studies also suggested that vitamin D status was associated with COPD risk.31 In the present study, we found that lower 25(OH)D concentrations could be related to higher COPD risk, providing further evidence for the association between 25(OH)D and COPD.

It is well-established that smoking is the most commonly encountered risk factor for COPD.32 Our findings indicated that 25(OH)D concentrations were inversely associated with COPD risk in both smokers and never-smokers. Extensive studies had reported smokers had a greater decline in FEV1, a higher incidence of abnormal lung function and a higher mortality rate from COPD.33 34 A nested case–control study found a relationship between vitamin D deficiency and COPD in smokers, but lacked evidence in never-smokers.35 Moreover, a cross-sectional study observed individuals exposed to secondhand smoking have increased risks of COPD.36 The current study showed that the association remained unchanged in never-smokers after adjustment for passive smoking. Furthermore, we also found that lower 25(OH)D concentrations might have a more significant effect on elevating the risk of COPD in current smokers. A meta-analysis including 24 studies with 11 340 participants suggested that smokers were likely to have lower 25(OH)D concentrations.37 Future clinical or laboratory studies are needed to explore the underlying mechanism.

Moreover, there is still debate on the optimal range of 25(OH)D. US Institute of Medicine considered 25(OH)D concentrations of 50 nmol/L (20 ng/mL) or higher to be adequate for bone and overall health,38 while the Endocrine Society stated that of 75 nmol/L (30 ng/mL) or higher provide increased health benefits.39 Our previous finding indicated that 45–60 nmol/L of 25(OH)D might be a potential interventional target for premature death.18 However, the direct evidence regarding optimal 25(OH)D on respiratory diseases was scarce. A recent Mendelian randomisation study suggested that 50–75 nmol/L of 25(OH)D may contribute to lower mortality of respiratory diseases.40 In our analysis, optimal concentrations of 50–60 nmol/L were observed to be associated with lower COPD risk.

For COPD survival, our results are consistent with previous studies indicating a negative association between pre-diagnostic 25(OH)D concentrations and COPD mortality.41 42 However, other longitudinal studies found no relationships between post-diagnostic 25(OH)D and COPD survival.12 19 These studies were restricted by the small sample size (n<500), and post-diagnostic 25(OH)D concentrations might lead to reverse causation.

Our analysis has several strengths. First, the large sample size and the morbidity and mortality data based on the NHS records provided sufficient power to detect associations. Second, serum 25(OH)D was assessed in a central laboratory using standardised, validated blood biochemistry methods with strict quality control, thereby minimising any measurement errors. Third, we comprehensively adjusted for broad demographics, lifestyle, and medical history. Several limitations also need to be acknowledged. First, reverse causality could not be excluded. However, all participants with baseline COPD were excluded from the analysis, and sensitivity analyses support the robustness of our findings. Second, a single measurement of 25(OH)D at recruitment might not represent long-term exposure concentrations. However, existing studies suggested a single measurement could adequately reflect vitamin D status over an extended period,43 and the moderate ICC revealed that time-dependent variation in vitamin D was unlikely to produce a substantial influence on our findings. Third, blood samples were collected from all participants at recruitment to minimise selection bias. However, 20% of participants were excluded, with 10% excluded due to missing available measurement for 25(OH)D and the remaining 10% excluded due to airway obstruction at baseline. Finally, since the majority of participants in this study were from the UK, the findings from this study might not apply to other populations.

In conclusion, we observed serum 25(OH)D concentrations were non-linearly associated with COPD risk in both smokers and never-smokers. Our study also indicated that higher concentrations of 25(OH)D were associated with improved survival of COPD. Whether lower concentrations of 25(OH)D are causal or contributory to COPD risk may spur future long-duration and large-scale RCTs.