No association of vitamin D binding protein gene polymorphisms with ulcerative colitis in Chinese patients

Object Vitamin D (VD) deficiency has been reported in patients with ulcerative colitis (UC), and polymorphism in the gene encoding the vitamin D binding protein (DBP) can affect the characteristics of DBP, thus affecting the level and function of VD in vivo . Previous studies have rarely reported on the potential relationship between DBP polymorphisms and UC. To investigate the associations between genetic variants in DBP genes and UC susceptibility in the Han Chinese population, in order to discern whether any differences exist between this population and those of other countries. Methods In this case-control study, the genotyping of DBP rs4588 and rs7041 polymorphisms was conducted using polymerase chain reaction (PCR)-ligase detection reactions, and the rs4588 and rs7041 genotypes were detected by PCR-restriction fragment length polymorphism. Results In our case-control cohort, no significant difference was observed in the UC risk for either of the two SNPs (rs4588 and rs7401) in the DBP genes ( P > 0.05). No association between UC susceptibility and the DBP gene haplotypes was found either. Conclusions Our results suggest that the two SNPs (rs4588 and rs7401) in the DBP genes may have no correlation with susceptibility to UC in the Chinese Han population. But interestingly, haplotype GC, which contains the rs4588 and rs7041 variants in the DBP gene, may affect the level of oxidative stress in UC patients, especially the level of MDA.


Introduction
Ulcerative colitis (UC) is a chronic, non-specific inflammatory disease of the colon and rectum with unknown aetiology, the incidence of which has been increasing worldwide over time [1]. The lesions, which are mainly localized in the colonic mucosa and submucosa, appear most commonly in the sigmoid colon and rectum. The pathogenesis of UC is multifactorial, involving genetic predisposition, environmental factors, dietary changes, immune disorders, and intestinal mucosal barrier dysfunction [2,3]. The main clinical manifestations of UC are diarrhoea, mucus-like bloody stools, abdominal pain, a prolonged course with repeated attacks, and this affects each patient's physical and mental health and quality of life [4].
Vitamin D (VD), which comprises a group of fat-soluble molecules, is converted into 1,25(OH) 2 D 3 in the body by two hydroxylation reactions, after which the converted form plays a biological function [5]. In addition to its well-known roles in calcium absorption and phosphorus metabolism, VD is also associated with the occurrence and development of UC [6]. In fact, VD deficiency has been reported in 60% of patients with inflammatory bowel disease [7]. This vitamin also plays a role in intestinal health mainly by resisting infection and regulating the immune response, including maintaining the balance of intestinal flora and the barrier function of the intestinal epithelium [8].
In the process of VD metabolism, the vitamin D binding protein (DBP) binds and transports VD and its metabolites to target organs, while also regulating the level of VD in the body [9][10][11]. Polymorphisms in the DBP gene can affect the characteristics of the protein it encodes, thus affecting the level and function of VD in vivo. DBP is an α2-glycosylated globulin[12] whose encoding gene, which is highly polymorphic and contains 13 exons and 12 introns, is located at positions 12 to 13 of the long arm of chromosome 4. The DBP is 458 amino acids long, and the two most functionally important single nucleotide polymorphisms (SNPs) in it are rs7041 and rs4588. Rs7041 is a missense mutation of GAT-GAG, which changes aspartic acid at position 416 to glutamic acid [9], and rs4588 is a missense mutation of ACG-AAG, which changes threonine at position 420 to lysine. There is linkage disequilibrium between these two loci, and their haplotype forms the most important three polymorphisms in the DBP gene, namely Gc1F, Gc1S, and Gc2 (Table 1) [13,14]. It is also known that different genotypes of the DBP and VD affect their binding affinities, with the affinity of VD metabolites to Gc1F being the strongest, and vice versa for Gc2. The stronger the affinity, the more VD exists in its combined form, a consequence of which is that free VD is reduced, and active VD entering the target tissue is reduced accordingly.
However, to the best of our knowledge, previous studies on VD have rarely reported on the potential relationship between DBP polymorphisms and UC. We herein conducted a case-control study to identify associations between DBP gene variants and UC susceptibility in the Han Chinese population to determine whether any differences exist between this population and those from other countries.

Study population
In this case-control study, we recruited 70 patients with UC and 84 healthy controls from December 2016 to December 2018 in Shanxi Province, China (The calculation formula of sample size obtained that 70 samples were needed for each group). The specific exclusion and inclusion criteria are shown in (Table 2). Informed consent was obtained from all participants before the study commenced. The protocol, which was approved by the Ethical Committee of Shanxi Medical University, China, is in accordance with the Declaration of Helsinki statements. The study questionnaire (which included demographic data and any factors affecting VD levels in vivo, such as milk consumption and outdoor sunlight exposure) was collected from the recruited volunteers, and the 5 ml of fasting venous blood collected from each volunteer into a heparinized tube was centrifuged. The plasma, red blood cells and white blood cells from each sample were collected and stored at −80°C for later use.
Genomic DNA extraction and genotype analysis Peripheral blood DNA was extracted strictly according to the instructions from the QIAamp DNA mini kit (QIAGEN, Germany), and the recovered samples were stored at −20°C for preservation. DBP rs4588 and rs7041 characterisation was conducted using polymerase chain reaction-ligase detection reactions (PCR-LDRs) with TaqMan genotyping assays on the 3730XL DNA analyser (ABI Co., USA). The sequences of the forward and reverse primers were 5'-GTT TTT CAG ACT GGC AGA GCG-3' and 5'-ACA CCA GGA AAA GCC TGT CAC-3', respectively, and the annealing temperature of the reactions was 60°C. The length of the amplified product was 259 bp.

Statistical analysis
All statistical analyses were performed using the SPSS 22.0 statistical package. The measurement data were expressed as mean±SD, and the counting data were expressed as frequencies (%). Analysis of variance (ANOVA) and the chi-square test were used for continuous variables and categorical variables, respectively. For continuous variables that did not satisfy the test for normality or homogeneity of variance, the Mann-Whitney U test in the nonparametric test was adopted and the results were expressed as median values and upper and lower quartiles [median(Q 25 -Q 75 )]. Hardy-Weinberg equilibrium was tested using the chi-square test and SHEsis software to analyse linkage imbalance and haplotypes. When the haplotype distribution frequency was less than 0.03, the software would automatically ignore the comparison between the two groups. Logistic regression analysis was used for the multivariate analysis, the results of which were considered statistically significant at P < 0.05 (two-tailed).

Study population and VD influencing factors
Our comparison of the baseline characteristics and VD influencing factors between the case and control groups is shown in Table 3 there was no significant differences in the values shown in Table 3 (P > 0.05).

Serum VD, DBP, inflammatory factors and oxidation indexes
The analysis results for serum VD, DBP, inflammatory factors and antioxidant indicators in the case and control groups are shown in Table 4. Among these results, the levels of the case group's serum DBP, TNF-α and IL-1β were higher than those of the control group (P < 0.05), whereas the serum 25(OH)D 3, 1,25(OH) 2 D 3, IL-23, MPO and SOD levels of the control group were higher than those of the case group (P < 0.05).
DBP gene sequencing results and Hardy-Weinberg equilibrium testing Sequence analysis of the PCR product from the rs4588 variant in the DBP gene revealed three sequence types (AA, CC and AC), with AA and CC being single peaks, and AC the overlapping peak (Figure 1 A, B and C). The gene sequence for rs7041 was found to be GG, TT, and GT (Figure 1 D, E and F).
The distributions results for the rs4588 and rs7041 genotypes are shown in Table 5. The genotypes for these two gene loci in the case and the control groups were evaluated using a Hardy-Weinberg equilibrium test, the results of which generated a value of P < 0.05, indicating that the selected population met the study requirements and was representative of the population as a whole (Table 5).

DBP polymorphisms and UC in the study groups
The genotype distributions and allele frequencies of the DBP gene variants rs4588 and rs7041 in the control and case groups are shown in Table 6, the differences of which were not statistically significant (P > 0.05).

Linkage disequilibrium and haplotype analyses
The linkage imbalance analysis showed that D' = 1.000 and r 2 = 0.162 ( Figure 2) between rs4588 and rs7041. We speculated that there is a chain disequilibrium reaction between rs4588 and rs7041. Haplotype analysis showed that there was no significant difference in frequency among the three main haplotypes between the case group and the control group (P > 0.05, Table 7).

Haplotype analysis results
The results of the analysis on VD, DBP, inflammatory factors and antioxidant indexes for the different haplotypes in the case and control groups are shown in Table 8. In the case group, the different haplotypes differed in their MDA distributions (P = 0.014), but there was no significant difference in the other indicators (P > 0.05).

Logistic regression analysis
Grouping was used as the dependent variable in the logistic stepwise regression analysis, with the cutoff value of α in = 0.05, α out = 0.10, a model likelihood ratio of 2 = 270.097, and statistical significance set at P < 0.001. Age, DBP, VD levels, haplotypes, serum inflammation markers and oxidative antioxidant levels were examined. Finally, a model with 6 independent variables showed statistical significance, and gender and 1,25(OH) 2 D 3 were found to be protective factors. The results from the model are shown in Table 9. Discussion UC, one of the main manifestations of inflammatory bowel disease (IBD), most commonly presents with abdominal pain, diarrhoea, and purulent mucus discharge with blood in the stools from patients. A few patients with this condition present with severe clinical disease of rapid inset, the progression of which develops quickly.
[15]The clinical signs can include the frequent passage of bloody stools, up to 30 times a day, as well as high fever, anaemia, nutritional disorders and weight loss. In addition to these signs, the present study also found that the BMI of patients in the case group was lower than that of the control group. In healthy people, the higher the BMI, the lower the VD level, and the lower the BMI, the higher the VD level, because of the fat dilution effect [16].
However, the morbid weight loss seen in patients with UC will reduce the VD level. One study found that more than half of its UC patients were low in VD [17]. Our investigations found that the levels of 25(OH)D 3 and 1,25(OH) 2 D 3 in patients with UC were lower than those in healthy people, but the opposite was found for DBP levels. The relationship between DBP and vitamin D concentration is unclear [18]. We hypothesized that the level of DBP needed for VD to function in the body is small when VD levels are low (as in the case group); consequently, the amount of free DBP is large when VD levels are low. When normal VD levels are present, as in the control group, more DBP is needed and there is less free DBP in the serum. Hence, people with UC will have higher DBP levels than healthy individuals, but lower VD levels.
Modern studies have argued that the pathogenesis of UC is multifactorial, and that imbalanced cytokine regulation is an important cause of this condition [19,20]. TNF-α is produced by a variety of immune cells, epidermal cells, endothelial cells and fibroblasts, and is involved in inflammatory and immune responses [21]. IL-1β, which is produced by activated macrophages as proproteins, is an important mediator of inflammatory responses and is involved in a variety of cellular activities, including cell proliferation, differentiation, and apoptosis. It has been shown that IL-1β is elevated in patients with UC and Crohn's disease and is also associated with disease severity [22]. It is well known that VD is widely involved in a variety of biological activities in the body, especially chronic diseases. Polymorphisms in VD pathway-related genes have been extensively studied in cancer, chronic inflammation, autoimmune diseases, and lipid metabolism disorders[9, 32, 33]. In recent years, some studies have highlighted the importance of VD in the pathogenesis of IBD [34,35]. DBP is an important component in the vitamin D signalling pathway and has multiple functions in the body, including transport of VD, regulation of inflammatory reactions in vivo and immune regulation [36].It plays a key role in signal transduction, and its GC coding gene polymorphism exists in more than 2000 SNP loci, and polymorphisms in these loci may affect the structure and concentration of the protein, leading to disease onset [10,37].
In the present study, we investigated the association between DBP variants and UC disease using a case-control design with 84 controls and 70 cases from the Chinese Han population. We found no significant association between the polymorphisms and risk of UC for genotype or haplotype. Our findings are consistent with a previous report by Luo et al [38]. Additionally, our logistic regression analysis further confirmed that the DBP gene haplotypes were independently correlated with UC after adjustment for other confounders.
A variety of factors affect the level and activity of VD during its activation, including the DBP, and most VD in the circulation is combined with the DBP. The biological activity of VD may depend mainly on this free component [12] . Therefore, the gene encoding DBP is a key candidate for the VD pathway and could play an important role in maintaining the overall level of VD and regulating the amount of free VD available. Mutations in the DBP gene cause changes in serum DBP levels, which are associated with changes in plasma VD concentrations [39][40][41]. A recent study reported that two SNPs in the DBP gene (rs7041 and rs4588) can produce three common variants affecting the amount of free VD metabolites. In addition, the same study also determined that SNPs in DBP (rs7041, rs4588 and rs2282679) were the genetic determinants of reduced 25(OH)D levels [42].
More than 120 DBP gene variants have been found through gene sequencing, and the existence of these variants affect the function of DBP, thus affecting the affinity between 25(OH)D 3 and DBP, making the ability of the DBP to bind VD differ significantly in the variants [43]. Previous studies have shown that DBP gene variants are associated with the risk of autoimmune diseases, such as thyroid autoimmune diseases, diabetes and asthma [44][45][46]. In a European case-control study, the DBP gene with a homozygous rs4588 SNP was significantly associated with the risk of Crohn's disease. Similarly, haplotype Gc2, which consists of rs4588 and rs7401 variants, appears to be a protective factor in the UC group [47]. There are few reports on DBP gene polymorphism and UC susceptibility in the Asian population. The association between genetic variation in the DBP and IBD has been reported mainly in a European cohort [47]. Our results suggest that there is no association between these variants and the risk of UC. As there is little evidence of a link between DBP variation and IBD susceptibility, multicentre studies will be needed to confirm these results.
In the present study, we found an interesting result. Specifically, in the patients with UC, the haplotypes are significantly associated with the MDA level, and the highest level of MDA is found when the haplotype is CG. Studies have shown that UC is often accompanied by a certain degree of oxidative stress injury, which stimulates the excessive generation of oxygen free radicals, leading to a decline in the activity of SOD, and inducing an increase in MDA concentration and aggravating the disease [48]. variants, may be a risk factor for the UC group.

Conclusion
In summary, our study is the first to reveal the relationship between the DBP gene and UC in the Han Chinese population. The results show that haplotype GC, which comprises rs4588 and rs7041 variants of the DBP gene, may affect the level of oxidative stress in UC patients, especially their MDA levels. However, the results need to be replicated and confirmed in a multi-centre study.

Study limitations
Our study has some limitations. The number of participants included was small, especially the patients. The source range of the population was also narrow, so its representativeness is not ideal. In addition, the limited economic conditions meant that we only examined two loci in the DBP gene, and the confounding factors were not

Availability of data and materials
All data generated or analysed during this study are included in this published article.

Consent for publication
Agree to published.  Tables   Table 1 Relationship