Association of the CYP24A1-rs2296241 polymorphism of the vitamin D catabolism enzyme with hormone-related cancer risk: a meta-analysis

Abstract

Background

The evidence for vitamin D reducing cancer risk is inconsistent, and it is not clear whether this reduction is related to variation in cytochrome P450 (CYP)24A1, the only enzyme known to degrade active vitamin D. We focused on evaluating the association of CYP24A1-rs2296241 polymorphism with hormone-related cancer risk by conducting a meta-analysis.

Methods

A systematic literature search was conducted in April 2014 (updated in December 2014) to identify eligible studies. A random-effects model was used to pool the odds ratio (OR).

Results

Eleven studies including 5,145 cases and 5,136 controls were considered for the allelic model, and eight studies of 3,959 cases and 3,560 controls were utilized for the additive, recessive, and dominant models. There was no significant association between CYP24A1-rs2296241 and hormone-related cancer risk in any of the models, yet substantial heterogeneity was observed. Subgroup analyses indicated that CYP24A1-rs2296241 variation reduced the prostate cancer risk in the additive (OR 0.91, 95% confidence interval 0.85–0.97) and recessive (OR 0.80, 95% confidence interval 0.67–0.95) models, with no evidence of heterogeneity.

Conclusion

This meta-analysis indicated that CYP24A1-rs2296241 polymorphism reduced the androgen-related prostate cancer risk in additive and recessive models. More genetic loci are needed to confirm the effect of CYP24A1 variation on the risk of prostate cancer.

Keywords:

• CYP24A1
• hormone-related cancer
• meta-analysis
• polymorphism
• vitamin D

Introduction

The scientific evidence linking vitamin D with carcinogenesis is increasing, but data suggesting that high levels of 25-hydroxyvitamin D3 [25(OH)D] reduce the risk of cancer are inconsistent. This could be because measured serum 25(OH)D levels may not correspond to vitamin D exposure or reflect tissue-specific levels. Before vitamin D can act as an anticancer agent, it must go through two steps of hydroxylation. Vitamin D is hydroxylated by 27-hydroxylase (cytochrome P450 [CYP]27A1 to 25(OH)D in the liver, followed by further metabolism by 1α-hydroxylase [CYP]27B1) to 1α,25-dihydroxyvitamin D3 (1,25(OH)₂D₃).¹ The level of the main circulating form of vitamin D, 25(OH)D, is a widely accepted biomarker of vitamin D status, and 1,25(OH)₂D₃ is the activity of one that influence metabolic pathway, cell functions, and expression of target genes.² The final step in vitamin D metabolism is degradation of 25(OH)D and 1,25(OH)₂D₃ to 24,25(OH)D₃ and 1,24,25(OH)₂D₃, respectively, which occurs via 25-hydroxyvitamin D 24-hydrolase (encoded by the CYP24A1 gene).³ CYP24A1 has been clearly established as the main enzyme responsible for the degradation of active vitamin D. It is evident that CYP24A1 works in balance with CYP27B1. Elevated CYP24A1 expression and a reduced rate of CYP24A1 gene silencing has been reported in specific tumors, including the pancreas,⁴ lung cancer,^5–8 prostate cancer,⁹ colon cancers,^10–12 oral cancer,¹³ and non-Hodgkin’s lymphoma.¹⁴

CYP24A1 is located on 20q13.2 of chromosome 20; a number of polymorphisms of CYP24A1 have been identified, and the list is growing rapidly.^8,9 Epidemiological studies have investigated the association between CYP24A1 variation and several hormone-related cancers, including thyroid carcinoma,¹⁵ prostate cancer,^16–19 and breast cancer.^20,21 However, the results have been inconsistent, possibly because of small sample sizes, low statistical power, and clinical heterogeneity. In addition, there are overall 45 genetic loci of CYP24A1 in these hormone-related cancers, and 42 of them have reported only once. Locus rs2762941 has been studied in two articles, while rs927650 has been studied in five articles but two of them lacked original data. Only rs2296241 was found in all the studies, so we assessed the association between the CYP24A1-rs2296241 variant and hormone-related cancer risk by conducting a meta-analysis.

Methods

Study eligibility and selection

A systematic literature search was conducted in April 2014 (updated in December 2014) via PubMed, using the terms “CYP24A1”, “Vitamin D” or “25(OH)D” or “1,25(OH)₂D₃”, “cancer” or “carcinoma” as keywords. No restrictions were imposed. All eligible original studies, review articles, and other relevant studies were searched manually.

The criteria for inclusion were: study subjects were enrolled in case-control studies; the clinical outcome was hormone-related cancers; the exposure of interest was the CYP24A1-rs2296241 polymorphism; numbers of cases and controls with G/A, GG/GA, and GG/AA were provided, and genotype frequencies in the controls departed from Hardy–Weinberg equilibrium; and the odds ratio (OR) and corresponding 95% confidence interval (CI) were reported, or data to calculate them were included.

Data extraction

Two authors independently extracted data and collected the following: author, duration of follow-up, country, type of cancer, total number of cases and controls, and the number of cases and controls with G/A, GG/GA, and GG/AA for CYP24A1-rs2296241. The quality of each study was independently assessed by two authors using the Newcastle-Ottawa Scale.²² The quality score was determined using three blocks of eight entries: selection of study population, contrast between case and control, and measurement of exposure factors. Total score ranged from 0 (worst) to 9 (best). A study was considered of high quality if the score was ≥5.

Statistical analysis

The meta-analysis was performed using allelic contrast (G versus A), recessive (GG versus GA+AA), dominant (GG+GA versus AA), and additive (GG versus AA) models. Point estimates of risk, and the OR and 95% CI were determined using a random-effects model. Hardy–Weinberg equilibrium was examined using a chi-square test. I² was used to evaluate for statistical heterogeneity.^23,24 Potential publication bias was assessed by examining funnel plots and using Begg’s test and Egger’s test.^25,26 All statistical analyses were performed with STATA version 10.0 (Stata Corporation, College Station, TX) and SAS version 9.1 (SAS Institute, Cary, NC, USA).

Results

We identified 97 potentially relevant articles in our initial search. Of these, 44 papers on laboratory studies were excluded after the first screening based on the abstract or title, leaving 53 articles for full-text review. We excluded 22 articles for not being case-control studies. Among the remaining 31 articles, five were excluded for only the survey of cancer tissues, three for the observation of tumor recurrence and death, ten for lacking gene locus data, four for not having pooled data, and two for not including hormone-related cancers. Finally, seven publications representing eleven studies were included in the meta-analysis. A flow chart showing the study selection is presented in Figure 1.

Figure 1 Flowchart of study selection.

The characteristics of the eligible studies are shown in Table 1. The allelic contrast model included all eleven studies from the seven publications. In the other three models, six publications representing eight studies and three studies of prostate cancer¹⁶ were excluded because they lacked detailed data for further evaluation. These studies were published between 2006 and 2014, and included two breast cancer,^20,21 seven prostate cancer,^16–19 and two thyroid carcinoma papers.¹⁵ Two studies were conducted in Germany, six in the USA, two in Canada, and one in Korea.

Table 1 Characteristics of case-control studies included in the meta-analysis

Study	Reference	Ethnicity	Country	Case/control	Case			Control			OR	HWE Y/N (P-value)
					AA	AG	GG	AA	AG	GG	(95% CI)
PTC	Penna-Martinez et al^Citation15	German	Germany (2006–2008)	205/302	57	101	47	89	151	62	0.922 (0.72–1.19)	Y (0.87)
FTC	Penna-Martinez et al^Citation15	German	Germany (2007–2009)	48/302	9	22	17	89	151	62	0.597 (0.39–0.93)	Y (0.87)
Breast cancer	Anderson et al^Citation20	Colored	Canada (2002–2003)	1,777/1,839	449	777	330	468	791	371	1.034 (0.94–1.14)	Y (0.28)
Breast cancer	McCullough et al^Citation21	Caucasian, African-American	USA (1998–2001)	500/500	141	254	99	139	253	98	1.003 (0.84–1.20)	Y (0.41)
Prostate cancer	Holt et al^Citation18	Caucasian or African-American	Canada (2002–2005)	711/718	170	356	166	183	371	151	0.921 (0.79–1.07)	Y (0.15)
Prostate cancer	Holt et al^Citation18	Caucasian or African-American	USA (2002–2005)	116/68	37	50	25	15	39	12	1.13 (0.74–1.75)	Y (0.1)
Prostate cancer	Holick et al^Citation17	Caucasian or African-American	USA (1993–1996)	630/565	152	285	134	157	275	107	0.88 (0.75–1.04)	Y (0.54)
Prostate cancer	Beuten et al^Citation16	non-Hispanic Caucasians	USA (2001–2009)	609/348	48.4 (G%)	51.6 (A%)		47 (G%)	53 (A%)		0.95 (0.73–1.24)	–
Prostate cancer	Beuten et al^Citation16	Hispanic Caucasian	USA (2001–2009)	195/514	55.4 (G%)	44.6 (A%)		55.9 (G%)	44.1 (A%)		1.02 (0.73–1.42)	–
Prostate cancer	Beuten et al^Citation16	African-American	USA (2001–2009)	82/109	54.9 (G%)	45.1 (A%)		49.7 (G%)	50.3 (A%)		0.81 (0.5–1.33)	–
Prostate cancer	Oh et al^Citation19	Korean	Korea (1995–2002)	272/173	64 (G%)	36 (A%)		55.8 (G%)	44.2 (A%)			Y (>0.05)

Abbreviations: PTC, papillary thyroid carcinoma; FTC, follicular thyroid carcinoma; HWE, Hardy–Weinberg equilibrium; CI, confidence interval; OR, odds ratio.

The quality assessment showed that the quality scores ranged from 5 to 7 with a median score of 6, suggesting that all studies were of high quality (Table S1).

We found no significant association between CYP24A1-rs2296241 polymorphism and hormone-related cancer risk in any of the models (additive, OR 0.94, 95% CI 0.86–1.02; recessive, OR 0.86, 95% CI 0.72–1.01; allelic contrast, OR 1.03, 95% CI 0.96–1.11; dominant, OR 0.91, 95% CI 0.79–1.04; Figures 25). However, substantial heterogeneity was observed (I²>48.4%). Therefore, hormone-related cancers were further divided into thyroid cancer, breast cancer, and prostate cancer, and a sensitivity analysis was performed (Table S2).

Figure 2 Forest plot showing the association between CYP24A1-rs2296241 polymorphism and hormone-related cancer risk in the additive model.

Note: Weights are from random effects analysis.
Abbreviations: CI, confidence interval; CYP, cytochrome P450; TSH, thyroid-stimulating hormone; FTC, follicular thyroid carcinoma; PTC, papillary thyroid carcinoma.

For thyroid cancer, including 353 cases and 305 controls, there was no statistically significance between the CYP24A1-rs2296241 polymorphism and cancer risk for any model (Figures 25), and heterogeneity was still present in all models (all I²>50%). For breast cancer, including 2,277 cases and 2,339 controls, there was no association in any of the models (Figures 2–5), and no evidence of heterogeneity was observed (all I²=0%).

However, the CYP24A1-rs2296241 polymorphism significantly reduced the risk of prostate cancer, including 1,729 cases and 1,524 controls, in the additive (GG versus AA, OR 0.91, 95% CI 0.85–0.97, Figure 2) and recessive (GG versus GA+AA, OR 0.80, 95% CI 0.67–0.95, Figure 3) models, while no association was significantly observed in allelic contrast and dominant models (Figures 3 and 4). Heterogeneity was not observed in any model (all I²<25%), except for the dominant model (I²=52.2%).

Figure 3 Forest plot showing the association between CYP24A1-rs2296241 polymorphism and hormone-related cancer risk in the recessive mode.

Note: Weights are from random effects analysis.
Abbreviations: CI, confidence interval; CYP, cytochrome P450; TSH, thyroid-stimulating hormone; FTC, follicular thyroid carcinoma; PTC, papillary thyroid carcinoma.

Figure 4 Forest plot showing the association between the CYP24A1-rs2296241 polymorphism and hormone-related cancer risk in the allelic contrast mode.

Note: Weights are from random effects analysis.
Abbreviations: CI, confidence interval; CYP, cytochrome P450; TSH, thyroid-stimulating hormone; FTC, follicular thyroid carcinoma; PTC, papillary thyroid carcinoma.

Figure 5 Forest plot showing the association between CYP24A1-rs2296241 polymorphism and hormone-related cancer risk in the dominant mode.

Note: Weights are from random effects analysis.
Abbreviations: CI, confidence interval; CYP, cytochrome P450; TSH, thyroid-stimulating hormone; FTC, follicular thyroid carcinoma; PTC, papillary thyroid carcinoma.

We did not find evidence of publication bias with either Begg’s test or Egger’s test in any of the models (all Begg’s tests P>0.05, all Egger’s tests P>0.05). Therefore, it is unlikely that publication bias had a significant influence on the observed association between CYP24A1-rs2296241 variation and reduced risk of prostate cancer in the additive model (GG versus AA) or recessive model (GG versus GA+AA).

Discussion

To our knowledge, until now, there has been no reported meta-analysis evaluating the association between CYP24A1 polymorphism and hormone-related cancer risk. In this meta-analysis, we found that variation in CYP24A1-rs2296241 was not significantly associated with the risk of hormone-related cancer in four models with evidence of substantial heterogeneity.

Substantial heterogeneity was not surprising, given the differences in tumor site and hormone type. Sensitivity analyses revealed a significant association between variation in CYP24A1-rs2296241 and a decreased risk of prostate cancer in the additive and recessive models. This hormonal influence is hardly surprising because the significance of androgens in the development of prostate cancer has been known for more than half century and androgens are major contributors to prostatic carcinogenesis.²⁷ Moreover, 1,25(OH)₂D₃, the main target degraded by CYP24A1 and expressed in most nucleated cells in local tissue, can also influence prostaglandin synthesis. Prostaglandin promotes carcinogenesis and facilitates cancer progression.²⁸ Even so, the mechanisms underlying the association between variation in CYP24A1 and the risk of cancer are unclear. CYP24A1-rs2296241 (conservation score 0.99) is a synonymous single nucleotide polymorphism (SNP) in exon 4. Although the SNP does not change the amino acid sequence, it could influence regulation and expression of the CYP24A1 gene and thereby alter protein levels. In addition, mutations in human CYP24A1 have been shown to modulate the regioselectivity of the enzyme.²⁹ Hence, CYP24A1 is deemed to be a 1,25(OH)₂D₃-binding protein first and a catabolic enzyme second.³⁰ This enzyme is rate-limiting for the amount of local vitamin D in cancer tissue, and elevated expression is more likely to lead to an adverse prognosis in cancer. Another explanation is that the true functional SNPs near the promoter region 5′ of exon 1 in CYP24A1, including rs2296241, rs4809960, rs2585428, and rs6022999, have a functional impact on vitamin D response element binding, and thereby could counteract the antitumorigenic effects of 1,25(OH)₂D₃.^11,31,32 Therefore, there are many claims that CYP24A1 has been identified as a proto-oncogene.^7,8 In addition, serum androgen levels were negatively correlated with vitamin D status in men, and androgens are required for growth of prostate cancer in synergy with vitamin D receptor activity.³³ Therefore, CYP24A1-rs2296241 variation may be a reduced risk factor, meanwhile prolonged androgen stimulation may contribute to prostatic cancer.

Although there was no association between variation and thyroid cancer, there was significant heterogeneity in the four models. There was a small overlap in the CI for different thyroid-stimulating hormone-related cancers, and this could have affected the heterogeneity. Other than excess thyroid-stimulating hormone, pregnancies terminated by spontaneous or induced abortions can contribute to an increase in risk of thyroid cancer and significant heterogeneity.²⁷

We know that breast cancer is not only affected by estrogen and progesterone but also by age at menarche, menopause, and first full-term pregnancy.²⁷ Our subgroup analysis revealed no heterogeneity or association of risk between breast cancer and variation in CYP24A1-rs2296241, suggesting it is not a significant risk factor for breast cancer.

Several limitations in the current meta-analysis should be addressed. First, most studies had insufficient numbers, except for one,²⁰ which may have attenuated the statistical power, particularly for subgroup analysis. Second, because of the lack of original data, our results were based on unadjusted estimates of OR without adjustment for an individual’s disease history. Third, some heterogeneity existed that may be ascribed to roles of hormone and methods used for genotyping. Hence, we used the random-effects model that generated wider CIs for all genetic models. Finally, only one genetic locus of CYP24A1 variation was analyzed, because it appeared in all eleven studies.

The number of CYP24A1 polymorphisms in the genomic databases currently stands at around 50 and is growing rapidly. Two CYP24A1 polymorphisms (rs34043203 and rs2762934) have been associated with increased breast cancer risk and one (rs1570669) with reduced breast cancer risk.³⁵ There were a significant association of the CYP24A1-rs2296241 variant with a decreased risk of oral cancer¹³ and head and neck cancer.³⁴ Significantly altered risks of recurrence/progression were observed in relation to genotype for two CYP24A1 SNPs (rs927650 and rs2762939), and five CYP24A1 SNPs (rs3787557, rs4809960, rs2296241, rs2585428, and rs6022999) significantly altered the risk of death from prostate cancer.³⁶ High levels of CYP24A1 are also found in prostate cancer cells.^10,37 If overexpression of CYP24A1 directly affects tumor proliferation, tumor-targeted treatment with CYP24A1-specific inhibitors or mutagens may be effective in slowing tumor growth.

Conclusion

In summary, this meta-analysis indicates that variation in CYP24A1-rs2296241 may be associated with a reduced risk of androgen-related prostate cancer in additive (GG versus AA) and recessive (GG versus GA+AA) models. However, more genetic loci are needed to confirm the effect of CYP24A1 on the risk of androgen-related cancer.

Acknowledgments

This work was supported by grants from the National Natural Scientific Funding of China (81072286, 81372979), Jiangsu Natural Science Foundation (BK2012619), Chinese Nutrition Society Nutrition Research Foundation-DSM Research Fund, and Collaborative Innovation Center of Radiation Medicine, Jiangsu Higher Education Institution.

Supplementary materials

Table S1 Newcastle-Ottawa scale evaluation criteria of case-control study

Disease	Reference	Study population				Groups			Exposure	Total scores
		a	b	c	d	e	f	g	h
PTC	Penna-Martinez et al^Citation15	0	1	1	1	0	1	1	1	6
FTC	Penna-Martinez et al^Citation15	0	1	1	1	0	1	1	1	6
Breast cancer	Anderson et al^Citation20	1	0	1	1	1	0	1	0	5
Breast cancer	McCullough et al^Citation21	0	1	1	1	1	1	1	0	6
Prostate cancer	Holt et al^Citation18	1	1	1	1	1	1	1	0	7
Prostate cancer	Holt et al^Citation18	1	1	1	1	1	1	1	0	7
Prostate cancer	Holick et al^Citation17	1	0	1	1	0	1	1	0	5
Prostate cancer	Beuten et al^Citation16	0	1	0	1	1	1	1	1	6
Prostate cancer	Beuten et al^Citation16	0	1	0	1	1	1	1	1	6
Prostate cancer	Beuten et al^Citation16	0	1	0	1	1	1	1	1	6
Prostate cancer	Oh et al^Citation19	1	1	1	1	1	1	1	1	7

Notes: a, cases to determine whether appropriate; b, representative cases; c, select of control; d, determination of control; e, comparability of cases and controls; f, determination of exposure factors; g, using the same method to determine the cases and controls of exposure factors; h, measurement of exposure factors, no response rate.

Abbreviations: PTC, papillary thyroid carcinoma; FTC, follicular thyroid carcinoma.

Table S2 Meta-analysis of the association between CYP24A1-rs2296241 polymorphism and cancer risk

Model	Group	Studies	Test of association		Test of heterogeneity
		(n)	OR	95% CI	I^2	P-value
Additive	Thyroid cancera	2	0.78	(0.51–1.18)	60.2	0.113
	Breast cancer	2	1.14	(0.97–1.11)	0	0.695
	Prostate cancer	4	0.91	(0.85–0.97)	11.4	0.336
	Overall	8	0.94	(0.86–1.02)	56.3	0.025
Recessive	Thyroid cancer	2	0.67	(0.37–1.22)	57.4	0.125
	Breast cancer	2	1.07	(0.92–1.24)	0	0.597
	Prostate cancer	4	0.80	(0.67–0.95)	0	0.923
	Overall	8	0.9	(0.77–1.05)	48.4	0.059
Dominant	Thyroid cancer	2	0.8	(0.51–1.25)	26.7	0.243
	Breast cancer	2	1.01	(0.88–1.15)	0	0.99
	Prostate cancer	4	1.92	(0.66–1.28)	52.2	0.099
	Overall	8	0.91	(0.79–1.04)	36.4	0.138
Allelic	Thyroid cancer	2	1.29	(0.87–1.37)	65.1	0.09
	Breast cancer	2	1	(0.90–1.11)	0	1
	Prostate cancer	7	0.98	(0.86–1.12)	13.4	0.329
	Overall	11	1.02	(0.94–1.11)	20.6	0.253

Note:

a Includes papillary thyroid carcinoma and follicular thyroid carcinoma.

Abbreviations: CYP, cytochrome P450; CI, confidence interval; OR, odds ratio.

Disclosure

The authors report no conflicts of interest in this work.