
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Departments of 1 Environmental Health, 2 Biostatistics, 3 Nutrition, and 4 Epidemiology, Harvard School of Public Health; 5 Department of Medicine, Massachusetts General Hospital, 6 Department of Biostatistics and Computational Biology, Dana-Farber Cancer Institute, 7 Thoracic Surgery Unit, Department of Surgery, Massachusetts General Hospital, and 8 Channing Laboratory, Brigham and Women's Hospital, Harvard Medical School, Boston, Massachusetts
Requests for reprints: Wei Zhou, Occupational Health Program, Harvard School of Public Health, 665 Huntington Avenue, Boston, MA 02115. Phone: 617-432-1641. Fax: 617-432-6981. E-mail: wzhou{at}hsph.harvard.edu
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
The human VDR gene contains 11 exons and spans
75 kb, and several polymorphisms are "functional" by in vitro studies. The Cdx-2 G>A polymorphism (rs11568820) is located within the Cdx-2-binding site in the VDR gene promoter, with the A allele being associated with significantly higher VDR transcriptional activity than the G allele (5, 6). The FokI C>T polymorphism (rs10735810) is in the translational start site of VDR and has distinct structural consequences for the VDR, with the T allele (also called "f") being less efficient in exerting 1,25(OH)2D effects than the C allele (also called "F"; ref. 7). The BsmI C>T polymorphism (rs1544410) is located in intron 8 at the 3' end of VDR and is in high linkage disequilibrium with several other reported polymorphisms, including ApaI, TaqI, and poly-A (8, 9), with the T allele (also called "B") being associated with increased VDR mRNA expression and increased serum levels of 1,25(OH)2D compared with the C allele (also called "b"; refs. 10, 11). Polymorphisms of VDR have been associated with metastasis, recurrence, treatment response, or prognosis of breast cancer (12), malignant melanoma (13), prostate cancer (14, 15), and colorectal cancer (16); however, the results are not consistent in different studies, with little reproducibility from study to study.
Our previous data suggested that high vitamin D levels (the joint effects of summer surgery season and high recent vitamin D intake) at the time of treatment initiation may be associated with improved survival of early-stage nonsmall cell lung cancer (NSCLC) patients (17). Based on the "functional" data of VDR polymorphisms, we hypothesized that the variant A allele of the Cdx-2 polymorphism, the wild-type C allele of the FokI polymorphism, and the variant T allele of the BsmI polymorphism, which are thought to be associated with higher VDR transcriptional activities or higher vitamin D levels, are associated with improved survival of NSCLC and have joint effects with vitamin D status at the time of treatment initiation. Furthermore, we hypothesized that the joint effects of the three polymorphisms are stronger than the individual effect. Because histologic difference has been reported for the effects of 1,25(OH)2D and VDR on lung cancer (18), we investigated the above a priori hypotheses in our cohort of early-stage NSCLC patients, overall and in different histologic cell types of NSCLC.
| Materials and Methods |
|---|
|
|
|---|
Among the 456 eligible subjects, 83 patients were excluded because they did not have adequate blood or DNA samples for genotyping, leaving a total of 373 subjects. There were no statistically significant differences in patient characteristics between those with and those without genotype data. Characteristics examined include age, gender, histologic cell type, stage, radiotherapy, and chemotherapy after surgery. There were also no differences in overall survival (OS) between the groups, neither. All patients had surgical resection as the initial treatment, including wedge (24%), lobectomy (62%), bilobectomy (2%), pneumonectomy (6%), sleeve lobectomy (3%), and lobectomy plus wedge (1%). Additionally, 30 (8%) patients received postoperative radiation and/or adjuvant chemotherapy. The study was approved by the Human Subjects Committee of MGH and Harvard School of Public Health, and informed consent was collected from each participant.
VDR Genotyping
DNA was extracted from peripheral blood samples using the Puregene DNA Isolation kit (Gentra Systems, Minneapolis, MN). The three VDR polymorphisms were genotyped by the 5' nuclease assay (Taqman) using the ABI Prism 7900HT Sequence Detection System (Applied Biosystems, Foster City, CA). The primers, probes, and reaction conditions were available upon request. Genotyping was done by laboratory personnel blinded to case-control status, and a random 5% of the samples were repeated to validate genotyping procedures. Two authors reviewed independently all genotyping results.
Outcome Data Collection
OS and recurrence-free survival (RFS) were the end points in this study. OS was calculated from the date of surgery to the date of last follow-up or death from any cause. Dates of death were obtained and cross-checked using at least one of the following four methods: (a) inpatient and outpatient medical records, (b) MGH tumor registry, (c) Social Security Death Index, and (d) confirmation with the patient's primary care physician and/or family. Patients who were not deceased were censored at the last date they were known to be alive based on date of last contact. This date was verified by methods (b) and/or (d) as described above. Median follow-up time for this cohort was computed among alive subjects. RFS was defined as the time from the date of surgery to the first date of recurrence of cancer or death from any cause or the date of last follow-up. Date of recurrence was obtained by reviewing the hospital and outpatient records of all patients. For those 14% of patients who had their primary follow-up outside of the MGH system, we contacted the primary physician to obtain follow-up information.
Statistical Analysis
Demographic, clinical, and histologic information was compared across different genotype groups using Pearson
2 tests (for categorical variables) and Kruskal-Wallis tests (for continuous variables), where appropriate. Hardy-Weinberg equilibrium of each polymorphism was tested using the
2 test, and detection of linkage disequilibrium between the three polymorphisms was based on Lewontin's D'. The associations between individual VDR polymorphisms (or the total number of "protective" alleles from the three polymorphisms) and OS and RFS were estimated using the method of Kaplan and Meier and assessed using the log-rank test. Cox proportional hazards models were used as our primary analyses, controlling for multiple possible covariates simultaneously, including age, gender, stage, smoking status, and surgery season, where appropriate, with each covariate, including genotype groups or number of "protective" alleles represented by indicator variables. In addition to the overall analysis, we also investigated the effects of VDR polymorphisms on NSCLC in different histologic cell types [i.e., adenocarcinoma and squamous cell carcinoma (SCC), respectively].
Haplotype frequencies and individual haplotypes were generated using the expectation-maximization algorithm, which reconstruct individual probabilities for individual phasing accuracy based on unphased genotype data, as well as estimates on the overall haplotype frequencies and their SEs (19-22). The associations between VDR haplotypes and survival were estimated using the "expectation substitution" approach (19, 20, 22), which treats expected haplotype scores (calculated under additive model) as observed covariates in a standard Cox proportional hazards model, instead of assigning each subject with the most likely haplotype pair. All reported Ps are from two-sided tests. Ps < 0.05 were considered statistically significant. All analyses were done using Statistical Analysis System software version 9 (SAS Institute, Cary, NC).
| Results |
|---|
|
|
|---|
Detailed demographic, clinical, and treatment information by different histologic cell types is presented in Table 1 . Compared with adenocarcinoma, SCC patients had lower frequencies of females (37% versus 54%), stage IA (37% versus 58%), and never smokers (1% versus 9%). Surgery seasons were similar between different histologic cell types. Because of the small sample sizes of bronchioloalveolar carcinoma and large cell carcinoma, we limited the genotype survival analysis to adenocarcinoma and SCC in the subgroup analysis.
|
Individual VDR Polymorphisms and OS and RFS
The results of log-rank test and Cox proportional hazards models for VDR polymorphisms are presented in Table 2
(for OS). No association was found between the VDR Cdx-2 polymorphism and OS in the whole population or in adenocarcinoma patients. For the Cdx-2 polymorphism in SCC, due to the small number of patients with A/A genotype (n = 5) and similar survival effects with the G/A genotype (Table 2), we combined the G/A and A/A genotype groups in the analysis. The VDR Cdx-2 polymorphism was associated with statistically significant better survival among SCC patients. The 5-year OS survival rates [95% confidence interval (95% CI)] of the Cdx-2 polymorphisms were 41% (28-53) for the G/G genotype and 55% (39-71) for the G/A+A/A genotypes, respectively (P = 0.04, log-rank test; Fig. 1A
). In the Cox proportional hazards model, the adjusted hazard ratio (AHR) was 0.56 (95% CI, 0.33-0.95) for the G/A+A/A versus G/G. Similar improved survival effects of the Cdx-2 polymorphism were found for RFS in SCC: the 5-year RFS survival rates were 34% (95% CI, 22-46) for the G/G genotype and 50% (95% CI, 34-66) for the G/A+A/A genotypes, respectively (P = 0.03, log-rank test; Fig. 1B), with the AHR of 0.57 (95% CI, 0.34-0.94) for G/A+A/A versus G/G. Further analysis showed that the interaction between Cdx-2 polymorphism (G/A+A/A versus G/G) and histologic cell types (SCC versus adenocarcinoma) was statistically significant (P = 0.02 for both OS and RFS) in the Cox proportional hazards model. Similar results were observed in the Caucasian-only analysis.
|
|
|
Combined VDR Polymorphisms and OS and RFS
We also analyzed the joint effects of the three VDR polymorphisms based on the "function" of each polymorphism and total number of "protective" alleles. For the purposes of this analysis, the Cdx-2 G/G, FokI T/T, and BsmI C/C genotypes, which are associated with lower VDR transcriptional activity, lower VDR mRNA levels, or lower circulating vitamin D levels, were assigned as having zero "protective" allele; conversely, the Cdx-2 A/A, FokI C/C, and BsmI T/T genotypes were assigned as having two "protective" alleles. A total of four groups was generated: group 1, the reference group, with zero or one "protective" allele; group 2, with two "protective" alleles; group 3, with three "protective" alleles; and group 4, with four or more "protective" alleles. Results of log-rank test and Cox proportional hazards models are shown in Table 4
. There was no statistically significant association between the joint polymorphisms and OS or RFS, overall or in adenocarcinoma patients. Among SCC patients, those with two or more "protective" alleles have statistically better survival when compared with the reference group in Cox proportional hazards models; however, no clear "dose-response" effect was observed. For OS, the AHRs were 0.20 (95% CI, 0.09-0.48), 0.40 (95% CI, 0.19-0.87), and 0.43 (95% CI, 0.19-0.97) for subjects with two, three, and four or more "protective" alleles (Ptrend = 0.71); for RFS, the corresponding AHRs were 0.23 (95% CI, 0.10-0.51), 0.43 (95% CI, 0.20-0.89), and 0.44 (95% CI, 0.20-0.97), respectively (Ptrend = 0.70). Similar results were observed in the Caucasian-only analysis.
|
|
| Discussion |
|---|
|
|
|---|
The associations between VDR and cancer metastases and survival are supported by animal models. Using cells isolated from VDR knockout and wild-type mice, both normal and transformed mammary cells derived from wild-type mice are growth inhibited by 1,25(OH)2D. However, cells derived from VDR knockout mice are completely unresponsive to 1,25(OH)2D (23). An inverse relationship between VDR levels and both colonic hyperproliferation and oxidative stress has been found in a mouse model (24). In humans, significant higher expression levels of VDR have been found in lung cancer tissues than in normal tissues (25, 26), and VDR activity may be enhanced by ß-catenin (27), which makes cells in the tumor mass more adherent to each other and reduces the likelihood of mobilization of large numbers of malignant cells into the lymphatic or blood circulation.
VDR is a highly polymorphic gene. Although functional significances between different genotypes of VDR polymorphisms have been reported based on in vitro studies, the associations between VDR polymorphisms and cancer survival are not consistent across studies, and there is no report on the association between VDR genotypes or haplotypes and lung cancer survival (28). The BsmI C>T polymorphism is the most widely studied VDR polymorphism in cancer prognosis. The T/T genotype may protect the development of metastases among breast cancer patient (29, 30) and may affect the prognosis of rectal cancer by influencing erbB-2 oncogene expression, although it alone did not influence survival (16). The BsmI T allele has also been associated with reduced acute graft-versus-host disease when present in the patient's genotype (31) and protected against recurrence of locally advanced prostate cancer (15). In a small study with 191 mostly Caucasian prostate cancer patients, the C allele of the FokI C>T polymorphism was associated with an increased risk of aggressive prostate cancer and could therefore be associated with worse prognosis (14). We did not observe the association between the individual BsmI or FokI polymorphism and NSCLC survival, whereas the joint effects of the three polymorphisms and haplotype analysis suggested better survival is associated with the "protective" alleles of the three polymorphisms, although no "dose-response" relationship was observed for the number of "protective" alleles (Table 4). Therefore, all of the three investigated VDR polymorphisms may contribute to better survival of NSCLC patients, with the Cdx-2 polymorphism having the strongest effect.
In this study, the association with VDR polymorphisms was observed among SCC patients whereas not among adenocarcinoma patients. In our previous analysis, we also observed a slightly stronger effect of surgery season and/or vitamin D intake on NSCLC survival among SCC patients than adenocarcinoma patients, although the difference was not statistically significant (data not shown). In vitro studies have suggested that vitamin D inhibits growth of a lung SCC cell line whereas not adenocarcinoma cell line, and the mRNA levels of VDR were much higher in the SCC cell line than the adenocarcinoma cell line (18). The mechanism for this histologic difference of VDR polymorphisms may be associated with the effect of retinoic acid, retinoic acid receptors (RAR), or retinoid X receptors (RXR). VDR is a member of the steroid/retinoid receptor superfamily of nuclear receptors. In the cell, 1,25(OH)2D binds to the VDR, and the VDR*1,25(OH)2D complex then interacts with the RXR to form a 1,25(OH)2D*VDR*RXR heterodimer complex, which then interacts with vitamin Dresponsive elements. RARs also function as heterodimers with RXR proteins (2). Because both retinol and vitamin D require RXR proteins for their actions, high doses of retinol may antagonize vitamin D actions (32, 33). In lung cancer, the expression levels of RARß gene have been found to be highly expressed in adenocarcinoma cell lines whereas not detectable in SCC cell lines (34). Retinoic acid has also shown its growth-inhibitory activity when applied to the adenocarcinoma cell line of human lung and stomach whereas is inactive on the esophageal SCC cell lines (35). Therefore, in adenocarcinoma patients, the effect of VDR may be overwhelmed by the highly expressed RAR levels, which may be one reason that we did not observe a significant effect of VDR polymorphism on survival among adenocarcinoma patients.
There are several limitations for this study. (a) Sample size. Although this is a case follow-up study with moderate sample size, the numbers are small in the SCC subgroup and for the gene-environment and gene-gene joint analyses. (b) Selection of candidate polymorphisms. We investigated three "functional" VDR polymorphisms based on in vitro data, and it is possible that these polymorphisms may not be truly functional. We did not investigate the other reported VDR polymorphisms, including ApaI, TaqI, and poly-A, which may introduce bias in the results, especially in the haplotype analysis. However, adding more polymorphisms in the results may also introduce multiple comparisons in the genotype analysis. Given the high linkage disequilibrium of these polymorphisms with the BsmI polymorphism, their association with NSCLC survival will be similar to the results of the BsmI polymorphism. (c) Missing data. In our population, 83 patients were excluded because they did not have adequate blood or DNA samples for genotyping. However, the distributions of demographic, histologic, treatment, and survival characteristics were very similar between patients with and without genotype information. (d) Survival data collection. In this population, recurrence data were collected retrospectively and patients were not on a prescribed surveillance schedule. We attempted to contact local physicians whenever patients were followed outside of the MGH system (14%). We collected the vital status data for each patient. However, we could not distinguish between the death from lung cancer and the death from other causes. Because the 5-year OS rates in this population were 65%, 53%, 34%, and 38%, respectively, for stages IA to IIB, the vast majority of these patients likely died from lung cancer. (e) We did not have complete information on either performance status or weight loss in our early-stage patients. However, in a random sample of 100 patients who had performance status data available, surgical resectability is a reasonable surrogate measure for good performance status. We also had very incomplete weight loss information, with most of the available information dichotomized into weight loss present or absent. Our preliminary data in a random sample of 40 stage IA patients and 10 stage II patients suggested that disease stage (in early-stage patients) may be crudely correlated with weight loss. In our models, we do adjust for disease stage. These factors are unlikely to be associated with VDR polymorphisms and will not introduce systematic bias in the results. (f) We lack information on the serum levels of vitamin D or other genes that are involved in vitamin D metabolism, which may limit the generalizability of our results.
In conclusion, this is the first study suggesting that VDR polymorphisms may be associated with improved survival in early-stage NSCLC patients. Our results suggest that the A allele (G/A+A/A genotypes) of the Cdx-2 polymorphism, the combined "protective" genotypes of the Cdx-2, FokI, and BsmI polymorphisms, and the A-C-T haplotype (Cdx-2-FokI-BsmI) are associated with better OS and RFS among SCC patients whereas not among adenocarcinoma patients. The results are consistent with our previous findings on surgery season and vitamin D intake, further suggesting that vitamin D may be associated with improved survival in early-stage NSCLC patients. These results need to be confirmed by independent prospective studies or clinical trials.
| Acknowledgments |
|---|
| Footnotes |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Received 1/12/06; revised 8/ 9/06; accepted 8/18/06.
| References |
|---|
|
|
|---|
,25-Dihydroxyvitamin D(3) is a preventive factor in the metastasis of lung cancer. Carcinogenesis 2005;26:42940.
heterodimers are specified by dinucleotide differences in the vitamin D-responsive elements of the osteocalcin and osteopontin genes. Mol Endocrinol 1996;10:144456.
,25-Dihydroxyvitamin D3 and all-trans-retinoic acid inhibit the growth of a lung cancer cell line. Anticancer Res 1996;16:26539.[Medline]This article has been cited by other articles:
![]() |
W. Janssens, A. Lehouck, C. Carremans, R. Bouillon, C. Mathieu, and M. Decramer Vitamin D Beyond Bones in Chronic Obstructive Pulmonary Disease: Time to Act Am. J. Respir. Crit. Care Med., April 15, 2009; 179(8): 630 - 636. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. S. Heist, W. Zhou, Z. Wang, G. Liu, D. Neuberg, L. Su, K. Asomaning, B. W. Hollis, T. J. Lynch, J. C. Wain, et al. Circulating 25-Hydroxyvitamin D, VDR Polymorphisms, and Survival in Advanced Non-Small-Cell Lung Cancer J. Clin. Oncol., December 1, 2008; 26(34): 5596 - 5602. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. M. Ochs-Balcom, M. S. Cicek, C. L. Thompson, T. C. Tucker, R. C. Elston, S. J.Plummer, G. Casey, and L. Li Association of vitamin D receptor gene variants, adiposity and colon cancer Carcinogenesis, September 1, 2008; 29(9): 1788 - 1793. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Zhou, R. S. Heist, and D. C. Christiani In Reply J. Clin. Oncol., December 1, 2007; 25(34): 5538 - 5539. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |