Next Article in Journal
Complexities of Prostate Cancer
Next Article in Special Issue
Vitamin D and Beta Cells in Type 1 Diabetes: A Systematic Review
Previous Article in Journal
Learning from TCR Signaling and Immunological Synapse Assembly to Build New Chimeric Antigen Receptors (CARs)
Previous Article in Special Issue
The Effects of Vitamin D Supplementation in COVID-19 Patients: A Systematic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Vitamin D Deficiency and COVID-19: A Biological Database Study on Pathways and Gene-Disease Associations

by
Ángela Alcalá-Santiago
1,2,3,
Miguel Rodríguez-Barranco
2,4,5,*,
Marta Rava
6,7,
María Ángeles Jiménez-Sousa
7,8,
Ángel Gil
2,3,9,10,
María José Sánchez
2,4,5,11,* and
Esther Molina-Montes
1,2,3,5
1
Department of Nutrition and Food Science, Faculty of Pharmacy, University of Granada, 18071 Granada, Spain
2
Instituto de Investigación Biosanitaria ibs.GRANADA, 18012 Granada, Spain
3
Institute of Nutrition and Food Technology (INYTA) ‘José Mataix’, Biomedical Research Centre, University of Granada, Avenida del Conocimiento s/n, 18071 Granada, Spain
4
Andalusian School of Public Health, Cuesta del Observatorio 4, 18012 Granada, Spain
5
CIBER de Epidemiología y Salud Pública (CIBERESP), 28029 Madrid, Spain
6
National Center of Epidemiology (CNE), Institute of Health Carlos III (ISCIII), 28029 Madrid, Spain
7
CIBER de Enfermedades Infecciosas (CIBERINFEC), 28029 Madrid, Spain
8
Unit of Viral Infection and Immunity, National Center for Microbiology (CNM), Institute of Health Carlos III (ISCIII), 28029 Madrid, Spain
9
Department of Biochemistry and Molecular Biology II, Faculty of Pharmacy, University of Granada, 18071 Granada, Spain
10
CIBER de Obesidad y Nutrición (CIBEROBN), 28029 Madrid, Spain
11
Department of Preventive Medicine and Public Health, Faculty of Medicine, University of Granada, 18011 Granada, Spain
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(22), 14256; https://doi.org/10.3390/ijms232214256
Submission received: 4 November 2022 / Revised: 11 November 2022 / Accepted: 14 November 2022 / Published: 17 November 2022
(This article belongs to the Special Issue The Role of Vitamin D in Human Health and Diseases 2.0)

Abstract

:
Vitamin D (VD) is a fat-soluble vitamin, and pivotal for maintaining health. Several genetic markers have been related to a deficient VD status; these markers could confer an increased risk to develop osteoporosis and other chronic diseases. A VD deficiency could also be a determinant of a severe COVID-19 disease. This study aimed to interrogate genetic/biological databases on the biological implications of a VD deficiency and its association with diseases, to further explore its link with COVID-19. The genetic variants of both a VD deficiency and COVID-19 were identified in the genome-wide association studies (GWAS) catalog and other sources. We conducted enrichment analyses (considering corrected p-values < 0.05 as statistically significant) of the pathways, and gene-disease associations using tools, such as FUMA, REVIGO, DAVID and DisGeNET, and databases, such as the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO). There were 26 and 46 genes associated with a VD deficiency and COVID-19, respectively. However, there were no genes shared between the two. Genes related to a VD deficiency were involved in the metabolism of carbohydrates, retinol, drugs and xenobiotics, and were associated with the metabolic syndrome and related factors (obesity, hypertension and diabetes mellitus), as well as with neoplasms. There were few enriched pathways and disease connections for the COVID-19-related genes, among which some of the aforementioned comorbidities were also present. In conclusion, genetic factors that influence the VD levels in the body are most prominently associated with nutritional and metabolic diseases. A VD deficiency in high-risk populations could be therefore relevant in a severe COVID-19, underlining the need to examine whether a VD supplementation could reduce the severity of this disease.

1. Introduction

Vitamin D (VD) is a fat-soluble vitamin with multiple functions, the most important is bone mineralization by regulating the levels of calcium and phosphorus in the bone matrix. A VD deficiency is therefore related to bone deformation (rickets) or pain (osteomalacia) [1,2]. This deficiency can be caused by an insufficient diet supply and a low sun exposure, but also by the inadequate absorption of VD from the diet, or abnormal conversion of VD to its bioactive metabolites [3]. It is well-known that there are inter-individual differences in the absorption of VD from the diet that depend on the molecular forms of VD (whether it is vitamin D3 from animal sources or vitamin D2 from plant-based foods) [4], the type of food (the content of fats, dietary fiber, etc.) [4,5], or the chemical modifications of VD in supplements and fortified foods [6]. Furthermore, other endogenous factors, such as age [7], obesity [8], diseases that imply gastrointestinal or renal disorders [9] and genetic variations [5,8] also define the individual’s VD status. Regarding the latter, it has been estimated that genetic variants (including single nucleotide polymorphisms, SNPs) associated with VD explain about 20% of the heritability to the serum levels of this vitamin [10]. Thus, there are several risk factors of VD deficiency, among which the genetic factors seem to be key. Indeed, there have been identified over 60 genetic loci that seem to influence the VD levels in the blood, i.e., the bioavailable part of VD [10,11]. Fewer genetic variants seem to be strictly related to a VD deficiency in the population, according to recent genome-wide association studies (GWAS) that sought the genetic determinants of this phenotype [12].
VD levels of less than 20 ng/mL in serum indicate a VD deficiency [13]. VD is measured in the form of 25-hydroxyvitamin D 25(OH)D, also known as calcidiol, which reflects VD produced endogenously, as well as that obtained from foods and supplements. Moreover, this regards the most active metabolite of VD [14]. While dietary reference intakes (DRIs) have been established for VD, to prevent its deficiency [15], there is still scant evidence to conclude that adequate intake levels of this vitamin prevent potential adverse health outcomes other than bone health [16]. Despite the absence of information, in recent years there have been plenty of studies on the beneficial role of VD in cancer [17], cardiovascular disease (CVD) [18], diabetes mellitus [19] and obesity [20], amongst other diseases. Moreover, low VD levels have been related to a higher prevalence of comorbidities in the clinical setting [21]. Together, these studies have shown that VD has many biological functions beyond regulating the skeletal homeostasis [22]. In fact, other well-recognized actions of VD are related to cell growth, inflammatory processes, and the immune response in the body [2]. Due to its interaction with the immune system, modulating both the innate and adaptive immune systems, it has become more evident that VD is also involved in infectious diseases triggered by viruses, such as the human immunodeficiency virus (HIV) and the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that causes the COVID-19 disease [23,24].
The COVID-19 disease has been of global concern since the proclamation of this pandemic by the WHO [25]. COVID-19 causes diseases including, but not restricted to, respiratory and non-respiratory diseases, e.g., neurological complications and cardiovascular diseases, and conditions related to hypertensive disorders, amongst others [26,27,28]. Its pathogenesis has been elucidated to a great extent and several risk factors of this disease have been proposed, including the A type blood group, the presence of comorbidities and some genetic factors [29,30]. For instance, the variant rs17047200 in the TLL-1 gene seems to confer a higher risk of a severe COVID-19 disease [31,32,33] Regarding the comorbidities, certain diseases in COVID-19 patients are extremely worrying since deleterious effects have been observed in multimorbid patients, if infected by the virus [34,35,36]. Indeed, the prevalence of comorbidities is relatively high among those who develop severe COVID-19 [37,38,39]. Interestingly, some studies have also shown that COVID-19 patients tend to have deficient VD levels and a worse disease progression [39,40,41].
Nowadays, there are various biological databases that constitute large repositories on genes and associated diseases, and cellular, molecular and metabolic pathways linking genes with diseases. These tools allow for conducting bioinformatics studies by integrating and analyzing data from different sources, thereby facilitating the interpretation and visualization of the results. To the best of our knowledge, an attempt has been made to apply such tools to investigate the connection between comorbidities in the COVID-19 disease [42]. However, the role of VD was not accounted for in this study.
Thus, the objective of the current study was to interrogate the biological databases on the potential gene-disease associations between genetic factors related to a VD deficiency and the comorbid conditions and to shed light on the underlying mechanisms between VD and the comorbidities in the COVID-19 disease. We hypothesized that a susceptibility to a VD deficiency and/or to COVID-19, could increase the risk to develop more severe and fatal COVID-19 outcomes, given that comorbidities are common to both.

2. Results

2.1. Description of Genetic Variants of a VD Deficiency and COVID-19

The 20 SNPs associated with a VD deficiency, at a p-value threshold of 5 × 10−8, were annotated to 26 genes. Some variants were allocated in several genes of the UDP Glucuronosyltransferase Family 1 and 2. The selected variants and their corresponding genes are listed in Table 1.
Some genes to highlight are GC (VD binding protein gene), CALCB (calcitonin related polypeptide beta gene), and LIPC (lipase C, hepatic type gene). Furthermore, the 26 genes were distributed among chromosomes 1, 2, 4, 8, 11, 12, 15, 19 and 20. Chromosomes 4 and 11 showed variants with a high degree of correlation (R2 > 0.8) between them, i.e., some genetic variants were in the linkage disequilibrium (LD). For instance, in chromosome 4 (Figure S1A) the rs2282679 and rs11723621 variants were in the LD. The variants rs4944958 and rs12785878, located in chromosome 11, were also in the LD (R2 = 0.987) (Figure S1B). Fourteen variants were intronic.
Concerning COVID-19, 57 genetic variants were obtained, which were annotated to 46 genes (Table S1). The genes ACE2 (angiotensin converting enzyme 2 gene) and genes of the major histocompatibility complex (HLA-DRB1) were of particular interest. All of the 46 genes were distributed across 17 chromosomes. Some variants in chromosomes 3, 9, 19 and 21 were in the LD. In particular, in chromosome 3 (Figure S2), eleven pairs of variants with an R2 greater than 0.8 were found (for example, rs71325088/rs13078854, rs10490770/rs13078854, rs10490770/rs71325088). In chromosome 9 (Figure S3A), chromosome 19 (Figure S3B) and chromosome 21 (Figure S3C), the variant pairs rs912805253/rs657152; rs2109069/rs12610495, and rs13050728/rs17860115, respectively, were in the LD, too. Forty-one variants were intronic.
Complete information for both sets of genetic variants, regarding the relative allele frequencies (RAF) and effect sizes of the associations, are shown in Tables S2 and S3.
Interestingly, there were no overlapping genes between a VD deficiency and COVID-19.

2.2. Functional Analyses and Enrichment Analyses

The resources used for the functional and enrichment analyses (FUMA, REVIGO and DAVID), derived information on gene expression patterns and shared molecular functions between the genes, regarding both cellular/molecular functions and metabolic pathways.
For the VD deficiency genes, the heat map illustrating the expression profiles by genes and tissues (Figure 1) showed that the genes SMARCA4, PDE3B, NADSYN1, ALDH1A2 and CYP2R1 were highly expressed in adipose tissue, whereas SMARCA4, NADSYN1 and, the CYP2R1 genes were expressed in tissues of the circulatory system, especially in arteries, and in tissues of the immune system. Furthermore, the genes GC, HAL, LIPC, SULT2A1 and UGT2B7, UGT1A8 and UGT1A9 showed a notable expression in the liver. In the lung, the genes SMARCA4 and NADSYN1 showed a relatively high expression. Likewise, for the COVID-19 genes, the expression profiles in the different tissues (Figure S4) showed that the HLA-DRB1, LAMP2 and PLSCR1 genes were more expressed in the adipose tissues; HLA-DRB1, LAMP2, PLEKHA4, PLSCR1, RAVER1, TES and THBS3 genes were more represented in tissues of the circulatory system, especially in arteries and HLA-DRB1, HLA-DQA1, IFNAR2, OAS1, RAVER1, TES and TYK2 genes were more represented in tissues of the immune system.
The results on the enrichment analyses performed on the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene ontology (GO) databases, for the genes associated with a VD deficiency, are shown in Figure 2. The enrichment analyses included all of the genes associated with a VD deficiency and showed that some metabolic pathways were overrepresented (Figure 2A). In particular, it turned out that those related to the metabolism of sugars and carbohydrates (e.g., the ascorbate and aldarate metabolism using KEGG) were enriched with respect to all background genes. In fact, the genes of the UDP Glucuronosyltransferase Family 1 and Family 2 explained 30% of the genes present in this pathway. These genes were also enriched in the retinol metabolism pathway, together with gene ALDH1A2, as well as in the pathways related to the cytochrome P450 system. A total number of 17 genes were also involved in 25 biological processes, according to the GO pathways (Figure 2B), including the uronic acid metabolic processes, and the metabolic processes related to xenobiotics, drugs, steroids and other hormones, as well as lipids and flavonoids. The molecular functions in MsigDB were related to these processes, too; for example, the glucuronosyltransferase activity emerged as the most important function (data not shown). The Reactome database showed an enrichment for glucuronidation, the VD metabolism, the metabolism of lipids and steroids, and biological oxidations (data not shown). Similarly, in WikiPathways we observed an enrichment for glucuronidation, the VD metabolism, the oxidation by cytochrome P450 and the pathways related to nuclear receptors and the transcription factors NRF2. Some of these processes were semantically similar to the GO terms in REVIGO, namely the metabolic processes of vitamins, steroids, flavonoids, lipids and xenobiotics (Figure S5A).
Regarding COVID-19 (Figure 3), we observed that the selected genes were involved in the immunological processes, such as antigen processing and presentation pathways, and immune regulation pathways (e.g., of natural killer cells, and cytokines). Among them, pathways related to autoimmune diseases and asthma were enriched in the KEGG database (Figure 3A). Enriched signaling pathways were also restricted to the immunological processes (Figure 3B). The HLA genes (the gene family of the major histocompatibility complex, MHC), together with the IL10RB and IFNA10 genes, mainly, were driving these pathways. Of note, no molecular function was enriched in the GO database. The pathway maps of COVID-19 genes showed interactions between these genes and cytokines, as well as processes that comprise the immune recognition against viral infections (Figure S6). In addition, in the MsigDB the biological processes and immunologic signatures that appeared significantly enriched, were also linked to pathways mediated by cytokines and other immunologic cells (for example, CD4 T cells and macrophages). There were two main semantically similar GO terms: negative regulation of the viral genomes and receptor signaling pathways (Figure S5B).

2.3. Gene-Disease Associations

The platforms and tools used for the gene-disease association analyses (DisGeNET and ClusterProfiler) produced gene-disease clusters, networks and heat maps by diseases and groups of diseases. Results on the variant-disease associations are not shown for the sake of clarity and because the results were fairly similar.
For a VD deficiency, there were 14 genes matching gene IDs in the DisGeNET database: CYP2R1, CYP24A1, NADSYN1, HAL, ALDH1A2, LIPC, SMARCA4, SULT2A1, BCAS1, UGT1A9, UGT1A6, UGT2B7, GC and PDE3B. The gene-disease networks of these genes (Figure S7), showed that there were strong interactions between the genes GC, ALDH1A2 and the Cytochrome P450 Family 2 genes, with several chronic diseases, including neurodegenerative diseases, cardiovascular diseases, musculoskeletal diseases, renal kidney diseases and various cancer types. Separate clusters were seen for the remaining genes. For instance, the gene LIPC formed a cluster of metabolic diseases. The gene-disease heat maps (Figure S8) and gene-disease heat maps by groups of diseases (Figure 4A) provided a deeper insight into these associations, revealing that genes of a VD deficiency have been mostly associated with diseases, such as diabetes and insulin resistance (Figures S7 and S8).
By the disease groups, a higher scoring (strength of association) was noted for neoplasms, congenital diseases (mostly driven by the NADSYN1 gene), male and female urogenital diseases and nutritional/metabolic diseases. Other diseases received a lower score. Importantly, the LIPC gene was more prominently associated with nutritional/metabolic diseases, with the genes GC, HAL and the Cytochrome P450 Family also being implicated.
For COVID-19, 29 genes were available in the DisGeNET database. Networks (Figure S9) and heatmaps of gene-disease associations, at the individual level (Figure S10) and by groups of diseases (Figure 4B) showed that: the five top gene-disease associations were congenital diseases, digestive system diseases, neoplasms, respiratory tract diseases and infections. Among the genes associated with respiratory tract diseases, the most significant ones were ATP11A and DPP9. Less importantly, though, were urogenital, nervous system and cardiovascular diseases. Associations with nutritional and metabolic diseases were also less important, and were mainly linked to the genes RMST, HK1, HLA-DQA1 and LAMP2. In diseases of the immune system, there were also highly relevant genes, such as HLA-DQA1, HLA-DRB1 and TYK. In cardiovascular diseases, the genes FBRSL1, IL10RB, ACE2 and ABO were those that were outstanding (Figure S10).
The results on the enrichment analysis for the genes associated with a VD deficiency are shown in Table 2. This analysis revealed the diseases in which the studied genes were more enriched, that is, more overrepresented. VD deficiency genes were firstly overrepresented in VD-related conditions (VD status and rickets). In addition, these genes were enriched for some metabolic diseases (for example, metabolic syndrome; corrected p-value = 2.933 × 10−5 and cardiovascular diseases (coronary heart disease; corrected p-value = 1.732 × 10−4).
The results on the enrichment analysis for the genes associated with COVID-19 are shown in Table 3. COVID-19-related genes were more enriched in conditions of the liver or respiratory tract (such as coughing and rales), and to a certain degree with autoimmune-related diseases.
Using the ClusterProfiler tool, the enrichment analysis was extended further to 50 diseases. The enrichment plots are shown in Figure 5. Several cancer types, metabolic syndrome and its components (such as obesity, diabetes, hypertension) emerged first. VD deficiency genes were highly represented in these diseases (Ratio > 0.5). Interestingly, some of these diseases (diabetes mellitus, coronary heart disease and cardiovascular diseases) were also significantly enriched (corrected p-value < 0.05, and Ratio > 0.3) for COVID-19-associated genes. Thus, there were common comorbid conditions enriched in genes associated with both a VD deficiency and COVID-19. Some viral infections and respiratory tract diseases were also apparent, the number of enriched genes was also high for these diseases in COVID-19.

3. Discussion

In this study, using several biological databases and tools, we have analyzed whether genes associated with a VD deficiency and the COVID-19 disease were interrelated, and whether molecular, cellular and metabolic pathways and diseases connected to these genes are shared between the two. Whereas there were no common genes, this study shows that the genes associated with a VD deficiency were mainly involved in the biological processes and pathways related to the metabolism of carbohydrates, retinol, steroids, drugs and xenobiotics. A VD deficiency is also more prominently associated with nutritional and metabolic diseases, all being comorbid conditions in the COVID-19 disease. Likely, genes associated with this illness tended to be linked to similar diseases, suggesting that a genetic predisposition to a VD deficiency might have a significant impact on COVID-19-related outcomes. Gene ontologies and pathways for COVID-19 were more predominantly related to natural killer cells and the production of cytokines, as well as to other biological processes implicated in the inflammatory and immune responses that are activated by this disease.
An earlier bioinformatics study on the potential genes associated with COVID-19 and comorbidities (hypertension, diabetes mellitus and coronary artery disease) found that there were eleven shared genes (TLR4, NLRP3, MBL2, IL6, IL1RN, IL1B, CX3CR1, CCR5, AGT, ACE and F2) [42]. These genes were enriched in nearly 300 biological pathways and 15 molecular functions, all of which are related to immune and inflammatory processes, which reinforces the importance of these mechanisms in the COVID-19 disease. Likewise, the authors showed that the three comorbidities were also connected to metabolic and immune system alterations. In our study, we did not consider genes associated with commonly present comorbidities in COVID-19, since our focus was put on a VD deficiency. Furthermore, it is important to highlight that, unlike this previous study, we kept in the analyses, only the genetic variants that met the p-value threshold of the GWAS studies, to retain the true positives, i.e., those robustly associated variants with the trait of interest. As generally known, to account for multiple testing in the GWAS studies, it is convenient to keep a stricter p-value threshold [43]. Despite this issue, our results also support that the SNPs associated with COVID-19 were implicated in the immune system function and inflammation-related diseases.
Among the best-established genetic variants in the COVID-19 disease, there are those of the gene ACE. This gene encodes the angiotensin-converting enzyme in the renin-angiotensin-aldosterone system (RAAS), which is involved in the blood pressure regulation and electrolyte balance, and possibly also in the pathophysiology of COVID-19 [44]. Genetic variants of the ACE gene affecting the treatment response of ACE inhibitors have been also described [45]. Moreover, exacerbation of symptoms in COVID-19 has been associated with a genetic predisposition to severe illness via this gene in multiple studies [46,47]. According to these studies, the alterations in the expression of the ACE gene lead to the accumulation of angiotensin, which seem to promote the infection of SARS-CoV-2 [46,48]. It is important to highlight that VD plays a central role in the RAAS, since it regulates the expression and production of renin [49]. The fact that VD and the RAAS are connected might explain why VD deficiency genes were enriched in diseases including but not limited to hypertension and coronary heart disease, in our study. Other COVID-19-related genes that were implicated in nutritional/metabolic diseases and other comorbid conditions (the ABO gene and others) seem to have no clear relationship with the metabolism and actions of VD [50]. The expression of the immune system-related SNPs in COVID-19 (the IFNA, IL10 and HLA genes), however, could be affected by the levels of VD. High levels of VD are known to suppress the B cell and T cell proliferation, the monocyte production of inflammatory cytokines (interleukins and tumor necrosis factor-alpha, TNF-α) and the expression of the MHC class II molecules [51]. Thus, given the importance of VD in the modulation of the immune function, it has been proposed that VD could have a positive impact on the COVID-19 disease [52,53]. In the last few years, various systematic reviews and meta-analyses have updated the evidence on the potential effects of a VD supplementation in patients diagnosed with COVID-19 in the clinical setting [54,55,56]. Furthermore, lower inflammatory markers have been observed in COVID-19 patients treated with VD [57]. However, not all studies support that VD supplementation is an effective treatment to reduce the risk of severe COVID-19 and its mortality. Therefore, there is still no solid evidence on this issue [56]. Nevertheless, a profound VD deficiency has been observed in patients with a severe COVID-19 illness [40,41], supporting that VD could reduce the harmful effects of the SarsCov-2 infection. As far as we are aware, whether these COVID-19 patients present a genetic susceptibility to a VD deficiency, has been not investigated.
Regarding a VD deficiency, the genes LIPC, HAL, GC and the Cytochrome P450 Family genes were among those which are more strongly associated with nutritional/metabolic diseases. The LIPC gene encodes the hepatic triglyceride lipase, which is expressed in the liver. It catalyzes the hydrolysis of triglycerides and phospholipids present in the circulating plasma. Its mutation increases the likelihood to develop dyslipidemia, hypertension and insulin resistance, and in turn, the metabolic syndrome [58]. Carriers of the risk alleles of the GC gene have been also associated with a higher risk of developing this syndrome and its components in various study populations [59,60]. A variation in the GC gene has also been associated with obesity [61]. As we have shown, this gene is prominently expressed in the adipose tissue, where VD can also regulate the biological processes including adipogenesis, inflammation, oxidative stress and metabolism in mature adipocytes [61]. Accordingly, low VD levels can lead to obesity, this is another component of the metabolic syndrome and an important cardiovascular disease risk factor [22,61]. As for the other genes, there is less evidence supporting their role in the metabolic syndrome or related components. Overall, however, all of these VD deficiency genes have not been linked to COVID-19, in any study. Thus, while there is no a common genetic background between VD and COVID-19, our results support that it is plausible that a VD deficiency causes comorbid conditions, posing a higher risk to develop a severe COVID-19 disease. In fact, as we have already mentioned, the prevalence of comorbidities in hospitalized patients with a severe COVID-19 is high [34]. VD could be key to treat the COVID-19 disease in these patients [5]. For instance, VD has been shown to improve the prognosis of COVID-19 in patients with renal kidney disease, which is also determined by the metabolic syndrome risk factors [62].
There are some limitations of this study. Due to a lack of matching gene IDs in some databases, because their associations are still uncertain in these genetic/biological data repositories or collections, we were not able to evaluate the connections with other potential pathways and diseases. Furthermore, we did not consider the gene-environment interactions or the potential confounders, such as smoking, in our analyses, which might have affected the results obtained. However, to our knowledge, this issue is not accounted for in any of these resources. Moreover, our search of genetic variants was restricted to a VD deficiency; those for which a VD deficiency was not the reported trait, were removed from the analyses. This way, the genetic variants that have been associated with the VD levels in blood in the GWAS studies, but not strictly with a VD deficiency, were not included. For example, the gene VDR, which encodes for the VD receptor (VDR), is activated by 1,25(OH)D and forms a heterodimer with the retinoid X receptor to regulate the transcription of other genes. The VDR has been found in virtually all cell types, which may explain its multiple and novel actions of VD on different tissues. Indeed evidence of the extraskeletal effects of VD, including xenobiotic detoxification, oxidative stress reduction, neuroprotective functions, antimicrobial defense, immunoregulation, anti-inflammatory/anticancer actions and cardiovascular benefits, is mediated by VDR [22]. Importantly, the expression of VDR depends on the levels of VD; the expression is only decreased when the VD levels are deficient [63]. Therefore, polymorphisms of the VDR gene were not included in this study. To check the extent to which this gene could influence our results, we carried out a secondary analysis involving all genes, including the VDR gene, where we observed that the genetic variants of the latter were overrepresented in the bone-related diseases (rheumatoid arthritis, rickets and osteoporosis), in various neoplasms and hypertensive disease. Furthermore, with a similar intensity to the GC gene, the gene VDR seemed to be involved in nutritional and metabolic diseases. Moreover, the same diseases came up in the disease-enrichment analyses and similar processes emerged in other functional analyses (Figure S11). Thus, not having considered the VDR genetic variants should not have materially affected our results.
Concerning the strengths of this study, we have used several resources and web-based tools to research the relationship between SNPs, genes, molecular functions, pathways and diseases. The use of several tools has enabled us to confirm the findings and to enhance the interpretations and our understanding on these relations. Finally, we evaluated whether proxy SNPs, i.e., SNPs in the LD among the selected genes and chromosomes, had any redundancy in our variant-disease analyses, possibly affecting multiple testing corrections in the enrichment analyses. As aforementioned, there were few SNP pairs in the LD; removing those SNPs did not alter the results. Furthermore, the enrichment and pathway analyses were undertaken on all genetic variants (and annotated genes), regardless of the population ancestry, since we included both those of European and non-European ancestry. It is unlikely that this issue has affected our results, given that the molecular functions, pathways and disease phenotypes are overall highly conserved [64]. Finally, we set the p-value threshold to 5 × 10−8, to keep only the true associations between the genetic variants with a VD deficiency and COVID-19 for the analyses.

4. Materials and Methods

4.1. Selection of Genetic Variants

We used two databases to identify SNPs related to a VD Deficiency and COVID-19: the GWAS catalog (https://www.ebi.ac.uk/gwas/ (accessed on 21 June 2022) and PhenoScanner (http://www.phenoscanner.medschl.cam.ac.uk/ (accessed on 21 June 2022)) database.
Both the GWAS catalog and PhenoScanner offer a high-quality curated collection of all published GWAS studies [65]. Specifically, the goal of these resources is to facilitate the cross-reference of genetic variants with a wide range of phenotypes. For example, PhenoScanner contains more than 65 billion associations and more than 150 million unique genetic variants [66,67]. We sought in both genomic databases, the SNPs associated with a VD deficiency or COVID-19 with a p-value lower than 5 × 10−8.
We extracted 20 SNPs (as of June 2022) associated with a VD deficiency, which came from seven GWAS studies. Of this set of SNPs, 17 SNPs were from populations with European ancestry and three SNPs were from Asian populations.
Regarding COVID-19, we applied the same steps. In this case, the number of selected variants was 57 (as of August 2022). Likewise, the majority of those SNPs were identified in the European GWAS studies.

4.2. Annotation of Genes

We used the GWAS catalog [65] and Ensembl [68,69] to annotate the SNPs previously identified to their genes. The information related to the variants of interest was extracted from the GWAS catalog. Indeed, the GWAS catalog database provides information about the genes to which the variants of interest are associated, the chromosome and its position, amongst the other data of interest. Ensemble (https://www.ensembl.org/index.html (accessed on 3 November 2022)) was used to complement this information via the variant effect predictor (VEP) application (https://www.ensembl.org/info/docs/tools/vep/index.html accessed on 22 June 2022) [69], and GeneCards (https://www.genecards.org/ (accessed on 22 June 2022)) was used to complete the name of each gene. These tools allowed us to retrieve information on the location of the variant and on the consequences of the variation in the DNA sequence (whether the variant is intronic, missense, etc.). The GRCh38.p13 genome was used as a reference for the gene annotations.
The same procedure was carried out in both cases, for a VD deficiency and for COVID-19.

4.3. Linkage Disequilibrium Analysis

We subjected all variants to the LD analysis to explore the existence of closely related genetic variants within a chromosome [70]. Given that the alleles of such variants are inherited together more commonly than expected, they are assumed to be correlated; the higher the correlation (R2), the closer the variants [71]. These variants could therefore imply a redundancy in the analysis, if not accounted for. We used the LD matrix tool provided by LDLINK (National Cancer Institute, NHI). (https://ldlink.nci.nih.gov/?tab=ldmatrix (accessed on 22 June 2022)) to identify the genetic variants in the LD. An interactive heat map matrix of the pairwise LD statistics for each chromosome in which more than one variant associated with a VD deficiency or COVID-19 was present, was obtained.

4.4. Identification of the Overlapping Genes

The related genes between a VD deficiency and COVID-19 were cross-analyzed, using custom functions developed using R software version 4.1.2 [72]. More precisely, two lists of genes and variants, one with the VD deficiency genes and the other with the COVID-19 genes were crossed to identify the common variants and genes. Venn diagrams were created to further visualize the potential overlapping genes.

4.5. Biological Database Studies

We performed the functional and enrichment analyses on the list of genes and variants selected in the previous steps, using several tools and resources. The type of analysis and the tools used to carry out these analyses are described below.

4.5.1. In Silico Functional Analyses and Enrichment Analyses

First, we run functional and enrichment analyses using DAVID, FUMA and REVIGO.
DAVID is a bioinformatics resource that enables the functional annotation and enrichment analysis of gene lists. This tool allows to agglomerate the gene identifiers from resources, such as Ensembl, Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) [73]. Through DAVID (https://david.ncifcrf.gov/tools.jsp (accessed on 23 June 2022)) we gathered information on the diseases that are related to the selected genes and variants, the functional annotations and pathways.
With FUMA, using the GENE2FUNC module (https://fuma.ctglab.nl/gene2func (accessed on 11 July 2022)), we obtained information on the expression patterns of the genes and shared molecular functions between the genes. In particular, we obtained the gene expression heat maps and results on the functional analyses for a given molecular, cellular or metabolic pathway, and even the presence in certain diseases. Additionally, the expression of the genes of interest in different tissues (from Gtex v6 data) [74] were extracted to compare the expression profiles. FUMA is another web-based platform that can be used to prioritize, visualize and interpret the GWAS results [75]. Its main objective is to generate functional biological information on the GWAS hits. It also provides gene-based, pathway and tissue enrichment results. The platform is linked to popular resources, such as GO, WikiPathways and MsigDB.
In the enrichment analysis, hypergeometric tests were performed to check if the genes of interest are overrepresented in a predefined gene set. In our case, the latter gene set, i.e., the so-called background genes, were the “universe”. The gene sets were compared with the aforementioned background genes regarding the pathways and diseases using different data repositories (MsigDB, WikiPathways, GO and KEGG). GO divides the functions of genes into three parts: cellular component, molecular function, and biological process. KEGG provides information on the metabolic pathways and also annotates enzymes that catalyze each step of these pathways. The Molecular Signatures Database (MsigDB) contains a collection of canonical pathways derived from KEGG, Reactome, WikiPathways and others. The immunologic signature gene set of MsigDB was also interrogated via FUMA. All of these databases were accessed in July 2022.
REVIGO is a web platform that uses GO terms from previously conducted enrichment analyses for further visualization and interpretation of the results. It uses a clustering algorithm, based on semantic similarity measures, giving several outputs to help in the interpretation [76].
All GO terms were obtained from the enrichment analyses conducted on the genes of interest, using ClusterProfiler package in the R software [77,78]. These analyses were performed separately for the VD deficiency and COVID-19 genes, and by combining the two. In addition to the GO terms, we obtained the p-values that indicated the degree of association between the gene and the GO term.
In REVIGO (http://revigo.irb.hr/FAQ.aspx#q01 (accessed on 15 July 2022)), each GO term, together with its p-value, was introduced, to obtain the different graphs for the interpretation of the results (multidimensional scaling plots, tree maps and tag clouds).

4.5.2. SNP-Disease and Gene-Disease Associations, including Enrichment Analyses

To explore the associations between the VD deficiency and COVID-19-associated genes, we used the DisGeNET database in R software (https://www.disgenet.org/disgenet2r (accessed on 3 November 2022)) from the disgenet2r package [79]. This tool allows to retrieve, visualize and expand the variant-disease and gene-disease associations.
The degree of association between a certain SNP and a disease, or a gene and a disease, is shown through a “score”. Indeed, the DisGeNET score ranks the gene-disease associations, according to their level of evidence calculated by an algorithm that considers the number and type of sources present in the database, and the number of publications that support the association [80,81].
The list of genes of interest was introduced as the HGNC symbols or Entrez gene identifiers. These results were represented as gene-disease heat maps; gene-class of diseases heat maps; or gene-disease networks. In this study, the DisGeNET curated database was used, which contains 9703 genes, 11,181 diseases and 84,038 associations (https://www.disgenet.org/dbinfo). This database was also accessed in 30 July 2022.
In addition, an enrichment analysis on the diseases in DisGeNET was performed for the same input of genes and the curated database available in DisGeNET. In the enrichment analysis, again, Fisher’s test was applied to all genes (universe; N = 9703) concerning the set of genes that were previously identified (VD deficiency and/or COVID-19). The p-values resulting from the multiple Fisher tests were corrected by the Benjamini–Hochberg (BH) method and for the false discovery rate using a FDR p-value of 0.05 [82]. The number of genes from the input list that were associated with the outcome reflects the BgRatio.

4.6. Workflow of the Analysis

Figure 6 shows our workflow diagram, illustrating the different steps and tools used to carry out our study. In summary, we started the analysis by extracting the variants associated with a Vitamin D (VD) deficiency and COVID-19, from the GWAS catalog and PhenoScanner database (p-value threshold 5 × 10−8) [65,67]. Subsequently, we used the Variant Effect Prediction tool (VEP, in Ensembl) to annotate these variants to genes [68]. The linkage disequilibrium (LD) of the selected variants was examined [70,71], as well as the existence of the overlapping variants or genes. Then, a study of the biological databases was carried out; firstly, an in silico functional analysis and enrichment analysis (using the platforms DAVID [73], FUMA [75] and REVIGO [76]) was conducted, and secondly, a study of the SNP-disease and gene-disease associations, including enrichment analyses (using disgenet2r package and ClusterProfiler) was carried out.

5. Conclusions

In the current study, we have made use of the genetic variants of a VD deficiency and COVID-19 disease to explore the molecular functions, metabolic pathways, and diseases relevant to both a VD deficiency and COVID-19. Using several complementary approaches, we provide insight into the link between a genetic susceptibility to low levels of VD and to the COVID-19 disease. Specifically, our results support the hypothesis that a genetically predicted VD deficiency might imply a higher risk of developing comorbid conditions, and in turn, a severe or fatal COVID-19. Therefore, the role of VD in this context should be addressed in future studies to evaluate whether the supplementation of VD in the population susceptible to a VD deficiency could reduce the COVID-19 severity. Furthermore, molecular and functional studies are needed to evidence the molecular mechanisms by which a VD deficiency and the COVID-19 disease are connected. This knowledge is also crucial to implement efficient strategies to prevent and control this disease.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232214256/s1. Figure S1: Linkage disequilibrium for vitamin D deficiency; Figure S2: Linkage disequilibrium for COVID-19 (chromosome 3); Figure S3: Linkage disequilibrium for COVID-19 (chromosome 9, 19 and 21); Figure S4: Expression heat map by genes and tissues (COVID-19); Figure S5: REVIGO for vitamin d deficiency and covid 19; Figure S6: Enrichment analysis of COVID-19 genes regarding GO terms in MSigDB on immunological signatures; Figure S7: Vitamin D deficiency gene-disease network; Figure S8: Vitamin D deficiency gene-disease heat map; Figure S9: COVID-19 gene-disease network; Figure S10: COVID-19 gene-disease heat map; Figure S11: Enrichment analysis of VD deficiency and the VDR gene; Table S1: Genes and variants associated with COVID-19; Table S2: GWAS catalog of vitamin D deficiency; Table S3: GWAS catalog of COVID-19.

Author Contributions

Conceptualization, E.M.-M., M.R. and M.Á.J.-S.; methodology, E.M.-M. and Á.A.-S.; formal analysis, Á.A.-S.; investigation, Á.A.-S. and E.M.-M.; data curation, Á.A.-S.; writing—original draft preparation, Á.A.-S. and E.M.-M.; writing—review and editing, Á.A.-S., M.R., M.R.-B., M.Á.J.-S., Á.G., M.J.S. and E.M.-M.; supervision, E.M.-M., M.R.-B. and M.J.S.; project administration, E.M.-M.; funding acquisition, E.M.-M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Project PECOVID-0200-2020, funded by Consejería de Salud de la Junta de Andalucía and co-funded by the European Regional Development Fund (ERDF-FEDER).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Links to resources and tools used in this study have been provided in the text. Datasets (list of genes) used are provided in the supplementary material.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of the data; in the writing of the manuscript; nor in the decision to publish the results.

References

  1. DeLuca, H.F. Metabolism of Vitamin D: Current Status. Am. J. Clin. Nutr. 1976, 29, 1258–1270. [Google Scholar] [CrossRef] [PubMed]
  2. Pike, J.W.; Christakos, S. Biology and Mechanisms of Action of the Vitamin D Hormone. Endocrinol. Metab. Clin. N. Am. 2017, 46, 815–843. [Google Scholar] [CrossRef] [PubMed]
  3. Kumar Maurya, V.; Aggarwal, M. Factors Influencing the Absorption of Vitamin D in GIT: An Overview. J. Food Sci. Technol. 2017, 54, 3753–3765. [Google Scholar] [CrossRef]
  4. Weber, F. Absorption Mechanisms for Fat-Soluble Vitamins and the Effect of Other Food Constituents. Prog. Clin. Biol. Res. 1981, 77, 119–135. [Google Scholar] [PubMed]
  5. Niramitmahapanya, S.; Harris, S.S.; Dawson-Hughes, B. Type of Dietary Fat Is Associated with the 25-Hydroxyvitamin D3 Increment in Response to Vitamin D Supplementation. J. Clin. Endocrinol. Metab. 2011, 96, 3170–3174. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Natri, A.-M.; Salo, P.; Vikstedt, T.; Palssa, A.; Huttunen, M.; Kärkkäinen, M.U.M.; Salovaara, H.; Piironen, V.; Jakobsen, J.; Lamberg-Allardt, C.J. Bread Fortified with Cholecalciferol Increases the Serum 25-Hydroxyvitamin D Concentration in Women as Effectively as a Cholecalciferol Supplement. J. Nutr. 2006, 136, 123–127. [Google Scholar] [CrossRef] [Green Version]
  7. Clemens, T.L.; Zhouf, X.-Y.; Myles, M.; Endres, D.; Lindsay, R. Serum Vitamin D2 and Vitamin D3 Metabolite Concentrations and Absorption of Vitamin D2 in Elderly Subjects. J. Clin. Endocrinol. Metab. 1986, 63, 656–660. [Google Scholar] [CrossRef]
  8. Wortsman, J.; Matsuoka, L.Y.; Chen, T.C.; Lu, Z.; Holick, M.F. Decreased Bioavailability of Vitamin D in Obesity. Am. J. Clin. Nutr. 2000, 72, 690–693. [Google Scholar] [CrossRef] [Green Version]
  9. Farraye, F.A.; Nimitphong, H.; Stucchi, A.; Dendrinos, K.; Boulanger, A.B.; Vijjeswarapu, A.; Tanennbaum, A.; Biancuzzo, R.; Chen, T.C.; Holick, M.F. Use of a Novel Vitamin D Bioavailability Test Demonstrates That Vitamin D Absorption Is Decreased in Patients with Quiescent Crohnʼs Disease. Inflamm. Bowel Dis. 2011, 17, 2116–2121. [Google Scholar] [CrossRef]
  10. Manousaki, D.; Mitchell, R.; Dudding, T.; Haworth, S.; Harroud, A.; Forgetta, V.; Shah, R.L.; Luan, J.; Langenberg, C.; Timpson, N.J.; et al. Genome-Wide Association Study for Vitamin D Levels Reveals 69 Independent Loci. Am. J. Hum. Genet. 2020, 106, 327–337. [Google Scholar] [CrossRef]
  11. Fu, L.; Yun, F.; Oczak, M.; Wong, B.Y.L.; Vieth, R.; Cole, D.E.C. Common Genetic Variants of the Vitamin D Binding Protein (DBP) Predict Differences in Response of Serum 25-Hydroxyvitamin D [25(OH)D] to Vitamin D Supplementation. Clin. Biochem. 2009, 42, 1174–1177. [Google Scholar] [CrossRef] [PubMed]
  12. Sepulveda-Villegas, M.; Elizondo-Montemayor, L.; Trevino, V. Identification and Analysis of 35 Genes Associated with Vitamin D Deficiency: A Systematic Review to Identify Genetic Variants. J. Steroid Biochem. Mol. Biol. 2020, 196, 105516. [Google Scholar] [CrossRef] [PubMed]
  13. Holick, M.F.; Chen, T.C. Vitamin D Deficiency: A Worldwide Problem with Health Consequences. Am. J. Clin. Nutr. 2008, 87, 1080S–1086S. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Peppelenbosch, M.P.; Fuhler, G.M.; Verhaar, A.P.; van der Laan, L.J.W. Action and Function of Vitamin D in Digestive Tract Physiology and Pathology. Curr. Med. Chem. 2017, 24, 928–936. [Google Scholar] [CrossRef]
  15. IOM (Institute of Medicine). Dietary Reference Intakes for Calcium and Vitamin D; Ross, A.C., Taylor, C.L., Yaktine, A.L., del Valle, H.B., Eds.; National Academies Press: Washington, DC, USA, 2011; ISBN 978-0-309-16394-1. [Google Scholar]
  16. Newberry, S.J.; Chung, M.; Shekelle, P.G.; Suttorp Booth, M.; Liu, J.L.; Ruelaz Maher, A.; Motala, A.; Cui, M.; Perry, T.; Shanman, R.; et al. Vitamin D and Calcium: A Systematic Review of Health Outcomes (Update). Evid. Rep. Technol. Assess. (Full Rep.) 2014, 217, 1–929. [Google Scholar] [CrossRef] [Green Version]
  17. Zhang, R.; Zhang, Y.; Liu, Z.; Pei, Y.; Xu, P.; Chong, W.; Hai, Y.; He, L.; He, Y.; Yu, J.; et al. Association between Vitamin D Supplementation and Cancer Mortality: A Systematic Review and Meta-Analysis. Cancers 2022, 14, 3717. [Google Scholar] [CrossRef]
  18. Pei, Y.-Y.; Zhang, Y.; Peng, X.-C.; Liu, Z.-R.; Xu, P.; Fang, F. Association of Vitamin D Supplementation with Cardiovascular Events: A Systematic Review and Meta-Analysis. Nutrients 2022, 14, 3158. [Google Scholar] [CrossRef]
  19. Lo, A.C.Q.; Lo, C.C.W. The Effect of Vitamin D Supplementation on Glycemic Control/Glucose Metabolism and Maternal-Neonatal Outcomes in Women with Established Gestational Diabetes Mellitus: An Updated Meta-Analysis. Clin. Nutr. 2022, 41, 2420–2423. [Google Scholar] [CrossRef]
  20. Dewansingh, P.; Reckman, G.A.R.; Mijlius, C.F.; Krijnen, W.P.; van der Schans, C.P.; Jager-Wittenaar, H.; van den Heuvel, E.G.H.M. Protein, Calcium, Vitamin D Intake and 25(OH)D Status in Normal Weight, Overweight, and Obese Older Adults: A Systematic Review and Meta-Analysis. Front. Nutr. 2021, 8, 718658. [Google Scholar] [CrossRef]
  21. Furuie, I.N.; Mauro, M.J.J.; Petruzziello, S.; Riechi, S.C.; Petterle, R.R.; Boguszewski, C.L.; Borba, V.Z.C. Two Threshold Levels of Vitamin D and the Prevalence of Comorbidities in Outpatients of a Tertiary Hospital. Osteoporos. Int. 2018, 29, 433–440. [Google Scholar] [CrossRef]
  22. Gil, A.; Plaza-Diaz, J.; Mesa, M.D. Vitamin D: Classic and Novel Actions. Ann. Nutr. Metab. 2018, 72, 87–95. [Google Scholar] [CrossRef] [PubMed]
  23. Wang, Y.; Huang, X.; Wu, Y.; Li, A.; Tian, Y.; Ren, M.; Li, Z.; Zhang, T.; Wu, H.; Wang, W. Increased Risk of Vitamin D Deficiency Among HIV-Infected Individuals: A Systematic Review and Meta-Analysis. Front. Nutr. 2021, 8, 722032. [Google Scholar] [CrossRef] [PubMed]
  24. Ghasemian, R.; Shamshirian, A.; Heydari, K.; Malekan, M.; Alizadeh-Navaei, R.; Ebrahimzadeh, M.A.; Ebrahimi Warkiani, M.; Jafarpour, H.; Razavi Bazaz, S.; Rezaei Shahmirzadi, A.; et al. The Role of Vitamin D in the Age of COVID-19: A Systematic Review and Meta-analysis. Int. J. Clin. Pract. 2021, 75, e14675. [Google Scholar] [CrossRef]
  25. Cucinotta, D.; Vanelli, M. WHO Declares COVID-19 a Pandemic. Acta Biomed. 2020, 91, 157–160. [Google Scholar] [CrossRef]
  26. Tossetta, G.; Fantone, S.; delli Muti, N.; Balercia, G.; Ciavattini, A.; Giannubilo, S.R.; Marzioni, D. Preeclampsia and Severe Acute Respiratory Syndrome Coronavirus 2 Infection: A Systematic Review. J. Hypertens. 2022, 40, 1629–1638. [Google Scholar] [CrossRef]
  27. Delli Muti, N.; Finocchi, F.; Tossetta, G.; Salvio, G.; Cutini, M.; Marzioni, D.; Balercia, G. Could SARS-CoV-2 Infection Affect Male Fertility and Sexuality? APMIS 2022, 130, 243–252. [Google Scholar] [CrossRef] [PubMed]
  28. Dale, L. Neurological Complications of COVID-19: A Review of the Literature. Cureus 2022, 14, e27633. [Google Scholar] [CrossRef]
  29. Ragab, D.; Salah Eldin, H.; Taeimah, M.; Khattab, R.; Salem, R. The COVID-19 Cytokine Storm; What We Know So Far. Front. Immunol. 2020, 11, 1446. [Google Scholar] [CrossRef]
  30. Clerbaux, L.-A.; Albertini, M.C.; Amigó, N.; Beronius, A.; Bezemer, G.F.G.; Coecke, S.; Daskalopoulos, E.P.; del Giudice, G.; Greco, D.; Grenga, L.; et al. Factors Modulating COVID-19: A Mechanistic Understanding Based on the Adverse Outcome Pathway Framework. J. Clin. Med. 2022, 11, 4464. [Google Scholar] [CrossRef]
  31. Grimaudo, S.; Amodio, E.; Pipitone, R.M.; Maida, C.M.; Pizzo, S.; Prestileo, T.; Tramuto, F.; Sardina, D.; Vitale, F.; Casuccio, A.; et al. PNPLA3 and TLL-1 Polymorphisms as Potential Predictors of Disease Severity in Patients With COVID-19. Front. Cell Dev. Biol. 2021, 9, 627914. [Google Scholar] [CrossRef]
  32. Gupta, K.; Kaur, G.; Pathak, T.; Banerjee, I. Systematic Review and Meta-Analysis of Human Genetic Variants Contributing to COVID-19 Susceptibility and Severity. Gene 2022, 844, 146790. [Google Scholar] [CrossRef]
  33. van der Made, C.I.; Netea, M.G.; van der Veerdonk, F.L.; Hoischen, A. Clinical Implications of Host Genetic Variation and Susceptibility to Severe or Critical COVID-19. Genome Med. 2022, 14, 96. [Google Scholar] [CrossRef] [PubMed]
  34. Ejaz, H.; Alsrhani, A.; Zafar, A.; Javed, H.; Junaid, K.; Abdalla, A.E.; Abosalif, K.O.A.; Ahmed, Z.; Younas, S. COVID-19 and Comorbidities: Deleterious Impact on Infected Patients. J. Infect. Public Health 2020, 13, 1833–1839. [Google Scholar] [CrossRef]
  35. Bigdelou, B.; Sepand, M.R.; Najafikhoshnoo, S.; Negrete, J.A.T.; Sharaf, M.; Ho, J.Q.; Sullivan, I.; Chauhan, P.; Etter, M.; Shekarian, T.; et al. COVID-19 and Preexisting Comorbidities: Risks, Synergies, and Clinical Outcomes. Front. Immunol. 2022, 13, 890517. [Google Scholar] [CrossRef]
  36. Tian, Y.; Qiu, X.; Wang, C.; Zhao, J.; Jiang, X.; Niu, W.; Huang, J.; Zhang, F. Cancer Associates with Risk and Severe Events of COVID-19: A Systematic Review and Meta-Analysis. Int. J. Cancer 2021, 148, 363–374. [Google Scholar] [CrossRef] [PubMed]
  37. Ticinesi, A.; Parise, A.; Cerundolo, N.; Nouvenne, A.; Prati, B.; Chiussi, G.; Guerra, A.; Meschi, T. Multimorbidity and Frailty Are the Key Characteristics of Patients Hospitalized with COVID-19 Breakthrough Infection during Delta Variant Predominance in Italy: A Retrospective Study. J. Clin. Med. 2022, 11, 5442. [Google Scholar] [CrossRef] [PubMed]
  38. Ahlström, B.; Frithiof, R.; Larsson, I.-M.; Strandberg, G.; Lipcsey, M.; Hultström, M. A Comparison of Impact of Comorbidities and Demographics on 60-Day Mortality in ICU Patients with COVID-19, Sepsis and Acute Respiratory Distress Syndrome. Sci. Rep. 2022, 12, 15703. [Google Scholar] [CrossRef]
  39. Mishra, P.; Parveen, R.; Bajpai, R.; Agarwal, N. Vitamin D Deficiency and Comorbidities as Risk Factors of COVID-19 Infection: A Systematic Review and Meta-Analysis. J. Prev. Med. Public Health 2022, 55, 321–333. [Google Scholar] [CrossRef]
  40. Cereda, E.; Bogliolo, L.; Klersy, C.; Lobascio, F.; Masi, S.; Crotti, S.; de Stefano, L.; Bruno, R.; Corsico, A.G.; di Sabatino, A.; et al. Vitamin D 25OH Deficiency in COVID-19 Patients Admitted to a Tertiary Referral Hospital. Clin. Nutr. 2021, 40, 2469–2472. [Google Scholar] [CrossRef]
  41. Jude, E.B.; Ling, S.F.; Allcock, R.; Yeap, B.X.Y.; Pappachan, J.M. Vitamin D Deficiency Is Associated With Higher Hospitalization Risk From COVID-19: A Retrospective Case-Control Study. J. Clin. Endocrinol. Metab. 2021, 106, e4708–e4715. [Google Scholar] [CrossRef]
  42. Feng, S.; Song, F.; Guo, W.; Tan, J.; Zhang, X.; Qiao, F.; Guo, J.; Zhang, L.; Jia, X. Potential Genes Associated with COVID-19 and Comorbidity. Int. J. Med. Sci. 2022, 19, 402–415. [Google Scholar] [CrossRef] [PubMed]
  43. Pe’er, I.; Yelensky, R.; Altshuler, D.; Daly, M.J. Estimation of the Multiple Testing Burden for Genomewide Association Studies of Nearly All Common Variants. Genet. Epidemiol. 2008, 32, 381–385. [Google Scholar] [CrossRef] [PubMed]
  44. Yamamoto, N.; Nishida, N.; Yamamoto, R.; Gojobori, T.; Shimotohno, K.; Mizokami, M.; Ariumi, Y. Angiotensin-Converting Enzyme (ACE) 1 Gene Polymorphism and Phenotypic Expression of COVID-19 Symptoms. Genes 2021, 12, 1572. [Google Scholar] [CrossRef] [PubMed]
  45. Arnett, D.K.; Davis, B.R.; Ford, C.E.; Boerwinkle, E.; Leiendecker-Foster, C.; Miller, M.B.; Black, H.; Eckfeldt, J.H. Pharmacogenetic Association of the Angiotensin-Converting Enzyme Insertion/Deletion Polymorphism on Blood Pressure and Cardiovascular Risk in Relation to Antihypertensive Treatment. Circulation 2005, 111, 3374–3383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Saengsiwaritt, W.; Jittikoon, J.; Chaikledkaew, U.; Udomsinprasert, W. Genetic Polymorphisms of ACE1, ACE2, and TMPRSS2 Associated with COVID-19 Severity: A Systematic Review with Meta-Analysis. Rev. Med. Virol. 2022, 32, e2323. [Google Scholar] [CrossRef]
  47. Yamamoto, N.; Yamamoto, R.; Ariumi, Y.; Mizokami, M.; Shimotohno, K.; Yoshikura, H. Does Genetic Predisposition Contribute to the Exacerbation of COVID-19 Symptoms in Individuals with Comorbidities and Explain the Huge Mortality Disparity between the East and the West? Int. J. Mol. Sci. 2021, 22, 5000. [Google Scholar] [CrossRef]
  48. Khoo, S.K.; Hurst, T.; Webb, M.J.; Dickie, G.; Kearsley, J.; Parsons, P.G.; Mackay, E.V. Cisplatin Chemotherapy of Ovarian Cancer: Is Short-Term In Vitro Chemosensitivity Predictive of Long-Term Patient Survival? Aust. N. Z. J. Obstet. Gynaecol. 1988, 28, 313–317. [Google Scholar] [CrossRef]
  49. Li, Y.C. Vitamin D Regulation of the Renin–Angiotensin System. J. Cell. Biochem. 2003, 88, 327–331. [Google Scholar] [CrossRef]
  50. Hashemi, S.M.A.; Thijssen, M.; Hosseini, S.Y.; Tabarraei, A.; Pourkarim, M.R.; Sarvari, J. Human Gene Polymorphisms and Their Possible Impact on the Clinical Outcome of SARS-CoV-2 Infection. Arch. Virol. 2021, 166, 2089–2108. [Google Scholar] [CrossRef]
  51. Aranow, C. Vitamin D and the Immune System. J. Investig. Med. 2011, 59, 881–886. [Google Scholar] [CrossRef]
  52. Malek Mahdavi, A. A Brief Review of Interplay between Vitamin D and Angiotensin-Converting Enzyme 2: Implications for a Potential Treatment for COVID-19. Rev. Med. Virol. 2020, 30, e2119. [Google Scholar] [CrossRef] [PubMed]
  53. Datta, P.K.; Liu, F.; Fischer, T.; Rappaport, J.; Qin, X. SARS-CoV-2 Pandemic and Research Gaps: Understanding SARS-CoV-2 Interaction with the ACE2 Receptor and Implications for Therapy. Theranostics 2020, 10, 7448–7464. [Google Scholar] [CrossRef] [PubMed]
  54. Migliorini, F.; Vaishya, R.; Eschweiler, J.; Oliva, F.; Hildebrand, F.; Maffulli, N. Vitamins C and D and COVID-19 Susceptibility, Severity and Progression: An Evidence Based Systematic Review. Medicina 2022, 58, 941. [Google Scholar] [CrossRef]
  55. D’Ecclesiis, O.; Gavioli, C.; Martinoli, C.; Raimondi, S.; Chiocca, S.; Miccolo, C.; Bossi, P.; Cortinovis, D.; Chiaradonna, F.; Palorini, R.; et al. Vitamin D and SARS-CoV2 Infection, Severity and Mortality: A Systematic Review and Meta-Analysis. PLoS ONE 2022, 17, e0268396. [Google Scholar] [CrossRef]
  56. Bassatne, A.; Basbous, M.; Chakhtoura, M.; el Zein, O.; Rahme, M.; El-Hajj Fuleihan, G. The Link between COVID-19 and VItamin D (VIVID): A Systematic Review and Meta-Analysis. Metabolism 2021, 119, 154753. [Google Scholar] [CrossRef] [PubMed]
  57. Sharif-Askari, F.S.; Hafezi, S.; Sharif-Askari, N.S.; Alsayed, H.A.H.; Mdkhana, B.; Selvakumar, B.; Temsah, M.-H.; Saddik, B.; al Anouti, F.; Halwani, R. Vitamin D Modulates Systemic Inflammation in Patients with Severe COVID-19. Life Sci. 2022, 307, 120909. [Google Scholar] [CrossRef]
  58. Kraja, A.T.; Vaidya, D.; Pankow, J.S.; Goodarzi, M.O.; Assimes, T.L.; Kullo, I.J.; Sovio, U.; Mathias, R.A.; Sun, Y.V.; Franceschini, N.; et al. A Bivariate Genome-Wide Approach to Metabolic Syndrome: STAMPEED Consortium. Diabetes 2011, 60, 1329–1339. [Google Scholar] [CrossRef] [Green Version]
  59. Hidalgo-Bravo, A.; Rivera-Paredez, B.; León-Reyes, G.; Patiño, N.; Castillejos-López, M.; Salmerón, J.; Velázquez-Cruz, R. Unravelling the Contribution of the Rs7041 and Rs4588 Polymorphisms of the GC Gene and Serum VDBP Levels for Developing Metabolic Syndrome in the Mexican Population. Int. J. Mol. Sci. 2022, 23, 10581. [Google Scholar] [CrossRef]
  60. Zhao, T.; Zhang, D.; Liu, Y.; Feng, M.; Xu, Z.; Huang, H.; Zhang, L.; Li, W.; Li, X. The Association Between GC Gene Polymorphisms and Metabolic Syndrome in Chinese Rural Population: A Case-Control Study. Diabetes Metab. Syndr. Obes. 2022, 15, 165–174. [Google Scholar] [CrossRef]
  61. Ruiz-Ojeda, F.J.; Anguita-Ruiz, A.; Leis, R.; Aguilera, C.M. Genetic Factors and Molecular Mechanisms of Vitamin D and Obesity Relationship. Ann. Nutr. Metab. 2018, 73, 89–99. [Google Scholar] [CrossRef]
  62. Hsieh, M.-C.; Hsiao, P.-J.; Liao, M.-T.; Hou, Y.-C.; Chang, Y.-C.; Chiang, W.-F.; Wu, K.-L.; Chan, J.-S.; Lu, K.-C. The Role of Vitamin D in SARS-CoV-2 Infection and Acute Kidney Injury. Int. J. Mol. Sci. 2022, 23, 7368. [Google Scholar] [CrossRef] [PubMed]
  63. Agrawal, T.; Gupta, G.K.; Agrawal, D.K. Vitamin D Deficiency Decreases the Expression of VDR and Prohibitin in the Lungs of Mice with Allergic Airway Inflammation. Exp. Mol. Pathol. 2012, 93, 74–81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Kere, J. Genetics of Complex Disorders. Biochem. Biophys. Res. Commun. 2010, 396, 143–146. [Google Scholar] [CrossRef] [PubMed]
  65. Buniello, A.; MacArthur, J.A.L.; Cerezo, M.; Harris, L.W.; Hayhurst, J.; Malangone, C.; McMahon, A.; Morales, J.; Mountjoy, E.; Sollis, E.; et al. The NHGRI-EBI GWAS Catalog of Published Genome-Wide Association Studies, Targeted Arrays and Summary Statistics 2019. Nucleic Acids Res. 2019, 47, D1005–D1012. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Kamat, M.A.; Blackshaw, J.A.; Young, R.; Surendran, P.; Burgess, S.; Danesh, J.; Butterworth, A.S.; Staley, J.R. PhenoScanner V2: An Expanded Tool for Searching Human Genotype-Phenotype Associations. Bioinformatics 2019, 35, 4851–4853. [Google Scholar] [CrossRef] [Green Version]
  67. Staley, J.R.; Blackshaw, J.; Kamat, M.A.; Ellis, S.; Surendran, P.; Sun, B.B.; Paul, D.S.; Freitag, D.; Burgess, S.; Danesh, J.; et al. PhenoScanner: A Database of Human Genotype-Phenotype Associations. Bioinformatics 2016, 32, 3207–3209. [Google Scholar] [CrossRef] [Green Version]
  68. Cunningham, F.; Allen, J.E.; Allen, J.; Alvarez-Jarreta, J.; Amode, M.R.; Armean, I.M.; Austine-Orimoloye, O.; Azov, A.G.; Barnes, I.; Bennett, R.; et al. Ensembl 2022. Nucleic Acids Res. 2022, 50, D988–D995. [Google Scholar] [CrossRef]
  69. Howe, K.L.; Achuthan, P.; Allen, J.; Allen, J.; Alvarez-Jarreta, J.; Amode, M.R.; Armean, I.M.; Azov, A.G.; Bennett, R.; Bhai, J.; et al. Ensembl 2021. Nucleic Acids Res. 2021, 49, D884–D891. [Google Scholar] [CrossRef]
  70. Borg, J.; Georgitsi, M.; Aleporou-Marinou, V.; Kollia, P.; Patrinos, G.P. Genetic Recombination as a Major Cause of Mutagenesis in the Human Globin Gene Clusters. Clin. Biochem. 2009, 42, 1839–1850. [Google Scholar] [CrossRef]
  71. Machiela, M.J.; Chanock, S.J. LDlink: A Web-Based Application for Exploring Population-Specific Haplotype Structure and Linking Correlated Alleles of Possible Functional Variants. Bioinformatics 2015, 31, 3555–3557. [Google Scholar] [CrossRef] [Green Version]
  72. R Core Team. R: A Language and Environment for Statistical Computing. 2020. Available online: https://www.scirp.org/%28S%28351jmbntvnsjt1aadkposzje%29%29/reference/referencespapers.aspx?referenceid=2903185 (accessed on 3 November 2022).
  73. Sherman, B.T.; Hao, M.; Qiu, J.; Jiao, X.; Baseler, M.W.; Lane, H.C.; Imamichi, T.; Chang, W. DAVID: A Web Server for Functional Enrichment Analysis and Functional Annotation of Gene Lists (2021 Update). Nucleic Acids Res. 2022, 50, W216–W221. [Google Scholar] [CrossRef] [PubMed]
  74. GTEx Consortium. Genetic Effects on Gene Expression across Human Tissues. Nature 2017, 550, 204–213. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Watanabe, K.; Taskesen, E.; van Bochoven, A.; Posthuma, D. Functional Mapping and Annotation of Genetic Associations with FUMA. Nat. Commun. 2017, 8, 1826. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Supek, F.; Bošnjak, M.; Škunca, N.; Šmuc, T. REVIGO Summarizes and Visualizes Long Lists of Gene Ontology Terms. PLoS ONE 2011, 6, e21800. [Google Scholar] [CrossRef] [Green Version]
  77. Yu, G.; Wang, L.-G.; Han, Y.; He, Q.-Y. ClusterProfiler: An R Package for Comparing Biological Themes Among Gene Clusters. OMICS 2012, 16, 284–287. [Google Scholar] [CrossRef]
  78. Wu, T.; Hu, E.; Xu, S.; Chen, M.; Guo, P.; Dai, Z.; Feng, T.; Zhou, L.; Tang, W.; Zhan, L.; et al. ClusterProfiler 4.0: A Universal Enrichment Tool for Interpreting Omics Data. Innovation 2021, 2, 100141. [Google Scholar] [CrossRef]
  79. Piñero, J.; Ramírez-Anguita, J.M.; Saüch-Pitarch, J.; Ronzano, F.; Centeno, E.; Sanz, F.; Furlong, L.I. The DisGeNET Knowledge Platform for Disease Genomics: 2019 Update. Nucleic Acids Res. 2020, 48, D845–D855. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Piñero, J.; Bravo, À.; Queralt-Rosinach, N.; Gutiérrez-Sacristán, A.; Deu-Pons, J.; Centeno, E.; García-García, J.; Sanz, F.; Furlong, L.I. DisGeNET: A Comprehensive Platform Integrating Information on Human Disease-Associated Genes and Variants. Nucleic Acids Res. 2017, 45, D833–D839. [Google Scholar] [CrossRef]
  81. Piñero, J.; Queralt-Rosinach, N.; Bravo, À.; Deu-Pons, J.; Bauer-Mehren, A.; Baron, M.; Sanz, F.; Furlong, L.I. DisGeNET: A Discovery Platform for the Dynamical Exploration of Human Diseases and Their Genes. Database 2015, 2015, bav028. [Google Scholar] [CrossRef] [Green Version]
  82. Glueck, D.H.; Mandel, J.; Karimpour-Fard, A.; Hunter, L.; Muller, K.E. Exact Calculations of Average Power for the Benjamini-Hochberg Procedure. Int. J. Biostat. 2008, 4, 11. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Heat map of the VD deficiency gene expression in different tissues. The tissue-specific expression patterns are based on GTEx v6 RNA−seq data. Cells in red represent a higher expression, compared to the cells in blue. Gene expression comparisons between tissues (horizontal comparison) within a gene (y axis) are comparable but not those of different genes within a tissue (vertical comparison). Thus, cells in red represent a higher expression of genes in a corresponding tissue, compared to other tissues, but do not represent a higher expression, compared to other genes. Source: FUMA.
Figure 1. Heat map of the VD deficiency gene expression in different tissues. The tissue-specific expression patterns are based on GTEx v6 RNA−seq data. Cells in red represent a higher expression, compared to the cells in blue. Gene expression comparisons between tissues (horizontal comparison) within a gene (y axis) are comparable but not those of different genes within a tissue (vertical comparison). Thus, cells in red represent a higher expression of genes in a corresponding tissue, compared to other tissues, but do not represent a higher expression, compared to other genes. Source: FUMA.
Ijms 23 14256 g001
Figure 2. Enrichment analysis of a VD deficiency, regarding the KEGG pathways (A) and GO in MSigDB on the biological processes (B). Pathways and processes that are overrepresented in the gene set of interest are shown. Furthermore, only genes present in a given pathway or process are given. The overrepresentation of the gene set of interest (input genes) regarding pathways and functions was tested using the hypergeometric test against the pre−defined gene sets (background genes), obtained from MsigDB or KEGG. Input genes that are overlapping in the pathway or process, the enrichment p−value and the proportion of the overlapping genes (input genes relative to the tested gene set) are shown. Source: FUMA.
Figure 2. Enrichment analysis of a VD deficiency, regarding the KEGG pathways (A) and GO in MSigDB on the biological processes (B). Pathways and processes that are overrepresented in the gene set of interest are shown. Furthermore, only genes present in a given pathway or process are given. The overrepresentation of the gene set of interest (input genes) regarding pathways and functions was tested using the hypergeometric test against the pre−defined gene sets (background genes), obtained from MsigDB or KEGG. Input genes that are overlapping in the pathway or process, the enrichment p−value and the proportion of the overlapping genes (input genes relative to the tested gene set) are shown. Source: FUMA.
Ijms 23 14256 g002
Figure 3. Enrichment analysis of COVID-19, regarding the KEGG pathways (A) and signaling pathways in the Reactome database (B). Pathways that are overrepresented in the gene set of interest are shown. Furthermore, only genes present in a given pathway are given. Overrepresentation of the gene set of interest (input genes), regarding the pathways was tested using the hypergeometric test against the pre−defined gene sets (background genes) obtained from Reactome or KEGG. Input genes that are overlapping in the pathway, the enrichment p-value and the proportion of the overlapping genes (input genes relative to the tested gene set) are shown. Source: FUMA.
Figure 3. Enrichment analysis of COVID-19, regarding the KEGG pathways (A) and signaling pathways in the Reactome database (B). Pathways that are overrepresented in the gene set of interest are shown. Furthermore, only genes present in a given pathway are given. Overrepresentation of the gene set of interest (input genes), regarding the pathways was tested using the hypergeometric test against the pre−defined gene sets (background genes) obtained from Reactome or KEGG. Input genes that are overlapping in the pathway, the enrichment p-value and the proportion of the overlapping genes (input genes relative to the tested gene set) are shown. Source: FUMA.
Ijms 23 14256 g003
Figure 4. Gene-disease heat map by groups of diseases for VD deficiency genes (A) and COVID-19 (B). The x-scale shows the list of genes (14 for VD and 29 for COVID-19); the y-axis shows the diseases grouped by the MeSH disease classes. The color scale is proportional to the percentage of diseases in each disease class. Source: DisGeNET.
Figure 4. Gene-disease heat map by groups of diseases for VD deficiency genes (A) and COVID-19 (B). The x-scale shows the list of genes (14 for VD and 29 for COVID-19); the y-axis shows the diseases grouped by the MeSH disease classes. The color scale is proportional to the percentage of diseases in each disease class. Source: DisGeNET.
Ijms 23 14256 g004
Figure 5. Enrichment Analysis of the genes associated with a VD deficiency (A) and COVID-19 (B). Representation of the existing association between the diseases and the enriched genes, based on the ratio (input genes relative/associated) and the FDR. The size of the dots refers to the number of genes enriched for that disease (the greater the number of genes, the larger the dot). The FDR is defined by a color scale, the closer it is to red, the greater the FDR value and the greater the association.
Figure 5. Enrichment Analysis of the genes associated with a VD deficiency (A) and COVID-19 (B). Representation of the existing association between the diseases and the enriched genes, based on the ratio (input genes relative/associated) and the FDR. The size of the dots refers to the number of genes enriched for that disease (the greater the number of genes, the larger the dot). The FDR is defined by a color scale, the closer it is to red, the greater the FDR value and the greater the association.
Ijms 23 14256 g005aIjms 23 14256 g005b
Figure 6. Workflow of the data analysis.
Figure 6. Workflow of the data analysis.
Ijms 23 14256 g006
Table 1. Genes and variants associated with a VD deficiency. Information on the genes to which these variants belong, the risk allele of each SNP, location in the genome (chromosome: chr) for the hg38 coordinates, and the p-value (existing degree of significance between that variant and the disease). Study accession numbers obtained from the GWAS catalog. Population: European (GCST000697, GCST90020244), Asian (GCST90101732).
Table 1. Genes and variants associated with a VD deficiency. Information on the genes to which these variants belong, the risk allele of each SNP, location in the genome (chromosome: chr) for the hg38 coordinates, and the p-value (existing degree of significance between that variant and the disease). Study accession numbers obtained from the GWAS catalog. Population: European (GCST000697, GCST90020244), Asian (GCST90101732).
VariantRisk AlleleGeneLocation
(Chr: Position)
p-ValueStudy
Accession
PubMed ID
rs3750297APADI11:172331813 × 10−10GCST90020244 34308111
rs12123821TFLG-AS11:1522066766 × 10−26
rs4845491C SMCP, LCE6A1:1528770937 × 10−10
rs3755322GUGT1A9, UGT1A5, UGT1A10, UGT1A7, UGT1A8, UGT1A62:2337131416 × 10−12
rs6600893C UGT2B11, UGT2B74:691131835 × 10−15
rs2282679GGC4:717426661 × 10−200
rs2205262CLINC005368:1159996595 × 10−11
rs7129781CCYP2R111:148908714 × 10−33
rs4944958ANADSYN111:714570279 × 10−143
rs964184GZPR111:1167782013 × 10−14
rs10859995THAL12:959819043 × 10−33
rs1800588TLIPC, ALDH1A215:584314762 × 10−10
rs55791371CSMARCA419:110774771 × 10−9
rs10426201GSULT2A119:478814928 × 10−20
rs17217119GBCAS1, CYP24A120:541260515 × 10−16
rs11723621GGC4:717496452 × 10−24GCST90101732 34852423
rs7041GGC4:717526172 × 10−9
rs11023332GPDE3B11:147625643 × 10−11
rs12785878GNADSYN111:714564032 × 10−27GCST000697 20541252
rs10741657GCALCB, CYP2R111:148933323 × 10−20
Table 2. Enrichment analysis of the genes associated with a VD deficiency obtained in disgenet2r (R package). The ratio column represents the number of genes from the input list that are associated with the disease. The BgRatio are the number of genes associated with the disease from the length of the universe, in the example, all genes in DisGeNET (the curated or background genes). The p-values resulting from the multiple Fisher tests are corrected for the false discovery rate (FDR) using the Benjamini–Hochberg (BH) method.
Table 2. Enrichment analysis of the genes associated with a VD deficiency obtained in disgenet2r (R package). The ratio column represents the number of genes from the input list that are associated with the disease. The BgRatio are the number of genes associated with the disease from the length of the universe, in the example, all genes in DisGeNET (the curated or background genes). The p-values resulting from the multiple Fisher tests are corrected for the false discovery rate (FDR) using the Benjamini–Hochberg (BH) method.
DescriptionFDRRatioBgRatio
Vitamin D3 measurement5.029 × 10−115/1414/21,666
Vitamin D measurement5.029 × 10−115/1414/21,666
Rickets2.933 × 10−54/1472/21,666
Metabolic Syndrome X2.933 × 10−58/141125/21,666
Serum albumin measurement3.233 × 10−56/14433/21,666
Glucose measurement3.689 × 10−54/1489/21,666
Elevated blood glucose level3.689 × 10−54/1489/21,666
High density lipoprotein cholesterol level quantitative trait locus 124.847 × 10−52/142/21,666
Coronary heart disease1.732 × 10−4 8/141576/21,666
Vitamin D Deficiency2.267 × 10−4 4/14153/21,666
Table 3. Enrichment analysis of the genes associated with COVID-19, obtained in disgenet2r (R package). The ratio column represents the number of genes from the input list that are associated with the disease. The BgRatio are the number of genes associated with the disease from the length of the universe, in the example, all genes in DisGeNET (the curated or background genes). The p-values resulting from the multiple Fisher tests are corrected for the false discovery rate (FDR) using the Benjamini–Hochberg (BH) method.
Table 3. Enrichment analysis of the genes associated with COVID-19, obtained in disgenet2r (R package). The ratio column represents the number of genes from the input list that are associated with the disease. The BgRatio are the number of genes associated with the disease from the length of the universe, in the example, all genes in DisGeNET (the curated or background genes). The p-values resulting from the multiple Fisher tests are corrected for the false discovery rate (FDR) using the Benjamini–Hochberg (BH) method.
DescriptionFDRRatioBgRatio
Primary biliary cirrhosis3.305 × 10−4 8/29478/21,666
Coughing4.286 × 10−4 6/29235/21,666
Reticulocyte count (procedure)4.286 × 10−4 6/29234/21,666
Eosinophilia-Myalgia Syndrome, L-Tryptophan-Related9.429 × 10−4 2/292/21,666
Liver diseases1.163 × 10−3 9/291019/21,666
Rales1.163 × 10−3 3/2923/21,666
Sarcoidosis, Pulmonary1.163 × 10−3 4/2981/21,666
Inclusion Body Myositis (disorder)1.163 × 10−3 4/2987/21,666
Juvenile Graves disease1.163 × 10−3 2/293/21,666
Podoconiosis1.163 × 10−3 2/293/21,666
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alcalá-Santiago, Á.; Rodríguez-Barranco, M.; Rava, M.; Jiménez-Sousa, M.Á.; Gil, Á.; Sánchez, M.J.; Molina-Montes, E. Vitamin D Deficiency and COVID-19: A Biological Database Study on Pathways and Gene-Disease Associations. Int. J. Mol. Sci. 2022, 23, 14256. https://doi.org/10.3390/ijms232214256

AMA Style

Alcalá-Santiago Á, Rodríguez-Barranco M, Rava M, Jiménez-Sousa MÁ, Gil Á, Sánchez MJ, Molina-Montes E. Vitamin D Deficiency and COVID-19: A Biological Database Study on Pathways and Gene-Disease Associations. International Journal of Molecular Sciences. 2022; 23(22):14256. https://doi.org/10.3390/ijms232214256

Chicago/Turabian Style

Alcalá-Santiago, Ángela, Miguel Rodríguez-Barranco, Marta Rava, María Ángeles Jiménez-Sousa, Ángel Gil, María José Sánchez, and Esther Molina-Montes. 2022. "Vitamin D Deficiency and COVID-19: A Biological Database Study on Pathways and Gene-Disease Associations" International Journal of Molecular Sciences 23, no. 22: 14256. https://doi.org/10.3390/ijms232214256

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop