Association between population vitamin D status and SARS-CoV-2 related serious-critical illness and deaths: An ecological integrative approach

Abstract

BACKGROUND

Vitamin D population status may have possible unappreciated consequences to the coronavirus disease 2019 (COVID-19) pandemic. Α significant association between vitamin D sufficiency and reduction in clinical severity and inpatient mortality from COVID-19 disease has recently been shown, while a recent study has claimed lower COVID-19 cases in European countries with a better vitamin D status. Low serum 25-hydroxyvitamin-D [25(OH)D] was identified as an independent risk factor for COVID-19 infection and hospitalization, and administration of 0.532 mg (21280 IU) of calcifediol or 25(OH)D, followed by 0.266 mg on days 3 and 7 and then weekly until discharge or intensive care unit admission significantly reduced the need for intensive care unit treatment.

AIM

To elucidate the role of vitamin D European population status in the COVID-19 pandemic, data from the Worldometer were analyzed.

METHODS

Linear regression explored the correlation between published representative-standardized population vitamin D concentrations and the number of total cases/million (M), recovered/M, deaths/M and serious-critically ill/M from COVID-19 for 26 European countries populated > 4 M (Worldometer). Life expectancy was analyzed with semi-parametric regression. Weighted analysis of variance/analysis of covariance evaluated serious-critical/M and deaths/M by the vitamin D population status: Deficient < 50, insufficient: 50-62.5, mildly insufficient > 62.5-75 and sufficient > 75 nmol/L, while controlling for life expectancy for deaths/M. Statistical analyses were performed in XLSTAT LIFE SCIENCE and R (SemiPar Library).

RESULTS

Linear regression found no correlation between population vitamin D concentrations and the total cases-recovered/M, but negative correlations predicting a reduction of 47%-64%-80% in serious-critical illnesses/M and of 61%-82%-102.4% in deaths/M further enhanced when adapting for life expectancy by 133-177-221% if 25(OH)D concentrations reach 100-125-150 nmol/L, sustained on August 15, 2020, indicating a truthful association. Weighted analysis of variance was performed to evaluate serious-critical/M (r² = 0.22) by the vitamin D population status and analysis of covariance the deaths/M (r² = 0.629) controlling for life expectancy (r² = 0.47). Serious-critical showed a decreasing trend (P < 0.001) from population status deficient (P < 0.001) to insufficient by 9.2% (P < 0.001), to mildly insufficient by 47.6% (P < 0.044) and to sufficient by 100% (reference, P < 0.001). For deaths/M the respective decreasing trend (P < 0.001) was 62.9% from deficient (P < 0.001) to insufficient (P < 0.001), 65.15% to mildly insufficient (P < 0.001) and 78.8% to sufficient (P = 0.041).

CONCLUSION

Achieving serum 25(OH)D 100-150 nmol/L (40-60 ng/mL) (upper tolerable daily doses followed by maintenance proposed doses not requiring medical supervision, Endocrine Society) may protect from serious-critical illness/death from COVID-19 disease.

Key Words: COVID-19; SARS-CoV-2; Vitamin D status; Vitamin D concentrations; 25-hydroxyvitamin-D; Immunity

Core Tip: To elucidate the role of vitamin D in the coronavirus disease 2019 (COVID-19) pandemic, we examined associations between published representative and standardized European population vitamin D data and the Worldometer COVID-19 data. Linear regression found no correlation between population vitamin D concentrations and the total cases-recovered/million (M), but negative correlations predicting a reduction of 47%-64%-80% in serious-critical illnesses/M and of 61%-82%-102.4% in deaths/M further enhanced when adapting for life expectancy by 133-177-221% if 25-hydroxyvitamin-D concentrations reach 100-125-150 nmol/L. Weighted analysis of variance/analysis of covariance showed a decreasing trend (P < 0.001) evaluating serious-critical/M (r² = 0.22) and the deaths/M (r² = 0.629) after controlling for life expectancy (r² = 0.47), by vitamin D population status, respectively.

INTRODUCTION

Vitamin D deficiency and insufficiency is a global health issue affecting probably many more than 1 billion children and adults worldwide, with institutionalized elderly being at higher risk of exhibiting lower 25-hydroxyvitamin-D [25(OH)D] blood concentrations. According to a systematic review of vitamin D status in populations worldwide, 37.3% of the studies reported 25(OH)D mean concentrations < 50 nmol/L in newborns and institutionalized elderly, who are at higher risk of exhibiting lower 25(OH)D concentrations[1]. Public health policy development is needed to reduce risk for potential health consequences of an inadequate vitamin D status[1], with consequences that should not be underestimated, especially now with this unprecedented pandemic of coronavirus disease 2019 (COVID-19)[2]. The initial universal lockdown for a period of 2-3 mo and the consequent repeated lockdowns along with the social distancing measures would further reduce the incidental solar vitamin D3 production, worsening the population’s vitamin D status[3]. Strong evidence supports the role of vitamin D particularly in preventing rickets and osteomalacia[4]. While circulating 25(OH)D concentrations below 30 nmol/L (12 ng/mL) are associated with an increased risk of rickets/osteomalacia, 25(OH)D concentrations between 50-125 nmol/L (20-50 ng/mL) appear to be safe and sufficient to promote skeletal health in the general population[5]. A serum 25(OH)D concentration of at least ≥ 50 nmol/L at the end of winter (10-20 nmol/L higher at the end of summer, to allow for seasonal decrease) is required for optimal musculoskeletal health[6]. Supplements of vitamin D in low doses together with calcium, alone or in combination with antiresorptive drugs may prevent hip or any type of fracture and have been evaluated in osteoporotic and osteopenic patients for primary as well as secondary prevention[7-9]. However, the role of vitamin D in innate and adaptive immunity remains rather underappreciated, with possible consequences and public health implications, leading to an increased risk for infectious diseases, autoimmune disorders and cancers[10]. Even if a recent randomized control trial (RCT) did not show lower incidence of invasive cancer in men ≥ 50 years or women ≥ 55 receiving 2000 IU of vitamin D3 daily up to 5 years[11], the study did report a statistically significant 25% reduced risk for cancer mortality. The study, however, had several limitations. Only 13% of the participants were vitamin D deficient [25(OH)D < 50 nmol/L], and 42%-45% of the participants were receiving a vitamin D supplement and multivitamins at inclusion. The participants, including the placebo group, were permitted to take up to 800 IU of vitamin D daily. This is the likely explanation why the mean baseline blood concentration of 25(OH)D was 74.5 nmol/L for the participants in this study[11]. The optimal 25(OH)D concentration is at least 75 nmol/L (30 ng/mL), which is what the mean baseline level was for the participants in the VITAL study. Secondary analyses from the VITAL study should be also considered as they indicate that the vitamin D dose was too low, since significant benefits were found for cancer incidence for those with body mass index (BMI) < 25 kg/m² and almost as significant for blacks. In fact, the authors speculated that the possible trial regimen–associated effects on cancer incidence among normal-weight participants and suggestive effects among black participants, which contrast with the null cardiovascular findings in these groups, may be explained by different vitamin D requirements for these outcomes. The Endocrine Society, which made its recommendations in 2011 for the treatment and prevention of vitamin D deficiency, concluded that to guarantee bone health, a blood level of 25(OH)D of at least 75 nmol/L (30 ng/mL) is required (https://www.endocrine.org/clinical-practice-guidelines/vitamin-d-deficiency)[12]. Beyond musculoskeletal health however, it has been found that vitamin D supplementation significantly reduced the risk of cancer death by 15% in a systematic review and meta-analysis of 52 trials with a total of 75454 participants[13], and it has been suggested that better health outcomes may occur in the range of 100-150 nmol/L[10]. The largest meta-analysis ever conducted of all studies published between January 1, 1966 and January 15, 2013 dealing with all-cause mortality related to serum 25(OH)D showed that 25(OH)D < 75 nmol/L was associated with higher all-cause mortality, its reduction being maintained with 25(OH)D ≥ 175 nmol/L (70 ng/mL), without a U-shaped curve as previously reported[10]. Achieving such concentrations with supplements and sensible sun exposure for a normal weight adult requires 2000–5000 IU daily intake of vitamin D2/D3, practically all year long except maybe during sunny vacations[14]. With vitamin D adequacy relying mainly (80%-90%) on sun exposure rather than on dietary sources (10%-20%), if not on supplementation, these doses should be adapted accordingly during lockdowns. It should also be recognized that sensible sun exposure has many additional health benefits not only in the immune system but also in improving the feeling of well-being[15]. At this time, neither the World Health Organization nor any other public health authority has issued any official advice or recommendation on vitamin D, or any other nutrients, to the best of our knowledge.

A quadratic relationship was found between vitamin D deficiency in countries affected by COVID-19 and the latitudes, implying a possible relation[16]. When mortality/ million (M) is plotted against latitude, all countries below 35 degrees North, above which people do not receive sufficient sunlight to retain adequate 25(OH)D concentrations during winter, have relatively lower mortality, implying a role for vitamin D status in outcomes from COVID-19[17]. Vitamin D is strongly affected by ozone variability, since ozone filters ultraviolet B, an important factor for vitamin D synthesis. A statistically significant link between ozone concentration and incidence of COVID-2019 disease in 34 countries was established[18]. Going back to the 1918-1919 influenza pandemic, substantial correlations were found for associations of July ultraviolet B dose in the United States with case fatality rates and rates of pneumonia[19]. Α significant association between vitamin D sufficiency and reduction in clinical severity and inpatient mortality was very recently shown[20]. Thus, to elucidate further the possible role of vitamin D population status in the COVID-19 pandemic, we examined the associations between published representative and standardized population vitamin D data on European population vitamin D status and the Worldometer COVID-19 data.

MATERIALS AND METHODS

Accessing data on European countries at the Worldometer, on June 19, 2020, we analyzed the 28 countries populated > 4 M (Table 1). For months, Swedish public health authorities have defended their controversial decision not to lock down the country in response to the global COVID-19 pandemic, with the country experiencing dramatic casualties. Thus, Sweden was excluded from analysis. The remaining 27 European countries adopted a defensive strategy during the current pandemic, even with delays and hesitations, as in the United Kingdom. Moldova was also excluded as no published vitamin D status data were found. For the remaining 26 countries, we used linear regression to explore the correlation between reported representative and standardized population vitamin D concentrations[21-28] and the number of total cases/M and recovered/M until June 19, 2020 as well as the deaths/M and the serious-critically ill/M from COVID-19 on that date (Table 1). Since mortality of COVID-19 disease has been shown to increase rapidly in respect to age, life expectancy (LE), an age-related index, was analyzed using a semi-parametric regression approach using Worldometer data. Weighted (https://doi.org/10.13094/SMIF-2015-00001) analysis of variance (ANOVA)/analysis of covariance (ANCOVA) was performed to evaluate serious-critical/M and deaths/M by the vitamin D population status - categorized as deficient (D) < 50, insufficient (IN) 50-62.5, mildly insufficient (MIN) > 62.5-75 and sufficient (S) > 75 nmol/L – while controlling for LE for deaths/M. To test whether these correlations withstand at another completely different momentum of this pandemic, which would be an indication of a truthful association, although still not a proof of causality, we also checked the above correlations and the differences between consecutive points of the same variables on August 15, 2020. All statistical analyses were performed in XLSTAT LIFE SCIENCE version April 1, 2020 (copyright Addinsoft 1995-2020) and R (R Core Team 2017), with the use of the SemiPar library.

Table 1 European coronavirus disease 2019 data from the Worldometer on June 19, 2020, compared to life expectancy and to available representative and standardized data on the European population vitamin D status[21-28].

	Country	Total cases/M	Total recovered	Serious critical	Deaths/M	Life expectancy in yr	Population 25(OH)D in nmol/L	Population, M
1	Russia	3899	324406	2300	54	72.99	39.7	145.93
2	Germany	2273	174400	396	107	81.88	50.1	83.77
3	United Kingdom	4447	N/A	379	626	81.77	47.4	67.87
4	France	2431	73887	752	454	83.13	60.0	65.26
5	Italy	3939	180544	168	571	84.01	45.0	60.46
6	Spain	6253	N/A	617	580	83.99	59.9	46.75
7	Ukraine	800	16033	343	23	72.50	29.0	43.74
8	Poland	827	15698	87	35	79.27	32.0	37.84
9	Romania	1216	16555	184	77	76.50	65.0	19.24
10	Netherlands	2885	N/A	57	355	82.78	64.7	17.13
11	Belgium	5219	16751	55	837	82.17	49.3	11.58
12	Czechia	968	7472	9	31	79.85	62.5	10.70
13	Greece	311	1374	10	18	82.80	54.3	10.42
14	Portugal	3772	24477	67	150	82.65	55.4	9.66
15	Sweden	5550	N/A	272	500	83.33	68.7	9.44
16	Hungary	422	2581	15	59	77.31	60.6	83.33
17	Belarus	6067	35275	92	36	75.20	72.0	9.00
18	Austria	1918	16141	7	76	82.05	51.7	8.73
19	Serbia	1454	11511	18	30	76.47	65.7	8.65
20	Switzerland	3,608	28900	17	226	84.25	46.0	6.94
21	Bulgaria	529	1941	13	27	75.49	38.7	5.79
22	Denmark	2139	11282	6	104	81.40	65.0	5.54
23	Finland	1287	6200	2	59	82.48	67.7	5.45
24	Slovakia	289	1447	0	5	78.00	81.5	5.41
25	Norway	1,609	8138	5	45	82.94	71.0	4.93
26	Ireland	5137	22698	28	347	82.81	56.4	4.10
27	Croatia	555	2142	0	26	79.02	46.9	4.03
28	Moldova	3249	7525	455	111	72.30	N/A	10.09

25(OH)D: 25-hydroxyvitamin-D; M: Million.

RESULTS

From the 26 European countries included in the analysis, populated 714.661 M in total, nine (54.17%, 387.15 M) had a vitamin D deficient status, eight an insufficient status (33.58%, 240.022 M), eight a mild insufficiency status (11.48%, 82.023 M) and only one country, Slovakia, a sufficient status (0.76%, 5.459 M). There was no correlation between the total cases/M nor the recovered/M and the European population vitamin D concentrations. Negative correlations were recognized regarding the total deaths/M (Figure 1A), predicting a reduction of deaths/M by 20% if the 25(OH)D concentration reaches 50 nmol/L (related to the number calculated at 25), by 40% at 75, by 61% at 100, by 82% at 125 and by 102.4% at 150 nmol/L and the serious-critical/M (Figure 1B), predicting a reduction of serious-critically ill/M by 16% if 25(OH)D concentration reaches 50 nmol/L (related to the number calculated at 25), by 31% at 75, by 47% at 100, by 64% at 125 and by 80% at 150 nmol/L.

Figure 1 Linear regression on June 19, 2020 related to available representative and standardized data on the European population vitamin D concentrations (x axis, nmol/L). A: Of the total deaths/million (M); B: Of the serious-critical cases/M.

Population vitamin D concentrations vs life expectancy exhibits a non-linear relationship (Figure 2): Higher life expectancy until approximately 77 years of age is characterized by better vitamin D concentrations, while practically reaching a plateau at 82 years, and then by a decline as expected in the elderly. There is a non-linear relationship between life expectancy and deaths/M with a dramatic increase in deaths/M after approximately 80 years (Figure 3). LE (i.e. age) seems to interfere with the effect of a better vitamin D concentration to the total number of deaths/M, rendering the vitamin D benefit even more important than the unadjusted one: A reduction in total deaths/M by 44% if 25(OH)D concentration reaches 50 nmol/L (related to the number calculated at 25), by 88% at 75, by 133% at 100, by 177% at 125 and by 221% at 150 nmol/L. The analytical form for the model on the deaths/M accounting for a potential non-linear effect of LE is year = -2675-4.111*vitamin D + f(LE), where f(.) is a non-linear smooth function of life expectancy. The P value for the term f(LE) was estimated via likelihood ratio test to be P = 0.042, indicating a statistically significant effect of life expectancy on deaths/M after adjusting for vitamin D concentration.

Figure 2 Population vitamin D concentrations vs life expectancy exhibits a non-linear relationship.

Figure 3 Non-linear relationship between life expectancy and deaths/million.

Weighted (https://doi.org/10.13094/SMIF-2015-00001) ANOVA was performed to evaluate serious-critical/M and ANCOVA for deaths/M by the population vitamin D status while controlling for LE. Given the r², about 22% of the variability of the dependent variable serious-critical/M could be explained by the population vitamin D status. A decreasing trend from population status D [β = 8.684, standard error (SE) = 2.196, 95% confidence interval (CI): 4.372/12.996, P < 0.001], IN (β = 7.883, SE = 2.205, 95%CI: 3.553/12.213, P < 0.001), MIN (β = 4.548, SE = 2.252, 95%CI: 0.126/8.169, P = 0.044) to S (LE mean 0.0, SE 2.181, 95%CI: -4.282/4.282, P < 0.001) was found with an average reduction of serious-critical/M of 9.2% from vitamin D status deficient to insufficient, of 47.6% from deficient to mildly insufficient and 100% from deficient to sufficient (reference, Figure 4). Regarding deaths/M (Figure 5), given the r², about 63% of the variability of the dependent variable deaths/M could be explained by the two variables, LE alone accounting for 47%. A decreasing trend from population status deficient (β = 150.375, SE = 8.859, 95%CI: 132.982/167.768, P < 0.001), insufficient (β = -72.514, SE = 10.336, 95%CI: -150.170/-55.866, P < 0.001), mildly insufficient (β = -80.518, SE = 12.556, 95%CI: -105.170/-55.866, P < 0.001) to sufficient (β = -129.122, SE = 62.915, 95%CI: -252.644/-5.599, P = 0.041) was found with an average reduction of deaths/M of 62.9% from vitamin D status deficient to insufficient, of 65.15% from deficient to mildly insufficient and 78.8% from deficient to sufficient.

Figure 4 Least square means of serious-critical/million for factor population vitamin D status. D: Deficiency; IN: Insufficiency; MIN: Mild insufficiency; S: Sufficiency.

Figure 5 Least square means of deaths/million for factor population vitamin D status. D: Deficiency; IN: Insufficiency; MIN: Mild insufficiency; S: Sufficiency.

On August 15, 2020, the above correlations were sustained and the differences between consecutive points for the two variables serious-critical/M and deaths/M in the two time points were correlated, not proving causality but suggesting a truthful association.

DISCUSSION

We explored any possible correlation between the population vitamin D status - influenced by various factors - and COVID-19 disease, in particular total cases, serious-critical illness and deaths. In contrast to a recently published study[29], we found no association between the vitamin D status of the European populations and the total confirmed cases/M of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections when we analyzed data from the Worldometer on June 19, 2020 on 26 European countries populated > 4 M. However, the negative correlations that we found between population vitamin D status and serious-critical/M and deaths/M show a clear tendency, even if they do not prove causality, namely after adjusting for LE, underlining the importance of an optimal vitamin D status especially in the elderly[30]. On August 15, 2020, at a completely different time point of this pandemic, before the second wave even had started, the above associations were sustained, suggesting a truthful correlation. Since the risk of COVID-19 disease increases rapidly with respect to age, an age-related index, such as LE, was found, as expected, to be a more important predictor of death rates. Thus, according to our results, a higher 25(OH)D concentration may protect from serious-critical illness and death from COVID-19 disease even more in the elderly but does not seem to prevent SARS-CoV-2 from spreading, in contrast to a recent study[29], which however reported also a negative correlation between the mean population vitamin D concentrations of 20 European countries and deaths/M from COVID-19 on April 8, 2020. Our findings also coincide with a recent study from Maghbooli et al[20] showing that vitamin D sufficiency [a serum 25(OH)D > 75 nmol/L (30 ng/mL)] reduced risk for adverse clinical outcomes in patients with COVID-19 infection: 6.3% of the patients who had a blood 25(OH)D concentration of at least 100 nmol/L (40 ng/mL) succumbed to the infection compared to 9.7% and 20% who died and had a circulating blood level above and below 75 nmol/L (30 ng/mL), respectively[20,31], suggesting that a blood level of at least 100 nmol/L (40 ng/mL) may be optimal for obtaining vitamin D’s immunomodulatory benefit.

Various parameters played a significant role in the spread of the current pandemic. Among them, air travel and direct connections with China and particularly Wuhan, where the epidemic started. Then, health policymaking with mass quarantine was instituted in most countries, influencing the course of the disease, but with no central coordination of the measures taken during the first wave of the pandemic, not even in the core of the European Union itself. Timing of the lockdowns, at least in the first wave, seemed to have been the main factor affecting the number of the cumulative deaths – although this has been strongly debated (https://thefatemperor.com/published-papers-and-data-on-lockdown-weak-efficacy-and-lockdown-huge-harms/), along with travel and border restrictions. Recent research emphasizes the importance of face masks while self-protection measures seem to be better implemented by populations with higher educational levels. Temperature also appears to have a small but statistically significant impact on the viral transmission rate as countries with daily average temperatures below 20 °C had a faster transmission rate. Most probably, genetic predisposition must have played a fundamental role in the susceptibility in SARS-CoV-2 infection[32,33]. The recent discovery of robust genetic signals relating to key host antiviral defense mechanisms and mediators of inflammatory organ damage in COVID-19 may lead to targeted treatment with existing drugs[33]. Most recent evidence show that angiotensin-I converting enzyme-2 (ACE2) expression and/or polymorphism could also influence both the individual susceptibility to SARS-CoV-2 infection and the outcome of the COVID-19 disease[34]. Thus, the integrity of our immune system and its ability to fight back with a coordinated way, keeping asymptomatic or within the subclinical spectrum most of the people infected and saving the lives of the severely infected, is a crucial factor. And there is significant evidence that vitamin D deficiency may compromise both innate and acquired immunity responses, leading to increased vulnerability to infections as to autoimmune responses and disorders[35].

The vitamin D status of a population is dependent on a variety of factors including supplementation and food fortification strategies, latitude of the country, season as well as on the local nutritional and sun exposure habits, especially in the non-institutionalized elderly[36]. The vitamin D status in the winter is even lower[1,37,38], with underappreciated consequences to the immune function[39,40]. Ideally, we should be able to analyze data on vitamin D status of the elderly in winter. Thus, a major limitation of our ecological approach is that we had to rely on published - but perhaps not always completely representative - data on the vitamin D status of the populations in Europe. However, data analyzed are based mainly on “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” recently published in the European Journal of Endocrinology[21] – presenting not only representative nationally or regionally as possible but also standardized population vitamin D concentrations -, a systematic review of vitamin D status in southern European countries[22], and a very important study applying the protocols developed by the National Institutes of Health-led international Vitamin D Standardization Program to serum 25(OH)D data from representative childhood/teenage and adult/older (we chose data from older adults) European populations, representing a sizable geographical footprint, to better quantify the prevalence of vitamin D deficiency in Europe[28]. Keeping in mind that the population vitamin D status reflects that of the elderly, which by default will be worse, we tried to analyze the most recently validated and representative data possible, whereas from the available data for each country we chose data from older adults in winter where provided, and in any case from Caucasian descent. Ideally, we should be able to analyze data on 25(OH)D concentrations of the patients as in an interesting recent report from Switzerland, which found significantly lower circulating 25(OH)D concentrations [27.75 nmol/L (11.1 ng/mL), P = 0.004] in polymerase chain reaction-positive for SARS-CoV-2 patients compared with negative patients [61.5 nmol/L (24.6 ng/mL)], even after stratifying patients according to age > 70 years[41]. Another important issue would be the differences in assessment mainly of the COVID-19 deaths in the various European countries. However, the World Health Organization had already issued the “International guidelines for certification and classification (coding) of COVID-19 as cause of death, April 20, 2020” 2 mo earlier to our analysis, allowing us to assume that they must had already been adopted by the European Countries responsible public health authorities. Furthermore, Worldometer.info mainly collects data from official reports, directly from governmental communication channels. An additional important limitation is the true evaluation of the number of affected subjects in the variable countries: Since not all patients infected with COVID-19 are symptomatic, the cases/M are dependent upon the percentage of the population tested and the consistency of the frequency of testing during the disease period evaluated, not to mention that several patients or carriers have been tested several times. Furthermore, the definition of case includes a carrier as well as a patient. Unfortunately, this limitation could not be overcome with the publicly available COVID-19 data at the time of our analysis. However, we had to report the absence of any correlation between total cases/M and the population vitamin D status in the sample we analyzed, in contrast to a recently published study with the opposite results[29]. Assessment of serious-critical cases in the European countries may also have been limited at some points by the shortcoming of intensive care unit (ICU) beds as well as the introduction of different drugs and “cocktail” treatments from country to country. Albeit, on June 19, 2020, the first wave of the pandemic in Europe was kind of winding down, and not particularly effective new or repurposed medication had at least been qualified at that point as such to change significantly the clinical course of the serious-critical patients, other than the accumulated experience of the health workers fighting on the frontline.

Independent researchers increasingly call for optimization of vitamin D status for enhanced immune protection against COVID-19 at least in older adults, hospital inpatients, nursing home residents and other vulnerable groups, extending this recommendation to the general population[42]. The elderly (> 65 years) have a higher risk for vitamin D deficiency due to decreased sun exposure and reduced ability for cutaneous synthesis[38], whereas aging exerts significant effects on all cells of the innate immune system[40], making vitamin D sufficiency even more valuable in this group. Early nutritional supplementation in non-critically ill patients hospitalized for COVID-19 has been implemented in hospital protocols providing 50000 UI/wk if 25(OH)D < 50 nmol/L and 25000 UI/wk if 25(OH)D < 75 nmol/L aiming at improved immunologic recovery with reduced levels of inflammation, immune activation, and increased immunity against pathogens[43].

The COVID-19 pandemic presents a puzzling challenge without specific treatment yet with timely administration being crucial for all current regimens on clinical trial or use. This is also the case for vitamin D, and this might be the reason why in a recent RCT, a single enteral dose of 540000 IU of vitamin D3 or matched placebo started late within 12 h after the decision to admit the critically ill (unrelated to COVID- 19) vitamin D deficient patient to an intensive care unit, had no benefit at a 90-d all-cause, all-location mortality[44]. Regarding vitamin D, we know that respiratory viruses downregulate vitamin D receptor expression in human bronchial epithelial cells, while improvement in vitamin D status increases antiviral defenses via cathelicidins and innate interferon pathways[45]. Vitamin D has a 12% overall protective effect against bacterial and viral acute respiratory tract infection, increased to 19% in those individuals on daily or weekly regimen compared to those on monthly boluses and up to 70% when vitamin D deficiency is corrected with daily supplementation[46]. Bioavailable 25(OH)D is inversely associated with illness severity in critically ill ICU patients associated with increased mortality and morbidity[47]. Calcitriol [1,25(OH)₂D₃] alleviates lipopolysaccharide induced acute lung injury and prevents the adult respiratory distress syndrome by minimizing the alveolar damage[48]. Vitamin D is also a negative endocrine regulator of the renin-angiotensin system. The mechanism for SARS-CoV-2 infection is the requisite binding of the virus to the membrane-bound form of ACE2 and internalization of the complex by the host cell. Recognition that ACE2 is the main host receptor by SARS-CoV-2 to infect human has prompted new therapeutic approaches to block the enzyme or reduce its expression to prevent cellular entry of SARS-CoV-2 in tissues expressing ACE2 (lung, heart, kidney, brain, and gut). Thus, it seems that both stimulation of the immune system and inhibition of renin-angiotensin system are mechanisms by which vitamin D may play a beneficial role in COVID-19 infection[49]. Vitamin D repletion in critical illness with a more aggressive dosing is showing similarly promising results with vitamin C repletion in septic shock[50] and may be able to prevent the cytokine storm that seems to be killing people rather than the virus itself[51]. C-reactive protein is a surrogate marker for unregulated inflammation and cytokine storm and is associated with vitamin D deficiency. Retrospective data and indirect evidence also show a possible role for vitamin D in reducing complications attributed to and the cytokine storm itself[52]. Moreover, recent research revealed that vitamin D receptor signaling in macrophages regulates a shift between proinflammatory and anti-inflammatory activation during ER stress-induced inflammation[53]. Thus, supplementation within recommended upper safety limits, for specific nutrients such as vitamins C and D, warrants optimal nutritional status to insure a well-functioning immune system protecting against viral infections[54].

Recent research demonstrated that low serum 25(OH)D was an independent risk factor for COVID-19 infection and hospitalization analyzing data from 14,000 members of Leumit Health Services in Israel[55]. A very recent pilot randomized clinical study demonstrated that administration of a relatively high dose (0.532 mg-21280 IU) of calcifediol or 25(OH)D, followed by 0.266 mg on days 3 and 7, and then weekly until discharge or ICU admission, significantly reduced the need for ICU treatment of patients requiring hospitalization due to proven COVID-19 disease[56]. In a single-center, retrospective cohort study concerning 489 patients, likely deficient vitamin D status was associated with increased COVID-19 risk[57].

Our analysis took place at two completely different time points during the beginning and the end of the first wave of this pandemic. We needed to confirm our first results at a completely different time point of the first wave. One could not attempt to extend this type of approach to the second wave or the third wave, which is now hitting Europe, first because the virus has significantly spread into the European populations. Secondly, after extended lockdowns and limited - if any - summer vacations and with no public health authority having officially advised supplementation with vitamin D even aiming to protect musculoskeletal health, one can hypothesize that the European population vitamin D status had to be worse, and this could also be one of the main reasons why the second and third waves appeared more deadly than the first, at least in several European countries, even in sunny countries as Greece. One of the main outcomes of our analysis though is that vitamin D does not prevent SARS-CoV-2 from spreading, while it may protect from serious-critical illness and death form COVID-19 disease, with significant and substantial protection being obtained at a 25(OH)D concentration of 100-150 nmol/L (40-60 ng/mL).

CONCLUSION

At this time and despite the ongoing debate on “The Big Vitamin D Mistake”[15], referring to the statistical error in the estimation of the Recommended Dietary Allowance of vitamin D discovered by Veugelers and Ekwaru[58] in 2014 and confirmed by Heaney et al[59]: About 4000 IU/d (3385) are needed to ensure 50 nmol/L in 97.5% of the population, about 6000 IU/d (6201) are needed to achieve the Endocrine Society’s recommendation of 75 nmol/L and about 9000 IU/d (9122) to reach 100 nmol/L, and even if the vitamin D deficiency pandemic is still being questioned[60], no one should confuse the global consensus on the minimum vitamin D doses needed to prevent nutritional rickets[61], with the doses needed to exert all of its extra-skeletal health benefits[62], particularly those related to our immune system. Apart from the known disagreement between the Endocrine Society and the Institute of Medicine (IOM) but also the discrepancy between the IOM and the Scientific Advisory Committee on Nutrition in Great Britain, two equally respectable government advisory committees, who after reviewing the same evidence, ended up with a twofold difference in target concentrations in serum 25(OH)D and similarly divergent conclusions for intakes of vitamin D[12], one can notice that differences concerning upper tolerable limits for vitamin D administration are limited. The more conservative IOM advises up to (upper tolerable limit) 1500 IU daily in infants < 1 year, 2500 IU in children 1-3 years, 3000 IU in children 4-8 years and up to 4000 IU for everybody after 9 years of age; where the Endocrine Society advices are up to 2000 IU for infants < 2 years, up to 4000 IU for children 1-18 years and up to 10000 IU for adults, adult pregnant and lactating women as well as the elderly, underlining that obese people may need up to two to three times more, as it may be needed to correct vitamin D deficiency or to treat specific conditions such as rickets, osteomalacia, hyperparathyroidism, malabsorption syndromes or if on medications affecting vitamin D’s metabolism. However, the doses that the Endocrine Society practice committee characterizes as not requiring medical supervision are practically identical to the IOM’s upper tolerable limits. Thus, supplementation with vitamin D within recommended safety limits, with doses that do not require prior measurement of the 25(OH)D concentration or medical supervision, apart from the already established protective role in bone mineral density[63], may also assure a well-functioning immune system[64].

In 2011, the Endocrine Society published the Endocrine Society Practice Guidelines on vitamin D, recommending how to treat and prevent vitamin D deficiency in children and adults. Based on the literature these recommendations were related to maximizing musculoskeletal health. However, in 2011 there was not enough scientific evidence for the Committee to recommend improvement in vitamin D status for reducing risk of many chronic illnesses or improving immune function. During the past decade, however, numerous studies have been conducted demonstrating that improvement in vitamin D status reduces risk for upper respiratory tract viral infections as well as having a wide variety of effects on both innate and acquired immunity[39,65]. A recent randomized controlled double-blind clinical trial assessed the impact of vitamin D supplementation on calcium metabolism and non-calcemic broad gene expression by relating them to the individual’s responsiveness to varying doses of vitamin D3[66]. Thirty healthy adults were randomized to receive 600, 4000 or 10000 IU/d of vitamin D3 for 6 mo. Circulating parathyroid hormone (PTH), 25(OH)D, calcium and peripheral white blood cells broad gene expression were evaluated. The investigators reported dose-dependent increase in circulating 25(OH)D concentrations, decreased PTH concentrations and no change in serum calcium levels. A plateau in circulating PTH levels was achieved at 16 wk in the 4000 and 10000 IU/d groups. There was a dose-dependent 25(OH)D alteration in broad gene expression with 162, 320 and 1289 genes up- or down-regulated in their white blood cells, respectively. Thus, improvement in vitamin D status does have a dramatic effect on immune cell activity. However, can it therefore be expected that everyone who improves their vitamin D status would experience the same genomic influences on their immune system if they raised their blood level of 25(OH)D to the same degree? Carlberg and Haq[67] gave daily 3200 IU of vitamin D3 to 71 prediabetic patients for 5 mo and found robust changes in total gene expression in peripheral blood mononuclear cells only in about half the subjects. Shirvani et al[66] observed in healthy adults who were vitamin D deficient and who received this same dose of vitamin D and raised their blood concentrations of 25(OH)D to the same degree, marked differences in the level of expression of the same genes. They reported that 60% of the healthy vitamin D deficient adults who received 10000 IU daily for 6 mo had a robust response in gene expression compared to the other 40% who had minimum to modest responses even though these subjects raised their blood concentrations of 25(OH)D in the same range of 60-90 ng/mL (150-225 nmol/L).

With all of this compelling information, it is reasonable for all responsible Public Health Authorities to consider advising their populations to enhance their immune system by improving their vitamin D status by encouraging sensible sun exposure and by taking vitamin D supplements (if not already on adequate supplementation or medically prohibited due to a vitamin D hypersensitivity disorder) at the doses which, as proposed by the Endocrine Society Guideline Committee in 2011, do not require previous laboratory testing nor medical supervision. To prevent nutritional rickets, daily doses of 400-1000 IU in infants, 600-1000 in children and 1500-2000 in teenagers (should be treated as adults) and adults, are needed. However, to achieve higher circulating concentrations of 25(OH)D at the range of 100-150 nmol/L (40-60 ng/ml), appearing according to our analysis to be necessary for substantially improving immune function and protect from COVID-19 disease, without any risk of toxicity[68], higher doses can be used. As mentioned above, the Endocrine Society Practice Guidelines recommends the safe upper limit for infants < 1 year is 2000 IU daily, children 1-18 years 4000 and adults (including elderly and adult pregnant-lactating women) 10000 IU, unless they are obese, requiring two to three times more. Thus, after a necessary initial repletion for up to 2 mo with these upper tolerable doses, the Endocrine Society’s Committee’s maintenance proposed doses, which can be safely given without medical supervision to prevent vitamin D deficiency and are practically identical with the IOM’s upper tolerable limits, can be continued: i.e. up to 1000 IU/d for infants aged < 6 m, 1500 for age 6 m - 1 year, 2500 for 1-3 years, 3000 for children 4-8 years and 4000 for children > 8 years, with adults, pregnant/lactating women and adolescents requiring a daily intake of 4000-5000 (8000-10000 if obese) to maintain circulating concentrations of 25(OH)D at the range of 100-150 nmol/L. For teenagers and adults on a weekly scheme, these doses translate to about 50000 or if obese 100,000 IU, this being equivalent to approximately 6000 IU daily and 12000 IU for obese, respectively.

These doses will achieve blood concentrations of 25(OH)D of at least 75 nmol/L (30 ng/mL) aiming at the preferred range of 100-150 nmol/L (40-60 ng/mL), without any risk of toxicity[68]. It has been estimated that once a blood concentration of 25(OH)D reaches 50 nmol/L (20 ng/mL) that for every 100 IU ingested, the blood concentration increases by approximately 0.6-1 ng/mL. A good example of this dosing was reported by Shirvani et al[66] who demonstrated that circulating concentrations of 25(OH)D were maintained in the range of 24.3 ± 4.1, 40.8 ± 3.8 and 78.6 ± 13.5 ng/mL, in vitamin D deficient adults who ingested 600, 4000 and 10000 IU daily for 6 mo. These data are supported by a population based Canadian study demonstrating that some adults taking up to 20000 IU daily for more than a year maintained a blood concentration of 25(OH)D in the range of 60-80 ng/mL without any evidence of toxicity[69]. This study also nicely demonstrated the effect of BMI on vitamin D status. The authors observed that those who had a BMI > 30 kg/m² needed to ingest 2.5 times more vitamin D to maintain the same blood level as a normal weight adult.

Achieving circulating concentrations of 25(OH)D in the range of 100-150 nmol/L (40-60 ng/mL) appears to optimize vitamin D’s effect on improving immune function, thereby substantially reducing the risk for serious-critical infections, particularly from SARS-CoV-2 according to our study, and possibly modulating the immune response, helping to prevent the dangerous cytokine storm often leading to COVID-19 related deaths. The COVID-19 pandemic is an unprecedented medical emergency for the modern world, and we may not possess the luxury, the time nor even the ethical argument to wait the definite results on RCTs while people are dying[70], while prospective well designed studies are needed to conclude on the impact of the vitamin D status on COVID-19 morbidity and mortality[71]. These trials are hopefully awaited, but before a medical emergency of this magnitude we need to remember that Evidence Based Medicine is not necessarily synonymous to RCTs. We do know that vitamin D enhances immune function. We know the extent of vitamin D deficiency, and we know that restrictions and lockdowns have probably worsened the populations’ vitamin D status. Thus, until then, decisions are taken based on and adapted to the best available evidence. And, as far as vitamin D, the evidence is there[51], justifying even the use of vitamin D as a possible adjuvant therapy for COVID-19 disease[72]. A preponderance of evidence does suggest that vitamin D deficiency increases mortality. Our findings predict a striking reduction of serious-critical illness and deaths from COVID-19 if 25(OH)D concentrations reach 100-150 nmol/L (40-60 ng/ml), and very recently SARS-CoV-2 positivity was found to be strongly and inversely associated with circulating 25(OH)D concentrations irrespective of latitudes, races/ethnicities, both sexes and age ranges[73]. Slovakia, at five deaths/M, having the lowest mortality rate in Europe from COVID-19 disease at the time of our analysis, a 125-fold lower than in the UK where official advice remains that 25(OH)D deficiency is < 25 nmol/L (https://www.nice.org.uk/advice/es28/evidence/evidence-review-pdf-8777674477), is a characteristic paradigm, being practically the only country in Europe with a 25(OH)D status meeting the Endocrine Society’s recommended level of sufficiency > 75 nmol/L (30 ng/mL).

From a public health perspective, given the established safety of even high doses, and the potential benefits in enhancing innate and adaptive immunity[74], mitigating also the inflammatory response[3], the recommendation of intensive supplementation with vitamin D as possible prophylaxis with safe doses that do not require prior measurement or medical supervision, must be seriously considered, especially now that the world is facing the third deadly wave of this pandemic, forcing populations into repeated new lockdowns without the broad availability of specific medications yet and while awaiting for vaccinations to be widely available and plausible.

There is no need to require a measurement of serum 25(OH)D before recommending treatment and/or supplementation with vitamin D. This is supported by the observation that ingesting 50000 IU of vitamin D every 2 wk for up to 6 years is not associated with any toxicity[75]. Furthermore, this study was conducted in a clinical setting and all patients were prescribed this vitamin D therapy without the knowledge of their baseline serum 25(OH)D concentration. After completion of the study, the baseline levels were measured. Some of the patient’s had a blood concentration of 25(OH)D of 125 nmol/L (50 ng/mL) and after being on 50000 IU of vitamin D once every 2 wk, their 25(OH)D concentration reached 200 nmol/L (80 ng/mL) without any evidence of toxicity[75].

There is essentially no vitamin D naturally occurring in the diet apart from oily fish, cod liver oil and sun-dried mushrooms. The modern way of life deprives us from sun exposure together with the warning to avoid all direct sun exposure by the national and international Dermatology Societies contributing to the worldwide vitamin D deficiency pandemic: Approximately 40% of the world’s population is vitamin D deficient, i.e. 25(OH)D < 50 nmol/L (20 ng/mL) and 60% or insufficient i.e. 50-79 nmol/L (20-29 ng/mL). Therefore, we also need to consider worldwide recommendations for vitamin D food fortification that is practiced in several countries including the United States, Canada, and Finland to name a few. Most other countries either do not encourage or forbid food fortification with vitamin D. Recently, in 2017, India implemented fortification of milk and cooking oil with vitamin D2 as a means of reducing vitamin D deficiency that is common in both children and adults in this sunny Asian subcontinent.

Vitamin D is safe, not toxic and inexpensive. In the “shade” of the modern way of life, the human body cannot produce enough vitamin D from sun exposure, as our hunter gatherer forefathers did and as Maasai herders and the Hazda continue to do. Vitamin D may improve and modulate immune response against SARS-CoV-2. With all the above data, the limitations and the perspectives discussed, the possible benefit in the fight against SARS-CoV-2 should the protection against COVID-19 serious-critical illnesses and death with vitamin D prove truthful, and this without any risk of toxicity, the gain for humanity as well global public health might be just invaluable.

ARTICLE HIGHLIGHTS

Research background

Recent studies have claimed lower coronavirus disease 2019 (COVID-19) cases in European countries with a better vitamin D status and a significant association between vitamin D sufficiency and reduction in clinical severity and inpatient mortality from COVID-19 disease. Low serum 25(OH)D was identified as an independent risk factor for COVID-19 infection and hospitalization, and administration of calcifediol or 25(OH)D significantly reduced the need for intensive care unit treatment.

Research motivation

Vitamin D population status may indeed have possible unappreciated consequences to the COVID-19 pandemic, a hypothesis that needed to be further elucidated.

Research objectives

Following an ecological integrative approach, we examined the associations between published representative and standardized European population vitamin D data and the Worldometer COVID-19 data at two completely different time points of the first wave of this pandemic. If any sustained correlations were to be found, they would be an indication of a truthful association, even though they could not prove causality.

Research methods

Using linear regression, we explored the correlation between published representative and standardized population vitamin D concentrations and the number of total cases/million (M), recovered/M, deaths/M and serious-critically ill/M from COVID-19 for 26 European countries populated > 4 M. Life expectancy (LE) was also analyzed with semi-parametric regression. Weighted analysis of variance/analysis of covariance evaluated serious-critical/M and deaths/M by the vitamin D population status: deficient < 50, insufficient: 50-62.5, mildly insufficient > 62.5-75 and sufficient > 75 nmol/L, while controlling for LE for deaths/M. Statistical analyses were performed in XLSTAT LIFE SCIENCE and R (SemiPar library).

Research results

No correlation was found between population vitamin D concentrations and the total cases-recovered/M, but negative correlations were depicted predicting a reduction of 47%-64%-80% in serious-critical illnesses/M and of 61%-82%-102.4% in deaths/M, further enhanced when adapting for LE by 133%-177%-221% if 25(OH)D concentrations reach 100-125-150 nmol/L. Weighted analysis of variance evaluated serious-critical/M (r² = 0.22) by the vitamin-D population status and analysis of covariance the deaths/M (r²= 0.629) while controlling for LE (r² = 0.47). Serious-critical showed a decreasing trend (P < 0.001) from population status deficient (P < 0.001) to insufficient by 9.2% (P < 0.001), to mildly insufficient by 47.6% (P = 0.044) and to sufficient by 100% (reference, P < 0.001). For deaths/M the respective decreasing trend (P < 0.001) was 62.9% from deficient to insufficient (P < 0.001), 65.15% to mildly insufficient (P < 0.001) and 78.8% to sufficient (P = 0.041).

Research conclusions

A higher 25(OH)D concentration may protect from serious-critical illness and death from COVID-19 disease - even more in the elderly - but does not seem to prevent severe acute respiratory syndrome coronavirus 2 from spreading.

Research perspectives

Considering the ongoing pandemic situation, the presented results are useful for public health systems to advise their populations to enhance their immune system by improving their vitamin D status. Specifically, achieving a serum 25(OH)D concentration of 100-150 nmol/L (40-60 ng/mL) with vitamin D2/D3 supplementation using the upper tolerable daily doses for up to 2 mo (infants < 1 year 2000 IU daily, children 1-18 years 4000 and adults including elderly and adult pregnant-lactating women 10000 IU, unless they are obese requiring 2-3 times more) followed by the maintenance proposed doses not requiring medical supervision, as proposed by the Endocrine Society and being practically identical with the Institute of Medicine’s upper tolerable limits (up to 1000 IU/d for infants aged < 6 mo, 1500 for age 6 mo - 1 year, 2500 for 1-3 years, 3000 for children 4-8 years and 4000 IU for children > 8 years, with adults, pregnant-lactating women and adolescents requiring a daily intake of 4000-5000 unless they are obese requiring two to three times more) may protect from serious-critical illness and death from COVID-19 disease.

ACKNOWLEDGEMENTS

We thank the experts in biostatistics Alexandros Gryparis and Arash Shirvani for guiding us in performing the statistical analysis.