Vitamin D and mortality: meta-analysis of individual participant data from a large consortium of cohort studies from Europe and the United States
BMJ 2014; 348 doi: https://doi.org/10.1136/bmj.g3656 (Published 17 June 2014) Cite this as: BMJ 2014;348:g3656- Ben Schöttker, postdoctoral scientist1,
- Rolf Jorde, professor23,
- Anne Peasey, postdoctoral scientist4,
- Barbara Thorand, senior researcher in epidemiology5,
- Eugène H J M Jansen, postdoctoral scientist6,
- Lisette de Groot, professor of nutrition and ageing7,
- Martinette Streppel, postdoctoral scientist7,
- Julian Gardiner, research associate4,
- José Manuèl Ordóñez-Mena, PhD candidate18,
- Laura Perna, postdoctoral scientist1,
- Tom Wilsgaard, professor9,
- Wolfgang Rathmann, senior researcher in epidemiology10,
- Edith Feskens, professor7,
- Ellen Kampman, professor7,
- Galatios Siganos, research associate9,
- Inger Njølstad, professor9,
- Ellisiv Bøgeberg Mathiesen, professor11,
- Růžena Kubínová, senior researcher12,
- Andrzej Pająk, professor13,
- Roman Topor-Madry, senior researcher13,
- Abdonas Tamosiunas, professor14,
- Maria Hughes, postdoctoral scientist15,
- Frank Kee, professor15,
- Martin Bobak, professor4,
- Antonia Trichopoulou, professor1617,
- Paolo Boffetta, professor1618,
- Hermann Brenner, professor1
- on behalf of the Consortium on Health and Ageing: Network of Cohorts in Europe and the United States (CHANCES)
- 1Division of Clinical Epidemiology and Aging Research, German Cancer Research Center, Im Neuenheimer Feld 581, 69120 Heidelberg, Germany
- 2Tromsø Endocrine Research Group, Department of Clinical Medicine, University of Tromsø (UiT) the Arctic University of Norway, 9037 Tromsø, Norway
- 3Division of Internal Medicine, University Hospital of North Norway, 9038 Tromsø, Norway
- 4Department Epidemiology and Public Health, University College London, London WC1E 6BT, UK
- 5Institute of Epidemiology II, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany
- 6Laboratory for Health Protection Research, National Institute for Public Health and the Environment, 3720 BA Bilthoven, Netherlands
- 7Division of Human Nutrition, Wageningen University, 6700 EV Wageningen, Netherlands
- 8Network Aging Research (NAR), University of Heidelberg, Heidelberg, Germany
- 9Epidemiology of Chronic Diseases Research Group, Department of community medicine, UiT the Arctic University of Norway, 9037 Tromsø, Norway
- 10German Diabetes Center, Institute for Biometrics and Epidemiology, Leibniz Center for Diabetes Research at Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
- 11Brain and Circulation Research Group, Department of clinical medicine, UiT the Arctic University of Norway, 9037 Tromsø, Norway
- 12National Institute of Public Health, Prague, Czech Republic
- 13Jagiellonian University Medical College, Faculty of Health Sciences, Krakow, Poland
- 14Institute of Cardiology of Lithuanian University of Health Sciences, Kaunas, Lithuania
- 15UKCRC Centre of Excellence for Public Health, Queens University Belfast, Northern Ireland, UK
- 16Hellenic Health Foundation, Kaisareias 13 & Alexandroupoleos, Athens 11527, Greece
- 17University of Athens, Medical School, Department of Hygiene, Epidemiology and Medical Statistics, Mikras Asias 75 st, Athens 11527, Greece
- 18Institute for Translational Epidemiology and Tisch Cancer Institute, Icahn School of Medicine at Mount Sinai, New York, NY, USA
- Correspondence to: B Schöttker b.schoettker{at}dkfz.de
- Accepted 19 May 2014
Abstract
Objective To investigate the association between serum 25-hydroxyvitamin D concentrations (25(OH)D) and mortality in a large consortium of cohort studies paying particular attention to potential age, sex, season, and country differences.
Design Meta-analysis of individual participant data of eight prospective cohort studies from Europe and the US.
Setting General population.
Participants 26 018 men and women aged 50-79 years
Main outcome measures All-cause, cardiovascular, and cancer mortality.
Results 25(OH)D concentrations varied strongly by season (higher in summer), country (higher in US and northern Europe) and sex (higher in men), but no consistent trend with age was observed. During follow-up, 6695 study participants died, among whom 2624 died of cardiovascular diseases and 2227 died of cancer. For each cohort and analysis, 25(OH)D quintiles were defined with cohort and subgroup specific cut-off values. Comparing bottom versus top quintiles resulted in a pooled risk ratio of 1.57 (95% CI 1.36 to 1.81) for all-cause mortality. Risk ratios for cardiovascular mortality were similar in magnitude to that for all-cause mortality in subjects both with and without a history of cardiovascular disease at baseline. With respect to cancer mortality, an association was only observed among subjects with a history of cancer (risk ratio, 1.70 (1.00 to 2.88)). Analyses using all quintiles suggest curvilinear, inverse, dose-response curves for the aforementioned relationships. No strong age, sex, season, or country specific differences were detected. Heterogeneity was low in most meta-analyses.
Conclusions Despite levels of 25(OH)D strongly varying with country, sex, and season, the association between 25(OH)D level and all-cause and cause-specific mortality was remarkably consistent. Results from a long term randomised controlled trial addressing longevity are being awaited before vitamin D supplementation can be recommended in most individuals with low 25(OH)D levels.
Introduction
Although vitamin D is obtained from diet and dietary supplements, the main source of vitamin D is its production in skin under the influence of solar ultraviolet B radiation. The most commonly measured vitamin D metabolite is serum 25-hydroxyvitamin D (25(OH)D) because of its greater half life (~3 weeks) and up to 1000-fold higher serum levels compared with the physiologically active metabolite 1,25-dihydroxyvitamin D (half life of a few hours).1 As ultraviolet B exposure varies with the latitude of regions and during the year, so mean 25(OH)D concentrations of populations also vary accordingly.2 3 Furthermore, women are generally more prone to low 25(OH)D concentrations than men (possibly because of a positive correlation of 25(OH)D with testosterone levels),4 and vitamin D deficiency is especially common among elderly people, who often have less sun exposure (because of reduced outdoor activity) and limited capacity of the skin to produce vitamin D metabolites.5 However, it is unclear if and to what extent these influences on production and maintenance of sufficient vitamin D metabolite concentrations affect the prognostic association of low 25(OH)D concentrations with mortality.
The objective of this meta-analysis of a consortium of eight prospective cohort studies from different European countries and the United States was to investigate the prognostic association of 25(OH)D measurements with all-cause, cardiovascular, and cancer mortality, paying particular attention to potential differences between countries, sexes, age groups, and seasons of blood draw.
Methods
This manuscript was written in adherence to the Strengthening The Reporting of Observational Studies in Epidemiology (STROBE) and the Meta-analysis Of Observational Studies in Epidemiology (MOOSE) statements.6 7
Study design and study population
This investigation is a meta-analysis of individual participant data of seven population based cohorts with measured 25(OH)D of the Consortium on Health and Ageing: Network of Cohorts in Europe and the United States (CHANCES, www.chancesfp7.eu) together with the third US National Health And Nutrition Examination Survey (NHANES III). The cohorts of the CHANCES consortium were chosen because all variables needed for this project were harmonised according to pre-defined rules. The harmonisation rules were discussed among the CHANCES partners until a consensus was reached. The NHANES survey was added to include a North American cohort and because it has an open access policy that enabled data access with subsequent conversion according to the CHANCES data harmonisation rules. The cohorts cover 16 European countries and the United States. The cohorts’ key characteristics are summarised in table 1⇓, and the study designs are described in the supplemental data and previous publications.2 8 9 10 11 12 All CHANCES cohorts and the NHANES are conducted in accordance with the declaration of Helsinki, and written informed consent was obtained from all study participants.
Inclusion and exclusion criteria
To make cohorts more comparable, analyses were restricted to study participants aged 50-79 years. Current smokers were excluded from the Tromsø study because the 25(OH)D assay used resulted in smokers having 15-20% higher 25(OH)D concentrations than non-smokers which was not reproducible with other assays.9 In addition, study participants with missing 25(OH)D measurement, missing values for covariates used in the main model, or participants lost to follow-up due to unknown reasons were excluded. The derivation of the final number of included study participants is shown for each cohort in table 1⇑.
Mortality ascertainment
All cohorts ascertained deaths by region, state, or country wide registries with a coverage of ≥98.5% and obtained information about the underlying cause of death for at least 97.6% of the deceased. The mean follow-up time ranged from 4.2 to 15.9 years (table 1⇑). All deaths coded with ICD-10 codes I00-I99 were considered cardiovascular deaths, and cancer deaths were defined by the ICD-10 codes C00-C97.
Measurement of 25(OH)D
Measurement of 25(OH)D was conducted in the full cohorts of ESTHER, Tromsø, and NHANES III, whereas in SENECA and MONICA/KORA random subsamples were measured, and in the three HAPIEE cohorts 25(OH)D was measured in a nested case-control design (including study participants who died from any cause as cases). Details of the sampling procedures, the design of the nested case-control study, and the manufacturers of the 25(OH)D assays are given in table 1⇑, and further details have been described previously.2 8 9 11 13 14
Covariate assessment
Sociodemographic, lifestyle, anthropometric, and disease history data were assessed by self administered questionnaires or in interviews. In addition to self reported information, some studies measured weight and height and validated the history of common chronic diseases by consulting medical records or registries (such as cancer registries) (supplementary table 1). If measured anthropometric data or validated diagnoses were available, these were used in the analysis, and self reported information was used only to fill missing information. Prevalent cardiovascular disease was defined as a history of myocardial infarction or stroke.
Statistical analyses
This meta-analysis of individual participant data followed a two step approach by analysing the cohorts individually first and conducting meta-analyses of the obtained effect estimates thereafter.
Differences in 25(OH)D concentrations according to baseline characteristics were assessed with the Wilcoxon rank-sum test. The longitudinal analyses with mortality outcomes were conducted using 25(OH)D quintiles because there is still no consensus about cut off values for 25(OH)D risk categories8 15 and the applied 25(OH)D assays in the cohorts were not commonly standardised. The cut off values for 25(OH)D quintiles were allowed to vary by cohort and by analysed subgroup to ensure that there was always 20% of the analysed study population in each quintile (see supplementary table 2 for applied cut off values).
Cox proportional hazards or logistic regression models (the latter for case-control studies nested in the HAPIEE cohorts) were used to estimate hazard ratios or odds ratios, respectively, with 95% confidence intervals for comparison of increasing 25(OH)D quintiles compared with the highest quintile. When the risk is low, logistic and Cox regression usually produce similar results and the odds ratio can be regarded as an approximation of the more accurate hazard ratio that also takes time to the event into account. We use the term “risk ratios” for both the effect estimates. Analyses for cardiovascular and cancer mortality were stratified by presence or absence of a history of the disease. For sensitivity analyses of all-cause mortality, data were restricted to subjects without a history of cardiovascular disease or cancer, respectively. For analyses on cause-specific mortality, competing risks analysis was employed to correct for the competing risk of death by other causes. Weights for subjects who died of competing causes of death were determined according to an extension of the Fine and Gray method.16
For all outcomes, three models were applied with adjustment for an increasing number of established risk factors for all-cause mortality and determinants of serum 25(OH)D concentrations. In model 1, adjustment was made for age, sex, and season of blood draw. Model 2 was the fully adjusted model with additional adjustment for education, body mass index, smoking, and physical activity. Model 3 additionally included diseases that could potentially be intermediates in the association of 25(OH)D and mortality (that is, a history of diabetes, hypertension, cardiovascular disease, and cancer). Model 2 results were considered the main results. Age was modelled linearly and the other variables categorically with dummy variables using the categories shown in table 2⇓. In addition, the following adjustments for cohort-specific characteristics were made in all three models: survey for MONICA/KORA, latitude of study centre for SENECA, and region (west, mid-west, north-east, south) as well as ethnicity (non-Hispanic white, non-Hispanic black, Mexican-American, other) for the NHANES III. Proportions of non-white study participants in the European studies were negligible.
For all meta-analyses, Mantel-Haenszel weighting and random effects models were used,17 taking the sample size of the cohorts and the possibility of statistical heterogeneity among the studies into account. The latter was examined with Cochrane’s Q test and the I2 statistic. In sensitivity analyses, the nested case-control studies and, separately, death in the first three years of follow up were excluded from the meta-analyses, but this did not change the overall findings and these results are therefore not shown. Dose-response relationships were assessed by estimation of a trend underlying the summary effect estimates of the 25(OH)D quintiles applying the best fitting equation of the form y=m x b (y, effect estimate; m, constant1, x=number of quintile (that is, 1, 2, 3, 4, or 5), b=constant2).
Subgroups to be analysed were chosen a priori with the aim to stratify for the important non-modifiable risk factors for vitamin D deficiency: sex, age, and season of blood draw.8 18 19 Cohort-specific subgroup analysis results were pooled in meta-analyses, and in a sensitivity analysis the pooling of results was restricted to those cohorts that could contribute data to all subgroup analyses (ESTHER, Tromsø and NHANES III). In a further sensitivity analysis, analyses were stratified by smoking status (current; yes or no) to examine whether excluding current smokers from the Tromsø study could have affected the results of the overall meta-analysis.
All statistical tests were two sided using an α level of 0.05. The meta-analyses were conducted with the statistical software Comprehensive Meta-Analysis 2.0 (Biostat, Englewood NJ, USA), the dose-response trend was assessed with Microsoft Excel 2010 (Redmond, Washington, USA), and all other analyses were conducted with SAS, version 9.2 (Cary, North Carolina, USA). The competing risks analysis was carried out with a SAS macro provided by G Heinze, Medical University of Vienna, Austria.20
Results
The baseline characteristics of the analytical samples from participating cohorts and the nested case-control studies are shown in supplementary tables 3 and 4, respectively. All cohorts included men and women and all age groups (50-79 years), with exception of the SENECA study, which only recruited in the age range 70-79 years. The median 25(OH)D concentrations by cohort varied between 24 and 62 nmol/L. Lower values in the Tromsø, MONICA/KORA, SENECA, and HAPIEE cohort studies can in part be explained by there being few or no serum samples being collected in summer or autumn. Interestingly, the medians of 25(OH)D concentrations in the nested case-control studies were consistently lower among cases than controls (supplementary table 4).
Table 2⇓ shows 25(OH)D concentrations by baseline characteristics of the study participants (HAPIEE cohorts are again shown separately for cases and controls). There was no clear trend in 25(OH)D by age, but the median 25(OH)D concentrations were consistently lower among women than men. As expected, there was a clear pattern of 25(OH)D concentrations measured in winter or spring being substantially lower than those measured in summer or autumn. In addition, the median 25(OH)D concentrations were, with few exceptions among the cases and controls of the HAPIEE cohorts, increasing with categories of education, lowest in obese individuals, lowest in current smokers, and higher among subjects who undertook vigorous physical activity. There were no systematic differences in 25(OH)D concentrations between subjects with and without history of chronic diseases (diabetes, hypertension, cardiovascular disease, or cancer).
The median follow-up time in individual cohorts varied between 4.2 and 15.8 years, during which 6695 study participants died, including 2624 deaths from cardiovascular diseases and 2227 deaths from cancer (table 1⇑). Compared with models adjusted for age, sex, and season (model 1), the risk ratios by 25(OH)D quintile (using the top quintile as reference) were attenuated for all outcomes by additional adjustment for conventional risk factors for premature mortality (model 2; supplementary tables 5-7). Further addition of common chronic diseases to the model (model 3) did not lead to further material attenuations of the risk ratios.
The forest plots of the meta-analyses of main results (model 2), comparing the bottom and the top 25(OH)D quintile are shown for all-cause mortality in fig 1⇓ and for cardiovascular and cancer mortality in fig 2⇓. A low degree of heterogeneity of effects in the individual cohort studies can be seen from the forest plots. The horizontal bars (the confidence intervals) overlap, and the boxes (the effect estimates) are mostly close together. The pooled effect estimate in the meta-analysis for all-cause mortality revealed a 1.6-fold higher mortality in the bottom quintile compared with the top quintile (risk ratio 1.57 (95% confidence interval 1.36 to 1.81), fig 1⇓). Furthermore, the association with cardiovascular mortality in subjects without a history of cardiovascular disease (risk ratio 1.41 (1.18 to 1.68)) and in subjects with a history of cardiovascular disease (risk ratio 1.65 (1.22 to 2.22)) and the association with cancer mortality in subjects with a history of cancer (1.70 (1.00 to 2.88)) were of similar strength and all statistically significant (fig 2⇓; detailed results shown in supplementary tables 6-9). The only exception was cancer mortality in subjects without a history of cancer, for whom there was no association (risk ratio 1.03 (0.89 to 1.20)).
The dose-response trends across all 25(OH)D quintiles show curvilinear inverse associations for all-cause mortality and cardiovascular mortality, and for cancer mortality in subjects with a history of cancer (fig 3⇓; detailed results shown in supplementary tables 5-9). For all-cause mortality, the risk ratio for bottom versus top 25(OH)D quintile was almost unchanged if subjects with a history of cardiovascular disease or of cancer were excluded (supplementary table 10-11).
In analyses stratified for sex, age, and season, the association of the bottom quintile of 25(OH)D concentration with the outcomes did not differ substantially in the subgroups (pooled results are shown in fig 4⇓, detailed results are shown in supplementary tables 12-20). Restricting the analyses to cohorts that can contribute to all subgroups produced similar results (not shown), as did stratification by smoking status or exclusion of events in the first three years of follow-up (not shown).
Heterogeneity statistics for the meta-analyses usually indicated a low, statistically non-significant level of heterogeneity. When statistically significant, heterogeneity could be attenuated to non-significant levels by removal of one outlier study (supplementary tables 5-20). In the numerous meta-analyses, six out of the eight studies were identified as the outlier between one and three times, which shows that being the outlier was rather randomly distributed among the studies and did not indicate incomparability of the study’s results.
Discussion
Principal findings
In this meta-analysis of eight prospective cohort studies from Europe and the United States, the lowest quintile of serum 25(OH)D concentration was associated with increased all-cause and cardiovascular mortality, with a curvilinear association between 25(OH)D concentration and these outcomes. An association with cancer mortality was only observed if subjects had a history of cancer. Results were consistent across study populations, sexes, age groups, and seasons of blood draw.
Comparison with other studies
This is the first meta-analysis of individual participant data on the relationship of 25(OH)D and mortality. Previous meta-analyses of published study results faced problems with hardly comparable vitamin D categories (including tertiles, quartiles, quintiles, or clinically defined categories). Nevertheless, the pooled effect estimates from this meta-analysis for all-cause mortality (risk ratio 1.57 (95% confidence interval 1.36 to 1.81)) and cardiovascular mortality (1.41 (1.18 to 1.68) in subjects without cardiovascular disease at baseline) comparing the lowest and highest quintile of 25(OH)D concentration agreed with the meta-analyses of Zittermann et al21 for all-cause mortality (risk ratio 1.41 (1.04 to 1.90)) and of Wang et al22 for cardiovascular mortality (1.42 (1.19 to 1.71)). The most recent meta-analysis, by Chowdhury et al,23 that recalculated published effect estimates to derive comparisons of the bottom versus the top tertile of vitamin D levels came to similar effect estimates for all-cause mortality (relative risk 1.35 (1.22 to 1.49)) and cardiovascular mortality (1.35 (1.13 to 1.61)). None of the studies included in the reviews was conducted in eastern European populations, which differ from western European populations with respect to lifestyle (such as alcohol consumption and nutrition) and socioeconomic characteristics.24 Our meta-analysis shows that the associations of low 25(OH)D concentrations with all-cause and cardiovascular mortality outcomes are as strong in eastern Europe as in western Europe, northern Europe, or the US.
With respect to cancer mortality, the recent meta-analyses of Chowdhury et al23 and Yin et al25 showed weak, albeit statistically significant, elevated pooled risk ratios (1.14 (1.01 to 1.29) and 1.20 (1.02 to 1.43), respectively, after change of reference category to highest 25(OH)D category in the meta-analysis of Yin et al). Our meta-analysis revealed a clear difference in the association of 25(OH)D and cancer mortality between subjects with a history of cancer (risk ratio 1.70 (1.00 to 2.88)) and those without (1.03 (0.89 to 1.20)). This difference might explain why the pooled effect estimates of the previous meta-analyses were weak and why a large degree of heterogeneity was seen among the included studies. Moreover, it speaks for an important role of vitamin D in cancer prognosis (as recently shown by others for colorectal and breast cancer prognosis26). However, we cannot exclude reverse causality, that is, that the cancer might have led to low 25(OH)D levels. Furthermore, our study with the endpoint cancer mortality cannot make assumptions about a potential role of vitamin D in early phases of the carcinogenic process. A recent systematic umbrella review of meta-analyses of observational studies assessing the association of 25(OH)D levels with site-specific cancer incidences concluded that there is evidence for an association of 25(OH)D levels with colorectal cancer whereas there is inconclusive or no evidence for an association of 25(OH)D with other cancer sites.27
The median 25(OH)D concentrations varied strongly between northern, central and eastern Europe and the US in the same seasons of blood collection. However, it is uncertain how much of these differences are related to different ultraviolet B exposure, vitamin D consumption in the diet, and vitamin D fortification habits in the countries or are just artefacts of the different 25(OH)D assays employed. The distributions of various covariates were similar across the studies in this meta-analysis, but differences in the season of sample collection could have led to some of the observed differences. Importantly, the 25(OH)D assays used in the eight cohort studies were not standardised to a common standard, and large differences in concentrations may simply be due to different assay standards and materials.28 29 However, a previous analysis of the SENECA study that used a common 25(OH)D assay for analysis in 12 European countries was in agreement with our finding that population means of 25(OH)D concentrations are higher in northern than in southern Europe (differences of ≥20 nmol/L were reported).2 The authors of the SENECA study mainly suspected country differences in vitamin D fortification of foods, such as margarine, as possible causes of the large regional differences.30 This could also explain the higher 25(OH)D levels in the US, where fortification of milk is the rule whereas it is the exception in Europe.31
In addition, this analysis in eight population based studies confirmed the observation of previous studies that 25(OH)D concentrations vary strongly by age, sex, season, education, obesity, physical activity, and smoking.8 18 19 Nevertheless, in our meta-analyses of mortality endpoints, the results were comparable despite cut-off values set at the quintiles varying by cohort (and therefore country), or by other factors such as age, sex, and season. Overall, this finding may support the view that low 25(OH)D concentrations might be a marker for a poor health status rather than a cause of premature mortality. The direction from a poor health state to vitamin D deficiency seems plausible because subjects with poor health typically spend less time outdoors because of reduced capacity for physical activity, limiting their production of provitamin D3 in the skin under solar ultraviolet B radiation.32
Implications for vitamin D supplementation recommendations
The question of causality could be best answered by randomised controlled trials. Two recent systematic reviews of such trials confirmed that vitamin D supply has an effect on mortality.23 33 However, the effects were weaker than suggested by observational studies, and vitamin D supplementation seemed to be effective only for the administration of vitamin D3 in subjects with low 25(OH)D levels at baseline.23 33 Mortality was only a secondary outcome in most of the trials, initially designed to study the effects of vitamin D and calcium supplementation on osteoporotic outcomes. New trials with a focus on non-skeletal outcomes are needed. Four large trials have started in 2012 and are described elsewhere.34 Among these, the VIDAL study, which is currently undertaking a feasibility study, will be particularly informative with respect to mortality outcomes because it will be specifically designed to assess these endpoints in older adults (65-84 year old citizens of Great Britain). First results of this trial are expected between 2017 and 2020. Until then, it seems reasonable to follow the US Institute of Medicine recommendations that evidence is not yet sufficient to recommend vitamin D3 supplementation for subjects with vitamin D insufficiency (defined as 30-50 nmol/L 25(OH)D) but that subjects with vitamin D deficiency (<30 nmol/L 25(OH)D) could profit from vitamin D3 supplementation by maintaining bone health and reduction of fracture risk.35
However, there is no consensus on the ideal cut-off values for vitamin D deficiency, because they are currently based on the prevention of osteoporotic outcomes, and other outcomes might also be relevant.8 15 In our view, all-cause mortality is an outcome with particularly high public health relevance and should therefore be given high priority. A recent analysis of the ESTHER cohort showed that the Institute of Medicine cut-offs for vitamin D deficiency and insufficiency are not only appropriate for osteoporotic outcomes but also for mortality prediction.8 The former ESTHER analysis had been adjusted for the main determinants of 25(OH)D levels and therefore corrected for variations in 25(OH)D levels according to season and sex. In this analysis, widely varying cut-off points specific for country, sex, and season produced similar results, which raises the question whether country, sex, and season might be considered for defining vitamin D deficiency. Future population based cohort studies with adequately standardised 25(OH)D measurements are essential to obtain sex- and season-specific cut-points for vitamin D deficiency. Caution should be used in generalising existing cut-offs to other regions with different latitudes.
Limitations and strengths
This meta-analysis of individual participant data has several strengths: the variety of cohorts from all over Europe and the US, the overall large size enabling subgroup analyses, almost complete registry-based follow-ups, and the common statistical analysis strategy. We did not use repeated 25(OH)D measurements during follow-up because this would not have had a strong influence on the results. Previous analyses of the Tromsø and ESTHER study showed that 25(OH)D concentrations are fairly stable even up to 14 years.8 36 The main limitation of this meta-analysis of prospective cohort studies is its observational nature. Despite adjustment for known potential confounders, we cannot rule out the possibility that the observed associations are confounded by other unmeasured factors, such as impairments of the immune system.32 In addition, residual confounding by variables that could have been measured with a higher degree of detail (such as physical activity) is possible but not a likely explanation for the observed strong associations between 25(OH)D concentrations and mortality outcomes.
Conclusion
In this meta-analysis of eight prospective cohort studies, the lowest 25(OH)D quintile was associated with increased all-cause mortality, cardiovascular mortality, and cancer mortality (in subjects with a history of cancer). The inverse association across quintiles was curvilinear, and the effects were remarkably consistent across countries, sexes, seasons of blood draw, and age groups despite 25(OH)D cut-off values varying according to these characteristics. To identify cut-off values for vitamin D deficiency, 25(OH)D variation by geographic region, sex, and season might need to be taken into account.
What is already known on this subject
Mean serum 25-hydroxyvitamin D (25(OH)D) concentrations vary by country, sex, age, and season of blood draw
Vitamin D deficiency has been linked to all-cause and cardiovascular mortality, whereas results on cancer mortality have been inconsistent
What this study adds
In this large consortium of eight cohort studies from Europe and the United States, the bottom 25(OH)D quintile was associated with increased all-cause and cardiovascular mortality and with cancer mortality in subjects with a history of cancer but not in subjects without a history of cancer. These relationships were compellingly consistent across countries, sexes, age groups, and seasons of blood draw, although cut-off values for the 25(OH)D quintiles varied among the analyses
In clinical practice, cut-off values for vitamin D deficiency might need to be made region-, sex-, and season-specific to identify those in the population with the relatively lowest 25(OH)D concentrations
Vitamin D may play an important role in cancer prognosis
Notes
Cite this as: BMJ 2014;348:g3656
Footnotes
Collaborators of the CHANCES group: Bernd Holleczek (Saarland Cancer Registry, Germany); Christa Meisinger (Institute of Epidemiology II, German Research Center for Environmental Health, Germany); Wolfgang Koenig (Department of Internal Medicine II-Cardiology, University of Ulm Medical Center, Germany); Nicole Jankovic (Division of Human Nutrition, Wageningen University, Netherlands); Kari Kuulasmaa (National Institute for Health and Welfare, Helsinki, Finland).
Contributors: All authors fulfil the criteria of authorship. Contributors participated in collection of data, analysis of data, coordination of the consortium, or obtaining of funding. Ben Schöttker is the guarantor of this work.
Funding: This analysis was part of the CHANCES project funded in the FP7 framework programme of DG-RESEARCH in the European Commission (Grant No 242244). The CHANCES project is coordinated by the Hellenic Health Foundation, Greece. Further funding sources of participating cohorts: The ESTHER study was funded by the Baden-Württemberg state Ministry of Science, Research and Arts (Stuttgart, Germany), the Federal Ministry of Education and Research (Berlin, Germany) and the Federal Ministry of Family Affairs, Senior Citizens, Women and Youth (Berlin, Germany). Measurements of 25(OH)D in men were conducted in the context of the German Cancer Aid project number 108250 and 108426. The Tromsø Study was performed by the University of Tromsø in cooperation with the National Health Screening Service. The study was funded by the Norwegian Research Council. The MONICA/KORA Augsburg cohort study was financed by the Helmholtz Zentrum München, German Research Center for Environmental Health which is funded by the German Federal Ministry of Education and Research and by the State of Bavaria and supported by grants from the German Federal Ministry of Education and Research. Additional funds for the present study were provided by the German Research Foundation (TH-784/2-1 and TH-784/2-2), the Federal Ministry of Health, the University of Ulm, the German Diabetes Center and the Ministry of Innovation, Science, Research, and Technology of the state North Rhine Westphalia. For this analysis, data of the MONICA/KORA Augsburg study were provided by the MORGAM Data Centre, at THL in Helsinki, Finland. CHANCES funding was used for data coding and transfer. The SENECA study was funded as a Concerted Action within the EURONUT programme of the EU. The HAPIEE study was funded by the Welcome Trust (064947 and 081081), the US National Institute on Ageing (R01 AG23522-01), and a grant from Mac Arthur Foundation.
The study’s funders played no role in study design; in the collection, analysis, and interpretation of data; in the writing of the report; and in the decision to submit the article for publication. The researchers were independent from the funders.
Competing interests: All authors have completed the ICMJE uniform disclosure form at www.icmje.org/coi_disclosure.pdf (available on request from the corresponding author) and declare: no support from any organisation for the submitted work; no financial relationships with any organisations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work.
Ethical approval: The included studies have been approved by local ethics committees: ESTHER: Medical Faculty of the University of Heidelberg and the Medical Association of Saarland; HAPIEE: University College London (Great Britain), National Institute of Public Health (Prague, Czech Republic), Jagiellonian University (Krakow, Poland) and Lithuanian University of Health Sciences (Kaunas, Lithuania); Tromsø Study: The Regional Committee for Medical and Health Research Ethics and the Data Inspectorate of Norway; MONICA/KORA: Local authorities and the Bavarian State Chamber of Physicians. SENECA: Local ethics approval by the SENECA participating centres (e.g. Wageningen University ethics committee).
Informed consent has been obtained from all participants included in the analysed studies and the studies are being conducted in accordance with the declaration of Helsinki.
Transparency: The lead author (the manuscript’s guarantor) affirms that the manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned have been explained.
Data sharing: No additional data available
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 3.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/3.0/.