Figures
Abstract
Background
Polymorphisms of genes encoding components of the vitamin D pathway including vitamin D receptor (VDR) and vitamin D binding protein (DBP) have been widely investigated because of the complex role played by vitamin D in cancer tumorogenesis. In this study, we investigated the association between VDR and DBP gene polymorphisms and HBV-related HCC risk in a Chinese population.
Methods
Study subjects were divided into three groups: 184 HBV patients with HCC, 296 HBV patients without HCC, and 180 healthy controls. The VDR rs2228570, and rs3782905 and the DBP rs7041 polymorphisms were genotyped using PCR-RFLP and the VDR rs11568820 polymorphism was genotyped by PCR-SSP, respectively. DNA sequencing was performed to validate the genotype results.
Results
We found that there were significant differences in the genotype and allele frequencies of the VDR rs2228570 and DBP rs7041 polymorphisms between HBV patients with HCC and healthy controls. The rs2228570 T allele was associated with a significant increased HBV-related HCC risk as compared with the C allele. The rs2228570 TT and TT/TC genotypes were correlated with a significant increased HBV-related HCC risk when compared with the wild-type CC homozygote. Similarly, the rs7041 G allele was associated with a significant increased HBV-related HCC risk as compared with the T allele. The rs7041 GG and GG/TG genotypes were correlated with a significant increased HBV-related HCC risk when compared with the wild-type TT homozygote. However, we did not observe any significant effect of VDR rs11568820, and rs3782905 polymorphisms on HBV-related HCC risk in this population. In haplotype analysis, we also did not find any significant differences in haplotype frequencies of the VDR gene between HBV patients with HCC and the healthy controls.
Conclusions
We conclude that the VDR rs2228570 and DBP rs7041 polymorphisms may contribute to increased susceptibility to HBV-related HCC in the Chinese population. Due to the marginal significance, further large and well-designed studies in diverse ethnic populations are needed to confirm our results.
Citation: Peng Q, Yang S, Lao X, Li R, Chen Z, Wang J, et al. (2014) Association of Single Nucleotide Polymorphisms in VDR and DBP Genes with HBV-Related Hepatocellular Carcinoma Risk in a Chinese Population. PLoS ONE 9(12): e116026. https://doi.org/10.1371/journal.pone.0116026
Editor: Balraj Mittal, Sanjay Gandhi Medical Institute, India
Received: August 17, 2014; Accepted: December 1, 2014; Published: December 26, 2014
Copyright: © 2014 Peng et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the paper.
Funding: The study was supported by the National Natural Science Foundation of China (No. 81260302) and the Guangxi education and innovation projects for graduate student in 2014 (No. YCBZ2014031). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Hepatocellular carcinoma (HCC) is the fifth most common malignancy and the third leading cause of cancer mortality worldwide [1]. The distribution of HCC is imbalance throughout the world, with particularly high incidence rates in Eastern Asia, and Sub-Saharan Africa, and China alone accounted for more than 50% of all HCC cases [1]. Major risk factors for HCC include infection with hepatitis B virus (HBV) and hepatitis C virus (HCV), dietary aflatoxin B1 exposure, and alcohol consumption [2], [3]. Overall, 75–85% of HCC patients are attributable to persistent viral infection of HBV, especially in developing countries such as China [1]. However, only a minority of lifelong chronic carriers of HBV eventually develops HCC, suggesting the importance of genetic susceptibility in HBV-related HCC.
Vitamin D is a fat soluble secosteroid which is involved in a wide variety of biological processes like bone metabolism, modulation of immune response, cell proliferation and cell differentiation [4]. Studies have reported that vitamin D implicated in inhibition of carcinogenesis by induction of differentiation, promotion of apoptosis, and inhibition of proliferation and angiogenesis [5]. Vitamin D receptor (VDR) is a crucial mediator for the cellular effects of vitamin D and additionally interacts with other cell-signaling pathways that influence cancer development [6]. It is an intracellular hormone receptor that specifically binds the biologically active form of vitamin D and interacts with specific nucleotide sequences (response elements) of target genes to produce a variety of biological effects. The VDR gene is located on chromosome 12q12–q14 and is highly polymorphic [7]. Many single-nucleotide polymorphisms (SNPs) have been identified and genetic variations in the VDR gene may phenotypically appear as interindividual variations in limiting rates of vitamin D synthesis in the skin, hydroxylation in the liver and kidney, transport, metabolism and degradation that would ultimately influence the anti-tumor effect of vitamin D. Moreover, the VDR gene variants have been reported to associate with increased risk of cancers including breast [8], prostate [9], and colorectal [10].
Vitamin D binding protein (DBP) is the primary transport protein of 25(OH)D in circulation, with approximately 88% of 25(OH)D being bound to DBP [11]. In addition to being a carrier for vitamin D metabolites, DBP also has anti-inflammatory and immunoregulatory functions, and has been identified playing a role in several chronic disease including cancers [12]. The gene for the DBP is also a logical candidate in the vitamin D pathway because of its major function in transporting vitamin D metabolites in the circulation. Variants in this gene have been shown to alter plasma concentrations of 25-OHD [13]. In addition, the variants in the DBP gene have been implicated in cancers such as breast [14] and prostate [15].
In light of the important biological function of VDR and DBP genetic polymorphisms in cancer development and progression, we hypothesized that genetic variants of the VDR and DBP genes were associated with increased susceptibility to HBV-related HCC. VDR rs11568820, rs2228570, and rs3782905 and DBP rs7041 polymorphisms are the most common studied SNPs which have been found modulate the expressions of vitamin D and play critical roles in carcinogenesis [12], [16], [17], [18]. Therefore, we performed a case–control study to investigate the association between VDR rs11568820, rs2228570, and rs3782905 and DBP rs7041 polymorphisms and HBV-related HCC susceptibility in a Chinese population.
Materials and Methods
Study population
This is a hospital-based case-control study performed in a total of 660 subjects, including 184 HBV patients with HCC, 296 HBV patients without HCC, and 180 healthy controls. All study subjects were consecutively recruited between July 2013 and November 2013 from the First Affiliated Hospital of Guangxi Medical University, Guangxi, China without any restrictions on age, gender, drinking and smoking. The HBV patients with HCC were diagnosed based on pathological findings combined with at least one positive HCC image on computed tomography (CT) or magnetic resonance imaging (MRI), sometimes combined with serum AFP analysis (>400 ng/ml). All patients with HBV infection selected for this study (including HCC and non-HCC patients) were further confirmed by being HBsAg (hepatitis B surface antigen) positive, HbcAb (hepatitis B virus core antibody) positive and HBeAg (hepatitis B e antigen) or HBeAb (hepatitis B e antibody) positive for at least 6 months. The presence of HCC was excluded from the non-HCC patients and healthy controls by magnetic resonance imaging (MRI), computed tomography (CT), ultrasonography, histology, and other laboratory tests. Moreover, patients with positive laboratory tests for hepatitis C virus (HCV; anti-HCV and/or HCV-RNA), human immunodeficiency virus (HIV), alcoholic liver disease, or suspected autoimmune diseases with antinuclear antibody titer greater than 1∶160 were excluded from this study. A sample of approximately 2–4 ml venous blood was collected from each participant. The following laboratory parameters were obtained for each participant at the time of whole-blood collection: serum total bilirubin (T-Bil), direct bilirubin (D-Bil), aspartate aminotransferase (AST), alanine aminotransferase (ALT), albumin (ALB), gamma-glutamyltransferase (GGT), alpha fetoprotein (AFP), and HBV-DNA. In addition, HBsAg, HBsAb, HbcAb, HBeAg, and HBeAb were also collected. All participants were exclusively enrolled from Guangxi district with written informed consent and the study was performed with the approval of the Ethics Committee of the First Affiliated Hospital of Guangxi Medical University.
Genomic DNA preparation and genotyping
Genomic DNA was extracted from 2 mL of peripheral blood using the standard phenol–chloroform method, dissolved in 300 µL of Tris–HCl buffer (10 mmol/L, pH 8.0) containing 1 mmol/L ethylenediaminetetraacetic acid (EDTA), and stored at −80°C until use. The SNPs of VDR rs2228570, and rs3782905 and DBP rs7041 were detected using polymerase chain reaction–restriction fragment length polymorphism (PCR-RFLP), and the VDR rs11568820 polymorphism was genotyped by polymerase chain reaction-sequence specific primer (PCR-SSP). Primer sequences, annealing temperatures, restriction enzymes, and product size in genotype analysis are shown in Table 1. To control the quality of genotyping, 50% of samples were repeated.
The gene fragments containing the polymorphic sites were amplified by PCR. PCR was performed in a total volume of 25 µL reaction mixture containing 1 U DNA Taq polymerase (Shanghai Biocolor, Shanghai, China), 1× buffer, 50 ng genomic DNA, 1.5 mmol/L MgCl2, 0.2 µmol/L primers, and 0.08 mmol/L dNTPs. The PCR cycles were as follows: 94°C for 5 min, 30 cycles of denaturing at 94 °C for 45 s, annealing at the indicated temperature for 45 s, extension at 72 °C for 45 s, and a final extension at 72 °C for 5 min in a Perkin-Elmer thermocycler (2700, Applied Biosystems, Foster City, CA, USA). After amplification, the PCR products were identified under ultraviolet light after electrophoresis for 15 min with 120 V in 2.5% agarose gel stained with ethidium bromide (EB). Then the PCR products were digested by allele-specific restriction enzymes overnight at 37°C. Cleaved DNA fragments were identified by ultraviolet light after electrophoresis in 2.5% agarose gel stained by EB (Figure 1–3).
Lane 10 shows marker; Lanes 1, 2 and 3 show rs2228570 TC genotype; lanes 4, 5 and 6 show rs2228570 TT genotype; lanes 7, 8 and 9 show rs2228570 CC genotypes.
Lane 10 shows marker; Lanes 1, 2 and 3 show rs3782905 GG genotype; lanes 4, 5 and 6 show rs3782905 CG genotype; lanes 7, 8 and 9 show rs3782905 CC genotypes.
Lane 10 shows marker; Lanes 1, 2 and 3 show rs7041 GG genotype; lanes 4, 5 and 6 show rs7041 TG genotype; lanes 7, 8 and 9 show rs7041 TT genotype.
VDR rs11568820 polymorphism was detected using forward primers designed with the last nucleotide complementary to the allelic variant substitution base as shown in Table 1. PCR products were illuminated by ultraviolet light after electrophoresis in 2.5% agarose gel stained by EB (Figure 4). The allelic type was determined according to the presence or absence of the desired length of PCR products (235 bp). In order to avoid pseudopositive or pseudonegative results in SSP-PCR, positive and negative controls were included each time. To confirm the genotyping results, 10% of PCR-amplified DNA samples in each group (totally 66 specimens) were examined by DNA sequencing, and the results were 100% concordant (Figure 5–8.).
Lane 10 shows marker; Lanes 1, 3 and 4 show rs11568820 AG genotype; Lanes 2, 6 and 7 show rs11568820 GG genotype; Lanes 5, 8 and 9 show rs11568820 AA genotype.
The arrows in –
show GG, AG, and AA genotypes, respectively
The arrows in –
show CC, TC, and TT genotypes, respectively
The arrows in –
show CC, GC, and GG genotypes, respectively
The arrows in –
show TT, TG, and GG genotypes, respectively
Statistical analysis
Continuous variables were expressed as mean±standard deviation (SD). Differences in demographic characteristics and laboratory parameters such as T-Bil, D-Bil, AST, ALT, ALB, GGT, AFP, and HBV-DNA were compared among groups with one-way ANOVA (for continuous variables) and χ2 test (for categorical variables). The Hardy–Weinberg equilibrium (HWE) was evaluated with a goodness of fit χ2-test with one degree of freedom to compare the observed genotype frequencies among the subjects with the expected genotype frequencies. Genotype and allele frequencies of VDR rs11568820, rs2228570, and rs3782905 and DBP rs7041 polymorphisms were compared among various groups using the χ2 test and fisher's exact test when appropriate, and the odds ratios (ORs) and 95% confidence intervals (CIs) were calculated by binary logistic regression adjusted for age, gender, smoking, and drinking to assess the relative risk conferred by a particular allele and genotype. To assess the effect of the genotypes carrying the variant allele we ran analysis under the dominant model. The linkage disequilibrium (LD) among the three polymorphisms in the VDR gene was quantified using Shi's standardized coefficient D′ [19]. The haplotypes and their frequencies were estimated based on a Bayesian algorithm using the Phase program [20]. All tests were two-sided with a significance level of P<0.05. All statistical analyses were performed using SPSS statistical software package version 13.0 (SPSS, Inc., Chicago, IL, USA).
Results
Demographic and clinical characteristics
The demographic and laboratory parameters of all subjects enrolled in this study were summarized in Table 2. The demographic data included age, gender, smoking, and drinking, and the laboratory parameters included T-Bil, D-Bil, AST, ALT, ALB, GGT, AFP, and HBV-DNA. Briefly, the healthy controls and the HBV patients without HCC were significantly younger than the HBV patients with HCC (P<0.001). The HBV patients with HCC had a higher proportion of smokers and drinkers compared with the other two groups (P<0.05). In addition, the HBV patients with HCC and the other two groups without HCC (healthy controls and HBV patients without HCC) had statistical significant different laboratory results for T-Bil, D-Bil, AST, ALT, ALB, GGT, AFP (P<0.001), and HBV-DNA (P<0.05).
Genotype and allele distributions of the VDR and DBP polymorphisms
The genotype and allele distributions of VDR and DBP polymorphisms in 184 HBV patients with HCC, 296 HBV patients without HCC, and 180 healthy controls were determined using PCR-RFLP and PCR-SSP methods. The genotype and allele distributions of VDR rs11568820, rs2228570, and rs3782905 and DBP rs7041 polymorphisms were shown in Table 3. The genotype frequencies for the 4 SNPs in healthy controls were consistent with the study of Asian populations (Han Chinese in Beijing) in the International HapMap project (http://hapmap.ncbi.nlm.nih.gov/cgi-perl/gbrowse/hapmap24_B36/). Furthermore, the genotype frequencies for the 4 SNPs were all consistent with HWE in the healthy controls (P = 0.238 for rs11568820, P = 0.351 for rs2228570, P = 0.650 for rs3782905, and P = 0.819 for rs7041, respectively).
HBV patients without HCC compare with healthy controls
Genotype and allele frequencies of the VDR rs11568820, rs2228570, and rs3782905 and DBP rs7041 polymorphisms between HBV patients without HCC and healthy controls are presented in Table 3. There were no significant differences in the genotype and allele frequencies of VDR rs11568820, rs2228570, and rs3782905 and DBP rs7041 polymorphisms between the HBV patients without HCC and healthy controls. Binary logistic regression analyses adjusted for age, gender, smoking and drinking also did not reveal any significant difference between the VDR rs11568820, rs2228570, and rs3782905 and DBP rs7041 polymorphisms and CHB risk.
HBV patients with HCC compare with HBV patients without HCC
Genotype and allele frequencies of the VDR rs11568820, rs2228570, and rs3782905 and DBP rs7041 polymorphisms between HBV patients with HCC and HBV patients without HCC are shown in Table 3. We also did not observe any significant differences in genotype and allele frequencies of the VDR rs11568820, rs2228570, and rs3782905 and DBP rs7041 polymorphisms between HBV patients with HCC and HBV patients without HCC. Similarly, binary logistic regression analyses adjusted for age, gender, smoking and drinking also did not reveal any significant difference between the VDR rs11568820, rs2228570, and rs3782905 and DBP rs7041 polymorphisms and HCC risk.
HBV patients with HCC compare with healthy controls
Genotype and allele frequencies of the VDR rs11568820, rs2228570, and rs3782905 and DBP rs7041 polymorphisms between HBV patients with HCC and healthy controls are summarized in Table 3. The results showed that there were no significant differences in the genotype and allele frequencies of the VDR rs11568820 and rs3782905 polymorphisms between the HBV patients with HCC and healthy controls. Binary logistic regression analyses adjusted for age, gender, smoking and drinking also did not reveal any significant difference between VDR rs11568820 and rs3782905 polymorphisms and HCC risk. However, with respect to VDR rs2228570 polymorphism, the frequencies of the CC, TC and TT genotypes were 29.4, 46.7 and 23.9% in healthy controls and were 21.7, 48.9 and 29.3% in HBV patients with HCC, respectively. The frequencies of the C and T alleles were 52.8 and 47.2% in healthy controls and were 46.2 and 53.8% in HBV patients with HCC, respectively. The results indicated that there were significant differences in the genotype and allele frequencies of the VDR rs2228570 polymorphism between HBV patients with HCC and healthy controls. Binary logistic regression analyses adjusted for age, gender, smoking and drinking showed that the rs2228570 T allele was associated with a significant increased HCC risk as compared with the C allele (adjusted OR = 1.49, 95%CI, 1.08–2.06, P = 0.022). The rs2228570 TT and the rs2228570 TT/TC genotypes were correlated with a significant increased HCC risk when compared with the wild-type CC homozygote (adjusted OR = 2.15, 95% CI, 1.13–4.08, P = 0.032 for TT, and adjusted OR = 1.82, 95% CI, 1.04–3.27, P = 0.039 for TT/TC). With respect to DBP rs7041 polymorphism, the frequencies of the TT, TG and GG genotypes were 52.2, 40.6 and 7.2% in healthy controls and were 47.3, 42.4 and 10.3% in HBV patients with HCC, respectively. The frequencies of the T and G alleles were 72.5 and 27.5% in healthy controls and were 68.5 and 31.5% in HBV patients with HCC, respectively. We observed that there were also significant differences in the genotype and allele frequencies of the DBP rs7041 polymorphism between HBV patients with HCC and healthy controls. Binary logistic regression analyses adjusted for age, gender, smoking and drinking revealed that the rs7041 G allele was associated with a significant increased HCC risk as compared with the T allele (adjusted OR = 1.50, 95% CI, 1.06–2.14, P = 0.034). The rs7041 GG and the rs7041 TG/GG genotypes were correlated with a significant increased HCC risk when compared with the wild-type CC homozygote (adjusted OR = 2.58, 95% CI, 1.12–5.92, P = 0.023 for TT, and adjusted OR = 1.74, 95% CI, 1.02–3.84, P = 0.039 for TG/GG).
Haplotype distributions of VDR gene in HBV patients with HCC and healthy controls
Haplotype analyses were performed in the VDR gene between HBV patients with HCC and healthy controls using the SHEsis software and the possible eight haplotype frequencies are shown in Table 4. Overall, strong linkage disequilibrium was observed between the rs2228570 C allele and the rs11568820 G allele (D' = 0.943) and rs3782905 C allele (D' = 0.962). Major GCC haplotype accounted for 25.7% and 22.5% of the eight haplotypes in both the HBV patients with HCC and healthy controls, respectively. However, the haplotype frequencies of the VDR gene in HBV patients with HCC were not significantly different from that in healthy controls (P>0.05).
Discussion
Hepatocellular carcinoma (HCC) is one of the most common causes of cancer mortality worldwide [1]. A better understanding of the mechanism involved in HCC is essential for us to develop new preventions and improve the existing therapies. It has been well accepted that HCC is a complex disease for which the underlying causes remain unclear. In this study, we investigated the VDR and DBP gene polymorphisms and determined whether the genetic factors are related to the occurrence of HBV-related HCC in a Chinese population. The results revealed that the VDR rs2228570 and DBP rs7041 polymorphisms were both significantly associated with HBV-related HCC. The rs2228570 T allele was associated with a significant increased HCC risk as compared with the C allele. The rs2228570 TT and the rs2228570 TT/TC genotypes were correlated with a significant increased HCC risk when compared with the wild-type CC homozygote. Similarly, the rs7041 G allele was associated with a significant increased HCC risk as compared with the T allele. The rs7041 GG and the rs7041 TG/GG genotypes were correlated with a significant increased HCC risk when compared with the wild-type CC homozygote. However, in haplotype analyses, we did not observe any significant association between VDR gene haplotypes and HBV-related HCC risk. Our results suggest that the VDR rs2228570 and DBP rs7041 polymorphisms may contribute to increased susceptibility to HBV-related HCC in the Chinese population. However, because the significant differences were only found in HBV patients with HCC versus healthy controls comparison but not in HBV patients with HCC versus HBV patients without HCC comparison, the results should be interpreted with caution.
Vitamin D is a steroid hormone with potent antitumor effects. Several investigations have demonstrated that 1, 25-dihydroxy vitamin D3 [1,25(OH)2D3], the active form of vitamin D, plays an essential role in the development of cancers by regulating the expression of tumor-related genes or mediating inhibition of growth of cancer cells [21], [22]. The VDR, a nuclear transcription factor that belongs to the steroid hormone receptor family, is an important regulator of the vitamin D pathway. It was reported that the VDR mediates gene expression of vitamin D [23] and most of the effects of active vitamin D in a wide variety of tissues [24]. In addition, in vitro studies have demonstrated that the VDR plays a critical role in regulating cell proliferation, differentiation and induction of apoptosis [25], [26]. A recent study by Beilfuss et al. [27] demonstrated that the polymorphisms in the VDR gene influenced VDR expressions and production and involved in the fibrotic processes in human hepatic stellate cells. The rs2228570 polymorphism is located in the 5′ end of the VDR gene. This polymorphism results in an alternative transcription initiation site, leading to a protein variant with three additional amino acids at the amino terminus [28]. The change in the VDR protein structure may lead to altered biological functions of vitamin D and involved in carcinogenesis [29]. Shafia et al. [30] reported that the VDR rs2228570 polymorphism was involved in the increased susceptibility to development and progression in multiple myeloma in an ethnic Kashmiri population, the homozygous (TT) genotype was associated with higher risk for developing multiple myeloma. Sarkissyan et al. [31] found that genetic variation of VDR rs2228570 influence CRC risk, particularly in African American cohorts. Similarly, Mohapatra et al. [4] reported that the VDR rs2228570 polymorphism was associated with an increased epithelial ovarian cancer risk in Indian population. In the present study, we found that the VDR rs2228570 polymorphism contributed to increased susceptibility to HBV-related HCC in Chinese population, which was consistent with the findings above.
DBP is a multifunctional protein known for its role in the transport of vitamin D metabolites. Almost all the 25(OH) D3 and active hormone 1, 25(OH)2D3 exist in the circulation bound to serum DBP [11]. It was reported that DBP regulate the half-life of 25(OH) D3 in the circulation through stabilization of the hormone and helps maintain serum vitamin D levels through reuptake by proximal tubule cells in the kidney [32]. In addition, DBP also has anti-inflammatory and immunoregulatory functions, and has been identified play important roles in several chronic disease including cancers [33]. The gene for the DBP is highly polymorphic, with more than 120 known variants and with different frequency distributions in diverse populations [34]. The rs7041 is one of the most common SNPs in the DBP gene which has been shown conversely associated with plasma concentrations of 25(OH)D3 [13], [35]. Zhou et al. [36] reported that the DBP rs7041 polymorphism was correlated with increased gastrointestinal cancer susceptibility in a Chinese population. Anderson et al. [14] demonstrated that the DBP rs7041 polymorphism interacts with Vitamin D exposure and contributes to increased breast cancer risk among a Caucasian Women in Canada. However, Poynter et al. [37] did not found a significant association between DBP rs7041 polymorphism and colorectal cancer in an American population. In this study, we observed that the DBP rs7041 polymorphism contributed to an increased HBV-related HCC risk in Chinese population. The inconsistent results may be attributed to the diverse physiological and pathophysiological effects of DBP rs7041 polymorphism on different ethnicities and different tissues.
Haplotypes are a group of closely relevant genetic markers present on one chromosome which tend to be inherited together and more frequent than expected by chance in a block pattern due to the presence of linkage disequilibrium [38]. In the present study, strong linkage disequilibrium was observed between the rs2228570 C allele and the rs11568820 G allele (D' = 0.943) and rs3782905 C allele (D' = 0.962) in the VDR gene, which was consistent with the results from Kaabachi et al. [39], suggesting that the three SNPs in the VDR gene may interact with each other and influence cancer risk. However, in the haplotype analysis, we did not observe any significant association between haplotypes and HBV-related HCC risk. The reason may be that the relatively small sample size in each haplotype has insufficient statistical power to detect a slight effect or may have generated a fluctuated risk estimate. Therefore, further studies with large sample size in diverse ethnicities are still needed to validate our findings.
Some limitations of this study need to be addressed. First of all, the study subjects were from hospitals and may not be representative of the general population. However, we applied rigorous epidemiological design by limiting factors of selecting subjects and used further statistical adjustment to minimize the potential biases. Second, the study population was restricted to Guangxi district, so it does not permit extrapolation of the results to other populations in other areas. Third, the sample size of this study is relatively small, especially for haplotype analysis, which may not have enough statistical power to explore the true association. Therefore, large population-based prospective studies with ethnically diverse populations are warranted to further elucidate the impact of VDR and DBP SNPs on HCC susceptibility.
In conclusion, our study suggested that the VDR rs2228570 and DBP rs7041 polymorphisms were associated with a significantly increased risk of HBV-related HCC in this Chinese population. This is the first report that the potentially functional SNP in the promoter region of VDR and DBP may contribute to HBV-related HCC susceptibility in the Chinese population. Due to the marginal significance, larger well-designed epidemiological studies with ethnically diverse populations and functional evaluations are warranted to confirm our findings.
Author Contributions
Conceived and designed the experiments: QP SL XQ. Performed the experiments: SY XL RL JW. Analyzed the data: RL ZC. Contributed reagents/materials/analysis tools: SL XQ. Wrote the paper: QP SY.
References
- 1. Jemal A, Bray F, Center MM, Ferlay J, Ward E, et al. (2011) Global cancer statistics. CA Cancer J Clin 61:69–90.
- 2. El-Serag HB, Rudolph KL (2007) Hepatocellular carcinoma: epidemiology and molecular carcinogenesis. Gastroenterology 132:2557–2576.
- 3. Kuper H, Tzonou A, Kaklamani E, Hsieh CC, Lagiou P, et al. (2000) Tobacco smoking, alcohol consumption and their interaction in the causation of hepatocellular carcinoma. Int J Cancer 85:498–502.
- 4. Mohapatra S, Saxena A, Gandhi G, Koner BC, Ray PC (2013) Vitamin D and VDR gene polymorphism (FokI) in epithelial ovarian cancer in Indian population. J Ovarian Res 6:37.
- 5. Bikle D (2009) Nonclassic actions of vitamin D. J Clin Endocrinol Metab 94:26–34.
- 6. Deeb KK, Trump DL, Johnson CS (2007) Vitamin D signalling pathways in cancer: potential for anticancer therapeutics. Nat Rev Cancer 7:684–700.
- 7. Zhou H, Xu C, Gu M (2009) Vitamin D receptor (VDR) gene polymorphisms and Graves' disease: a meta-analysis. Clin Endocrinol (Oxf) 70:938–945.
- 8. Dalessandri KM, Miike R, Wiencke JK, Farren G, Pugh TW, et al. (2012) Vitamin D receptor polymorphisms and breast cancer risk in a high-incidence population: a pilot study. J Am Coll Surg 215:652–657.
- 9. Shui IM, Mucci LA, Kraft P, Tamimi RM, Lindstrom S, et al. (2012) Vitamin D-related genetic variation, plasma vitamin D, and risk of lethal prostate cancer: a prospective nested case-control study. J Natl Cancer Inst 104:690–699.
- 10. Rasool S, Kadla SA, Rasool V, Qazi F, Khan T, et al. (2014) Role of the VDR Bsm I and Apa I polymorphisms in the risk of colorectal cancer in Kashmir. Oncol Res Treat 37:345–349.
- 11. Powe CE, Ricciardi C, Berg AH, Erdenesanaa D, Collerone G, et al. (2011) Vitamin D-binding protein modifies the vitamin D-bone mineral density relationship. J Bone Miner Res 26:1609–1616.
- 12. Anic GM, Weinstein SJ, Mondul AM, Mannisto S, Albanes D (2014) Serum vitamin D, vitamin D binding protein, and risk of colorectal cancer. PLoS One 9:e102966.
- 13. Sinotte M, Diorio C, Berube S, Pollak M, Brisson J (2009) Genetic polymorphisms of the vitamin D binding protein and plasma concentrations of 25-hydroxyvitamin D in premenopausal women. Am J Clin Nutr 89:634–640.
- 14. Anderson LN, Cotterchio M, Cole DE, Knight JA (2011) Vitamin D-related genetic variants, interactions with vitamin D exposure, and breast cancer risk among Caucasian women in Ontario. Cancer Epidemiol Biomarkers Prev 20:1708–1717.
- 15. Kidd LC, Paltoo DN, Wang S, Chen W, Akereyeni F, et al. (2005) Sequence variation within the 5' regulatory regions of the vitamin D binding protein and receptor genes and prostate cancer risk. Prostate 64:272–282.
- 16. Bai Y, Yu Y, Yu B, Ge J, Ji J, et al. (2009) Association of vitamin D receptor polymorphisms with the risk of prostate cancer in the Han population of Southern China. BMC Med Genet 10:125.
- 17. Zhao XZ, Yang BH, Yu GH, Liu SZ, Yuan ZY (2014) Polymorphisms in the vitamin D receptor (VDR) genes and skin cancer risk in European population: a meta-analysis. Arch Dermatol Res 306:545–553.
- 18. Kaabachi W, Kaabachi S, Rafrafi A, Amor AB, Tizaoui K, et al. (2014) Association of vitamin D receptor FokI and ApaI polymorphisms with lung cancer risk in Tunisian population. Mol Biol Rep 41:6545–6553.
- 19. Shi YY, He L (2005) SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Res 15:97–98.
- 20. Stephens M, Smith NJ, Donnelly P (2001) A new statistical method for haplotype reconstruction from population data. Am J Hum Genet 68:978–989.
- 21. Chen Y, Liu W, Sun T, Huang Y, Wang Y, et al. (2013) 1,25-Dihydroxyvitamin D promotes negative feedback regulation of TLR signaling via targeting microRNA-155-SOCS1 in macrophages. J Immunol 190:3687–3695.
- 22. Fu B, Wang H, Wang J, Barouhas I, Liu W, et al. (2013) Epigenetic regulation of BMP2 by 1,25-dihydroxyvitamin D3 through DNA methylation and histone modification. PLoS One 8:e61423.
- 23. Walters MR (1992) Newly identified actions of the vitamin D endocrine system. Endocr Rev 13:719–764.
- 24. Baker AR, McDonnell DP, Hughes M, Crisp TM, Mangelsdorf DJ, et al. (1988) Cloning and expression of full-length cDNA encoding human vitamin D receptor. Proc Natl Acad Sci U S A 85:3294–3298.
- 25. Slattery ML (2007) Vitamin D receptor gene (VDR) associations with cancer. Nutr Rev 65:S102–104.
- 26. Gu H, Wang X, Zheng L, Tang W, Dong C, et al. (2014) Vitamin D receptor gene polymorphisms and esophageal cancer risk in a Chinese population: a negative study. Med Oncol 31:827.
- 27. Beilfuss A, Sowa JP, Sydor S, Beste M, Bechmann LP, et al. (2014) Vitamin D counteracts fibrogenic TGF-beta signalling in human hepatic stellate cells both receptor-dependently and independently. Gut
- 28. Uitterlinden AG, Fang Y, Van Meurs JB, Pols HA, Van Leeuwen JP (2004) Genetics and biology of vitamin D receptor polymorphisms. Gene 338:143–156.
- 29. Arai H, Miyamoto K, Taketani Y, Yamamoto H, Iemori Y, et al. (1997) A vitamin D receptor gene polymorphism in the translation initiation codon: effect on protein activity and relation to bone mineral density in Japanese women. J Bone Miner Res 12:915–921.
- 30. Shafia S, Qasim I, Aziz SA, Bhat IA, Nisar S, et al. (2013) Role of vitamin D receptor (VDR) polymorphisms in susceptibility to multiple myeloma in ethnic Kashmiri population. Blood Cells Mol Dis 51:56–60.
- 31. Sarkissyan M, Wu Y, Chen Z, Mishra DK, Sarkissyan S, et al. (2014) Vitamin D receptor FokI gene polymorphisms may be associated with colorectal cancer among African American and Hispanic participants. Cancer 120:1387–1393.
- 32. Nykjaer A, Dragun D, Walther D, Vorum H, Jacobsen C, et al. (1999) An endocytic pathway essential for renal uptake and activation of the steroid 25-(OH) vitamin D3. Cell 96:507–515.
- 33. Chishimba L, Thickett DR, Stockley RA, Wood AM (2010) The vitamin D axis in the lung: a key role for vitamin D-binding protein. Thorax 65:456–462.
- 34. Pibiri F, Kittles RA, Sandler RS, Keku TO, Kupfer SS, et al. (2014) Genetic variation in vitamin D-related genes and risk of colorectal cancer in African Americans. Cancer Causes Control 25:561–570.
- 35. Wang W, Ingles SA, Torres-Mejia G, Stern MC, Stanczyk FZ, et al. (2014) Genetic variants and non-genetic factors predict circulating vitamin D levels in Hispanic and non-Hispanic White women: the Breast Cancer Health Disparities Study. Int J Mol Epidemiol Genet 5:31–46.
- 36. Zhou L, Zhang X, Chen X, Liu L, Lu C, et al. (2012) GC Glu416Asp and Thr420Lys polymorphisms contribute to gastrointestinal cancer susceptibility in a Chinese population. Int J Clin Exp Med 5:72–79.
- 37. Poynter JN, Jacobs ET, Figueiredo JC, Lee WH, Conti DV, et al. (2010) Genetic variation in the vitamin D receptor (VDR) and the vitamin D-binding protein (GC) and risk for colorectal cancer: results from the Colon Cancer Family Registry. Cancer Epidemiol Biomarkers Prev 19:525–536.
- 38. Seng KC, Seng CK (2008) The success of the genome-wide association approach: a brief story of a long struggle. Eur J Hum Genet 16:554–564.
- 39. Kaabachi W, Kaabachi S, Rafrafi A, Amor AB, Tizaoui K, et al. (2014) Association of vitamin D receptor FokI and ApaI polymorphisms with lung cancer risk in Tunisian population. Mol Biol Rep