Association between CYP1A1 Ile462Val Polymorphism and Oral Squamous Cell Carcinoma Susceptibility: Evidence from 13 Investigations.

CYP1A1 Ile462Val polymorphism might play a key role in pathogenesis of oral squamous cell carcinoma (OSCC). Many case-control studies have investigated the association between CYP1A1 Ile462Val polymorphism and OSCC susceptibility. However, the conclusions are inconsistent. To aim a convincible conclusion, we carried out a meta-analysis to systematically evaluate the association of CYP1A1 Ile462Val polymorphism with OSCC susceptibility. We searched Pubmed, Web of Science, Ovid and Embase databases for available publications. The odds ratio (OR) with the corresponding 95% confidence interval (95% CI) was carried out to estimate the association. A total of 13 papers including 1468 cases and 2183 controls were included, a significant increased OSCC risk was observed in recessive model (OR=1.64, 95% CI=1.08-2.49), but not other genetic models. Our results suggest that the homozygous variant of CYP1A1 Ile462Val might be a risk factor of OSCC.

Introduction

Oral cancer is one of the most commonly occurred cancers worldwide and causes a considerable problem to human health being a high mortality rates and disfigurement 1, 2. About 90% of oral malignant neoplasms are oral squamous cell carcinomas (OSCC), followed by adenocarcinoma and, rarely, other types 3. Despite advances in treatment for OSCC, the 5-year survival rate remains undesirable 4-6. Hence, investigating the risk factors and developing the early diagnosis for treatment and prevention of OSCC are urgently needed.

Epidemiological studies showed that OSCC have been associated with high tobacco abuse and alcohol consumption 7-11. However, not all individuals who have smoking and alcohol habits develop this fatal disease, suggesting that individual genetic factor may also be involved in disease etiology. The research results of human genome project (HGP) have demonstrated that different individuals are all 99.9% the same genes and the tiny genetic difference sequences are mainly gene single nucleotide polymorphisms (SNPs). Therefore, inter-individual differences in expression of SNPs might contribute to variability in risk towards various types of cancers including OSCC. Currently, the published evidences showed that there were significant associations of gene polymorphisms with the susceptibility of many cancers, such as Glutathione S-transferase (GSTs) and Cytochrome P4501 A1 (CYP1A1) gene polymorphisms with lung of squamous cell carcinoma, polymorphism of 8q24 rsl3281615 with risk of breast cancer, CYP1A1 and GSTs gene polymorphisms with head and neck cancer 12-19. However, the associations of CYP1A1 Ile462Val variant with OSCC risks are inconsistent. For example, some studies 20-22 reported that CYP1A1*2C polymorphism was a risk factor of OSCC, however, some studies 23, 24 opposed it.

CYP1A1 is a member of the CYP family that participates in the metabolism of xenobiotics and endogenous compounds, encoding for the aryl hydrocarbon hydrolase (AHH), which is involved in activation of polycyclic aromatic hydrocarbon (PAHs) and aromatic amines and is expressed in oral tissue 25. Because CYP1A1 is able to activate carcinogenic PAHs and its expression and function are affected by gene polymorphisms, people pay more attention to the association of cancer and CYP1A1. According to the published studies, the CYP1A1 have several single nucleotide polymorphisms, which may alter the activities of their enzymes and increase carcinogen activation and yield to carcinogenicity. The CYP1A1*2C or CYP1A1 exon7 or Ile462Val is one of the most common polymorphisms, which is a transition from valine to isoleucine at codon 462 in the exon7 region resulting in the introduction of an BsrDI restriction site and associating with increase in enzyme activity and hence affecting the risks of carcinoma 26, 27. Exon7 restriction site polymorphism brings about three genotypes, wild type (Val/Val), heterozygous variant (Val/Ile) and homozygous variant (Ile/Ile) 28.

Considering the significance of CYP1A1 Ile462Val polymorphism in the occurrence and development of malignancies including OSCC, we systematically evaluated the association between CYP1A1 Ile462Val polymorphism and OSCC susceptibility through a meta-analysis.

Materials and Methods

Search strategy

We searched the following databases: Pubmed, Web of Science, Ovid and Embase databases without language limitations, and the last research was updated on August 8, 2014. The search process was designed to find primarily all relevant articles and the search strategies were listed as follows: (1) Cytochrome P450 1A1 or P450 1A1 or CYP1A1 or CYP1A1∗2C or exon7 or Ile462Val; and (2) genotype or polymorphism or allele or variant; and (3) oral squamous cell carcinoma or OSCC or mouth neoplasm or oral cancer or oral carcinoma or oral tumor. Searched results were in dependently screened by two authors according to the titles, abstracts and types of articles, and irrelevant papers were dropped out. Manual review of the references cited in the selected articles was undertaken to obtain articles that might have been missed in the search process. Then we downloaded the relevant papers and further screened to identify potentially eligible studies. If essential data were not provided in the original articles, every effort was made to contact the authors.

Inclusion/exclusion criteria

All relevant case-control studies were screened, irrespective of languages. In the meta-analysis, the following criteria were set and reviewed by two independent authors (Xiao-Lei Yang and Shang Xie): (1) studies should be concerned about the relation of CYP1A1 Ile462Val polymorphism with OSCC risk, and cases of identification of OSCC were histologically confirmed; (2) each trial should be an observational study (case-control or cohort), and the subjects should be human beings; (3) papers must offer the size of the sample, and the genetic distribution or the original information that can help infer the relevant results; (4) when multiple publications from a particular research group reported data from overlapping samples, the study reporting the largest or latest dataset was included.

Exclusion criteria: (1) review articles or editorials or meta-analysis; (2) case reports, or lack of case-control studies; (3) studies that estimated the risk of secondary tumors, prognosis, recurrence, or response to treatment. For conflicting evaluation, an agreement was reached following a discussion. If a consensus could not be attained, another author was invited to solve the dispute and then a final result was made by the majority of the votes. After rigorous searching, we viewed all papers in accordance with the criteria defined above for further analysis.

Data extraction

All data was independently reviewed and extracted from the included papers by two investigators. Differences between reviewers were solved by discussion, or through consultation if necessary. The following characteristic were collected from each study, ethnicity, country, sample size, control source, matching contents, the Hardy-Weinberg Equilibrium (HWE), and genetic distribution of cases and controls. If there are studies deviating from HWE, we will increase the analysis of subgroup based on the HWE. When the data were not clear nor presented by the author in the publication, contacting them for further details were attempted.

Quality assessment

We performed the Newcastle-Ottawa Scale (NOS) 29 to estimate the methodological quality of included studies. The NOS system categorizes into three dimensions including selection, comparability, and ascertainment of outcome. A star system is used to assess the quality of all selected studies. The scores of NOS ranges from zero (the lowest) to nine (the highest) stars, with more stars indicating a better methodological quality. The assessment was performed independently by two investigators and the inconformity was solved by a discussion, or consultation if necessary.

Statistical analysis

All the data management and analysis for this meta-analysis were performed with STATA 11.0 software (Stata corporation, College Station, TX). The odds ratios (ORs) with corresponding 95% confidence intervals (CIs) were used to estimate the association between CYP1A1 Ile462Val polymorphism and OSCC risk. In order to calculate heterogeneity of studies, the Chi-Square test was used and significance was set at P value less than 0.05 level 30. The inconsistency index I2 was calculated to assess the variation caused by heterogeneity. If the P value of the heterogeneity test was >0.10 and I2<25%, the fixed-effects model was performed to calculate the combined OR, which assumed the same homogeneity of effect size across all studies. If the P value of the heterogeneity test was <0.10 and I2>25%, it showed that the between-study heterogeneity was statistically significant, and a random effect model was used to estimate the pooled OR. The funnel plot was used to test the underlying publication bias, and the funnel plot asymmetry was estimated by Egger's linear regression 31. Sensitivity analyses were performed to identify the influence of the individual studies on the combined OR. In the analysis, we excluded each study to assess if stability between the remaining studies was reached.

Trial sequential analysis (TSA)

According to the Cochrane Handbook, meta-analyses and systematic reviews are considered to be the best available evidence if all eligible trials are included. However, 'the best available evidence' might not be always equal to 'sufficient evidence'. Based on this issue, we applied the TSA to estimate the power of the current conclusions 32-34. The methods of using TSA were based on the 'User manual for Trial Sequential Analysis (TSA)'.

Results

Characteristics of included studies

As shown in Figure , a total of 332 studies were retrieved by the literature search. Of which, 290 studies were excluded as being irrelevant to CYP1A1 Ile462Val or OSCC or genetic polymorphism or not human studies. And one potential eligible article was obtained by screening the references of reviews. After more detailed assessment for the left 43 potential eligible articles, one article was excluded for animal study; two papers were reviews. Following this, 12 papers were concerned about CYP1A1*2A or other genes, but not CYP1A1 Ile462Val. Besides, six papers were excluded because of failing to provide sufficient data. In addition, seven studies were excluded because the cases just were diagnosed oral cancer, were not confirmed the identification of OSCC. And two studies were excluded as they were concerned of prognosis and only contained the cases, lacking of controls. Finally, 13 papers 20-24, 35-42 were conformed to the inclusion criteria and they were eligible for the meta-analysis of CYP1A1 Ile462Val polymorphism.

We established a database concerning of the information extracted from each included paper. Summaries of these studies were presented in Table which included the first author, ethnicity, country, number and characteristics of cases and controls, and other relevant information. Of the 13 studies included in this meta-analysis, eight studies were performed in Asians, four in Caucasians, and one in mixed race. The number of cases and controls in the papers included in this meta-analysis for varied from 31 to 190, 60 to 241, respectively. The frequency of the CYP1A1 Ile462Val homozygous variant allele (Val/Val) in cases group varied from 0-13.4%, and the controls group were present in 0-6.2%.

The results of quality assessment

According to NOS system, of all observational studies are awarded a maximum of four stars in selection, two stars in comparability, and three stars in exposure. In this meta-analysis, the assessment results for the selected studies ranges from six to eight stars (Table ), indicating all of the selected studies are moderate-high qualities in methodology.

Test of heterogeneity and quantitative synthesis

We performed a heterogeneity analysis of the dominant model (Val/Val + Val/Ile vs. Ile/Ile), recessive model (Val/Val vs. Val/Ile + Ile/Ile), homozygous model (Val/Val vs. Ile/Ile), and allele comparing (Val vs. Ile). Owing to no heterogeneity was observed in recessive model (I2=22.0%, PQ- test=0.248), fixed-effect models were used to synthesized the data (Figure ). The overall results in recessive model suggested that individuals carrying homozygous variant (Val/Val) might have an increased risk of OSCC when compared with those who carried the homozygous Ile allele (OR=1.64, 95% CI=1.08-2.49).

As for the dominant model (I2=72.0%, PQ- test=0.000), homozygous model (I2=54.7%, PQ- test=0.024), random-effect models were used to synthesize the data, respectively. However, results in dominant model, homozygous model failed to obvious association of CYP1A1 Ile462Val with risks of OSCC. In order to further explore the observed heterogeneity in dominant model and additive model, subgroup analyses were performed by ethnicity, source of control, and HWE and the details were shown in Table .

As for allele comparing model, obvious heterogeneity were observed (I2 =66.4%; PQ- test=0.002), but the result had no statistical significance (P=0.308). The details of allele comparing model were also shown in the Table .

Publication bias

The Begg's funnel plot was used to evaluate the possible publication bias. Besides, the Begg's test and Egger's linear regression were used to the quantitative evaluation of the symmetry of the meta-analysis funnel plots and the results were as follows: (1) dominant model: Begg's Test, P=0.502 and Egger's linear regression test: t=0.49, P=0.634; (2) recessive model: Begg's Test, P=0.230 and Egger's linear regression test: t=0.26, P=0.803; (3) homozygous model: Begg's Test, P=0.230 and Egger's linear regression test: t=0.89, P=0.416; (4) allele comparing model: Begg's Test, P=0.371 and Egger's linear regression test: t=-1.42, P=0.192. The data of all four models indicated that there were not significant publication biases for all them.

Sensitivity analysis

In order to evaluate the stability of the results and reveal the influence of each study on the pooled ORs, sensitivity analysis were performed by excluding each case-control study respectively. When all of the four models were performed sensitivity analysis, the estimates in these four models changed between lower CI limits and upper CI limits, suggesting the results in this meta-analysis were stable.

Trial sequential analysis

Thirteen trials (3651 subjects) were used to investigate the association between CYP1A1 Ile462Val polymorphism and OSCC risk. Using the recessive model (including 10 trials with 2547 subjects) as an example, we performed the TSA and found that the required information size to demonstrate clear conclusions was 5174 subjects (Figure ). The cumulative z-curve does not cross the trial monitoring boundary before reaching the required power, which indicates that the cumulative evidence is insufficient and further trials are necessary (Figure ). The results of other models were not shown as the study methods were similar.

Discussion

Oral cancer is the cancer of mouth, including squamous cell carcinoma, adenocarcinoma, verrucous carcinoma, and so on. Different histopathologic types of cancers might have different genetic susceptibilities, for example, CYP1A1 Ile462Val polymorphism is a risk factor of squamous cell carcinoma of lung, but the associations vary in different histological types of lung cancer 17, 43. Therefore, it is more reasonable to assess the association of gene polymorphisms with risks of OSCC, oral adenocarcinoma and other types, separately.

To our knowledge, this is the first meta-analysis to assess the association between CYP1A1 Ile462Val polymorphism and OSCC risk. Although there were two previous meta-analyses 44, 45 concerning of CYP1A1 Ile462Val polymorphism and risks of oral cancer, the results did not involve in single histopathologic type and therefore they failed to be on behalf of the association of CYP1A1 Ile462Val polymorphism with OSCC risk. For instance, Singh, et al. 46 found no association between CYP1A1 Ile462Val polymorphism and oral cancer risk. However, the cases included in this research are OSCC and verrucous carcinoma. Thus, the results might fail to be on behalf of the association of OSCC and CYP1A1 Ile462Val polymorphism. Based on what mentioned above and for obtaining a powerful conclusion concerning about risks of OSCC and CYP1A1 Ile462Val polymorphism, we performed this systematical meta-analysis to evaluate it.

In the present meta-analysis, for overall data, the results of recessive model showed that the homozygous Val allele mutation (Val/Val) might play a critical role in the occurrence and development of OSCC (OR=1.64, 95% CI=1.08-2.49), which not only suggested that the same homogeneity of effect size was across all relevant studies but also were inconsistent with the previous meta-analyses 44, 45. As the results shown, the individuals who carried homozygous mutation (Val/Val) compared with those who carried wild genotype (Ile/Ile) and heterozygous variant (Val/Ile) might have an increased susceptibility of OSCC, implying that the OSCC might have different susceptibility from other histological types of oral cancer. As for the other three models (homozygous model, dominant model and allele comparing model), they failed to show obvious association between CYP1A1 Ile462Val polymorphism and OSCC risk, and the heterogeneities among studies were observed in them, respectively.

In order to explore the source of heterogeneity, we performed the subgroup analyses. After the subgroup analyses by ethnicity, source of control, HWE, the heterogeneities were not removed, indicating that other factors, such as age, gender, country, lifestyle, social status, smoking and alcohol habits might also yield to heterogeneities. The confused results may be resulted from few studies, few subjects, or other potential factors, and the more credible conclusions need the future studies containing large sample sizes and well-designed criteria to confirm.

Inevitably, several disadvantages should be noticed. First, the selected papers only involved in data of Asia, America, Europe, but without the related data about Africa and Australia. People come from different places may share different genotypes, thus, applicability of the results are limited. Second, in the current meta-analysis, only 13 articles with 1468 cases and 2183 controls were included. Therefore, the study might only have limited statistical power. Third, we performed a subgroup analysis by ethnicity, source of control, HWE, but the other factors, such as gender, age, were not performed for data limitations. Finally, heterogeneity was existed, which might weaken the reliability of conclusions. Based on the limitations mentioned above, the results should be considered with caution.

Overall, in spite of these limitations, the results of this analysis suggest that the homozygous variant might be a risk factor of OSCC. And future studies focusing on CYP1A1 Ile462Val polymorphism containing large sample sizes and well-designed criteria are necessary to make the conclusions more credible.

Table 1

Characteristics of the studies evaluating the association of CYP1A1 Ile 462 Val polymorphism with the susceptibility of OSCC

Year	First author	Ethnicity	Country	cases				controls				Control source	Matching	NOS	HWE
				Val/Val	Val/Ile	Ile/Ile	All	Val/Val	Val/Ile	Ile/Ile	All
1997	Park	Caucasian	USA	23*		108	131	10*		121	131	Healthy	Age, gender	8	Yes
1999	Katoh	Asian	Japan	6	36	50	92	8	50	89	147	Non-cancer	---	6	Yes
1999	Morita	Asian	Japan	2	6	23	31	6	54	104	164	Healthy	---	7	Yes
2000	Sato	Asian	Japan	19	55	68	142	6	46	90	142	Non-cancer	Age, gender	8	Yes
2002	Hahn	Caucasian	Germany	0	4	90	94	0	6	86	92	Healthy	---	6	Yes
2002	Kao	Asian	China	8	84	14	106	2	78	66	146	Non-cancer	---	6	Yes
2003	Gronau	Caucasian	Germany	0	18	55	73	1	35	100	136	Non-cancer	Age, smoking &alcohol habits	8	Yes
2006	Sugimura	Asian	Japan	8	33	81	112	15	88	138	241	Non-cancer	---	6	Yes
2006	Leichsenring	Mixed	Brazil	2	18	52	72	0	16	44	60	Healthy	Age, gender, smoking habits	8	Yes
2008	Buch	Caucasian	USA	28*		162	190	65*		345	410	Healthy	Age, gender	8	NA
2008	Sam	Asian	India	3	35	149	187	4	36	180	220	Non-cancer	Age, gender	8	Yes
2009	Amtha	Asian	Indonesia	36*		45	81	85*		77	162	Non-cancer	Age, gender	8	NA
2012	Balaji	Asian	India	3	34	120	157	5	26	101	132	Non-cancer	Age	6	Yes

NA: not available; HWE: Hardy-Weinberg Equilibrium; *: Val/Val + Val/Ile

Note: control source, Non-cancer means patients are in hospital without cancers.

Table 2

Meta-analysis of the association between CYP1A1 Ile 462 Val polymorphism and OSCC risk

Variables	Val/Val+Val/Ile vs. Ile/Ile					Val/Val vs. Val/Ile+Ile/Ile					Val/Val vs. Ile/Ile					Val vs. Ile
	OR(95% CI)	I2(%)	P (Q- test)	P		OR(95% CI)	I2 (%)	P(Q- test)	P		OR(95% CI)	I2 (%)	P(Q- test)	P		OR(95% CI)	I2 (%)	P(Q- test)	P
Total	1.17 (0.86-1.61)	72	<0.001	0.321		1.64 (1.08-2.49)	22.0	0.248	0.019		1.74 (0.84-3.60)	54.7	0.024	0.133		1.16 (0.88-1.52)	66.4	0.002	0.308
Ethnicity
Caucasian	1.12 (0.66-1.88)	50.4	0.109	0.676		0.61 (0.02-15.28)	---	---	0.766		0.60 (0.02-15.07)	---	---	0.758		1.22 (0.87-1.70)	0.0	0.652	0.537
Asian	1.21 (0.78-1.89)	80.9	<0.001	0.398		1.63 (1.07-2.50)	36.9	0.147	0.024		1.75 (0.79-3.92)	64.5	0.01	0.169		0.84 (0.49-1.45)	75.5	0.0	0.254
Mixed	1.06 (0.49-2.29)	---	---	0.887		4.29 (0.20-91.12)	---	---	0.35		4.24 (0.20-90.62)	---	---	0.355		1.17 (0.59-2.35)	---	---	0.654
Control Source
Healthy	1.04 (0.65-1.67)	46.3	0.114	0.876		2.38 (0.59-9.65)	38.6	0.62	0.225		1.97 (0.44-8.22)	0.0	0.554	0.386		0.91 (0.57-1.45)	0.0	0.603	0.682
Non-cancer	1.25 (0.82-1.91)	79.9	0.000	0.297		1.59 (1.03-2.46)	0.0	0.134	0.037		1.69 (0.72-4.00)	65.3	0.008	0.231		1.23 (0.89-1.70)	74.9	0.001	0.214
HWE
Yes	1.27 (0.87-1.83)	73.9	0.000			1.64 (1.08-2.49)	22.0	0.248	0.019		1.74 (0.84-3.60)	54.7	0.024	0.133		1.16 (0.88-1.52)	66.4	0.002	0.308
NA	0.83 (0.58-1.18)	0.0	0.521			---	---	---	---		---	---	---	---		---	---	---	---

I2: the variation in OR attributable to heterogeneity; P(Q-test)>0.05, heterogeneity was not statistically significant; P>0.05, no statistical significant