Vitamin D and socioeconomic deprivation mediate COVID-19 ethnic health disparities ================================================================================== * Leonardo Mariño-Ramírez * Maria Ahmad * Lavanya Rishishwar * Shashwat Deepali Nagar * Kara K. Lee * Emily T. Norris * I. King Jordan ## Abstract Ethnic minorities in developed countries suffer a disproportionately high burden of COVID-19 morbidity and mortality, and COVID-19 ethnic disparities have been attributed to social determinants of health. Vitamin D has been proposed as a modifiable risk factor that could mitigate COVID-19 health disparities. We investigated the relationship between vitamin D and COVID-19 susceptibility and severity using the UK Biobank, a large progressive cohort study of the United Kingdom population. Structural equation modelling was used to evaluate the ability of vitamin D, socioeconomic deprivation, and other known risk factors to mediate COVID-19 ethnic health disparities. Asian ethnicity is associated with higher COVID-19 susceptibility, compared to the majority White population, and Asian and Black ethnicity are both associated with higher COVID-19 severity. Socioeconomic deprivation mediates all three ethnic disparities and shows the highest overall signal of mediation for any COVID-19 risk factor. Vitamin supplements, including vitamin D, mediate the Asian disparity in COVID-19 susceptibility, and serum 25-hydroxyvitamin D (calcifediol) levels mediate Asian and Black COVID-19 severity disparities. Several measures of overall health also mediate COVID-19 ethnic disparities, underscoring the importance of comorbidities. Our results support ethnic minorities’ use of vitamin D as both a prophylactic and a supplemental therapeutic for COVID-19. ## Introduction The COVID-19 pandemic has had a disproportionate impact on ethnic minorities in developed countries1-4. In the United Kingdom (UK), Asian and Black minority groups have higher rates of COVID-19 diagnosis and mortality than the majority White population5. In the United States (US), American Indian, Black, and Latino populations have higher rates of COVID-19 cases, hospitalization, and death compared to White and Asian groups6. Ethnic disparities in COVID-19 have been attributed to a number of different factors, including a high burden of comorbidities, over-crowded living and working conditions, and reduced access to healthcare. Many of the risk factors associated with COVID-19 ethnic disparities are linked to structural barriers resulting from socioeconomic deprivation, segregation, and racial discrimination7. The vitamin D hypothesis holds that low levels of circulating vitamin D among ethnic minority patients may also contribute to COVID-19 disparities8. Vitamin D plays an important role in a number of immune related processes, including defense against viral infections9. Biosynthesis of vitamin D depends on exposure to sunlight and ultraviolet irradiation of 7-dehydrocholesterol, which is converted to vitamin D3 (cholecalciferol) in the skin and then 25-hydroxyvitamin D (calcifediol) in the liver10. Melanin absorbs ultraviolet radiation, thereby reducing the amount of irradiation available for vitamin D synthesis in the skin. Accordingly, darker skinned individuals have reduced vitamin D biosynthesis, and lower levels of circulating 25-hydroxyvitamin D, than lighter skinned individuals11. This is particularly true for ethnic minorities living in high latitude countries like the UK12. Studies on the role of vitamin D in immune defense against COVID-19 have yielded conflicting results. Retrospective cohort studies in the US found that vitamin D deficiency was associated with COVID-19 infection, with the strongest effect seen for Black patients13,14. A randomized control trial of hospitalized COVID-19 patients in Spain showed that treatment with a high dose of 25-hydroxyvitamin D significantly reduced admission to the intensive care unit15. A commentary in the British Medical Journal pointed to more than 40 such studies claiming that vitamin D mitigates COVID-1916. On the other hand, the UK National Institute for Health and Care Excellence surveyed the existing literature and concluded that there is insufficient evidence to support the use of vitamin D supplements for the treatment or prevention of COVID-1917. A pair of studies on the UK Biobank, a large progressive cohort study of the UK population, found that circulating levels of 25-hydroxyvitamin D are not significantly associated with the risk of COVID-19 disease or mortality18,19. The objective of this study was to investigate the relationship between vitamin D and ethnic health disparities using the UK Biobank. The previous studies of the UK Biobank cohort relied on multivariable modelling to evaluate the association of vitamin D with COVID-19 outcomes when all other variables, including patient ethnicity, were held constant18,19. The multivariable modelling approach does not directly address which risk factors may explain ethnic differences in COVID-19 outcomes. For this study, we used structural equation modelling (SEM) to evaluate the ability of vitamin D, socioeconomic deprivation, and other known risk factors to mediate COVID-19 ethnic health disparities. The SEM approach allows for the formulation and testing of causal hypotheses with observational data as well as the characterization of potential intervening mechanisms. Our study cohort includes an order of magnitude greater sample size than the previous UK Biobank studies allowing for greater statistical resolution on the question of COVID-19 disparities. ## Methods ### Study cohort and participant attributes The study cohort was obtained from the UK Biobank, a prospective cohort study intended to investigate the effects of genetic, environmental, and lifestyle factors on human health and disease20. The UK Biobank provides deep phenotypic, clinical, and genetic data for more than 500,000 participants enrolled between 2006 and 201021. We extracted the following information for UK Biobank participants: (1) age (Field 21022 Age at recruitment), (2) sex (Field 31: Sex), (3) ethnicity (Field 21000: ethnic background), (4) body mass index (Field 21001: Body mass index (BMI)), (5) smoking status (Field 20116: Smoking status), (6) diastolic blood pressure (Field 4079: Diastolic blood pressure, automated reading), (7) systolic blood pressure (Field 4080: Systolic blood pressure, automated reading), (8) Townsend deprivation index (Field 189: Townsend deprivation index at recruitment), (9) overall health rating (Field 2178: Overall health rating), (10) long-standing illness (Field 2188: Long-standing illness, disability, or infirmity), (11) supplements (Field 6155: Vitamin and mineral supplements), (12) vitamin D (Field 30890: Vitamin D), and (13) skin color (Field 1717: Skin colour). UK biobank participants’ serum 25-hydroxyvitamin D (calcifediol) levels were measured in nmol/L units via chemiluminescent immunoassay analysis (CLIA) using the DiaSorin Liaison XL platform as previously reported (see Field 30890: Vitamin D description)21,22. Participants’ ages ranged from 40 to 69 at enrollment, and age at enrollment was split into decades for analysis. The categories for self-reported skin color (SRSC), which are Very fair, Fair, Light olive, Dark olive, Brown, and Black, were collapsed into four categories for analysis: Fair, Olive, Brown, and Black. Participant self-reported supplement use was considered for vitamins A, B, C, D, E, V9, and multivitamins. The effects of supplement use were modelled individually and jointly. Socioeconomic deprivation (SED) was measured by the Townsend deprivation index (TDI), a widely used measure of SED, which is known to be associated with poor health outcomes23. The TDI incorporates four factors into its score: (1) unemployment (as a percentage of those aged 16 and economically active), (2) non-car ownership (as a percentage of all households, (3) non-home ownership (as a percentage of all households), and (4) household overcrowding24. UK Biobank TDI values range from -6 to 11, where negative values in indicate less deprivation and higher values indicate more deprivation, and the raw values were split into deciles for analysis using the thresholds provided by the UK Biobank. ### Participant Ethnicity UK Biobank participants self-identified as belonging to one of six ethnic groups at enrollment – Asian, Black, Chinese, Mixed, White, or Other – each of which is made up of a specific set of ethnic backgrounds. The UK Biobank ethnic group designations are based on the UK National Healthcare Service (NHS) categories, which were adopted from the UK Census. Since the UK Census makes no distinction between race and ethnicity, the UK Biobank ethnic groups used here may correspond to racial groups used in the US. The three largest ethnic groups in the UK Biobank – Asian, Black, and White – were used for analysis here. The Asian ethnic group consists of participants with Indian, Pakistani, Bangladeshi, or Any other Asian ethnic background. The Asian group does not include Chinese ethnicity, which has its own group. The Black ethnic group consists of participants with Caribbean, African, and Any other Black ethnic background. The White ethnic group consists of participants with British, Irish, and Any other White ethnic background. ### COVID-19 phenotypes and prevalence Data on COVID-19 outcomes were taken from the UK Biobank covid19_result table (accessed December 20, 2020). COVID-19 susceptibility is defined as whether or not a participant tested positive for SARS-CoV-2, where cases are participants with positive test results and controls are participants with negative test results. SARS-CoV-2 PCR tests were performed on participant nose and throat swabs as previously described25. COVID-19 severity is defined as whether or not a participant with a SARS-CoV-2 positive test result was hospitalized, where cases are hospitalized participants with positive test results and controls are non-hospitalized participants with positive test results. Unadjusted prevalence values for COVID-19 susceptibility and severity were calculated as the number of cases divided by the number of cases and controls. ### Statistical analyses and visualization All statistical analyses were performed using the R statistical language v4.0.426. Bar graphs and heat maps were generated using the R ‘ggplot’ v3.3.3 library. The effects of ethnicity on COVID-19 susceptibility and severity (cases and controls) were modeled with age and sex as covariates using multivariable logistic regression with the ‘glm’ function in R: COVID-19 ∼ Ethnicity + Age + Sex. Forest plots showing the logistic regression coefficient point estimates (β-values), 95% confidence intervals, and P-values were generated using the ‘viz-forest’ function from the ‘meta-viz’ v0.3.1 package in R27. Structural equation modelling (SEM) and mediation analyses were performed using the ‘sem’ function from the R statistical package ‘Lavaan’28. SEM was used to evaluate the ability of vitamin D, socioeconomic deprivation, and other known COVID-19 risk factors to mediate ethnic disparities in COVID-19 susceptibility and severity. Barron and Kenny’s stepwise method for mediation was used to evaluate the mediation effect of each risk factor independently, and then all risk factors that showed significant effects when considered independently were carried forward to a final multivariable structural equation model29. The ‘sem’ function was used to formulate each structural equation model and to calculate model coefficient values (β-values) and their P-values. For the final multivariable structural equation models, the mediation effects and their P-values were calculated using the ‘sem’ function. Mediation effects (indirect, direct, and total) are expressed as effect sizes (β-values), P-values, and percent mediation values. Structural equation path diagrams were generated using the web-based diagram tool ‘Lucidchart’. Additional details on the SEM approach used here are provided in the Supplementary Material. ## Results ### COVID-19 ethnic disparities We characterized the prevalence of COVID-19 susceptibility and severity for the three largest ethnic groups in the UK Biobank: Asian, Black, and White. Demographic and ethnic characteristics for the study cohort are summarized in Table 1. Numbers are shown for the total cohort and for SARS-CoV-2 test positive cases and hospitalized cases, which were used to define the COVID-19 susceptibility and severity phenotypes. 35,339 study participants were tested for SARS-CoV-2, resulting in 6,049 positive cases, 1,630 of whom were hospitalized. View this table: [Table 1.](http://medrxiv.org/content/early/2021/09/22/2021.09.20.21263865/T1) Table 1. Study cohort characteristics. Numbers and the percent share for each cohort category – total cohort, SARS-CoV-2 positive cases, and COVID-19 hospitalized cases – are shown for age, ethnicity, and sex. The prevalence of COVID-19 susceptibility and severity were compared across ethnic groups, stratified by sex and age (Figure 1A). Males showed higher overall prevalence for COVID-19 susceptibility and severity compared to females. The prevalence of COVID-19 susceptibility decreases with age, whereas severity increases with age. The Asian ethnic group shows the highest prevalence for COVID-19 susceptibility, followed by the Black and White groups. The Black ethnic group shows the highest prevalence for COVID-19 severity, followed by the Asian and White groups. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/22/2021.09.20.21263865/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2021/09/22/2021.09.20.21263865/F1) Figure 1. Ethnic disparities in COVID-19 susceptibility and severity. (A) Ethnic group prevalence values for COVID-19 susceptibility and severity: Asian (red), Black (blue), and White (orange). Prevalence values are stratified by sex and age brackets as shown. (B) Forest plots showing regression coefficient effect size estimates (β-values), 95% confidence intervals, and P-values for COVID-19 susceptibility and severity multivariable logistic regression models. Multivariable logistic regression was used to model the joint effects of age, sex, and ethnicity on COVID-19 outcomes (Figure 1B). Age, sex, and ethnicity are all significantly associated with COVID-19 susceptibility and severity. Age shows the strongest effects on both susceptibility and severity, followed by sex and ethnicity. Asian ethnicity is significantly associated with higher susceptibility compared to the White reference group, whereas the Black and White groups do not show a significant difference in susceptibility. The Black and Asian groups both show significantly higher severity than the White reference group. ### Skin color, vitamin D, and socioeconomic deprivation We explored the plausibility of the vitamin D hypothesis for COVID-19 disparities by looking at differences skin color and serum vitamin D levels among study participants from the Asian, Black, and White ethnic groups. Study participants from the Black ethnic group reported mainly black and brown skin color, and Asian participants reported primarily brown skin color, followed by olive and fair (Figure 2A). Participants from the White group reported mostly fair skin color, followed by olive. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/22/2021.09.20.21263865/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2021/09/22/2021.09.20.21263865/F2) Figure 2. Ethnic differences in skin color, vitamin D levels, and socioeconomic deprivation. (A) Heatmap showing relative frequency (density) of participant self-reported skin color for Asian, Black, and White ethnic groups. (B) Mean ±95% confidence interval vitamin D levels for Asian (red), Black (blue), and White (orange) ethnic groups. (C) Mean ±95% confidence interval Townsend deprivation index (TDI) deciles for ethnic groups. Participants from the White ethnic group showed the highest mean serum vitamin D levels, followed by the Black and Asian groups (Figure 2B). Differences in vitamin D levels between the White and other ethnic groups were much higher than the difference seen between the Black and Asian groups. We also compared socioeconomic deprivation (SED) among ethnic groups using the composite Townsend deprivation index (TDI). The Black ethnic group shows the highest mean TDI decile followed by the Asian and White groups (Figure 2C). ### Testing the vitamin D hypothesis We used structural equation modelling (SEM) to test the vitamin D hypothesis by evaluating the ability of vitamin D, and other known COVID-19 risk factors, to mediate the observed disparities in COVID-19 susceptibility and severity. In addition to vitamin D, we considered skin color, body mass index (BMI), systolic and diastolic blood pressure, long-standing illness (LSI), overall health ratio (OHR), smoking status, and supplement use as potential mediators of COVID-19 ethnic disparities. We were particularly interested in studying the role of socioeconomic deprivation, measured by the TDI, as a potential mediator of COVID-19 disparities. Age and sex were treated as exogenous variables with a direct effect on COVID-19 susceptibility and severity. We evaluated the mediation potential for each risk factor independently, and all risk factors that showed significant effects when considered alone were carried forward to a final multivariable structural equation model. This was done for each of the three significant COVID-19 ethnic disparities observed here: Susceptibility (Asian) and Severity (Asian and Black) (Figure 3). Each step in the SEM procedure is illustrated in Supplementary Figure 1 and the results for each step are shown in Supplementary Figures 2-5. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/22/2021.09.20.21263865/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2021/09/22/2021.09.20.21263865/F3) Figure 3. Mediation of COVID-19 ethnic disparities. Structural equation models are shown with ethnicity predictor variables (purple boxes), mediator variables (red boxes), exogenous predictor variables (gray boxes), and outcome variables (blue boxes). Arrows indicate the direction of predictor-to-outcome relationships. Direct effects of ethnicity on COVID-19 outcomes are shown along with the indirect (mediated) effects of ethnicity. Effect size estimates (β-values) are shown for each predictor-to-outcome relationship, and effect size statistical significance levels are indicated as *0.01 (2020). 26. Team, R. C. R: A language and environment for statistical computing., (R Foundation for Statistical Computing, 2013). 27. Kossmeier, M., Tran, U. S. & Voracek, M. Visualizing Meta-Analytic Data with R Package Metaviz. (R Foundation for Statistical Computing, 2019). 28. Rosseel, Y. Lavaan: An R package for structural equation modeling and more. Journal of Statistical Software 48, 1–36 (2012). 29. Baron, R. M. & Kenny, D. A. The moderator-mediator variable distinction in social psychological research: conceptual, strategic, and statistical considerations. J Pers Soc Psychol 51, 1173–1182, doi:10.1037//0022-3514.51.6.1173 (1986). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1037/0022-3514.51.6.1173&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3806354&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F22%2F2021.09.20.21263865.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1986F285400010&link_type=ISI)