Among 1017 eligible participants, there were 283 deaths, 119 of which were due to colorectal cancer. The median time of follow-up of participants who were still alive was 116 months (range 41–238 months). The median predicted 25(OH)D level was 27.18

ng

ml
−1 in NHS and 29.18

ng

ml
−1 in HPFS. This pattern is consistent with what we saw in our previous analysis of pre-diagnosis plasma 25(OH)D levels and colorectal cancer survival, in which circulating 25(OH)D concentrations were also slightly higher in men than in women: median plasma 25(OH)D level was 27.1

ng

ml
−1 in men and 23.9 and 25.7

ng

ml
−1 for two laboratory runs in women (
Ng et al, 2008). The difference in mean post-diagnosis predicted 25(OH)D levels between the highest and lowest deciles in this study cohort was 10

ng

ml
−1, which is similar to the difference in mean actual circulating 25(OH)D levels across extreme deciles of predicted 25(OH)D in the original and validation cohorts.
Baseline characteristics according to quintiles of post-diagnosis predicted 25(OH)D levels are shown in . Overall, participants with higher predicted 25(OH)D levels were more likely to be of white race, have a lower BMI, report higher physical activity, have higher calcium intake, and were more likely to be diagnosed with colorectal cancer in the summer or autumn. Other prognostic characteristics did not differ significantly between quintiles.
| Table 2Baseline characteristics of cohort according to quintile of predicted 25(OH)D (n=1017) |
Higher post-diagnosis predicted 25(OH)D levels were associated with a significant reduction in colorectal cancer-specific and overall mortality (, respectively, and ). This relationship remained largely unchanged after adjusting for other predictors of cancer survival (). Compared with patients with post-diagnosis 25(OH)D scores in the lowest quintile, those in the highest quintile had an adjusted HR of 0.50 (95% CI, 0.26–0.95; P trend=0.02) for cancer-specific mortality and 0.62 (95% CI, 0.42–0.93; P trend=0.002) for overall mortality.
| Table 3Age-adjusted and multivariate hazard ratios of death according to quintile of predicted 25(OH)D in the entire study cohort (n=1017) |
We considered the possibility that the 25(OH)D score may be acting as a surrogate for the causal factor, such as BMI or physical activity, which are both in the prediction equation. We therefore repeated our analyses after adjusting for BMI and physical activity. When BMI was included, the significant relationship between post-diagnosis 25(OH)D score and cancer-specific and overall mortality did not change. The adjusted HR was 0.51 (95% CI, 0.26–0.99; P trend=0.04) for cancer-specific mortality and 0.62 (95% CI, 0.41–0.94; P trend=0.005) for overall mortality, comparing extreme quintiles. Similarly, when we controlled for physical activity, we continued to observe a benefit for higher post-diagnosis predicted 25(OH)D levels, with an adjusted HR of 0.44 (95% CI, 0.22–0.87; P trend=0.01) for cancer-specific mortality and 0.67 (95% CI, 0.44–1.00; P trend=0.02) for overall mortality. When both BMI and physical activity were included in the model, the adjusted HR was 0.45 (95% CI, 0.22–0.91; P trend=0.02) for cancer-specific mortality and 0.66 (95% CI, 0.43–1.03; P trend=0.03) for overall mortality.
Moreover, when both BMI and physical activity were included in our model, the remaining components of the post-diagnosis 25(OH)D score that were not ‘accounted for' were region of residence, race, and dietary and supplemental vitamin D intake. We therefore explored the impact of each of these individual variables on mortality in our study cohort, and found that patients who reported higher total vitamin D intake showed a trend towards lower risk of death (P trend=0.08). Compared with those in the lowest quintile, patients in the highest quintile of vitamin D intake had an adjusted HR of 0.72 (95% CI, 0.49–1.04) for overall mortality.
We also adjusted for calcium intake in our models and found a persistent significant association between post-diagnosis predicted 25(OH)D levels and survival. The adjusted HR was 0.44 (95% CI, 0.23–0.87; P trend=0.008) for cancer-specific mortality and 0.56 (95% CI, 0.37–0.84; P trend=0.0005) for overall mortality, comparing extreme quintiles. Of note, the addition of race and season of diagnosis to the multivariable model also did not change the adjusted HRs for cancer-specific and overall mortality. When race was added to the model, the adjusted HR comparing extreme quintiles was 0.53 (95% CI, 0.28–1.02; P trend=0.04) for cancer-specific mortality and 0.65 (95% CI, 0.43–0.97; P trend=0.004) for overall mortality. When season of diagnosis was included, the adjusted HR was 0.49 (95% CI, 0.26–0.94; P trend=0.02) for cancer-specific mortality and 0.63 (95% CI, 0.42–0.93; P trend=0.004) for overall mortality.
Given that lower levels of post-diagnosis predicted 25(OH)D could reflect the presence of occult cancer or other major illness, we excluded patients who died within 6 months of their post-diagnosis assessment. We continued to observe significant reductions in the risk of cancer-specific and overall mortality with increasing post-diagnosis 25(OH)D scores, with participants in the highest quintile having an adjusted HR of 0.50 (95% CI, 0.25–0.98; P trend=0.03) for cancer-specific mortality and 0.63 (95% CI, 0.42–0.95; P trend=0.003) for overall mortality.
In a separate analysis, pre-diagnosis 25(OH)D scores were calculated for colorectal cancer patients with available information (n=1955). Higher pre-diagnosis 25(OH)D scores were found to be associated with a decrease in cancer-specific and overall mortality (P trend=0.03 and 0.01, respectively). When we adjusted for pre-diagnosis predicted 25(OH)D levels as well as other predictors of cancer survival in our model, higher post-diagnosis predicted 25(OH)D levels were still associated with a significant reduction in both cancer-specific (P trend=0.02) and overall (P trend=0.008) mortality, whereas the effect of pre-diagnosis predicted 25(OH)D level was no longer significant.
We examined the influence of post-diagnosis 25(OH)D scores across the other predictors of cancer mortality (). Stratified analyses of cancer-specific and overall survival showed no significant interactions; the inverse relationship between post-diagnosis predicted 25(OH)D levels and mortality remained largely unchanged across most subgroups. Of note, there was a trend towards a greater impact of higher post-diagnosis 25(OH)D scores on overall survival in patients diagnosed in the winter or spring compared with those diagnosed in summer or autumn (P interaction=0.06).