1. Introduction
Cardiovascular diseases (CVDs) are a group of disorders of the heart and blood vessels. They are classified into coronary heart disease, cerebrovascular disease, peripheral arteriopathies, rheumatic heart disease, congenital heart disease, deep vein thrombosis, and pulmonary embolism. According to the World Health Organization, CVD is the leading cause of death worldwide. It is estimated that in 2015 (the last year for which data have been published) 17.7 million people died from CVDs, representing 31% of all registered global deaths. Of these deaths, more than 80% take place in low- and middle-income countries, affecting men and women almost equally [1]; hence the utmost importance of identifying the risk factors involved in the development of this disease, given the enormous social and economic implications [2].
Although the etiology of CVD has not been clearly established, it has been shown that the causes of these disorders are multi-factorial, due to the combination of environmental risk factors, such as smoking, lack of physical activity, eating habits, high blood pressure, type 2 diabetes, and dyslipidemias, with genetic predisposition [3,4,5]. The search for genes that predispose to CVD has led to the identification of human variations of deoxyribonucleic acid (DNA), evaluation of the risk profile and adoption of preventive or therapeutic measures [4]. The vitamin D receptor (VDR) gene and its single nucleotide polymorphisms (SNPs) have received special attention due to their association with a cardiometabolic risk profile. The precise mechanism underlying its influence on pathogenesis is still unclear and may be the result of different factors [6,7,8]. Firstly, the VDR gene is found in vascular smooth muscle cells and endothelial cells, potentially affecting their growth and proliferation [9]. In addition, VDR activation induces an increase dependent on nitric oxide concentration in endothelial cells, and enhances the angiogenic properties of endothelial progenitor cells [10,11]. Secondly, vitamin D could regulate immune cells, inhibiting the release of pro-inflammatory cytokines and increasing the release of anti-inflammatory cytokines, thus playing a role in blood vessel protection [12]. Thirdly, vitamin D is an important regulator of the renin-angiotensin-system (RAS). Vitamin D deficiency would involve activation of the renin gene, producing an increase in angiotensin II, which can lead to hypertension and ventricular hypertrophy [13,14]. Additionally, it would prompt an increase in the production of reactive oxygen species (ROS) and the activation of G proteins such as Rho A, resulting in the inhibition of the pathways necessary for intracellular glucose transport and, therefore, the development of insulin resistance and the onset of metabolic syndrome [15]. Finally, vitamin D is associated with an atherogenic lipid profile that includes increased serum LDL and decreased HDL levels. The biologically active form of vitamin D suppresses foam cell formation, decreases cholesterol uptake by macrophages, and induces LDL autophagy through gene regulation through VDR. However, short-term replenishment of 25-hydroxyvitamin D levels does not improve the lipid profile in humans [16,17]. In short, vitamin D can lower blood pressure values and have anti-inflammatory, anti-proliferative, anti-hypertrophic, anti-fibrotic, anti-diabetic, and anti-thrombotic effects, beneficially modulating classic cardiovascular risk factors [18].
Vitamin D can be obtained from 7-dehydrocholesterol after exposure to sunlight and through food [19,20,21]. Vitamin D is hydroxylated to 25-hydroxyvitamin D3 (25(OH)D3) in the liver and 1α-hydroxylated in the kidney, forming the active hormone 1,25-dihydroxyvitamin D3 (1.25(OH)2D3) [22,23]. This active metabolite binds to the VDR ligand-binding domain (LBD), forming a heterodimer with the retinoid X receptor (RXR) that binds to vitamin D response elements (VDRE) in the promoter region of target genes modulating transcriptional activation [24,25]. The VDR is expressed almost ubiquitously regulating approximately 3% of the genome; more than 900 genes participate in many physiological processes [26,27].
The gene encoding the nuclear vitamin D receptor is a large gene more than 100 Kb in length [28], found in the long arm of chromosome 12 (locus 12q13.1) [29]. It has at least five promoter regions, eight coding exons and at least six non-coding exons that are alternately spliced [30]. It encodes a protein that contains 427 amino acids [26] belonging to the family of steroid receptors for retinoic acid, thyroid hormone, sex hormones, and adrenal steroids [31]. To date, 470 SNPs have been identified at the VDR locus [32], but there are five common polymorphisms that have been widely studied due to their effects on various physiological and pathological phenotypes; Cdx2 (rs11568820), FokI (rs2228570), TaqI (rs731236), BsmI (rs1544410) and ApaI (rs7975232) [33,34].
Recently, the importance of these polymorphisms and their haplotypes has been increasingly recognized as more studies have linked them to different diseases [35]. Some SNPs play a key role in modifying 1.25(OH)2D3 uptake, due to their capacity to modify vitamin D function [36]. However, the exact molecular mechanism explaining the association between VDR polymorphisms and serum levels of 25(OH)2D3 remains unknown [37]. The VDR SNP FokI (rs2228570, exon 2, C > T, formerly known as rs10735810) is characterized by the presence of two ATG start codons separated by six nucleotides, modifying the length and functional activity of the protein. Thus, FokI is the only VDR polymorphism with functional impact, as it involves the loss of a transcription start site. It is not in linkage disequilibrium with other SNPs, so the associations with the VDR genotype FokI are considered independent markers of the VDR gene [29]. The C allele results in a shorter variant of the VDR protein (424 amino acids) and a longer variant associated with the T allele (427 amino acids) [38]. Although no significant differences have been reported in ligand affinity, DNA binding or transactivation activity between the two allelic forms of FokI, the shorter variant shows greater receptor activity than the longer variant as it appears to interact more efficiently with the transcription factor TFIIB [39,40]. The SNPs BsmI (rs1544410, intron 8, G > A), TaqI (rs731236, exon 9, T > C) and ApaI (rs7975232, intron 8, C > A) are found near the 3’UTR region of the VDR gene [6,27,35]. This region contains many polymorphisms and therefore a strong linkage disequilibrium could explain the associations observed with Bsm-Apa-Taq. These polymorphisms are involved in the regulation of VDR expression, notably through regulation of messenger ribonucleic acid (mRNA) stability [41,42]. Therefore, they can influence the expression of the protein without altering its structure or function [43]. The polymorphism Cdx2 (rs11568820, exon 1, G > A) is located in the promoter region of the 5’ end of the VDR gene. The VDR specifically interacts with functional enhancer elements in the VDR gene in the small intestine, regulating its gene expression and therefore calcium absorption [44,45]. The transcriptional activity of the VDR promoter with the G allele is lower compared to the A allele, decreasing intestinal absorption of calcium [44] and potentially influencing central obesity [27].
Numerous studies have been carried out under this conceptual framework to investigate the association of these polymorphisms with the risk of developing CVD. However, the results obtained are contradictory and further studies are required in different populations to obtain more information about the influence of these SNPs on the susceptibility of developing CVD. Based on the foregoing, we carried out this study to evaluate the association of these polymorphisms in the VDR gene (BsmI (rs1544410), TaqI (rs731236), ApaI (rs7975232), FokI (rs2228570) and Cdx2 (rs11568820)) and CVD risk in a Caucasian population.
2. Results
2.1. Patients Characteristics
A total of 246 cases with CVD and 246 controls of Caucasian origin from Southern Spain were included in the study. Controls were matched to cases on age and sex (1:1). Their clinical, socio-demographic, and pathological characteristics are detailed in Table 1. The median age of the patients was 70 (63.25, 76) years. The group of cases consisted of 132 women (53.7%) and 114 men (46.3%). Most of the cases were non-smokers, non-drinkers, overweight, not dyslipidemic, hypertensive and did not have diabetes. As regards cardiovascular disease, 84 (34.14%) suffered from cardiac arrhythmias, 15 (6.10%) presented cardiomyopathy, 15 (6.10%) cerebrovascular disease, 27 (10.98%) heart failure, 46 (18.70%) valvular heart disease and 59 (23.98%) vascular peripheral disease.
The control group consisted of 142 women (142/246; 57.7%) and 104 men (104/246; 42.3%), with a mean age of 69.5 (63, 76) years. Most of the controls were non-smokers, non-drinkers, had a healthy weight, not dyslipidemic, did not have hypertensive and did not have diabetes.
Cases and controls showed significant differences in terms of BMI (p = 0.0003) and diabetes (p < 0.001; OR = 10.53; 95% CI = 1.32–3.03; diabetes vs. non-diabetes).
2.2. Genotype Distribution
The genotype frequency for the control group coincided with the expected values according to the Hardy-Weinberg equilibrium model (Table S1). The coefficients D’ and r2 for the evaluation of linkage disequilibrium are described in Table S2 and Figure 1. All SNPs had minor allele frequencies above 1%, so none of them had to be excluded from the analysis (Table S3). The estimated haplotype frequencies are presented in Table S4.
2.3. Influence of Genetic Polymorphisms on CVD Risk
Bivariate analysis was performed according to the following models: genotypic, additive, allelic, dominant, and recessive. VDR FokI (rs2228570) was the only polymorphism to show an association with CVD risk in the genotypic (pχ² = 0.0063), additive (pχ² = 0.0267), allelic (pχ² = 0.0230) and recessive (pχ² = 0.0015) models (Table S5). This significant association was also found after making adjustments for the Bonferroni test, in which both the genotypic and recessive models showed that patients with the TT genotype presented a higher risk of developing CVD (pBonferroni-corrected = 0.0305; OR = 2.35; 95% CI = 1.33–4.28; TT vs. CC and pBonferroni-corrected =0.0075; OR = 2.39; 95% CI = 1.39–4.24; TT vs. C; Table 2, respectively). The logistic regression analysis adjusted for BMI and diabetes revealed that the TT genotype was associated with a higher risk of CVD in both the genotypic model (p = 0.0430; OR = 2.30; 95% CI = 1.06–5.37; TT vs. CC; Table 3) and the recessive model (p = 0.0099; OR = 2.71; 95% CI = 1.31–6.07; TT vs. C; Table 3). However, no association was found in either the genotypic logistic regression model adjusted for BMI and diabetes (p = 0.2286; OR = 0.72; 95% CI = 0.43–1.22; TC vs. CC; Table 3), the dominant model (p = 0.8648 OR = 0.95; 95% CI = 0.58–157; T vs. CC; Table 3) or the additive model (p = 0.2151; OR = 1.24; 95% CI = 0.88–1.76; T vs. C; Table 3). The other polymorphisms analyzed did not show statistically significant associations with the development of CVD in any of the models studied (Table S5). Haplotype analysis revealed that the GAC haplotype (p = 0.047; OR = 0.34; 95% CI = 0.12–0.98) was associated with an increased risk of CVD (Table 4).
3. Discussion
Cardiovascular disease, which involves the heart, brain, and peripheral circulation, represents a major health problem throughout the world. It has multiple genetic and environmental components that contribute to the observed phenotype [46]. The role played by vitamin D in cardiovascular function is receiving more and more attention and, in particular, the effect of the respective polymorphisms in the VDR gene. Most studies describing VDR SNPs have focused on five polymorphisms: BsmI (rs1544410), TaqI (rs731236), ApaI (rs7975232), FokI (rs2228570) and Cdx2 (rs11568820) [32,33].
The VDR FokI polymorphism (rs2228750) is one of the most researched. It is the only VDR polymorphism with functional impact, as it involves the loss of a transcription start site [29]. Previous studies have reported contradictory results, so the relationship between this polymorphism and the development of CVD is not clear [47,48,49,50,51,52,53,54]. In our study, the logistic regression analysis adjusted for body mass index and diabetes revealed that the VDR FokI TT genotype showed an association with increased risk of CVD in the genotypic model (p < 0.05; Table 3) and in the recessive model (p < 0.05; Table 3). Several studies are consistent with our results, for example one study conducted in a Caucasian population (from Poland) (58 cases/142 controls), in which statistically significant differences were found in relation to the T allele of the VDR FokI polymorphism, identifying it as a risk factor for heart failure and hypertension (p = 0.03) [55]. Another study in a Caucasian population (from Egypt) (50 cases/50 controls) also reported an association between the VDR FokI polymorphism and the presence of congenital heart disease defects. In particular, patients carrying the VDR FokI-T allele were at greater risk of suffering congenital heart disease defects [(p = 0.006; OR = 3; 95% CI = 1–8; TC vs. CC); (p = 0.03; OR = 10; 95% CI = 1–97; TT vs. CC); (p = 0.001; OR = 3; 95% CI = 2–6; T vs. C)] [50]. Another study in a Caucasian population (from Egypt) (244 cases/138 controls) shows that allele T frequency is higher in patients with CVD (group A and B) which means VDR FokI-T allele increase CVD risk [(p < 0.001; OR = 2.397; 95% CI = 1.518–3.786; C vs. T); (p < 0.001; OR = 2.298; 95% CI = 1.443–3.660; C vs. T) [56]. Furthermore, the study carried out by Nakhl et al. (2019) in a Caucasian population (Mediterranean region) with 50 patients reported that the presence of VDR FokI polymorphisms is higher compared to the rest of the Caucasian population [46]. Moreover, the results of our study are in accord with others conducted in an Asian population. For example, the study of Hao et al. (2019) in a Chinese population (145 cases/90 controls), in which statistically significant differences were found between cases and controls in relation to the T allele of the VDR FokI polymorphism, identifying it as a risk factor for the disease [(p < 0.0001; OR = 2.45; 95% CI = 1.60–3.77; T vs. C); (p < 0.0001; OR = 3.42; 95% CI = 1.94–6.07; TC vs. CC); (p = 0.02; OR = 4.73; 95% CI = 1.35–16.27; TT vs. CC)] [47]. Likewise, another study carried out in an Asian population (from China) (860 cases/862 controls) showed that individuals carrying the VDR FokI-T allele were at a higher risk of developing coronary artery disease (CAD) [(p < 0.001; OR = 1.47; 95% CI = 1.22–1.85; CT vs. CC); (p < 0.001; OR = 2.12; 95% CI = 1.43–2.88; TT vs. CC); (p < 0.001; OR = 1.68; 95% CI = 1.26–2.17; TC + TT vs. CC)] [48]. Similarly, a study carried out in an Asian population (from India) (292 cases/219 controls) reported a statistically significant association between the T genotype of the VDR FokI polymorphism and susceptibility to developing ischemic cerebrovascular accident (CVA) (p = 0.01; OR = 1.89; 95% CI = 1.17–3.07; TT vs. CC) [49]. Furthermore, according to the meta-analysis conducted by Tabaei et al. (2021) in Caucasian and Asian populations (3741 cases/1998 controls), the T allele of the VDR FokI SNP was significantly associated with the CAD risk in all genetic models including; dominant model (p = 0.02; OR = 1.47, 95% CI 1.06–2.03), recessive model (p < 0.001; OR = 1.39, 95% CI 1.14–1.71), allelic model (p = 0.01; OR = 1.38, 95% CI 1.08–1.76, p = 0.01), TT vs. CC model (p = 0.02; OR = 1.58, 95% CI 1.07–2.34), and TC vs. TT model (p = 0.04; OR = 1.38, 95% CI 1.00–1.92), having a strong significant association for the Asian population but not the Caucasian [57].
However, a study carried out in a Caucasian population (from Turkey) (54 cases/58 controls) evaluated the influence of the VDR FokI polymorphism on the risk of arteriovenous fistula (AVF), but found no statistically significant association (p = 0.168; OR = 1.683; 95% CI = 0.803–3.531) [52]. Another study carried out in a Caucasian population (United States) (205 cases/206 controls) also evaluated the influence of the VDR FokI polymorphism on the risk of congestive heart failure (HF) in patients with hypertension, without finding any statistically significant association in any of the analyzed models [(p > 0.05 OR = 0.82; 95% CI = 0.47–1.42; TC vs. CC); (p > 0.05; OR = 0.85; 95% CI = 0.39–1.88; TT vs. CC); (p > 0.05; OR = 0.82; 95% CI = 0.49–1.39; CT + TT vs. CC)] [53]. Finally, a meta-analysis composed of six studies in Caucasian and Asian populations (2725 cases/1017 controls) evaluated the influence of the VDR FokI polymorphism on the risk of CAD, but did not observe any statistically significant association in any of the models analyzed [(p > 0.05; OR = 0.93; 95% CI = 0.76–1.10; CT + TT vs. CC); (p> 0.05; OR = 1.14; 95% CI = 0.86–1.43; TT vs. CT + CC); (p > 0.05; OR = 1.00; 95% CI = 0.88–1.12; T vs. C); (p > 0.05; OR = 0.88; 95% CI = 0.71–1.06; CT vs. CC); (p > 0.05; OR = 1.06; 95% CI = 0.77–1.35; TT vs. CC)] [54].
According to the results shown in the majority of the studies, including ours, the VDR FokI SNP presents statistical association with CVD in Asian and Caucasian (from Europe) populations. However, this statistical significance does not occur in US Caucasian populations. It seems to indicate that ethnicity may influences this correlation.
Otherwise, BsmI (rs1544410), TaqI (rs731236) and ApaI (rs7975232) do not result in a structural change of the VDR protein. However, these polymorphisms can modify the expression of the VDR gene and increase susceptibility to developing CVD [29,43]. In our study, the VDR BsmI polymorphism did not show an association with susceptibility to developing CVD in any of the models analyzed (Table 2). Various studies in Caucasian populations are consistent with our results, such as the studies conducted by Ewida et al. (2021) [58] and Nakhl et al. (2019), in which there were no association between those polymorphisms and CVD [46]. However, according to the results of the study carried out by Abouzid et al. (2021) they suggest that VDR BsmI-AA genotype is protective of CVD (p = 0.04) [53]. Moreover, according to the meta-analysis by Tabaei et al. (2021) in Caucasian and Asian populations (6169 cases/2392 controls), VDR ApaI-CC genotype was significantly associated with CAD risk in all genetic models, including: dominant model (p = 0.002; OR = 1.20, 95% CI 1.07–1.34), allelic model (p = 0.004; OR = 1.12, 95% CI 1.04–1.21), CC vs. AA model (p = 0.01; OR = 1.22, 95% CI 1.04–1.44), and AC vs. AA model (p = 0.007; OR = 1.18, 95% CI 1.05–1.33), other than the recessive model (p = 0.17; OR = 1.10, 95% CI 0.96–1.27) [57].
Finally, the Cdx2 polymorphism (rs11568820) is located in the promoter region of the 5′ end of the VDR gene regulating gene expression [44]. So far, no study has evaluated the impact of the VDR Cdx2 polymorphism on the susceptibility for developing CVD. In our study, no significant association was found between the VDR Cdx2 polymorphism (rs11568820) and the risk of developing CVD in any of the models analyzed (Table 2).
The main limitation of this study was the limited sample size compared to other studies, especially in terms of cases, which could have prevented the detection of certain associations. Another limitation of our study is that serum levels of 25(OH)D are not available as it is a retrospective study. Thus, it has not been possible to correlate the levels of active vitamin D and the genotypes. However, despite the limited sample size, after adjusting for the Bonferroni test to avoid false-positive associations, the effect of VDR FokI remained. The strengths of our study were the homogeneity of the sample, especially in terms of the cases, which consisted of only University Hospital Virgen de las Nieves patients diagnosed by the same team, and also from the same geographical area, thus increasing uniformity. In addition, the controls chosen were older than the cases to reduce potential selection bias.
4. Materials and Methods
4.1. Study Design
A retrospective case-control study was carried out.
4.2. Study Subjects
This included 246 patients with CVD and 246 controls of Caucasian origin from Southern Spain, with a case/control ratio 1. The cases were recruited at the University Hospital Virgen de las Nieves, Granada, Spain, from March 2013 to October 2022. Controls were individuals over 18 years of age who had resided in the same geographic area, with no personal history of CVD, recruited at the same hospital. This sample size has been chosen based on previous studies and the prevalence of these gene polymorphisms in a Spanish Caucasian population.
This case-control study was approved by the Ethics Committee of the Andalusian Health System (SAS) and carried out in accordance with the Declaration of Helsinki (code: 0957-N-21). All subjects participating in the study signed a written informed consent for the extraction of saliva or blood samples and their donation to the biobank. The samples were coded and treated confidentially.
4.3. Socio-Demographic and Clinical Variables
The socio-demographic data collected and included in the study were sex, age, smoking, alcoholic habit, body mass index (BMI), dyslipidemia, hypertension, diabetes, and cardiovascular disease suffered by the patients. Individuals were classified as non-smokers if they had never smoked or smoked < 100 cigarettes in their life, as ex-smokers if they had smoked ≥ 100 cigarettes in their life but were not currently smoking, and as active smokers if they had smoked ≥ 100 cigarettes in their life and were currently smoking. The individuals were classified according to standard drinking units (SDU) as non-drinkers if they were abstainers or did not consume alcohol regularly, as active drinkers if their alcohol consumption >4 SDU/day in men and >2.5 SDU/day in women, and as ex-drinkers if their alcohol consumption was >4 SDU units/day in men and >2.5 SDU/day in women, but they were not currently drinking [59]. For BMI, following Spanish Society for the Study of Obesity criteria, individuals were classified as underweight (BMI < 18.5), healthy weight (18.5 < BMI < 24.9), overweight (25 < BMI < 29.9), class I obese (30 < BMI < 34.9), class II obese (35 < BMI < 39.9) and class III obese (BMI > 40) [60]. The following criteria are based on the guidelines created by the working group of the European Society of Cardiology and the European Society of Atherosclerosis. For dyslipidemia, individuals were classified as dyslipidemic [triglycerides > 150 mg/dL; high-density lipoprotein (HDL) < 40 mg/dL in men and HDL < 48 mg/dL in women; low-density lipoprotein (LDL) > 100 mg/dL; total cholesterol > 200 mg/dL)] and not dyslipidemic. For hypertension, individuals were classified as hypertensive [systolic blood pressure (SBP) > 140 mmHg and diastolic blood pressure (DBP) > 90 mmHg] and non-hypertensive [61]. For diabetes, patients were classified according to American Diabetes Association Criteria into diabetics [steroid diabetes (glucocorticoid-induced hyperglycemia), type I diabetes (fasting blood glucose > 126 mg/dL produced by an autoimmune reaction against insulin-producing cells) and type II diabetes (fasting blood glucose > 126 mg/dL caused by insulin resistance)] and non-diabetics [58]. Both patients with type I and II diabetes have been grouped in the group of cases and controls with diabetes, without excluding any patient. As regards cardiovascular disease, individuals were classified into those with arterial thromboembolism, atrial fibrillation, heart failure, valvular heart disease or venous thromboembolism.
4.4. Genetic Variables
4.4.1. DNA Isolation
The DNA samples, isolated from blood or saliva, were obtained from the University Hospital Virgen de las Nieves Biobank, which belongs to the Biobank of the Andalusian Public Health Systems. Blood samples were collected in BD Vacutainer® tubes with EDTA K3 as anticoagulant (3 mL). Saliva samples were collected in 50-mL BD FalconTM conical tubes (BD, Plymouth, UK). DNA extraction was performed using the DNA QIAamp DNA Mini extraction kit (Qiagen GmbH, Hilden, Germany), in accordance with the specifications provided by the manufacturer, from DNA purification from blood or saliva, and stored at −40°C. DNA concentration and purity were measured using a UV NanoDrop 2000 spectrophotometerTM with the absorbance ratio at 280/260 and 280/230.
4.4.2. Genotyping and Quality Control
The genetic polymorphisms BsmI (rs1544410), TaqI (rs731236), ApaI (rs7975232), FokI (rs2228570) and Cdx2 (rs11568820) were determined by means of real-time polymerase chain reaction (PCR) for allelic discrimination using TaqMan® probes (ABI Applied Biosystems, 7300 Real-Time PCR System), in accordance with the manufacturer’s instructions. The ID assays used for the polymorphisms were as follows: C___2404008_10 for TaqI (rs731236), C__28977635_10 for ApaI (rs7975232), C__12060045_20 for FokI (rs2228570) and C___2880808_10 for Cdx2 (rs11568820). The BsmI polymorphism (rs1544410) was analyzed using an assay customized by ThermoFisher Scientific (Waltham, MA, USA) coded as AN324M4. The criteria for SNPs quality control were: (1) missing genotype rate per SNP < 0.05; (2) minor allele frequency > 0.01; (3) p value > 0.05 in Hardy-Weinberg equilibrium test; (4) missing genotype rate between cases and control < 0.05.
4.4.3. Statistical Analysis
Cases and controls were matched by age and gender with 1:1 propensity score matching method. Quantitative data were expressed as the results (±standard deviation) for variables with normal distribution and medians or percentiles (25 and 75) for variables with non-normal distribution. The Shapiro-Wilks test was used to verify normality.
Hardy-Weinberg equilibrium and pairwise haplotype frequencies were estimated, and Lewontin’s D prime (D′) and the linkage disequilibrium coefficient (r2) were calculated. The bivariate association analysis between CVD risk and polymorphisms was performed with multiple models (genotypic, additive, allelic, dominant, and recessive) using Pearson’s Chi-square test and Fisher’s exact test, providing Odds ratio (OR) values and the corresponding 95% confidence interval (95% CI). The models were defined as follows: allelic (D vs. d), dominant ((DD, Dd) vs. dd), recessive (DD vs. (Dd, dd), genotypic (DD vs. dd, and Dd vs. dd) and additive, with D being the minor allele and d, the major allele. The Bonferroni correction was used for multiple comparisons. Unconditional multiple logistic regression models (genotypic, dominant, and recessive) were considered to determine the influence of possible confounding variables on CVD risk. All tests were two-tailed with a significant level of p < 0.05 and were performed using free access software for the PLINK whole genome association analysis toolset and the statistical program R 3.2.2 [62,63]. Linkage disequilibrium was performed with Haploview 4.2 software and haplotype analysis with SNPStats, a web tool for the analysis of association studies [64,65].
5. Conclusions
The VDR FokI polymorphism (rs2228570) was significantly associated with the development of CVD. According to it, this SNP could be used as a risk biomarker for the disease mentioned. Otherwise, no influence was found of the polymorphisms VDR BsmI (rs1544410), TaqI (rs731236), ApaI (rs7975232) and Cdx2 (rs11568820) on the risk of developing CVD in our patients. Further studies should be done in order to understand and to confirm the role of VDR polymorphisms in the CVD.
P.G.R. and C.P.R. designed and conceived the experiments; P.G.R. performed the experiments; P.G.R. and J.M.G.N. collected samples; P.G.R. and L.E.P.L. analyzed the data; J.M.G.N. and S.R.T. contributed materials and analysis tools; P.G.R. prepared the original draft; C.P.R., J.M.G.N., L.E.P.L. and M.d.C.R.T. reviewed and edited the draft; P.G.R., C.P.R., and S.R.T. reviewed the analysis and interpretation; A.J.M. supervised funding acquisition. Every author participated in reviewing the manuscript and improving the intellectual content. All authors have read and agreed to the published version of the manuscript.
The study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics and Research Committee of the Sistema Andaluz de Salud (code: 0957-N-21).
Informed consent was obtained from all subjects involved in the study.
The results of this study are part of the doctoral thesis presented by Paula González Rojo at the University of Granada.
The authors declare no conflict of interest.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Enlarge this image.Figure 1. Chromosomal location for (vitamin D receptor) VDR, and linkage disequilibrium (LD).
Clinico-pathologic characteristics of cardiology cases and controls.
| Cases | Controls | χ² | p-Value | Reference | OR | 95% CI | |||
|---|---|---|---|---|---|---|---|---|---|
| N | n (%) | N | n (%) | ||||||
| Gender | 246 | 246 | |||||||
| Male | 114 (46.3) | 104 (42.3) | 0.8237 | 0.3641 | |||||
| Female | 132 (53.7) | 142 (57.7) | |||||||
| Age | 246 | 70 (63.25, 76) | 246 | 69.5 (63, 76) | |||||
| Tobacco consumption | 233 | 203 | |||||||
| Current-smokers | 41 (17.6) | 27 (13.3) | 1.9502 | 0.3772 | |||||
| Former-smokers | 70 (30.0) | 70 (34.5) | |||||||
| Non-smokers | 122 (52.4) | 106 (52.2) | |||||||
| Alcohol consumption | 220 | 173 | |||||||
| Current-drinkers | 46 (20.9) | 27 (15.6) | 2.0823 | 0.3531 | |||||
| Former-drinkers | 8 (3.6) | 5 (2.9) | |||||||
| Non-drinkers | 166 (75.5) | 141 (81.5) | |||||||
| Body Mass Index | 182 | 29.63 ± 5.14 | 114 | 28.03 ± 4.89 | |||||
| Underweight | 1 (0.5) | 1 (0.9) | 0.0066 a | Normal (healthy weight) | 1.21 | 0.05–31.44 | |||
| Normal (healthy weight) | 33 (18.1) | 40 (35.1) | 1 | ||||||
| Overweight | 68 (37.4) | 30 (26.3) | 2.74 | 1.47–5.20 | |||||
| Obese Class I (Moderately obese) | 46 (25.3) | 32 (28.1) | 1.74 | 0.92–3.34 | |||||
| Obese Class II (Severely obese) | 29 (15.9) | 10 (8.8) | 3.51 | 1.53–8.57 | |||||
| Obese Class III (Very severely obese) | 5 (2.7) | 1 (0.9) | 6.06 | 0.92–119.19 | |||||
| Dyslipidemia | 246 | 246 | |||||||
| No | 158 (64.2) | 155 (63.0) | 0.0790 | 0.7786 | |||||
| Yes | 88 (35.8) | 91 (37.0) | |||||||
| Hypertension | 246 | 246 | |||||||
| No | 100 (40.6) | 115 (46.8) | 1.8588 | 0.1728 | |||||
| Yes | 146 (48) | 131 (53.2) | |||||||
| Diabetes | 246 | 246 | |||||||
| No | 170 (69.1) | 201 (81.7) | 10.5320 | 0.0012 a | No | 1 | |||
| Yes | 76 (30.9) | 45 (18.3) | 1.99 | 1.32–3.06 | |||||
| Disease | 246 | ||||||||
| Cardiac Arrhythmias | 84 (34.15) | ||||||||
| Cardiomyopathy | 15 (6.10) | ||||||||
| Cerebrovascular disease | 15 (6.10) | ||||||||
| Heart failure | 27 (10.98) | ||||||||
| Heart valve disease | 46 (18.70) | ||||||||
| Peripheral vascular disease | 59 (23.98) | ||||||||
Shade means the value is significant. a p-value for t test.
Influence of VDR FokI (rs2228570) gene polymorphism on risk of cardiovascular disease.
| Models | Genotype | Cases [n (%)] | Controls [n (%)] | p-Value a | Adjusted |
OR c | 95% CI |
|---|---|---|---|---|---|---|---|
| Genotypic | TT | 45 (18.3) | 21 (8.5) | 0.0061 | 0.0305 | 2.35 | 1.33–4.28 |
| CT | 99 (40.2) | 113 (45.9) | 0.96 | 0.66–1.41 | |||
| CC | 102 (41.5) | 112 (45.5) | 1 | ||||
| Dominant | T | 144 (58.5) | 134 (54.5) | 0.3631 | 1 | ||
| CC | 102 (41.5) | 112 (45.5) | |||||
| Recessive | TT | 45 (18.3) | 21 (8.5) | 0.0022 | 0.0075 | 2.39 | 1.39–4.24 |
| C | 201 (81.7) | 225 (91.5) | |||||
| Allelic | T | 189 (38.41) | 155 (31.50) | 0.0230 | 0.1151 | ||
| C | 303 (61.59) | 337 (68.50) | |||||
| Additive | - | - | - | 0.02673 d | 0.1336 |
CI: confidence interval; OR: odds ratio. a p-value for χ2-test. b p-value for Bonferroni correction. c Unadjusted or crude ORs.d p-value for t test. Shade means the value is significant.
Influence of clinical characteristic and VDR FokI (rs2228570) gene polymorphism on risk of cardiovascular disease.
| Genotypic | Dominant | Recessive | Additive | ||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| TT vs. CC | CT vs. CC | T vs. CC | TT vs. C | T vs. C | |||||||||||
| p-Value | OR | 95% CI | p-Value | OR | 95% CI | p-Value | OR | 95% CI | p-Value | OR | 95% CI | p-Value | OR | 95% CI | |
| Body Mass Index | 0.0265 | 1.06 | 1.01–1.12 | 0.0265 | 1.06 | 1.01–1.12 | 0.0458 | 1.05 | 1.00–1.11 | 0.0413 | 1.05 | 1.00–1.11 | 0.0634 | 1.05 | 0.99–1.10 |
| Diabetes | |||||||||||||||
| Yes | 0.0359 | 1.91 | 1.06–3.55 | 0.0359 | 1.91 | 1.06–3.55 | 0.0468 | 1.83 | 1.02–3.38 | 0.0316 | 1.94 | 1.07–3.59 | 0.0372 | 1.89 | 1.05–3.50 |
| VDR FokI (rs2228570) | 0.0430 | 2.30 | 1.06–5.37 | 0.2286 | 0.72 | 0.43–1.22 | 0.8648 | 0.95 | 0.58–1.57 | 0.0099 | 2.71 | 1.31–6.07 | 0.2151 | 1.24 | 0.88–1.76 |
Shade means the value is significant.
Haplotype association with response.
| Rs1544410 | Rs7975232 | Rs731236 | Freq | OR (95% CI) | p-Value | |
|---|---|---|---|---|---|---|
| 1 | G | C | T | 0.4613 | 1.00 | --- |
| 2 | A | A | C | 0.3685 | 0.92 (0.70–1.22) | 0.570 |
| 3 | G | A | T | 0.1079 | 0.81 (0.52–1.24) | 0.330 |
| 4 | A | A | T | 0.0251 | 1.11 (0.51–2.38) | 0.800 |
| 5 | G | A | C | 0.0189 | 0.34 (0.12–0.98) | 0.047 |
| 6 | A | C | T | 0.0114 | 1.22 (0.41–3.65) | 0.730 |
| Rare | * | * | * | 0.0069 | 0.54 (0.12–2.49) | 0.430 |
Freq: haplotype frequency. Shade means the value is significant. * reference to rare haplotypes as there is no symbol to identify the group.
Supplementary Materials
The following are available online at
References
1. WHO. Cardiovascular Diseases; World Health Organization: Ginebra, Switzerland, 2017; Available online: https://www.who.int/news-room/fact-sheets/detail/cardiovascular-diseases-(cvds) (accessed on 3 June 2022).
2. Wang, W.; Jiang, B.; Sun, H.; Ru, X.; Sun, D.; Wang, L.; Wang, L.; Jiang, Y.; Li, Y.; Wang, Y. et al. Prevalence, incidence, and mortality of stroke in China: Results from a nationwide population-based survey of 480687 adults. Circulation; 2017; 135, pp. 759-771. [DOI: https://dx.doi.org/10.1161/CIRCULATIONAHA.116.025250] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28052979]
3. Della-Morte, D.; Guadagni, F.; Palmirotta, R.; Testa, G.; Caso, V.; Paciaroni, M.; Abete, P.; Rengo, F.; Ferroni, P.; Sacco, R.L. et al. Genetics of ischemic stroke, stroke-related risk factors, stroke precursors and treatments. Pharmacogenomics; 2012; 13, pp. 595-613. [DOI: https://dx.doi.org/10.2217/pgs.12.14] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22462751]
4. Abbate, R.; Sticchi, E.; Fatini, C. Genetics of cardiovascular disease. Clin. Cases Min. Bone Metab.; 2008; 5, pp. 63-66.
5. Reiner, Z.; Catapano, A.L.; De Backer, G.; Graham, I.; Taskinen, M.R.; Wiklund, O.; Agewall, S.; Alegria, E.; Chapman, M.J.; Durrington, P. et al. ESC/EAS Guidelines for the management of dyslipidaemias: The task force for the management of dyslipidaemias of the European Society of Cardiology (ESC) and the European Atherosclerosis Society (EAS). Eur. Heart J.; 2011; 32, pp. 1769-1818.
6. Abu el Maaty, M.A.; Hassanein, S.I.; Sleem, H.M.; Gad, M.Z. Vitamin D receptor gene polymorphisms (TaqI and ApaI) in relation to 25-hydroxyvitamin D levels and coronary artery disease incidence. J. Recept. Signal Transduct. Res.; 2015; 35, pp. 391-396. [DOI: https://dx.doi.org/10.3109/10799893.2014.959593]
7. He, L.; Wang, M. Association of vitamin d receptor-a gene polymorphisms with coronary heart disease in Han Chinese. Int. J. Clin. and Exp. Med.; 2015; 8, pp. 6224-6229.
8. Wang, E.W.-L.; Pang, M.Y.; Siu, P.M.; Lai, C.K.; Woo, J.; Collins, A.R.; Benzie, I.F. Vitamin D status and cardiometabolic risk factors in young adults in Hong Kong: Associations and implications. Asia Pac. J. Clin. Nutr.; 2018; 27, pp. 231-237.
9. Somjen, D.; Weisman, Y.; Kohen, F.; Gayer, B.; Limor, R.; Sharon, O.; Jaccard, N.; Knoll, E.; Stern, N. 25-hydroxyvitamin D3-1alpha-hydroxylase is expressed in human vascular smooth muscle cells and is upregulated by parathyroid hormone and estrogenic compounds. Circulation; 2005; 111, pp. 1666-1671. [DOI: https://dx.doi.org/10.1161/01.CIR.0000160353.27927.70]
10. Grundmann, M.; Haidar, M.; Placzko, S.; Niendorf, R.; Darashchonak, N.; Hubel, C.A.; von Versen-Höynck, F. Vitamin D improves the angiogenic properties of endothelial progenitor cells. Am. J. Physiol. Cell Physiol.; 2012; 303, pp. 954-962. [DOI: https://dx.doi.org/10.1152/ajpcell.00030.2012]
11. Molinari, C.; Uberti, F.; Grossini, E.; Vacca, G.; Carda, S.; Invernizzi, M.; Cisari, C. 1α,25-Dihydroxycholecalciferol Induces Nitric Oxide Production in Cultured Endothelial Cells. Cell. Physiol. Biochem.; 2011; 27, pp. 661-668. [DOI: https://dx.doi.org/10.1159/000330075]
12. Tenderich, G.; Berthorld, H.K.; Stehle, P.; Koerfer, R.; Schleithoff, S.S.; Zittermann, A. Vitamin D supplementation improves cytokine profiles in patients with congestive heart failure: A double-blind, randomized, placebo-controlled trial. Am. J. Clin. Nutr.; 2006; 83, pp. 754-759.
13. Li, Y.C.; Kong, J.; Wei, M.; Chen, Z.-F.; Liu, S.Q.; Cao, L.-P. 1,25-Dihydroxyvitamin D(3) is a negative endocrine regulator of the renin-angiotensin system. J. Clin. Investig.; 2002; 110, pp. 229-238. [DOI: https://dx.doi.org/10.1172/JCI0215219] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12122115]
14. Li, Y.C.; Qiao, G.; Uskokovic, M.; Xiang, W.; Zheng, W.; Kong, J. Vitamin D: A negative endocrine regulator of the renin-angiotensin system and blood pressure. J. Steroid Biochem. Mol. Biol.; 2004; 89–90, pp. 387-392. [DOI: https://dx.doi.org/10.1016/j.jsbmb.2004.03.004]
15. Rammos, G.; Tseke, P.; Ziakka, S. Vitamin D, the renin-angiotensin system, and insulin resistance. Int. Urol. Nephrol.; 2008; 40, pp. 419-426. [DOI: https://dx.doi.org/10.1007/s11255-007-9244-4]
16. Oh, J.; Riek, A.E.; Darwech, I.; Funai, K.; Shao, J.; Chin, K.; Sierra, O.L.; Carmeliet, G.; Ostlund, R.E., Jr.; Bernal-Mizrachi, C. Deletion of macrophage Vitamin D receptor promotes insulin resistance and monocyte cholesterol transport to accelerate atherosclerosis in mice. Cell Rep.; 2015; 10, pp. 1872-1886. [DOI: https://dx.doi.org/10.1016/j.celrep.2015.02.043] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25801026]
17. Kumar, S.; Nanduri, R.; Bhagyaraj, E.; Kalra, R.; Ahuja, N.; Chacko, A.P.; Tiwari, D.; Sethi, K.; Saini, A.; Chandra, V. et al. Vitamin D3-VDR-PTPN6 axis mediated autophagy contributes to the inhibition of macrophage foam cell formation. Autophagy; 2020; 17, pp. 2273-2289. [DOI: https://dx.doi.org/10.1080/15548627.2020.1822088] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32917126]
18. Mozos, I.; Marginean, O. Links between vitamin D deficiency and cardiovascular diseases. Biomed. Res. Int.; 2015; 2015, 109275. [DOI: https://dx.doi.org/10.1155/2015/109275] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26000280]
19. Pineda Lancheros, L.E.; Pérez Ramírez, C.; Sánchez Martín, A.; Gálvez Navas, J.M.; Martínez Martínez, F.; Ramírez Tortosa, M.d.C.; Jiménez Morales, A. Impact of Genetic Polymorphisms on the Metabolic Pathway of Vitamin D and Survival in Non-Small Cell Lung Cancer. Nutrients; 2021; 13, 3783. [DOI: https://dx.doi.org/10.3390/nu13113783]
20. Wacker, M.; Holick, M.F. Sunlight and vitamin D: A global perspective for health. Dermato-Endocrinology; 2013; 5, pp. 51-108. [DOI: https://dx.doi.org/10.4161/derm.24494]
21. Borel, P.; Caillaud, D.; Cano, N. Vitamin D bioavailability: State of the art. Crit. Rev. Food Sci. Nutr.; 2013; 55, pp. 1193-1205. [DOI: https://dx.doi.org/10.1080/10408398.2012.688897]
22. Pilz, S.; Verheyen, N.; Grübler, M.R.; Tomaschitz, A.; März, W. Vitamin D and cardiovascular disease prevention. Nat. Rev. Cardiol.; 2016; 13, 404.
23. Christakos, S.; Dhawan, P.; Verstuyf, A.; Verlinden, L.; Carmeliet, G. Vitamin D: Metabolism, molecular mechanism of action, and pleiotropic effects. Physiol. Rev.; 2016; 96, pp. 365-408. [DOI: https://dx.doi.org/10.1152/physrev.00014.2015] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26681795]
24. Haussler, M.R.; Whitfield, G.K.; Kaneko, I.; Haussler, C.A.; Hsieh, D.; Hsieh, J.-C.; Jurutka, P.W. Molecular Mechanisms of Vitamin D Action. Calcif. Tissue Int.; 2013; 92, pp. 77-98. [DOI: https://dx.doi.org/10.1007/s00223-012-9619-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22782502]
25. Rachez, C.; Lemon, B.D.; Suldan, Z.; Bromleigh, V.; Gamble, M.J.; Näär, A.M.; Erdjument-Bromage, H.; Tempst, P.; Freedman, L.P. Ligand-dependent transcription activation by nuclear receptors requires the DRIP complex. Nature; 1999; 398, pp. 824-828. [DOI: https://dx.doi.org/10.1038/19783]
26. Bouillon, R.; Carmeliet, G.; Verlinden, L.; van Etten, E.; Verstuyf, A.; Luderer, H.F.; Lieben, L.; Mathieu, C.; Demay, M. Vitamin D and human health: Lessons from vitamin D receptor null mice. Endocr. Rev.; 2008; 29, pp. 726-776. [DOI: https://dx.doi.org/10.1210/er.2008-0004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18694980]
27. Sangkaew, B.; Nuinoon, M.; Jeenduang, N. Association of vitamin D receptor gene polymorphisms with serum 25(OH)D levels and metabolic syndrome in Thai population. Gene; 2018; 659, pp. 59-66. [DOI: https://dx.doi.org/10.1016/j.gene.2018.03.047] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29555202]
28. Vaidya, A.; Sun, B.; Forman, J.P.; Hopkins, P.N.; Brown, N.J.; Kolatkar, N.S.; Williams, G.H.; Williams, J.S. The Fok1 vitamin D receptor gene polymorphism is associated with plasma renin activity in Caucasians. Clin. Endocrinol.; 2011; 74, pp. 783-790. [DOI: https://dx.doi.org/10.1111/j.1365-2265.2011.03991.x]
29. Uitterlinden, A.G.; Fang, Y.; Van Meurs, J.B.; Pols, H.A.; Van Leeuwen, J.P. Genetics and biology of vitamin D receptor polymorphisms. Gene; 2004; 338, pp. 143-156. [DOI: https://dx.doi.org/10.1016/j.gene.2004.05.014]
30. Hong, Y.J.; Kang, E.S.; Ji, M.J.; Choi, H.J.; Oh, T.; Koong, S.-S.; Jeon, H.J. Association between Bsm1 Polymorphism in Vitamin D Receptor Gene and Diabetic Retinopathy of Type 2 Diabetes in Korean Population. Endocrinol. Metab.; 2015; 30, pp. 469-474. [DOI: https://dx.doi.org/10.3803/EnM.2015.30.4.469]
31. Margolis, R.N.; Christakos, S. The nuclear receptor superfamily of steroid hormones and vitamin D gene regulation. An update. Ann. N. Y. Acad. Sci.; 2010; 1192, pp. 208-214. [DOI: https://dx.doi.org/10.1111/j.1749-6632.2009.05227.x]
32. Sygitowicz, G.; Pera, L.; Sitkiewicz, D. Vitamin D receptor (VDR) polymorphism and the risk of cardiovascular events. Kardiol Pol.; 2014; 72, pp. 64-66. [DOI: https://dx.doi.org/10.5603/KP.2014.0005] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24469749]
33. Serrano, J.C.E.; De Lorenzo, D.; Cassanye, A.; Martín-Gari, M.; Espinel, A.; Delgado, M.A.; Pamplona, R.; Portero-Otin, M. Vitamin D receptor BsmI polymorphism modulates soy intake and 25-hydroxyvitamin D supplementation benefits in cardiovascular disease risk factors profile. Genes Nutr.; 2013; 8, pp. 561-569. [DOI: https://dx.doi.org/10.1007/s12263-013-0356-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24158768]
34. Zhao, Y.; Liao, S.; He, J.; Jin, Y.; Fu, H.; Chen, X.; Fan, X.; Xu, H.; Liu, X.; Jin, J. et al. Association of vitamin D receptor gene polymorphisms with metabolic syndrome: A case-control design of population-based cross-sectional study in North China. Lipids Health Dis.; 2014; 13, 129. [DOI: https://dx.doi.org/10.1186/1476-511X-13-129] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25106919]
35. Dorsch, M.P.; Nemerovski, C.W.; Ellingrod, V.L.; Cowger, J.A.; Dyke, D.B.; Koelling, T.M.; Wu, A.H.; Aaronson, K.D.; Simpson, R.U.; Bleske, B.E. Vitamin D Receptor Genetics on Extracellular Matrix Biomarkers and Hemodynamics in Systolic Heart Failure. J. Cardiovasc. Pharmacol. Ther.; 2014; 19, pp. 439-445. [DOI: https://dx.doi.org/10.1177/1074248413517747]
36. Mahjoubi, I.; Kallel, A.; Sbaï, M.H.; Ftouhi, B.; Ben Halima, M.; Jemaa, Z.; Moncef, F.; Hedia, S.; Riadh, J.; Naziha, K. Lack of association between FokI polymorphism in vitamin D receptor gene (VDR) & type 2 diabetes mellitus in the Tunisian population. Indian J. Med. Res.; 2016; 144, pp. 46-51.
37. Hajj, A.; Chedid, R.; Chouery, E.; Megarbané, A.; Gannagé-Yared, M.-H. Relationship between vitamin D receptor gene polymorphisms, cardiovascular risk factors and adiponectin in a healthy young population. Pharmacogenomics; 2016; 17, pp. 1675-1686. [DOI: https://dx.doi.org/10.2217/pgs-2016-0045]
38. Bhakat, R.; Chandra, L.; Saxena, A.; Sarda, A.K.; Krishnamurthy, K.; Yadav, P. Evaluation of metabolic syndrome and vitamin D receptor gene polymorphism in male factor infertility. Indian J. Clin. Biochem.; 2017; 32, pp. 468-472. [DOI: https://dx.doi.org/10.1007/s12291-017-0632-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29062180]
39. Vuolo, L.; Di Somma, C.; Faggiano, A.; Colao, A. Vitamin D and cancer. Front. Endocrinol.; 2012; 3, pp. 1-13. [DOI: https://dx.doi.org/10.3389/fendo.2012.00058]
40. Jurutka, P.W.; Remus, L.S.; Whitfield, G.K.; Thompson, P.D.; Hsieh, J.C.; Zitzer, H.; Tavakkoli, P.; Galligan, M.A.; Dang, H.T.; Haussler, C.A. The polymorphic N terminus in human vitamin D receptor isoforms influences transcriptional activity by modulating interaction with transcription factor IIB. Mol. Endocrinol.; 2000; 14, pp. 401-420. [DOI: https://dx.doi.org/10.1210/mend.14.3.0435]
41. Decker, C.J.; Parker, R. Diversity of cytoplasmic functions for the 3′ untranslated region of eukaryotic transcripts. Curr Opin Cell Biol.; 1995; 7, pp. 386-392. [DOI: https://dx.doi.org/10.1016/0955-0674(95)80094-8]
42. Morrison, N.A.; Qi, J.C.; Tokita, A.; Kelly, P.J.; Crofts, L.; Nguyen, T.V.; Sambrook, P.N.; Eisman, J.A. Prediction of bone density from vitamin D receptor alleles. Nature; 1994; 367, pp. 284-287. [DOI: https://dx.doi.org/10.1038/367284a0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8161378]
43. Györffy, B.; Vásárhelyi, B.; Krikovszky, D.; Madácsy, L.; Tordai, A.; Tulassay, T.; Szabó, A. Gender-specific association of vitamin D receptor polymorphism combinations with type 1 diabetes mellitus. Eur. J. Endocrinol.; 2002; 147, pp. 803-808. [DOI: https://dx.doi.org/10.1530/eje.0.1470803] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12457456]
44. Arai, H.; Miyamoto, K.I.; Yoshida, M.; Yamamoto, H.; Taketani, Y.; Morita, K.; Kubota, M.; Yoshida, S.; Ikeda, M.; Watabe, F. et al. The Polymorphism in the caudal-related homeodomain protein Cdx-2 binding element in the human vitamin D receptor gene. J. Bone Min. Res.; 2001; 16, pp. 1256-1264. [DOI: https://dx.doi.org/10.1359/jbmr.2001.16.7.1256]
45. Yamamoto, H.; Miyamoto, K.-I.; Li, B.; Taketani, Y.; Kitano, M.; Inoue, Y.; Morita, K.; Pike, J.W.; Takeda, E. The Caudal-Related Homeodomain Protein Cdx-2 Regulates Vitamin D Receptor Gene Expression in the Small Intestine. J. Bone Miner. Res.; 1999; 14, pp. 240-247. [DOI: https://dx.doi.org/10.1359/jbmr.1999.14.2.240] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9933478]
46. Nakhl, S.; Sleilaty, G.; Chouery, E.; Salem, N.; Chaine, R.; Farès, N. FokI vitamin D receptor gene polymorphism and serum 25-hydroxyvitamin D in patients with cardiovascular risk. Arch. Med. Sci. Atheroscler. Dis.; 2019; 4, pp. e298-e303. [DOI: https://dx.doi.org/10.5114/amsad.2019.91437]
47. Hao, Y.; Chen, Y. Vitamin D levels and vitamin D receptor variants are associated with chronic heart failure in Chinese patients. J. Clin. Lab. Anal.; 2019; 33, e22847. [DOI: https://dx.doi.org/10.1002/jcla.22847]
48. Jun, M.; Xue-Qiang, G.; Jia, L.; Yang-Jing, X.; Cheng, Z.; Ge, J. Interactions between vitamin D receptor (VDR) gene and Interleukin-6 gene and environment factors on coronary heart disease risk in a Chinese Han population. Oncotarget; 2017; 8, pp. 78419-78428. [DOI: https://dx.doi.org/10.18632/oncotarget.19472]
49. Prabhakar, P.; Majumdar, V.; Kulkarni, G.B.; Christopher, R. Genetic variants of vitamin D receptor and susceptibility to ischemic stroke. Biochem. Biophys. Res. Commun.; 2015; 456, pp. 631-636. [DOI: https://dx.doi.org/10.1016/j.bbrc.2014.12.007]
50. Mokhtar, W.A.; Fawzy, A.; Allam, R.M.; Amer, R.M.; Hamed, M.S. Maternal vitamin D level and vitamin D receptor gene polymorphism as a risk factor for congenital heart diseases in offspring; An Egyptian case-control study. Genes Dis.; 2018; 6, pp. 193-200. [DOI: https://dx.doi.org/10.1016/j.gendis.2018.08.001]
51. Lu, S.; Guo, S.; Hu, F.; Guo, Y.; Yan, L.; Ma, W.; Wang, Y.; Wei, Y.; Zhang, Z.; Wang, Z. The associations between the polymorphisms of vitamin D receptor and coronary artery disease: A systematic review and meta-analysis. Medicine; 2016; 95, e3467. [DOI: https://dx.doi.org/10.1097/MD.0000000000003467]
52. Huzmeli, C.; Bagci, G.; Candan, F.; Bagci, B.; Akkaya, L.; Kayatas, M. Association of vitamin D receptor gene TaqI, FokI and ApaI variants with arteriovenous fistula failure in hemodialysis patients. J. Vasc. Access; 2018; 19, pp. 303-310. [DOI: https://dx.doi.org/10.1177/1129729817752860] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29544394]
53. Wilke, R.A.; Simpson, R.U.; Mukesh, B.N.; Bhupathi, S.V.; Dart, R.A.; Ghebranious, N.R.; McCarty, C.A. Genetic variation in CYP27B1 is associated with congestive heart failure in patients with hypertension. Pharmacogenomics; 2009; 10, pp. 1789-1797. [DOI: https://dx.doi.org/10.2217/pgs.09.101] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19891555]
54. Alizadeh, S.; Djafarian, K.; Alizadeh, H.; Mohseni, R.; Shab-Bidar, S. Common variants of vitamin D receptor gene polymorphisms and susceptibility to coronary artery disease: A systematic review and meta-analysis. J. Nutr. Nutr.; 2017; 10, pp. 9-18. [DOI: https://dx.doi.org/10.1159/000455914] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28351026]
55. Abouzid, M.; Kruszyna, M.; Burchardt, P.; Kruszyna, L.; Główka, F.K.; Karazniewicz-Łada, M. Vitamin D Receptor Gene Polymorphism and Vitamin D Status in Population of Patients with Cardiovascular Disease—A Preliminary Study. Nutrients; 2021; 13, 3117. [DOI: https://dx.doi.org/10.3390/nu13093117]
56. Ewida, S.M.; Salem, A.; Shaker, Y.M.; Samy, N.; Yassen, I.; Hassan-Mohamed, R. Vitamin D levels and vitamin D receptor genetic variants in Egyptian cardiovascular disease patients with and without diabetes. Egypt. J. Med. Hum. Genet.; 2021; 22, 55. [DOI: https://dx.doi.org/10.1186/s43042-021-00174-9]
57. Tabaei, S.; Motallebnezhad, M.; Tabaee, S.S. Vitamin D Receptor (VDR) Gene Polymorphisms and Risk of Coronary Artery Disease (CAD): Systematic Review and Meta-analysis. Biochem. Genetics; 2021; 59, pp. 813-836. [DOI: https://dx.doi.org/10.1007/s10528-021-10038-x]
58. The American Diabetes Association (ADA). Classification and diagnosis of diabetes: Standards of medical care in diabetes—2019. Diabetes Care; 2019; 42, (Suppl. S1), pp. S13-S28. [DOI: https://dx.doi.org/10.2337/dc19-S002]
59. Ochoa Mangado, E.; Madoz-Gúrpide, A.; Vicente Muelas, N. Diagnóstico y tratamiento de la dependencia de alcohol. Med. Segur. Trab.; 2009; 55, pp. 26-40. [DOI: https://dx.doi.org/10.4321/S0465-546X2009000100003]
60. Salas-Salvadó, J.; Rubio, M.A.; Barbany, M.; Moreno, B. Consenso SEEDO 2007 para la evaluación del sobrepeso y la obesidad y el establecimiento de criterios de intervención terapéutica. Med. Clin.; 2007; 128, pp. 184-196. [DOI: https://dx.doi.org/10.1016/S0025-7753(07)72531-9]
61. Catapano, A.L.; Graham, I.; De Backer, G.; Wiklund, O.; Chapman, M.J.; Drexel, H.; Hoes, A.; Jennings, C.; Landmesser, U.; Pedersen, T. Guía ESC/EAS 2016 sobre el tratamiento de las dislipemias. Rev. Esp. Cardiol.; 2017; 70, pp. 1151-11564. [DOI: https://dx.doi.org/10.1016/j.recesp.2016.11.052]
62. R Development Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020; Available online: https://www.R-project.org/ (accessed on 6 May 2022).
63. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.; Bender, D.; Maller, J.; Sklar, P.; de Bakker, P.I.; Daly, M.J. et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet.; 2007; 81, pp. 559-575. [DOI: https://dx.doi.org/10.1086/519795] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17701901]
64. Solé, X.; Guinó, E.; Valls, J.; Iniesta, R.; Moreno, V. SNPStats: A web tool for the analysis of association studies. Bioinformatics; 2006; 22, pp. 1928-1929. [DOI: https://dx.doi.org/10.1093/bioinformatics/btl268] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16720584]
65. Barrett, J.C.; Fry, B.; Maller, J.; Daly, M.J. Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics; 2005; 21, pp. 263-265. [DOI: https://dx.doi.org/10.1093/bioinformatics/bth457] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15297300]
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Cristina Pérez Ramírez 2 ; Gálvez Navas, José María 2
; Pineda Lancheros, Laura Elena 1
; Susana Rojo Tolosa 1
; María del Carmen Ramírez Tortosa 3
; Alberto Jiménez Morales 1 1 Pharmacogenetics Unit, Pharmacy Service, University Hospital Virgen de las Nieves, 18014 Granada, Spain
2 Pharmacogenetics Unit, Pharmacy Service, University Hospital Virgen de las Nieves, 18014 Granada, Spain; Department of Biochemistry and Molecular Biology II, Faculty of Pharmacy, Campus Universitario de Cartuja, University of Granada, 18071 Granada, Spain
3 Department of Biochemistry and Molecular Biology II, Faculty of Pharmacy, Campus Universitario de Cartuja, University of Granada, 18071 Granada, Spain