Genetic polymorphisms of CYP1A1 and risk of leukemia: a meta-analysis.

The associations between CYP1A1 polymorphisms and risk of leukemia have been studied extensively, but the results have been inconsistent. Therefore, in this study, we performed a meta-analysis to clarify associations of three CYP1A1 polymorphisms (T3801C, A2455G, and C4887A) with the risks of acute lymphoblastic leukemia (ALL), acute myeloid leukemia (AML), and chronic myeloid leukemia (CML). Medline, EMBASE, and China National Knowledge Infrastructure databases were searched to collect relevant studies published up to April 20, 2015. The extracted data were analyzed statistically, and pooled odds ratios with 95% confidence intervals were calculated to quantify the associations. Overall, 26 publications were included. Finally, T3801C was associated with an increased risk of AML in Asians under the dominant model. For A2455G, the risk of ALL was increased among Caucasians in the recessive model and the allele-contrast model; A2455G was also associated with an increased risk of CML among Caucasians under the recessive model, dominant model, and allele-contrast model. For C4887A, few of the included studies produced data. In conclusion, the results suggest that Asians carrying the T3801C C allele might have an increased risk of AML and that Caucasians with the A2455G GG genotype might have an increased risk of ALL. Further investigations are needed to confirm these associations.

Background

Leukemia is a malignant tumor of the hematopoietic system.1 There are four common types of leukemia: acute lymphoblastic leukemia (ALL), acute myeloid leukemia (AML), chronic lymphoblastic leukemia, and chronic myeloid leukemia (CML). Environmental exposures and genetic susceptibility play roles in the etiology of leukemia.2 For environmental exposures, benzene, ionizing radiation, and cytotoxic therapy are some of the proposed causes of acute leukemia; for genetic susceptibility, the most common type is the single nucleotide polymorphism (SNP) that consists of a variation at a single base pair. Depending on where it is located, SNPs can interfere with a gene’s function, affecting metabolic pathways and further affecting the course of the disease and its progress. SNPs in xenobiotic system, DNA repair system, and cell regulation have been identified as risk factors in childhood leukemia.3,4

CYP1A1 and leukemia

Cytochrome P450 (CYP) enzyme catalyzes the phase I metabolism reaction. The CYP1A1 gene (CYP1A1) is a member of the CYP family that participates in the metabolism of xenobiotics and endogenous compounds, particularly polycyclic aromatic hydrocarbons such as benzo[a]pyrene.5 Genetic polymorphisms have been reported for CYP1A1 that may alter the function of enzyme and thus influence the ability of enzymes to metabolize the chemical carcinogens and mutagens,6,7 which may influence how susceptible individuals are to contracting leukemia.8

SNPs in CYP1A1

The most commonly studied SNP in CYP1A1 is the T3801C polymorphism (also referred to as CYP1A1*2A, m1, or rs4646903). This SNP was localized to chromosome 15q229 and is characterized by a T to C mutation at nucleotide 3,801 in the 3′-flanking region of CYP1A1. Some studies have indicated that the T3801C polymorphism can alter the level of gene expression or the stability of messenger RNA (mRNA), resulting in the enzyme exhibiting highly inducible activity.10 Another common genetic mutant site is A2455G (also referred to as CYP1A1*2C, m2 allele, or rs1048943), which has been widely examined with regard to cancer susceptibility.11–13 This SNP is characterized by an A to G mutation at nucleotide 2,455 in the region of CYP1A1.14 The C4887A mutation (CYP1A1*4 allele) is also an important mutant site, which results in the replacement of Thr by Asn in codon 461, near the site of the A2455G mutation.15 Previous research has shown that this polymorphism can alter the level of gene expression or the stability of mRNA, resulting in the enzyme exhibiting highly inducible activity. In other words, certain variant genotypes of CYP1A1 that change the enzymatic activity appear to crucially influence the risk of leukemia.16

Aim of this study

Previous studies of the relationships between CYP1A1 polymorphisms and risk of leukemia have generated controversial results, and so whether CYP1A1 polymorphisms are a risk factor for leukemia remains uncertain. In this study, we performed a quantitative meta-analysis with the aim of generating more reliable results by applying an analysis with greater statistical power.

Materials and methods

Search strategy

We searched Medline (PubMed), EMBASE, and CNKI (China National Knowledge Infrastructure) databases for publications written in English or Chinese up to April 20, 2015. The search was based on the following Boolean combination of keywords: “CYP1A1” (“Cytochrome P450 1A1”) AND [“ALL” (“acute lymphocytic leukemia” OR “acute lymphoblastic leukemia”) OR “AML” (“acute myeloid leukemia” OR “acute myeloblastic leukemia” OR “acute myelocytic leukemia”) OR “CML” (“chronic myeloid leukemia” OR “chronic myelocytic leukemia” OR “chronic myeloblastic leukemia”) OR “leukemia” (“leucocythemia”)]. All searched studies were retrieved, and their bibliographies were checked for other relevant publications. Review articles and bibliographies of other relevant studies identified were searched manually to find additional eligible studies.

Inclusion criteria

Published studies were included if they fulfilled the following criteria:

Use of a case-control design.

Providing sufficient data on the distribution of CYP1A1 polymorphisms in leukemia and in controls or sufficient information for such data to be calculated.

Focusing on the associations between CYP1A1 polymorphisms and the risk of AML, ALL, and/or CML.

Data extraction

Two investigators independently carried out an eligibility evaluation and data abstraction for each potentially eligible study. Any disagreements were further discussed and resolved by consensus. The following data were abstracted from each of the included studies: name of the first author, year of publication, race distribution of the study population, genotype distributions of the cases and controls, and numbers of cases and controls. We also checked that the distribution of the genotypes in controls conformed with Hardy–Weinberg equilibrium (HWE).

Statistical analysis

The meta-analysis examined the overall associations under the recessive model (CC versus TC + TT for T3801C, GG versus AG + AA for A2455G, and CC versus AC + AA for C4887A), the dominant model (CC + TC versus TT for T3801C, GG + AG versus AA for A2455G, and CA + AA versus CC for C4887A), and the allele-contrast model (C allele versus T allele for T3801C, G allele versus A allele for A2455G, and A allele versus C for C4887A).

The raw data for genotype distribution were used to calculate the study-specific estimates of odds ratio (OR) and 95% confidence interval (CI). Review Manager (version 5.3) software (The Nordic Cochrane Centre, Copenhagen, Denmark) was used to implement the meta-analysis. The presence of heterogeneity was assessed using Cochran’s Q statistic and quantified using the I2 statistic, where a value above 50% indicates the presence of a very high degree of heterogeneity.17

ORs and 95% CIs were calculated using the Z-test to assess the associations between the CYP1A1 SNPs and risk of leukemia. The threshold for statistical significance was set at P<0.05. The fixed effects model was used when no significant clinical or statistical heterogeneity was present (I2≤50%),18 whereas the random effects model was used when there was substantial heterogeneity (I2>50%).

When substantial heterogeneity was present, sensitivity analysis was performed by excluding individual studies; outlying studies were identified and excluded, and the I2 estimates for these different sets of studies were examined.

Potential publication bias was estimated by constructing funnel plots. If most of the data appeared at the top of a funnel plot and were distributed approximately symmetrically, this was considered to indicate the absence of obvious publication bias and vice versa.19 Funnel plots were not constructed when there were too few analyzed studies.

Results

Overview of the characteristics of the included studies

The flow chart for the study selection process is depicted in Figure 1. We selected 217 articles identified by the Boolean combination of the 5 keywords listed in the “Search strategy” section. Checking for duplicates resulted in the removal of 118 articles. Of the remaining 99 articles, 45 did not report studies of the associations between CYP1A1 SNPs and leukemia, 6 did not focus on leukemia, 10 were review articles, and 11 did not provide enough data, and so all of these articles were also discarded. Moreover, the genotypes in controls did not conform with HWE in the study of Chen et al.20 Therefore, ultimately only 26 studies related to the relationships between CYP1A1 SNPs and the risk of leukemia remained in this meta-analysis, among which 13, 5, and 2 related to ALL, AML, and CML, respectively, 5 articles related to both ALL and AML, and 1 article related to all three disease types.

Figure 1

Flow chart of study selection.

Abbreviations: ALL, acute lymphoblastic leukemia; AML, acute myeloid leukemia; CML, chronic myeloid leukemia; HWE, Hardy–Weinberg equilibrium; SNP, single nucleotide polymorphism.

Study characteristics

The basic data of the 26 articles included in the meta-analysis were extracted and are listed in Table 1. In the study of Majumdar et al21 the distributions of the A2455G and C4887A genotypes in controls did not conform with HWE (although the T3801C genotypes did), and so this study was excluded when we analyzed the roles of the CYP1A1 A2455G and C4887A polymorphisms. Overall, the meta-analysis included the following studies:

Table 1

Study characteristics

Study	Ethnicity	Cases	SNP	Number of cases	Number of controls	PHWE in controls
Agha et al29	Caucasian	ALL	T3801C, C4887A	186	200	0.0577
Aydin-Sayitoglu et al30	Caucasian	ALL	T3801C	249	140	0.3057
		AML	T3801C	249	140	0.3057
Bolufer et al31	Caucasian	ALL	T3801C	443	454	0.151
		AML	T3801C	443	454	0.151
Bonaventure et al32	Caucasian	ALL	T3801C	493	549	0.2309
		AML	T3801C	493	549	0.2309
Krajinovic et al27	Caucasian	ALL	T3801C	177	304	0.3941
			C4887A			0.1355
Ouerhani et al33	Caucasian	ALL	T3801C	100	106	0.5731
Swinney et al12	Caucasian	ALL	T3801C, A2455G	258	646	0.1272
Yamaguti et al13	Caucasian	ALL	T3801C	99	99	0.1701
			A2455G			0.9485
Selvin49	Caucasian	ALL	A2455G	175	175	0.9347
Gao et al34	Asian	ALL	T3801C	78	112	0.1723
			A2455G			0.1357
Joseph et al35	Asian	ALL	T3801C	118	118	0.5939
			A2455G			0.141
Lee et al36	Asian	ALL	T3801C	164	164	0.4015
			A2455G			0.4514
Suneetha et al37	Asian	ALL	T3801C	92	150	0.806
Balta et al38	Mixed	ALL	T3801C	144	185	0.0907
		AML	T3801C	144	185	0.0907
Canalle et al39	Mixed	ALL	T3801C	113	221	0.083
Clavel et al40	Mixed	ALL	T3801C	219	105	0.1861
		AML	T3801C	219	105	0.1861
Gallegos-Arreola et al41	Mixed	ALL	T3801C	210	228	0.7728
Razmkhah et al43	Mixed	ALL	A2455G	85	94	0.4353
		AML	A2455G	85	94	0.9543
		CML	A2455G	85	94	0.9543
Gallegos-Arreola et al48	Mixed	ALL	A2455G	136	136	0.315
D’Alo et al44	Caucasian	AML	T3801C	193	273	0.1604
			A2455G			0.4279
			C4887A			0.7257
Jiang et al45	Asian	AML	T3801C	98	120	0.1062
Kim et al46	Asian	AML	T3801C	415	1,700	0.6805
			A2455G			0.0688
Majumdar et al21	Asian	AML	T3801C	110	126	0.8415
Yamaguti et al47	Mixed	AML	T3801C	133	133	0.1381
			A2455G			0.3451
Al-Achkar et al42	Caucasian	CML	A2455G	126	172	0.3253
Taspinar et al22	Caucasian	CML	A2455G	107	135	0.5889

Abbreviations: ALL, acute lymphoblastic leukemia; AML, acute myeloid leukemia; CML, chronic myeloid leukemia; PHWE, probability of adherence to the Hardy–Weinberg equilibrium; SNP, single nucleotide polymorphism.

Three studies involving 272 cases and 402 controls related to the role of the A2455G polymorphism in the risk of CML.

Four studies involving 846 cases and 2,201 controls related to the role of the A2455G polymorphism in the risk of AML.

Ten studies involving 1,330 cases and 3,688 controls related to the role of the T3801C polymorphism in the risk of AML.

One study involving 193 cases and 273 controls related to the role of the C4887A polymorphism in the risk of AML.

Eight studies involving 1,734 cases and 2,194 controls related to the role of the A2455G polymorphism in the risk of ALL.

Sixteen studies involving 2,299 cases and 3,209 controls related to the role of the T3801C polymorphism in the risk of ALL.

One study involving 177 cases and 295 controls related to the role of the C4887A polymorphism in the risk of ALL.

Meta-analysis results

The main results of the meta-analysis are listed in Table 2.

Table 2

ORs and heterogeneity results for CYP1A1 meta-analysis

SNP	Cancer	Ethnicity	n	Genetic model	I2	Contrast model	P, ORs (95% CIs)
T3801C	ALL	Caucasian	8	Recessive	26%	Fixed effects	0.28, 1.31 (0.81, 2.12)
				Dominant	61%	Random effects	0.10, 1.30 (0.95, 1.79)
				Allele contrast	32%	Fixed effects	0.05, 1.18 (1.00, 1.39)
		Asian	4	Recessive	63%	Random effects	0.73, 1.14 (0.54, 2.39)
				Dominant	68%	Random effects	0.06, 1.64 (0.98, 2.73)
				Allele contrast	80%	Random effects	0.19, 1.35 (0.86, 2.14)
		Mix	4	Recessive	76%	Random effects	0.21, 2.68 (0.57, 12.62)
				Dominant	82%	Random effects	0.32, 1.34 (0.75, 2.39)
				Allele contrast	92%	Random effects	0.34, 1.42 (0.69, 2.91)
A2455G	ALL	Caucasian	4	Recessive	0%	Fixed effects	0.002, 2.23 (1.36, 3.68)
				Dominant	39%	Fixed effects	0.12, 1.22 (0.95, 1.58)
				Allele contrast	41%	Fixed effects	0.010, 1.31 (1.07, 1.61)
		Asian	3	Recessive	11%	Fixed effects	0.14, 1.57 (0.87, 2.83)
				Dominant	72%	Random effects	0.29, 1.40 (0.76, 2.59)
				Allele contrast	0%	Fixed effects	0.89, 1.02 (0.80, 1.29)
		Mix	1	Recessive		Fixed effects
				Dominant		Fixed effects	0.94, 1.03 (0.45, 2.34)
				Allele contrast		Fixed effects	0.94, 1.03 (0.47, 2.26)
C4887A	ALL	Caucasian	1	Recessive		Fixed effects	0.61, 1.67 (0.23, 11.99)
				Dominant		Fixed effects	0.14, 0.57 (0.27, 1.21)
				Allele contrast		Fixed effects	0.22, 0.65 (0.33, 1.30)
T3801C	AML	Caucasian	4	Recessive	0%	Fixed effects	0.26, 1.74 (0.66, 4.55)
				Dominant	47%	Fixed effects	0.16, 0.82 (0.62, 1.08)
				Allele contrast	59%	Random effects	0.63, 0.90 (0.59, 1.37)
		Asian	3	Recessive	0%	Fixed effects	0.60, 0.93 (0.71, 1.22)
				Dominant	0%	Fixed effects	0.03, 1.26 (1.02, 1.55)
				Allele contrast	0%	Fixed effects	0.23, 1.09 (0.95, 1.26)
		Mix	3	Recessive	0%	Fixed effects	0.72, 1.19 (0.46, 3.05)
				Dominant	74%	Random effects	0.65, 1.24 (0.48, 3.22)
				Allele contrast	72%	Random effects	0.81, 1.11 (0.48, 2.55)
A2455G	AML	Caucasian	1	Recessive
				Dominant		Fixed effects	0.47, 0.78 (0.39, 1.53)
				Allele contrast		Fixed effects	0.48, 0.78 (0.40, 1.53)
		Asian	1	Recessive		Fixed effects	0.39, 0.79 (0.47, 1.35)
				Dominant		Fixed effects	0.42, 1.09 (0.88, 1.36)
				Allele contrast		Fixed effects	0.72, 1.03 (0.87, 1.23)
		Mix	2	Recessive	0%	Fixed effects	0.50, 1.64 (0.39, 6.92)
				Dominant	0%	Fixed effects	0.0009, 1.96 (1.32, 2.90)
				Allele contrast	0%	Fixed effects	0.002, 1.72 (1.22, 2.43)
C4887A	AML	Caucasian	1	Recessive		Fixed effects	0.39, 2.85 (0.26, 31.63)
				Dominant		Fixed effects	0.005, 2.16 (1.27, 3.69)
				Allele contrast		Fixed effects	0.004, 2.08 (1.26, 3.44)
A2455G	CML	Caucasian	2	Recessive	17%	Fixed effects	0.004, 19.61 (2.62, 146.76)
				Dominant	69%	Random effects	0.0010, 3.93 (1.74, 8.87)
				Allele contrast	1%	Fixed effects	P,0.00001, 3.92 (2.61, 5.90)
		Mix	1	Recessive		Fixed effects	0.89, 0.80 (0.03, 20.00)
				Dominant		Random effects	0.79, 0.88 (0.34, 2.29)
				Allele contrast		Fixed effects	0.70, 0.84 (0.34, 2.07)

Abbreviations: ALL, acute lymphoblastic leukemia; AML, acute myeloid leukemia; CI, confidence interval; CML, chronic myeloid leukemia; OR, odds ratio; SNP, single nucleotide polymorphism.

CYP1A1 T3801C was not associated with ALL risk in Caucasians or Asians

In the meta-analysis of 16 studies that involved data from 2,299 cases and 3,209 controls and investigated the effect of the T3801C polymorphism on the risk of ALL, a high degree of heterogeneity was observed under the recessive model (I2=76%, P<0.00001), dominant model (I2=68%, P<0.00001), and allele-contrast model (I2=80%, P<0.00001). A sensitivity analysis was conducted since I2 exceeded 50% for all three models. However, removing any particular article from the analysis did not reduce the degree of heterogeneity, and so the random effects model was used. This identified an association between the T3801C SNP and the risk of ALL under the dominant model (P=0.006, OR =1.39, 95% CI =1.10–1.76) and the allele-contrast model (P=0.03, OR =1.31, 95% CI =1.02–1.67), but not under the recessive model (P=0.27, OR =1.43, 95% CI =0.76–2.70). Moreover, the funnel plot did not reveal any obvious publication bias.

Different results can be obtained when including different studies due to population differences in the T3801C mutation frequency or linkage disequilibrium block. Subgroup analysis according to race was conducted, which revealed no marked associations in the three genetic comparisons among Asians (recessive model: P=0.73, OR =1.14, 95% CI =0.54–2.39; dominant model P=0.06, OR =1.64, 95% CI =0.98–2.73; and allele-contrast model: P=0.19, OR =1.35, 95% CI =0.86–2.14) and Caucasians (recessive model: P=0.28, OR =1.31, 95% CI =0.81–2.12; dominant model: P=0.10, OR =1.30, 95% CI =0.95–1.79; and allele-contrast model: P=0.05, OR =1.18, 95% CI =1.00–1.39) (Figure 2).

Figure 2

Meta-analysis of the association between CYP1A1 T3801C gene polymorphism and ALL risk under three models: (A) recessive, (B) dominant, and (C) allele contrast.

Abbreviations: ALL, acute lymphoblastic leukemia; CI, confidence interval; M–H, Mantel–Haenszel type; OR, odds ratio; SE, standard error.

CYP1A1 A2455G was associated with ALL risk in Caucasians

In the meta-analysis of eight studies that involved the data from 1,734 patients and 2,194 controls and investigated the effect of the A2455G polymorphism on the risk of ALL, the forest plot revealed no heterogeneity under the recessive model (I2=0%, P=0.49), dominant model (I2=44%, P=0.09), and allele-contrast model (I2=20%, P=0.27). The fixed effects model was therefore used for these three contrast models, which revealed an association between the A2455G polymorphism and the risk of ALL under the recessive model (P=0.0007, OR =1.93, 95% CI =1.32–2.82), the dominant model (P=0.02, OR =1.25, 95% CI =1.03–1.52), and the allele-contrast model (P=0.04, OR =1.17, 95% CI =1.01–1.36). The funnel plot did not reveal any obvious publication bias.

In subgroup analysis according to race, the risk of leukemia was higher in Caucasians under the recessive model (P=0.002, OR =2.23, 95% CI =1.36–3.68) and the allele-contrast model (P=0.010, OR =1.31, 95% CI =1.07–1.61), but not under the dominant model (P=0.12, OR =1.22, 95% CI =0.95–1.58) (Figure 3).

Figure 3

Meta-analysis of the association between CYP1A1 A2455G gene polymorphism and ALL risk under three models: (A) recessive, (B) dominant, and (C) allele contrast.

Abbreviations: CML, chronic myeloid leukemia; CI, confidence interval; M–H, Mantel–Haenszel type; OR, odds ratio; ALL, acute lymphoblastic leukemia.

CYP1A1 C4887A was not associated with ALL risk in Caucasians

In the meta-analysis of the single study that involved 177 cases and 295 controls (all of whom were Caucasians) and investigated the effect of the C4887A polymorphism on the risk of ALL, no association was found between the CYP1A1 C4887A SNP and the risk of ALL under any contrast model when using the fixed effects model.

CYP1A1 T3801C was associated with AML risk in Asians

In the meta-analysis of ten studies that involved 1,330 cases and 3,688 controls and investigated the effect of the T3801C polymorphism on the risk of AML, a high degree of heterogeneity was observed under the recessive model (I2=62%, P=0.010), dominant model (I2=65%, P=0.002), and allele-contrast model (I2=72%, P=0.0002). A sensitivity analysis was therefore conducted, which revealed that the study of Majumdar et al21 was an outlier. Removing their data from the analysis reduced the degree of heterogeneity. Under the allele-contrast model, the degree of heterogeneity was still high, and so a random effects model was applied. No association was found between the T3801C SNP and the risk of AML under the recessive model (P=0.91, OR =0.98, 95% CI =0.76–1.27), the dominant model (P=0.07, OR =1.15, 95% CI =0.99–1.34), or the allele-contrast model (P=0.71, OR =1.05, 95% CI =0.82–1.34). The funnel plot did not reveal any obvious publication bias.

However, in subgroup analysis according to race, the risk was increased among Asians under the dominant model (P=0.03, OR =1.26, 95% CI =1.02–1.55), but not under the recessive model (P=0.60, OR =0.93, 95% CI =0.71–1.22) or the allele-contrast model (P=0.23, OR =1.09, 95% CI =0.95–1.26) (Figure 4).

Figure 4

Meta-analysis of the association between CYP1A1 T3801C gene polymorphism and AML risk under three models: (A) recessive, (B) dominant, and (C) allele contrast.

Abbreviations: AML, acute myeloid leukemia; CI, confidence interval; M–H, Mantel–Haenszel type; OR, odds ratio; SE, standard error.

CYP1A1 A2455G was not associated with AML risk in Caucasians or Asians

The meta-analysis of four studies that involved 846 patients and 2,201 controls – one study involved Caucasians, one study involved Asians, and two studies involved mixed races – found no association between the A2455G SNP and the risk of AML under any contrast model.

CYP1A1 C4887A was associated with AML risk in Caucasians

In the meta-analysis of the single study that involved 193 cases and 273 controls (all of whom were Caucasians), when using the fixed effects model there was a significant association between the C4887A SNP and the risk of AML under the dominant model (P=0.005, OR =2.16, 95% CI =1.27–3.69) and the allele-contrast model (P=0.004, OR =2.08, 95% CI =1.26–3.44). However, this relationship needs to be confirmed in further studies.

CYP1A1 A2455G was associated with CML risk in Caucasians

In the meta-analysis of the three studies that involved 272 patients and 402 controls and investigated the effect of the A2455G polymorphism on the risk of CML, the forest plot revealed no heterogeneity under the recessive model (I2=43%, P=0.17), and so the fixed effects model was applied. There was a high degree of heterogeneity under the dominant model (I2=81%, P=0.006) and the allele-contrast model (I2=80%, P=0.006), and so a sensitivity analysis was conducted. However, removing any particular study from the analysis did not reduce the degree of heterogeneity, and so the random effects model was used for these two contrast models. An association was found between the A2455G SNP and the risk of CML under the recessive model (P=0.001, OR =9.91, 95% CI =2.43–40.36) and the allele-contrast model (P=0.04, OR =2.56, 95% CI =1.05–6.23).

In subgroup analysis according to race, the risk of leukemia was increased in Caucasians under the recessive model (P=0.004, OR =19.61, 95% CI =2.62–146.76), dominant model (P=0.0010, OR =3.93, 95% CI =1.74–8.87), and allele-contrast model (P<0.00001, OR =3.92, 95% CI =2.61–5.90) (Figure 5).

Figure 5

Meta-analysis of the association between CYP1A1 A2455G gene polymorphism and CML risk under three models: (A) recessive, (B) dominant, and (C) allele contrast.

Abbreviations: ALL, acute lymphoblastic leukemia; CI, confidence interval; M–H, Mantel–Haenszel type; OR, odds ratio; SE, standard error; CML, chronic myeloid leukemia.

Discussion

This study identified 26 articles describing research into the associations between the CYP1A1 T3801C, CYP1A1 A2455G, and CYP1A1 C4887A polymorphisms and the risks of ALL, AML, and CML. The conclusions drawn using the data obtained in these studies were not consistent. This situation prompted the use of meta-analysis in this study – which is a common method applied in systematic reviews – to integrate these contradictions, and to obtain more accurate results. All of the finally included studies were of high quality and had a rigorous scientific design, accurate data, and clear results. The main conclusion is that the presence of the CYP1A1 T3801C polymorphism increased the risks of AML among Asians, while the CYP1A1 A2455G polymorphism increased the risks of ALL among Caucasians.

Few studies have investigated the CYP1A1 polymorphism as a possible risk factor for CML.22 The present meta-analysis found an association between the A2455G SNP and the risk of CML among Caucasians. However, this result was based on a relatively small number of cases where the statistical power was weak, and so, this conclusion should be treated with caution. The accumulation of disadvantageous environmental factors might be required to develop chronic disease, and so, investigating environment – genetic interactions – would be meaningful for this kind of chronic disease. Although, in this study, we endeavored to extract relevant information regarding smoking and alcohol consumption from the primary literature, insufficient data were obtained. Further investigations concerning the interactions between smoking, alcohol consumption, and gene variations with CML are required.

This meta-analysis included studying the associations between three SNPs (T3801C, A2455G, and C4887A) and the risk of AML. The results showed that the CYP1A1 T3801C polymorphism might not be significantly correlated with the risk of AML. Possible racial variations in presentation, treatment patterns, and survival with respect to AML might exist,23 which made subgroup analysis according to race necessary. We found a significant increased risk among Asians under the dominant model, but not under the recessive or allele-contrast model. Compared with the meta-analysis of Zhuo et al24 who also researched the association between T3801C and the risk of AML, the dominant model for Asians utilized a fixed effects in the present meta-analysis, which makes our conclusions more convincing. For A2455G, the subgroup analyses found no association among either Asians or Caucasians with the risk of AML. The meta-analysis of Zhuo et al25 also investigated the association between A2455G and the risk of AML, but our conclusions are more convincing due to the inclusion of a larger number of articles. In addition, we believe that the article of Zhuo et al contained some data errors.25 The research data of Majumdar et al21 defined the allele *2B as an m1 (T3801C) and m2 (A2455G) combined locus, and so the data for *2B cannot represent the data for A2455G. Moreover, the raw data for A2455G genotypes in controls did not conform with HWE, and so our meta-analysis excluded that group of data. Zhuo et al also included the combined mutation data, which undermines the credibility of their meta-analysis.25 Our meta-analysis included only one study of the C4887A SNP and revealed an association between this SNP and the risk of AML under the dominant and allele-contrast models. However, the sample was too small to reliably characterize the intrinsic link between the CYP1A1 C4887A polymorphism and risk of AML. Further investigations are needed to confirm this conclusion.

Our meta-analysis also studied the associations between three SNPs (T3801C, A2455G, and C4887A) and the risk of ALL. The results showed that the CYP1A1 T3801C polymorphism might have a strong correlation with an increased risk of ALL. However, in the subgroup analyses stratified by race, the data did not reveal such an association among either Asians or Caucasians when they were analyzed separately. Our analysis included more studies than that of Zhuo et al26 which also researched the association between T3801C and the risk of ALL, and so our conclusions are more plausible. The results for the CYP1A1 A2455G polymorphism indicated that this is significantly correlated with the risk of ALL among Caucasians. Again, our analysis included more studies than that of Zhuo et al25 which also researched the association between A2455G and ALL. In addition, we believe that the article of Zhuo et al contained some data errors.25 Krajinovic et al27 found that CYP1A1*2B represents both m1 (T3801C) and m2 (A2455G), and so, the data for *2B cannot represent the data for A2455G alone. Moreover, Zhuo et al included combined mutation data, which would produce incorrect results.25 More importantly, our study only extracted data for a single race from the study of Swinney et al12 which makes our race subgroup analyses more credible. Only one study related to the CYP1A1 C4887A polymorphism was included in our meta-analysis. The results showed that this polymorphism might not be correlated with the risk of ALL. Also, the subgroup analyses found no association among either Asians or Caucasians when they were analyzed separately. The sample was too small to reliably characterize the intrinsic link between the CYP1A1 C4887A polymorphism and the risk of ALL, and so further studies are needed to clarify this.

There are two implications of this study that are particularly significant: 1) the CYP1A1 A2455G polymorphism can be used as a screening index for the risk of ALL and CML, while the CYP1A1 T3801C polymorphism can be used as a screening index for the risk of AML; and 2) genetic screening can be used to identify individuals with the susceptible genotype who could take early defense measures to reduce their risk of leukemia.

Three limitations of this meta-analysis should be considered when interpreting its findings. First, our findings only related to the relationships between specific genes and diseases, whereas leukemia is influenced by multiple factors, and so, future studies need to consider age, race, family history, and environmental and lifestyle factors when attempting to draw more accurate conclusions. Second, a relatively small number of eligible studies were included in this meta-analysis, which could limit the generalizability of the results. This also restricted the statistical power of the subgroup analyses. As is well known, the inclusion of studies with very small samples may lead to overestimations of the associations.28 Third, only published studies were included in the present meta-analysis. A certain degree of publication bias is always present due to the possibility of nonsignificant or negative findings not being published.

In conclusion, this meta-analysis suggests that Asians who carry the T3801C C allele might have an increased risk of AML, while Caucasians with the A2455G GG genotype might have an increased risk of ALL. In addition, we found an association between the CYP1A1 A2455G polymorphism and risk of CML among Caucasians. Future large studies of these sites are recommended in order to confirm the findings of this study and also to eventually identify the underlying mechanisms.