Abstract
Background: Previous studies rarely evaluated the associations between vitamin D–binding protein and free vitamin D with colorectal cancer risk. We assessed these biomarkers and total 25-hydroxyvitamin D in relation to colorectal cancer risk in a sample of African Americans.
Methods: Cases comprised 224 African American participants of the Southern Community Cohort Study diagnosed with incident colorectal cancer. Controls (N = 440) were selected through incidence density sampling and matched to cases on age, sex, and race. Conditional logistic regression was used to calculate odds ratios (OR) and 95% confidence intervals (CI) for associations between biomarker levels and colorectal cancer risk.
Results: Vitamin D was inversely associated with colorectal cancer risk where the OR per-SD increase in total and free 25-hydroxyvitamin D were 0.82 (95% CI, 0.66–1.02) and 0.82 (95% CI, 0.66–1.01), respectively. Associations were most apparent among cases diagnosed >3 years after blood draw: ORs for the highest tertile versus the lowest were 0.69 (95% CI, 0.21–0.93) for total 25-hydroxyvitamin D and 0.71 (95% CI, 0.53–0.97) for free 25-hydroxyvitamin D. Inverse associations were seen in strata defined by sex, BMI, and anatomic site, although not all findings were statistically significant. Vitamin D–binding protein was not associated with colorectal cancer risk.
Conclusions: Our findings suggest that total and free 25-hydroxyvitamin D may be inversely associated with colorectal cancer risk among African Americans.
Impact: These findings highlight a potential role for vitamin D in colorectal cancer prevention in African Americans. Cancer Epidemiol Biomarkers Prev; 26(8); 1242–7. ©2017 AACR.
Introduction
Previous meta-analyses of prospective epidemiologic studies report an inverse association between total circulating 25-hydroxyvitamin D, the clinical measure of vitamin D status, and colorectal cancer risk (1, 2). Most previous studies have focused exclusively on total 25-hydroxyvitamin D or vitamin D3 in association with cancer risk, potentially missing important associations between bioavailable 25-hydroxyvitamin D, free 25-hydroxyvitamin D, and colorectal cancer risk. Bioavailable 25-hydroxyvitamin D is defined as the portion of 25-hydroxyvitamin D not bound by vitamin D–binding protein, and it accounts for approximately 10% to 15% of total circulating 25-hydroxyvitamin D. Free 25-hydroxyvitamin D constitutes less than 1% of total 25-hydroxyvitamin D in circulation and is completely unbound from albumin and vitamin D–binding protein (3–5). Bioavailable and free 25-hydroxyvitamin D are hypothesized to have higher biologic potential to act on cells than bound 25-hydroxyvitamin D (4, 6, 7). However, recent studies have found mixed results for the associations between these biomarkers and colorectal cancer risk. One study reports strong inverse associations between colorectal cancer risk and free and bioavailable 25-hydroxyvitamin D (8), whereas other studies find null or modest positive associations between the biomarkers and colorectal cancer risk (9–11). In addition, no study of bioavailable or free 25-hydroxyvitamin D has been conducted in a sample with large numbers of African Americans, whose total 25-hydroxyvitamin D levels tend to be substantially lower than in people of European descent. Herein, we evaluate the associations of circulating vitamin D–binding protein, total and free 25-hydroxyvitamin D, and colorectal cancer risk in a study sample of African Americans.
Materials and Methods
Study population
Data for the current analysis arise from a nested case–control study within the Southern Community Cohort Study (SCCS), a previously described prospective cohort study designed to investigate health disparities (12, 13). The SCCS enrolled more than 85,000 participants from 2002 to 2009 from 12 states in the southeastern United States Eligible participants were ages 40 to 79 at enrollment, English-speaking, and provided a blood sample to the study for biomarker measurement. Participants provided information on lifestyle factors, demographics, and personal medical history, primarily via in-person interviews. Cases for the current study were identified through linkages with state cancer registries through April 1, 2014, as diagnosed with an incident colorectal cancer after study enrollment. Colorectal cancer was defined by International Classification of Diseases-Oncology-3 codes C180-189, C199, and C209. We identified controls through incidence density sampling of the cohort among subjects who were free of any cancer diagnosis except skin cancer at the time of the matching case's diagnosis. We individually matched cases and controls by age, sex, and race at a ratio of two controls for each case. Because of small sample size of white participants (68 cases), analysis was restricted to African American participants for a total of 224 colorectal cancer cases and 440 controls. The SCCS was approved by institutional review boards at Vanderbilt University Medical Center and Meharry Medical College. All participants provided written informed consent.
Biomarker assessment
Circulating vitamin D and vitamin D–binding protein levels were measured at Heartland Assays, Inc. (Ames, Iowa). The method for quantitative determination of 25-hydroxyvitamin D was an FDA approved direct, competitive chemiluminescence immunoassay using the DiaSorin LIAISON 25-OH Vitamin D Total assay (14). This assay is co-specific for 25-hydroxyvitamin D3 and 25-hydroxyvitamin D2. The vitamin D–binding protein measurements were completed using the polyclonal Human Vitamin D–Binding Protein ELISA kit manufactured by Genway Biotech Inc. A proxy for free circulating 25-hydroxyvitamin D was calculated as the molar ratio of 25-hydroxyvitamin D to vitamin D–binding protein multiplied by 105 (15).
Statistical analysis
Frequency distributions of participant characteristics were tabulated by case–control status. Because of possible differences in sun exposure by season of blood sample draw, we created variables for 25-hydroxyvitamin D and free 25-hydroxyvitamin D that accounted for calendar week of sample collection. In separate statistical models, we regressed total and free 25-hydroxyvitamin D against calendar week using the periodic function described by Gail and colleagues (16). Seasonal-adjusted total and free 25-hydroxyvitamin values were created by adding the residuals from the periodic function model to the mean values of total and free 25-hydroxyvitamin D derived from the periodic function models. Season of participant recruitment did not vary by case status. Conditional logistic regression was used to calculate odds ratios (OR) and 95% confidence intervals (CI) for the associations between biomarker levels and colorectal cancer risk. Conditional regression statistical models included the following matched variables and potential confounders: age (5 year groupings), sex, body mass index (BMI; <18.5, 18.5–24.9, 25.0–29.9, 30.0–34.9, ≥35.0 kg/m2), education (<9 years, 9–11 years, high school, some college, college graduate and beyond), smoking (never, former, and current), physical activity (tertiles), alcohol intake (women: none, 0< drink/day ≤1, >1 drink/day; men: none, 0< drinks/day ≤2, >2 drinks/day), history of colorectal cancer screening (yes, no/unknown), and family history of colorectal cancer diagnosis in a first degree relative (yes, no/unknown). Missing covariate data (typically for <2%–4% of participants) were set to sex-specific medians. We calculated associations between vitamin D biomarkers and colorectal cancer risk with the biomarkers categorized in tertiles based on the distribution among controls and by per standard deviation change. P values for trend tests were calculated by treating ordinal vitamin D biomarker variables as continuous in the statistical model. Additional analyses stratified by length of time between blood sample donation and case diagnosis were conducted to define a potential biologically relevant exposure period and to assess the hypothesis that colorectal cancers diagnosed soon after study recruitment might not be attributable to baseline vitamin D status. Associations between vitamin D biomarkers and colorectal cancer risk were further evaluated using conditional logistic regression statistical models among subgroups defined by sex, BMI, and anatomic site (colon/rectum). Potential interactions between vitamin D–binding protein and vitamin D levels were evaluated by completing likelihood ratio tests comparing statistical models with and without the addition of cross-product terms. Statistical analyses were performed using SAS statistical software (version 9.4; SAS Institute Inc.).
Results
Cases and controls were similar for matched factors as well as for BMI and alcohol use (Table 1). In comparison to controls, cases tended to have less education, and lower household income. Cases were less likely to have been screened for colorectal cancer and were less likely to be current smokers. Cases and controls had similar median values for vitamin D–binding protein (301.4 and 300.6 μg/mL), total 25-hydroxyvitamin D (12.5 and 13.0 ng/mL), and free 25-hydroxyvitamin D (5.8 and 6.4). African Americans had lower 25-hydroxyvitamin D biomarker values than their white counterparts; the median value for total 25-hydroxyvitamin D and free 25-hydroxyvitamin D among white controls (N = 122, not included in primary analysis) were 18.6 ng/mL and 9.0. Vitamin D–binding protein levels were similar for both races (300.6 μg/mL for African Americans and 299.3 μg/mL for whites). Total and free 25-hydroxyvitamin D were highly correlated with each other (Pearson correlation coefficient = 0.96). Vitamin D–binding protein was not correlated with total 25-hydroxyvitamin D (Pearson correlation coefficient = 0.01) and was negatively correlated with free 25-hydroxyvitamin D (Pearson correlation coefficient = −0.23).
Selected characteristics of African American participants at baseline by case status, SCCS, 2002–2009
Baseline characteristics, N (%) . | Controls (N = 440) . | Cases (N = 224) . |
---|---|---|
Age, years, median (IQR) | 55 (13) | 55 (13) |
Male sex | 203 (46.1) | 100 (44.6) |
Education | ||
<High school | 163 (37.1) | 104 (46.4) |
High school | 147 (33.4) | 68 (30.4) |
>High school | 120 (27.3) | 47 (21.0) |
Household income, $ | ||
<15,000 | 275 (62.5) | 148 (66.1) |
15,000–49,999 | 136 (30.9) | 67 (29.9) |
≥50,000 | 16 (3.6) | 2 (0.9) |
Colorectal cancer screening, ever | 121 (27.5) | 55 (24.6) |
Family history of colorectal cancer | 22 (5.0) | 12 (5.4) |
BMI, kg/m2, median (IQR) | 28.3 (9.9) | 29.4 (10.0) |
Current smoker | 184 (41.8) | 75 (33.5) |
Non- and moderate alcohol consumers | 337 (76.6) | 175 (78.1) |
Physical activity, MET-h/d, median (IQR) | 16.3 (22.7) | 15.9 (21.0) |
Sedentary time, hours, median (IQR) | 8.0 (6.2) | 8.0 (6.3) |
Total energy, kcal, median (IQR) | 2,212 (1748) | 2,261 (1826) |
Meat, servings/day, median (IQR) | 2 (2) | 2 (1) |
Fruits and vegetable, servings/day, median (IQR) | 3 (2) | 3 (2) |
Folate, micrograms/day, median (IQR) | 432 (373) | 471 (388) |
Vitamin D supplement use | 0 (0) | 0 (0) |
25-hydroxyvitamin D, ng/mL, median (IQR) | 13.0 (10.2) | 12.5 (10.7) |
Vitamin D–binding protein, μg/mL, median (IQR) | 300.6 (56.9) | 301.4 (62.2) |
Baseline characteristics, N (%) . | Controls (N = 440) . | Cases (N = 224) . |
---|---|---|
Age, years, median (IQR) | 55 (13) | 55 (13) |
Male sex | 203 (46.1) | 100 (44.6) |
Education | ||
<High school | 163 (37.1) | 104 (46.4) |
High school | 147 (33.4) | 68 (30.4) |
>High school | 120 (27.3) | 47 (21.0) |
Household income, $ | ||
<15,000 | 275 (62.5) | 148 (66.1) |
15,000–49,999 | 136 (30.9) | 67 (29.9) |
≥50,000 | 16 (3.6) | 2 (0.9) |
Colorectal cancer screening, ever | 121 (27.5) | 55 (24.6) |
Family history of colorectal cancer | 22 (5.0) | 12 (5.4) |
BMI, kg/m2, median (IQR) | 28.3 (9.9) | 29.4 (10.0) |
Current smoker | 184 (41.8) | 75 (33.5) |
Non- and moderate alcohol consumers | 337 (76.6) | 175 (78.1) |
Physical activity, MET-h/d, median (IQR) | 16.3 (22.7) | 15.9 (21.0) |
Sedentary time, hours, median (IQR) | 8.0 (6.2) | 8.0 (6.3) |
Total energy, kcal, median (IQR) | 2,212 (1748) | 2,261 (1826) |
Meat, servings/day, median (IQR) | 2 (2) | 2 (1) |
Fruits and vegetable, servings/day, median (IQR) | 3 (2) | 3 (2) |
Folate, micrograms/day, median (IQR) | 432 (373) | 471 (388) |
Vitamin D supplement use | 0 (0) | 0 (0) |
25-hydroxyvitamin D, ng/mL, median (IQR) | 13.0 (10.2) | 12.5 (10.7) |
Vitamin D–binding protein, μg/mL, median (IQR) | 300.6 (56.9) | 301.4 (62.2) |
NOTE: Subjects with data missing for the characteristic of interest are not included in this analysis.
Abbreviation: IQR, interquartile range.
ORs for the associations between total and free 25-hydroxyvitamin D with colorectal risk indicated inverse associations, although all confidence intervals crossed unity (Table 2). A null association was observed between vitamin D–binding protein and colorectal cancer risk. Results were equivalent from minimally and fully adjusted statistical models.
Associations between vitamin D biomarkers and colorectal cancer risk, African American participants of the SCCS
Biomarker . | Controls . | Cases . | ORa (95% CI) . | Ptrend . | ORb (95% CI) . | Ptrend . |
---|---|---|---|---|---|---|
25(OH)D, ng/mL | ||||||
Tertile 1: ≤10.49 | 102 | 71 | 1 (Ref) | 1 (Ref) | ||
Tertile 2: 10.50–16.08 | 103 | 47 | 0.65 (0.41–1.04) | 0.62 (0.38–1.01) | ||
Tertile 3: >16.08 | 107 | 57 | 0.77 (0.48–1.21) | 0.73 (0.45–1.17) | ||
Per standard deviation increase | 0.87 (0.71–1.06) | 0.17 | 0.82 (0.66–1.02) | 0.07 | ||
VDBP, μg/mL | ||||||
Tertile 1: ≤280.20 | 146 | 74 | 1 (Ref) | 1 (Ref) | ||
Tertile 2: 280.21–320.00 | 145 | 69 | 0.93 (0.62–1.40) | 0.95 (0.63–1.45) | ||
Tertile 3: >320.00 | 149 | 81 | 1.10 (0.74–1.64) | 1.08 (0.72–1.63) | ||
Per standard deviation increase | 1.04 (0.88–1.23) | 0.66 | 1.03 (0.87–1.23) | 0.70 | ||
Free 25(OH)Dc | ||||||
Tertile 1: ≤5.02 | 103 | 60 | 1 (Ref) | 1 (Ref) | ||
Tertile 2: 5.03–7.77 | 102 | 57 | 0.92 (0.57–1.46) | 0.85 (0.53–1.39) | ||
Tertile 3: >7.77 | 107 | 58 | 0.90 (0.56–1.44) | 0.83 (0.51–1.36) | ||
Per standard deviation increase | 0.87 (0.71–1.06) | 0.16 | 0.82 (0.66–1.01) | 0.07 |
Biomarker . | Controls . | Cases . | ORa (95% CI) . | Ptrend . | ORb (95% CI) . | Ptrend . |
---|---|---|---|---|---|---|
25(OH)D, ng/mL | ||||||
Tertile 1: ≤10.49 | 102 | 71 | 1 (Ref) | 1 (Ref) | ||
Tertile 2: 10.50–16.08 | 103 | 47 | 0.65 (0.41–1.04) | 0.62 (0.38–1.01) | ||
Tertile 3: >16.08 | 107 | 57 | 0.77 (0.48–1.21) | 0.73 (0.45–1.17) | ||
Per standard deviation increase | 0.87 (0.71–1.06) | 0.17 | 0.82 (0.66–1.02) | 0.07 | ||
VDBP, μg/mL | ||||||
Tertile 1: ≤280.20 | 146 | 74 | 1 (Ref) | 1 (Ref) | ||
Tertile 2: 280.21–320.00 | 145 | 69 | 0.93 (0.62–1.40) | 0.95 (0.63–1.45) | ||
Tertile 3: >320.00 | 149 | 81 | 1.10 (0.74–1.64) | 1.08 (0.72–1.63) | ||
Per standard deviation increase | 1.04 (0.88–1.23) | 0.66 | 1.03 (0.87–1.23) | 0.70 | ||
Free 25(OH)Dc | ||||||
Tertile 1: ≤5.02 | 103 | 60 | 1 (Ref) | 1 (Ref) | ||
Tertile 2: 5.03–7.77 | 102 | 57 | 0.92 (0.57–1.46) | 0.85 (0.53–1.39) | ||
Tertile 3: >7.77 | 107 | 58 | 0.90 (0.56–1.44) | 0.83 (0.51–1.36) | ||
Per standard deviation increase | 0.87 (0.71–1.06) | 0.16 | 0.82 (0.66–1.01) | 0.07 |
Abbreviations: 25(OH)D, 25-hydroxyvitamin D; VDBP, vitamin D–binding protein; Ref, reference.
aAnalyses result from conditional logistic regression models where cases and controls are matched on age, race, and sex. Analyses are adjusted for calendar week of sample collection, and body mass index.
bIn addition, adjusted for education, smoking, physical activity, alcohol intake, history of colorectal cancer screening, and family history of colorectal cancer.
cFree 25(OH)D biomarker is calculated as 25(OH)D:VDBP molar ratio (×103) and is a proxy for free 25-hydroxyvitamin D status.
We carried out stratified analyses by the time between blood sample donation and case diagnosis to evaluate the possible influence of reverse causation and determine a biologically relevant window of exposure (Table 3, Supplementary Table). In analyses among cases diagnosed with colorectal cancer more than 3 years after blood draw and their matched controls, we found inverse associations between total and free 25-hydroxyvitamin D with colorectal cancer risk. Among participants with greater than 3 years between blood draw and follow-up, a per standard deviation increase in 25-hydroxyvitamin D resulted in a 31% (95% CI, 7%–79%) decreased risk of colorectal cancer. Further analyses among these cases and their controls showed that the inverse association was more apparent in women than men, and in obese (BMI ≥ 30 kg/m2) more than non obese participants, although interaction tests were not statistically significant (Fig. 1). The inverse association was somewhat stronger for colon cancer than rectal cancer (Pinteraction = 0.70). Analyses restricted to cases diagnosed within 3 years of blood draw and their controls produced positive but non-statistically significant associations between tertiles of total and free 25-hydroxyvitamin D and colorectal cancer risk; among this subgroup, when total and free 25-hydroxyvitamin D were modeled in per standard deviation increments, we observed null associations with colorectal cancer risk (Table 3). Results were similar, in analyses stratified by ≤ or > 2 years of time between blood draw and diagnosis (Supplementary Table). Vitamin D–binding protein was not associated with colorectal cancer risk in any of the stratified analyses (Table 3). In addition, we did not observe evidence of interaction between vitamin D–binding protein levels and either total or free 25-hydroxyvitamin D in risk of colorectal cancer (all Pinteraction > 0.05).
Associations between vitamin D biomarkers and colorectal cancer risk by time since blood draw, African American participants of the SCCS
. | Time between blood draw and cancer diagnosis . | |||||
---|---|---|---|---|---|---|
. | ≤3 years . | >3 years . | ||||
Biomarker . | Cases . | ORa (95% CI) . | Ptrend . | Cases . | ORa (95% CI) . | Ptrend . |
25(OH)D, ng/mL | ||||||
Tertile 1: ≤10.49 | 22 | 1 (Ref) | 49 | 1 (Ref) | ||
Tertile 2: 10.50–16.08 | 28 | 1.22 (0.62–2.38) | 19 | 0.36 (0.19–0.69) | ||
Tertile 3: >16.08 | 32 | 1.45 (0.75–2.79) | 25 | 0.45 (0.25–0.83) | ||
Per SD increase | 0.98 (0.76–1.28) | 0.90 | 0.69 (0.21–0.93) | 0.02 | ||
VDBP, μg/mL | ||||||
Tertile 1: ≤280.20 | 24 | 1 (Ref) | 50 | 1 (Ref) | ||
Tertile 2: 280.21–320.00 | 25 | 1.08 (0.57–2.03) | 44 | 0.91 (0.56–1.48) | ||
Tertile 3: >320.00 | 38 | 1.55 (0.86–2.82) | 43 | 0.85 (0.52–1.38) | ||
Per SD increase | 1.26 (1.00–1.60) | 0.05 | 0.89 (0.72–1.10) | 0.27 | ||
Free 25(OH)Db | ||||||
Tertile 1: ≤5.02 | 20 | 1 (Ref) | 40 | 1 (Ref) | ||
Tertile 2: 5.03–7.77 | 29 | 1.38 (0.71–2.71) | 28 | 0.60 (0.33–1.09) | ||
Tertile 3: >7.77 | 33 | 1.64 (0.83–3.23) | 25 | 0.51 (0.27–0.95) | ||
Per SD increase | 0.96 (0.74–1.25) | 0.76 | 0.71 (0.53–0.97) | 0.03 |
. | Time between blood draw and cancer diagnosis . | |||||
---|---|---|---|---|---|---|
. | ≤3 years . | >3 years . | ||||
Biomarker . | Cases . | ORa (95% CI) . | Ptrend . | Cases . | ORa (95% CI) . | Ptrend . |
25(OH)D, ng/mL | ||||||
Tertile 1: ≤10.49 | 22 | 1 (Ref) | 49 | 1 (Ref) | ||
Tertile 2: 10.50–16.08 | 28 | 1.22 (0.62–2.38) | 19 | 0.36 (0.19–0.69) | ||
Tertile 3: >16.08 | 32 | 1.45 (0.75–2.79) | 25 | 0.45 (0.25–0.83) | ||
Per SD increase | 0.98 (0.76–1.28) | 0.90 | 0.69 (0.21–0.93) | 0.02 | ||
VDBP, μg/mL | ||||||
Tertile 1: ≤280.20 | 24 | 1 (Ref) | 50 | 1 (Ref) | ||
Tertile 2: 280.21–320.00 | 25 | 1.08 (0.57–2.03) | 44 | 0.91 (0.56–1.48) | ||
Tertile 3: >320.00 | 38 | 1.55 (0.86–2.82) | 43 | 0.85 (0.52–1.38) | ||
Per SD increase | 1.26 (1.00–1.60) | 0.05 | 0.89 (0.72–1.10) | 0.27 | ||
Free 25(OH)Db | ||||||
Tertile 1: ≤5.02 | 20 | 1 (Ref) | 40 | 1 (Ref) | ||
Tertile 2: 5.03–7.77 | 29 | 1.38 (0.71–2.71) | 28 | 0.60 (0.33–1.09) | ||
Tertile 3: >7.77 | 33 | 1.64 (0.83–3.23) | 25 | 0.51 (0.27–0.95) | ||
Per SD increase | 0.96 (0.74–1.25) | 0.76 | 0.71 (0.53–0.97) | 0.03 |
Abbreviations: 25(OH)D, 25-hydroxyvitamin D; VDBP, vitamin D–binding protein; Ref, reference; SD, standard deviation.
aAnalyses result from conditional logistic regression models where cases and controls are matched on age, race, and sex. Analyses are adjusted for calendar week of sample collection, body mass index, education, smoking, physical activity, alcohol intake, history of colorectal cancer screening, and family history of colorectal cancer.
bFree vitamin D biomarker is calculated as 25(OH)D:VDBP molar ratio (x103) and is a proxy for free 25-hydroxyvitamin D status.
Associations between free vitamin D and colorectal cancer risk by selected participant characteristics. Odds ratios are presented as a per standard deviation increase in level of free vitamin D and result from conditional logistic regression models where cases and controls are matched on race, age, and sex. Odds ratios and 95% confidence intervals are adjusted for calendar week of sample collection, body mass index, education, smoking, physical activity, alcohol intake, history of colorectal cancer screening, and family history of colorectal cancer. Analyses include African American cases with greater than 3 years between blood draw and diagnosis and their matched controls. Free vitamin D is calculated as 25(OH)D:VDBP molar ratio (×103) and is a proxy for free 25-hydroxyvitamin D status. a, P values for interaction were calculated by inclusion of cross-product terms for free vitamin D and the variable of interest. We were unable to calculate a P value for interaction between sex and free vitamin D status because sex was used as a matching factor.
Associations between free vitamin D and colorectal cancer risk by selected participant characteristics. Odds ratios are presented as a per standard deviation increase in level of free vitamin D and result from conditional logistic regression models where cases and controls are matched on race, age, and sex. Odds ratios and 95% confidence intervals are adjusted for calendar week of sample collection, body mass index, education, smoking, physical activity, alcohol intake, history of colorectal cancer screening, and family history of colorectal cancer. Analyses include African American cases with greater than 3 years between blood draw and diagnosis and their matched controls. Free vitamin D is calculated as 25(OH)D:VDBP molar ratio (×103) and is a proxy for free 25-hydroxyvitamin D status. a, P values for interaction were calculated by inclusion of cross-product terms for free vitamin D and the variable of interest. We were unable to calculate a P value for interaction between sex and free vitamin D status because sex was used as a matching factor.
Discussion
We find inverse associations between total and free 25-hydroxyvitamin D and colorectal cancer risk, with associations limited to cases diagnosed more than three years after their blood draw. Our study is the largest to date of the association of vitamin D and colorectal cancer risk in African Americans. Our results are consistent with potential anticarcinogenic effects of vitamin D among African Americans and raise the possibility that lower total 25-hydroxyvitamin D levels contribute to the higher risk of colorectal cancer among African Americans than other racial groups. Other lines of evidence support an inverse association between vitamin D and overall cancer risk. Animal and cancer cell line studies show the presence of vitamin D is associated with decreased cellular inflammation and proliferation, and increased cellular differentiation and apoptosis (17–21). In addition, the vitamin D receptor is expressed in normal colon and rectal cells (22), as is the vitamin D activating enzyme CYP27B1 (23), providing evidence that vitamin D has biologic function in these cell types.
Yet, studies in humans yield conflicting results. Secondary analyses from randomized controlled trials (RCT) observe null associations between vitamin D supplementation and colorectal cancer risk (18, 24). In addition, a 2015 publication describing a RCT of daily vitamin D supplementation of 1000 IU for 3 to 5 years reports no association with risk of recurrent colorectal adenoma (25). However, these RCTs have specific limitations, most notably a small number of colorectal cancers or adenomas, limited treatment adherence, short-term follow-up, and strict eligibility criteria (limited generalizability). In contrast to RCTs, observational epidemiologic evidence supports an inverse association between total circulating 25-hydroxyvitamin D and colorectal cancer risk. A 2011 meta-analysis of prospective epidemiologic studies computed an OR for colorectal cancer of 0.66 for participants in the highest versus lowest quartile of 25-hydroxyvitamin D levels (2). More recent publications that described data from the Women's Health Initiative and the Prostate, Lung, Colorectal and Ovarian Cancer Screening trials also report 25-hydroxyvitamin D to be inversely association with colorectal cancer risk (9, 26). Similar to these previous epidemiologic studies, we observe a decreased risk for participants in the highest versus lowest tertile of total 25-hydroxyvitamin D, and a consistent inverse association when total 25-hydroxyvitamin D is modeled in per standard deviation increments.
Several recent studies have investigated associations between colorectal cancer risk and vitamin D–binding protein, free, and bioavailable 25-hydroxyvitamin D. Our results confirm the null association between colorectal cancer risk and vitamin D–binding protein reported in four recent publications (8–11), and find an inverse association with free 25-hydroxyvitamin D more than 3 years after blood draw, which is consistent with results from one of the previous studies (9). One previous study (8), but not the others (9–11), reports inverse associations between total, free, and bioavailable 25-hydroxyvitamin D with colorectal cancer risk that were most apparent in participants with values of vitamin D–binding protein below the median. We do not find evidence of interaction between vitamin D–binding protein and total or free 24-hydroxyvitamin D. Anic and colleagues (10) report positive rather than inverse associations between total and free 25-hydroxyvitamin D with colorectal cancer risk using data from the Alpha-Tocopherol, Beta-Carotene Prevention (ATBC) trial. Of note, the ATBC study only includes male smokers. Cigarette smoking may affect vitamin D metabolism; one investigation finds a polycyclic aromatic hydrocarbon created during cigarette smoking inactivates the hormonal form of vitamin D in cell lines, which may have implications for the association between vitamin D and colorectal cancer risk in smokers (27).
Among participants with more than three years between blood draw and diagnosis, we evaluated whether the associations between vitamin D biomarkers and colorectal cancer risk varied by sex, BMI, or anatomic site. The reduced risk ORs were primarily seen in women and obese participants. Previous epidemiologic studies either find no difference in association by sex (28), or similar to our study, a more prominent association in women (9, 29). In contrast to our study, Wu and colleagues (29) observed a stronger association between total 25-hydroxyvitamin D and colorectal cancer risk in participants with BMI below the median value for study participants. Previous studies report mixed findings for whether the association between total circulating 25-hydroxyvitamin D and cancer risk varies by colon or rectal site. A meta-analysis reports a stronger inverse association in rectal than colon tumors; however, the interaction is not statistically significant at P < 0.05 (2), and the largest individual study included in the meta-analysis reports stronger inverse associations with colon than with rectal tumors (28).
In the SCCS, among participants with three or fewer years between blood draw and diagnosis there are null associations between total and free 25-hydroxyvitamin D with colorectal cancer risk. The presence of undiagnosed disease among cases diagnosed soon after baseline may have affected vitamin D biomarker status in these individuals. It is conceivable that participants with colorectal cancer may have changed their behavior in activities known to affect vitamin D levels, such as reduced physical activity, greater time indoors and lower sun exposure, shortly before their diagnosis. The change in activities may have led to the null association. In our study, more than three years between blood draw and diagnosis may represent a relevant biological window of exposure in the associations between circulating total and free 25-hydroxyvitamin D and colorectal cancer risk. The natural history of colorectal cancer suggests that it make take upwards of 10 years for an adenomatous polyp to develop to colorectal cancer (30). If our 25-hydroxyvitamin D biomarker categorization of high vs. low status correlates with long term exposure, then it is plausible that the association between 25-hydroxyvitamin D and colorectal cancer risk may be most evident in participants with longer follow-up time and more time for the process of carcinogenesis to develop.
Strengths of our study include the prospective study design, and the inclusion of vitamin D–binding protein, total, and free vitamin D in evaluating the association of vitamin D and colorectal cancer risk. Our study is one of the largest evaluations of vitamin D biomarkers and colorectal cancer risk in African Americans. Similar to the few previous studies that have included African American participants, we observe the inverse association between 25-hydroxyvitamin D biomarkers and colorectal cancer risk in African Americans (9, 31). Our study has limitations such that the relatively short follow-up time and limited sample size in subgroups left us without sufficient power to make definitive conclusions about the associations of vitamin D biomarkers with cancer risk in subgroups of white participants and those diagnosed with rectal cancer. Unfortunately, we did not have albumin measured for all participants consequently we were unable to assess the relationship between bioavailable vitamin D and colorectal cancer risk in this analysis, and we chose to use a proxy for free 25-hydroxyvitamin D instead of directly measuring it. However, previous studies have noted the high correlation between measured and calculated free 25-hydroxyvitamin D (32), which gives us confidence that our calculated proxy for free 25-hydroxyvitamin D is similar to actual circulating levels of the biomarker. In addition, because this is an observational study there is the potential for residual or uncontrolled confounding factors to influence our results, and we are unable to make causal inferences about the associations between vitamin D biomarkers and colorectal cancer risk.
Conclusions
Our prospective study finds inverse associations between total and free 25-hydroxyvitamin D and colorectal cancer risk among African American participants with longer follow up. The current study adds to the limited literature on vitamin D and cancer among African Americans. Importantly, there is a high correlation between total and free circulating 25-hydroxyvitamin D, and consequently similar associations between both measures and risk of colorectal cancer. The high correlation implies that measurement of total circulating 25-hydroxyvitamin D may be the appropriate vitamin D biomarker to measure in epidemiologic studies. The relationship between vitamin D and colorectal cancer risk has been understudied in African Americans. The present investigation provides one of the few opportunities to evaluate the influence of 25-hydroxyvitamin D measures on colorectal cancer risk at lower ranges than were examined in previous studies.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Disclaimer
The contents of this publication are solely the responsibility of the authors and do not necessarily represent the official views of the CDC or the Mississippi Cancer Registry. Cancer data for Southern Community Cohort Study cancer cases from West Virginia have been provided by the West Virginia Cancer Registry. The opinions expressed are those of the authors and do not necessarily represent those of the West Virginia Cancer Registry.
Authors' Contributions
Conception and design: S. Warren Andersen, X.-O. Shu, W. Zheng
Development of methodology: S. Warren Andersen
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): Q. Cai, M.D. Steinwandel, W.J. Blot
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): S. Warren Andersen, N.K. Khankari, W. Zheng
Writing, review, and/or revision of the manuscript: S. Warren Andersen, X.-O. Shu, Q. Cai, P.W. Jurutka, W.J. Blot, W. Zheng
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): S. Warren Andersen
Acknowledgments
We thank Ms. Kimberly A. Kreth and Ms. Nan Kennedy, Division of Epidemiology, Vanderbilt University Medical Center, for their assistance in preparing the article. Data on Southern Community Cohort Study cancer cases used in this publication were provided by the Alabama Statewide Cancer Registry; Kentucky Cancer Registry, Lexington, KY; Tennessee Department of Health, Office of Cancer Surveillance; Florida Cancer Data System; North Carolina Central Cancer Registry, North Carolina Division of Public Health; Georgia Comprehensive Cancer Registry; Louisiana Tumor Registry; Mississippi Cancer Registry; South Carolina Central Cancer Registry; Virginia Department of Health, Virginia Cancer Registry; Arkansas Department of Health, Cancer Registry, Little Rock, AR. Data on Southern Community Cohort Study cancer cases from Mississippi were collected by the Mississippi Cancer Registry which participates in the National Program of Cancer Registries of the Centers for Disease Control and Prevention.
Grant Support
The Southern Community Cohort Study is funded by grant R01 CA92447 (principal investigators: W.J. Blot and W. Zheng) from the National Cancer Institute at the NIH, including special allocations from the American Recovery and Reinvestment Act (3R01 CA09244708S1). S. Warren Andersen is supported by K99 CA207848 and the Vanderbilt Molecular and Genetic Epidemiology of Cancer training program (U.S. NIH grant R25 CA160056: awarded to X.-O. Shu). N.K. Khankari is supported by the Vanderbilt Molecular and Genetic Epidemiology of Cancer training program (NIH grant R25 CA160056: awarded to X.-O. Shu). The Arkansas Central Cancer Registry is fully funded by a grant from the National Program of Cancer Registries, Centers for Disease Control and Prevention.
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.