Association of Serum 25-Hydroxyvitamin D Concentrations with All-Cause and Cause-Specific Mortality Among Adult Patients with Existing Cardiovascular Disease

Abstract

Background:

Vitamin D insufficiency and deficiency are common in patients with cardiovascular disease (CVD). We aimed to prospectively examine the associations of serum 25-hydroxyvitamin D [25(OH)D] concentrations with all-cause and cause-specific mortality among adult patients with existing CVD.

Methods:

We included 37,080 patients with CVD from the UK Biobank study, a prospective cohort of half a million participants aged 40–69 years. We defined patients with CVD as those who suffered coronary heart disease, atrial fibrillation, heart failure, or stroke. The associations of serum 25(OH)D concentration with all-cause and cause-specific mortality were examined by using multivariable Cox regression models and competing risk analyses.

Results:

Among 37,080 patients with CVD at baseline, 57.5% were subjected to vitamin D deficiency (i.e., 25[OH]D <50 nmol/L). During a median follow-up of 8.9 years, 3920 total deaths occurred, including 1411 deaths from CVD, 1497 deaths from cancer, 368 deaths from respiratory disease, and 644 other-cause deaths. We observed non-linear inverse associations for all-cause, cancer, respiratory disease, and other-cause mortality (P-nonlinearity <0.001) and approximately linear inverse associations for CVD mortality (P-nonlinearity = 0.069). Among CVD patients with vitamin D deficiency, per 10 nmol/L increment in serum 25(OH)D concentrations was associated with an 11% reduced risk for all-cause mortality and 10% reduced risk for CVD mortality.

Conclusion:

Among patients with existing CVD, increasing levels in serum 25(OH)D were independently associated with a decreased risk of all-cause and cause-specific mortality. These findings suggest that elevated serum 25(OH)D concentration benefits CVD patients with vitamin D deficiency.

Background

Cardiovascular disease (CVD) remains the leading cause of mortality and poses significant health burdens globally.¹ The past few decades have witnessed an escalating research interest in the potential role of vitamin D in CVD prevention.^{2, 3} Observational studies have linked low 25‑hydroxyvitamin D (25[OH]D) serum levels with increased risk of subsequent CVD.^3–7 Despite the accumulating mechanism- and population-based evidence of vitamin D on preventing CVD, recent randomized controlled trials (RCTs) failed to establish the cardiovascular benefits of vitamin D supplementation in the general population or those without a previous history of CVD.^8–10 Although results from these trials are not encouraging, vitamin D status remains an important worldwide public-health concern.¹¹ Vitamin D insufficiency and deficiency are common in patients with CVD.^{12, 13} Existing evidence from RCTs indicates that vitamin D supplementation exetys beneficial effects on coronary artery disease (CAD)¹⁴, left ventricular structure and function in patients with chronic heart failure (CHF)¹⁵, and inflammatory milieu in patients with CHF.¹⁶

Evidence from observational studies is limited and inconsistent in this regard.^17–24 Prior cohort studies among patients with CVD¹⁷, HF^18–20, suspected coronary artery disease²¹, or suspected stable angina pectoris²² found that plasma 25(OH)D concentrations are inversely associated with all-cause and cardiovascular mortality. However, two cohort studies did not find a significant association between vitamin D levels and all-cause mortality or secondary cardiovascular event incidence.^{23, 24} Meanwhile, previous studies are subjected to several limitations, including a relatively small sample size, ignoring dietary factors, physical activity, season of vitamin D assessment, and CVD duration and other confounding factors, which might limit the interpretation of the results. Improved understanding of the association between vitamin D status and mortality may shed new light on the potential of vitamin supplementation to provide clinical benefits for patients with CVD in further RCTs. To fill the knowledge gap, we aimed to prospectively examine the associations of serum 25(OH)D concentrations with all-cause and cause-specific mortality among adult patients with existing CVD in a large cohort.

Methods

Study population

The UK Biobank is a prospective cohort of half a million participants, aged 40–69 years, recruited between 2006 and 2010 throughout the UK.²⁵ These participants completed extensive touch screen questionnaires, physical and functional measures, and collection of biological samples. Follow-up was conducted through linkages to routinely available national datasets. Ethical approval was obtained by the National Health Service National Research Ethics Service (11/NW/0382) and renewed by the North West–Haydock Research Ethics Committee (16/NW/0274). All participants provided informed written consent. The background information about UK Biobank and details is available on the website http://www.ukbiobank.ac.uk.

Definitions of baseline CVD were constructed using the 9th revision of the International Statistical Classification of Diseases (ICD-9), ICD-10, and Office of Population, Censuses and Surveys-4 codes as well as self-reported data fields with choice-, disease- or procedure-specific codes. We defined patients with CVD as those who suffered coronary heart disease (CHD), atrial fibrillation (AF), HF, or stroke. Detailed definition is described in Additional file 1: Table S1. In the present study, we excluded participants without CVD at baseline (n=461,643), leaving 40,850 participants with pre-existing CVD. Participants with missing data on serum 25(OH)D (n=3760) and pregnancy (n=10) were further excluded from analyses. Finally, 37,080 patients with CVD were included (Additional file 1: Figure S1).

Assessment of vitamin D status

Vitamin D status was classified according to the serum levels of 25(OH)D. The Endocrine Society Clinical Practice Guidelines define vitamin D severe deficiency as serum 25(OH)D levels <25.0 nmol/L, moderate deficiency as 25.0–49.9 nmol/L, insufficient as 50.0–74.9 nmol/L, and sufficient as ≥75.0 nmol/L.²⁶ In the UK biobank, analysis of serum 25(OH)D concentrations utilized immunoassay analyzers (DiaSorin Liaison XL Analyzer, made in Diasorin S.p.A) by a direct competitive chemiluminescent immunoassay method, with an analytical range of 10–375 nmol/L.

Assessment of outcomes

In the UK Biobank study, mortality data of each participant were obtained by linkages to National Health Service (NHS) datasets, including the NHS Digital (for England and Wales) and the NHS Central Register (for Scotland). The date of death and the causes of death were provided and coded using the ICD-10 system. In addition to the all-cause mortality, primary cause of death was extracted from the UK Biobank Cause of Death Registry, including cardiovascular mortality (ICD-10: I00-I99), cancer mortality (ICD-10: C00-C97), respiratory disease mortality (ICD-10: J00-J99), and other-cause mortality (mortality excluding cardiovascular, cancer, and respiratory disease).

Covariates

Information on age, sex, ethnicity, education, household income, and lifestyle behaviors was acquired using touch screen questionnaires. Participants were labelled as never, former, and current smokers according to the summarized smoking status (Field ID: 20116). Daily average alcohol consumption was described depending on the drinking frequency and the number of drink-equivalents/day. Physical activity was classified in accordance with 2018 Physical Activity Guidelines for Americans as inactive (those with no documented leisure time physical activity), insufficient (<150 min/week of moderate activity and <75 min/week of vigorous activity), and active (≥150 min/week of moderate activity and/or ≥75 min/week of vigorous activity).²⁷ We constructed a healthy diet score with reference to the dietary priorities for cardiometabolic health recommended by American Heart Association.²⁸ Definitions of each component of a healthy diet score are shown in Additional file 1: Table S2. The scale of the healthy diet score ranged from 0 to 10, and a higher score equates to a much healthier dietary pattern. Adherence to a healthy diet was defined as participants who had at least 5 scores of healthy diet components. Socioeconomic status was indicated by Townsend deprivation index scores, and higher Townsend scores equate to higher levels of socioeconomic deprivation.²⁹ Details regarding health conditions and drug use were ascertained by touch screen questionnaires, face-to-face interviews, and linkage to electronic health records. A trained nurse measured the height, weight, and blood pressure during the initial assessment. Serum was collected through venipuncture, and biomarkers including glycated hemoglobinA1c (HbA1c), lipids, and C-reactive protein (CRP) were measured. The estimated glomerular filtration rate (creatinine-cystatin C equation, eGFRcr-cys) was calculated from serum creatinine and cystatin C.³⁰

Statistical analysis

Multivariable Cox regression models to the estimate hazard ratios (HRs) and 95% confidence intervals (CIs) were used for the associations of the serum 25(OH)D concentration with all-cause mortality. Competing risk analyses were conducted using the cause-specific hazard function model to estimate hazards for CVD, cancer, respiratory disease, and other-cause mortality.^{31, 32} The time to events was calculated from the date of the blood sample collection to the death or the censoring date (14 February 2018), whichever came first. Participants with severe deficiency of vitamin D (serum 25[OH]D <25.0 nmol/L) were selected as the reference group. Models were successively adjusted for age (continuous), sex (male, female), ethnicity (White, mixed, Asian, Black, Chinese, others), education (college or university, vocational qualification, upper secondary, lower secondary, others), Townsend deprivation index (in quintiles), household income (<18,000; 18,000–30,999; 31,000–51,999; 52,000–100,000; >100000£), smoking status (never smoker, former smoker, current smoker), alcohol consumption (0, 0.1–4.9, 5.0–14.9, 15.0–19.9, 20.0–29.9, ≥30.0 g/day), physical activity (inactive, insufficient, active), healthy diet score (in quintiles), BMI (<18.5, 18.5–22.9, 23.0–24.9, 25.0–29.9, 30.0–34.9, ≥35.0 kg/m²), eGFR_cr-cys (<30.0, 30.0–60.0, 60.0–90.0, ≥90.0 mL min^-1 per 1.73 m²), CRP (in quintiles), anti-hypertensive medication use, cholesterol lowering medication use, diabetes medication use (none, only oral medication, only insulin, or insulin and oral medication), history of cancer, diabetes, hypertension, and CVD duration (<1.0, 1.0–4.9, 5.0–9.9, ≥10.0 years). Missing values of covariates were treated as dummy variables. The dose–response curves presenting the hazard of serum 25(OH)D were fitted by using the restricted cubic spline model with four knots (rms, hmisc, lattice, and survival packages in R software).

Several sensitivity analyses were performed to test the robustness of our results. We performed subgroup analyses across age, sex, BMI, smoking status, physical activity, dietary supplement use, antihypertensive treatment, cholesterol lowering medication, and CVD duration. The joint test was used to obtain a P value for interaction for examining the statistical significance of the difference between subgroups. Considering the influence of seasonal fluctuations and time spent outdoors on circulating 25(OH)D concentrations, we further adjusted for the month of blood collection (January through December, categorical), time spent outdoors in summer (continuous), and time spent outdoors in winter (continuous) in Cox regression models. We further adjusted for blood pressure, HbA1c, lipids, and dietary supplement use, including vitamin D supplements, multivitamin supplements, mineral supplesments, fish oil, and glucosamine. We also excluded patients with diagnosed thyroid and parathyroid diseases. In view of the unavailable data on blood parathyroid hormone levels in UK Biobank, we adjusted for the serum levels of calcium and phosphate.³³ We also considered energy intake and dietary vitamin D intake of 12,505 patients who completed 24h dietary recalls between April 2009 and June 2012. To minimize the potential reverse causation bias, we excluded patients who died within 2 years and re-examined the association between serum 25(OH)D and mortality.

Statistical analyses were performed between 1 November, 2020 and 16 March, 2021. SAS version 9.4 (SAS Institute, USA) and R software (The R Foundation, http://www.r-project.org, version 4.0.2) were utilized for analyses and plotting with a two-sided significance threshold of P <0.05.

Results

Baseline characteristics

A total of 37,080 patients with CVD (mean age, 61.4 years [SD, 6.4 years]; 12,668 [34.2%] females) were included in the present analysis. Table 1 shows the baseline characteristics according to serum 25(OH)D levels. Among these patients, 5772 (15.6%); 15,559 (42.0%); 11,451 (30.9%); and 4298 (11.6%) were in the vitamin D status of severe deficiency, moderate deficiency, insufficient, and sufficient, respectively. Compared with patients with vitamin D deficiency, those with higher serum 25(OH)D concentrations were more likely to be older, male, non-smoker, and to have lower levels of socioeconomic deprivation, drink more alcohol, adhere to a healthy diet pattern, exercise more, use more dietary supplements, have a lower BMI, and have lower prevalence of hypertension and diabetes but a higher prevalence of cancer.

The adjusted means (95% CIs) of cardiometabolic markers were presented according to serum 25(OH)D concentrations (Table 2). Higher levels of serum 25(OH)D were significantly associated with lower levels of triglycerides, total cholesterol, low-density lipoprotein direct, lipoprotein(a), apolipoprotein-A, apolipoprotein-B, systolic blood pressure, diastolic blood pressure, HbA1c, CRP, and eGFRcr-cys (all P_trend <0.05). However, levels of high-density lipoprotein cholesterol did not change significantly across serum 25(OH)D categories (P_trend=0.75).

Serum 25(OH)D and mortality

During 320,387 person-years of follow-up (median 8.9 years), 3920 total deaths occurred, including 1411 deaths from CVD, 1497 deaths from cancer, 368 deaths from respiratory disease, and 644 other-cause deaths. In Cox regression analyses, compared with patients in severe deficiency of vitamin D (serum 25[OH]D <25 nmol/L), those with moderate deficiency, insufficient, or sufficient vitamin D status showed a decreased hazard for all-cause and cause-specific mortality (Table 3). These associations remained robust after stepwise adjustment for confounders. In the fully adjusted model, the HRs and 95% CIs from the lowest to the highest serum 25(OH)D categories (<25.0, 25.0–49.9, 50.0–74.9, and ≥75.0 nmol/L) were 1.00 (reference), 0.80 (0.74, 0.87), 0.69 (0.63, 0.76), and 0.69 (0.61, 0.78), respectively, for all-cause mortality; 1.00 (reference), 0.79 (0.68, 0.90), 0.69 (0.59, 0.81), and 0.65 (0.53, 0.81), respectively, for CVD mortality; 1.00 (reference), 0.89 (0.77, 1.03), 0.77 (0.66, 0.90), and 0.80 (0.66, 0.98), respectively, for cancer mortality; 1.00 (reference), 0.73 (0.56, 0.94), 0.49 (0.36, 0.67), and 0.58 (0.38, 0.88), respectively, for respiratory disease mortality; and 1.00 (reference), 0.73 (0.62, 0.90), 0.62 (0.51, 0.87), and 0.64 (0.46, 0.88), respectively, for other-cause mortality (Model 3 in Table 3).

We further fitted smoothing splines to present the dose–response relationship between serum 25(OH)D concentrations and risk of mortality. All-cause and CVD mortality decreased with increasing serum 25(OH)D concentrations and then reached a plateau at around 50 nmol/L 25(OH)D (Figure 1). The test of nonlinearity was statistically significant (P-nonlinearity <.001) for all-cause mortality but not significant for CVD mortality (P-nonlinearity = .069). Similar non-linear trends were observed for cancer, respiratory disease, and other-cause mortality (Additional file 1: Figure S2). Among patients with vitamin D deficiency (serum 25(OH)D <50.0 nmol/L), per 10 nmol/L increment in serum 25(OH)D concentrations was associated with a 11% (HR, 0.89 [95%CI: 0.85, 0.92]) reduced risk for all-cause mortality, 10% (HR, 0.90 [95%CI: 0.85, 0.96]) reduced risk for CVD mortality, 9% (HR, 0.91 [95%CI: 0.85, 0.97]) reduced risk for cancer mortality, 15% (HR, 0.85 [95%CI: 0.75, 0.96]) reduced risk for respiratory disease mortality, and a 15% (HR, 0.85 [95%CI: 0.77, 0.93]) reduced risk for other-cause mortality (Table 3). However, among patients with 25(OH)D concentrations ≥50.0 nmol/L, increment in 25(OH)D concentrations was not significantly associated with a decreased risk for all-cause and cause-specific mortality.

We also investigated the associations of vitamin D status with all-cause, CVD, cancer, respiratory disease, and other-cause mortality of various CVD types, including CHD, AF, HF, and stroke (Table 4). Among patients with stroke, the HRs and 95% CIs from the lowest to the highest serum 25(OH)D categories were 1.00 (reference), 0.77 (0.65, 0.90), 0.65 (0.54, 0.79), and 0.55 (0.41, 0.72), respectively, for all-cause mortality and 1.00 (reference), 0.80 (0.60, 1.07), 0.64 (0.45, 0.91), and 0.36 (0.21, 0.64), respectively, for CVD mortality.

Sensitivity analyses

In stratified analyses, the associations of serum 25(OH)D concentration with all-cause and CVD mortality were robust across the strata of age, sex, BMI, smoking status, physical activity, dietary supplement use, antihypertensive treatment, cholesterol lowering medication, and CVD duration (Additional file 1: Table S3). We observed that serum 25(OH)D concentrations increased with prolonged time spent outdoors, and the seasonal changes in 25(OH)D concentration showed a peak during the summer months and a trough during the winter months (Additional file 1: Figure S3). Additional adjustment of the confounders of the month of blood collection and time spend outdoors did not change the associations between serum 25(OH)D and mortality risk (Additional file 1:Table S4, Model 1). Further sensitivity analyses with additional adjustment of dietary supplement use, mean arterial pressure, HbA1c, triglycerides, and high-density lipoprotein cholesterol produced similar results (Additional file 1:Table S4, Models 2 and 3). These associations did not change appreciably with further adjustment of serum calcium and phosphate after excluding participants with a history of thyroid or parathyroid diseases (Additional file 1:Table S5). In addition, we excluded participants who died within 2 years of follow-up and re-examined the associations. Results did not alter the significance of the associations between 25(OH)D concentration and outcomes (Additional file 1:Table S6 ). Similar results were also obtained with further adjustment of energy intake and dietary vitamin D intake (Additional file 1:Table S7).

Discussion

In this large prospective cohort study, near 60% of patients with CVD were subjected to vitamin D deficiency (serum 25[OH]D < 50 nmol/L). Among patients with existing CVD, increasing levels in serum 25(OH)D were independently associated with a decreased risk of all-cause and cause-specific mortality. Such associations presented non-linear trends for all-cause, cancer, respiratory disease, and other-cause mortality, and a linear trend for cardiovascular mortality. Our findings underscore the importance of optimizing vitamin D status for patients with existing CVD, especially for those subjected to vitamin D deficiency.

Although numerous observational studies have examined the association of vitamin D with all-cause and cause-specific mortality, most studies focused on the general population and deliberately excluded patients with known CVD.³⁴ Existing evidence for the long-term association between vitamin D status and adverse outcomes in CVD patients is inconsistent and insufficient.^17–24 One cohort study including 1125 German patients with stable CHD found no significant association of serum 25 (OH)D levels with secondary cardiovascular event incidence and all-cause mortality.²⁴ In another longitudinal study with 946 stable CHD patients in the San Francisco Bay Area, Welles et al. found that 25(OH)D levels under 50 nmol/L remain independently associated with cardiovascular events, but the independent association is no longer present (adjusted HR, 1.11 [95% CI: 0.85, 1.44]) after further adjustment for potential biological mediators.²³ Three small-sample studies found that vitamin D deficiency is associated with an increased risk of ischemic heart disease events among 244 patients with prior CVD (including ischemic heart disease, peripheral vascular disease, and stroke).¹⁷ In addition, a low 25(OH)D concentration is associated with a poor prognosis in patients with HF (n = 148 and n = 548).^{19, 20} A clinic-based study in Israel found that vitamin D deficiency is more prevalent in HF patients (n = 3009) than in the control group (n = 46,825).¹⁸ HF patients with severe deficiency of vitamin D (25[OH]D < 25.0 nmol/L) have higher risk of mortality than those with serum 25(OH)D ≥ 75.0 nmol/L (adjusted HR, 1.61 [95% CI: 1.08, 2.41]). The authors observed that vitamin D supplementation is associated with reduced mortality in HF patients with an adjusted HR, 0.68 [95% CI: 0.54–0.85]). However, the clinic-based design and relatively short follow-up duration (median 518 days) might limit the generalizability of these findings to the broader population of CVD patients and long-term outcomes. Similar limitations also remained in the hospital-based single-center study comprising 3316 patients of White ethnicity referred for coronary angiography.²¹ In another study among 4114 white patients suspected of having stable angina pectoris, Eirik et al. found that plasma 25(OH)D concentrations are nonlinearly (U-shaped) associated with all-cause mortality; near 70 nmol/L 25(OH)D corresponds to the lowest mortality risk, and their analysis suggested increased all-cause mortality at concentrations > 100 nmol/L.²² However, our findings were not consistent with this trend. All-cause mortality decreased with increasing serum 25(OH)D concentration and then reached a plateau at around 50 nmol/L 25(OH)D. A small amount could potentially decrease the risk of cardiovascular mortality with higher 25(OH)D concentrations (> 100 nmol/L) in our analyses. This continuing downward trend of all-cause and cardiovascular mortality was more pronounced among patients with stroke. Further clinical trials are warranted to determine the benefits of vitamin D supplementation among patients with existing CVD and subtypes, and the threshold range of 25(OH)D is applicable to vitamin D supplementation.

To the best of our knowledge, this dataset is the largest to examine the associations of serum 25(OH)D concentrations with all-cause and cause-specific mortality among adult patients with existing CVD. Compared with previous studies, the present study considered much more comprehensive confounders, conducted numerous sensitivity analyses, and included various CVD subtypes. Our results suggested that CVD patients with vitamin D deficiency (25[OH]D < 50 nmol/L) were more likely to benefit from optimizing vitamin D status than those without vitamin D deficiency. An interesting finding was that the analysis of patients with stroke showed impressive reductions in CVD mortality for the group with sufficient vitamin D (25[OH]D ≥ 75.0 nmol/L), with an adjusted HR of 0.36 (95% CI: 0.21, 0.64). This finding implied that stroke patients might derive additional cardiovascular benefits from maintaining more adequate vitamin D status than patients with heart disease. Although insufficient evidence has been proposed in this regard, our findings provided a novel perspective of vitamin D status and cardiovascular morbidity and mortality in patients with stroke. Nowadays, several RCTs are underway to assess the effect of vitamin D supplementation on patients with heart failure (i.e., NCT03416361, NCT03289637), idiopathic cardiomyopathies (i.e., NCT02517814), and coronary artery disease (i.e., NCT01570309, NCT02996721). In the future, these clinical trials may further demonstrate our findings.

The strengths of our study include the large number of participants and the wide range of serum 25(OH)D levels. The high-quality data from UK Biobank favor the definition and classification of various CVD subtypes and allow us to examine the long-term association of vitamin D status with all-cause and cause-specific mortality. Associations were robust after adjustment for a wide range of potential confounding factors, including socioeconomic status, dietary and lifestyle factors, comorbidities, and CVD duration. In addition, the robustness of these findings was confirmed by comprehensive sensitivity analyses. Several potential limitations need to be considered. First, while our research cannot make a causal inference because of the observational nature, it provides evidence to support the importance of conducting future studies using a more rigorous study design, including RCTs. Second, although we considered a wide range of covariates and performed adequate sensitivity analyses, residual confounders and potential bias may still be present. Third, since the UK Biobank cohort included the largely white population, the results from this study may not be directly generalizable to other populations. Fourth, we were unable to access the relationship between dynamic changes in 25(OH)D concentrations and mortality because the 25(OH)D measurements were not repeated.

Conclusion

In the UK Biobank, nearly 60% of CVD patients were subjected to vitamin D deficiency. Among patients with existing CVD, increasing levels in serum 25(OH)D were independently associated with a decreased risk of all-cause and cause-specific mortality in a non-linear and dose-response manner. Compared with CVD patients with serum 25(OH)D ≥ 50 nmol/L, those subjected to vitamin D deficiency benefited more from an increment in serum 25(OH)D. Our findings provided novel clues awaiting further validation in clinical trials.

Abbreviations

CVD: Cardiovascular disease

25[OH]D: 25‑hydroxyvitamin D

CAD: coronary artery disease

CHF: chronic heart failure

RCTs: randomized controlled trials

CHD: coronary heart disease

AF: atrial fibrillation

HbA1c: glycated hemoglobinA1c

CRP: C-reactive protein

HRs: hazard ratios

CIs: confidence intervals

Declarations

Ethics approval and consent to participate: Ethical approval was obtained by the National Health Service National Research Ethics Service (11/NW/0382) and renewed by the North West–Haydock Research Ethics Committee (16/NW/0274). All participants provided informed written consent.

Consent for publication/Data sharing: Data from UK Biobank are available to all researchers upon making an application. This research has been conducted using the UK Biobank Resource under Application 63454.

Availability of data and materials: All data generated or analysed during this study are included in this published article [and its supplementary information files].

Competing interests: The authors declare that they have no competing interests.

Founding: L.C. is supported by the National Key Research and Development Program of China (2020YFC2006300) and the National Natural Science Foundation of China (82003461). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Authors' contributions: L.C. conceived and designed the study. L.D. and M.L. analyzed and interpreted data. L.D. drafted the manuscript, and L.C. revised it critically for important intellectual content. All authors provided final approval of the version to be published. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted.

Acknowledgements: We thank all UK Biobank participants and staff.

References

1. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet (London, England). 2018;392:1736-1788.
2. Pilz S, Verheyen N, Grübler MR, Tomaschitz A and März W. Vitamin D and cardiovascular disease prevention. Nat Rev Cardiol. 2016;13:404-417.
3. Michos ED, Cainzos-Achirica M, Heravi AS and Appel LJ. Vitamin D, Calcium Supplements, and Implications for Cardiovascular Health: JACC Focus Seminar. J Am Coll Cardiol. 2021;77:437-449.
4. Chin K, Zhao D, Tibuakuu M, Martin SS, Ndumele CE, Florido R, Windham BG, Guallar E, Lutsey PL and Michos ED. Physical Activity, Vitamin D, and Incident Atherosclerotic Cardiovascular Disease in Whites and Blacks: The ARIC Study. J Clin Endocrinol Metab. 2017;102:1227-1236.
5. Lutsey PL, Michos ED, Misialek JR, Pankow JS, Loehr L, Selvin E, Reis JP, Gross M, Eckfeldt JH and Folsom AR. Race and Vitamin D Binding Protein Gene Polymorphisms Modify the Association of 25-Hydroxyvitamin D and Incident Heart Failure: The ARIC (Atherosclerosis Risk in Communities) Study. JACC Heart Fail. 2015;3:347-356.
6. Robinson-Cohen C, Hoofnagle AN, Ix JH, Sachs MC, Tracy RP, Siscovick DS, Kestenbaum BR and de Boer IH. Racial differences in the association of serum 25-hydroxyvitamin D concentration with coronary heart disease events. JAMA. 2013;310:179-188.
7. Schneider ALC, Lutsey PL, Selvin E, Mosley TH, Sharrett AR, Carson KA, Post WS, Pankow JS, Folsom AR, Gottesman RF and Michos ED. Vitamin D, vitamin D binding protein gene polymorphisms, race and risk of incident stroke: the Atherosclerosis Risk in Communities (ARIC) study. Eur J Neurol. 2015;22:1220-1227.
8. Manson JE, Cook NR, Lee IM, Christen W, Bassuk SS, Mora S, Gibson H, Gordon D, Copeland T, D'Agostino D, Friedenberg G, Ridge C, Bubes V, Giovannucci EL, Willett WC and Buring JE. Vitamin D Supplements and Prevention of Cancer and Cardiovascular Disease. The New England journal of medicine. 2019;380:33-44.
9. Barbarawi M, Kheiri B, Zayed Y, Barbarawi O, Dhillon H, Swaid B, Yelangi A, Sundus S, Bachuwa G, Alkotob ML and Manson JE. Vitamin D Supplementation and Cardiovascular Disease Risks in More Than 83 000 Individuals in 21 Randomized Clinical Trials: A Meta-analysis. JAMA Cardiol. 2019;4:765-776.
10. Autier P, Mullie P, Macacu A, Dragomir M, Boniol M, Coppens K, Pizot C and Boniol M. Effect of vitamin D supplementation on non-skeletal disorders: a systematic review of meta-analyses and randomised trials. Lancet Diabetes Endocrinol. 2017;5.
11. Holick MF and Chen TC. Vitamin D deficiency: a worldwide problem with health consequences. Am J Clin Nutr. 2008;87:1080S-1086S.
12. Kim DH, Sabour S, Sagar UN, Adams S and Whellan DJ. Prevalence of hypovitaminosis D in cardiovascular diseases (from the National Health and Nutrition Examination Survey 2001 to 2004). Am J Cardiol. 2008;102:1540-1544.
13. Lee JH, Gadi R, Spertus JA, Tang F and O'Keefe JH. Prevalence of vitamin D deficiency in patients with acute myocardial infarction. Am J Cardiol. 2011;107:1636-1638.
14. Wu Z, Wang T, Zhu S and Li L. Effects of vitamin D supplementation as an adjuvant therapy in coronary artery disease patients. Scand Cardiovasc J. 2016;50.
15. Witte KK, Byrom R, Gierula J, Paton MF, Jamil HA, Lowry JE, Gillott RG, Barnes SA, Chumun H, Kearney LC, Greenwood JP, Plein S, Law GR, Pavitt S, Barth JH, Cubbon RM and Kearney MT. Effects of Vitamin D on Cardiac Function in Patients With Chronic HF: The VINDICATE Study. J Am Coll Cardiol. 2016;67:2593-2603.
16. Schleithoff SS, Zittermann A, Tenderich G, Berthold HK, Stehle P and Koerfer R. Vitamin D supplementation improves cytokine profiles in patients with congestive heart failure: a double-blind, randomized, placebo-controlled trial. Am J Clin Nutr. 2006;83:754-759.
17. Wasson LT, Shimbo D, Rubin MR, Shaffer JA, Schwartz JE and Davidson KW. Is vitamin D deficiency a risk factor for ischemic heart disease in patients with established cardiovascular disease? 10-year follow-up of the Nova Scotia Health Survey. Int J Cardiol. 2011;148:387-389.
18. Gotsman I, Shauer A, Zwas DR, Hellman Y, Keren A, Lotan C and Admon D. Vitamin D deficiency is a predictor of reduced survival in patients with heart failure; vitamin D supplementation improves outcome. Eur J Heart Fail. 2012;14:357-366.
19. Liu LCY, Voors AA, van Veldhuisen DJ, van der Veer E, Belonje AM, Szymanski MK, Silljé HHW, van Gilst WH, Jaarsma T and de Boer RA. Vitamin D status and outcomes in heart failure patients. Eur J Heart Fail. 2011;13:619-625.
20. Schierbeck LL, Jensen TS, Bang U, Jensen G, Køber L and Jensen J-EB. Parathyroid hormone and vitamin D--markers for cardiovascular and all cause mortality in heart failure. Eur J Heart Fail. 2011;13:626-632.
21. Dobnig H, Pilz S, Scharnagl H, Renner W, Seelhorst U, Wellnitz B, Kinkeldei J, Boehm BO, Weihrauch G and Maerz W. Independent association of low serum 25-hydroxyvitamin d and 1,25-dihydroxyvitamin d levels with all-cause and cardiovascular mortality. Arch Intern Med. 2008;168:1340-1349.
22. Degerud E, Nygård O, de Vogel S, Hoff R, Svingen GFT, Pedersen ER, Nilsen DWT, Nordrehaug JE, Midttun Ø, Ueland PM and Dierkes J. Plasma 25-Hydroxyvitamin D and Mortality in Patients With Suspected Stable Angina Pectoris. J Clin Endocrinol Metab. 2018;103:1161-1170.
23. Welles CC, Whooley MA, Karumanchi SA, Hod T, Thadhani R, Berg AH, Ix JH and Mukamal KJ. Vitamin D deficiency and cardiovascular events in patients with coronary heart disease: data from the Heart and Soul Study. Am J Epidemiol. 2014;179:1279-1287.
24. Grandi NC, Breitling LP, Vossen CY, Hahmann H, Wüsten B, März W, Rothenbacher D and Brenner H. Serum vitamin D and risk of secondary cardiovascular disease events in patients with stable coronary heart disease. Am Heart J. 2010;159:1044-1051.
25. Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, Downey P, Elliott P, Green J, Landray M, Liu B, Matthews P, Ong G, Pell J, Silman A, Young A, Sprosen T, Peakman T and Collins R. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015;12:e1001779.
26. Holick MF, Binkley NC, Bischoff-Ferrari HA, Gordon CM, Hanley DA, Heaney RP, Murad MH and Weaver CM. Evaluation, treatment, and prevention of vitamin D deficiency: an Endocrine Society clinical practice guideline. J Clin Endocrinol Metab. 2011;96:1911-1930.
27. Piercy KL, Troiano RP, Ballard RM, Carlson SA, Fulton JE, Galuska DA, George SM and Olson RD. The Physical Activity Guidelines for Americans. JAMA. 2018;320:2020-2028.
28. Mozaffarian D. Dietary and Policy Priorities for Cardiovascular Disease, Diabetes, and Obesity: A Comprehensive Review. Circulation. 2016;133:187-225.
29. Foster HME, Celis-Morales CA, Nicholl BI, Petermann-Rocha F, Pell JP, Gill JMR, O'Donnell CA and Mair FS. The effect of socioeconomic deprivation on the association between an extended measurement of unhealthy lifestyle factors and health outcomes: a prospective analysis of the UK Biobank cohort. Lancet Public Health. 2018;3:e576-e585.
30. Inker LA, Schmid CH, Tighiouart H, Eckfeldt JH, Feldman HI, Greene T, Kusek JW, Manzi J, Van Lente F, Zhang YL, Coresh J and Levey AS. Estimating glomerular filtration rate from serum creatinine and cystatin C. The New England journal of medicine. 2012;367:20-29.
31. Austin PC, Lee DS and Fine JP. Introduction to the Analysis of Survival Data in the Presence of Competing Risks. Circulation. 2016;133:601-609.
32. Lau B, Cole SR and Gange SJ. Competing risk regression models for epidemiologic data. Am J Epidemiol. 2009;170:244-256.
33. Khundmiri SJ, Murray RD and Lederer E. PTH and Vitamin D. Compr Physiol. 2016;6:561-601.
34. Chowdhury R, Kunutsor S, Vitezova A, Oliver-Williams C, Chowdhury S, Kiefte-de-Jong JC, Khan H, Baena CP, Prabhakaran D, Hoshen MB, Feldman BS, Pan A, Johnson L, Crowe F, Hu FB and Franco OH. Vitamin D and risk of cause specific death: systematic review and meta-analysis of observational cohort and randomised intervention studies. BMJ. 2014;348:g1903.

Tables

Table 1. Baseline characteristics of 37080 patients with cardiovascular disease in UK Biobank

Data were expressed as the mean (SD), median (IQR) or n (%).

*Adherance to a healthy diet was defined as participants who had at least 5 score of healthy diet components.

Abbreviations: BMI, body mass index; eGFRcr-cys, estimated glomerular filtration rate (creatinine–cystatin C equation); CVD, cardiovascular disease; IQR, interquartile range; SD, standard deviation.

Table 2. Adjusted means of cardiometabolic markers according to serum 25(OH)D concentrations

Data were expressed as the adjusted mean (95% CI) and adjusted for age (continuous), sex (male, female), ethnicity (White, mixed, Asian, Black, Chinese, others, or unknown), education (college or university, vocational qualification, upper secondary, lower secondary, others, or unknown), Townsend deprivation index (in quintiles), household income (<18000, 18000–30999, 31000–51999, 52000–100000, >100000 £, or unknown), smoking status (never smoker, former smoker, current smoker, or unknown), alcohol consumption (0, 0.1–4.9, 5.0–14.9, 15.0–19.9, 20.0–29.9, ≥30.0 g/day, or unknown), physical activity (inactive group, insufficient group, active group, or unknown), healthy diet score (in quintiles), BMI (<18.5, 18.5–22.9, 23.0–24.9, 25.0–29.9, 30.0–34.9, or ≥35.0 kg/m²), antihypertensive medication use (yes, no), cholesterol lowering medication use (yes, no), diabetes medication use (none, only oral medication, only insulin, or insulin and oral medication), history of cancer, diabetes, hypertension, and duration of CVD (<1.0, 1.0–4.9, 5.0–9.9, or ≥10.0 years).

Abbreviations: eGFRcr-cys, estimated glomerular filtration rate (creatinine–cystatin C equation); HbA1C, glycated hemoglobinA1c.

Table 3. Associations of serum 25(OH)D concentrations with all-cause and cause-specific mortality

Model 1: adjusted for age (continuous), sex (male, female), and ethnicity (White, mixed, Asian, Black, Chinese, others, or unknown)

Model 2: adjusted for model 1 plus education (college or university, vocational qualification, upper secondary, lower secondary, others, or unknown), Townsend deprivation index (in quintiles), household income (<18000, 18000–30999, 31000–51999, 52000–100000, >100000 £, or unknown), smoking status (never smoker, former smoker, current smoker, or unknown), alcohol consumption (0, 0.1–4.9, 5.0–14.9, 15.0–19.9, 20.0–29.9, ≥30.0 g/day, or unknown), physical activity (inactive group, insufficient group, active group, or unknown), healthy diet score (in quintiles), and BMI (<18.5, 18.5–22.9, 23.0–24.9, 25.0–29.9, 30.0–34.9, or ≥35.0 kg/m²).

Model 3: adjusted for model 2 plus eGFRcr-cys (<30.0, 30.0–60.0, 60.0–90.0, ≥90.0 ml min^-1 per 1.73 m²), C-reactive protein (in quintiles), antihypertensive medication use (yes, no), cholesterol lowering medication use (yes, no), diabetes medication use (none, only oral medication, only insulin, or insulin and oral medication), history of cancer, diabetes, hypertension, and duration of CVD (<1.0, 1.0–4.9, 5.0–9.9, or ≥10.0 years).

Abbreviations: BMI, body mass index; CVD, cardiovascular disease; eGFRcr-cys, estimated glomerular filtration rate (creatinine–cystatin C equation).

Table 4. Hazard ratios for all-cause and cause-specific mortality by serum 25(OH)D among patients with CVD subtypes

Data are presented as hazard ratio (95% confidence interval) and adjusted for  age (continuous), sex (male, female), ethnicity (White, mixed, Asian, Black, Chinese, others, or unknown), education (college or university, vocational qualification, upper secondary, lower secondary, others, or unknown), Townsend deprivation index (in quintiles), household income (<18000, 18000–30999, 31000–51999, 52000–100000, >100000 £, or unknown), smoking status (never smoker, former smoker, current smoker, or unknown), alcohol consumption (0, 0.1–4.9, 5.0–14.9, 15.0–19.9, 20.0–29.9, ≥30.0 g/day, or unknown), physical activity (inactive group, insufficient group, active group, or unknown), healthy diet score (in quintiles), BMI (<18.5, 18.5–22.9, 23.0–24.9, 25.0–29.9, 30.0–34.9, or ≥35.0 kg/m2), eGFRcr-cys (<30.0, 30.0–60.0, 60.0–90.0, ≥90.0 ml min^-1 per 1.73 m²), C-reactive protein (in quintiles), antihypertensive medication use (yes, no), cholesterol lowering medication use (yes, no), diabetes medication use (none, only oral medication, only insulin, or insulin and oral medication), history of cancer, diabetes, hypertension, and duration of CVD (<1.0, 1.0–4.9, 5.0–9.9, or ≥10.0 years).

Abbreviations: BMI, body mass index; CVD, cardiovascular disease; eGFRcr-cys, estimated glomerular filtration rate (creatinine–cystatin C equation).