Next Article in Journal
Effects of Five Amino Acids (Serine, Alanine, Glutamate, Aspartate, and Tyrosine) on Mental Health in Healthy Office Workers: A Randomized, Double-Blind, Placebo-Controlled Exploratory Trial
Next Article in Special Issue
Association of Vitamin D in Different Trimester with Hemoglobin during Pregnancy
Previous Article in Journal
Revisional Surgery of One Anastomosis Gastric Bypass for Severe Protein–Energy Malnutrition
Previous Article in Special Issue
Associations between Maternal Selenium Status and Cord Serum Vitamin D Levels: A Birth Cohort Study in Wuhan, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Association of Vitamin D and Its Pathway Genes’ Polymorphisms with Hypertensive Disorders of Pregnancy: A Prospective Cohort Study

1
Department of Public Health, and Department of Anesthesiology, Second Affiliated Hospital of Zhejiang University School of Medicine, Hangzhou 310058, China
2
Department of Epidemiology & Health Statistics, School of Public Health, School of Medicine, Zhejiang University, Hangzhou 310058, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 23 April 2022 / Revised: 27 May 2022 / Accepted: 2 June 2022 / Published: 6 June 2022

Abstract

:
Objective: We aimed to explore the effect of single nucleotide polymorphism (SNP) in the genes of the vitamin D (VitD) metabolic pathway and its interaction with VitD level during pregnancy on the development of hypertensive disorders of pregnancy (HDP). Methods: The study was conducted in the Zhoushan Maternal and Child Health Care Hospital, China, from August 2011 to May 2018. The SNPs in VitD metabolic pathway-related genes were genotyped. Plasma 25-hydroxyvitamin vitamin D (25(OH)D) levels was measured at first (T1), second (T2), and third (T3) trimesters. The information of systolic blood pressure (SBP) and diastolic blood pressure (DBP), and the diagnosis of HDP were extracted from the electronic medical record system. Multivariable linear and logistic regression models and crossover analysis were applied. Results: The prospective cohort study included 3699 pregnant women, of which 105 (2.85%) were diagnosed with HDP. After adjusting for potential confounders, VitD deficiency at T2, as well as the change of 25(OH)D level between T1 and T2, were negatively associated with DBP at T2 and T3, but not HDP. Polymorphisms in CYP24A1, GC, and LRP2 genes were associated with blood pressure and HDP. In addition, VitD interacted with CYP24A1, GC, and VDR genes’ polymorphisms on blood pressure. Furthermore, participants with polymorphisms in CYP24A1-rs2248137, LRP2-rs2389557, and LRP2-rs4667591 and who had VitD deficiency at T2 showed an increased risk of HDP. Conclusions: The individual and interactive association between VitD deficiency during pregnancy and SNPs in the genes of the VitD metabolic pathway on blood pressure and HDP were identified.

1. Introduction

Hypertensive disorders of pregnancy (HDP), including gestational hypertension, preeclampsia, eclampsia, pregnancy complicated with chronic hypertension, and chronic hypertension complicated with preeclampsia [1], accounted for nearly 18% of all maternal deaths worldwide [2]. Its increasing prevalence and related risks for maternal and child health as well as cardiovascular diseases later in life has garnered great attention in the field of public health [3,4]. The risk factors for HDP are advanced age, primipara, multiple pregnancy, family history of hypertension, high pre-pregnancy body mass index (BMI), and high basal blood pressure [5].
Approximately 5% to 7% of pregnancies are complicated by preeclampsia [6]. While the cause of preeclampsia is not fully discerned, previous studies have suggested that abnormal placentation and angiogenesis were central to the pathogenesis of this syndrome [6]. In recent years, growing evidence of the association between maternal hypovitaminosis D and increased risk of HDP has been suggested [7,8]. Compared to non-pregnant state, there are significant changes in vitamin D (VitD) metabolism during pregnancy, and the serum levels of VitD binding protein (VDBP) [9], as well as the active form, 1,25-dihydroxyvitamin (1,25(OH)2D) [10], increased notably. It is believed that not only the kidneys but also the placenta and decidua produce and secret 1,25(OH)2D during pregnancy [11]. Moreover, VitD receptors and related metabolic enzymes have been discovered in the placenta and decidua [12], indicating a potential role for VitD in implantation and placental function, outside of its well-established role in skeletal health [13].
To date, trial evidence appears insufficient to lean towards a protective effect of VitD supplementation during pregnancy against the risk of preeclampsia owing to small sample size or low study quality [14,15]. In addition, findings from observational studies in regard to the association between maternal VitD status and HDP are discrepant due to the large heterogeneity between study designs, lack of adherence to standardized outcome definitions, and different gestational weeks of VitD detection [8,16]. On the other hand, genetic variants in the VitD metabolic pathway have also been shown to participate in the pathogenesis of blood pressure increase and preeclampsia [8,17], which suggests a possible interaction between VitD and its pathway gene variants for HDP. The concentration or effect of VitD can be highly regulated due to the variation of key protein expression or activity. 25(OH)D is the main circulating metabolism and is considered the biological marker of VitD status. The main metabolic enzymes involved in the synthesis, transport, reabsorption, and inactivation of VitD include 25-hydroxylase (CYP3A4), 1-hydroxylase (CYP27B1), vitamin D-binding protein (GC), 24-hydroxylas and metaling (LRP2), and 24-hydroxylase (CYP24A1). Moreover, VitD receptor (VDR) regulates VitD metabolism through binding 1,25(OH)2D [18].
So far, most studies have only focused on the relationship between VitD status during pregnancy or gene variation in the VitD metabolic pathway and HDP, without considering the possible interaction between them. This study aimed to explore the association of VitD status in three trimesters of pregnancy with the risk of HDP, and to explore the interactive effect between maternal VitD level and genetic variants in the VitD metabolic pathways (GC, CYP24A1, CYP3A4, CYP27B1, LRP2, VDR) on gestational blood pressure and HDP.

2. Materials and Methods

2.1. Study Design and Participants

The Zhoushan Pregnant Women Cohort (ZPWC) is an ongoing prospective cohort, conducted in Zhoushan Maternal and Child Health Care Hospital, Zhejiang. This study was based on the data of ZPWC from August 2011 to May 2018. We recruited pregnant women aged between 18 and 45 years at their first prenatal visit. A more detailed description of the inclusion and exclusion criteria can be seen in a previous study [19]. In addition, pregnant women without extreme/missing information of blood pressure and who measured plasma 25(OH)D levels in the first, second, or third trimester were included in the study. In addition, because gestational hypertension (GH), preeclampsia, and eclampsia are different from pregnancy complicated with chronic hypertension and chronic hypertension complicated with preeclampsia in pathogenesis and clinical treatment, pregnant women with chronic hypertension before pregnancy were also excluded [1]. Informed consent was obtained from all participants before the investigation. The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of Zhoushan Maternal and Child Health Care Hospital on 9 January 2011 (Ethical Approval Code: 2011-05).

2.2. Collection of Data and Blood Sample

The interviewers conducted face-to-face questionnaire surveys. Socio-demographic characteristics, lifestyle, and health behavior in the first (T1: 8th–14th gestational week), second (T2: 24th–28th gestational week), and third (T3: 32nd–36th gestational week) trimester, and 42nd day postpartum were collected. At each visit, professional nurses and inspectors were responsible for drawing and centrifuging fasting venous blood samples under 4 °C and separating the plasma and white blood cells, which were then stored at −80 °C until use.

2.3. Measurement of 25(OH)D Concentrations

Plasma 25(OH)D2 and 25(OH)D3 concentrations (reported in ng/mL) were measured by Liquid chromatography–tandem mass spectrometry (API 3200MD (Applied Bio-systems/MDS Sciex, Framingham, MA, USA)). The lowest sensitivity of 25(OH)D2 and 25(OH)D3 was 2 ng/mL and 5 ng/mL, respectively. The intra-assay and inter-assay coefficient variance were 1.47–7.24% and 4.48–6.74% for 25(OH)D2 and 2.50–7.59% and 4.44–6.76% for 25(OH)D3, respectively [19]. The 25(OH)D concentrations were the sum of 25(OH)D2 and 25(OH)D3.

2.4. Data Extraction

According to the guidelines of pregnant women prenatal health care, the first check-up and registration was conducted on the 8th–12th gestational week. After registration, 12 check-ups at 16, 20, 24, 28, 30, 32, 34, 36, 37, 38, 39, and 40 weeks of pregnancy were followed, along with a birth check every three days until delivery was performed after the 40th week, and a postpartum visit on the 42th day after delivery. The information including height, gestational age, and follow-up information (e.g., weight, systolic blood pressure (SBP), diastolic blood pressure (DBP), etc.), socio-demographic characteristics (e.g., age, education level, etc.), reproductive history (e.g., gravidity, parity, threatened abortion, and fetal malformation, etc.), history of present diseases (e.g., diabetes, etc.), pregnant complications (such as gestational diabetes mellitus, preeclampsia, and kidney disease, etc.), intrapartum complications (e.g., fetal distress, placenta previa, and placental abruption, etc.), was extracted from an electronic medical recorder system (EMRS).

2.5. Covariates Assessment

According to Endocrine Society Clinical Practice Guidelines, we defined plasma 25(OH)D < 20 ng/mL (50 nmol/L) as VitD deficiency [20], and 25(OH)D concentrations ≥ 20 ng/mL as VitD non-deficiency. The change of 25(OH)D level during pregnancy is defined as a difference of 25(OH)D level between three trimesters. The following parameters were also defined: Pre-pregnancy body mass index (BMI) = weight (kg)/height2 (m2), gestational weight gain (continuous) = the weight on the day of VitD test at T1, T2, or T3, the pre-pregnancy weight, educational level (senior high school and below, college and above), gravity (1, ≥2, missing), parity (0, ≥1, missing), basal blood pressure (the level of blood pressure at the first prenatal examination or early pregnancy, continuous), the seasons of blood pressure measurement (divided as followed: spring (March to May), summer (June to August), fall (September to November), and winter (December to February) based on the sunshine intensity and duration in different months [21]).

2.6. HDP Definition

In perinatal care, SBP and DBP would be routinely measured [22]; we extracted the data from EMRS. In a sitting position, blood pressure measurement was performed from the right hand with a standard mercury sphygmomanometer. GH onset was defined as SBP ≥ 140 and/or DBP ≥ 90 mm Hg after the 20th gestational week (according to last menstruation date and B-ultrasound) in at least two consecutive examinations [23]. On the basis of GH, urinary protein ≥ +1 on a dipstick was defined as preeclampsia [1]. Eclampsia was defined as the presence of new-onset grand mal seizures in a woman with preeclampsia [24]. GH, preeclampsia, and eclampsia were combined as the group of HDP in later analysis.

2.7. SNP Selection and Genotyping

VitD-related SNP were selected if they met any one of the following conditions [25,26]: (1) SNPs positively associated with 25(OH)D concentration reported in the literature, and the minimum allele frequency (MAF) ≥ 10%; (2) SNPs displayed in the functional region in the NCBI database: exon region, intron splicing point, 5′end and 3′end regulatory regions, and MAF ≥ 10%; (3) HapMap Chinese database, including gene regions, SNPs within 1500 bp at the 5′end and 3′end; (4) selected by HaploView, the conditions are: MAF ≥ 10%; R2 ≥ 0.8. Finally, a total of 34 SNPs in the VitD metabolic pathway were selected (CYP27B1: rs10877012, CYP3A4: rs2242480, rs4646437, LRP2: rs4667591, rs10210408, rs2228171, rs7600336, rs2544381, rs2544390, rs2389557, GC: rs16846876, rs12512631, rs17467825, rs2070741, rs2282679, rs3755967, rs2298850, rs4588, rs7041, rs222020, rs1155563, rs2298849, VDR: rs2228570, rs7975232, rs11568820, rs2238136, rs2853559, rs4334089, rs10783219, CYP24A1: rs6013897, rs2762934, rs2209314, rs6127118, rs2248137).
The conventional phenol–chloroform extraction method was used to extract DNA from the peripheral blood leukocytes, which was then stored in TE-buffer at −80 °C. DNA was diluted to 10 ng/μL using a Nanodrop® ND-1000 Spectrophotometer (Thermo Fisher Scientific Inc., Wilmington, NC, USA) for SNP analysis. A Sequenom MassARRAY iPLEX Gold platform (Sequenom, San Diego, CA, USA) was used for SNP genotyping. The call rate of these SNPs was over 98%, which conformed to the Hardy–Weinberg equilibrium.

2.8. Statistical Analysis

The characteristics between HDP and non-HDP groups were compared by t-test for continuous variables and by chi-squared test for categorical variables. Latent mixture modeling (PROC TRAJ) was used to identify subgroups that shared similar VitD patterns. Model fit was assessed using the Bayesian Information Criterion. We initiated a model with three trajectories, and then compared the BIC to that with two. The model with three trajectories identified fit best [27] (Figure S1). Restricted cubic spline (RCS) analyses were used to characterize the dose-response association and explore the potential linear or nonlinear relationship of 25(OH)D level in three trimesters, the change of 25(OH)D level during pregnancy with blood pressure in three trimesters, and HDP. Multivariable adjusted analyses with three knots were used. Test result for nonlinearity was checked first. If the test for nonlinearity was not significant, test result for overall association and linearity was checked, with a significant result indicating a linear association [28]. Multivariate adjusted RCS analysis showed that there was no nonlinear association of 25(OH)D level in three trimesters, the change of 25(OH)D level during pregnancy with blood pressure, and HDP during pregnancy (Pfor non-linear > 0.05) (Figures S2–S6). The Hardy–Weinberg equilibrium (HWE) of genotyped SNPs was tested using the χ2 test.
A multiple linear regression model and a multivariate logistic regression model, combined with a crossover analysis method were utilized to explore the association between VitD and its metabolic pathway-related gene variants as well as their interactions with SBP, DBP, and HDP. The generalized linear model was used to analyze the relationship of the change of 25(OH)D level during pregnancy with SBP and DBP, and the multivariate logistic regression model was used to analyze the association between the changes in 25(OH)D levels and the trajectory of VitD during pregnancy with HDP. Models were adjusted for the following potential confounders: pre-pregnancy BMI, maternal age, gestational weight gain, gestational week, educational level, parity, basal blood pressure, and the seasons of blood pressure measurement.
β (se) for linear regression, ORs, and corresponding 95% CIs for logistic regression were calculated, respectively. All test results were considered statistically significant at a value of p < 0.05. RCS analyses were performed using R software (version 3.6.3); the other analyses were performed using SAS (version 9.4, SAS Institute, Cary, NC, USA).

3. Results

3.1. Subject Characteristics

The demographic characteristics of participants with HDP or non-HDP were compared and are shown in Table 1. The prospective cohort study included 3699 pregnant women, of which 105 (2.85%) were diagnosed with HDP. The mean age was 29.30 ± 3.95 years for HDP participants and 28.67 ± 3.64 years for non-HDP participants. Compared with non-HDP participants, HDP women had higher pre-pregnancy BMI (21.16 ± 2.91 kg/m2 vs. 23.62 ± 4.05 kg/m2, p < 0.0001). The SBP and DBP levels in three trimesters were higher in HDP than non-HDP. However, VitD deficiency in three trimesters, educational level, gravity, and parity were not significantly different between the two groups. The characteristics of participants in the SNP analysis are shown in Supplementary Table S1. Pregnant women with HDP had higher pre-pregnancy BMI than the non-HDP group. There was no significant difference in weight gain and 25(OH)D level in three trimesters, educational level, gravity, and parity between the two groups (Table S1).

3.2. The Association between 25(OH)D in Three Trimesters and HDP

After being adjusted for potential confounders, 25(OH)D level at T1 was negatively associated with SBP (β (se) = −0.05 (0.02), p = 0.0287) and DBP (β (se) = −0.05 (0.02), p = 0.0190) at T1. In addition, 25(OH)D level at T2 was negatively associated with DBP at T2 and T3, respectively (β (se) = −0.10 (0.02), p < 0.0001, β (se) = −0.07 (0.02), p = 0.0003) (Table 2). The association between VitD deficiency in three trimesters with SBP and DBP were consistent with the above results (Table S2). For each 1 ng/mL increase in 25(OH)D changes between T1 and T2, DBP at T2 and T3 decreased by 0.11 (se = 0.02) mmHg and 0.11 (se = 0.02) mmHg, respectively (p < 0.0001). (Table 3). Three subgroups of participants with data of 25(OH)D levels in three trimesters were identified by latent mixture modeling. Compared with women whose 25(OH)D levels remained low from T1 to T3, women whose 25(OH)D levels gradually increased at T2 and T3 or whose 25(OH)D levels remained high during pregnancy had lower DBP at T3 (β (se) = −1.13 (0.46), p = 0.0137, β (se) = −1.74 (0.74), p = 0.0195) (Table 4). However, there was no significant association between 25(OH)D levels, VitD deficiency in three trimesters, the change of 25(OH)D levels, or the VitD trajectory during pregnancy with HDP (Tables S3–S6).

3.3. The Association between SNP and HDP

The association of each SNP genotype with SBP and DBP at T1, T2, and T3 are shown in Tables S7–S9, respectively. Polymorphisms in CYP24A1-rs2248137 was significantly associated with higher SBP at T1 and DBP at T2 and T3. Polymorphisms in CYP24A1-rs2762934 were significantly associated with higher DBP at T1 and SBP at T2. Polymorphisms in LRP2-rs4667591 were significantly associated with higher SBP at T1 and DBP at T3. Polymorphisms in GC-rs2070741, rs222020, and rs2298849 were associated with higher SBP at T2. Polymorphisms in LRP2-rs2544390 were associated with higher DBP at T3. Furthermore, polymorphisms in CYP24A1-rs2248137, CYP24A1-rs2762934, CYP24A1-rs6127118, and GC-rs2070741 were associated a higher risk of HDP (Table 5). However, there was no significant association between other genes’ polymorphisms and HDP.

3.4. The Interaction between Single SNP and VitD Deficiency in Three Trimesters on the Risk of HDP

Results of the crossover analysis are shown in Tables S10–S13. Polymorphisms of seven SNPs (rs16846876, rs2282679, rs17467825, rs2298849, rs2298850, rs3755967, and rs4588) in GC gene and VitD deficiency at T2 might exert interactions on DBP at T2. In addition, VDR-rs2228570 and VitD deficiency at T2 might exert interaction on SBP at T2. Furthermore, women with mutations in CYP24A1-rs2248137, LRP2-rs2389557, and LRP2-rs4667591 and had VitD deficiency at T2 showed increased risk of HDP (Table 6).

4. Discussion

In the present study, 25(OH)D level at T2, as well as 25(OH)D change between T1 and T2, were significantly inversely associated with DBP at T2 and T3. However, significant associations between maternal VitD deficiency in any trimesters and HDP were not observed. Polymorphism in CYP24A1, GC, and LRP2 was associated with blood pressure, and polymorphism in CYP24A1 and GC was associated with increased risk of HDP. Furthermore, interactive effects between VitD deficiency and polymorphisms in CYP24A1, GC, and VDR genes on blood pressure were identified. Women with polymorphisms in CYP24A1 and LRP2 genes and had VitD deficiency at T2 showed a higher risk of HDP.
Previous findings on the association between VitD level during pregnancy and HDP were not consistent. A prospective observational study conducted in southern China found that there were no significant differences in the risk of HDP among women with different levels of VitD at 16–20-week gestation [29]. A case-control study conducted in Iran found that pregnant women with VitD deficiency (25(OH)D < 20 ng/mL) had higher blood pressure and increased risk of preeclampsia than those with VitD insufficiency (25(OH)D: 20~30 ng/mL) [8]. The prospective Swedish GraviD cohort study, including 1413 pregnant women, found that 25(OH)D was positively associated with T1 blood pressure [16]; however, both 25(OH)D level at T3 and change in 25(OH)D level from T1 to T3 were significantly and negatively associated with preeclampsia, but not with the risk of GH [30]. Another nested case-control study carried out among Australian pregnant women found that higher levels of VitD (25(OH)D > 75 nmol/L) in early pregnancy (10–14 weeks) could prevent the occurrence of early-onset preeclampsia (p = 0.09); however, women with low levels of 25(OH)D (<37.5 ng/mL) in the first trimester of pregnancy had a tendency toward reduced risk of preeclampsia (p = 0.07) [31]. Conflicting data for an association of VitD during pregnancy with HDP results from a number of sources, including large heterogeneity between study designs, different ethnicities, different subtypes of HDP included in the analysis, variable quality of measurement for 25(OH)D, and inconsistent definition of VitD status [32]. On the other hand, studies have shown that the gene variation of key enzymes in VitD synthesis, transport and metabolism pathway would also affect the levels and effects of 25(OH)D and 1,25(OH)2D [25,33]. Furthermore, genetic mutations in the VitD metabolic pathway were also associated with increased risk of HDP [8].
The active form of VitD (1,25(OH)2D) needs to bind to VDR to exert its biological function. Relevant studies related to genetic variants in the VitD metabolic pathway with HDP were mainly focus on three SNPs (rs2228570, rs731236, and rs1544410) of VDR gene. Rezavand et al. [8] found that, compared with VDR-rs2228570 TC and TT + TC genotypes, the SBP and DBP of CC genotype were higher, and the risk of preeclampsia increased by 1.72 times. However, no association was found between VDR-rs731236, VDR-rs1544410, and preeclampsia. Knabl et al. [34] also reported that there was a strong association between the polymorphisms in rs10735810 and rs1544410 of VDR and the risk of GH. The polymorphisms in rs10735810 affect plasma renin activity and may be associated with a reduced risk of GH [34]. In this study, VitD deficiency at T2 interacted with the variants of VDR-rs2238136 on DBP and VDR-rs2228570 on SBP at T2.
The CYP24A1 gene is located in 20q13-2, which is mainly expressed in the kidney and encodes the catabolic enzymes of 1,25(OH)2D and 25(OH)D [35]. Evidence relating to the association between CYP24A1 gene polymorphism and susceptibility to hypertension, especially among pregnant women, is scare. A case-control study among the Chinese Han population found that CYP24A1-rs56229249 significantly decreased the hypertension risk in homozygote and recessive models [36]. In addition, rs2762940 was related to hypertension risk in men, and rs56229249 was a protective factor against hypertension in women [36]. The comprehensive genetic association study in the Women’s Genome Health Study (WGHS) found that CYP24A1-rs2296241 showed significant associations with SBP, DBP, mean arterial pressure, and pulse pressure [37]. In this study, we found that gene variants in CYP24A1-rs2248137, CYP24A1-rs2762934, and CYP24A1-rs6127118 were associated with increased risk of HDP. Furthermore, CYP24A1-rs6013897 interacted with VitD deficiency at T2 on HDP. On the other hand, LRP2 is located on 2q24-q31, which is a member of the low-density lipoprotein receptor family and encodes megalin protein. In the kidney, megalin and cubilin combine together with hydroxylate 25(OH)D3 into 1,25(OH)2D3 [38]. Studies regarding the association between LRP2 genes and VitD with the risk for HDP are still lacking. This study found that the mutations of LRP2-rs2389557 and LRP2-rs4667591 and VitD deficiency at T2 had a combined effect on the risk of HDP.
The GC gene encodes VitD binding protein (VDBP) [39], which is the major transporter of VitD. About 85% to 90% of 25(OH)D is bound to VDBP in circulation [40]. VDBP can aggravate or enhance various biological processes during pregnancy, such as immune regulation, glucose metabolism, and blood pressure regulation [39]. The GC-1 subtype was more common in pregnant women with preeclampsia than in those without preeclampsia, which was considered as a potential early detection genetic marker for women at risk of preeclampsia [41]. In HIV endemic areas of South Africa, compared with women with normal blood pressure, two SNPs of GC gene (rs4588 and rs7041) were more common in pregnant women with preeclampsia, and were not related to HIV status [42]. Furthermore, GC-rs4588 polymorphism was associated with early-onset (<34 weeks) and late-onset (≥34 weeks of pregnancy) preeclampsia, while GC-rs7041 was associated with early-onset eclampsia [42]. A nested case-control study of 170 American women from Massachusetts tracked the levels of VDBP and 25(OH)D throughout pregnancy to examine whether these biomarkers were associated with blood pressure or the risk of preeclampsia, but found no significant correlation of VDBP or 25(OH)D levels with preeclampsia [43]. At present, the combined effect of GC gene polymorphism and VitD during pregnancy on HDP is not clear. A study focused on preterm birth found that rs7041 variants interacted with VitD at T2 on the gestational week of delivery and preterm birth [44]. Our study found that the variant of GC-rs2070741 was associated with higher SBP at T2 and increased risk of HDP. Mutations at GC rs16846876, rs2282679, rs17467825, rs2298849, rs2298850, rs3755967, and rs4588 interacted with VitD deficiency at T2 on higher DBP at T2.
To our knowledge, this is the first prospective cohort study exploring the association between VitD in three trimesters and VitD pathway gene variants as well as their interactions on SBP, DBP, and the risk of HDP. However, limitations could not be neglected. First of all, some subjects had a lack of 25(OH)D data at T2 and T3, and therefore selection bias might exist. However, subgroup analysis of pregnant women with VitD detected at T1 and T2 showed that the results were almost consistent with the results in the whole study population. Secondly, as the prevalence of HDP in this study was relatively low (2.84%), the association between VitD and different HDP subtypes (GH, preeclampsia, eclampsia) could not be explored. However, studies have shown that, although these subtypes can appear alone, they are progressive manifestations of a single process and share common etiology [45,46]. Lastly, the relatively single ethnic population of this study may also limit the extrapolation of findings.

5. Conclusions

This study found that the level of 25(OH)D at T1 and T2 was negatively correlated with DBP at T2. In addition, polymorphisms in VitD metabolic pathway genes, including CYP24A1 and GC, increased the risk of HDP. Furthermore, gene variants in CYP24A1 and LRP2 and VitD deficiency at T2 showed combined effect on the risk of HDP, but the specific mechanism remains to be further investigated. The results of this study provide a scientific basis for the clinical detection of VitD during pregnancy and the supplementation of VitD during pregnancy.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/nu14112355/s1, Figure S1. Trajectory of 25(OH)D level during pregnancy. Figure S2. Dose-response relationships of 25(OH)D levels at T1 and T2 with blood pressure at T1 and T2. Figure S3. Dose-response relationships of 25(OH)D levels at T1, T2, and T3 with blood pressure at T3. Figure S4. Dose-response relationships of 25(OH)D levels at T1, T2, and T3 with HDP. Figure S5. Dose-response relationships of the change of 25(OH)D levels during pregnancy with blood pressure at T2 and T3. Figure S6. Dose-response relationships between the change of 25(OH)D levels during pregnancy and HDP. Table S1. Baseline characteristics of pregnant women in SNP analysis. Table S2. The relationship of VitD deficiency in three trimesters with blood pressure. Table S3. The association between 25(OH)D levels in three trimesters and HDP. Table S4. The relationship between VitD deficiency in three trimesters with HDP. Table S5. The association between the change of 25(OH)D levels during pregnancy and HDP. Table S6. The association between the trajectory of VitD during pregnancy and HDP. Table S7. The association of single SNP with SBP and DBP at T1. Table S8. The association of single SNP with SBP and DBP at T2. Table S9. The association of single SNP with SBP and DBP at T3. Table S10. The association of single SNP and VitD at T1 with blood pressure at T1. Table S11. The association between single SNP and VitD at T2 with blood pressure at T2. Table S12. The association between single SNP and VitD at T3 with blood pressure at T3. Table S13. The association between single SNP and VitD at T2 with blood pressure at T3.

Author Contributions

Conceptualization, Y.Y.; methodology, M.M.; validation, S.S.; resources, H.C.; formal analysis, Z.P.; investigation, X.A.; writing—review and editing, M.M.; visualization, H.Z., P.C. and Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Chinese National Natural Science Foundation (81973055), the National Key Research and Development Programme of China (No. 2021YFC2701901), and Major research and development projects of the Zhejiang Science and Technology Department (2018C03010), the Key Laboratory of Intelligent Preventive Medicine of Zhejiang Province (2020E10004) and the Leading Innovative and Entrepreneur Team Introduction Program of Zhejiang (2019R01007).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of Zhoushan Maternal and Child Health Care Hospital on 9 January 2011 (Ethical Approval Code: 2011-05).

Informed Consent Statement

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

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available because they contain information that could compromise the privacy of research participants.

Acknowledgments

We thank all the participants who took part in this study. We acknowledge the support of staff in Zhoushan Maternal and Child Care Hospital, who conducted and managed the cohort.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. American College of Obstetricians and Gynecologists’ Task Force on Hypertension in Pregnancy. Hypertension in Pregnancy. Obstet. Gynecol. 2013, 122, 1122–1131. [Google Scholar] [CrossRef]
  2. Khan, K.S.; Wojdyla, D.; Say, L.; Gülmezoglu, A.M.; Van Look, P.F. WHO analysis of causes of maternal death: A systematic review. Lancet 2006, 367, 1066–1074. [Google Scholar] [CrossRef]
  3. Say, L.; Chou, D.; Gemmill, A.; Tunçalp, Ö.; Moller, A.-B.; Daniels, J.; Gülmezoglu, A.M.; Temmerman, M.; Alkema, L. Global causes of maternal death: A WHO systematic analysis. Lancet Glob. Health 2014, 2, E323–E333. [Google Scholar] [CrossRef] [Green Version]
  4. Theilen, L.H.; Meeks, H.; Fraser, A.; Esplin, M.S.; Smith, K.R.; Varner, M.W. Long-term mortality risk and life expectancy following recurrent hypertensive disease of pregnancy. Am. J. Obstet. Gynecol. 2018, 219, 107.e1–107.e6. [Google Scholar] [CrossRef] [Green Version]
  5. Visintin, C.; Mugglestone, M.A.; Almerie, M.Q.; Nherera, L.M.; James, D.; Walkinshaw, S. On behalf of the Guideline Development Group Management of hypertensive disorders during pregnancy: Summary of NICE guidance. BMJ 2010, 341, c2207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Rana, S.; Lemoine, E.; Granger, J.P.; Karumanchi, S.A. Preeclampsia: Pathophysiology, Challenges, and Perspectives. Circ. Res. 2019, 124, 1094–1112. [Google Scholar] [CrossRef] [PubMed]
  7. Adela, R.; Borkar, R.M.; Mishra, N.; Bhandi, M.M.; Vishwakarma, G.; Varma, B.A.; Ragampeta, S.; Banerjee, S.K. Lower Serum Vitamin D Metabolite Levels in Relation to Circulating Cytokines/Chemokines and Metabolic Hormones in Pregnant Women with Hypertensive Disorders. Front. Immunol. 2017, 8, 273. [Google Scholar] [CrossRef] [Green Version]
  8. Rezavand, N.; Tabarok, S.; Rahimi, Z.; Vaisi-Raygani, A.; Mohammadi, E.; Rahimi, Z. The effect of VDR gene polymorphisms and vitamin D level on blood pressure, risk of preeclampsia, gestational age, and body mass index. J. Cell. Biochem. 2018, 120, 6441–6448. [Google Scholar] [CrossRef]
  9. Zhang, J.Y.; Lucey, A.J.; Horgan, R.; Kenny, L.C.; Kiely, M. Impact of pregnancy on vitamin D status: A longitudinal study. Br. J. Nutr. 2014, 112, 1081–1087. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Ganguly, A.; Tamblyn, J.A.; Finn-Sell, S.; Chan, S.; Westwood, M.; Gupta, J.; Kilby, M.D.; Gross, S.; Hewison, M. Vitamin D, the placenta and early pregnancy: Effects on trophoblast function. J. Endocrinol. 2018, 236, R93–R103. [Google Scholar] [CrossRef]
  11. Weisman, Y.; Harell, A.; Edelstein, S.; David, M.; Spirer, Z.; Golander, A. 1α, 25-Dihydroxyvitamin D3, and 24,25-dihydroxyvitamin D3 in vitro synthesis by human decidua and placenta. Nature 1979, 281, 317–319. [Google Scholar] [CrossRef]
  12. Zehnder, D.; Evans, K.N.; Kilby, M.D.; Bulmer, J.N.; Innes, B.A.; Stewart, P.M.; Hewison, M. The Ontogeny of 25-Hydroxyvitamin D3 1α-Hydroxylase Expression in Human Placenta and Decidua. Am. J. Pathol. 2002, 161, 105–114. [Google Scholar] [CrossRef]
  13. Karras, S.N.; Wagner, C.L.; Castracane, V.D. Understanding vitamin D metabolism in pregnancy: From physiology to pathophysiology and clinical outcomes. Metabolism 2018, 86, 112–123. [Google Scholar] [CrossRef] [PubMed]
  14. De-Regil, L.M.; Palacios, C.; Lombardo, L.K.; Pena-Rosas, J.P. Vitamin D supplementation for women during pregnancy. Cochrane Database Syst. Rev. 2016, 1, CD008873. [Google Scholar]
  15. Palacios, C.; Kostiuk, L.K.; Pena-Rosas, J.P. Vitamin D supplementation for women during pregnancy. Cochrane Database Syst. Rev. 2019, 7, CD008873. [Google Scholar] [CrossRef]
  16. Barebring, L.; O’Connell, M.; Winkvist, A.; Johannsson, G.; Augustin, H. Serum cortisol and vitamin D status are independently associated with blood pressure in pregnancy. J. Steroid Biochem. Mol. Biol. 2019, 189, 259–264. [Google Scholar] [CrossRef]
  17. Magnus, M.C.; Miliku, K.; Bauer, A.; Engel, S.M.; Felix, J.F.; Jaddoe, V.W.V.; ALawlor, D.; London, S.J.; Magnus, P.; McGinnis, R.; et al. Vitamin D and risk of pregnancy related hypertensive disorders: Mendelian randomisation study. BMJ 2018, 361, k2167. [Google Scholar] [CrossRef] [Green Version]
  18. Barry, E.L.; Rees, J.R.; Peacock, J.L.; Mott, L.A.; Amos, C.I.; Bostick, R.M.; Figueiredo, J.C.; Ahnen, D.J.; Bresalier, R.S.; Burke, C.A.; et al. Genetic variants in CYP2R1, CYP24A1, and VDR modify the efficacy of vitamin D3 supplementation for increasing serum 25-hydroxyvitamin D levels in a randomized controlled trial. J. Clin. Endocrinol. Metab. 2014, 99, E2133–E2137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Shao, B.; Mo, M.; Xin, X.; Jiang, W.; Wu, J.; Huang, M.; Wang, S.; Muyiduli, X.; Si, S.; Shen, Y.; et al. The interaction between prepregnancy BMI and gestational vitamin D deficiency on the risk of gestational diabetes mellitus subtypes with elevated fasting blood glucose. Clin. Nutr. 2020, 39, 2265–2273. [Google Scholar] [CrossRef]
  20. Holick, M.F.; Binkley, N.C.; Bischoff-Ferrari, H.A.; Gordon, C.M.; Hanley, D.A.; Heaney, R.P.; Murad, M.H.; Weaver, C.M. Endocrine Society. Evaluation, Treatment, and Prevention of Vitamin D Deficiency: An Endocrine Society Clinical Practice Guideline. J. Clin. Endocrinol. Metab. 2011, 96, 1911–1930. [Google Scholar] [CrossRef] [Green Version]
  21. Thorne, H.C.; Jones, K.H.; Peters, S.P.; Archer, S.N.; Dijk, D.-J. Daily and Seasonal Variation in the Spectral Composition of Light Exposure in Humans. Chrono Int. 2009, 26, 854–866. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Zhao, P.; Diao, Y.; You, L.; Wu, S.; Yang, L.; Liu, Y. The influence of basic public health service project on maternal health services: An interrupted time series study. BMC Public Health 2019, 19, 824–831. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Brown, M.A.; Magee, L.A.; Kenny, L.C.; Karumanchi, S.A.; McCarthy, F.P.; Saito, S.; Hall, D.R.; Warren, C.E.; Adoyi, G.; Ishaku, S.; et al. Hypertensive Disorders of Pregnancy: ISSHP Classification, Diagnosis, and Management Recommendations for International Practice. Hypertension 2018, 72, 24–43. [Google Scholar] [CrossRef] [Green Version]
  24. American College of Obstetricians and Gynecologists. Practice ACoO. ACOG practice bulletin. Diagnosis and management of preeclampsia and eclampsia. International journal of gynaecology and obstetrics: The official organ of the International Federation of Gynaecology and Obstetrics. Obstet. Gynecol. 2002, 77, 67–75. [Google Scholar]
  25. Shao, B.; Jiang, S.; Muyiduli, X.; Wang, S.; Mo, M.; Li, M.; Wang, Z.; Yu, Y. Vitamin D pathway gene polymorphisms influenced vitamin D level among pregnant women. Clin. Nutr. 2018, 37, 2230–2237. [Google Scholar] [CrossRef]
  26. Wu, J.; Shao, B.; Xin, X.; Luo, W.; Mo, M.; Jiang, W.; Si, S.; Wang, S.; Shen, Y.; Yu, Y. Association of vitamin D pathway gene polymorphisms with vitamin D level during pregnancy was modified by season and vitamin D supplement. Clin. Nutr. 2020, 40, 3650–3660. [Google Scholar] [CrossRef]
  27. Li, W.; Jin, C.; Vaidya, A.; Wu, Y.; Rexrode, K.; Zheng, X.; Gurol, M.E.; Ma, C.; Wu, S.; Gao, X. Blood Pressure Trajectories and the Risk of Intracerebral Hemorrhage and Cerebral Infarction: A Prospective Study. Hypertension 2017, 70, 508–514. [Google Scholar] [CrossRef] [PubMed]
  28. Desquilbet, L.; Mariotti, F. Dose-response analyses using restricted cubic spline functions in public health research. Stat. Med. 2010, 29, 1037–1057. [Google Scholar] [CrossRef] [PubMed]
  29. Zhou, J.; Su, L.; Liu, M.; Liu, Y.; Cao, X.; Wang, Z.; Xiao, H. Associations between 25-hydroxyvitamin D levels and pregnancy outcomes: A prospective observational study in southern China. Eur. J. Clin. Nutr. 2014, 68, 925–930. [Google Scholar] [CrossRef] [Green Version]
  30. Bärebring, L.; Bullarbo, M.; Glantz, A.; Agelii, M.L.; Jagner; Ellis, J.; Hulthen, L.; Schoenmakers, I.; Augustin, H. Preeclampsia and Blood Pressure Trajectory during Pregnancy in Relation to Vitamin D Status. PLoS ONE 2016, 11, e0152198. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Schneuer, F.J.; Roberts, C.L.; Guilbert, C.; Simpson, J.M.; Algert, C.S.; Khambalia, A.Z.; Tasevski, V.; Ashton, A.W.; Morris, J.M.; Nassar, N. Effects of maternal serum 25-hydroxyvitamin D concentrations in the first trimester on subsequent pregnancy outcomes in an Australian population. Am. J. Clin. Nutr. 2014, 99, 287–295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. O’Callaghan, K.M.; Kiely, M. Systematic Review of Vitamin D and Hypertensive Disorders of Pregnancy. Nutrients 2018, 10, 294. [Google Scholar] [CrossRef] [Green Version]
  33. Wang, T.J.; Zhang, F.; Richards, J.B.; Kestenbaum, B.; van Meurs, J.B.; Berry, D.; Kiel, D.P.; Streeten, E.A.; Ohlsson, C.; Koller, D.L.; et al. Common genetic determinants of vitamin D insufficiency: A genome-wide association study. Lancet 2010, 376, 180–188. [Google Scholar] [CrossRef] [Green Version]
  34. Knabl, J.; Vattai, A.; Ye, Y.; Jueckstock, J.; Hutter, S.; Kainer, F.; Mahner, S.; Jeschke, U. Role of Placental VDR Expression and Function in Common Late Pregnancy Disorders. Int. J. Mol. Sci. 2017, 18, 2340. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Sakaki, T.; Kagawa, N.; Yamamoto, K.; Inouye, K. Metabolism of vitamin D3 by cytochromes P450. Front. Biosci. 2005, 10, 119–134. [Google Scholar]
  36. Bao, Q.; Wang, D.; Zhang, Y.; Bao, L.; Jia, H. The Impact of CYP24A1 Polymorphisms on Hypertension Susceptibility. Kidney Blood Press. Res. 2020, 45, 28–37. [Google Scholar] [CrossRef]
  37. Wang, L.; Chu, A.; Buring, J.E.; Ridker, P.M.; Chasman, D.I.; Sesso, H.D. Common Genetic Variations in the Vitamin D Pathway in Relation to Blood Pressure. Am. J. Hypertens. 2014, 27, 1387–1395. [Google Scholar] [CrossRef] [Green Version]
  38. Marzolo, M.-P.; Farfán, P. New Insights into the Roles of Megalin/LRP2 and the Regulation of its Functional Expression. Biol. Res. 2011, 44, 89–105. [Google Scholar] [CrossRef] [Green Version]
  39. Norman, A.W. 1α,25(OH)2 Vitamin D3Vitamin D Nuclear Receptor (VDR) and Plasma Vitamin D-Binding Protein (DBP) Structures and Ligand Shape Preferences for Genomic and Rapid Biological Responses; Elsevier Inc.: Amsterdam, The Netherlands, 2008. [Google Scholar]
  40. Bikle, D.D.; Malmstroem, S.; Schwartz, J. Current Controversies: Are Free Vitamin Metabolite Levels a More Accurate Assessment of Vitamin D Status than Total Levels? Endocrinol. Metab. Clin. N. Am. 2017, 46, 901–918. [Google Scholar] [CrossRef]
  41. Mekbeb, T. The association of serum proteins with preeclampsia. Ethiop. Med. J. 1990, 28, 1–30. [Google Scholar]
  42. Naidoo, Y.; Moodley, J.; Ramsuran, V.; Naicker, T. Polymorphisms within vitamin D binding protein gene within a Preeclamptic South African population. Hypertens. Pregnancy 2019, 38, 260–267. [Google Scholar] [CrossRef]
  43. Powe, C.E.; Seely, E.W.; Rana, S.; Bhan, I.; Ecker, J.; Karumanchi, S.A.; Thadhani, R. First trimester vitamin D, vitamin D binding protein, and subsequent preeclampsia. Hypertension 2010, 56, 758–763. [Google Scholar] [CrossRef] [Green Version]
  44. Wang, S.; Xin, X.; Luo, W.; Mo, M.; Si, S.; Shao, B.; Shen, Y.; Cheng, H.; Yu, Y. Association of vitamin D and gene variants in the vitamin D metabolic pathway with preterm birth. Nutrition 2021, 89, 111349. [Google Scholar] [CrossRef] [PubMed]
  45. Phoswa, W.N. Dopamine in the Pathophysiology of Preeclampsia and Gestational Hypertension: Monoamine Oxidase (MAO) and Catechol-O-methyl Transferase (COMT) as Possible Mechanisms. Oxidative Med. Cell. Longev. 2019, 2019, 1–8. [Google Scholar] [CrossRef] [PubMed]
  46. Mammaro, A.; Carrara, S.; Cavaliere, A.; Ermito, S.; Dinatale, A.; Pappalardo, E.M.; Militello, M.; Pedata, R. Hypertensive disorders of pregnancy. J. Prenat. Med. 2009, 3, 1–5. [Google Scholar] [PubMed]
Table 1. Baseline characteristics of pregnant women.
Table 1. Baseline characteristics of pregnant women.
VariablesNon-HDP (N = 3594)HDP (N = 105)p
Mean ± SD
Age, years28.67 ± 3.6429.30 ± 3.950.0811
Pre-pregnancy BMI, kg/m221.16 ± 2.9123.62 ± 4.05<0.0001
T1 (N = 3302)
Weight gain, kg 0.01 ± 0.170.02 ± 0.140.3841
SBP, mmHg103.53 ± 9.31112.21 ± 10.58<0.0001
DBP, mmHg68.35 ± 6.6774.48 ± 5.99<0.0001
25(OH)D, ng/mL 17.85 ± 8.3817.70 ± 7.100.8629
T2 (N = 2479)
Weight gain, kg 5.61 ± 3.826.25 ± 4.500.1971
SBP, mmHg107.17 ± 9.24116.60 ± 14.44<0.0001
DBP, mmHg69.13 ± 7.8276.61 ± 8.84<0.0001
25(OH)D, ng/mL 23.28 ± 10.3822.91 ± 9.600.7827
T3 (N = 1549)
Weight gain, kg 11.91 ± 3.7311.20 ± 4.690.2181
SBP, mmHg108.86 ± 9.72123.30 ± 15.68<0.0001
DBP, mmHg70.99 ± 7.3882.43 ± 8.44<0.0001
25(OH)D, ng/mL 26.53 ± 11.2826.20 ± 11.010.8468
N (%)
VitD deficiency at T1 a2176 (67.85)66 (69.47)0.7386
VitD deficiency at T2 b1067 (44.15)30 (48.39)0.5067
VitD deficiency at T3 c476 (31.63)13 (29.55)0.7696
Educational level 0.1263
≤High school957 (26.63)35 (33.33)
>High school2637 (73.37)70 (66.67)
Gravity 0.5389
11652 (45.97)43 (40.95)
≥21822 (50.70)59 (56.19)
Unknown120 (3.34)3 (2.86)
Parity 0.4887
02015 (56.07)65 (61.90)
≥1771 (21.45)20 (19.05)
Unknown808 (22.48)20 (19.05)
Abbreviations: HDP, hypertensive disorders in pregnancy; BMI, body mass index; SBP, systolic blood pressure; DBP, diastolic blood pressure; VitD, vitamin D. a N = 3302, b N = 2479, c N = 1549.
Table 2. Association between 25(OH)D levels in three trimesters with blood pressure.
Table 2. Association between 25(OH)D levels in three trimesters with blood pressure.
Trimesters of 25(OH)DNSBP, mmHgDBP, mmHg
β (se)pβ (se)p
Blood pressure at T1 (N = 3302)
T13302−0.02 (0.02)0.31980.02 (0.01)0.2056
Blood pressure at T2 (N = 2479)
T12125−0.05 (0.02)0.0287−0.05 (0.02)0.0190
T224790.03 (0.02)0.1675−0.10 (0.02)<0.0001
Blood pressure at T3 (N = 1549)
T113280.03 (0.03)0.2361−0.02 (0.02)0.3259
T213900.04 (0.03)0.1214−0.07 (0.02)0.0003
T315490.04 (0.02)0.0541−0.02 (0.02)0.2486
Abbreviations: SBP, systolic blood pressure, DBP, diastolic blood pressure. Adjusted for pre-pregnancy BMI, maternal age, gestational weight gain, gestational week, educational level, parity, basal blood pressure, and the seasons of blood pressure measurement.
Table 3. The association between the change of 25(OH)D levels during pregnancy and blood pressure at T2 and T3.
Table 3. The association between the change of 25(OH)D levels during pregnancy and blood pressure at T2 and T3.
The Change of Trimesters NThe Change of 25(OH)D Levels, ng/mL *SBP, mmHgDBP, mmHg
β (se)pβ (se)p
Blood pressure at T2 (N = 2479)
Between T1 and T221253.50 (84.59)0.03 (0.02)0.1217−0.11 (0.02)<0.0001
Blood pressure at T3 (N = 1549)
Between T1 and T212122.40 (81.27)0.03 (0.03)0.3142−0.11 (0.02)<0.0001
Between T1 and T313286.59 (98.02)0.06 (0.03)0.0294−0.02 (0.02)0.3516
Between T2 and T313903.40 (87.23)0.02 (0.03)0.36620.04 (0.02)0.0405
Abbreviations: SBP, systolic blood pressure; DBP, diastolic blood pressure. * Presented as the median (range). Adjusted for pre-pregnancy BMI, maternal age, gestational weight gain, gestational week, educational level, parity, basal blood pressure, the seasons of blood pressure measurement, and 25(OH)D level at T1.
Table 4. The association between the trajectory of VitD during pregnancy and blood pressure at T3.
Table 4. The association between the trajectory of VitD during pregnancy and blood pressure at T3.
Trajectory of VitDN (%)SBP, mmHgDBP, mmHg
β (se)pβ (se)p
Subgroup 1621 (51.24)Ref Ref
Subgroup 2469 (38.70)0.48 (0.60)0.4216−1.13 (0.46)0.0137
Subgroup 3122 (10.07)1.58 (0.98)0.1052−1.74 (0.74)0.0195
Abbreviations: SBP, systolic blood pressure; DBP, diastolic blood pressure; VitD, vitamin D. Adjusted for pre-pregnancy BMI, maternal age, gestational weight gain, gestational week, educational level, parity, basal blood pressure, and the seasons of blood pressure measurement.
Table 5. The relationship between single SNP and HDP *.
Table 5. The relationship between single SNP and HDP *.
SNPGenotypesNCase (%)Crude ModelAdjusted Model *
OR (95%CI)pOR (95%CI)p
CYP24A1
rs2209314TT94128 (3.0)Ref Ref
CT130938 (2.9)0.97 (0.59–1.60)0.91980.99 (0.60–1.64)0.9648
CC4434 (0.9)0.30 (0.10–0.85)0.0240.30 (0.10–0.87)0.026
rs2248137GG93418 (1.9)Ref Ref
GC45320 (4.4)2.35 (1.23–4.49)0.00962.62 (1.32–5.21)0.0059
CC64319 (3.0)1.55 (0.81–2.98)0.18851.80 (0.92–3.53)0.0869
rs2762934GG59922 (3.7)Ref Ref
GA1195 (4.2)1.15 (0.43–3.10)0.78191.07 (0.38–3.00)0.9051
AA71 (14.3)4.37 (0.50–37.88)0.18069.98 (1.06–94.04)0.0444
rs6013897TT52920 (3.8)Ref Ref
AT1728 (4.7)1.24 (0.54–2.87)0.61311.26 (0.53–3.02)0.5964
AA221 (4.5)1.21 (0.16–9.46)0.85461.53 (0.19–12.57)0.6905
rs6127118GG96324 (2.5)Ref Ref
AG161641 (2.5)1.02 (0.61–1.70)0.94390.96 (0.57–1.61)0.8736
AA1197 (5.9)2.45 (1.03–5.80)0.04262.38 (0.98–5.77)0.0542
CYP27B1
rs10877012TT112527 (2.4)Ref Ref
GT120432 (2.7)1.11 (0.66–1.87)0.69251.14 (0.67–1.93)0.6354
GG36012 (3.3)1.40 (0.70–2.80)0.33661.61 (0.80–3.25)0.1856
CYP3A4
rs2242480CC152941 (2.7)Ref Ref
CT100426 (2.6)0.96 (0.59–1.59)0.88790.96 (0.58–1.60)0.8885
TT1595 (3.1)1.18 (0.46–3.03)0.73291.37 (0.53–3.56)0.5164
rs4646437GG53022 (4.2)Ref Ref
AG1827 (3.8)0.92 (0.39–2.20)0.85760.83 (0.34–2.04)0.6868
GC
rs1155563TT95129 (3.0)Ref Ref
TC129031 (2.4)0.78 (0.47–1.31)0.34990.76 (0.45–1.29)0.3091
CC45011 (2.4)0.80 (0.39–1.61)0.52620.73 (0.36–1.50)0.3886
rs12512631TT46321 (4.5)Ref Ref
CT2417 (2.9)0.63 (0.26–1.50)0.29730.77 (0.31–1.89)0.5683
CC181 (5.6)1.24 (0.16–9.75)0.83931.07 (0.13–8.84)0.9466
rs16846876AA127438 (3.0)Ref Ref
AT113325 (2.2)0.73 (0.44–1.22)0.23570.69 (0.41–1.16)0.1632
TT2929 (3.1)1.03 (0.49–2.16)0.92830.89 (0.41–1.93)0.7745
rs17467825AA125437 (3.0)Ref Ref
GA115029 (2.5)0.85 (0.52–1.39)0.52080.80 (0.48–1.32)0.3792
GG2996 (2.0)0.67 (0.28–1.61)0.3750.54 (0.22–1.33)0.181
rs2070741TT49517 (3.4)Ref Ref
GT2138 (3.8)1.10 (0.47–2.58)0.83171.10 (0.46–2.66)0.8251
GG183 (16.7)5.62 (1.49–21.28)0.0114.77 (1.12–20.21)0.0341
rs222020TT27010 (3.7)Ref Ref
CT34412 (3.5)0.94 (0.40–2.21)0.88670.85 (0.35–2.07)0.7284
CC1106 (5.5)1.50 (0.53–4.23)0.44361.42 (0.49–4.15)0.5169
rs2282679TT125436 (2.9)Ref Ref
GT114130 (2.6)0.91 (0.56–1.49)0.71850.87 (0.52–1.43)0.5726
GG3066 (2.0)0.68 (0.28–1.62)0.38080.55 (0.22–1.35)0.1913
rs2298849AA112027 (2.4)Ref Ref
GA121938 (3.1)1.30 (0.79–2.15)0.30031.32 (0.79–2.20)0.2835
GG3667 (1.9)0.79 (0.34–1.83)0.58090.82 (0.35–1.90)0.6376
rs2298850GG122935 (2.8)Ref Ref
CG115129 (2.5)0.88 (0.54–1.45)0.6210.84 (0.50–1.39)0.4891
CC3056 (2.0)0.68 (0.29–1.64)0.39670.55 (0.22–1.38)0.2051
rs3755967CC125036 (2.9)Ref Ref
CT114930 (2.6)0.90 (0.55–1.48)0.68750.86 (0.52–1.42)0.5515
TT3066 (2.0)0.67 (0.28–1.62)0.37680.54 (0.22–1.35)0.189
rs4588GG124136 (2.9)Ref Ref
GT114629 (2.5)0.87 (0.53–1.43)0.57890.83 (0.50–1.37)0.4576
TT3066 (2.0)0.67 (0.28–1.60)0.36790.54 (0.22–1.34)0.1819
rs7041AA144539 (2.7)Ref Ref
CA106128 (2.6)0.98 (0.60–1.60)0.92681.03 (0.62–1.70)0.9026
CC1954 (2.1)0.76 (0.27–2.14)0.59680.87 (0.30–2.48)0.7934
LRP2
rs10210408CC89522 (2.5)Ref Ref
TC132540 (3.0)1.24 (0.73–2.09)0.43221.16 (0.68–1.99)0.5816
TT48610 (2.1)0.83 (0.39–1.78)0.63780.81 (0.38–1.74)0.5919
rs2228171TT93224 (2.6)Ref Ref
CT39415 (3.8)1.50 (0.78–2.89)0.22791.48 (0.75–2.94)0.2607
CC3018 (2.7)1.03 (0.46–2.32)0.93751.14 (0.50–2.61)0.7534
rs2389557AA1948 (4.1)Ref Ref
GA36312 (3.3)0.79 (0.32–1.98)0.62180.78 (0.31–2.00)0.6094
GG1669 (5.4)1.33 (0.50–3.54)0.56391.24 (0.45–3.42)0.6722
rs2544381GG38820 (5.2)Ref Ref
CG2847 (2.5)0.46 (0.19–1.12)0.08620.48 (0.19–1.18)0.108
CC521 (1.9)0.36 (0.05–2.75)0.32490.37 (0.05–2.92)0.349
rs2544390CC1999 (4.5)Ref Ref
CT37014 (3.8)0.83 (0.35–1.95)0.66980.70 (0.29–1.71)0.4308
TT1545 (3.2)0.71 (0.23–2.16)0.54420.68 (0.22–2.15)0.5131
rs4667591TT24510 (4.1)Ref Ref
GT34711 (3.2)0.77 (0.32–1.84)0.55580.82 (0.33–2.03)0.671
GG1328 (6.1)1.52 (0.58–3.94)0.39291.68 (0.62–4.56)0.3081
rs7600336CC2368 (3.4)Ref Ref
TC34214 (4.1)1.22 (0.50–2.95)0.66431.10 (0.44–2.75)0.8428
TT1486 (4.1)1.20 (0.41–3.54)0.73571.30 (0.43–3.93)0.6374
VDR
rs10783219AA99826 (2.6)Ref Ref
TA127536 (2.8)1.09 (0.65–1.81)0.75121.03 (0.61–1.74)0.8982
TT4289 (2.1)0.80 (0.37–1.73)0.57530.76 (0.35–1.66)0.4869
rs11568820CC2196 (2.7)Ref Ref
TC35017 (4.9)1.81 (0.70–4.67)0.21821.93 (0.72–5.14)0.19
TT1535 (3.3)1.20 (0.36–4.00)0.76751.48 (0.43–5.14)0.5328
rs2228570GG21210 (4.7)Ref Ref
GA36415 (4.1)0.87 (0.38–1.97)0.73510.86 (0.37–2.00)0.7204
AA1483 (2.0)0.42 (0.11–1.55)0.1910.41 (0.11–1.56)0.1905
rs2238136CC48319 (3.9)Ref Ref
TC21810 (4.6)1.17 (0.54–2.57)0.68791.16 (0.51–2.60)0.7259
rs2853559GG31916 (5.0)Ref Ref
GA31310 (3.2)0.62 (0.28–1.40)0.25310.68 (0.29–1.56)0.3599
AA872 (2.3)0.45 (0.10–1.98)0.28750.40 (0.09–1.84)0.2418
rs4334089GG2277 (3.1)Ref Ref
AG35015 (4.3)1.41 (0.56–3.51)0.46371.45 (0.56–3.75)0.4397
AA1486 (4.1)1.33 (0.44–4.03)0.6171.54 (0.49–4.85)0.4563
rs7975232CC38211 (2.9)Ref Ref
CA28515 (5.3)1.87 (0.85–4.14)0.1211.84 (0.81–4.20)0.1458
AA603 (5.0)1.78 (0.48–6.56)0.38941.77 (0.46–6.78)0.4031
Abbreviations: VitD, vitamin D; HDP, hypertensive disorders in pregnancy. * Adjusted for pre-pregnancy BMI, maternal age, gestational weight gain, educational level, parity, and basal blood pressure.
Table 6. The interaction between SNPs and VitD in three trimesters on the risk of HDP.
Table 6. The interaction between SNPs and VitD in three trimesters on the risk of HDP.
SNPGenotypesVitDT1T2T3
NOR (95%CI)NOR (95%CI)NOR (95%CI)
CYP24A1
rs2209314CC/CT≥20545Ref724Ref557Ref
TT≥203161.60 (0.66–3.87)3831.41 (0.65–3.08)2941.86 (0.80–4.28)
CC/CT<2010261.18 (0.57–2.45)5411.13 (0.53–2.40)2601.20 (0.43–3.29)
TT<205601.52 (0.70–3.32)3091.41 (0.62–3.18)1451.04 (0.28–3.79)
Pinteraction 0.8774 0.7996 0.5739
rs2248137GG≥20357Ref451Ref339Ref
GC≥201172.50 (0.81–7.73)1114.45 (1.25–15.79) *1063.81 (1.11–13.07) *
CC≥202201.04 (0.32–3.31)2693.90 (1.33–11.41) *1952.22 (0.66–7.47)
GG<204800.77 (0.28–2.09)2832.11 (0.67–6.63)1311.17 (0.22–6.23)
GC<203171.80 (0.70–4.67)1335.42 (1.71–17.23) *631.05 (0.12–9.51)
CC<203771.65 (0.64–4.26)1711.17 (0.22–6.04)902.47 (0.56–10.81)
Pinteraction 0.5195 0.1442 0.8045
rs6127118GG≥20288Ref405Ref316Ref
AG/AA≥205722.01 (0.66–6.10)7020.83 (0.38–1.81)5372.97 (0.99–8.88)
GG<205822.29 (0.77–6.87)2910.77 (0.28–2.12)1331.86 (0.40–8.58)
AG/AA<2010101.63 (0.56–4.77)5581.12 (0.52–2.42)2712.18 (0.62–7.65)
Pinteraction 0.855 0.2332 0.5262
CYP27B1
rs10877012TT≥20363Ref445Ref331Ref
GT≥203721.14 (0.40–3.23)5101.24 (0.53–2.93)4051.12 (0.42–2.94)
GG≥201212.67 (0.86–8.29)1511.84 (0.60–5.61)1152.56 (0.83–7.85)
TT<206611.28 (0.52–3.15)3611.08 (0.41–2.84)1680.84 (0.21–3.27)
GT<207251.41 (0.59–3.41)3701.26 (0.51–3.12)1801.32 (0.41–4.21)
GG<202011.79 (0.58–5.49)1142.80 (0.97–8.10)542.26 (0.45–11.45)
Pinteraction 0.6139 0.4639 0.9865
CYP3A4
rs2242480CC≥20490Ref644Ref504Ref
CT≥203140.89 (0.34–2.31)4041.16 (0.51–2.64)2931.08 (0.44–2.68)
TT≥20541.74 (0.37–8.17)591.57 (0.35–7.15)522.47 (0.66–9.19)
CC<208991.04 (0.51–2.13)4801.53 (0.75–3.13)2471.25 (0.48–3.22)
CT<206041.16 (0.54–2.47)3180.53 (0.17–1.61)1420.51 (0.11–2.33)
TT<20871.52 (0.41–5.62)513.13 (0.85–11.47)173.30 (0.39–28.05)
Pinteraction 0.7911 0.696 0.5117
GC
rs1155563TT≥20331Ref440Ref326Ref
TC≥204010.81 (0.31–2.11)4981.04 (0.44–2.44)3920.45 (0.17–1.15)
CC≥201270.81 (0.21–3.09)1671.31 (0.44–3.95)1330.87 (0.27–2.80)
TT<205311.26 (0.55–2.88)2721.70 (0.70–4.15)1230.21 (0.03–1.65)
TC<207660.78 (0.34–1.80)4031.02 (0.40–2.55)1990.92 (0.33–2.52)
CC<202900.89 (0.33–2.42)1720.87 (0.26–2.89)830.94 (0.25–3.54)
Pinteraction 0.9456 0.5655 0.2934
rs16846876AA≥20437Ref557Ref419Ref
AT≥203661.22 (0.48–3.07)4501.18 (0.52–2.70)3530.34 (0.12–0.99)
TT≥20581.31 (0.27–6.38)1021.46 (0.43–4.97)811.83 (0.62–5.37)
AA<207141.65 (0.75–3.60)3791.71 (0.77–3.79)1680.51 (0.14–1.84)
AT<206720.87 (0.37–2.07)3621.00 (0.40–2.51)1821.00 (0.38–2.69)
TT<202061.12 (0.37–3.36)1080.72 (0.15–3.38)550.48 (0.06–3.85)
Pinteraction 0.4024 0.3482 0.8437
rs17467825AA≥20447Ref578Ref430Ref
GA≥203520.79 (0.31–1.99)4371.06 (0.48–2.34)3370.41 (0.15–1.07)
GG≥20620.49 (0.06–3.95)940.59 (0.12–3.00)880.56 (0.13–2.33)
AA<206921.15 (0.55–2.38)3471.55 (0.71–3.40)1520.36 (0.08–1.61)
GA<206910.91 (0.43–1.93)3830.84 (0.35–2.04)1910.85 (0.32–2.25)
GG<202130.65 (0.21–1.98)1210.81 (0.22–2.96)620.76 (0.17–3.48)
Pinteraction 0.8118 0.8253 0.2029
rs2282679TT≥20449Ref578Ref429Ref
GT≥203490.81 (0.32–2.04)4341.27 (0.57–2.80)3360.42 (0.16–1.10)
GG≥20630.47 (0.06–3.80)980.62 (0.12–3.11)890.56 (0.13–2.34)
TT<206881.10 (0.53–2.29)3451.67 (0.75–3.71)1530.36 (0.08–1.60)
GT<206890.97 (0.46–2.04)3810.92 (0.37–2.25)1880.87 (0.33–2.30)
GG<202170.65 (0.21–1.97)1230.88 (0.24–3.22)630.76 (0.17–3.48)
Pinteraction 0.7444 0.6979 0.2075
rs2298849AA≥20332Ref430Ref346Ref
GA≥204113.12 (1.02–9.59) *5110.92 (0.41–2.05)3860.93 (0.38–2.27)
GG≥201201.36 (0.24–7.60)1700.43 (0.09–1.95)1241.23 (0.37–4.10)
AA<206922.29 (0.77–6.85)3780.72 (0.29–1.82)1731.01 (0.33–3.11)
GA<207002.33 (0.78–6.95)3621.28 (0.57–2.87)1781.25 (0.41–3.82)
GG<202041.94 (0.51–7.40)1111.01 (0.28–3.70)54
Pinteraction 0.5918 0.1481 0.4554
rs2298850GG≥20441Ref564Ref416Ref
CG≥203540.77 (0.31–1.95)4371.34 (0.60–3.01)3410.41 (0.15–1.10)
CC≥20630.46 (0.06–3.71)990.66 (0.13–3.33)890.56 (0.13–2.37)
GG<206721.04 (0.50–2.20)3391.83 (0.81–4.12)1500.38 (0.09–1.70)
CG<206940.90 (0.42–1.91)3850.85 (0.33–2.20)1910.75 (0.26–2.12)
CC<202150.64 (0.21–1.96)1220.94 (0.25–3.50)630.79 (0.17–3.65)
Pinteraction 0.759 0.5603 0.2452
rs3755967CC≥20448Ref575Ref427Ref
CT≥203510.80 (0.32–2.01)4381.25 (0.57–2.76)3400.41 (0.16–1.08)
TT≥20630.47 (0.06–3.78)980.61 (0.12–3.10)890.56 (0.13–2.32)
CC<206851.09 (0.52–2.29)3451.66 (0.75–3.69)1530.36 (0.08–1.58)
CT<206950.96 (0.46–2.02)3830.91 (0.37–2.24)1890.86 (0.33–2.27)
TT<202170.64 (0.21–1.96)1230.87 (0.24–3.20)630.75 (0.16–3.45)
Pinteraction 0.7354 0.7072 0.2031
rs4588GG≥20448Ref571Ref422Ref
GT≥203480.81 (0.32–2.05)4391.23 (0.56–2.72)3390.40 (0.15–1.07)
TT≥20640.45 (0.06–3.57)980.60 (0.12–3.05)900.54 (0.13–2.24)
GG<206761.11 (0.53–2.32)3411.66 (0.75–3.70)1490.37 (0.08–1.65)
GT<206960.91 (0.43–1.93)3820.80 (0.31–2.05)1920.70 (0.25–1.98)
TT<202150.65 (0.21–1.99)1240.84 (0.23–3.11)630.75 (0.16–3.46)
Pinteraction 0.8164 0.6081 0.2578
rs7041AA≥20423Ref573Ref454Ref
CA≥203582.46 (0.96–6.32)4530.92 (0.42–2.04)3370.65 (0.27–1.58)
CC≥20820.97 (0.12–8.11)840.55 (0.07–4.25)630.50 (0.06–3.91)
AA<208922.03 (0.87–4.75)4781.15 (0.55–2.40)2260.80 (0.30–2.13)
CA<206121.21 (0.46–3.14)3161.04 (0.43–2.50)1580.61 (0.17–2.18)
CC<20882.18 (0.54–8.77)5621
Pinteraction 0.2275 0.8097 0.8926
LRP2
rs10210408CC≥20289Ref354Ref287Ref
TC≥204181.75 (0.60–5.10)5581.45 (0.58–3.62)4150.69 (0.29–1.62)
TT≥201551.56 (0.40–6.00)1991.25 (0.39–4.03)1540.30 (0.07–1.42)
CC<205261.83 (0.66–5.12)2941.10 (0.36–3.33)1270.41 (0.09–1.94)
TC<207821.64 (0.61–4.40)3991.96 (0.79–4.89)2140.86 (0.32–2.31)
TT<202901.16 (0.34–3.91)1580.64 (0.13–3.15)640.44 (0.06–3.55)
Pinteraction 0.497 0.9331 0.3864
rs2228171TT≥20291Ref374Ref292Ref
CT≥20920.97 (0.19–4.86)892.33 (0.65–8.38)630.98 (0.21–4.70)
CC≥20841.61 (0.39–6.59)1040.91 (0.19–4.49)740.37 (0.05–3.00)
TT<205481.21 (0.48–3.07)2941.06 (0.36–3.16)1400.60 (0.16–2.29)
CT<202861.68 (0.62–4.56)962.12 (0.59–7.56)60
CC<201941.26 (0.38–4.14)1112.37 (0.67–8.41)591.22 (0.25–5.97)
Pinteraction 0.9483 0.4054 0.5358
rs2389557AA/GA≥20113Ref102Ref87Ref
GG≥20372.15 (0.18–25.61)351.66 (0.22–12.36)262.11 (0.21–20.73)
AA/GA<204332.50 (0.55–11.48)1561.04 (0.20–5.46)1000.55 (0.07–4.57)
GG<201253.57 (0.70–18.10)337.09 (1.21–41.47)18
Pinteraction 0.7002 0.8493 0.2569
rs4667591TT/GT≥20120Ref116Ref94Ref
GG≥20317.98 (0.64–98.90)216.84 (0.76–61.89)190.90 (0.04–21.95)
TT/GT<204584.84 (0.63–37.13)1461.74 (0.35–8.58)930.15 (0.01–2.10)
GG<201007.64 (0.88–66.52)447.44 (1.11–49.79) *261.60 (0.13–19.78)
Pinteraction 0.8345 0.7612 0.3993
VDR
rs10783219AA≥20309Ref431Ref336Ref
TA≥204121.15 (0.44–3.05)5230.80 (0.36–1.77)3971.66 (0.67–4.12)
TT≥201400.96 (0.24–3.87)1560.41 (0.09–1.90)1200.77 (0.16–3.86)
AA<206021.28 (0.52–3.15)2970.78 (0.30–2.02)1451.44 (0.41–5.02)
TA<207391.15 (0.47–2.78)4191.00 (0.45–2.21)1940.71 (0.18–2.78)
TT<202530.90 (0.29–2.80)1351.04 (0.33–3.31)652.65 (0.66–10.67)
Pinteraction 0.8452 0.1733 0.945
rs2238136CC≥20111Ref99Ref83Ref
TC/TT≥20404.94 (0.42–58.60)386.78 (0.84–54.44)300.77 (0.06–9.56)
CC<203635.26 (0.68–40.67)1264.10 (0.66–25.32)770.44 (0.06–3.28)
TC/TT<201963.82 (0.45–32.40)652.57 (0.29–22.94)42
Pinteraction 0.4588 0.0586 0.9648
Abbreviations: VitD, vitamin D; HDP, hypertensive disorders in pregnancy. Adjusted for pre-pregnancy BMI, maternal age, educational level, parity, basal blood pressure. * p < 0.05.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Si, S.; Mo, M.; Cheng, H.; Peng, Z.; Alifu, X.; Zhou, H.; Chi, P.; Zhuang, Y.; Yu, Y. The Association of Vitamin D and Its Pathway Genes’ Polymorphisms with Hypertensive Disorders of Pregnancy: A Prospective Cohort Study. Nutrients 2022, 14, 2355. https://0-doi-org.brum.beds.ac.uk/10.3390/nu14112355

AMA Style

Si S, Mo M, Cheng H, Peng Z, Alifu X, Zhou H, Chi P, Zhuang Y, Yu Y. The Association of Vitamin D and Its Pathway Genes’ Polymorphisms with Hypertensive Disorders of Pregnancy: A Prospective Cohort Study. Nutrients. 2022; 14(11):2355. https://0-doi-org.brum.beds.ac.uk/10.3390/nu14112355

Chicago/Turabian Style

Si, Shuting, Minjia Mo, Haoyue Cheng, Zhicheng Peng, Xialidan Alifu, Haibo Zhou, Peihan Chi, Yan Zhuang, and Yunxian Yu. 2022. "The Association of Vitamin D and Its Pathway Genes’ Polymorphisms with Hypertensive Disorders of Pregnancy: A Prospective Cohort Study" Nutrients 14, no. 11: 2355. https://0-doi-org.brum.beds.ac.uk/10.3390/nu14112355

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