Next Article in Journal
Controlled Release of Therapeutics from Thermoresponsive Nanogels: A Thermal Magnetic Resonance Feasibility Study
Previous Article in Journal
PARP Inhibitors in the Treatment of Early Breast Cancer: The Step Beyond?
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 12
Issue 6
10.3390/cancers12061379
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

The Role of Vitamin D Receptor Gene Polymorphisms in Colorectal Cancer Risk

by
Ippokratis Messaritakis
1,*,
Asimina Koulouridi
1,
Maria Sfakianaki
1,
Konstantinos Vogiatzoglou
1,
Nikolaos Gouvas
2,
Elias Athanasakis
3,
John Tsiaoussis
4,
Evangelos Xynos
5,
Dimitriοs Mavroudis
1,6,
Maria Tzardi
7 and
John Souglakos
1,6
1
Laboratory of Translational Oncology, Medical School, University of Crete, 70013 Heraklion, Greece
2
Medical School, University of Cyprus, 20537 Nicosia, Cyprus
3
Department of General Surgery, Heraklion University Hospital, 71110 Heraklion, Greece
4
Department of Anatomy, School of Medicine, University of Crete, 70013 Heraklion, Greece
5
Department of Surgery, Creta Interclinic Hospital of Heraklion, 71305 Heraklion, Greece
6
Department of Medical Oncology, University General Hospital of Heraklion, 71110 Heraklion, Greece
7
Laboratory of Pathology, University General Hospital of Heraklion, 70013 Heraklion, Greece
*
Author to whom correspondence should be addressed.
Cancers 2020, 12(6), 1379; https://0-doi-org.brum.beds.ac.uk/10.3390/cancers12061379
Submission received: 28 April 2020 / Revised: 25 May 2020 / Accepted: 25 May 2020 / Published: 27 May 2020
Browse Figures
Versions Notes

Abstract

:
Vitamin D deficiency has been associated with increased colorectal cancer (CRC) incidence risk and mortality. Vitamin D mediates its action through the binding of the vitamin D receptor (VDR), and polymorphisms of the VDR might explain these inverse associations. The aim of the study was the investigation of the relevance of rs731236; Thermus aquaticus I (TaqI), rs7975232; Acetobacter pasteurianus sub. pasteurianus I (ApaI), rs2228570; Flavobacterium okeanokoites I (FokI) and rs1544410, Bacillus stearothermophilus I (BsmI) polymorphisms of the VDR gene to colorectal carcinogenesis (CRC) and progression. Peripheral blood was obtained from 397 patients with early operable stage II/III (n = 202) and stage IV (n = 195) CRC. Moreover, samples from 100 healthy donors and 40 patients with adenomatous polyps were also included as control groups. Genotyping in the samples from patients and controls was performed using polymerase chain reaction-restriction fragment length polymorphisms (PCR-RFLP). A significant association was revealed between all four polymorphisms and cancer. Individuals with homozygous mutant (tt, aa, ff or bb) genotypes were more susceptible to the disease (p < 0.001). All of the mutant genotypes detected were also significantly associated with stage IV (p < 0.001), leading to significantly decreased survival (p < 0.001). Moreover, all four polymorphisms were significantly associated with KRAS (Kirsten ras oncogene) mutations and Toll-like receptor (TLR2, TLR4 and TLR9) genetic variants. In multivariate analysis, tt, aa and ff genotypes emerged as independent factors associated with decreased overall survival (OS) (p = 0.001, p < 0.001 and p = 0.001, respectively). The detection of higher frequencies of the VDR polymorphisms in CRC patients highlights the role of these polymorphisms in cancer development and progression.

1. Introduction

Colorectal cancer (CRC) is the third leading cause of cancer both in men and women [1]. CRC is of major importance to public health, with increased mortality rates worldwide and accounting for 9% of all cancers [1]. CRC is a multifactorial disease and involves complex interactions between environmental and genetic factors [2,3]. However, understanding all of the mechanisms and the associations between these factors is not an easy task. In the past, it has been hypothesized that higher incidence rates of CRC in areas with low sunlight exposure might be attributable to lower levels of vitamin D [4]. Indeed, it has been reported that vitamin D deficiency is a common phenomenon in Saudi Arabia, especially among women [5]. Nowadays, a large number of candidates’ genes have been identified as responsible for their potential role in tumorigenesis.
The importance of vitamin D for bone health is well established, but lately its role beyond the skeletal system has gained scientific attention [6]. Vitamin D regulates cellular differentiation and proliferation in normal and malignant tissues, regulates proliferation, apoptosis and cell adhesion in tumor cells and modifies tumor angiogenesis, invasion and metastasis, along with reducing oxidative DNA damage [7]. Vitamin D deficiency has been associated with various cancer types [8,9], whereas, increased vitamin D serum levels play a role in decreased colorectal adenoma risk [9]. Vitamin D mediates its action by binding to the vitamin D receptor (VDR), a member of the nuclear receptor superfamily, and is expressed on various cell types, including colorectal epithelial cells, enabling the transactivation of target genes [10,11]. Thus, the VDR gene has been implicated in CRC. Over 60 single nucleotide polymorphisms (SNPs) of the VDR gene, located in the promoter region in exons 2–9, both in their proximity and in the 3’-UTR (3’-untranslated) region, have been studied in relation to cancer occurrence and prognosis [12,13]. However, only a few of them are potentially functional and affect the expression of the VDR gene in relation to CRC risk. These include TaqI (rs731236; Thermus aquaticus I), located in exon 9 [9,14], ApaI (rs7975232; Acetobacter pasteurianus sub. pasteurianus I) and BsmI (rs1544410, Bacillus stearothermophilus I), located in the intron between exons 8 and 9 [14,15,16] and FokI (rs2228570; Flavobacterium okeanokoites I), located in exon 2 [15,17].
To this end, we aimed to investigate four single nucleotide polymorphisms (SNPs)—TaqI, ApaI, FokI and BsmI—of the VDR gene within patients with sporadic CRC, for the first time in the Greek population, and to evaluate their association with the risk of cancer development and progression. The selection of these SNPs was based on the common VDR SNPs sites examined in other populations in previous genetic epidemiological studies. In addition, the correlation of the expression of these molecules and previously genotyped Toll-like receptor (TLR) variants (TLR2: 196-to-174 del; TLR4: Asp299Gly and Thr399Ile; TLR9: T1237C and T1486C) in the same patients’ and controls’ samples [18] was also investigated.

2. Results

2.1. Patients’ Demographics and Molecular Characteristics

From 09/2003 to 11/2013, 397 patients were recruited in the study and presented newly diagnosed CRC and histologically documented disease. The patients’ characteristics are listed in Table 1. The median age of the patients was 65 years, 246 (62.0%) were males, 202 (50.9%) were of stage II/III, 372 (93.7%) had PS-ECOG (Performance Status according to the Eastern Cooperative Oncology Group) 0–1, 205 (47.4%) had a high tumor grade and 279 (70.3%) had a colon/sigmoid tumor location. Moreover, 230 (58.7%) patients had early operable stage II/III disease, 64 (27.8%) of whom relapsed, whereas 223 (56.2%) patients had stage IV disease, 197 (88.3%) of whom relapsed, as shown in Table 1 and Table S1.
The VDR gene TaqI t allele (silent T → C transition in exon 9), ApaI a allele (T → G transition in intron 8), FokI f allele (C → T transition at the junction of intron 1 and exon 2) and BsmI b allele (G → A transition in intron 8) genotypes and allele frequencies were investigated in all 397 patients and are shown in Table 1 and Table S1; whereas KRAS (Kirsten ras oncogene) mutations were investigated in 245 patients due to there being no sample availability, shown in Table 1 and Table S1.

2.2. Analysis of the VDR Gene Polymorphisms

The VDR gene TaqI, ApaI, FokI and BsmI polymorphisms amplification products were expected to be 740 bp, 740 bp, 265 bp and 825 bp, respectively. The PCR products were digested by the TaqI, ApaI, FokI and BsmI enzymes, respectively. Following electrophoresis, 143 (36.0%), 129 (32.5%) and 125 (36.5%) patients presented the homozygous mutant (tt), the heterozygous (Tt) and the wild type (TT) genotypes, respectively, for TaqI polymorphisms, as shown in Figure 1A, Table 1 and Table S1. Similarly, 130 (32.7%), 122 (30.7%) and 145 (36.5%) patients presented the homozygous mutant (aa), the heterozygous (Aa) and the wild type (AA) genotypes, respectively, for ApaI polymorphisms, as shown in Figure 1B, Table 1 and Table S1. Moreover, 128 (32.2%), 147 (37.0%) and 122 (30.7%) patients presented the homozygous mutant (ff), the heterozygous (Ff) and the wild type (FF) genotypes, respectively, for FokI polymorphisms, as shown in Figure 1C, Table 1 and Table S1. Finally, 131 (33.0%), 191 (48.1%) and 75 (18.9%) patients presented the homozygous mutant (bb), the heterozygous (Bb) and the wild type (BB) genotypes, respectively, for BsmI polymorphisms, as shown in Figure 1D, Table 1 and Table S1.
The results showed that the allelic frequencies of all four polymorphisms were significantly associated with the patient group compared to the control groups (TaqI, p > 0.001; ApaI, p > 0.001; FokI, p > 0.001; BsmI, p > 0.001), thus highlighting the role of these polymorphisms in the disease. More specifically, it was shown that both the healthy donors and the adenomatous polyps controls presented mainly the wild type and the heterozygous genotypes, and this was observed in all four (TaqI, ApaI, FokI and BsmI) polymorphisms, as shown in Table 2.

2.3. Association of VDR Variants and Disease Stage

Table 2 shows the association observed between the VDR polymorphisms and disease stage of the patients. The tt genotype was more prevalent in stage IV patients, whereas the Tt and TT genotypes were mostly seen in stage II/III patients (62.4% vs. 12.0%, 28.6% vs. 36.1% and 9.0% vs. 51.9%, respectively; p < 0.001), as shown Table 2 and Table S1. The aa genotype was more prevalent in stage IV patients, whereas the Aa and AA genotypes were also mostly prevalent in stage II/III patients (61.4% vs. 6.7%, 25.4% vs. 35.6% and 13.2% vs. 57.7%, respectively; p < 0.001), as shown in Table 2 and Table S1. Moreover, the ff genotype was also more prevalent in stage IV patients, whereas the Ff and the FF genotypes were mostly met in stage II/III patients (56.6% vs. 10.1%, 30.7% vs. 42.8% and 12.7% vs. 47.1%, respectively; p < 0.001), as shown in Table 2 and Table S1. Similarly, the bb genotype was more prevalent in stage IV patients, whereas the Bb and BB genotypes were mostly prevalent in stage II/III patients (49.7% vs. 17.8%, 39.7% vs. 55.8% and 10.6% vs. 26.4%, respectively; p < 0.001), as shown in Table 2 and Table S1.

2.4. Correlation of VDR and Toll-Like Receptor Variants

The correlation between the different VDR and TLR variants was analyzed and is presented in Table 3. When analyzing the whole group of patients, a statistically significant coexistence was observed both between all the different VDR combinations and between the VDR and TLR genotype combinations, shown in Table 3 and Table S1. When patients were analyzed according to their disease stage, it was observed that stage II/III patients had a significant coexistence of only: TaqI-FokI, p = 0.005; TaqI-BsmI, p < 0.001; ApaI-BsmI, p = 0.014; FokI-BsmI, p < 0.001; FokI-Asp299Gly, p = 0.035; FokI-Thr399Ile, p = 0.035, as shown in Table 3 and Table S1. Similarly, in stage IV patients, a significant coexistence was observed in: TaqI-ApaI, p < 0.001; TaqI-BsmI, p < 0.001; ApaI-BsmI, p < 0.001; FokI-ApaI, p = 0.001; FokI-BsmI, p < 0.001; TaqI-T1237C, p = 0.042; TaqI-T1486C, p = 0.042; FokI-Asp299Gly, p = 0.037; FokI-Thr399Ile, p = 0.037, as shown in Table 3 and Table S1.

2.5. Association of TaqI, ApaI, FokI and BsmI Variants and KRAS Status

The VDR variants and KRAS statuses of the CRC patients presented a significant association, as shown in Table 4. The tt genotype was more prevalent in patients with a mutant KRAS status, whereas the Tt and the TT alleles were mostly seen in KRAS wild type patients (53.8% vs. 39.0%, 24.0% vs. 34.8% and 22.1% vs. 26.2%, respectively, p = 0.041). The aa genotype was also more frequent in KRAS mutants, whereas Aa and AA genotypes were more frequent in patients with KRAS wild type status (55.8% vs. 36.9%, 19.2% vs. 34.8% and 25.0% vs. 28.4%, respectively; p = 0.007). Similarly, the ff genotype was more prevalent in KRAS mutants, whereas the Ff and FF genotypes were more frequent in KRAS wild type patients (59.6% vs. 27.7%, 25.0% vs. 42.6% and 15.4% vs. 29.8%, respectively; p < 0.001). Finally, the bb genotype was also more frequent in KRAS mutants, whereas Bb and BB genotypes were more frequent in patients with KRAS wild type status (65.4% vs. 27.7%, 24.0% vs. 55.3% and 10.6% vs. 17.0%, respectively; p < 0.001), as shown in Table 4 and Table S1.

2.6. Association of VDR Variants and Clinical Outcome

Sixty-four (27.8%) adjuvant and 197 (49.7%) metastatic patients presented a disease progression following their adjuvant and first-line treatment, respectively, as shown in Table S1. The median disease-free survival (DFS) was 19 months (95% confidence interval (CI): 15.5–22.5) and the median overall survival (OS) was 155 months (95% CI: 59.1–250.9), respectively for stage II/III patients. According to the presence of different VDR genotypes, no significant differences were observed in DFS, whereas only a significant shorter OS in patients with the aa genotype (p < 0.001) was observed, as shown in Figure 2A.
For the case of stage IV patients, the median progression-free survival (PFS) was 8 months (95% CI: 7.1–8.9) and the median OS was 31 months (95% CI: 25.2–36.8), respectively, as shown in Table S1. Again, there was no difference in PFS, according to the presence of different VDR genotypes, whereas only a significant decrease in OS in patients with TaqI homozygous mutant or heterozygous alleles (p = 0.037) was observed, as shown in Figure 2B. The analysis of all patients presented a median OS of 75 months (95% CI: 56.8–93.2) prevailing a significantly decreased OS in patients with TaqI, ApaI, FokI and/or BsmI homozygous mutant alleles (p < 0.001, p < 0.001, p < 0.001 and p < 0.001, respectively), as shown in Figure 2C–F and Table S1.

2.7. Univariate and Multivariate Analysis Cox Regression Analysis

Univariate analysis revealed that PS (ECOG) and TLR2 polymorphisms were significantly associated with a shorter PFS and PS (ECOG) tumor grade, and all TLR and VDR polymorphisms were significantly associated with shorter OS, as shown in Table 5. In multivariate analysis, adjusting for these factors, PS (ECOG) and TLR2 polymorphisms emerged as independent factors associated with decreased PFS (HR: 1.6, 95% CI: 1.0–2.6, p = 0.04 and HR: 2.2, 95% CI: 1.2–4.0, p = 0.013, respectively). Moreover, PS (ECOG), tumor grade, TLR2 196-to-174 del, TaqI, ApaI and FokI variants emerged as independent factors associated with decreased OS (HR: 3.7, 95% CI: 2.4–5.9, p < 0.001; HR: 1.8, 95% CI: 0.7–5.1, p < 0.001; HR: 2.1, 95% CI: 1.3–3.4, p = 0.003; HR: 1.3, 95% CI: 1.0–1.5, p = 0.029; HR: 1.6, 95% CI: 1.3–2.0, p < 0.001; HR: 1.4, 95% CI: 1.1–1.8, p = 0.005), as shown in Table 5.

3. Discussion

Undoubtedly, the risk of sporadic CRC has been linked to environmental and genetic factors [2,19] and more recently to gut microbiota. The microbiota influences both human health and disease by affecting the development of the host immune system and by maintaining homeostasis to influence diseases and allergies that cannot simply be parsed into strict pathogenesis and commensalism [20,21]. Chronic infection and inflammation are the most important epigenetic factors contributing to tumorigenesis and tumor progression [22]. Moreover, vitamin D deficiency has been implicated, among other diseases and metabolic syndromes, in human malignancies [23,24]. Since vitamin D mediates its action by binding to the VDR gene, the suboptimal responsiveness of the VDR can be manifested as vitamin D deficiency. Furthermore, the VDR gene, as part of innate immunity, is responsible for the prevention and elimination of infection and the determination of the gut microbiome [25,26,27]. Thus, polymorphisms in the human VDR gene may play an important role in the structure of the gut microbiome. Interestingly, it has been reported previously that VDR conditional knockout (vdrΔIEC) in the intestinal epithelial or low intestinal VDR protein levels may lead to dysbiosis [25] and reduced autophagy, accompanied by a reduction in ATG16L1 (autophagy-related 16 like 1), an inflammatory bowel disease risk gene [27]; whereas, the absence of intestinal VDR leads to a susceptibility to colon cancer via reducing JAK/STAT (Janus kinases/signal transducer and activator of transcription proteins) signaling, which is a pathway with a critical role in intestinal and microbial homeostasis and in dampening inflammatory responses [27]. Therefore, the vitamin D/VDR pathway may significantly influence homeostasis, signaling between the microbiota and host in intestinal inflammation and tumorigenesis [27]. It is worth mentioning that, despite the fact that the VDR gene polymorphisms (TaqI, ApaI, FokI and BsmI) are considered to be non-functional, they might be linked to other functional polymorphisms elsewhere in the VDR gene, thus participating in a more complex gene network, enhancing or inhibiting the expression of VDR target genes. Such VDR gene polymorphisms are likely to affect transcriptional regulation, mRNA stability or protein translational efficiency, thus affecting the structure and functioning of VDR protein [9,28,29,30]. It has also been suggested that VDR gene polymorphisms define differential transcriptional VDR activity or mRNA stability in vitro [31]. Moreover, vitamin D is considered as an immunomodulatory molecule, which can modulate cytokine responses through T cells, thus representing an important link between TLR activation and innate immunity against microorganisms [32,33]. Others have described the TLR activation of human macrophages upregulating VDR gene expression, leading to the induction of cathelicidin and the consequent killing of intracellular Mycobacterium tuberculosis [34]. Moreover, the authors further analyzed Cyp27B1 (cytochrome P450 family 27 subfamily B member 1), which catalyzes the conversion of the inactive provitamin D3 hormone into its active form, and it was found to be significantly upregulated. Furthermore, the VDR and Cyp27B1 in monocytes and macrophages were also found to be upregulated. Thus, the authors concluded that the TLR induces the upregulation of the VDR and Cyp27B1 gene expressions in such cell types [34]. They also performed functional tests on the VDR and it was demonstrated that the VDR is functional in primary human monocytes and, when activated, it triggers the stimulation of antimicrobial peptides. Finally, the authors demonstrate that TLR activation in monocytes might lead to the activation of a microbicidal pathway, which depends on the production and action of vitamin D through VDR. Overall, the authors provided a potential explanation of how vitamin D may act as a key link between TLR activation and antimicrobial responses in innate immunity [34]. The aim of the current study was to evaluate the detection of VDR (TaqI, ApaI, FokI and BsmI) polymorphisms in adjuvant and metastatic CRC patients.
A number of studies on VDR polymorphisms have been performed in patients with sporadic CRC. Some studies reported contradictory results between various VDR genetic variants and CRC [35,36], or even no association [13,37], probably due to limitations such as the recruitment of patients with different characteristics, such as ethnicity, among these studies and the small sample size. However, many of these studies presented significant associations [35,38]. Indeed, in a meta-analysis conducted by Serrano et al., [39] the authors presented a significant increased risk for CRC in the presence of tt and aa genotypes, with strong frequency variations present among different ethnic groups. Moreover, Gandini et al., [24] reviewed 79 studies with more than 52,000 cancer cases, including CRC patients, and 62,000 controls, and demonstrated that BsmI, FokI and TaqI polymorphisms are associated with CRC. Additional studies reported that the bb genotype of the BsmI polymorphism are more susceptible to CRC, whereas BB, Bb, TtFf and TTFf genotypes are significantly associated with a decreased risk of CRC [14,35,40,41]. Our results are in accordance with these findings. In fact, we demonstrate a higher frequency of the tt, aa, ff and bb genotypes in CRC patients compared to the control groups, thus highlighting the role of these polymorphisms in colorectal carcinogenesis. Moreover, higher frequencies of the tt, aa, ff and bb genotypes were detected in metastatic CRC patients compared to stage II/III patients, emphasizing the role of these polymorphisms in CRC progression and in patients’ overall survival.
The Ras/MAPK (Ras/mitogen activated protein kinases) pathway and its continuous activation, due to mutations presented in codon 12 of the KRAS gene, plays an important role in treatment resistance in patients with various carcinomas, including CRC [42]. To this end, we also aimed to associate the frequency of VDR polymorphisms in patients with different KRAS statuses. Despite the fact that approximately 40% of the enrolled patients were not evaluated for their KRAS statuses, a significant association was demonstrated between the tt, aa, ff and bb genotypes with the KRAS mutant patients.
The TLR pathway increases the risk of colitis-associated CRC due to commensal gut microbiota [43,44]. TLRs play an important role in immunity and are expressed in various cell types, including tumor cells [45,46]. Since TLR polymorphisms have been associated with changes in susceptibility to many diseases, including cancers [47] and TLRs promote the survival of cancer cells [48], we also correlated the coexistence of previously genotyped TLR (TLR2, TLR4 and TLR9) polymorphisms [18] with the VDR polymorphisms. It has been previously reported that human TLRs are considered to be regulators in vitamin D/VDR signaling. In humans, when a pathogen is detected by TLRs, gene expressions of VDR and Cyp27B1 are induced [34,49].
Herein, we demonstrated that both TLR [18] and VDR polymorphisms are associated with an increased risk of CRC development and progression, with an impact on patients’ survival, and also demonstrated a significant correlation between TLR and VDR gene polymorphisms.
To our knowledge, this is the first time VDR gene polymorphisms have been investigated in the Greek population to evaluate their role in CRC risk and patients’ survival. The strength of the present study includes careful clinical and epidemiological data collection in combination with genotyping. Such data collection allowed us to investigate gene-to-gene or other interactions. Moreover, the relatively large sample size and the 10-year follow-up of the enrolled patients provided enough strength to discriminate significant interactions that have an impact on CRC risk and the survival of the patients. An additional strength is that the enrolled cases and controls were from the same ethnicity and were age and gender matched. A limitation of the study is the lack of vitamin D measurements in the plasma or sera of the enrolled patients and controls, due to the retrospective nature of the study. However, despite this limitation, the results of this study provide significant information regarding the association between TaqI, ApaI and FokI polymorphisms and susceptibility to CRC and its progression.

4. Materials and Methods

4.1. Patients’ Population

In total, 397 patients with colon adenocarcinoma were enrolled in the study, between 2003 and 2013, from the Department of Medical Oncology, University Hospital of Heraklion.

4.2. Ethics Approval and Consent to Participate

The study has been approved by the Ethics Committee/Institutional Review Board of the University Hospital of Heraklion (Number 7302/19-8-2009), and all patients signed a written informed consent form for their participation. All of the procedures performed were in accordance with the ethical standards of the institutional and/or national research committee and the 1964 Helsinki declaration, and its later amendments or comparable ethical standards

4.3. Blood and Tissue Samples from Control Groups

In parallel to patients’ samples, blood samples and formalin-fixed paraffin embedded (FFPE) tissues were also obtained from 100 healthy blood donors and from 40 patients with colon adenomas in the absence of CRC disease, respectively, which were used as controls in the study.

4.4. Genomic DNA Extraction

Peripheral blood mononuclear cells (PBMC) from all individuals (patients and healthy donors) were obtained using the Ficoll–Hypaque density gradient (d = 1077 g/mL; Sigma-Aldrich, GmbH, Darmstadt, Germany), as described previously [18]. Representative formalin fixed paraffin embedded (FFPE) specimens from the primary tumor were examined by an experienced pathologist and the appropriate area for microdissection was defined. Microdissection and malignant cells collection were performed using a piezoelectric microdissector (Eppendorf, Hamburg, Germany), as described previously [50].
The DNA extraction of all samples was performed using the MasterPure™ Complete DNA and RNA Purification Kit (Epicenter, Madison, WI, USA), according to the manufacturer’s instructions. NanoDrop ND-1000, version 3.3 (ThermoFisher Scientific, Waltham, MA, USA) was used for DNA quantification.

4.5. VDR and TLR Genotyping

For the genotyping of the SNPs at the TaqI, ApaI, FokI and BsmI positions of the VDR gene, polymerase chain reaction (PCR) and restriction fragment length polymorphism (RFLP) methods were used. The sequences of the primers used for the PCR amplification of the fragments are provided in Table 6. Allele types, SNP reference numbers and PCR conditions for all the analyzed polymorphisms are shown in Table 7.
TaqI (Minotech Biotechnology, IMBB-FORTH, Heraklion, Greece), ApaI (ThermoFisher Scientific, MA, USA), FokI (ThermoFisher Scientific) and BsmI (ThermoFisher Scientific) restriction enzymes were used to digest the amplified products of the VDR gene, according to manufacturer’s instructions. Briefly, 10 μL of each related PCR product was mixed with 1 μL of each restriction enzyme and 2 μL of 10× buffers. Diethyl pyrocarbonate (DEPC) treated water was added to a final volume of 20 μL (for TaqI) or 30 μL (for ApaI, FokI and BsmI). After incubation at 65 °C for 15 min (TaqI) or at 37 °C for 5 min (ApaI, FokI and BsmI), the restriction fragments were separated by electrophoresis on a 2% agarose gel, stained with Sybr Safe DNA Gel Stain (ThermoFisher Scientific), and were visualized with the AlphaImager ultraviolet transilluminator (Alpha Innotech Corp., San Leandro, CA, USA). The usual nomenclature for restriction fragment length polymorphism alleles was used in this study [51,52]. The lowercase (t, a, f, b) alleles represent the presence of the restriction site and the uppercase alleles (T, A, F, B) represent the absence of the restriction site.
Accordingly, TLR genotyping was performed as previously described by our group [18]. In brief, the determination of TLR2 196-to-174 Ins/Del polymorphism was performed by PCR, whereas the determinations of TLR4 and TLR9 were performed by PCR-RFLP. The PCR and PCR-RFLP conditions and primers sets that were used have been previously reported by our group [18].

4.6. KRAS Mutational Analysis

KRAS mutational analysis was performed by Sanger sequencing after the PCR amplification of KRAS exon 2. The PCR conditions and the primers set that was used have been previously reported by our group [53].

4.7. Study Design and Statistics

The current study is a retrospective, single institution study which aimed to investigate the VDR gene polymorphisms in CRC patients before the initiation of any treatment. Disease-free survival (DFS), progression-free survival (PFS) and overall survival (OS) were calculated as previously described [18]. In brief, DFS was calculated from the date of surgery to the date of disease recurrence, PFS was calculated from the date of diagnosis to documented disease progression or death from any cause, and OS was calculated from the date of diagnosis to the date of death, from any cause. Laboratory analysis was carried out blind to clinical data, and statistical analysis was based on contingency tables, including the calculations of hazard ratios (HR) and 95% CI, as previously described [18]. Statistical significance was set at p = 0.05.

5. Conclusions

In conclusion, the results of the present study highlight the significant role of VDR polymorphisms in carcinogenesis, disease progression and patients’ survival. Our data also showed a correlation between TLR and VDR expression and an increased impact on patients’ survival. Based on the present results, therapies targeting the activity of VDRs, including the modulation of the TLR/VDR pathways, might provide new approaches to the management of CRC.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6694/12/6/1379/s1, Table S1: Research raw data.

Author Contributions

Conceptualization, I.M. and J.S.; methodology, I.M., K.V. and M.S.; validation, I.M., M.S. and M.T.; data curation I.M., A.K., N.G., E.A., J.T., E.X., D.M., M.T. and J.S.; formal analysis, I.M.; supervision, I.M. and J.S.; writing—original draft preparation, I.M.; writing—review and editing, I.M and J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was fully supported by research grants from the Hellenic Society of Medical Oncology (HeSMO) and the Gastrointestinal Cancer Study Group (GIC-SG). The funders had no role in the design of the study, in the collection, analyses or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer Statistics, 2017. CA Cancer J. Clin. 2017, 67, 7–30. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Aran, V.; Victorino, A.P.; Thuler, L.C.; Ferreira, C.G. Colorectal Cancer: Epidemiology, Disease Mechanisms and Interventions to Reduce Onset and Mortality. Clin. Colorectal Cancer 2016, 15, 195–203. [Google Scholar] [CrossRef] [PubMed]
  3. Valle, L. Genetic predisposition to colorectal cancer: where we stand and future perspectives. World J. Gastroenterol. 2014, 20, 9828–9849. [Google Scholar] [CrossRef] [PubMed]
  4. Garland, C.F.; Garland, F.C. Do sunlight and vitamin D reduce the likelihood of colon cancer? Int. J. Epidemiol. 1980, 9, 227–231. [Google Scholar] [CrossRef]
  5. Alsanea, N.; Abduljabbar, A.S.; Alhomoud, S.; Ashari, L.H.; Hibbert, D.; Bazarbashi, S. Colorectal cancer in Saudi Arabia: incidence, survival, demographics and implications for national policies. Ann. Saudi Med. 2015, 35, 196–202. [Google Scholar] [CrossRef]
  6. Theodoratou, E.; Tzoulaki, I.; Zgaga, L.; Ioannidis, J.P. Vitamin D and multiple health outcomes: umbrella review of systematic reviews and meta-analyses of observational studies and randomised trials. BMJ 2014, 348, g2035. [Google Scholar] [CrossRef] [Green Version]
  7. Krajewski, W.; Dziegala, M.; Kolodziej, A.; Dembowski, J.; Zdrojowy, R. Vitamin D and urological cancers. Cent. Eur. J. Urol. 2016, 69, 139–147. [Google Scholar] [CrossRef] [Green Version]
  8. Lefkowitz, E.S.; Garland, C.F. Sunlight, vitamin D, and ovarian cancer mortality rates in US women. Int. J. Epidemiol. 1994, 23, 1133–1136. [Google Scholar] [CrossRef]
  9. Uitterlinden, A.G.; Fang, Y.; Van Meurs, J.B.; Pols, H.A.; Van Leeuwen, J.P. Genetics and biology of vitamin D receptor polymorphisms. Gene 2004, 338, 143–156. [Google Scholar] [CrossRef] [Green Version]
  10. Zhang, K.; Song, L. Association between vitamin D receptor gene polymorphisms and breast cancer risk: a meta-analysis of 39 studies. PLoS ONE 2014, 9, e96125. [Google Scholar] [CrossRef]
  11. Bouillon, R.; Van Schoor, N.M.; Gielen, E.; Boonen, S.; Mathieu, C.; Vanderschueren, D.; Lips, P. Optimal vitamin D status: A critical analysis on the basis of evidence-based medicine. J. Clin. Endocrinol. Metab. 2013, 98, E1283–E1304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Fuhrman, B.J.; Freedman, D.M.; Bhatti, P.; Doody, M.M.; Fu, Y.P.; Chang, S.C.; Linet, M.S.; Sigurdson, A.J. Sunlight, polymorphisms of vitamin D-related genes and risk of breast cancer. Anticancer Res. 2013, 33, 543–551. [Google Scholar] [PubMed]
  13. Perna, L.; Hoffmeister, M.; Schottker, B.; Arndt, V.; Haug, U.; Holleczek, B.; Burwinkel, B.; Ordonez-Mena, J.M.; Brenner, H. Vitamin D receptor polymorphism and colorectal cancer-specific and all-cause mortality. Cancer Epidemiol. 2013, 37, 905–907. [Google Scholar] [CrossRef] [PubMed]
  14. Laczmanska, I.; Laczmanski, L.; Bebenek, M.; Karpinski, P.; Czemarmazowicz, H.; Ramsey, D.; Milewicz, A.; Sasiadek, M.M. Vitamin D receptor gene polymorphisms in relation to the risk of colorectal cancer in the Polish population. Tumour Biol. 2014, 35, 12397–12401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Sarkissyan, M.; Wu, Y.; Chen, Z.; Mishra, D.K.; Sarkissyan, S.; Giannikopoulos, I.; Vadgama, J.V. Vitamin D receptor FokI gene polymorphisms may be associated with colorectal cancer among African American and Hispanic participants. Cancer 2014, 120, 1387–1393. [Google Scholar] [CrossRef] [Green Version]
  16. Bai, Y.H.; Lu, H.; Hong, D.; Lin, C.C.; Yu, Z.; Chen, B.C. Vitamin D receptor gene polymorphisms and colorectal cancer risk: a systematic meta-analysis. World J. Gastroenterol. 2012, 18, 1672–1679. [Google Scholar] [CrossRef]
  17. Jenab, M.; McKay, J.; Bueno-de-Mesquita, H.B.; van Duijnhoven, F.J.B.; Ferrari, P.; Slimani, N.; Jansen, E.H.J.M.; Pischon, T.; Rinaldi, S.; Tjonneland, A.; et al. Vitamin D receptor and calcium sensing receptor polymorphisms and the risk of colorectal cancer in European populations. Cancer Epidemiol. Biomarkers Prev. 2009, 18, 2485–2491. [Google Scholar] [CrossRef] [Green Version]
  18. Messaritakis, I.; Stogiannitsi, M.; Koulouridi, A.; Sfakianaki, M.; Voutsina, A.; Sotiriou, A.; Athanasakis, E.; Xynos, E.; Mavroudis, D.; Tzardi, M.; et al. Evaluation of the detection of Toll-like receptors (TLRs) in cancer development and progression in patients with colorectal cancer. PLoS ONE 2018, 13, e0197327. [Google Scholar] [CrossRef]
  19. De la Chapelle, A. Genetic predisposition to colorectal cancer. Nat. Rev. Cancer 2004, 4, 769–780. [Google Scholar] [CrossRef]
  20. Cho, I.; Blaser, M.J. The human microbiome: at the interface of health and disease. Nat. Rev. Genet. 2012, 13, 260–270. [Google Scholar] [CrossRef] [Green Version]
  21. Mandal, P. Molecular mechanistic pathway of colorectal carcinogenesis associated with intestinal microbiota. Anaerobe 2018, 49, 63–70. [Google Scholar] [CrossRef] [PubMed]
  22. Balkwill, F.; Coussens, L.M. Cancer: an inflammatory link. Nature 2004, 431, 405–406. [Google Scholar] [CrossRef] [PubMed]
  23. Holick, M.F. Vitamin D deficiency. N. Engl. J. Med. 2007, 357, 266–281. [Google Scholar] [CrossRef] [PubMed]
  24. Gandini, S.; Gnagnarella, P.; Serrano, D.; Pasquali, E.; Raimondi, S. Vitamin D receptor polymorphisms and cancer. Adv. Exp. Med. Biol. 2014, 810, 69–105. [Google Scholar] [PubMed]
  25. Wu, S.; Zhang, Y.G.; Lu, R.; Xia, Y.; Zhou, D.; Petrof, E.O.; Claud, E.C.; Chen, D.; Chang, E.B.; Carmeliet, G.; et al. Intestinal epithelial vitamin D receptor deletion leads to defective autophagy in colitis. Gut 2015, 64, 1082–1094. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, J.; Thingholm, L.B.; Skieceviciene, J.; Rausch, P.; Kummen, M.; Hov, J.R.; Degenhardt, F.; Heinsen, F.A.; Ruhlemann, M.C.; Szymczak, S.; et al. Genome-wide association analysis identifies variation in vitamin D receptor and other host factors influencing the gut microbiota. Nat. Genet. 2016, 48, 1396–1406. [Google Scholar] [CrossRef]
  27. Sun, J. The Role of Vitamin D and Vitamin D Receptors in Colon Cancer. Clin. Transl. Gastroenterol. 2017, 8, e103. [Google Scholar] [CrossRef]
  28. Papadopoulou, A.; Kouis, P.; Middleton, N.; Kolokotroni, O.; Karpathios, T.; Nicolaidou, P.; Yiallouros, P.K. Association of vitamin D receptor gene polymorphisms and vitamin D levels with asthma and atopy in Cypriot adolescents: a case-control study. Multidiscip. Respir. Med. 2015, 10, 26. [Google Scholar] [CrossRef] [Green Version]
  29. Bellamy, R.; Ruwende, C.; Corrah, T.; McAdam, K.P.; Thursz, M.; Whittle, H.C.; Hill, A.V. Tuberculosis and chronic hepatitis B virus infection in Africans and variation in the vitamin D receptor gene. J. Infect. Dis. 1999, 179, 721–724. [Google Scholar] [CrossRef] [Green Version]
  30. Uitterlinden, A.G.; Fang, Y.; van Meurs, J.B.; van Leeuwen, H.; Pols, H.A. Vitamin D receptor gene polymorphisms in relation to Vitamin D related disease states. J. Steroid Biochem. Mol. Biol. 2004, 89-90, 187–193. [Google Scholar] [CrossRef]
  31. Morrison, N.A.; Qi, J.C.; Tokita, A.; Kelly, P.J.; Crofts, L.; Nguyen, T.V.; Sambrook, P.N.; Eisman, J.A. Prediction of bone density from vitamin D receptor alleles. Nature 1994, 367, 284–287. [Google Scholar] [CrossRef] [PubMed]
  32. Arji, N.; Busson, M.; Iraqi, G.; Bourkadi, J.E.; Benjouad, A.; Bouayad, A.; Mariaselvam, C.; Salah, S.; Fortier, C.; Amokrane, K.; et al. Genetic diversity of TLR2, TLR4, and VDR loci and pulmonary tuberculosis in Moroccan patients. J. Infect. Dev. Ctries. 2014, 8, 430–440. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Meyer, V.; Saccone, D.S.; Tugizimana, F.; Asani, F.F.; Jeffery, T.J.; Bornman, L. Methylation of the Vitamin D Receptor (VDR) Gene, Together with Genetic Variation, Race, and Environment Influence the Signaling Efficacy of the Toll-Like Receptor 2/1-VDR Pathway. Front. Immunol. 2017, 8, 1048. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Liu, P.T.; Stenger, S.; Li, H.; Wenzel, L.; Tan, B.H.; Krutzik, S.R.; Ochoa, M.T.; Schauber, J.; Wu, K.; Meinken, C.; et al. Toll-like receptor triggering of a vitamin D-mediated human antimicrobial response. Science 2006, 311, 1770–1773. [Google Scholar] [CrossRef] [PubMed]
  35. Li, C.; Li, Y.; Gao, L.B.; Wang, Y.Y.; Zhou, B.; Lv, M.L.; Lu, H.M.; Zhang, L. Vitamin D receptor gene polymorphisms and the risk of colorectal cancer in a Chinese population. Dig. Dis. Sci. 2009, 54, 634–639. [Google Scholar] [CrossRef] [PubMed]
  36. Yaylim-Eraltan, I.; Arzu Ergen, H.; Arikan, S.; Okay, E.; Ozturk, O.; Bayrak, S.; Isbir, T. Investigation of the VDR gene polymorphisms association with susceptibility to colorectal cancer. Cell Biochem. Funct. 2007, 25, 731–737. [Google Scholar] [CrossRef]
  37. Theodoratou, E.; Farrington, S.M.; Tenesa, A.; McNeill, G.; Cetnarskyj, R.; Barnetson, R.A.; Porteous, M.E.; Dunlop, M.G.; Campbell, H. Modification of the inverse association between dietary vitamin D intake and colorectal cancer risk by a FokI variant supports a chemoprotective action of Vitamin D intake mediated through VDR binding. Int. J. Cancer 2008, 123, 2170–2179. [Google Scholar] [CrossRef]
  38. Vidigal, V.M.; Silva, T.D.; de Oliveira, J.; Pimenta, C.A.M.; Felipe, A.V.; Forones, N.M. Genetic polymorphisms of vitamin D receptor (VDR), CYP27B1 and CYP24A1 genes and the risk of colorectal cancer. Int. J. Biol. Markers 2017, 32, e224–e230. [Google Scholar] [CrossRef]
  39. Serrano, D.; Gnagnarella, P.; Raimondi, S.; Gandini, S. Meta-analysis on vitamin D receptor and cancer risk: focus on the role of TaqI, ApaI, and Cdx2 polymorphisms. Eur. J. Cancer Prev. 2016, 25, 85–96. [Google Scholar] [CrossRef]
  40. Kadiyska, T.; Yakulov, T.; Kaneva, R.; Nedin, D.; Alexandrova, A.; Gegova, A.; Savov, A.; Mitev, V.; Kremensky, I. Vitamin D and estrogen receptor gene polymorphisms and the risk of colorectal cancer in Bulgaria. Int. J. Colorectal Dis. 2007, 22, 395–400. [Google Scholar] [CrossRef]
  41. Alkhayal, K.A.; Awadalia, Z.H.; Vaali-Mohammed, M.A.; Al Obeed, O.A.; Al Wesaimer, A.; Halwani, R.; Zubaidi, A.M.; Khan, Z.; Abdulla, M.H. Association of Vitamin D Receptor Gene Polymorphisms with Colorectal Cancer in a Saudi Arabian Population. PLoS ONE 2016, 11, e0155236. [Google Scholar] [CrossRef] [PubMed]
  42. Rosa, R.; Melisi, D.; Damiano, V.; Bianco, R.; Garofalo, S.; Gelardi, T.; Agrawal, S.; Di Nicolantonio, F.; Scarpa, A.; Bardelli, A.; et al. Toll-like receptor 9 agonist IMO cooperates with cetuximab in K-ras mutant colorectal and pancreatic cancers. Clin. Cancer Res. 2011, 17, 6531–6541. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Kinugasa, T.; Akagi, Y. Status of colitis-associated cancer in ulcerative colitis. World J. Gastrointest. Oncol. 2016, 8, 351–357. [Google Scholar] [CrossRef] [PubMed]
  44. Niedzielska, I.; Niedzielski, Z.; Tkacz, M.; Orawczyk, T.; Ziaja, K.; Starzewski, J.; Mazurek, U.; Markowski, J. Toll-like receptors and the tendency of normal mucous membrane to transform to polyp or colorectal cancer. J. Physiol. Pharmacol. 2009, 60, 65–71. [Google Scholar] [CrossRef] [PubMed]
  45. Krishnan, J.; Selvarajoo, K.; Tsuchiya, M.; Lee, G.; Choi, S. Toll-like receptor signal transduction. Exp. Mol. Med. 2007, 39, 421–438. [Google Scholar] [CrossRef]
  46. Yu, L.; Chen, S. Toll-like receptors expressed in tumor cells: targets for therapy. Cancer Immunol. Immunother. 2008, 57, 1271–1278. [Google Scholar] [CrossRef]
  47. Medvedev, A.E. Toll-like receptor polymorphisms, inflammatory and infectious diseases, allergies, and cancer. J. Interferon Cytokine Res. 2013, 33, 467–484. [Google Scholar] [CrossRef] [Green Version]
  48. Pradere, J.P.; Dapito, D.H.; Schwabe, R.F. The Yin and Yang of Toll-like receptors in cancer. Oncogene 2014, 33, 3485–3495. [Google Scholar] [CrossRef] [Green Version]
  49. Schauber, J.; Dorschner, R.A.; Coda, A.B.; Buchau, A.S.; Liu, P.T.; Kiken, D.; Helfrich, Y.R.; Kang, S.; Elalieh, H.Z.; Steinmeyer, A.; et al. Injury enhances TLR2 function and antimicrobial peptide expression through a vitamin D-dependent mechanism. J. Clin. Investig. 2007, 117, 803–811. [Google Scholar] [CrossRef] [Green Version]
  50. Harsch, M.; Bendrat, K.; Hofmeier, G.; Branscheid, D.; Niendorf, A. A new method for histological microdissection utilizing an ultrasonically oscillating needle: demonstrated by differential mRNA expression in human lung carcinoma tissue. Am. J. Pathol. 2001, 158, 1985–1990. [Google Scholar] [CrossRef]
  51. Wang, Q.; Xi, B.; Reilly, K.H.; Liu, M.; Fu, M. Quantitative assessment of the associations between four polymorphisms (FokI, ApaI, BsmI, TaqI) of vitamin D receptor gene and risk of diabetes mellitus. Mol. Biol. Rep. 2012, 39, 9405–9414. [Google Scholar] [CrossRef] [PubMed]
  52. Messaritakis, I.; Samonis, G.; Dimopoulou, D.; Maraki, S.; Papadakis, J.A.; Daraki, V.; Fragaki, M.; Choulaki, C.; Andrianaki, A.M.; Kofteridis, D.P. Staphylococcus aureus nasal carriage might be associated with vitamin D receptor polymorphisms in type 2 diabetes. Clin. Microbiol. Infect. 2014, 20, 920–925. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Souglakos, J.; Philips, J.; Wang, R.; Marwah, S.; Silver, M.; Tzardi, M.; Silver, J.; Ogino, S.; Hooshmand, S.; Kwak, E.; et al. Prognostic and predictive value of common mutations for treatment response and survival in patients with metastatic colorectal cancer. Br. J. Cancer 2009, 101, 465–472. [Google Scholar] [CrossRef] [PubMed]
Cancers 12 01379 g001 550
Figure 1. PCR-RFLP representative agarose gels of different SNPs of the (A) TaqI (TT, Tt, tt), (B) ApaI (AA, Aa, aa), (C) FokI (FF, Ff, ff) and (D) BsmI (BB, Bb, bb), respectively.
Figure 1. PCR-RFLP representative agarose gels of different SNPs of the (A) TaqI (TT, Tt, tt), (B) ApaI (AA, Aa, aa), (C) FokI (FF, Ff, ff) and (D) BsmI (BB, Bb, bb), respectively.
Cancers 12 01379 g001
Cancers 12 01379 g002 550
Figure 2. Overall survival according to the detection of (A) ApaI polymorphisms and (B) TaqI polymorphisms in stage II/III patients, and (C) TaqI polymorphisms, (D) ApaI polymorphisms, (E) FokI polymorphisms and (F) BsmI polymorphisms in stage IV patients.
Figure 2. Overall survival according to the detection of (A) ApaI polymorphisms and (B) TaqI polymorphisms in stage II/III patients, and (C) TaqI polymorphisms, (D) ApaI polymorphisms, (E) FokI polymorphisms and (F) BsmI polymorphisms in stage IV patients.
Cancers 12 01379 g002
Table
Table 1. Patients’ demographics and characteristics.
Table 1. Patients’ demographics and characteristics.
CharacteristicsFrequency (n = 397)%
Age (range)65 (18–88)
<7026065.5
≥7013734.5
Gender
Male24662
Female15138
Stage
IIA–IIIC20250.9
IV19549.1
Location
Colon/Sigmoid27970.3
Rectum11829.7
PS (ECOG)
0–137293.7
≥2256.3
Surgery
Yes34787.8
No5012.2
Radiotherapy
Yes8020.4
No31379.6
Adjuvant Treatment
Yes23058.7
No16241.4
First Line Treatment
Yes22356.2
No17446.8
Grade
High20547.4
Low22852.6
TaqI (T to C)
Wt (TT)12531.5
Hetero (Tt)12932.5
Homo (tt)14336
ApaI (T to G)
Wt (AA)14536.5
Hetero (Aa)12230.7
Homo (aa)13032.7
FokI (C to T)
Wt (FF)12230.7
Hetero (Ff)14737
Homo (ff)12832.2
BsmI (T to C)
Wt (BB)7518.9
Hetero (Bb)19148.1
Homo (bb)13133
KRAS
Mutant10442.4
Wild type14157.6
ND152
Table
Table 2. Association of VDR polymorphisms between patients and control groups (healthy donors and adenomatous polyps) and between different patients’ stages (II/III and IV).
Table 2. Association of VDR polymorphisms between patients and control groups (healthy donors and adenomatous polyps) and between different patients’ stages (II/III and IV).
SNPGenotypePatient No (%)Healthy No (%)Adenomatous Polyps No (%)p-ValueStage II/III No (%)Stage IV No (%)p-Value
TaqIWt (TT)125 (31.5)71 (71)24 (60.0)<0.001108 (51.9%)17 (9.0%)<0.001
Hetero (Tt)129 (32.5)26 (26.0)11 (27.5)75 (36.1%)54 (28.6%)
Homo (tt)143 (36.0)3 (3.0)5 (12.5)25 (12.0%)118 (62.4%)
ApaIWt (AA)145 (36.5)56 (56.0)22 (55.0)<0.001120 (57.7%)25 (13.2%)<0.001
Hetero (Aa)122 (30.7)40 (40.0)17 (42.5)74 (35.6%)48 (25.4%)
Homo (aa)130 (32.7)4 (4.0)1 (2.5)14 (6.7%)116 (61.4%)
FokIWt (FF)122 (30.7)55 (55.0)21 (52.5)<0.00198 (47.1%)24 (12.7%)<0.001
Hetero (Ff)147 (37.0)40 (40.0)16 (40.0)89 (42.8%)58 (30.7%)
Homo (ff)128 (32.2)5 (5.0)3 (7.5)21 (10.1%)107 (56.6%)
BsmIWt (BB)75 (18.9)55 (55.0)13 (32.5)<0.00155 (26.4%)20 (10.6%)<0.001
Hetero (Bb)191 (48.1)43 (43.0)24 (60.0)116 (55.8%)75 (39.7%)
Homo (bb)131 (33.0)2 (2.0)3 (7.5)37 (17.8%)94 (49.7%)
Table
Table 3. Correlation of various VDR and TLR polymorphisms. All values presented correspond to p-values.
Table 3. Correlation of various VDR and TLR polymorphisms. All values presented correspond to p-values.
All patients TaqIApaIFokIBsmITLR4 (Asp299Gly)TLR4 (Thr399Ile)TLR9 (T1237C)TLR9 (T1486C)TLR2 (196-to-174 del)
TaqI-0.0000.0000.0000.0000.0000.0000.0000.000
ApaI0.000-0.0000.0000.0000.0000.0000.0000.000
FokI0.0000.000-0.0000.0000.0000.0020.0020.000
BsmI0.0000.0000.000-0.0040.0040.0260.0260.009
Adjuvant TaqIApaIFokIBsmITLR4 (Asp299Gly)TLR4 (Thr399Ile)TLR9 (T1237C)TLR9 (T1486C)TLR2 (196-to-174 del)
TaqI-0.3010.0050.0000.8030.8030.3880.3880.346
ApaI0.301-0.5880.0140.7330.7330.9440.9440.148
FokI0.0050.588-0.0000.0350.0350.8080.8080.272
BsmI0.0000.0140.000-0.9660.9660.9050.9050.266
Metastatic TaqIApaIFokIBsmITLR4 (Asp299Gly)TLR4 (Thr399Ile)TLR9 (T1237C)TLR9 (T1486C)TLR2 (196-to-174 del)
TaqI-0.0000.2030.0000.8740.8740.0420.042-
ApaI0.000-0.0010.0000.2950.2950.0230.023-
FokI0.2030.001-0.0000.0370.0370.9730.973-
BsmI0.0000.0000.000-0.7030.7030.8140.814-
Table
Table 4. Association of VDR gene polymorphisms and KRAS status.
Table 4. Association of VDR gene polymorphisms and KRAS status.
KRAS
SNPGenotypeWt (%)Mutant (%)p-Value
TaqIWt (TT)37 (26.2)23 (22.1)0.041
Hetero (Tt)49 (34.8)25 (24.0)
Homo (tt)55 (39.0)56 (53.8)
ApaIWt (AA)40 (28.4)26 (25.0)0.007
Hetero (Aa)49 (34.8)20 (19.2)
Homo (aa)52 (36.9)58 (55.8)
FokIWt (FF)42 (29.8)16 (15.4)<0.001
Hetero (Ff)60 (42.6)26 (25.0)
Homo (ff)39 (27.7)62 (59.6)
BsmIWt (BB)24 (17.0)11 (10.6)<0.001
Hetero (Bb)78 (55.3)25 (24.0)
Homo (bb)39 (27.7)68 (65.4)
Table
Table 5. Univariate and multivariate Cox Regression analysis.
Table 5. Univariate and multivariate Cox Regression analysis.
Univariate AnalysisMultivariate Analysis
PFSOSPFSOS
HR (95% CI)p-ValueHR (95% CI)p-ValueHR (95% CI)p-ValueHR (95% CI)p-Value
Gender (Male vs. Female)1.0 (0.8–1.4)0.9860.9 (0.7–1.2)0.411----
Histology (Colon/Sigmoid vs. Rectum)1.2 (0.9–1.6)0.2671.2 (0.9–1.7)0.172----
PS ECOG (≥2 vs. 0–1)2.6 (1.6–4.2)<0.0014.4 (2.8–6.9)<0.0011.6 (1.0–2.6)0.043.7 (2.4–5.9)<0.001
Age (≥70 vs. <70)1.1 (0.8–1.4)0.6521.3 (1.0–1.8)0.06----
Grade (High vs. Low)1.4 (0.6–3.4)0.4934.4 (1.8–11.0)0.001--1.8 (0.7–5.1)<0.001
KRAS1.0 (0.8–1.4)0.7871.2 (0.8–1.6)0.362----
TLR2 196-to-174 del (homozygous mutant vs. others)2.1 (1.2–3.9)0.0163.7 (2.4–5.9)<0.0012.2 (1.2–4.1)0.0132.1 (1.3–3.4)0.003
TLR4-Asp299Gly (homozygous mutant vs. others)1.2 (0.9–1.7)0.2592.3 (1.7–3.3)<0.001--1.4 (1.0–1.9)0.095
TLR4-Thr399Ile (homozygous mutant vs. others)1.2 (0.9–1.7)0.2592.3 (1.7–3.3)<0.001--1.4 (1.0–1.9)0.095
TLR9-T1237C (homozygous mutant vs. others)1.0 (0.7–1.4)0.9281.7 (1.3–2.3)<0.001--1.1 (0.8–1.6)0.336
TLR9-T1486C (homozygous mutant vs. others)1.0 (0.7–1.4)0.9281.7 (1.3–2.3)<0.001--1.1 (0.8–1.6)0.336
TaqI (homozygous mutant vs. others)1.2 (0.9–1.2)0.7832.2 (1.8–2.7)<0.001--1.3 (1.0–1.6)0.029
ApaI (homozygous mutant vs. others)1.1 (0.9–1.3)0.462.6 (2.1–3.1)<0.001--1.6 (1.3–2.0)<0.001
FokI (homozygous mutant vs. others)1.1 (0.9–1.2)0.7682.2 (1.9–2.6)<0.001--1.4 (1.1–1.8)0.005
BsmI (homozygous mutant vs. others)0.1 (0.9–1.3)0.6161.7 1.4–2.2)<0.001--1.0 (0.7–1.6)0.644
Table
Table 6. PCR primers designed to amplify fragments harboring the VDR single nucleotide polymorphisms.
Table 6. PCR primers designed to amplify fragments harboring the VDR single nucleotide polymorphisms.
SNPPrimerSequenceFragment Size
TaqIForward5′-CAG AGC ATG GAC AGG GAG CAA-3′740 bp uncleaved
Reverse3′-GCA ACT CCT CAT GGC TGA GGT CTC-5′495 bp, 245 bp or 290 bp, 245 bp, 205 bp
ApaIForward5′-CAG AGC ATG GAC AGG GAG CAA-3′740 bp uncleaved
Reverse3′-GCA ACT CCT CAT GGC TGA GGT CTC-5′530 bp, 210 bp
FokIForward5′-AGC TGG CCC TGG CAC TGA CTC TGC TCT-3′265 bp uncleaved
Reverse3′-ATG GAA ACA CCT TGC TTC TTC TCC CTC-5′196 bp, 69 bp
BsmIForward5′-CAA CCA AGA CTA CAA GTA CCG CGT CAG TGA-3′825 bp uncleaved
Reverse3′-AAC CAG CGG GAA GAG GTC AAG GG-5′650 bp, 175 bp
Table
Table 7. Allele types, single nucleotide polymorphism (SNP) reference numbers and PCR conditions.
Table 7. Allele types, single nucleotide polymorphism (SNP) reference numbers and PCR conditions.
SNPAllele TypeRef. NumberPCR Conditions
TaqIt allele: silent T → C transition in exon 9rs731236Initial heating at 94 °C for 2 min, followed by 30 cycles of denaturing (at 94 °C for 30 s), annealing (at 57 °C for 30 s) and chain extension (at 72 °C for 1 min), followed by a final extension step at 72 °C for 7 min
ApaIa allele: T → G transition in intron 8rs7975232Initial heating at 94 °C for 2 min, followed by 30 cycles of denaturing (at 94 °C for 30 s), annealing (at 57 °C for 30 s) and chain extension (at 72 °C for 1 min), followed by a final extension step at 72 °C for 7 min
FokIf allele: C → T transition at the junction of intron 1 and exon 2rs2228570Initial heating at 95 °C for 2 min, followed by 35 cycles of denaturing (at 94 °C for 30 s), annealing (at 58 °C for 30 s) and chain extension (at 72 °C for 1 min), followed by a final extension step at 72 °C for 7 min
BsmIb allele: G → A transitionrs1544410Initial heating at 95 °C for 2 min, followed by 35 cycles of denaturing (at 94 °C for 30 s), annealing (at 60 °C for 30 s) and chain extension (at 72 °C for 1 min), followed by a final extension step at 72 °C for 7 min

Share and Cite

MDPI and ACS Style

Messaritakis, I.; Koulouridi, A.; Sfakianaki, M.; Vogiatzoglou, K.; Gouvas, N.; Athanasakis, E.; Tsiaoussis, J.; Xynos, E.; Mavroudis, D.; Tzardi, M.; et al. The Role of Vitamin D Receptor Gene Polymorphisms in Colorectal Cancer Risk. Cancers 2020, 12, 1379. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers12061379

AMA Style

Messaritakis I, Koulouridi A, Sfakianaki M, Vogiatzoglou K, Gouvas N, Athanasakis E, Tsiaoussis J, Xynos E, Mavroudis D, Tzardi M, et al. The Role of Vitamin D Receptor Gene Polymorphisms in Colorectal Cancer Risk. Cancers. 2020; 12(6):1379. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers12061379

Chicago/Turabian Style

Messaritakis, Ippokratis, Asimina Koulouridi, Maria Sfakianaki, Konstantinos Vogiatzoglou, Nikolaos Gouvas, Elias Athanasakis, John Tsiaoussis, Evangelos Xynos, Dimitriοs Mavroudis, Maria Tzardi, and et al. 2020. "The Role of Vitamin D Receptor Gene Polymorphisms in Colorectal Cancer Risk" Cancers 12, no. 6: 1379. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers12061379

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