Next Article in Journal
Phylogenomics-Based Reconstruction and Molecular Evolutionary Histories of Brassica Photoreceptor Gene Families
Previous Article in Journal
Gene Variants Related to Cardiovascular and Pulmonary Diseases May Correlate with Severe Outcome of COVID-19
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 15
10.3390/ijms23158686
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Vitamin D-Related Single Nucleotide Polymorphisms as Risk Biomarker of Cardiovascular Disease

by
Paula González Rojo
1,
Cristina Pérez Ramírez
1,2,*,
José María Gálvez Navas
1,2,†,
Laura Elena Pineda Lancheros
1,†,
Susana Rojo Tolosa
1,†,
María del Carmen Ramírez Tortosa
2,‡ and
Alberto Jiménez Morales
1,‡
1
Pharmacogenetics Unit, Pharmacy Service, University Hospital Virgen de las Nieves, 18014 Granada, Spain
2
Department of Biochemistry and Molecular Biology II, Faculty of Pharmacy, Campus Universitario de Cartuja, University of Granada, 18071 Granada, Spain
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to the work.
‡
These authors contributed equally to the work.
Int. J. Mol. Sci. 2022, 23(15), 8686; https://doi.org/10.3390/ijms23158686
Submission received: 14 July 2022 / Revised: 1 August 2022 / Accepted: 1 August 2022 / Published: 4 August 2022
(This article belongs to the Section Molecular Pathology, Diagnostics, and Therapeutics)
Browse Figure
Versions Notes

Abstract

:
Cardiovascular diseases (CVDs) are a group of disorders of the heart and blood vessels. In addition to environmental risk factors, genetic predisposition increases the risk; this includes alterations in the vitamin D receptor gene (VDR). These alterations play a key role in modifying vitamin D uptake, being able to modify its function and increasing susceptibility to cardiovascular disorders. The aim of this study was to evaluate the association of polymorphisms in the VDR gene and risk of CVD in a Caucasian population. A retrospective case-control study was conducted comprising 246 CVD patients and 246 controls of Caucasian origin from Southern Spain. 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. The logistic regression analysis adjusted for body mass index 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) and the recessive model (p = 0.0099; OR = 2.71; 95% CI = 1.31–6.07; TT vs. C). Haplotype analysis revealed that the haplotype GAC (p = 0.047; OR = 0.34; 95% CI = 0.12–0.98) was associated with increased risk of CVD. The VDR polymorphisms FokI (rs2228570) was significantly associated with the development of CVD. No influence was observed of the VDR polymorphisms BsmI (rs1544410), TaqI (rs731236), ApaI (rs7975232) and Cdx2 (rs11568820) on the risk of developing CVD in the patients studied.

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.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms23158686/s1.

Author Contributions

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.

Funding

This research was funded by the ERDF funds (EU) from the Instituto de Salud Carlos III (PT13/0010/0039) by cofounding grants the Virgen de las Nieves University Hospital Biobank.

Institutional Review Board Statement

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 Statement

Informed consent was obtained from all subjects involved in the study.

Acknowledgments

The results of this study are part of the doctoral thesis presented by Paula González Rojo at the University of Granada.

Conflicts of Interest

The authors declare no conflict of interest.

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, 759–771. [Google Scholar] [CrossRef] [PubMed]
  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, 595–613. [Google Scholar] [CrossRef] [PubMed]
  4. Abbate, R.; Sticchi, E.; Fatini, C. Genetics of cardiovascular disease. Clin. Cases Min. Bone Metab. 2008, 5, 63–66. [Google Scholar]
  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, 1769–1818. [Google Scholar]
  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, 391–396. [Google Scholar] [CrossRef]
  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, 6224–6229. [Google Scholar]
  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, 231–237. [Google Scholar]
  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, 1666–1671. [Google Scholar] [CrossRef] [Green Version]
  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, 954–962. [Google Scholar] [CrossRef]
  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, 661–668. [Google Scholar] [CrossRef]
  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, 754–759. [Google Scholar]
  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, 229–238. [Google Scholar] [CrossRef] [PubMed]
  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, 387–392. [Google Scholar] [CrossRef]
  15. Rammos, G.; Tseke, P.; Ziakka, S. Vitamin D, the renin-angiotensin system, and insulin resistance. Int. Urol. Nephrol. 2008, 40, 419–426. [Google Scholar] [CrossRef]
  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, 1872–1886. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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, 2273–2289. [Google Scholar] [CrossRef] [PubMed]
  18. Mozos, I.; Marginean, O. Links between vitamin D deficiency and cardiovascular diseases. Biomed. Res. Int. 2015, 2015, 109275. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef]
  20. Wacker, M.; Holick, M.F. Sunlight and vitamin D: A global perspective for health. Dermato-Endocrinology 2013, 5, 51–108. [Google Scholar] [CrossRef] [Green Version]
  21. Borel, P.; Caillaud, D.; Cano, N. Vitamin D bioavailability: State of the art. Crit. Rev. Food Sci. Nutr. 2013, 55, 1193–1205. [Google Scholar] [CrossRef]
  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. [Google Scholar]
  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, 365–408. [Google Scholar] [CrossRef] [PubMed]
  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, 77–98. [Google Scholar] [CrossRef] [PubMed]
  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, 824–828. [Google Scholar] [CrossRef]
  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, 726–776. [Google Scholar] [CrossRef] [PubMed]
  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, 59–66. [Google Scholar] [CrossRef] [PubMed]
  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, 783–790. [Google Scholar] [CrossRef]
  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, 143–156. [Google Scholar] [CrossRef] [Green Version]
  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, 469–474. [Google Scholar] [CrossRef] [Green Version]
  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, 208–214. [Google Scholar] [CrossRef]
  32. Sygitowicz, G.; Pera, L.; Sitkiewicz, D. Vitamin D receptor (VDR) polymorphism and the risk of cardiovascular events. Kardiol Pol. 2014, 72, 64–66. [Google Scholar] [CrossRef] [PubMed]
  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, 561–569. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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, 439–445. [Google Scholar] [CrossRef] [Green Version]
  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, 46–51. [Google Scholar]
  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, 1675–1686. [Google Scholar] [CrossRef] [Green Version]
  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, 468–472. [Google Scholar] [CrossRef] [PubMed]
  39. Vuolo, L.; Di Somma, C.; Faggiano, A.; Colao, A. Vitamin D and cancer. Front. Endocrinol. 2012, 3, 1–13. [Google Scholar] [CrossRef] [Green Version]
  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, 401–420. [Google Scholar] [CrossRef]
  41. Decker, C.J.; Parker, R. Diversity of cytoplasmic functions for the 3′ untranslated region of eukaryotic transcripts. Curr Opin Cell Biol. 1995, 7, 386–392. [Google Scholar] [CrossRef]
  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, 284–287. [Google Scholar] [CrossRef] [PubMed]
  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, 803–808. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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, 1256–1264. [Google Scholar] [CrossRef]
  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, 240–247. [Google Scholar] [CrossRef] [PubMed]
  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, e298–e303. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef] [Green Version]
  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, 78419–78428. [Google Scholar] [CrossRef] [Green Version]
  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, 631–636. [Google Scholar] [CrossRef]
  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, 193–200. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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, 303–310. [Google Scholar] [CrossRef] [PubMed]
  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, 1789–1797. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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, 9–18. [Google Scholar] [CrossRef] [PubMed]
  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. [Google Scholar] [CrossRef]
  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. [Google Scholar] [CrossRef]
  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, 813–836. [Google Scholar] [CrossRef]
  58. The American Diabetes Association (ADA). Classification and diagnosis of diabetes: Standards of medical care in diabetes—2019. Diabetes Care 2019, 42 (Suppl. S1), S13–S28. [Google Scholar] [CrossRef] [Green Version]
  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, 26–40. [Google Scholar] [CrossRef] [Green Version]
  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, 184–196. [Google Scholar] [CrossRef]
  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, 1151–11564. [Google Scholar] [CrossRef]
  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, 559–575. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Solé, X.; Guinó, E.; Valls, J.; Iniesta, R.; Moreno, V. SNPStats: A web tool for the analysis of association studies. Bioinformatics 2006, 22, 1928–1929. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Barrett, J.C.; Fry, B.; Maller, J.; Daly, M.J. Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics 2005, 21, 263–265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ijms 23 08686 g001 550
Figure 1. Chromosomal location for (vitamin D receptor) VDR, and linkage disequilibrium (LD).
Figure 1. Chromosomal location for (vitamin D receptor) VDR, and linkage disequilibrium (LD).
Ijms 23 08686 g001
Table
Table 1. Clinico-pathologic characteristics of cardiology cases and controls.
Table 1. Clinico-pathologic characteristics of cardiology cases and controls.
CasesControlsχ²p-ValueReferenceOR95% CI
Nn (%)Nn (%)
Gender246 246
Male 114 (46.3) 104 (42.3)0.82370.3641
Female 132 (53.7) 142 (57.7)
Age24670 (63.25, 76)24669.5 (63, 76)
Tobacco consumption233 203
Current-smokers 41 (17.6) 27 (13.3)1.95020.3772
Former-smokers 70 (30.0) 70 (34.5)
Non-smokers 122 (52.4) 106 (52.2)
Alcohol consumption220 173
Current-drinkers 46 (20.9) 27 (15.6)2.08230.3531
Former-drinkers 8 (3.6) 5 (2.9)
Non-drinkers 166 (75.5) 141 (81.5)
Body Mass Index18229.63 ± 5.1411428.03 ± 4.89
Underweight 1 (0.5) 1 (0.9) 0.0066 aNormal (healthy weight)1.210.05–31.44
Normal (healthy weight) 33 (18.1) 40 (35.1)1
Overweight 68 (37.4) 30 (26.3)2.741.47–5.20
Obese Class I (Moderately obese) 46 (25.3) 32 (28.1)1.740.92–3.34
Obese Class II (Severely obese) 29 (15.9) 10 (8.8)3.511.53–8.57
Obese Class III (Very severely obese) 5 (2.7) 1 (0.9)6.060.92–119.19
Dyslipidemia246 246
No 158 (64.2) 155 (63.0)0.07900.7786
Yes 88 (35.8) 91 (37.0)
Hypertension246 246
No 100 (40.6) 115 (46.8)1.85880.1728
Yes 146 (48) 131 (53.2)
Diabetes246 246
No 170 (69.1) 201 (81.7)10.53200.0012 aNo1
Yes 76 (30.9) 45 (18.3)1.991.32–3.06
Disease246
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.
Table
Table 2. Influence of VDR FokI (rs2228570) gene polymorphism on risk of cardiovascular disease.
Table 2. Influence of VDR FokI (rs2228570) gene polymorphism on risk of cardiovascular disease.
ModelsGenotypeCases [n (%)]Controls [n (%)]p-Value aAdjusted
p-Value b		OR c95% CI
GenotypicTT45 (18.3)21 (8.5)0.00610.03052.351.33–4.28
CT99 (40.2)113 (45.9)0.960.66–1.41
CC102 (41.5)112 (45.5)1
DominantT144 (58.5)134 (54.5)0.36311
CC102 (41.5)112 (45.5)
RecessiveTT45 (18.3)21 (8.5)0.00220.00752.391.39–4.24
C201 (81.7)225 (91.5)
AllelicT189 (38.41)155 (31.50)0.02300.1151
C303 (61.59)337 (68.50)
Additive---0.02673 d0.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.
Table
Table 3. Influence of clinical characteristic and VDR FokI (rs2228570) gene polymorphism on risk of cardiovascular disease.
Table 3. Influence of clinical characteristic and VDR FokI (rs2228570) gene polymorphism on risk of cardiovascular disease.
GenotypicDominantRecessiveAdditive
TT vs. CCCT vs. CCT vs. CCTT vs. CT vs. C
p-ValueOR95% CIp-ValueOR95% CIp-ValueOR95% CIp-ValueOR95% CIp-ValueOR95% CI
Body Mass Index0.02651.061.01–1.120.02651.061.01–1.120.04581.051.00–1.110.04131.051.00–1.110.06341.050.99–1.10
Diabetes
Yes0.03591.911.06–3.550.03591.911.06–3.550.04681.831.02–3.380.03161.941.07–3.590.03721.891.05–3.50
VDR FokI (rs2228570)0.04302.301.06–5.370.22860.720.43–1.220.86480.950.58–1.570.00992.711.31–6.070.21511.240.88–1.76
Shade means the value is significant.
Table
Table 4. Haplotype association with response.
Table 4. Haplotype association with response.
Rs1544410Rs7975232Rs731236FreqOR (95% CI)p-Value
1GCT0.46131.00---
2AAC0.36850.92 (0.70–1.22)0.570
3GAT0.10790.81 (0.52–1.24)0.330
4AAT0.02511.11 (0.51–2.38)0.800
5GAC0.01890.34 (0.12–0.98)0.047
6ACT0.01141.22 (0.41–3.65)0.730
Rare***0.00690.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.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

González Rojo, P.; Pérez Ramírez, C.; Gálvez Navas, J.M.; Pineda Lancheros, L.E.; Rojo Tolosa, S.; Ramírez Tortosa, M.d.C.; Jiménez Morales, A. Vitamin D-Related Single Nucleotide Polymorphisms as Risk Biomarker of Cardiovascular Disease. Int. J. Mol. Sci. 2022, 23, 8686. https://doi.org/10.3390/ijms23158686

AMA Style

González Rojo P, Pérez Ramírez C, Gálvez Navas JM, Pineda Lancheros LE, Rojo Tolosa S, Ramírez Tortosa MdC, Jiménez Morales A. Vitamin D-Related Single Nucleotide Polymorphisms as Risk Biomarker of Cardiovascular Disease. International Journal of Molecular Sciences. 2022; 23(15):8686. https://doi.org/10.3390/ijms23158686

Chicago/Turabian Style

González Rojo, Paula, Cristina Pérez Ramírez, José María Gálvez Navas, Laura Elena Pineda Lancheros, Susana Rojo Tolosa, María del Carmen Ramírez Tortosa, and Alberto Jiménez Morales. 2022. "Vitamin D-Related Single Nucleotide Polymorphisms as Risk Biomarker of Cardiovascular Disease" International Journal of Molecular Sciences 23, no. 15: 8686. https://doi.org/10.3390/ijms23158686

APA Style

González Rojo, P., Pérez Ramírez, C., Gálvez Navas, J. M., Pineda Lancheros, L. E., Rojo Tolosa, S., Ramírez Tortosa, M. d. C., & Jiménez Morales, A. (2022). Vitamin D-Related Single Nucleotide Polymorphisms as Risk Biomarker of Cardiovascular Disease. International Journal of Molecular Sciences, 23(15), 8686. https://doi.org/10.3390/ijms23158686

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