key: cord-0892972-hxfdobzk authors: Infante, Marco; Buoso, Andrea; Pieri, Massimo; Lupisella, Santina; Nuccetelli, Marzia; Bernardini, Sergio; Fabbri, Andrea; Iannetta, Marco; Andreoni, Massimo; Colizzi, Vittorio; Morello, Maria title: Low Vitamin D Status at Admission as a Risk Factor for Poor Survival in Hospitalized Patients With COVID-19: An Italian Retrospective Study date: 2021-02-18 journal: Journal of the American College of Nutrition DOI: 10.1080/07315724.2021.1877580 sha: aebd11b6c0cb82a2e5fc201f10dc79e2dcc1bb36 doc_id: 892972 cord_uid: hxfdobzk OBJECTIVE: Preliminary findings suggest a relationship between lower serum 25-hydroxyvitamin D [25(OH)D] levels and incidence and severity of COVID-19. The aim of this study was to evaluate the relationship between vitamin D status at admission and different markers of inflammation, coagulation, and sepsis in hospitalized patients with COVID-19. METHOD: We conducted a retrospective study on 137 consecutive patients with SARS-CoV-2 infection and available data on serum 25(OH)D levels, who were admitted to our Institution between March 1 and April 30, 2020. Patients were divided into two groups: survivors (n = 78; 57%) and non-survivors (n = 59; 43%). RESULTS: At admission, all patients showed hypovitaminosis D. Median total serum 25(OH)D levels at admission were significantly higher in survivors than non-survivors (12 ng/mL vs 8 ng/mL; p < 0.01). Non-survivors exhibited significantly higher median levels of white blood cell (WBC) count, neutrophil-to-lymphocyte count ratio (NLR), high-sensitivity C-reactive protein (hsCRP), ferritin, interleukin 6 (IL-6), D-dimer, fibrinogen, and procalcitonin (PCT) compared to survivors at three different time points during hospitalization. In a multivariate analysis performed by a logistic regression model, serum 25(OH)D levels were significantly inversely associated with risk of COVID-19-related in-hospital mortality (odds ratio, 0.91; 95% confidence interval, 0.85–0.98; p = 0.01). According to receiver operating characteristic curve analysis, hsCRP, NLR, ferritin, and D-dimer were the best predictive biomarkers for poor prognosis of COVID-19, whereas IL-6, PCT, fibrinogen, 25(OH)D, WBC count, and tumor necrosis factor alpha (TNF-α) may serve as supportive biomarkers for worse clinical course of the disease. CONCLUSIONS: We found a markedly high prevalence (100%) of hypovitaminosis D in patients admitted to hospital with COVID-19, suggesting a possible role of low vitamin D status in increasing the risk of SARS-CoV-2 infection and subsequent hospitalization. The inverse association between serum 25(OH)D levels and risk of in-hospital mortality observed in our cohort suggests that a lower vitamin D status upon admission may represent a modifiable and independent risk factor for poor prognosis in COVID-19. The coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first reported in the city of Wuhan (China) in December 2019 and was declared a global pandemic by the World Health Organization (WHO) on March 11, 2020 . Over recent months, COVID-19 has placed unprecedented strain on healthcare systems worldwide, posing serious threats to global health. Italy was the first Western country to be hit by the COVID-19 outbreak. Older age and underlying comorbid conditions (such as cardiovascular disease, hypertension, obesity, and diabetes mellitus) have emerged as major risk factors for mortality related to COVID-19 (1) (2) (3) (4) (5) (6) (7) . A dysregulated immune response resulting in the so-called "cytokine release syndrome" (also known as "cytokine storm") has been shown to play a critical role in the pathophysiology of the most severe cases of COVID-19 (8, 9) . Patients with severe manifestations of COVID-19 exhibit significantly increased circulating levels of C-reactive protein and several pro-inflammatory cytokines and chemokines, including tumor necrosis factor alpha (TNF-a), interleukin-1 beta (IL-1b), interleukin 6 (IL-6), and interferon-gamma (IFN-c) (10) (11) (12) . In turn, increased circulating levels of such cytokines result in a systemic hyperinflammatory state characterized by increased activity of CD8þ cytotoxic T cells, augmented differentiation of T helper (Th) 17 cells, and reduced activity of regulatory T cells (13) . These abnormal immune responses can lead to acute respiratory distress syndrome and multiorgan failure (13) (14) (15) ). An increase in neutrophil counts and a marked reduction of peripheral lymphocyte counts (primarily CD4þ and CD8þ T cells) have also been reported in severe cases of COVID-19, with the degree of lymphopenia correlating with disease severity (11, 16) . Severe cases of COVID-19 also exhibit abnormal coagulation results, which consist of significantly elevated concentrations of fibrinogen, D-dimer, and other fibrin degradation products, along with significantly longer prothrombin time and activated partial thromboplastin time (11, 16, 17) . These findings suggest the development of overt disseminated intravascular coagulation (17) . Over the last decade, several mechanistic studies have shown that vitamin D exerts anti-inflammatory and immunomodulatory properties, beyond its well-established role in the regulation of calcium and bone homeostasis (18) . Vitamin D has been shown to play a pivotal role in the regulation of both innate and adaptive immune responses, promoting antiviral effector mechanisms, reducing the expression of several pro-inflammatory cytokines, and favoring tolerogenic responses (18) (19) (20) (21) (22) . Preclinical evidence supports that calcitriol (the active metabolite of vitamin D, also referred to as 1,25-dihydroxyvitamin D3) exerts various effects on both innate and adaptive immune systems, resulting in induction of anti-inflammatory pathways and immune tolerance. With regard to innate immunity, calcitriol is able to induce the transcription of antimicrobial peptides (e.g., cathelicidin and defensin b2) in several human cell lines (keratinocytes, myeloid cells, monocytes/ macrophages, and neutrophils) (23) (24) (25) (26) . Also, calcitriol (a) promotes the differentiation of monocytes/macrophages and enhances their chemotactic and phagocytic capacity (27, 28) ; (b) inhibits the synthesis of pro-inflammatory cytokines (including IL-6 and TNF-a) by monocytes and macrophages (29) ; (c) reduces macrophage surface expression of major histocompatibility complex-class II molecules, thus decreasing the macrophage antigen presentation and T cell stimulatory ability (30, 31) ; (d) promotes the shift of macrophage polarization from M1 phenotype (pro-inflammatory or "classically activated" macrophages) toward M2 phenotype (anti-inflammatory or "alternatively activated" macrophages) (32) ; and (e) modulates the differentiation and function of dendritic cells, rendering them more tolerogenic and reducing their antigen-presenting capacity (33) (34) (35) (36) (37) . With regard to adaptive immunity, calcitriol up-regulates regulatory T cells (38) and promotes the shift of T cells from an "effector" toward a "regulatory" and anti-inflammatory phenotype by reducing Th1 and Th17 cell differentiation and favoring Th2 cell differentiation (39) (40) (41) . Additionally, immune cells are both vitamin D targets and local producers of vitamin D (42) . Indeed, functional vitamin D receptor has been identified in almost all immune cells, including neutrophils, T cells, and antigen-presenting cells (macrophages and dendritic cells) (43) (44) (45) (46) as well as in human airway epithelial cells (47) . In addition, several immune cells (e.g., macrophages, dendritic cells, and T-and B-lymphocytes) have been found to express the vitamin D-activating enzymes 25-and 1a-hydroxylase (30, (48) (49) (50) (51) . In vitro studies also suggest that vitamin D plays an important role in local "respiratory homeostasis" either by inducing the expression of antimicrobial peptides or by directly affecting the replication of respiratory viruses (47) . Vitamin D deficiency represents a global pandemic afflicting more than 1 billion individuals across all age groups worldwide (52) . Moreover, there is an overlap between risk factors for vitamin D deficiency and severe COVID-19 (such as Black or Asian ethnic origin, older age, and obesity) (53) . Hence, over the last few months, several researchers have suggested vitamin D deficiency as an independent risk factor for COVID-19 infection and adverse outcomes in the context of established disease (53) (54) (55) (56) . Similarly, there has been a growing interest in a potential role for vitamin D as an adjuvant immunomodulatory agent able to prevent SARS-CoV-2 infection or counteract the development of the cytokine release syndrome and improve outcomes in the setting of COVID-19 (53, 54, 57) . Therefore, we conducted a retrospective cohort study among patients with COVID-19 admitted to our Institution during the Italian COVID-19 outbreak, comparing the levels of inflammatory markers at admission between survivors and nonsurvivors with confirmed SARS-CoV-2 infection. Our study primarily aimed to measure serum 25-hydroxyvitamin D [25(OH)D] levels upon inpatient admission, in order to evaluate the relationship between vitamin D status and different markers of inflammation, coagulation, and sepsis in this population. We conducted a retrospective study including patients with confirmed SARS-CoV-2 infection who were consecutively admitted to our Institution (Tor Vergata University Hospital-PTV, Rome, Italy) between March 1 and April 30, 2020. Serum 25(OH)D levels were measured at admission in all patients, as per our institutional protocol. The study was reviewed and approved by the Ethics Committee of University of Rome Tor Vergata (Registration Number: 141/ 20, July 23, 2020). At admission, all patients provided written informed consent to anonymous data collection and analysis for research purposes. Patients were divided into two groups according to the outcomes of survival and death, namely survivors and non-survivors. For each group, we considered the length of hospital stay. In addition, the length of stay in the intensive care unit (ICU) was considered as a clinical marker of disease severity (need for invasive mechanical ventilation). The medical records of patients were independently reviewed by two members of our research team. Epidemiological, clinical, radiological, and laboratory data were collected through electronic medical records (ModulabV R ) and recorded in an anonymous inpatient COVID-19 database. Initial diagnosis of COVID-19 was made by an infectious disease specialist based on clinical symptoms (cough, fever, dyspnea, and/or anosmia) and imaging tests (chest X-ray and/or computed tomography) indicative of acute respiratory tract infection and COVID-19 pneumonia. Laboratory confirmation of SARS-CoV-2 infection was made from nasopharyngeal swab samples obtained upon hospital admission and analyzed through real-time reverse transcription polymerase chain reaction (RT-PCR) for 2019-nCoV RNA extraction according to the manufacturer's instructions (RT-PCR kit Seegene AllplexTM 2019-nCoV Assay, Seegene, Seoul, South Korea). Hematology and biochemical parameters were measured on blood, serum, and plasma samples collected upon admission to the emergency department, infectious disease unit, or ICU. White blood cell (WBC) count was determined by using automated hematological analyzer (Dasit-Sysmex, Milan, Italy). We also determined the neutrophil-to-lymphocyte ratio (NLR) as a marker of systemic inflammation (58) (normal values range between 0.78 and 3.53) (59) . Serum levels of high-sensitivity C-reactive protein (hsCRP; reference range 0-5 mg/L) were measured by using an immunoturbidimetric method (Abbott Diagnostics, Milan, Italy). Serum levels of IL-6 (reference range: 0-50 pg/ mL) were measured using chemiluminescence method (IMMULITE 2000 instrument, Siemens, Milan, Italy). Serum levels of TNF-a (reference range: 0-12.4 pg/mL) were measured using the enzyme-linked immunosorbent assay (ELISA) technique (DRG, International Instruments GmbH, Marburg, Germany). Serum levels of ferritin (reference range: 21.81-274.66 ng/mL), procalcitonin (PCT; reference range: 0.01-0.50 ng/mL) were measured using the chemiluminescence method (Architect Instrument, Abbott, Milan, Italy). Total serum 25(OH)D was measured by electrochemiluminescence (Abbott Architect Instrument, Milan, Italy), with the limit of quantitative value at 2.2 ng/mL at 20% coefficient variation. Plasma fibrinogen concentrations (reference range: 200-400 mg/dL) were measured using the Clauss method (ACL-TOP instrumentation, Werfen, Milan, Italy). Plasma D-dimer levels (reference range: 0-500 ng/mL) were measured by ACL-TOP instrumentation (Werfen, Milan, Italy). Hematological and biochemical parameters were measured at three different time points: (1) time of hospital admission (T1), (2) midpoint of hospitalization (T2), and (3) 1 day before discharge or death (T3) for survivors and non-survivors, respectively. Total serum 25(OH)D levels were only measured at admission (T1). Descriptive statistics such as frequency, percentage, mean and standard deviation (SD), median, and percentiles were calculated. Both the histogram and the Kolmogorov-Smirnov test of normality (p value < 0.05) were used to check whether the data were normally distributed. In the presence of a normal distribution of data, parametric tests were used such as analysis of variance with Bonferroni post hoc test in the case of more than two variables, or t test in the case of two variables. Non-parametric tests, such as Kruskal-Wallis test (variables with more than two categories) and Mann-Whitney U test (variables with two categories), were used to test differences between different groups. A p value of less than 0.05 was considered statistically significant in all statistical analyses. The comparison between the percentages was performed using the v 2 test. Multivariate analysis performed by a logistic regression model was used to determine the independent association of serum 25(OH)D levels (expressed as continuous variable) and risk of COVID-19-related in-hospital mortality. We initially included in the logistic regression model the following covariates: (1) continuous variables: 25(OH)D, age, body mass index (BMI), WBC count, NLR, hsCRP, fibrinogen, Ddimer, IL-6, TNF-a, ferritin, and PCT and (2) categorical variables: sex, hypertension and cardiovascular disease (CVD), diabetes mellitus, obesity (expressed as a BMI value of !30 kg/m 2 ), and malignancy (active malignancy or history of previous malignancy). After backward elimination process, WBC count, hsCRP, IL-6, TNF-a, ferritin, PCT, BMI (continuous variable), hypertension and CVD, diabetes mellitus, and malignancy were excluded from the model because of the p value > 0.1. p values < 0.05 were considered statistically significant. All statistical analyses were performed using MedCalc Version 18.2.18 (MedCalc Software Ltd, Ostend, Belgium). Receiver operating characteristic (ROC) curve analysis was performed to assess sensitivity and specificity of different potential predictive biomarkers for worse clinical course and adverse outcomes of COVID-19. The area under the curve (AUC) was used to assess the prognostic accuracy. AUC values ranging from 0.5 to 1.0 indicate the prognostic marker discriminatory ability for COVID-19 morbidity and mortality. A total of 137 consecutive patients admitted to our Institution (between March 1 and April 30, 2020) were enrolled into this retrospective single-center study. All patients had laboratoryconfirmed cases of COVID-19, were Caucasian, and resided in the Lazio region. Participant demographics, admission characteristics, and length of stay are shown in Table 1 . Table 2 lists the prevalence of major comorbidities in our study population. Participants included 89 males (65%) and 48 females (35%). Patients were divided into two groups according to the outcomes of survival and death, namely survivors (n ¼ 78; 57%) and nonsurvivors (n ¼ 59; 43%). Mean age in the survivor group was significantly lower than that in the non-survivor group (65 ± 13 vs 70 ± 14 years, respectively; p ¼ 0.01). The survivor group included 43 males (55%) and 35 females (45%), while the non-survivor group included 46 males (78%) and 13 females (22%). The percentage of males was significantly higher in non-survivors compared to survivors (78% vs 55%, p < 0.005), whereas the percentage of females was significantly lower in non-survivors compared to survivors (22% vs 45%, p < 0.005). Median values of BMI between survivors and non-survivors were comparable (Table 1) . However, the percentage of obese patients (defined as a BMI value of !30 kg/m 2 ) was significantly higher in non-survivors compared to survivors (29% vs 15%; n ¼ 17 vs 12; p ¼ 0.007) ( Table 2 ). Among survivors, there were 27 patients (35%) with hypertension and CVD, 8 patients (10%) with diabetes mellitus, and 7 patients (9%) with active or previous malignancy. Among non-survivors, there were 25 patients (42%) with hypertension and CVD, 6 patients (10%) with diabetes mellitus, and 9 patients (15%) with active or previous malignancy. There was no statistically significant difference in the percentage of hypertension and CVD, diabetes mellitus, and active or previous malignancy between survivors and non-survivors ( Table 2) . All patients received the same standard care for the treatment of COVID-19 (as per our institutional protocol) consisting of a combination therapy with dexamethasone plus hydroxychloroquine and lopinavir/ritonavir administered shortly after the admission. None of the patients reported vitamin D supplementation prior to hospital admission. The mean length of hospital stay was 30 ± 18 days and 15 ± 10 days for survivors and non-survivors, respectively (p < 0.001). With regard to length of stay in the ICU, survivors and non-survivors spent a mean time in the ICU of 3 ± 7 days and 8 ± 8 days, respectively (p < 0.001) ( Table 1) . Hematological and biochemical parameters and pro-inflammatory cytokines were measured at three different time Length of in-hospital stay (days) Length of stay in ICU (days) points: (1) time of hospital admission (T1), (2) midpoint of hospitalization (T2), and (3) 1 day before discharge or death (T3) for survivors and non-survivors, respectively. Total serum 25(OH)D levels were only measured at admission (T1). Baseline values at the time of hospital admission (T1) and changes in hematological and biochemical parameters and pro-inflammatory cytokines at different time points during hospitalization (T2 and T3) are shown in Figure 1 (expressed as median values) and Figure 2 (expressed as median values, with the addition of SD and interquartile ranges) and listed in Supplementary Table S1 . Values of hematological and biochemical parameters in survivors (SU) and non-survivors (NSU) were compared at each time point, namely: T1-SU vs T1-NSU; T2-SU vs T2-NSU; and T3-SU vs T3-NSU. Kruskal-Wallis non-parametric test was used to test differences between groups at different time points. In survivors, median WBC count tended to increase from T1 (7 Â 10 3 /mL) to T2 (8 Â 10 3 /mL) and normalized at T3 (6 Â 10 3 /mL). In non-survivors, median WBC count was higher than survivors at T1 (8 Â 10 3 /mL), increased at T2 (13 Â 10 3 /mL), and remained above the normal range at T3 (15.5 Â 10 3 /mL). In survivors, median NLR was 4.1 at T1, 7.5 at T2, and 2.5 at T3. In non-survivors, median NLR was markedly high at T1 (NLR: 12) and sharply increased at T2 (NLR: 21) and T3 (NLR: 30). The difference in median WBC count and NLR between survivors and non-survivors was statistically significant at all time points (T1-SU vs T1-NSU, T2-SU vs T2-NSU, T3-SU vs T3-NSU; p < 0.05). Median hsCRP values were above the normal range at all time points in survivors: T1 (41 mg/L), T2 (24 mg/L), and T3 (3 mg/L). However, median hsCRP values were markedly higher in non-survivors at all time points, from T1 (128 mg/ L) to T2 (122 mg/L) and T3 (122 mg/L). The difference in hsCRP values between survivors and non-survivors was statistically significant at all time points (T1-SU vs T1-NSU, T2-SU vs T2-NSU, T3-SU vs T3-NSU; p < 0.05). In survivors, median fibrinogen levels were above the normal range at T1 (560 mg/dL), whereas they decreased and normalized at T2 (370 mg/dL) and T3 (310 mg/dL). In non-survivors, median fibrinogen levels were increased at all time points: T1 (660 mg/dL), T2 (600 mg/dL), and T3 (610 mg/dL). The difference in fibrinogen levels between survivors and non-survivors was statistically significant at all time points (T1-SU vs T1-NSU, T2-SU vs T2-NSU, T3-SU vs T3-NSU; p < 0.05). In survivors, median D-dimer levels were above the normal range at all time points, although they tended to decrease at T3: T1 (890 ng/mL), T2 (936 ng/mL), and T3 (699 ng/mL). In non-survivors, median D-dimer levels were abnormally high at all time points: T1 (1728.5 ng/mL), T2 (2164 ng/mL), and T3 (1945.5 ng/mL). The difference in D-dimer levels between survivors and non-survivors was statistically significant at all time points (T1-SU vs T1-NSU, T2-SU vs T2-NSU, T3-SU vs T3-NSU; p < 0.05). In survivors, median IL-6 levels were normal at T1 (29 pg/ mL), whereas they increased at T2 (81.5 pg/mL) and normalized at T3 (10 pg/mL). On the contrary, in non-survivors, median IL-6 levels were slightly increased at T1 (64 pg/mL), whereas they dramatically increased at T2 (202 pg/mL) and T3 (430 pg/mL). The difference in IL-6 levels between survivors and non-survivors was statistically significant at all time points (T1-SU vs T1-NSU, T2-SU vs T2-NSU, T3-SU vs T3-NSU; p < 0.05). In both survivors and non-survivors, TNF-a levels showed a similar Gaussian distribution. At T1, median TNF-a levels were slightly increased in both survivors and non-survivors: T1-SU (14.5 pg/mL) and T1-NSU (21 pg/mL). At T2, median TNF-a levels further increased in both survivors and nonsurvivors: T2-SU (23 pg/mL) and T2-NSU (35 pg/mL). At T3, median TNF-a levels decreased and normalized in survivors (T3-SU: 12 pg/mL), whereas they decreased without normalizing in non-survivors (T3-NSU: 21 pg/mL). Although non-survivors exhibited higher TNF-a levels at all time points compared to survivors, statistical significance between survivors and non-survivors was only observed at T1 (T1-SU vs T1-NSU; p < 0.05) and T3 (T3-SU vs T3-NSU; p < 0.05). In survivors, median ferritin levels were above the normal range at all time points in both survivors and non-survivors. In survivors, median ferritin levels were as follows: T1 (550 ng/mL), T2 (568.5 ng/mL), and T3 (501 ng/mL). In non-survivors, median ferritin levels were markedly elevated at T1 (1234 ng/mL) and T2 (1158 ng/mL), while they sharply increased at T3 (3342.5 ng/mL). The difference in ferritin levels between survivors and non-survivors was statistically significant at all time points (T1-SU vs T1-NSU, T2-SU vs T2-NSU, T3-SU vs T3-NSU; p < 0.05). In survivors, median PCT levels remained within the normal range at all time points: T1 (0.13 ng/mL), T2 (0.04 ng/mL), and T3 (0.04 ng/mL). In non-survivors, median PCT values were normal at T1 (0.47 ng/mL), whereas they were mildly elevated at T2 (0.6 ng/mL) and T3 (1.1 ng/mL). The difference in PCT levels between survivors and non-survivors was Total serum 25(OH)D levels represent the most reliable biomarker of vitamin D status (60) . At admission, all 137 patients showed hypovitaminosis D, defined as serum 25(OH)D levels <30 ng/mL according to the Endocrine Society guidelines on evaluation, treatment and prevention of vitamin D deficiency (61) . In our cohort, median total serum 25(OH)D levels at admission were significantly higher in survivors than non-survivors (12 ng/mL [25th-75th percentile: 7-15] vs 8 ng/mL [25th-75th percentile: 5-14]; p < 0.01) (Figure 3 ). We also evaluated 25(OH)D level as continuous variable through a logistic regression analysis to determine the independent association of 25(OH)D and in-hospital mortality. In our final logistic regression model after adjusting for major confounders (age, sex, obesity, NLR, fibrinogen, D-dimer), there was a significant inverse association between serum 25(OH)D levels and risk of COVID-19-related in-hospital mortality (odds ratio [OR], 0.91; 95% confidence interval [CI] 0.85-0.98; p ¼ 0.01; Table S2 ). Both survivors (n ¼ 78) and non-survivors (n ¼ 59) were further stratified into different subgroups according to their serum 25(OH)D status at admission (T1), as follows: Group 1: survivors with severe vitamin D deficiency, defined as serum 25(OH)D levels <10 ng/mL (n ¼ 31; 40% of all survivors). Group 2: survivors with mild to moderate vitamin D deficiency, defined as serum 25(OH)D levels between 10 and 19.9 ng/mL (n ¼ 32; 41% of all survivors). Group 3: survivors with vitamin D insufficiency, defined as serum 25(OH)D levels between 20 and 29.9 ng/mL (n ¼ 15; 19% of all survivors). Group 4: non-survivors with severe vitamin D deficiency, defined as serum 25(OH)D levels <10 ng/mL (n ¼ 38; 64% of all non-survivors). Group 5: non-survivors with mild to moderate vitamin D deficiency, defined as serum 25(OH)D levels between 10 and 19.9 ng/mL (n ¼ 17; 29% of all non-survivors). Group 6: non-survivors with vitamin D insufficiency, defined as serum 25(OH)D levels between 20 and 29.9 ng/mL (n ¼ 4; 7% of all non-survivors) (Supplementary Table S3) . Figure 4) . Moreover, the percentage of patients was significantly higher in the survivor subgroup with mild to moderate vitamin D deficiency compared to the non-survivor subgroup with mild to moderate vitamin D deficiency (group 2 vs group 5: 41% vs 29%, respectively; p < 0.05), as well as in the survivor subgroup with vitamin D insufficiency compared to the nonsurvivor subgroup with vitamin D insufficiency (group 3 vs group 6: 19% vs 7%, respectively; p < 0.05). Conversely, there was a trend toward a lower percentage of patients in the survivor subgroup with severe vitamin D deficiency compared to the non-survivor subgroup with severe vitamin D deficiency (group 1 vs group 4: 40% vs 64%, respectively; p ¼ 0.06) (Supplementary Table S3, Figure 4) . Table S4 ). Among non-survivors, mean Table S4 ). Among survivors, length of hospital stay decreased within groups that showed higher 25(OH)D levels, although statistical significance between different survivor subgroups was not reached. Overall, length of hospital stay was greater in all survivor groups compared to non-survivor groups stratified by 25(OH)D status, although statistical significance was observed only between survivor and non-survivor subgroups with mild to moderate and severe vitamin D deficiency (group 1 vs group 4, p < 0.05; group 2 vs group 5, p < 0.05) (Supplementary Table S4 ). Among survivors, mean length of stay in the ICU was 4 ± 9 days in group 1 [25(OH)D < 10 ng/mL], 3 ± 7 days in group 2 [25(OH)D between 10 and 19.9 ng/mL], and 0.6 ± 2 days in group 3 [25(OH)D between 20 and 29.9 ng/ mL] (Supplementary Table S4 Table S4 ). The difference in mean length of stay in the ICU was statistically significant between all survivor and non-survivor groups stratified by 25(OH)D status (group 1 vs group 4, group 2 vs group 5, group 3 vs group 6; p < 0.05) (Supplementary Table S4 ). Values of hematological parameters, biochemical parameters, and pro-inflammatory cytokines measured at the time of hospital admission (T1) in survivor and non-survivor subgroups stratified by serum 25(OH)D status are presented in Supplementary Appendix 1 and Supplementary Table S5 . With regard to the analysis examining the relationship between 25(OH)D status and hematological and biochemical parameters and pro-inflammatory cytokines, we only took into account measurements performed at T1 (time of hospital admission) because this was the same time point at which data on serum 25(OH)D levels were available. ROC curve analysis and AUC were performed to identify reliable predictive biomarkers for worse clinical course and adverse outcomes (expressed as mortality) of COVID-19. The AUC provides an overall measure of accuracy of different analytes. Based on AUC values of different analytes at admission, we observed that hsCRP (AUC ¼ 0.719), NLR (AUC ¼ 0.718), ferritin (AUC ¼ 0.709), and D-dimer (AUC¼ 0.708) at admission were the best predictive biomarkers for adverse outcome of COVID-19 (death). On the other hand, IL-6 (AUC ¼ 0.665), PCT (AUC ¼ 0.647), fibrinogen (AUC ¼ 0.634), 25(OH)D (AUC ¼ 0.617), WBC count (AUC ¼ 0.614), and TNF-a (AUC ¼ 0.570) could be considered as supportive biomarkers for disease severity and mortality ( Figure 5 , Table 3 ). Our study suggests an involvement of low vitamin D status as a potential risk factor for SARS-CoV-2 infection and COVID-19-related hospitalization. First, we found a markedly high prevalence of hypovitaminosis D (100%), which was present in all 137 patients upon admission, in accordance with findings observed in another study (62) . These findings are also in line with those observed in a recent large, real-word, population-based study, which identified an independent and significant association between a low 25(OH)D status (<30 ng/mL) and the increased likelihood of COVID-19 infection (63) . In a univariate analysis from the aforementioned study, low 25(OH)D levels were also significantly associated with an increased likelihood of hospitalization due to COVID-19 (63) . A single-center retrospective cohort study further confirmed that likely deficient vitamin D status was associated with increased risk for COVID-19 (defined by a positive PCR test result) (64) . In addition, a smaller retrospective cohort study found significantly lower serum 25(OH)D levels in SARS-CoV-2 PCR-positive patients compared to negative patients (65) . Accordingly, a retrospective observational analysis conducted in the US among 191,779 patients with SARS-CoV-2 results performed from mid-March through mid-June 2020 showed that the SARS-CoV-2 positivity rate was higher in patients with deficient 25(OH)D values (<20 ng/mL) compared to patients with adequate values (30-34 ng/mL) and those with values !55 ng/mL (66) . In our cohort, median total serum 25(OH)D levels at admission were significantly higher in survivors than nonsurvivors (12 vs 8 ng/mL). In non-survivors, median total serum 25(OH)D levels at admission were therefore indicative of severe vitamin D deficiency (<10 ng/mL). In the multivariate analysis performed by a logistic regression model, serum 25(OH)D levels were significantly inversely associated with risk of COVID-19-related in-hospital mortality (OR, 0.91; 95% CI, 0.85-0.98; p ¼ 0.01; Table S2 ), independent of age, sex, BMI, markers of inflammation, coagulation and sepsis, and presence of major comorbidities (hypertension and CVD, diabetes mellitus, obesity, and active malignancy or history of previous malignancy). Length of stay in the ICU was considered as a clinical marker of disease severity due to the need for invasive mechanical ventilation. Non-survivors showed a significantly shorter length of in-hospital stay than survivors (mean: 15 vs 30 days), as well as a significantly longer length of stay in the ICU compared to survivors (mean: 8 vs 3 days) ( Table 1) . With regard to subgroups stratified by 25(OH)D status, we observed that in survivors the length of hospital stay tended to be shorter in patients who had higher levels of 25(OH)D (Supplementary Table S4 ), suggesting that a higher vitamin D level may favor a faster recovery. Accordingly, a recent prospective cohort study conducted in hospitalized patients with and without COVID-19 aged !65 years found a significantly higher incidence of noninvasive ventilation support and high dependency unit admission among patients with vitamin D deficiency in the COVID-19-positive group. Importantly, study participants were considered vitamin D-replete in presence of 25(OH)D levels >30 nmol/L (corresponding to >12 ng/mL) (67) . In keeping with our findings, another Italian retrospective observational study conducted among 42 hospitalized patients with acute respiratory failure due to COVID-19 found a high prevalence of hypovitaminosis D (81%), and a survival analysis showed that patients with severe vitamin D deficiency [defined as 25(OH)D levels <10 ng/mL] exhibited a significantly higher mortality risk after 10 days of hospitalization compared to those with 25(OH)D levels !10 ng/mL (50% vs 5%, p ¼ 0.019) (62) . Another retrospective study conducted among 185 hospitalized patients with COVID-19 found that vitamin D deficiency [defined as serum total 25(OH)D levels <12 ng/mL] at the time of admission was associated with higher risk of invasive mechanical ventilation and death (68) . De Smet et al. (69) recently documented that male patients with COVID-19 exhibited progressively lower 25(OH)D levels with advancing disease radiologic stage (assessed by chest computed tomography), and vitamin D deficiency (<20 ng/mL) at admission was significantly and independently associated with COVID-19 mortality. Thus, a lower vitamin D status may partly account for adverse clinical outcomes and mortality in hospitalized patients with COVID-19. In particular, our results-together with the findings from the aforementioned studies (62, 67, 68) -may suggest the existence of a baseline 25(OH)D threshold falling within a given range (probably between 8 and 12 ng/mL), which might predict poor clinical outcomes in hospitalized patients with COVID-19. However, large prospective studies are needed to confirm this hypothesis. In our cohort, non-survivors also exhibited significantly higher levels of markers of inflammation, coagulation, and sepsis (WBC count, NLR, hsCRP, ferritin, IL-6, D-dimer, fibrinogen, and PCT) compared to survivors at all time points (T1, T2, T3) during hospitalization. Non-survivors also showed significantly higher levels of the pro-inflammatory cytokine TNF-a compared to survivors at all time points (T1, T2, T3), although statistical significance was reached only at T1 and T3. Similar findings have recently been confirmed in other studies comparing asymptomatic patients or those with mild to moderate COVID-19 and patients with severe COVID-19 (70, 71) , as well as in studies comparing hospitalized COVID-19-negative and COVID-19positive patients (67) . A lower vitamin D status may have contributed, at least in part, to exacerbating systemic inflammation and cytokine release syndrome among non-survivors in our cohort. This explanation may be reasonable in light of the anti-inflammatory and immunomodulatory properties of vitamin D (18, 19, 21) . A recent Spanish retrospective case-control study (72) found that mean serum 25(OH)D levels were significantly lower in 216 hospitalized patients with COVID-19 compared to 197 sex-matched population-based controls (14 ± 7 vs 20.9 ± 7.4 ng/mL, respectively). Patients with COVID-19 also showed a higher prevalence of vitamin D deficiencydefined as serum 25(OH)D < 20 ng/mL-compared to controls (82% vs 47%). Furthermore, 25(OH)D levels significantly and inversely correlated with serum ferritin and Ddimer levels (72) . Similarly, a Chinese case-control study (73) found that median serum 25(OH)D levels were significantly lower in hospitalized patients with COVID-19 compared to healthy controls (55.6 vs 72 nmol/L; 22.2 vs 28.8 ng/ mL, respectively). In addition, authors found significantly higher rates of vitamin D deficiency-defined as serum 25(OH)D < 20 ng/mL-in COVID-19 cases (41.9%) compared to healthy controls (11.1%). Among COVID-19 cases, median serum 25(OH)D levels were significantly lower in severe/critical cases (38.2 nmol/L; 15.2 ng/mL) than in mild/ moderate cases (56.6 nmol/L; 22.6 ng/mL). Interestingly, ROC curve analysis identified a serum 25(OH)D level of 41.19 nmol/L (16.47 ng/mL) as a potential threshold to predict risk of SARS-CoV-2 infection and disease severity (73) . In our cohort, the percentage of obese patients was significantly higher among non-survivors compared to survivors, further supporting a critical role of obesity as a risk factor for poor prognosis of COVID-19, as it has recently been demonstrated by other groups (6, 7) . However, obese patients frequently exhibit hypovitaminosis D due to a number of possible reasons, including volumetric dilution of vitamin D in the large fat mass of such individuals (74), blunted catecholamine-induced release of vitamin D3 and 25(OH)D3 from adipose tissue (75) , as well as altered activity and expression of vitamin D-metabolizing enzymes in the liver and extrahepatic tissues (75) (76) (77) . Thus, hypovitaminosis D may be one of the drivers of poor prognosis in obese patients with COVID-19. Overall, findings from our study further suggest that lower vitamin D levels may partly contribute to worsening the clinical course and prognosis in hospitalized patients with COVID-19. Importantly, preliminary findings from intervention studies (including randomized placebo-controlled trials) conducted among hospitalized patients with COVID-19 suggest that vitamin D supplementation (in different doses and formulations) may promote viral clearance and reduce the severity of the disease (78) (79) (80) . Cross-sectional clinical studies have shown that lower vitamin D levels are significantly associated with acute respiratory tract infections (81) (82) (83) . A British cohort study showed that each 10-nmol/L (4-ng/mL) increase in serum 25(OH)D levels was associated with a 7% lower risk of respiratory infection, after adjustment for adiposity, lifestyle, and socioeconomic factors (84) . A meta-analysis by Martineau et al. (85) demonstrated that vitamin D supplementation protected against acute respiratory tract infections, and protective properties of vitamin D were stronger in participants with baseline serum 25(OH)D levels of <25 nmol/L (corresponding to <10 ng/ mL, indicative of severe vitamin D deficiency) compared to those with baseline levels of !25 nmol/L (!10 ng/mL). Brenner et al. (86) recently assessed the association between serum 25(OH)D levels and mortality from respiratory diseases over a 15-year follow-up period in a cohort of 9548 adults aged 50 to 75 years. Authors found that participants with vitamin D insufficiency and deficiency exhibited strongly increased respiratory mortality compared to those with sufficient vitamin D status (86) . Of note, a serum 25(OH)D level of approximately !40 ng/mL may provide protection against acute viral respiratory infections, as it has been showed in a prospective cohort study conducted in 198 healthy adults (87) . Therefore, attainment and maintenance of target serum 25(OH)D levels of 40 to 60 ng/mL may be critical in order to achieve the immunomodulatory properties of vitamin D in vivo and to effectively prevent acute respiratory tract infections, including SARS-CoV-2 infection (19, 54, 57) . According to our ROC curve analysis, hsCRP (AUC ¼ 0.719), NLR (AUC ¼ 0.718), ferritin (AUC ¼0.709), and D-dimer (AUC¼ 0.708) at admission were the best predictive biomarkers for adverse outcomes of COVID-19. NLR is a well-known marker of systemic inflammation (58) , which has previously been suggested as a useful prognostic marker for hospital mortality in patients with acute exacerbation of chronic obstructive pulmonary disease (88) . Importantly, NLR has recently been shown to effectively predict mortality in hospitalized patients with COVID-19 (89) . According to our ROC curve analysis, IL-6 (AUC ¼ 0.665), PCT (AUC ¼ 0.647), fibrinogen (AUC ¼ 0.634), 25(OH)D (AUC¼ 0.617), WBC count (AUC ¼0.614), and TNF-a (AUC ¼ 0.570) may be considered as supportive biomarkers for COVID-19 mortality. The measurement of these analytes may also serve as a panel of supportive biomarkers for worse clinical course of COVID-19. Similarly, other studies showed that TNF-a (90) and PCT (91) upon admission are reliable predictors of disease severity and death in hospitalized patients with COVID-19. Finally, in our cohort, mean age of the survivor group was significantly lower than that of the non-survivor group (65 ± 13 vs 70 ± 14 years, respectively). Older age has emerged as one of the main risk factors for severe COVID-19 since the start of the pandemic (92) . Intriguingly, older adults also represent individuals at high risk for vitamin D deficiency, as skin storage depots of 7-dehydrocholesterol (the precursor of vitamin D) and the human skin capacity to synthesize cholecalciferol upon sunlight exposure both decrease with age (93). In particular, age older than 70 years is an important risk factor for impaired vitamin D status (94) . Hence, aging-associated vitamin D deficiency may represent one of the drivers of COVID-19-related mortality in older adults by potentially triggering systemic inflammatory responses and endothelial dysfunction (95) . We acknowledge that our findings must be interpreted with caution due to a series of major limitations of the present study, including the retrospective database design, the small sample size, the lack of information on chronic lung disease, as well as the lack of a healthy control group. Moreover, as a single-center study, we cannot generalize our results to other settings. Conversely, the main strength of the study is the assessment of vitamin D status of participants upon admission. Our results could also be analyzed in future meta-analyses of observational studies assessing the vitamin D status in patients with COVID-19. More detailed data regarding the pharmacological treatments employed for the management of COVID-19 in our patients will be further assessed. Our study found a high prevalence of hypovitaminosis D in patients admitted to hospital with COVID-19, suggesting a possible role of low vitamin D status in increasing the risk of SARS-CoV-2 infection and subsequent hospitalization. We also found that non-survivors exhibited significantly lower vitamin D levels at admission compared to survivors, as well as higher levels of markers of inflammation, coagulation, and sepsis. Serum 25(OH)D levels were significantly inversely associated with the risk of COVID-19-related inhospital mortality, independent of age, sex, markers of inflammation, coagulation and sepsis, and major comorbidities. We therefore suggest that a lower vitamin D status upon admission may be a modifiable risk factor and early predictive marker for adverse outcomes and mortality in hospitalized patients with COVID-19. However, larger prospective studies should be conducted to further confirm the existence of a causal relationship between hypovitaminosis D and SARS-CoV-2 infection and COVID-19 severity. Moreover, large multicenter, randomized, double-blind controlled trials are needed to address whether vitamin D supplementation can effectively reduce the risk of SARS-CoV-2 infection, as well as COVID-19-related hospitalization, morbidity, and mortality. Arterial hypertension and the risk of severity and mortality of COVID-19 Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China: summary of a report of 72 314 cases from the Chinese Center for Disease Control and Prevention The Northwell COVID-19 Research Consortium, et al. Presenting characteristics, comorbidities, and outcomes among 5700 patients hospitalized with COVID-19 in the New York City Area Obesity a risk factor for severe COVID-19 infection: multiple potential mechanisms Coronavirus infections and type 2 diabetes-shared pathways with therapeutic implications How important is obesity as a risk factor for respiratory failure, intensive care admission and death in hospitalised COVID-19 patients Obesity as a predictor for a poor prognosis of COVID-19: a systematic review SARS-CoV-2 and viral sepsis: observations and hypotheses Cytokine storms: understanding COVID-19 Clinical features of patients infected with 2019 novel coronavirus in Wuhan Longitudinal characteristics of lymphocyte responses and cytokine profiles in the peripheral blood of SARS-CoV-2 infected patients Hydroxychloroquine in the COVID-19 pandemic era: in pursuit of a rational use for prophylaxis of SARS-CoV-2 infection Cytokine release syndrome in severe COVID-19 The pathogenesis and treatment of the 'Cytokine Storm' in COVID-19 COVID-19 and multiorgan response Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study Abnormal coagulation parameters are associated with poor prognosis in patients with novel coronavirus pneumonia Vitamin D: not just the bone. Evidence for beneficial pleiotropic extraskeletal effects Editorial -Vitamin D status: a key modulator of innate immunity and natural defense from acute viral respiratory infections The role of vitamin D and omega-3 PUFAs in islet transplantation Immunologic effects of vitamin D on human health and disease Influence of vitamin D on islet autoimmunity and beta-cell function in type 1 diabetes Cutting edge: 1,25-dihydroxyvitamin D3 is a direct inducer of antimicrobial peptide gene expression Human cathelicidin antimicrobial peptide (CAMP) gene is a direct target of the vitamin D receptor and is strongly up-regulated in myeloid cells by 1,25-dihydroxyvitamin D3 The vitamin D-antimicrobial peptide pathway and its role in protection against infection Toll-like receptor triggering of a vitamin D-mediated human antimicrobial response Vitamin D and its analogs as regulators of immune activation and antigen presentation 1,25-Dihydroxyvitamin D3 exerts opposing effects to IL-4 on MHC class-II antigen expression, accessory activity, and phagocytosis of human monocytes Vitamin D inhibits monocyte/macrophage proinflammatory cytokine production by targeting MAPK phosphatase-1 Identification and immune regulation of 25-hydroxyvitamin D-1-alpha-hydroxylase in murine macrophages 25-Dihydroxyvitamin D3 curtails the inflammatory and T cell stimulatory capacity of macrophages through an IL-10-dependent mechanism 1,25-Dihydroxyvitamin D 3 promotes high glucose-induced M1 macrophage switching to M2 via the VDR-PPARc signaling pathway Vitamin D3 affects differentiation, maturation, and function of human monocyte-derived dendritic cells 1 Alpha,25-dihydroxyvitamin D3 inhibits differentiation, maturation, activation, and survival of dendritic cells leading to impaired alloreactive T cell activation Suppressive effect of 1alpha,25-dihydroxyvitamin D3 on type I IFN-mediated monocyte differentiation into dendritic cells: impairment of functional activities and chemotaxis 1,25-Dihydroxyvitamin D3 alters murine dendritic cell behaviour in vitro and in vivo 1,25-Dihydroxyvitamin D3 restrains CD4þ T cell priming ability of CD11c þ dendritic cells by upregulating expression of CD31 1,25-Dihydroxyvitamin D3 and IL-2 combine to inhibit T cell production of inflammatory cytokines and promote development of regulatory T cells expressing CTLA-4 and FoxP3 25-dihydroxyvitamin D3 induces an autoantigen-specific T-helper 1/Thelper 2 immune shift in NOD mice immunized with GAD65 (p524-543) 25-Dihydroxyvitamin d3 has a direct effect on naive CD4(þ) T cells to enhance the development of Th2 cells Vitamin D action: lessons from VDR and Cyp27b1 null mice Vitamin D: modulator of the immune system Vitamin D and immune function 1,25-dihydroxyvitamin D3 receptors in human leukocytes Human neutrophils express messenger RNA of vitamin D receptor and respond to 1alpha,25-dihydroxyvitamin D3 Vitamin D metabolism and signaling in the immune system Vitamin D modulation of innate immune responses to respiratory viral infections Identification of 25-hydroxyvitamin D3 1alpha-hydroxylase gene expression in macrophages DCs metabolize sunlight-induced vitamin D3 to 'program' T cell attraction to the epidermal chemokine CCL27 Immune regulation of 25-hydroxyvitamin-D3-1alpha-hydroxylase in human monocytes Modulatory effects of 1,25-dihydroxyvitamin D3 on human B cell differentiation The vitamin D deficiency pandemic: approaches for diagnosis, treatment and prevention Vitamin D for COVID-19: a case to answer? Evidence that vitamin D supplementation could reduce risk of influenza and COVID-19 infections and deaths Tripartite combination of candidate pandemic mitigation agents: Vitamin D, quercetin, and estradiol manifest properties of medicinal agents for targeted mitigation of the COVID-19 pandemic defined by genomics-guided tracing of SARS-CoV-2 targets in human cells Evidence regarding vitamin D and risk of COVID-19 and its severity Perspective: improving vitamin D status in the management of COVID-19 Prognostic significance of neutrophil-to-lymphocyte ratio in patients with sepsis: a prospective observational study What is the normal value of the neutrophil-to-lymphocyte ratio? 25-Hydroxyvitamin D as a biomarker of vitamin D status and its modeling to inform strategies for prevention of vitamin D deficiency within the population Evaluation, treatment, and prevention of vitamin D deficiency: an Endocrine Society clinical practice guideline Vitamin D deficiency as a predictor of poor prognosis in patients with acute respiratory failure due to COVID-19 Frenkel-Morgenstern M. Low plasma 25(OH) vitamin D level is associated with increased risk of COVID-19 infection: an Israeli population-based study Association of vitamin D status and other clinical characteristics with COVID-19 test results Cant u M. 25-Hydroxyvitamin D concentrations are lower in patients with positive PCR for SARS-CoV-2 SARS-CoV-2 positivity rates associated with circulating 25-hydroxyvitamin D levels Vitamin D status and outcomes for hospitalised older patients with COVID-19 Vitamin d deficiency and outcome of COVID-19 patients Serum 25(OH)D level on hospital admission associated with COVID-19 stage and mortality Analysis of vitamin D level among asymptomatic and critically ill COVID-19 patients and its correlation with inflammatory markers 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 Vitamin D status in hospitalized patients with SARS-CoV-2 infection Does serum vitamin D level affect COVID-19 infection and its severity?-A case-control study Volumetric dilution, rather than sequestration best explains the low vitamin D status of obesity Impaired release of vitamin D in dysfunctional adipose tissue: new cues on vitamin D supplementation in obesity Obesity represses CYP2R1, the vitamin D 25-hydroxylase, in the liver and extrahepatic tissues Obesity decreases hepatic 25-hydroxylase activity causing low serum 25-hydroxyvitamin D Effect of calcifediol treatment and best available therapy versus best available therapy on intensive care unit admission and mortality among patients hospitalized for COVID-19: a pilot randomized clinical study Short term, high-dose vitamin D supplementation for COVID-19 disease: a randomised, placebo-controlled, study (SHADE study) Vitamin D supplementation associated to better survival in hospitalized frail elderly COVID-19 patients: the GERIA-COVID quasi-experimental study An association of serum vitamin D concentrations < 40 nmol/L with acute respiratory tract infection in young Finnish men Association between serum 25-hydroxyvitamin D level and upper respiratory tract infection in the Third National Health and Nutrition Examination Survey On the epidemiology of influenza Hypp€ onen E. Vitamin D status has a linear association with seasonal infections and lung function in British adults Vitamin D supplementation to prevent acute respiratory infections: individual participant data metaanalysis Vitamin D insufficiency and deficiency and mortality from respiratory diseases in a cohort of older adults: potential for limiting the death toll during and beyond the COVID-19 pandemic? Serum 25-hydroxyvitamin d and the incidence of acute viral respiratory tract infections in healthy adults Prognostic role of neutrophil-lymphocyte ratio and platelet-lymphocyte ratio for hospital mortality in patients with AECOPD Leucocyte subsets effectively predict the clinical outcome of patients with COVID-19 pneumonia: a retrospective case-control study An inflammatory cytokine signature predicts COVID-19 severity and survival Procalcitonin levels in COVID-19 patients SARS-CoV-2 and COVID-19 in older adults: what we may expect regarding pathogenesis, immune responses, and outcomes Aging decreases the capacity of human skin to produce vitamin D3 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 Vitamin D: a novel player in endothelial function and dysfunction We acknowledge all medical staff involved in the diagnosis and treatment of patients with COVID-19 at Tor Vergata University Hospital-PTV. MI and MM designed the research project, wrote the paper, supervised the project, and equally contributed to the manuscript. AB collected and retrieved clinical data, analyzed results, and contributed to the research project. MP performed statistical analysis. SL and MN performed, collected, and retrieved biochemical data. SB, AF, MIa, MA, and VC supervised the research project and reviewed the manuscript. All authors edited the manuscript. No honorarium, grant, or other forms of payment were received by authors to write this manuscript. No author has potential conflicts of interest to disclose. This study did not receive any funding for the conduct of the research and/or preparation of the article.