Predicted 25(OH)D Score and Colorectal Cancer Risk According to Vitamin D Receptor Expression

Vitamin D has long been hypothesized to be associated with a lower risk of colorectal cancer (1–3). A 10 ng/mL increment in blood 25-hydroxyvitamin D [25(OH)D] level was associated with 26% lower risk of colorectal cancer [95% confidence interval (CI), 0.63–0.89] in a meta-analysis of nine prospective studies (4). However, colorectal cancer is a heterogeneous disease in which each tumor evolves through distinct carcinogenic pathways acquiring a unique set of genetic and epigenetic aberrations over time (5). The inhibitory property of vitamin D on colorectal carcinogenesis may differ according to specific tumoral features of the colon (6).

Numerous studies support that vitamin D receptor may mediate the anticarcinogenic effect of vitamin D (1, 7). The 1,25-dihydroxycholecalciferol [1,25(OH)₂D], a hydroxylated form of 25(OH)D, binds to the vitamin D receptor (VDR). Then, activated VDR heterodimerizes with the retinoid X receptor alpha (RXRA; refs. 1, 7). This complex translocates to the cell nucleus and binds to the vitamin D response element to regulate the transcription of a number of genes that control cellular proliferation, differentiation, angiogenesis, inflammation, and apoptosis (1). No previous studies have prospectively investigated the role of vitamin D according to VDR expression level in colorectal tumors. Thus, we hypothesized that a higher vitamin D status was associated with a lower risk of colorectal cancer with high VDR expression level but not with a risk of cancer with low VDR expression level.

To test our hypothesis, we prospectively examined the association of a long-term predicted 25(OH)D score (8) with the risk of colorectal cancer according to the level of VDR expression in the Nurses' Health Study (NHS) and Health Professionals' Follow-up Study (HPFS). The predicted 25(OH)D score comprehensively takes into account both endogenous and exogenous sources of plasma 25(OH)D; cumulatively averaged 25(OH)D scores, in particular, could well represent long-term levels of plasma 25(OH)D status (8). In the secondary analyses, we evaluated the association between the predicted 25(OH)D score and the risk of colorectal cancer as defined by the mutation status of KRAS, BRAF, and PIK3CA, which previously was found to interact with VDR (9–13).

The NHS enrolled 121,700 U.S. female registered nurses ages 30 to 55 from 11 U.S. states in 1976. The HPFS enrolled 51,529 U.S. male health professionals ages 40 to 75 from all 50 states in 1986. The details for both cohorts have been described previously (14). In brief, NHS and HPFS participants completed questionnaires inquiring about lifestyle factors and chronic disease history at the initiation of each cohort study. This information has been updated via biennial questionnaires thereafter.

In this analysis, the baseline year was 1986, when all the determinants of plasma 25(OH)D were first collected to compute the predicted 25(OH)D score. We excluded participants with a missing predicted 25(OH)D score and those previously diagnosed with cancer (except nonmelanoma skin cancer) or ulcerative colitis. After exclusion, we included 96,239 women from the NHS and 44,197 men from the HPFS. The Institutional review board at the Brigham and Women's Hospital (Boston, MA), and the Harvard School of Public Health (Boston, MA) approved the NHS and the HPFS studies, respectively. All participants provided informed consent at enrollment.

The derivation and validation of the predicted 25(OH)D score in NHS and HPFS have been previously described (8). In brief, significant clinical predictors of circulating plasma 25(OH)D, including race, UV-B flux, dietary vitamin D, supplementary vitamin D, body mass index, physical activity, alcohol intake (NHS only), and postmenopausal hormone use (NHS only), were regressed on a plasma 25(OH)D concentration among the population of NHS (N = 2,079) and HPFS (N = 911) who were disease-free and had their blood measured for plasma 25(OH)D. The predicted 25(OH)D score was computed by summing the value of each significant determinant that was weighted by their individual associations with plasma 25(OH)D.

The validity of the predicted 25(OH)D score was assessed by using the plasma 25(OH)D level, a gold standard indicator of vitamin D status (8) because it reflects vitamin D from both dietary and solar sources (15) with a three-week half-life (16). In the independent subpopulation of NHS and HPFS whose plasma 25(OH)D was not used for the creation of a predicted 25(OH)D score, the correlation coefficients between predicted 25(OH)D score and plasma 25(OH)D, adjusted for laboratory batch, age, and season of blood draw, were 0.33 for NHS and 0.30 for HPFS, respectively. The plasma 25(OH)D concentration increased with the 25(OH)D score; the mean differences of the actual 25(OH)D level between extreme deciles of the predicted 25(OH)D score were 8.7 ng/mg for both NHS and HPFS. In addition, the predicted 25(OH)D score per 10 ng/mL yielded an OR of 0.78 for colorectal cancer risk, which is similar to the risk estimate associated with measured plasma 25(OH)D (OR, 0.82 per 10 ng/mL increment; refs. 8, 17).

The participants reported newly diagnosed colorectal cancer on each follow-up questionnaire. Among non-respondents, we searched the National Death Index to identify the cause of fatality (18). Through these methods we identified approximately 96% to 97% of cases (19), and we requested permission to review the medical records and pathology reports. Investigators, blinded to the information of the participants, confirmed the reported cases and extracted information on histologic type, stage, and anatomic location of cancer. Among 2,604 colorectal cancer cases identified during the study follow-up, we collected 1,059 tumor specimens for this study. Among those, we obtained data on tissue microarray immunohistochemistry of VDR expression (N = 743) and mutation status of KRAS (N = 1,045), BRAF (N = 1,044), and PIK3CA (N = 970).

For this study, we measured the nuclear VDR expression level in colorectal cancer tissue because VDR modulates genomic transcription in cell nuclei (1, 7) and the nuclear VDR level, not the cytoplasmic VDR level, has been a prognostic marker of improved survival (20, 21). We retrieved formalin-fixed, paraffin-embedded colorectal cancer tissue blocks from hospitals throughout the United States where patients with colorectal cancer had undergone surgical resection. Tissue microarray blocks were constructed as previously described (22). Colorectal cancer cases with VDR expression are limited to those with available tissue microarray for the immunohistochemistry. The pathologists were blinded to patients' information on lifestyle or medical histories. The demographic characteristics of cases with available nuclear VDR expression were similar to those without nuclear VDR expression data [mean age, 56.9 vs. 57.5 years; women, 64 vs. 56%; White, 93% vs. 93%; mean BMI (kg/m²), 25.9 vs. 25.9; mean physical activity (metabolic equivalent score per week), 14.9 vs.16.3; mean pack years of smoking, 4.8 vs. 4.6; current multivitamin use, 39% vs. 38%; family history of colorectal cancer, 11% vs. 13%; mean red meat intake, 1.2 vs. 1.2 servings/day, mean calcium intake 884.2 vs. 864.8 g/day, mean folate intake 419.0 vs. 421.0 μg/day; P ≥ 0.05 for all comparisons].

For VDR immunostaining, deparaffinized tissue sections were heated using a pressure cooker in a microwave for 15 minutes in Antigen Retrieval Citra Solution, pH 6 (BioGenex Laboratories). Tissue sections were incubated with Dual Endogenous Enzyme Block (DAKO), then Serum Free Protein Block (DAKO), each for 15 minutes. Slides were incubated at 4°C overnight with a primary antibody against VDR (1:500, rabbit polyclonal; Novus Biologicals, NBP1-19478), diluted in Da Vinci Green Diluent (Biocare Medical). Envision anti-rabbit HRP-labeled polymer (DAKO) was applied to the sections for 30 minutes followed by visualization using the chromogen 3,3-diaminobenzidine (DAKO), and hematoxylin counterstain; VDR (N204) detects endogenous levels of VDR protein and synthetic peptide corresponding to the residues surrounding Asparagine 204 of human VDR.

The level of nuclear VDR expression in the tumor tissue was assessed using a semiquantitative immunoreactivity scoring (SIS) system (23). The staining intensity was scored as 1 (no immunostaining) to 4 (strong). The percentage of immunoreactive nuclear cells was rated from 0% to 100%. The SIS score was calculated by multiplying the scores for expression intensities with the percentage of positive cells resulting in the score variation from 0 to 400. We defined VDR expression as positive (SIS score ≥ 180) and negative (SIS score < 180); the cutoff for dichotomization of nuclear VDR expression level was chosen at the nuclear VDR level where survival length among colorectal cancer survivors significantly differs. In addition, we classified colorectal tumor into three subgroups using the tertile of VDR expression level in secondary analyses.

DNA was extracted from paraffin-embedded tissues. PCR and pyrosequencing were performed for KRAS (codons 12, 13, 61, and 146; refs. 24, 25), BRAF (codon 600; ref. 26), and PIK3CA (exons 9 and 20; refs. 27), as described previously.

Lifestyle and other information (body mass index, physical activity, smoking status, aspirin use, multivitamin use, family history of colorectal cancer in first-degree relatives, and history of endoscopy) was self-reported on biennial questionnaires. Dietary information was assessed via validated semiquantitative food frequency questionnaires with approximately 130 food items every 4 years in the NHS (28) and the HPFS (29). Daily nutrient intake was calculated by multiplying the frequency response of each specified food item by the nutrient content of the specified portion sizes then summing these products for all food items.

We categorized predicted 25(OH)D into quintiles separately within each cohort. To investigate long-term 25(OH)D exposure and minimize the influence of measurement error, we used the cumulatively averaged predicted 25(OH)D score as our exposure. Sensitivity analyses using a simple updated and baseline predicted 25(OH)D were conducted.

We calculated HRs and 95% confidence intervals (95% CI) separately for incidence of overall and subtypes of colorectal cancer using duplication method Cox proportional cause-specific hazards regression for competing risk data (30). Cases without information on relevant molecular subtypes were censored at the time of colorectal cancer diagnosis. Person-years of follow-up were calculated from the date of baseline questionnaire return to the date of diagnosis of colorectal cancer, date of death, loss to follow-up, or the end of follow-up (2008 for NHS; 2008 for HPFS), whichever came first. We included age (in months) and the year of questionnaire cycle as stratification variables. In multivariate analyses, we included potential confounding variables. We updated covariates in each questionnaire cycle to take into account potential changes over time. A missing indicator for missing responses of each covariate was created, if applicable. We tested for trend using the Wald test of the continuous variable set to the median values of quintiles of the 25(OH)D score. We pooled participants from the NHS and the HPFS and included a cohort indicator as stratification variable in the model. Before pooling the cohorts, we tested between-studies heterogeneities using Q statistic (31, 32) and found no significant heterogeneities (P_{heterogeneity} between studies > 0.19 for all analyses). To test the significance of differential association by the tumor molecular characteristic, we conducted the likelihood ratio test comparing the model fit that allows separate associations by different molecular subtypes (i.e., VDR-positive vs. VDR-negative cancers) to the model fit that assumed a common effect. With 80% power, the ratio of relative risk (RR) for detecting significant heterogeneity by tumor subtypes were 1.30 for the analysis of VDR expression, 1.15 for KRAS mutation, 1.65 for BRAF mutation, and 1.59 for PIK3CA mutation.

In addition, to evaluate the nonlinearity of the association of the predicted 25(OH)D score with the risk of overall and molecular subtypes of colorectal cancer, we compared the model fit between the model with the linear term and cubic spline terms and the model without spline terms (33–35); the log-likelihood ratio test comparing those two models was statistically significant suggesting a nonlinear association. Therefore, we considered the predicted 25(OH)D score as a categorical variable in our analyses.

We used the SAS software (SAS Institute, Inc., Version 9). All tests were two-sided and a P < 0.05 was considered statistically significant.

Among 140,418 participants followed up from 1986 through 2008, we documented 1,059 incident colorectal cancer cases (41% of all colorectal cancer cases) with tumor molecular data on the level of VDR expression (N = 743) and mutation status for KRAS (N = 1,045), BRAF (N = 1,044), and PIK3CA (N = 970). Table 1 summarizes the age-standardized characteristics of the study population at the median follow-up time according to the level of the predicted 25(OH)D score. As expected, participants with higher predicted 25(OH)D scores were more likely to live in the Southern states, have a lower body mass index, report higher physical activity, have higher multivitamin intake, and consume more vitamin D, calcium, and folate. The differences in the median across the extreme quintiles of the predicted 25(OH)D score were 7.8 ng/mL in men and 10.4 ng/mL in women, which is equivalent to 22–29 μg/day increase in vitamin D intake (ref. 36; Table 1).

Because the results from men and women did not significantly differ (Supplementary Tables S1 and S2), here, we present the results from the pooled analyses. We first examined the association of the predicted 25(OH)D score with the risk of overall colorectal cancer with available molecular data in our cohorts (Table 2). The predicted 25(OH)D score was inversely associated with the risk of colorectal cancer. In the pooled cohort, the multivariate HR (95% CI) was 0.52 (0.42–0.64) comparing the highest with the lowest quintile of the predicted 25(OH)D score (P_trend < 0.001). Because some individual components of the predicted 25(OH)D score are associated with colorectal cancer, we evaluated whether the observed association between the predicted 25(OH)D score and colorectal cancer incidence was confounded by those individual components (8). We found that the inclusion of body mass index, physical activity, and alcohol consumption into the model yielded similar results to our original main model (Table 2).

We further evaluated the association of the predicted 25(OH)D score with the risk of colorectal cancer according to VDR expression level (Table 2). We found that the predicted 25(OH)D score had a statistically significant inverse association with the risk of both VDR-positive and VDR-negative colorectal tumors. In the pooled cohort, multivariate HRs (95% CI) were 0.48 (0.30–0.78) for VDR-negative and 0.56 (0.42–0.75) for VDR-positive colorectal cancer comparing the highest with the lowest quintile of the predicted 25(OH)D score. The difference in the associations for VDR-positive versus VDR-negative tumors was not statistically significant (P_{heterogeneity} for subtypes = 0.75). In addition, we grouped the colorectal tumors using the tertile of VDR expression level to evaluate whether the inverse association of the predicted 25(OH)D score with colorectal cancer risk differed by VDR expression level with a linear trend; we found no evidence of linear trend for heterogeneity. Pooled multivariate HRs (95% CI) of predicted 25(OH)D comparing extreme quintiles were 0.47 (0.29–0.74), 0.63 (0.43–0.93), and 0.51 (0.33–0.78) for colorectal tumors with low, medium, and high VDR levels, respectively (P_{heterogeneity} for subtypes = 0.67).

We also examined whether the association of the predicted 25(OH)D score with the risk of colorectal cancer differed according to the mutation status of KRAS, BRAF, or PIK3CA (Table 3). A higher predicted 25(OH)D score was inversely associated with colorectal cancer risk, regardless of KRAS, BRAF, or PIK3CA status (P_{heterogeneity} for subtypes ≥ 0.22). Pooled multivariate HRs (95% CI) were 0.46 (0.35–0.60) for KRAS wild-type versus 0.70 (0.50–0.98) for KRAS-mutant, 0.52 (0.41–0.65) for BRAF wild-type versus 0.58 (0.35–0.94) for BRAF-mutant, and 0.51 (0.40–0.65) for PIK3CA wild-type versus 0.54 (0.32–0.90) for PIK3CA-mutant (P_trend ≤ 0.01). Given the interrelationship of key regulatory molecules that may obscure an association, we further evaluated the association of the predicted 25(OH)D score with colorectal cancer risk according to the combination of VDR with KRAS or PIK3CA mutation status (9, 10, 12, 13). However, we did not observe substantial or statistically significant differences in risk estimates according to any of these categories of tumor subtypes (data not shown).

We conducted sensitivity analyses excluding the participants who were diagnosed within the first 4 years of follow-up to minimize the influence of alterations in lifestyle due to prediagnostic symptoms or occult cancer. These results were not materially different from our original findings (data not shown).

In this large prospective study, the inverse association between the predicted 25(OH)D score and colorectal cancer incidence did not significantly differ according to VDR expression level. In addition, we did not observe significant differences in the association by mutation status of KRAS, BRAF, and PIK3CA. To the best of our knowledge, our study is the first to examine the association of vitamin D with the risk of colorectal cancer according to these tumor molecular features. Our results suggest that the apparent benefit of vitamin D on the risk of colorectal cancer may be uniform regardless of the heterogeneous carcinogenic pathways analyzed in this study.

In the current literature, several lines of mechanistic evidence have suggested that vitamin D may have a stronger inverse association with the risk of colorectal tumor with high level of VDR expression. Colon epithelial and cancer cells express both VDR (37) and 1-α-hydroxylase (CYP27B1; ref. 38), which converts 25(OH)D into the 1,25(OH)₂D that binds to VDR. The recent genome-wide association study identified numerous VDR-binding sites in the colon (39), thus suggesting that the colon may be more likely to be influenced by VDR signaling than are other organs. Furthermore, experimental studies observed that the growth arrest induced by vitamin D disappeared in VDR knockout cell lines (40, 41).

However, we found no evidence of differential association of the predicted 25(OH)D score with the risk of colorectal cancer according to VDR expression level. Furthermore, we also did not observe significantly varied association of the predicted 25(OH)D score with colorectal tumor defined by mutation status of KRAS, BRAF, and PIK3CA that are downstream of EGFR signaling (42). Previously, some studies suggested the anticancer effect of vitamin D in the oncogenic pathways involving VDR and RAS mitogen-activated protein kinase (MAPK; refs. 9, 12, 13) or PI3K-AKT pathways (9). Despite these experimental data suggesting the potential variation of association between vitamin D and colorectal cancer according to tumor features, there was no empirical evidence from prospective population-based studies.

The consistent and significant inverse association of the predicted 25(OH)D score with the risk of all subtypes of colorectal cancer we observed suggests that vitamin D may not act through a single mechanism to inhibit carcinogenesis. Supporting this speculation, experimental studies have demonstrated that vitamin D reduces proliferation, inflammation, and angiogenesis, stimulates differentiation and apoptosis, and enhances the immune system (1, 2). Recent observations that suggest the interaction of vitamin D with IGF1 (43) and the influence of vitamin D on the regulation of microRNAs (44) and chromatin epigenetic alteration (44, 45) add to the mechanistic complexity of the action of vitamin D. Confirmation from a future large study is warranted given that our study is the first epidemiologic study that used VDR expression level in colorectal tumor tissue. Nonetheless, the consistent inverse associations we observed may reflect the pleiotropic antitumoral biologic function of vitamin D.

Our result of the inverse association between the predicted 25(OH)D score and the risk of colorectal cancer is aligned with current literature. Previous results have consistently supported the benefit of vitamin D in preventing colorectal cancer (46). The largest nested case–control study that measured 25(OH)D level among 1,248 cases and 1,248 controls across 10 Western European countries reported the RR for colorectal cancer of 0.60 (95% CI, 0.46–0.80) comparing the highest to the lowest quintile of plasma 25(OH)D (47). The meta-analyses of 9 studies (7 cohort studies and 2 nested case–control studies) with 2,767 cases and 3,948 controls reported that the summary RR was 0.67 (95% CI, 0.54–0.80) comparing the highest versus lowest categories of blood 25(OH)D levels (4).

Our study has several strengths. Prospectively collected data eliminated recall bias. Bias that could have occurred in the selection of controls in a case–control design was avoided in our prospective study design. The follow-up rate of the participant in NHS and HPFS is high over 90%. With a long-term follow-up of 24 years, we were able to examine several potential latency periods relating vitamin D to cancer risk. We were able to evaluate several tumor markers concurrently through comprehensive assessment of tumor pathologic and molecular characteristics. Our integrative molecular pathologic epidemiology approach aimed to examine potential links between exposures and molecular signatures of disease, and obtain mechanistic insights on biologic influences of exposures on molecular pathways to the disease (48–50). Tumor molecular analyses have become increasingly common in research and clinical practice (51–54). By cumulatively updating the predicted 25(OH)D score, we could examine average long-term vitamin D status compared to using a single measurement of actual 25(OH)D (half-life 3 weeks; ref. 16). The predicted 25(OH)D score in our analyses accounts for both diet and sun exposure, capturing the variation of circulating plasma 25(OH)D comprehensively (8). In addition, detailed information on lifestyle, diet, and other risk factors allowed us to finely adjust for potential confounding factors for colorectal cancer.

Our study has several limitations. Regarding the use of the predicted 25(OH)D score, the predicted 25(OH)D level may be confounded by its predictors; however, our results did not change materially when we added individual components of the score into the main model. Although the predicted 25(OH)D score is a validated predictor of the plasma 25(OH)D level in the NHS and HPFS, the misclassification of the actual level of 25(OH)D is still possible due to factors that were not considered at the derivation of the predicted score. For example, a recently large pooled GWAS study identified polymorphisms that significantly predict circulating 25(OH)D levels (55). The measurement error might be nondifferential attenuating the associations observed.

In addition, we could not retrieve the tumor marker information for all of our colorectal cancer cases. However, the risk estimates for incident colorectal cancer using cancer cases with tumor marker information were similar to those for incident cancer in the entire population (56). Because of limited cases, our study had limited power to detect differential associations of colorectal cancer risk by molecular subtypes (5, 48). Another limitation is that our measurement of VDR expression level may not reflect changes in VDR levels as a tumor progresses (37). There is currently no standardized method to assess VDR expression level in colorectal tumors or a consensus on cutoff to classify VDR overexpression. However, results did not change when we alternatively divided colorectal cancer into two groups using the median value of tumoral VDR expression level (data not shown), and we did not observe any trend for the association of the predicted 25(OH)D score with the risk of colorectal cancer when we further classified the colorectal tumor by tertile categories of tumoral VDR expression level. In addition, we cannot rule out whether our result is driven by residual confounding. For example, VDR signaling needs not only vitamin D as a ligand but also the retinoid X receptor (57) and protein complexes (58) to be activated and to unwind the chromosomal constraint (57). However, the distribution of those multiple components may not be differential with respect to the level of either vitamin D or VDR expression in a colorectal tumor.

In conclusion, we observed that prediction of a high predicted 25(OH)D score was significantly associated with the lower risk of colorectal cancer regardless of VDR expression level or molecular features of colorectal tumors. Our observation supports a broad influence of vitamin D on colorectal carcinogenesis that may involve multiple biologic pathways including host immunity. As vitamin D status is determined by many modifiable lifestyle factors, our results may further indicate that diet and lifestyle associated with high level of vitamin D level can be recommended to prevent colorectal cancer. Future studies with large numbers of cases are warranted to replicate our observations.

No potential conflicts of interest were disclosed.

The content is solely the responsibility of the authors and does not necessarily represent the official views of NIH or ASCO. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the article.

Conception and design: Z.R. Qian, C.S. Fuchs, A.T. Chan, K. Ng, E. Cho, S. Ogino

Development of methodology: Z.R. Qian, C.S. Fuchs, E.L. Giovannucci, K. Ng, E. Cho, S. Ogino

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): Z.R. Qian, M. Yamauchi, K. Inamura, S.A. Kim, K. Mima, C.S. Fuchs, A.T. Chan, S. Ogino

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): S. Jung, K.A. Bertrand, K.C. Fitzgerald, Y. Sukawa, X. Zhang, M. Wang, S.A. Smith-Warner, K. Wu, C.S. Fuchs, A.T. Chan, E.L. Giovannucci, K. Ng, S. Ogino, R. Nishihara

Writing, review, and/or revision of the manuscript: S. Jung, Z.R. Qian, K.A. Bertrand, X. Zhang, S.A. Smith-Warner, K. Wu, C.S. Fuchs, A.T. Chan, E.L. Giovannucci, K. Ng, E. Cho, S. Ogino, R. Nishihara

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): K.A. Bertrand, K. Mima, A.T. Chan, K. Ng, S. Ogino, R. Nishihara

Study supervision: K. Wu, C.S. Fuchs, A.T. Chan, S. Ogino, R. Nishihara

The authors thank hospitals and pathology departments throughout the United States for generously providing them with tissue specimens. In addition, the authors thank the participants and staff of the NHS and the HPFS for their valuable contributions as well as the following state cancer registries for their help: AL, AZ, AR, CA, CO, CT, DE, FL, GA, ID, IL, IN, IA, KY, LA, ME, MD, MA, MI, NE, NH, NJ, NY, NC, ND, OH, OK, OR, PA, RI, SC, TN, TX, VA, WA, WY.

This work was supported by NIH grants (P01 CA87969 to S.E. Hankinson; P01 CA55075 to W.C. Willett; UM1 CA167552 to W.C. Willett; P50 CA127003 to C.S. Fuchs; R01 CA137178 to A.T. Chan; R01 CA151993 to S. Ogino; R01 CA136950 to E. Cho; R25 CA098566 to M.J. Stampfer, and K07 CA148894 to K. Ng), by grants from the American Society of Clinical Oncology (ASCO) Conquer Cancer Foundation (Career Development Award to K. Ng), Bennett Family Fund and the Entertainment Industry Foundation through National Colorectal Cancer Research Alliance and by postdoctoral fellowships from Asan Medical Center (to S.A. Kim), Japan Society for Promotion of Sciences (to K. Inamura), Takashi Tsuruo Memorial Fund (to K. Inamura), and the Uehara Memorial Foundation (to K. Mima).

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.