Population based mean Vitamin D levels in 19 European Countries & COVID-19 Mortality
========================================================================================

* Amar Ahmad
* Christian Heumann
* Raghib Ali
* Tim Oliver

## ABSTRACT

**Objectives** Reports early in the epidemic linking low mean national Vitamin D level with increased COVID-19 death, and until recently little research on the impact of Vitamin D deficiency on severity of COVID-19, led to this update of the initial report studying mortality up to the end of January 2021.

**Design and Setting** Coronavirus pandemic data for 19 European countries were downloaded from Our World in Data, which was last updated on January 24, 2021. Data from March 21, 2020 to January 22, 2021 were included in the statistical analysis. Vitamin-D (25)-HD mean data were collected by literature review. Poisson mixed-effect model was used to model the data.

**Results** European countries with Vitamin-D (25)-HD mean less than or equal to 50 have higher COVID-19 death rates as compared with European countries with Vitamin-D (25)-HD mean greater than 50, relative risk of 2.155 (95% CI: 1.068 – 4.347, p-value = 0.032). A statistically significant negative moderate Spearman rank correlation was observed between Vitamin-D (25)-HD mean and the number of COVID-19 deaths for each 14-day period during the COVID-19 pandemic time period.

**Conclusions** The observation of the significantly lower COVID-19 mortality rates in countries with lowest annual sun exposure but highest mean Vitamin-D (25)-HD levels provides support for the use of food fortification. The need to consider re-configuring vaccine strategy due to emergence of large number of COVID-19 variants provides an opportunity to undertake such therapeutic randomized control trials.

Key words
*   Population
*   Vitamin D
*   Europe
*   COVID-19
*   Mortality
*   Poisson mixed effect model

## INTRODUCTION

The COVID-19 pandemic has created a major global health problem because of its high number of patients requiring intensive care and death rate particularly in patients over 65-70 yrs of age. (1-3) The impact has been unequal globally, primarily due to the younger populations in, for example, Africa where the majority of the population are under 65 (4) not getting ill, except for South Africa (5). Non-pharmaceutical interventions used to control transmission, particularly universal lockdowns have been the primary means of reducing case incidence and this has caused serious economic harm in many countries (6-8).

Vitamin-D has a potential role in infections, which is implied from its effect on the innate and adaptive immune responses (9) and the fact that respiratory infections tend to disappear during summer. Furthermore, there is existing evidence for an association between vitamin-D deficiency and a risk of influenza infection and when the individuals have Vitamin D deficiency, they have reduced infection when entered into supplementation randomized trials (10, 11), though dose and scheduling is far from clear.

In the 12 months since the onset of the COVID-19 pandemic, no major trials of Vitamin D have been completed. There has however been increasing interest with 5 reports demonstrating a strong impact of Vitamin D deficiency on the behavior of COVID-19 illness (12-16) and limited data from three small randomized trial reporting improved clinical outcome in patients receiving replacement therapy (17) that has been followed by posting on-line of a 930 cluster randomized trial with Calciferol reporting a 64% reducation in deaths.

Population-based mean levels of vitamin-D were shown in an analysis duting the first 2 months of the pandemic to correlate inversely with COVID-19 mortality (18). The aim of this paper is to update these analyses after 12 months of the pandemic through statistical analysis focusing on the association between vitamin D (25)-HD average and COVID-19 mortality rates based on public record of the number of COVID-19 deaths in 2-week periods from 21/02/2020 – 22/01/2021 in 19 European countries.

## STUDY DESIGN

Coronavirus pandemic data for 19 European countries were downloaded from “Our World in Data” (see [https://ourworldindata.org/coronavirus](https://ourworldindata.org/coronavirus)), which was last updated on January 24, 2021 (19). Data from March 21, 2020 to January 22, 2021 were included in the statistical analysis. Data on vitamin D deficiency were collected via a literature review (18, 20, 21).

## STATISTICAL ANALYSIS

Last observation carried forward was used to impute missing values, whenever data were not available (22). The vitamin D mean value (18) was missing for Greece which was (single) imputed by the median vitamin D value of all countries.

The crude mortality rate (CMR) was computed using data from January 22, 2021 by dividing the total number of COVID-19 deaths by the corresponding population per 100,000. Spearman rank correlation was estimated between vitamin D averages variable and the number of COVID-19 total deaths and the crude mortality rate respectively. A jackknife empirical 95% confidence interval for Spearman’s correlation was computed (23).

Then the vitamin D values were binarised into 0 (vitamin D ≤ 50) versus 1 (vitamin D > 50). Wilcoxon rank-sum test was used to compare the total number of COVID-19 deaths and CMR between countries who have vitamin D averages of less than or equals to 50 versus countries that have vitamin D averages greater than 50 respectively.

Univariate and bivariate generalised linear regression models with Quasi-Poisson distribution were performed using the total number of COVID-19 deaths by January 22nd 2021 as an outcome. The predictors were the categorical vitamin D variable and the proportion of age 70+ for each country. The population variable was added as an offset in the model. Relative risks (RRs) and 95% confidence intervals (95% CI) were computed with corresponding z- and p-values.

A time-point variable of 22 time-periods of 14 days was created between 01/03/2020 and 22/01/2021 at each time-point. A Spearman rank correlation was used to evaluate the association between the number of COVID-19 deaths at each time-point and the average of vitamin D of the 19 European countries. To compute an approximate 95% confidence interval (using the normality approximation) for the estimated Spearman rank correlation, a leave-one-time-period-out cross-validation was performed at each time-point.

A Poisson mixed effects regression model (24) was fitted with the total number of new COVID-19 deaths as an outcome. The percentage of age 70+ and binarized vitamin D were used as fixed effect predictors. The time-periods and the 19 European countries variable were used as random effect variables with (04/04/2020 to 17/04/2020) and the UK as reference groups respectively. The populations variable was added as an offset in the model. Relative risks (RRs) and 95% confidence intervals (95% CI) were computed with corresponding z- and p-values.

All applied statistical tests were two-sided, p-value < 0.05 were considered statistically significant. Statistical analyses were performed in R version 4.0.2 (25)(25)(25).

## RESULTS

The data used in this statistical analysis are shown in Table 1 in descending order of the countries level of Vitamin D level. It can be seen visually that the crude mortality rate (CMR) increases with decreasing mean Vitamin D level. A statistically significant moderate negative Spearman’s ρ correlation was observed between the total number of COVID-19 deaths and the average of vitamin D, ρ = -0.516 (95% CI: -0.860 – -0.168) as well as between COVID-19 crude mortality rates and the average of vitamin D, ρ = -0.430 (95% CI: -0.805 – -0.081) (Figure 1). A similar result was seen in the numbers of Table 2. This shows the Spearman’s rho correlation analysis between the number of new COVID-19 deaths at each of the 22 time periods and the same reported mean vitamin D used on each occasion. All the observed negative Spearman’s rho correlation values were statistically significant.

View this table:
[Table 1:](http://medrxiv.org/content/early/2021/10/25/2021.03.11.21253361/T1)

Table 1: 
Data used in the statistical analysis. Data sorted by vitamin-D values

View this table:
[Table 2:](http://medrxiv.org/content/early/2021/10/25/2021.03.11.21253361/T2)

Table 2: 
Spearman correlation (95% CI) between vitamin-D and the new number of COVID-19 deaths at different time-points

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/10/25/2021.03.11.21253361/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2021/10/25/2021.03.11.21253361/F1)

Figure 1: Distribution of the total number of COVID-19 deaths with the vitamin-D group

Figure 1 shows the distribution of the total number of COVID-19 deaths and CMR within the vitamin D group (≤50 versus >50) respectively. A statistically significant difference was observed between the vitamin D groups (≤50 versus >50), Wilcoxon sum-rank test p-value of 0.036 and 0.012 respectively.

Table 3 shows the results of the univariate and bivariate generalised linear regression model with family quasi-Poisson using the total number of COVID-19 deaths from 22/03/2020 as an outcome. The countries with vitamin D average ≤ 50 have higher COVID-19 death rates as compared with countries with vitamin D average > 50, RR of 1.642 (95% CI 1.274 – 2.118, p-value = 0.02) univariately. This result holds after adjusting for the population age structure (ie percentage of the population with age 70+), RR of 1.663 (95% CI: 1.293 – 2.140, p-value < 0.0001). This result is similar to the result from the fitted Poisson mixed-effects models (Table 4), where countries with vitamin D average ≤ 50 have higher COVID-19 death rates as compared with countries with vitamin D average > 50 before and after adjusting for population age structure with RR of 2.197 (95% CI: 1.131 – 4.271, p-value = 0.02) and 2.155 (95% CI: 1.068 – 4.347, p-value = 0.032) respectively. The percentage of age 70+ was not statistically significant in either a univariate and bivariate regression models (Table 3 and 4).

View this table:
[Table 3:](http://medrxiv.org/content/early/2021/10/25/2021.03.11.21253361/T3)

Table 3: 
Univariate and bivariate quasi-Poisson regression models. The percentage of age 70+ and binarized vitamin D variables were added together

View this table:
[Table 4:](http://medrxiv.org/content/early/2021/10/25/2021.03.11.21253361/T4)

Table 4: 
Poisson mixed effect regression models. Univariate models, where fitted with single fixed effect variable. The percentage of age 70+ and binarized vitamin D variables were added together as fix-effect predictors.

## CONCLUSION

Despite regularly assuring the UK population that the pandemic would lessen during our summer months, there has been no funding for clinical trials of Vitamin D supplements and negative reviews mainly based on others data by PHE, SACN, NICE ([https://www.nice.org.uk/guidance/ng187](https://www.nice.org.uk/guidance/ng187)). This negative view focuses attention onto an on-going controversy dating back to the 1950s, i.e. whether food should be supplemented with Vitamin D. In the early years after the 2nd world war in the UK enthusiasm for supplementation was so great that children died from overdose hypercalcaemia. As a result, the UK has been more reluctant to sanction such supplementation while it has been used increasingly used in the Scandinavian countries (26). This could explain why most of these countries, despite having less sunshine than the UK, have lower RRs. The increasing data to support this interpretation started with Meltzer et al studying Vitamin D levels in 489 attendees in an urban medical centre in Chicago (13). They demonstrated a higher frequency of acquisition of COVID-19 cases in the Vitamin D deficient group (21.9%) than in the sufficient group (12.2%). This was subsequently confirmed by a much larger study (14) and two smaller studies (12, 15). The Boston group also reported data suggesting a link between Vitamin D deficiency and increased severity of disease (16).

In a recent systematic review and meta-analysis, Pereira et al (27) show, that vitamin D insufficiency increased COVID-19 mortality (OR = 1.82, 95% CI = 1.06–2.58).

Using UK Biobank data of 502,624 participants aged 37–73 years between 2006 and 2010, Hastie et al (28, 29) show that COVID-19 participants who had COVID-19 have lower vitamin-D level with median of 43.8 (IQR: 28.7-61.6) as compared with those who had no COVID-19 median of 47.2 (IQR: 32.7–62.7), Wilcoxon’s p-value < 0.01. However, this result was not statistically significant after adjustment for confounders, though this did not take allowance for the fact that some of the confounders were also associated with Vitamin D deficiency.

In a cohort of 185 patients at the Medical University Hospital Heidelberg-Germany, vitamin-D deficiency was associated with higher risk of invasive mechanical ventilation and deaths after adjusting for age, gender, and comorbidities, HR of 6.12 (95% CI: 2.79–13.42, p < 0.001) and 14.73 (95% CI 4.16–52.19, p < 0.001) (30, 31).

Karahan and Katkat (32) show that vitamin-D insufficiency was present in 93.1% of the patients with severe-critical COVID-19, and that vitamin-D 25(OH) mean was significantly lower in patients with severe-critical COVID-19 compared with moderate COVID-19.

Despite all of these strong indications that Vitamin D could play an important role in the control of the COVID-19 pandemic, formally published papers have less than 500 patients recruited to randomized trials of Vitamin D(27). There has been one prospective observational study on 410 patients from India which at first sight was considered as a failure of Vitamin D supplementation ([https://www.researchsquare.com/article/rs-129238/v1](https://www.researchsquare.com/article/rs-129238/v1), currently undergoing major revision). As it gave Vitamin D to two thirds of those with Vitamin D deficiency who were also younger than those who weren’t, the inverse of what is expected from other studies and overall there was a very low mortality rate, there was no difference in survival between those with and those without Vitamin D deficiency.

No data has been published from the UK about the effect of Vitamin D replacement. This is despite a continuous debate since April 2020 about Vitamin D deficiency contributing to excess deaths in the BAME community ([https://www.gov.uk/government/publications/covid-19-understanding-the-impact-on-bame-communities](https://www.gov.uk/government/publications/covid-19-understanding-the-impact-on-bame-communities)). Throughout the year, the BAME community are known to have 20% less circulating vitamin D (33) and even in summer more than 30% have severe deficiency compared to less than 6% in white European population (34) and probably have done so for many of the years since birth. In addition, there are two studies, one a randomized trial, that have clearly demonstrated that they need 10 times more Vitamin D than recommended in the SACN report (35, 36).

Clearly, accelerating the role out of Vaccines around the world should not diminish. However as our data could not adjust for well established COVID-19 confounders such as Diabetes, Obesity, social deprivation and poverty, the increasing availability of simple finger-prick techniques ([www.vitamindtest.org.uk](http://www.vitamindtest.org.uk)) that speed up measuring of Vitamin D should make better selection for such therapy trials easier to recruit. Given reports published on line but still under review (see [https://www.imperial.ac.uk/media/imperial-college/institute-of-global-health-innovation/REACT-2-round-5-preprint.pdf](https://www.imperial.ac.uk/media/imperial-college/institute-of-global-health-innovation/REACT-2-round-5-preprint.pdf) table 4) that the over 70’s population have a 62% failure rate to produce IgG antibodies >21 days after first vaccination compared to 24% in 50-69 year olds and 7% in <50 year olds and the older age groups have a well recognized risk of lower Vitamin D levels, there could be a strong case to evaluate Vitamin D supplements in such patients.

## Conclusion

The data from this statistical analysis shows a strong and statistically significant association between the Vitamin D deficiency and the total number of COVID-19 deaths in the 19 European countries included in this statistical analysis. The new vitamin D techniques for easier detection of deficiency using the finger prick technology should enable better selection of patients to benefit from treatment and prove it with .appropriately selected patients in treatment trials and reduce poor immune response in vaccine recipients

## Data Availability

Data available through Amar Ahmad email amar.ahmad@qmul.ac.uk

## Footnotes

*   The data from this statistical analysis shows a strong and statistically significant association between the Vitamin D deficiency and the total number of COVID19 deaths in the 19 European countries included in this statistical analysis. The new vitamin D techniques for easier detection of deficiency using the finger prick technology should enable better selection of patients to benefit from treatment and prove it with .appropriately selected patients in treatment trials and reduce poor immune response in vaccine recipients

*   Received March 11, 2021.
*   Revision received October 25, 2021.
*   Accepted October 25, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## REFERENCES

1.  1.Chen N, Zhou M, Dong X, Qu J, Gong F, Han Y, et al. Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study. Lancet. 2020;395(10223):507–13.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30211-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F25%2F2021.03.11.21253361.atom) 

2.  2.Sun K, Chen J, Viboud C. Early epidemiological analysis of the coronavirus disease 2019 outbreak based on crowdsourced data: a population-level observational study. Lancet Digit Health. 2020;2(4):e201–e8.
    
    

3.  3.Li Q, Guan X, Wu P, Wang X, Zhou L, Tong Y, et al. Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia. N Engl J Med. 2020;382(13):1199–207.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001316&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F25%2F2021.03.11.21253361.atom) 

4.  4.Massinga Loembe M, Tshangela A, Salyer SJ, Varma JK, Ouma AEO, Nkengasong JN. COVID-19 in Africa: the spread and response. Nat Med. 2020;26(7):999–1003.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0961-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F25%2F2021.03.11.21253361.atom) 

5.  5.Dyer O. Covid-19: Excess deaths point to hidden toll in South Africa as cases surge. BMJ. 2020;370:m3038.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzNzAvanVsMzBfOS9tMzAzOCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEwLzI1LzIwMjEuMDMuMTEuMjEyNTMzNjEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

6.  6.Strang P, Furst P, Schultz T. Excess deaths from COVID-19 correlate with age and socio-economic status. A database study in the Stockholm region. Ups J Med Sci. 2020;125(4):297–304.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F25%2F2021.03.11.21253361.atom) 

7.  7.Mele M, Magazzino C. Pollution, economic growth, and COVID-19 deaths in India: a machine learning evidence. Environ Sci Pollut Res Int. 2020.
    
    

8.  8.Adams-Prassl A, Cloyne J, Costa Dias M, Parey M, Ziliak JP. The COVID-19 Economic Crisis. Fisc Stud. 2020;41(3):489–92.
    
    

9.  9.Gruber-Bzura BM. Vitamin D and Influenza-Prevention or Therapy? Int J Mol Sci. 2018;19(8).
    
    

10. 10.Martineau AR, Jolliffe DA, Hooper RL, Greenberg L, Aloia JF, Bergman P, et al. Vitamin D supplementation to prevent acute respiratory tract infections: systematic review and meta-analysis of individual participant data. BMJ. 2017;356:i6583.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzNTYvZmViMTVfNC9pNjU4MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEwLzI1LzIwMjEuMDMuMTEuMjEyNTMzNjEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

11. 11.Jolliffe DA, Camargo CA, Sluyter JD, Aglipay M, Aloia JF, Ganmaa D, et al. Vitamin D supplementation to prevent acute respiratory infections: systematic review and meta-analysis of aggregate data from randomised controlled trials. medRxiv. 2020.
    
    

12. 12.Panagiotou G, Tee SA, Ihsan Y, Athar W, Marchitelli G, Kelly D, et al. Low serum 25-hydroxyvitamin D (25[OH]D) levels in patients hospitalized with COVID-19 are associated with greater disease severity. Clin Endocrinol (Oxf). 2020.
    
    

13. 13.Meltzer DO, Best TJ, Zhang H, Vokes T, Arora V, Solway J. Association of Vitamin D Status and Other Clinical Characteristics With COVID-19 Test Results. JAMA Netw Open. 2020;3(9):e2019722.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F25%2F2021.03.11.21253361.atom) 

14. 14.Kaufman HW, Niles JK, Kroll MH, Bi C, Holick MF. SARS-CoV-2 positivity rates associated with circulating 25-hydroxyvitamin D levels. PLoS One. 2020;15(9):e0239252.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0239252&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32941512&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F25%2F2021.03.11.21253361.atom) 

15. 15.Merzon E, Tworowski D, Gorohovski A, Vinker S, Golan Cohen A, Green I, et al. Low plasma 25(OH) vitamin D level is associated with increased risk of COVID-19 infection: an Israeli population-based study. FEBS J. 2020.
    
    

16. 16.Maghbooli Z, Sahraian MA, Ebrahimi M, Pazoki M, Kafan S, Tabriz HM, et al. Vitamin D sufficiency, a serum 25-hydroxyvitamin D at least 30 ng/mL reduced risk for adverse clinical outcomes in patients with COVID-19 infection. PLoS One. 2020;15(9):e0239799.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0239799&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F25%2F2021.03.11.21253361.atom) 

17. 17.Nikniaz L, Akbarzadeh M A.,, Hosseinifard H, Hosseini M. Vitamin D supplementation on mortality rate and clinical outcomes of COVID-19 patients: A systematic review and meta-analysis [www.medrxiv.org/content/10.1101/2021.01.04.21249219v1.full.pdf](http://www.medrxiv.org/content/10.1101/2021.01.04.21249219v1.full.pdf). medrxiv [Internet]. 2021.
    
    

18. 18.Ilie PC, Stefanescu S, Smith L. The role of vitamin D in the prevention of coronavirus disease 2019 infection and mortality. Aging Clin Exp Res. 2020;32(7):1195–8.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F25%2F2021.03.11.21253361.atom) 

19. 19.Hasell J, Mathieu E, Beltekian D, Macdonald B, Giattino C, Ortiz-Ospina E, et al. A cross-country database of COVID-19 testing. Sci Data. 2020;7(1):345.
    
    

20. 20.Singh S, Kaur R, Singh RK. Revisiting the role of vitamin D levels in the prevention of COVID-19 infection and mortality in European countries post infections peak. Aging Clin Exp Res. 2020;32(8):1609–12.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s40520-020-01619-8&link_type=DOI) 

21. 21.Lips P, Cashman KD, Lamberg-Allardt C, Bischoff-Ferrari HA, Obermayer-Pietsch B, Bianchi ML, et al. Current vitamin D status in European and Middle East countries and strategies to prevent vitamin D deficiency: a position statement of the European Calcified Tissue Society. Eur J Endocrinol. 2019;180(4):P23–P54.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1530/EJE-18-0736&link_type=DOI) 

22. 22.Vickers AJ, Altman DG. Statistics notes: missing outcomes in randomised trials. BMJ. 2013;346:f3438.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzNDYvanVuMDZfMi9mMzQzOCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzEwLzI1LzIwMjEuMDMuMTEuMjEyNTMzNjEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

23. 23.Wang R, Peng J. Jackknife Empirical likelihood intervals for Spearman’s rho. Peng, L.. North American Actuarial Journal. 2011;15(4):475–86.
    
    

24. 24.Bates D, Mächler M, Bolker BS W., “Fitting Linear Mixed-Effects Models Using lme4.” Journal of Statistical Software. 2015;67(1):1–48.
    
    

25. 25.RCT. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.2020.
    
    

26. 26.Spiro A, Buttriss JL. Vitamin D: An overview of vitamin D status and intake in Europe. Nutr Bull. 2014;39(4):322–50.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/nbu.12108&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25635171&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F25%2F2021.03.11.21253361.atom) 

27. 27.Pereira M, Dantas Damascena A, Galvao Azevedo LM, de Almeida Oliveira T, da Mota Santana J. Vitamin D deficiency aggravates COVID-19: systematic review and meta-analysis. Crit Rev Food Sci Nutr. 2020:1–9.
    
    

28. 28.Hastie CE, Mackay DF, Ho F, Celis-Morales CA, Katikireddi SV, Niedzwiedz CL, et al. Vitamin D concentrations and COVID-19 infection in UK Biobank. Diabetes Metab Syndr. 2020;14(4):561–5.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F25%2F2021.03.11.21253361.atom) 

29. 29.Hastie CE, Mackay DF, Ho F, Celis-Morales CA, Katikireddi SV, Niedzwiedz CL, et al. Corrigendum to “Vitamin D concentrations and COVID-19 infection in UK Biobank” [Diabetes Metabol Syndr: Clin Res Rev 2020 14 (4) 561-5]. Diabetes Metab Syndr. 2020;14(5):1315–6.
    
    

30. 30.Radujkovic A, Hippchen T, Tiwari-Heckler S, Dreher S, Boxberger M, Merle U. Vitamin D Deficiency and Outcome of COVID-19 Patients. Nutrients. 2020;12(9).
    
    

31. 31.Radujkovic A, Merle U. Reply to: “Vitamin D Insufficiency May Account for Almost Nine of Ten COVID-19 Deaths: Time to Act. Comment on: Vitamin D Deficiency and Outcome of COVID-19 Patients. Nutrients 2020, 12, 2757”. Nutrients. 2020;12(12).
    
    

32. 32.Karahan S, Katkat F. Impact of Serum 25(OH) Vitamin D Level on Mortality in Patients with COVID-19 in Turkey. J Nutr Health Aging. 2021;25(2):189–96.
    
    

33. 33.Raisi-Estabragh Z, McCracken C, Bethell MS, Cooper J, Cooper C, Caulfield MJ, et al. Greater risk of severe COVID-19 in Black, Asian and Minority Ethnic populations is not explained by cardiometabolic, socioeconomic or behavioural factors, or by 25(OH)-vitamin D status: study of 1326 cases from the UK Biobank. J Public Health (Oxf). 2020;42(3):451–60.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/pubmed/fdaa095&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F25%2F2021.03.11.21253361.atom) 

34. 34.Sutherland JP, Zhou A, Leach MJ, Hypponen E. Differences and determinants of vitamin D deficiency among UK biobank participants: A cross-ethnic and socioeconomic study. Clin Nutr. 2020.
    
    

35. 35.Hollis BW, Johnson D, Hulsey TC, Ebeling M, Wagner CL. Vitamin D supplementation during pregnancy: double-blind, randomized clinical trial of safety and effectiveness. J Bone Miner Res. 2011;26(10):2341–57.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jbmr.463&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21706518&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F10%2F25%2F2021.03.11.21253361.atom) 

36. 36.Kiely M, Cashman KD. Summary Outcomes of the ODIN Project on Food Fortification for Vitamin D Deficiency Prevention. Int J Environ Res Public Health. 2018;15(11).