Citation: Noel Pabalan, Neetu Singh, Eloisa Singian, Caio Parente Barbosa, Bianca Bianco, Hamdi Jarjanazi. Associations of CYP1A1 gene polymorphisms and risk of breast cancer in Indian women: a meta-analysis[J]. AIMS Genetics, 2015, 2(4): 250-262. doi: 10.3934/genet.2015.4.250
[1] | Tú Nguyen-Dumont, Jenna Stewart, Ingrid Winship, Melissa C. Southey . Rare genetic variants: making the connection with breast cancer susceptibility. AIMS Genetics, 2015, 2(4): 281-292. doi: 10.3934/genet.2015.4.281 |
[2] | Mohammad Hashemi, Fatemeh Bizhani, Hiva Danesh, Behzad Narouie, Mehdi Sotoudeh, Mohammad Hadi Radfar, Mehdi Honarkar Ramezani, Gholamreza Bahari, Mohsen Taheri, Saeid Ghavami . MiR-608 rs4919510 C>G polymorphism increased the risk of bladder cancer in an Iranian population. AIMS Genetics, 2016, 3(4): 212-218. doi: 10.3934/genet.2016.4.212 |
[3] | Ke-Sheng Wang, Xuefeng Liu, Daniel Owusu, Yue Pan, Changchun Xie . Polymorphisms in the ANKS1B gene are associated with cancer, obesity and type 2 diabetes. AIMS Genetics, 2015, 2(3): 192-203. doi: 10.3934/genet.2015.3.192 |
[4] | Khalid Eltahir Khalid . Vitamin D receptor gene polymorphisms in Sudanese children with type 1 diabetes. AIMS Genetics, 2016, 3(3): 167-176. doi: 10.3934/genet.2016.3.167 |
[5] | Alena V Malakhova, Olga I Rudko, Vladimir V Sobolev, Artemii V Tretiakov, Elena A Naumova, Zarema G Kokaeva, Julia E Azimova, Eugene A Klimov . PDE4B gene polymorphism in Russian patients with panic disorder. AIMS Genetics, 2019, 6(3): 55-63. doi: 10.3934/genet.2019.3.55 |
[6] | Agnès Collet, Julien Tarabeux, Elodie Girard, Catherine Dubois DEnghien, Lisa Golmard, Vivien Deshaies, Alban Lermine, Anthony Laugé, Virginie Moncoutier, Cédrick Lefol, Florence Copigny, Catherine Dehainault, Henrique Tenreiro, Christophe Guy, Khadija Abidallah, Catherine Barbaroux, Etienne Rouleau, Nicolas Servant, Antoine De Pauw, Dominique Stoppa-Lyonnet, Claude Houdayer . Pros and cons of HaloPlex enrichment in cancer predisposition genetic diagnosis. AIMS Genetics, 2015, 2(4): 263-280. doi: 10.3934/genet.2015.4.263 |
[7] | Carina Ladeira, Susana Viegas, Mário Pádua, Elisabete Carolino, Manuel C. Gomes, Miguel Brito . Relation between DNA damage measured by comet assay and OGG1 Ser326Cys polymorphism in antineoplastic drugs biomonitoring. AIMS Genetics, 2015, 2(3): 204-218. doi: 10.3934/genet.2015.3.204 |
[8] | Michael J. Murray . Drosophila models of metastasis. AIMS Genetics, 2015, 2(1): 25-53. doi: 10.3934/genet.2015.1.25 |
[9] | Tin Tin Su . Non-autonomous consequences of cell death and other perks of being metazoan. AIMS Genetics, 2015, 2(1): 54-69. doi: 10.3934/genet.2015.1.54 |
[10] | Jennifer K. Redig, Gameil T. Fouad, Darcie Babcock, Benjamin Reshey, Eleanor Feingold, Roger H. Reeves, Cheryl L. Maslen . Allelic Interaction between CRELD1 and VEGFA in the Pathogenesis of Cardiac Atrioventricular Septal Defects. AIMS Genetics, 2014, 1(1): 1-19. doi: 10.3934/genet.2014.1.1 |
Breast cancer is a leading cause of death among women in India [1]. The multifactorial nature of this disease is marked by the synergy of interaction among various environmental and genetic factors [2,3]. One genetic factor is the CYP1A1 protein, encoded by CYP1A1 gene, an interesting candidate that might influence susceptibility to breast cancer risk for reasons mentioned by Chen et al. [4]. CYP1A1 gene polymorphisms have been extensively studied, especially in relation to breast cancer susceptibility [4,5,6,7,8,9,10,11,12,13]. The CYP1A1 gene, located at 15q22-q24, comprises seven exons and six introns spanning 5,810 base pairs [14]. Of the four common polymorphisms of the CYP1A1 gene identified, we examine three in the Indian population: M1 (rs4646903: T3801C, giving rise to a MspI restriction site in the 3’-noncoding region [15], M2 (rs1048943: A2455G), resulting in an amino acid change at codon 462 from isoleucine to valine within the heme-binding domain of exon 7 [16]; and M4 (rs1799814: C2453A), resulting in an amino acid substitution of threonine with asparagine at codon 461 [17]. Of the three polymorphisms, M4 is very rare.
Though the functional significance of variant CYP1A1 genotypes is unclear [4], studies of CYP1A1 in cultured human lymphocytes showed significantly elevated levels of inducible enzyme activity among M2 genotypes compared with the wild-type genotype [18]. M2 alleles appear to be associated with CYP1A1 inducibility at the level of transcription followed by 3-fold elevation in aryl hydrocarbon hydroxylase enzyme activity [19]. The M1 allele was also reported to be more readily inducible than the CYP1A1 wild-type allele [20,21]. A number of meta-analyses have appeared addressing associations of CYP1A1 polymorphisms with breast cancer. To our knowledge, no meta-analysis, including the four we chose to compare our study with, has focused on a population of a single country in Asia. Yet three of them examined CYP1A1 associations with breast cancer in Asian populations using subgroup analysis. Yao et al. [6] and Chen et al. [4] did this among Asians in general and east-Asians, respectively. Still these populations were a heterogeneous mix of different Asian nationalities. The Sergentanis et al. [5] examination of the Chinese population was more homogeneous. Still none of these meta-analyses have addressed CYP1A1 associations with breast cancer in the Indian population. Given that breast cancer is an overriding epidemiological issue in India [9] and that reports do not necessarily agree with each other, we conducted a meta-analysis of published case-control studies to examine the effect of CYP1A1 genetic polymorphisms on breast cancer risk among Indian women.
Two strategies were used in searching MEDLINE using PubMed and Embase. In the first, search terms were broad such as “DNA polymorphism”, “breast cancer” and “xenobiotic”. In a separate second search, search terms were more specific: “CYP1A1” and “cytochrome P450 1A1”. The broad and specific search terms were combined with “breast cancer” and “Indian” to find association studies as of 14 April 2015. Studies were eligible if they had genotypic data with a case-control design restricted to English. Figure 1 outlines our study selection process in a flowchart following PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analyses) guidelines [22]. The two search strategies yielded a total of 31 citations, 23 of which were screened for not fulfilling our inclusion criteria. Excluding a duplicate article, full texts of the remaining seven were obtained which were judged as eligible; as such, they were included in the meta-analysis. We manually searched the reference sections of each of these studies for additional articles. Syamala et al. [7] provided sporadic and familial data thus were treated as separate studies. In all, this meta-analysis had a total number of eight studies from seven publications [7,8,9,10,11,12,13].
Two investigators independently extracted data and reached consensus on all the items. The following information was obtained from each publication: first author’s name, published year, city and region in India, state of controls, matching criteria, sample source, genotype data, number of cases and controls. We also calculated frequencies of the variant allele as well as deviations of the controls from the Hardy-Weinberg equilibrium (HWE).
Quality of the studies was evaluated with use of the Newcastle-Ottawa Scale (NOS) [23]. Each study was assessed based on three broad perspectives: selection, comparability, and exposure with a score ranging from 0 to 9. Study quality scores for high, medium and poor were ≥ 7, 4-6 and < 4, respectively.
Risks (odds ratios [ORs] and 95% confidenceintervals) of breast cancer in three CYP1A1 polymorphisms (M1, M2, and M4) were estimated for each study. Pooled ORs were calculated for the following genetic models using variant (var) and wild-type (wt) genotypes: (i) homozygous (var-var, genotypes compared with the wt-wt), (ii) recessive (var-var vs. wt-var + wt-wt) (iii) dominant (var-var + wt-var vs. wt-wt) and (iv) co-dominant genetic model (var vs. wt). To compare effects on the same baseline, we used raw data for genotype frequencies to calculate study-specific estimates of the OR. Pooled ORs were obtained using either the fixed [24] (in absence of heterogeneity) or random [25] (in its presence) effects models. Heterogeneity between studies was estimated using the I2-based Q test [26] and quantified with the I2 statistic which measures degree of inconsistency among studies [27]. Given the low power of the heterogeneity test [28], significance threshold was set at p = 0.10. Sources of heterogeneity were examined with subgroup analysis [27]. Sensitivity analysis, which involved omitting one study at a time and recalculating the pooled OR, was also used to test for robustness of the summary effects. Availability of analyzable data from the South and North Indian subpopulations from the included studies allowed this type of geographical subgrouping. Supported by precedent literature [29,30], the probability of differential risk associations between the South and North Indian subpopulations warranted testing for the presence of interactions; wherein the generated p values were subjected to the Bonferroni correction. Data were analyzed using Review Manager 5.3 (Copenhagen: Nordic Cochrane Centre, Cochrane Collaboration, 2014) and SigmaStat 3.5 (Systat Software, San Jose, CA). Significance was set at a p-value of ≤ 0.05 throughout except in heterogeneity estimation. Publication bias was not investigated because of the low sensitivity of the qualitative and quantitative tests when the number of studies is lower than ten [31].
Table 1 shows features of the included studies, published between 2005 and 2011. Five studies had Southern subjects [7,8,11,12,13] compared to two from the North [9,10]. There were more hospital-based [7,10,11,12] than population-based [8,9,13] studies. Two studies made no mention of matching [11,13]. Methodological quality of the studies indicates that they were moderate with most studies having a score of 6-7.
First author | Year | Indian city | Indian region | State of controls | Source of controls | Matching criteria | NOS |
Singh N | 2007 | Lucknow | North | normal | Hospital | menopause, diet | 6 |
Singh V | 2007 | Lucknow | North | healthy | Population | residence | 7 |
Surekha | 2009 | Hyderabad | South | healthy | Population | age | 7 |
Chacko | 2005 | Kerala | South | healthy | Hospital | age, menopause | 6 |
Syamala* | 2010 | Thiruvanan | South | cancer-free | Hospital | residence, gender | 6 |
Naushad | 2011 | Hyderabad | South | healthy | Hospital | no mention | 4 |
Kiruthiga | 2011 | Madurai | South | no mention | Population | no mention | 5 |
* Separate familial and sporadic data were considered as two studies; NOS: Newcastle-Ottawa Score. |
Table S1 summarizes genotypic features of the component studies in the meta-analysis. In the M1 polymorphism association studies (1151 cases/1527 controls), populations of South and North India were examined in four [7,11,12,13] and two [9,10] papers, respectively. These Indian geographical regions were likewise examined in four [8,9,10,11] and two [12,13] papers on the M2 polymorphism (1,039 cases/1022 controls). The two studies [9,10] and a single one [11] on the M4 polymorphism (631 cases/611 controls) had North and South Indian subjects, respectively. Control subjects deviated from the HWE in the following: (i) two studies [9,12] in M1; and (ii) one study for M2 [8] and two for M4 [10,11]. For both M1 and M2, the mean frequency of the minor allele was greater from the Northern populations, though this was significant only for M1 (M1: North x = 0.32 ± 0.08 standard deviation (SD) versus South x = 0.18 ± 0.05 SD, p = 0.03; M2: North x = 0.18 ± 0.06 SD versus South x = 0.16 ± 0.11 SD, p = 0.81).
Table S2 shows genotypic frequencies in the three CYP1A1 polymorphisms in premenopausal (M1: 279 cases/309 controls; M2: 277 cases/309 controls) and postmenopausal (M1: 230 cases/251 controls; M2: 232 cases/251 controls) women. The M1 and M2 polymorphisms had three studies each [9,10,12]. In the M4 polymorphism, premenopausal (168 cases/197 controls) and postmenopausal (118 cases/139 controls) subjects were all North Indians on account of two studies [9,10], one of which deviated from the HWE [10]. In the geography and menopausal status analyses, the probability of differential risk associations between North and South India as well as premenopausal and postmenopausal women warranted testing for presence of interactions.
In Table 2, all three polymorphisms showed overall increased risk across all comparisons. Significance was observed in the M1 (OR 1.61-3.43, p = 0.007-0.04) and M4 (OR 2.02-3.92, p = 0.02-0.04) polymorphisms. The significant homozygous and recessive M4 effects (OR 3.81-3.92, p = 0.02) were obtained in the absence of heterogeneity (pheterogeneity = 0.79-0.84, I2 = 0%). All other effects in all the polymorphisms were heterogeneous, particularly high for M1 (pheterogeneity < 0.00001, I2 = 85-96%). Figure 2 exemplifies the significant effects of M1 under conditions of high heterogeneity.
Test of association | Test of heterogeneity | ||||||
Genetic model | N | OR | 95% CI | pA | pB | I2 (%) | |
M1 | Homozygous | 7 | 2.31 | 0.98-5.46 | 0.06 | <0.0001 | 88 |
Recessive | 7 | 2.02 | 0.95-4.29 | 0.07 | <0.0001 | 85 | |
Dominant | 7 | 3.43 | 1.41-8.35 | 0.007 | <0.0001 | 96 | |
Codominant | 7 | 1.61 | 1.03-2.52 | 0.04 | <0.0001 | 92 | |
M2 | Homozygous | 6 | 2.53 | 0.65-9.95 | 0.18 | 0.003 | 72 |
Recessive | 6 | 2.39 | 0.68-8.42 | 0.18 | 0.009 | 68 | |
Dominant | 6 | 1.50 | 0.83-2.74 | 0.18 | <0.0001 | 89 | |
Codominant | 6 | 1.51 | 0.89-2.58 | 0.13 | <0.0001 | 90 | |
M4 | Homozygous | 2 | 3.92 | 1.25-12.33 | 0.02 | 0.79 | 0 |
Recessive | 2 | 3.81 | 1.22-11.96 | 0.02 | 0.84 | 0 | |
Dominant | 2 | 1.96 | 0.87-4.43 | 0.11 | 0.003 | 83 | |
Codominant | 2 | 2.02 | 1.03-3.95 | 0.04 | 0.01 | 78 | |
N: number of studies; OR: odds ratio; CI: confidence interval; pA and pB: p values for association and heterogeneity with significance set at < 0.05 and < 0.10, respectively; Values in bold indicate significant tests of association; R: Random-effects model; F: Fixed-effects model. |
Figure 2 Forest plot of association of the CYP1A1 M1 polymorphism on breast cancer in the overall analysis in the dominant genetic model of the overall analysis. Diamond denotes the pooled odds ratio (OR). Squares indicate the OR in each study, with square sizes directly proportional to the weight contribution (%) of each study. Horizontal lines represent 95% confidence intervals (CI). The Z test for overall effect indicates significance if p < 0.05. The chi-square test indicates heterogeneity (p < 0.10) warranting use of the random-effects model.
Table 3 shows subgroup associationsstratified by geography in the M1 polymorphism. Mainly, the North Indians were protected (OR 0.64-0.77, p = 0.03-0.55), but not their Southern counterparts, who were at significant risk (OR 2.20-4.34, p < 0.0001). The tests of interaction highlight these contrasting effects, with significance generated in all (p = 0.00002-0.02) but the recessive model. Table 4 compares increased risk effects between the two Indian populations in the M2 polymorphism which did not materially differ given the non-significance of the p values of interaction (p > 1).
South India | North India | ||||||||||
n = 5 | n = 2 | ||||||||||
Test of association | Test of heterogeneity | Test of association | Test of heterogeneity | ||||||||
OR | 95% CI | pA | pB | I2 | OR | 95% CI | pA | pB | I2 | pi | |
Homozygous | 4.34 | 3.07-6.12 | <0.0001 | 0.20 | 33 | 0.64 | 0.22-1.93 | 0.43 | 0.03 | 78 | 0.019 |
Recessive | 3.55 | 2.54-4.95 | <0.0001 | 0.50 | 0 | 0.69 | 0.20-2.34 | 0.55 | 0.01 | 83 | 0.23 |
Dominant | 3.73 | 0.88-15.71 | 0.007 | <0.0001 | 97 | 2.81 | 2.16-3.65 | <0.0001 | 0.63 | 0 | >1 |
Codominant | 1.17 | 1.00-1.37 | 0.05 | <0.0001 | 86 | 2.72 | 2.23-3.31 | <0.0001 | 0.65 | 0 | >1 |
Premenopausal | Postmenopausal | ||||||||||
n = 3 | n = 3 | ||||||||||
Test of association | Test of heterogeneity | Test of association | Test of heterogeneity | ||||||||
OR | 95% CI | pA | pB | I2 | OR | 95% CI | pA | pB | I2 | pi | |
Homozygous | 1.87 | 0.63-5.53 | 0.26 | 0.04 | 68 | 1.35 | 0.53-3.41 | 0.53 | 0.09 | 58 | >1 |
Recessive | 1.67 | 0.65-4.30 | 0.29 | 0.07 | 62 | 1.54 | 0.92-2.60 | 0.10 | 0.19 | 40 | >1 |
Dominant | 1.30 | 0.92-1.81 | 0.13 | 0.23 | 31 | 1.03 | 0.56-1.86 | 0.94 | 0.09 | 58 | >1 |
Codominant | 1.31 | 0.87-1.97 | 0.20 | 0.08 | 60 | 1.08 | 0.67-1.75 | 0.74 | 0.06 | 65 | >1 |
OR: odds ratio; CI: confidence interval; pA and pB: p values for association and heterogeneity with significance set at < 0.05 and < 0.10, respectively; pi: p values for tests of interaction after Bonferroni correction; Values in bold are significant in the tests of association and tests of interaction |
South India | North India | ||||||||||
n = 5 | n = 2 | ||||||||||
Test of association | Test of heterogeneity | Test of association | Test of heterogeneity | ||||||||
OR | 95% CI | pA | pB | I2 | OR | 95% CI | pA | pB | I2 | pi | |
Homozygous | 3.72 | 0.55-25.44 | 0.18 | 0.006 | 81 | 1.36 | 0.06-31.06 | 0.85 | 0.009 | 79 | >1 |
Recessive | 3.34 | 0.60-18.59 | 0.17 | 0.02 | 76 | 1.26 | 0.06-24.67 | 0.88 | 0.01 | 77 | >1 |
Dominant | 2.37 | 0.63-8.92 | 0.20 | <0.0001 | 92 | 1.17 | 0.54-2.55 | 0.69 | <0.0001 | 89 | >1 |
Codominant | 2.81 | 0.75-10.51 | 0.12 | <0.0001 | 94 | 1.10 | 0.58-2.09 | 0.76 | 0.0001 | 89 | >1 |
Premenopausal | Postmenopausal | ||||||||||
n = 3 | n = 3 | ||||||||||
Test of association | Test of heterogeneity | Test of association | Test of heterogeneity | ||||||||
OR | 95% CI | pA | pB | I2 | OR | 95% CI | pA | pB | I2 | pi | |
Homozygous | 1.69 | 0.20-14.35 | 0.63 | 0.03 | 73 | 1.79 | 0.89-3.59 | 0.10 | 0.30 | 17 | >1 |
Recessive | 1.74 | 0.25-12.11 | 0.58 | 0.04 | 69 | 1.62 | 0.85-3.11 | 0.14 | 0.35 | 6 | >1 |
Dominant | 0.90 | 0.42-1.94 | 0.79 | 0.009 | 79 | 1.43 | 0.96-2.11 | 0.08 | 0.23 | 33 | >1 |
Codominant | 0.99 | 0.46-2.14 | 0.98 | 0.0007 | 86 | 1.36 | 1.01-1.83 | 0.04 | 0.11 | 54 | >1 |
OR: odds ratio; CI: confidence interval; pA and pB: p values for association and heterogeneity with significance set at < 0.05 and < 0.10, respectively; pi: p values for tests of interaction after Bonferroni correction; Values in bold are significant in the tests of association. |
In all three polymorphisms, differential effects between premenopausal and postmenopausal women were not significant (Tests for Interaction: p > 1) where both subgroups generally exhibited increased risk effects (Tables 3-4). These effects were significant in the codominant model of M2 and M4 among postmenopausal women (Tables 4 and 5), suggesting this subgroup to be at particular risk. Among premenopausal women, M2 effects indicated codominant null and protective dominant effects suggesting variability of the results (Table 4).
Premenopausal | Postmenopausal | ||||||||||
n = 2 | n = 2 | ||||||||||
OR | 95% CI | pA | pB | I2 | OR | 95% CI | pA | pB | I2 | pi | |
Homozygous | 3.70 | 0.73-18.63 | 0.11 | 0.67 | 0 | 4.77 | 0.93-24.41 | 0.06 | 0.96 | 0 | >1 |
Recessive | 3.63 | 0.72-18.30 | 0.12 | 0.74 | 0 | 4.59 | 0.90-23.44 | 0.07 | 0.99 | 0 | >1 |
Dominant | 1.16 | 0.70-1.90 | 0.56 | 0.13 | 56 | 1.61 | 0.84-3.11 | 0.15 | 0.30 | 6 | >1 |
Codominant | 1.26 | 0.81-1.95 | 0.30 | 0.23 | 31 | 1.83 | 1.04-3.24 | 0.04 | 0.23 | 30 | >1 |
OR: odds ratio; CI: confidence interval; pA and pB: p values for association and heterogeneity with significance set at < 0.05 and < 0.10, respectively; pi: p values for tests of interaction after Bonferroni correction; Values in bold are significant in the tests of association. |
Direction of effects of the pooled ORs in M1 and M2 polymorphisms remained unchanged at increased risk with sensitivity treatment indicating that they were not materially altered suggesting robustness of our findings.
With a combined sample size of 5981 (2821 cases and 3160 controls) in three CYP1A1 polymorphisms, this meta-analysis points to evidence of overall associations between the M1, M2 and M4 polymorphisms and breast cancer among Indian women. However, this association was limited by between study heterogeneity. Its sources were examined with subgroup analysis. Stratifying the studies according to geography revealed the North India subgroup to be more heterogeneous than their Southern counterparts.
Three meta-analyses have addressed associations of the CYP1A1 polymorphisms in breast cancer with mixed populations [4,5,6]. The following polymorphisms were investigated by Chen et al. [4] (M1 and M2), Yao et al. [6] (M1) and Sergentanis et al. [5] (M1, M2, M3 and M4).
Our meta-analysis differs from these three in that we examined the M1, M2 and M4 polymorphisms in the Indian population only enabling discernment of differences in effects between North and South Indians. The Indian population is the second largest in the world with diverse social, cultural, linguistic and biological features spread across a wide topography [32].
A fourth meta-analysis that focused on the M1 polymorphism addressed its effects in breast cancer among South Indians as a subgroup of four comparisons [33]. These effects indicated significant increased risk for up to 4-fold in all the genetic models they used, which was similar to our findings from five studies. The difference between that meta-analysis and ours was that our statistical comparison between the North and South Indian populations resulted in significant differential effects, protecting the former and placing the latter at risk.
The combination of population stratification leading to multiple comparisons and borderline significance of some associations open the possibility of false-positive findings [34]. Avoiding this problem warrants that the subpopulation be defined in terms of marker-allele frequencies and factors that influence disease such as geography. In this meta-analysis, we addressed population stratification in terms of differential allele frequencies between two distinct regions of India, which intensified significance.
Still, our study has been limited by the following: (i) given the multiplicity of comparisons for different genetic models, geography and menopausal status subgroups, and the unavoidable flexibility of choosing and defining the correlates, associations may have been detected by chance alone; (ii) low sample sizes of the component studies which translate to weak statistical power. But then, the point of meta-analysis is to combine studies resulting in increase of statistical power; (iii) that most of the studies were hospital-based may preclude extrapolation of our findings to the general population; (iv) heterogeneity of most of our findings, (v) deviation of some studies from the HWE which may have biased summary outputs and point to methodological weaknesses, such as biased selection of subjects, genotyping errors population stratification [35]. However, this was addressed with sensitivity analysis. (vi) Our study was based on unadjusted data, as the genotype information stratified for the main confounding variables was unavailable in the original papers. In addition, addressing of the confounding factors varied across the studies.
Yet, despite these weaknesses, the strengths of our meta-analysis include: (i) statistical homogeneity of the M4 homozygous and recessive effects imply combinability of the component studies; (ii) consistency of overall findings across the genetic models; (iii) because we confined our meta-analysis to a single country, we were able to delineate differential summary effects between populations of North and South India; (iv) controls were uniformly defined (cancer-free, healthy), minimizing non-differential misclassification bias.
The findings we report here highlight two points. The first is that individual epidemiological studies may lack statistical power or may have other interfering factors which prevent the unmasking of overall associations. Likewise, geographically local associations with specific alleles may exist which are due to linkage disequilibrium with other disease modifying alleles and are not reproduced in other populations. The second point is a caution with regard to drawing conclusions from individual studies. We report results that were not apparent from an examination of individual studies. We find that meta-analysis is a useful tool to examine broad trends in CYP1A1 associations with breast cancer in Indian populations, and may avoid possible misleading conclusions based on only single-population studies.
The pathways of carcinogen metabolism are complex, mediated by the activities of multiple genes in concert with CYP1A1 such as GSTM1 [12]. It is conceivable that breast cancer risk related to any one locus will be small because gene-gene as well as gene-environment interactions are likely to operate. Individual studies in our meta-analysis have addressed gene-gene [7,12] and gene-environment [8,10,13]interaction and provided data which we excluded for concern with multiple testing issues.
To our knowledge, this is the firstmeta-analysis that details the association of CYP1A1 polymorphisms with breast cancer among Indian women. This focus on one country allowed us to examine differences with use of geography and menopausal subgroups. Because our meta-analysis summary effects are relevant only to this demography, we provide an epidemiological profile of Indian women’s susceptibility to breast cancer on account of the CYP1A1 polymorphisms. Additional well-designed studies based on sample sizes commensurate with detection of small genotypic risks should allow more definitive conclusions about the association of CYP1A1 polymorphisms and breast cancer among Indian women.
The authors declare no conflict of interest.
First author | Year | Indian Region | Case | Control | Case | Control | Case | Control | maf | HWE | |
M1 (1151/1527)* | TT | TC | CC | ||||||||
1 | Singh N | 2007 | North | 94 | 94 | 35 | 53 | 17 | 15 | 0.26 | 0.07 |
2 | Singh V | 2007 | North | 70 | 93 | 60 | 80 | 12 | 43 | 0.38 | 0.001 |
3 | Chacko | 2005 | South | 62 | 90 | 39 | 18 | 11 | 4 | 0.12 | 0.02 |
4 | Syamala s | 2010 | South | 93 | 243 | 87 | 106 | 39 | 18 | 0.19 | 0.22 |
5 | Syamala f | 2010 | South | 55 | 243 | 56 | 106 | 29 | 18 | 0.19 | 0.22 |
6 | Naushad | 2011 | South | 168 | 140 | 125 | 98 | 49 | 15 | 0.25 | 0.69 |
7 | Kiruthiga | 2011 | South | 30 | 36 | 17 | 12 | 3 | 2 | 0.16 | 0.45 |
M2 (1039/1022)* | AA | AG | GG | ||||||||
1 | Singh N | 2007 | North | 98 | 119 | 36 | 41 | 11 | 2 | 0.14 | 0.46 |
2 | Singh V | 2007 | North | 104 | 119 | 38 | 68 | 0 | 9 | 0.22 | 0.86 |
3 | Surekha | 2009 | South | 78 | 122 | 169 | 127 | 1 | 0 | 0.26 | <0.001 |
4 | Chacko | 2005 | South | 67 | 90 | 34 | 19 | 11 | 3 | 0.11 | 0.13 |
5 | Naushad | 2011 | South | 221 | 144 | 108 | 99 | 13 | 10 | 0.24 | 0.16 |
6 | Kiruthiga | 2011 | South | 34 | 48 | 4 | 2 | 12 | 0 | 0.02 | 0.89 |
M4 (631/611)* | CC | CA | AA | ||||||||
1 | Singh N | 2007 | North | 129 | 145 | 12 | 15 | 6 | 2 | 0.06 | 0.04 |
2 | Singh V | 2007 | North | 93 | 43 | 43 | 52 | 6 | 2 | 0.14 | 0.24 |
3 | Naushad | 2011 | South | 274 | 240 | 68 | 13 | 0 | 0 | 0.03 | 0.03 |
* CYP1A1 polymorphism (cases / controls); s: sporadic; f: familial; maf: minor allele frequency; HWE: Hardy-Weinberg Equilibrium p values. |
Premenopausal | Postmenopausal | |||||||||||||||||
First author | Year | Indian Region | Case | Control | Case | Control | Case | Control | maf | HWE | Case | Control | Case | Control | Case | Control | maf | HWE |
M1 | TT | TC | CC | TT | TC | CC | ||||||||||||
N (cases/controls) | (279/309) | (230/251) | ||||||||||||||||
Singh N | 2007 | North | 40 | 57 | 17 | 21 | 5 | 3 | 0.17 | 0.55 | 50 | 30 | 18 | 21 | 13 | 8 | 0.31 | 0.18 |
Singh V | 2007 | North | 53 | 58 | 43 | 35 | 9 | 13 | 0.31 | 0.35 | 17 | 35 | 17 | 35 | 3 | 10 | 0.34 | 0.79 |
Chacko | 2005 | South | 29 | 44 | 58 | 58 | 25 | 0 | 0.35 | 0.14 | 33 | 46 | 54 | 54 | 25 | 12 | 0.35 | 0.51 |
M2 | AA | AG | GG | AA | AG | GG | ||||||||||||
N (cases/controls) | (277/309) | (232/251) | ||||||||||||||||
Singh N | 2007 | North | 39 | 51 | 16 | 29 | 5 | 1 | 0.19 | 0.16 | 56 | 49 | 23 | 9 | 4 | 1 | 0.09 | 0.45 |
Singh V | 2007 | North | 80 | 70 | 25 | 41 | 0 | 5 | 0.22 | 0.74 | 24 | 49 | 13 | 27 | 0 | 4 | 0.22 | 0.91 |
Chacko | 2005 | South | 31 | 45 | 57 | 58 | 23 | 9 | 0.34 | 0.10 | 36 | 45 | 57 | 58 | 23 | 9 | 0.34 | 0.10 |
M4 | CC | CA | AA | CC | CA | AA | ||||||||||||
N (cases/controls) | (168/197) | (118/139) | ||||||||||||||||
Singh N | 2007 | North | 57 | 70 | 3 | 10 | 2 | 1 | 0.07 | 0.37 | 66 | 54 | 9 | 4 | 6 | 1 | 0.05 | 0.02 |
Singh V | 2007 | North | 67 | 83 | 34 | 32 | 4 | 1 | 0.15 | 0.27 | 36 | 59 | 9 | 20 | 2 | 1 | 0.14 | 0.63 |
maf: minor allele frequency; HWE: Hardy-Weinberg Equilibrium p values. |
[1] |
Dikshit R, Gupta PC, Ramasundarahettige C, et al. (2012) Cancer mortality in India: a nationally representative survey. Lancet 379: 1807-1816. doi: 10.1016/S0140-6736(12)60358-4
![]() |
[2] |
Vargo-Gogola T, Rosen JM (2007) Modelling breast cancer: one size does not fit all. Nat Rev Cancer 7: 659-672. doi: 10.1038/nrc2193
![]() |
[3] |
Nickels S, Truong T, Hein R, et al. (2013) Evidence of gene-environment interactions between common breast cancer susceptibility loci and established environmental risk factors. PLoS Genet 9: e1003284. doi: 10.1371/journal.pgen.1003284
![]() |
[4] |
Chen C, Huang Y, Li Y, et al. (2007) Cytochrome P450 1A1 (CYP1A1) T3801C and A2455G polymorphisms in breast cancer risk: a meta-analysis. J Hum Genet 52: 423-435. doi: 10.1007/s10038-007-0131-8
![]() |
[5] |
Sergentanis TN, Economopoulos KP (2010) Four polymorphisms in cytochrome P450 1A1 (CYP1A1) gene and breast cancer risk: a meta-analysis. Breast Cancer Res Treat 122: 459-469. doi: 10.1007/s10549-009-0694-5
![]() |
[6] |
Yao L, Yu X, Yu L (2010) Lack of significant association between CYP1A1 T3801C polymorphism and breast cancer risk: a meta-analysis involving 25,087 subjects. Breast Cancer Res Treat 122: 503-507. doi: 10.1007/s10549-009-0717-2
![]() |
[7] |
Syamala VS, Syamala V, Sheeja VR, et al. (2010) Possible risk modification by polymorphisms of estrogen metabolizing genes in familial breast cancer susceptibility in an Indian population. Cancer Invest 28: 304-311. doi: 10.3109/07357900902744494
![]() |
[8] |
Surekha D, Sailaja K, Rao DN, et al. (2009) Association of CYP1A1*2 polymorphisms with breast cancer risk: a case control study. Indian J Med Sci 63: 13-20. doi: 10.4103/0019-5359.49077
![]() |
[9] |
Singh V, Rastogi N, Sinha A, et al. (2007) A study on the association of cytochrome-P450 1A1 polymorphism and breast cancer risk in north Indian women. Breast Cancer Res Treat 101: 73-81. doi: 10.1007/s10549-006-9264-2
![]() |
[10] |
Singh N, Mitra AK, Garg VK, et al. (2007) Association of CYP1A1 polymorphisms with breast cancer in North Indian women. Oncol Res 16: 587-597. doi: 10.3727/000000007783629972
![]() |
[11] |
Naushad SM, Reddy CA, Rupasree Y, et al. (2011) Cross-talk between one-carbon metabolism and xenobiotic metabolism: implications on oxidative DNA damage and susceptibility to breast cancer. Cell Biochem Biophys 61: 715-723. doi: 10.1007/s12013-011-9245-x
![]() |
[12] |
Chacko P, Joseph T, Mathew BS, et al. (2005) Role of xenobiotic metabolizing gene polymorphisms in breast cancer susceptibility and treatment outcome. Mutat Res 581: 153-163. doi: 10.1016/j.mrgentox.2004.11.018
![]() |
[13] | Kiruthiga PV, Kannan MR, Saraswathi C, et al. (2011) CYP1A1 gene polymorphisms: lack of association with breast cancer susceptibility in the southern region (Madurai) of India. Asian Pac J Cancer Prev 12: 2133-2138. |
[14] |
Masson LF, Sharp L, Cotton SC, et al. (2005) Cytochrome P-450 1A1 gene polymorphisms and risk of breast cancer: a HuGE review. Am J Epidemiol 161: 901-915. doi: 10.1093/aje/kwi121
![]() |
[15] |
Kawajiri K, Nakachi K, Imai K, et al. (1990) Identification of genetically high risk individuals to lung cancer by DNA polymorphisms of the cytochrome P450IA1 gene. FEBS Lett 263: 131-133. doi: 10.1016/0014-5793(90)80721-T
![]() |
[16] | Hayashi SI, Watanabe J, Nakachi K, et al. (1991) PCR detection of an A/G polymorphism within exon 7 of the CYP1A1 gene. Nucleic Acids Res 19: 4797. |
[17] | Cascorbi I, Brockmoller J, Roots I (1996) A C4887A polymorphism in exon 7 of human CYP1A1: population frequency, mutation linkages, and impact on lung cancer susceptibility. Cancer Res 56: 4965-4969. |
[18] | Li Y, Millikan RC, Bell DA, et al. (2005) Polychlorinated biphenyls, cytochrome P450 1A1 (CYP1A1) polymorphisms, and breast cancer risk among African American women and white women in North Carolina: a population-based case-control study. Breast Cancer Res 7: R12-18. |
[19] |
Crofts F, Taioli E, Trachman J, et al. (1994) Functional significance of different human CYP1A1 genotypes. Carcinogenesis 15: 2961-2963. doi: 10.1093/carcin/15.12.2961
![]() |
[20] |
Kiyohara C, Hirohata T, Inutsuka S (1996) The relationship between aryl hydrocarbon hydroxylase and polymorphisms of the CYP1A1 gene. Jpn J Cancer Res 87: 18-24. doi: 10.1111/j.1349-7006.1996.tb00194.x
![]() |
[21] |
Li Y, Millikan RC, Bell DA, et al. (2004) Cigarette smoking, cytochrome P4501A1 polymorphisms, and breast cancer among African-American and white women. Breast Cancer Res 6: R460-473. doi: 10.1186/bcr814
![]() |
[22] |
Moher D, Liberati A, Tetzlaff J, et al. (2009) Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. Ann Intern Med 151: 264-269, W264. doi: 10.7326/0003-4819-151-4-200908180-00135
![]() |
[23] | Wells S PJ, Welch V. The newcastle-ottawa scale (NOS) for assessing the quality of nonrandomised studies in meta-analyses. Ottawa Health Research Institute, 2011. Available from: www.ohri.ca/programs/clinical_epidemiology/oxford/asp. |
[24] | Mantel N, Haenszel W (1959) Statistical aspects of the analysis of data from retrospective studies of disease. J Natl Cancer Inst 22: 719-748. |
[25] |
DerSimonian R, Laird N (1986) Meta-analysis in clinical trials. Control Clin Trials 7: 177-188. doi: 10.1016/0197-2456(86)90046-2
![]() |
[26] |
Lau J, Ioannidis JP, Schmid CH (1997) Quantitative synthesis in systematic reviews. Ann Intern Med 127: 820-826. doi: 10.7326/0003-4819-127-9-199711010-00008
![]() |
[27] |
Higgins JP, Thompson SG (2002) Quantifying heterogeneity in a meta-analysis. Stat Med 21: 1539-1558. doi: 10.1002/sim.1186
![]() |
[28] |
Higgins JP, Thompson SG, Deeks JJ, et al. (2003) Measuring inconsistency in meta-analyses. Bmj 327: 557-560. doi: 10.1136/bmj.327.7414.557
![]() |
[29] | Gautham M, Shyamprasad KM, Singh R, et al. (2014) Informal rural healthcare providers in North and South India. Health Policy Plan 29 Suppl 1: i20-29. |
[30] |
Ravindran RD, Vashist P, Gupta SK, et al. (2011) Prevalence and risk factors for vitamin C deficiency in north and south India: a two centre population based study in people aged 60 years and over. PLoS One 6: e28588. doi: 10.1371/journal.pone.0028588
![]() |
[31] |
Ioannidis JP, Trikalinos TA (2007) The appropriateness of asymmetry tests for publication bias in meta-analyses: a large survey. Cmaj 176: 1091-1096. doi: 10.1503/cmaj.060410
![]() |
[32] | Gadgil MaG, R. (1992) The fissure land: An ecological history of India; Press OU, editor. New Delhi: Oxford University Press. |
[33] |
He XF, Wei W, Liu ZZ, et al. (2014) Association between the CYP1A1 T3801C polymorphism and risk of cancer: evidence from 268 case-control studies. Gene 534: 324-344. doi: 10.1016/j.gene.2013.10.025
![]() |
[34] |
Wacholder S, Chanock S, Garcia-Closas M, et al. (2004) Assessing the probability that a positive report is false: an approach for molecular epidemiology studies. J Natl Cancer Inst 96: 434-442. doi: 10.1093/jnci/djh075
![]() |
[35] |
Thakkinstian A, McElduff P, D'Este C, et al. (2005) A method for meta-analysis of molecular association studies. Stat Med 24: 1291-1306. doi: 10.1002/sim.2010
![]() |
1. | Salih Q Ibrahem, Hussien Q Ahmed, Khalida M Amin, Genetic Variations in Cytochrome P450 1A1 and 1B1 Genes in a Cohort of Patients from Iraq Diagnosed with Breast Cancer, 2021, 15, 1178-2234, 10.1177/11782234211050727 | |
2. | Tajamul Hussain, Salman Alrokayan, Upadhyay Upasna, Manickam Pavithrakumari, Jaganathan Jayapriya, Vijay Kumar Kutala, Shaik Mohammad Naushad, Meta-analysis of genetic polymorphisms in xenobiotic metabolizing enzymes and their association with breast cancer risk, 2018, 97, 0022-1333, 523, 10.1007/s12041-018-0946-4 |
First author | Year | Indian city | Indian region | State of controls | Source of controls | Matching criteria | NOS |
Singh N | 2007 | Lucknow | North | normal | Hospital | menopause, diet | 6 |
Singh V | 2007 | Lucknow | North | healthy | Population | residence | 7 |
Surekha | 2009 | Hyderabad | South | healthy | Population | age | 7 |
Chacko | 2005 | Kerala | South | healthy | Hospital | age, menopause | 6 |
Syamala* | 2010 | Thiruvanan | South | cancer-free | Hospital | residence, gender | 6 |
Naushad | 2011 | Hyderabad | South | healthy | Hospital | no mention | 4 |
Kiruthiga | 2011 | Madurai | South | no mention | Population | no mention | 5 |
* Separate familial and sporadic data were considered as two studies; NOS: Newcastle-Ottawa Score. |
Test of association | Test of heterogeneity | ||||||
Genetic model | N | OR | 95% CI | pA | pB | I2 (%) | |
M1 | Homozygous | 7 | 2.31 | 0.98-5.46 | 0.06 | <0.0001 | 88 |
Recessive | 7 | 2.02 | 0.95-4.29 | 0.07 | <0.0001 | 85 | |
Dominant | 7 | 3.43 | 1.41-8.35 | 0.007 | <0.0001 | 96 | |
Codominant | 7 | 1.61 | 1.03-2.52 | 0.04 | <0.0001 | 92 | |
M2 | Homozygous | 6 | 2.53 | 0.65-9.95 | 0.18 | 0.003 | 72 |
Recessive | 6 | 2.39 | 0.68-8.42 | 0.18 | 0.009 | 68 | |
Dominant | 6 | 1.50 | 0.83-2.74 | 0.18 | <0.0001 | 89 | |
Codominant | 6 | 1.51 | 0.89-2.58 | 0.13 | <0.0001 | 90 | |
M4 | Homozygous | 2 | 3.92 | 1.25-12.33 | 0.02 | 0.79 | 0 |
Recessive | 2 | 3.81 | 1.22-11.96 | 0.02 | 0.84 | 0 | |
Dominant | 2 | 1.96 | 0.87-4.43 | 0.11 | 0.003 | 83 | |
Codominant | 2 | 2.02 | 1.03-3.95 | 0.04 | 0.01 | 78 | |
N: number of studies; OR: odds ratio; CI: confidence interval; pA and pB: p values for association and heterogeneity with significance set at < 0.05 and < 0.10, respectively; Values in bold indicate significant tests of association; R: Random-effects model; F: Fixed-effects model. |
South India | North India | ||||||||||
n = 5 | n = 2 | ||||||||||
Test of association | Test of heterogeneity | Test of association | Test of heterogeneity | ||||||||
OR | 95% CI | pA | pB | I2 | OR | 95% CI | pA | pB | I2 | pi | |
Homozygous | 4.34 | 3.07-6.12 | <0.0001 | 0.20 | 33 | 0.64 | 0.22-1.93 | 0.43 | 0.03 | 78 | 0.019 |
Recessive | 3.55 | 2.54-4.95 | <0.0001 | 0.50 | 0 | 0.69 | 0.20-2.34 | 0.55 | 0.01 | 83 | 0.23 |
Dominant | 3.73 | 0.88-15.71 | 0.007 | <0.0001 | 97 | 2.81 | 2.16-3.65 | <0.0001 | 0.63 | 0 | >1 |
Codominant | 1.17 | 1.00-1.37 | 0.05 | <0.0001 | 86 | 2.72 | 2.23-3.31 | <0.0001 | 0.65 | 0 | >1 |
Premenopausal | Postmenopausal | ||||||||||
n = 3 | n = 3 | ||||||||||
Test of association | Test of heterogeneity | Test of association | Test of heterogeneity | ||||||||
OR | 95% CI | pA | pB | I2 | OR | 95% CI | pA | pB | I2 | pi | |
Homozygous | 1.87 | 0.63-5.53 | 0.26 | 0.04 | 68 | 1.35 | 0.53-3.41 | 0.53 | 0.09 | 58 | >1 |
Recessive | 1.67 | 0.65-4.30 | 0.29 | 0.07 | 62 | 1.54 | 0.92-2.60 | 0.10 | 0.19 | 40 | >1 |
Dominant | 1.30 | 0.92-1.81 | 0.13 | 0.23 | 31 | 1.03 | 0.56-1.86 | 0.94 | 0.09 | 58 | >1 |
Codominant | 1.31 | 0.87-1.97 | 0.20 | 0.08 | 60 | 1.08 | 0.67-1.75 | 0.74 | 0.06 | 65 | >1 |
OR: odds ratio; CI: confidence interval; pA and pB: p values for association and heterogeneity with significance set at < 0.05 and < 0.10, respectively; pi: p values for tests of interaction after Bonferroni correction; Values in bold are significant in the tests of association and tests of interaction |
South India | North India | ||||||||||
n = 5 | n = 2 | ||||||||||
Test of association | Test of heterogeneity | Test of association | Test of heterogeneity | ||||||||
OR | 95% CI | pA | pB | I2 | OR | 95% CI | pA | pB | I2 | pi | |
Homozygous | 3.72 | 0.55-25.44 | 0.18 | 0.006 | 81 | 1.36 | 0.06-31.06 | 0.85 | 0.009 | 79 | >1 |
Recessive | 3.34 | 0.60-18.59 | 0.17 | 0.02 | 76 | 1.26 | 0.06-24.67 | 0.88 | 0.01 | 77 | >1 |
Dominant | 2.37 | 0.63-8.92 | 0.20 | <0.0001 | 92 | 1.17 | 0.54-2.55 | 0.69 | <0.0001 | 89 | >1 |
Codominant | 2.81 | 0.75-10.51 | 0.12 | <0.0001 | 94 | 1.10 | 0.58-2.09 | 0.76 | 0.0001 | 89 | >1 |
Premenopausal | Postmenopausal | ||||||||||
n = 3 | n = 3 | ||||||||||
Test of association | Test of heterogeneity | Test of association | Test of heterogeneity | ||||||||
OR | 95% CI | pA | pB | I2 | OR | 95% CI | pA | pB | I2 | pi | |
Homozygous | 1.69 | 0.20-14.35 | 0.63 | 0.03 | 73 | 1.79 | 0.89-3.59 | 0.10 | 0.30 | 17 | >1 |
Recessive | 1.74 | 0.25-12.11 | 0.58 | 0.04 | 69 | 1.62 | 0.85-3.11 | 0.14 | 0.35 | 6 | >1 |
Dominant | 0.90 | 0.42-1.94 | 0.79 | 0.009 | 79 | 1.43 | 0.96-2.11 | 0.08 | 0.23 | 33 | >1 |
Codominant | 0.99 | 0.46-2.14 | 0.98 | 0.0007 | 86 | 1.36 | 1.01-1.83 | 0.04 | 0.11 | 54 | >1 |
OR: odds ratio; CI: confidence interval; pA and pB: p values for association and heterogeneity with significance set at < 0.05 and < 0.10, respectively; pi: p values for tests of interaction after Bonferroni correction; Values in bold are significant in the tests of association. |
Premenopausal | Postmenopausal | ||||||||||
n = 2 | n = 2 | ||||||||||
OR | 95% CI | pA | pB | I2 | OR | 95% CI | pA | pB | I2 | pi | |
Homozygous | 3.70 | 0.73-18.63 | 0.11 | 0.67 | 0 | 4.77 | 0.93-24.41 | 0.06 | 0.96 | 0 | >1 |
Recessive | 3.63 | 0.72-18.30 | 0.12 | 0.74 | 0 | 4.59 | 0.90-23.44 | 0.07 | 0.99 | 0 | >1 |
Dominant | 1.16 | 0.70-1.90 | 0.56 | 0.13 | 56 | 1.61 | 0.84-3.11 | 0.15 | 0.30 | 6 | >1 |
Codominant | 1.26 | 0.81-1.95 | 0.30 | 0.23 | 31 | 1.83 | 1.04-3.24 | 0.04 | 0.23 | 30 | >1 |
OR: odds ratio; CI: confidence interval; pA and pB: p values for association and heterogeneity with significance set at < 0.05 and < 0.10, respectively; pi: p values for tests of interaction after Bonferroni correction; Values in bold are significant in the tests of association. |
First author | Year | Indian Region | Case | Control | Case | Control | Case | Control | maf | HWE | |
M1 (1151/1527)* | TT | TC | CC | ||||||||
1 | Singh N | 2007 | North | 94 | 94 | 35 | 53 | 17 | 15 | 0.26 | 0.07 |
2 | Singh V | 2007 | North | 70 | 93 | 60 | 80 | 12 | 43 | 0.38 | 0.001 |
3 | Chacko | 2005 | South | 62 | 90 | 39 | 18 | 11 | 4 | 0.12 | 0.02 |
4 | Syamala s | 2010 | South | 93 | 243 | 87 | 106 | 39 | 18 | 0.19 | 0.22 |
5 | Syamala f | 2010 | South | 55 | 243 | 56 | 106 | 29 | 18 | 0.19 | 0.22 |
6 | Naushad | 2011 | South | 168 | 140 | 125 | 98 | 49 | 15 | 0.25 | 0.69 |
7 | Kiruthiga | 2011 | South | 30 | 36 | 17 | 12 | 3 | 2 | 0.16 | 0.45 |
M2 (1039/1022)* | AA | AG | GG | ||||||||
1 | Singh N | 2007 | North | 98 | 119 | 36 | 41 | 11 | 2 | 0.14 | 0.46 |
2 | Singh V | 2007 | North | 104 | 119 | 38 | 68 | 0 | 9 | 0.22 | 0.86 |
3 | Surekha | 2009 | South | 78 | 122 | 169 | 127 | 1 | 0 | 0.26 | <0.001 |
4 | Chacko | 2005 | South | 67 | 90 | 34 | 19 | 11 | 3 | 0.11 | 0.13 |
5 | Naushad | 2011 | South | 221 | 144 | 108 | 99 | 13 | 10 | 0.24 | 0.16 |
6 | Kiruthiga | 2011 | South | 34 | 48 | 4 | 2 | 12 | 0 | 0.02 | 0.89 |
M4 (631/611)* | CC | CA | AA | ||||||||
1 | Singh N | 2007 | North | 129 | 145 | 12 | 15 | 6 | 2 | 0.06 | 0.04 |
2 | Singh V | 2007 | North | 93 | 43 | 43 | 52 | 6 | 2 | 0.14 | 0.24 |
3 | Naushad | 2011 | South | 274 | 240 | 68 | 13 | 0 | 0 | 0.03 | 0.03 |
* CYP1A1 polymorphism (cases / controls); s: sporadic; f: familial; maf: minor allele frequency; HWE: Hardy-Weinberg Equilibrium p values. |
Premenopausal | Postmenopausal | |||||||||||||||||
First author | Year | Indian Region | Case | Control | Case | Control | Case | Control | maf | HWE | Case | Control | Case | Control | Case | Control | maf | HWE |
M1 | TT | TC | CC | TT | TC | CC | ||||||||||||
N (cases/controls) | (279/309) | (230/251) | ||||||||||||||||
Singh N | 2007 | North | 40 | 57 | 17 | 21 | 5 | 3 | 0.17 | 0.55 | 50 | 30 | 18 | 21 | 13 | 8 | 0.31 | 0.18 |
Singh V | 2007 | North | 53 | 58 | 43 | 35 | 9 | 13 | 0.31 | 0.35 | 17 | 35 | 17 | 35 | 3 | 10 | 0.34 | 0.79 |
Chacko | 2005 | South | 29 | 44 | 58 | 58 | 25 | 0 | 0.35 | 0.14 | 33 | 46 | 54 | 54 | 25 | 12 | 0.35 | 0.51 |
M2 | AA | AG | GG | AA | AG | GG | ||||||||||||
N (cases/controls) | (277/309) | (232/251) | ||||||||||||||||
Singh N | 2007 | North | 39 | 51 | 16 | 29 | 5 | 1 | 0.19 | 0.16 | 56 | 49 | 23 | 9 | 4 | 1 | 0.09 | 0.45 |
Singh V | 2007 | North | 80 | 70 | 25 | 41 | 0 | 5 | 0.22 | 0.74 | 24 | 49 | 13 | 27 | 0 | 4 | 0.22 | 0.91 |
Chacko | 2005 | South | 31 | 45 | 57 | 58 | 23 | 9 | 0.34 | 0.10 | 36 | 45 | 57 | 58 | 23 | 9 | 0.34 | 0.10 |
M4 | CC | CA | AA | CC | CA | AA | ||||||||||||
N (cases/controls) | (168/197) | (118/139) | ||||||||||||||||
Singh N | 2007 | North | 57 | 70 | 3 | 10 | 2 | 1 | 0.07 | 0.37 | 66 | 54 | 9 | 4 | 6 | 1 | 0.05 | 0.02 |
Singh V | 2007 | North | 67 | 83 | 34 | 32 | 4 | 1 | 0.15 | 0.27 | 36 | 59 | 9 | 20 | 2 | 1 | 0.14 | 0.63 |
maf: minor allele frequency; HWE: Hardy-Weinberg Equilibrium p values. |
First author | Year | Indian city | Indian region | State of controls | Source of controls | Matching criteria | NOS |
Singh N | 2007 | Lucknow | North | normal | Hospital | menopause, diet | 6 |
Singh V | 2007 | Lucknow | North | healthy | Population | residence | 7 |
Surekha | 2009 | Hyderabad | South | healthy | Population | age | 7 |
Chacko | 2005 | Kerala | South | healthy | Hospital | age, menopause | 6 |
Syamala* | 2010 | Thiruvanan | South | cancer-free | Hospital | residence, gender | 6 |
Naushad | 2011 | Hyderabad | South | healthy | Hospital | no mention | 4 |
Kiruthiga | 2011 | Madurai | South | no mention | Population | no mention | 5 |
* Separate familial and sporadic data were considered as two studies; NOS: Newcastle-Ottawa Score. |
Test of association | Test of heterogeneity | ||||||
Genetic model | N | OR | 95% CI | pA | pB | I2 (%) | |
M1 | Homozygous | 7 | 2.31 | 0.98-5.46 | 0.06 | <0.0001 | 88 |
Recessive | 7 | 2.02 | 0.95-4.29 | 0.07 | <0.0001 | 85 | |
Dominant | 7 | 3.43 | 1.41-8.35 | 0.007 | <0.0001 | 96 | |
Codominant | 7 | 1.61 | 1.03-2.52 | 0.04 | <0.0001 | 92 | |
M2 | Homozygous | 6 | 2.53 | 0.65-9.95 | 0.18 | 0.003 | 72 |
Recessive | 6 | 2.39 | 0.68-8.42 | 0.18 | 0.009 | 68 | |
Dominant | 6 | 1.50 | 0.83-2.74 | 0.18 | <0.0001 | 89 | |
Codominant | 6 | 1.51 | 0.89-2.58 | 0.13 | <0.0001 | 90 | |
M4 | Homozygous | 2 | 3.92 | 1.25-12.33 | 0.02 | 0.79 | 0 |
Recessive | 2 | 3.81 | 1.22-11.96 | 0.02 | 0.84 | 0 | |
Dominant | 2 | 1.96 | 0.87-4.43 | 0.11 | 0.003 | 83 | |
Codominant | 2 | 2.02 | 1.03-3.95 | 0.04 | 0.01 | 78 | |
N: number of studies; OR: odds ratio; CI: confidence interval; pA and pB: p values for association and heterogeneity with significance set at < 0.05 and < 0.10, respectively; Values in bold indicate significant tests of association; R: Random-effects model; F: Fixed-effects model. |
South India | North India | ||||||||||
n = 5 | n = 2 | ||||||||||
Test of association | Test of heterogeneity | Test of association | Test of heterogeneity | ||||||||
OR | 95% CI | pA | pB | I2 | OR | 95% CI | pA | pB | I2 | pi | |
Homozygous | 4.34 | 3.07-6.12 | <0.0001 | 0.20 | 33 | 0.64 | 0.22-1.93 | 0.43 | 0.03 | 78 | 0.019 |
Recessive | 3.55 | 2.54-4.95 | <0.0001 | 0.50 | 0 | 0.69 | 0.20-2.34 | 0.55 | 0.01 | 83 | 0.23 |
Dominant | 3.73 | 0.88-15.71 | 0.007 | <0.0001 | 97 | 2.81 | 2.16-3.65 | <0.0001 | 0.63 | 0 | >1 |
Codominant | 1.17 | 1.00-1.37 | 0.05 | <0.0001 | 86 | 2.72 | 2.23-3.31 | <0.0001 | 0.65 | 0 | >1 |
Premenopausal | Postmenopausal | ||||||||||
n = 3 | n = 3 | ||||||||||
Test of association | Test of heterogeneity | Test of association | Test of heterogeneity | ||||||||
OR | 95% CI | pA | pB | I2 | OR | 95% CI | pA | pB | I2 | pi | |
Homozygous | 1.87 | 0.63-5.53 | 0.26 | 0.04 | 68 | 1.35 | 0.53-3.41 | 0.53 | 0.09 | 58 | >1 |
Recessive | 1.67 | 0.65-4.30 | 0.29 | 0.07 | 62 | 1.54 | 0.92-2.60 | 0.10 | 0.19 | 40 | >1 |
Dominant | 1.30 | 0.92-1.81 | 0.13 | 0.23 | 31 | 1.03 | 0.56-1.86 | 0.94 | 0.09 | 58 | >1 |
Codominant | 1.31 | 0.87-1.97 | 0.20 | 0.08 | 60 | 1.08 | 0.67-1.75 | 0.74 | 0.06 | 65 | >1 |
OR: odds ratio; CI: confidence interval; pA and pB: p values for association and heterogeneity with significance set at < 0.05 and < 0.10, respectively; pi: p values for tests of interaction after Bonferroni correction; Values in bold are significant in the tests of association and tests of interaction |
South India | North India | ||||||||||
n = 5 | n = 2 | ||||||||||
Test of association | Test of heterogeneity | Test of association | Test of heterogeneity | ||||||||
OR | 95% CI | pA | pB | I2 | OR | 95% CI | pA | pB | I2 | pi | |
Homozygous | 3.72 | 0.55-25.44 | 0.18 | 0.006 | 81 | 1.36 | 0.06-31.06 | 0.85 | 0.009 | 79 | >1 |
Recessive | 3.34 | 0.60-18.59 | 0.17 | 0.02 | 76 | 1.26 | 0.06-24.67 | 0.88 | 0.01 | 77 | >1 |
Dominant | 2.37 | 0.63-8.92 | 0.20 | <0.0001 | 92 | 1.17 | 0.54-2.55 | 0.69 | <0.0001 | 89 | >1 |
Codominant | 2.81 | 0.75-10.51 | 0.12 | <0.0001 | 94 | 1.10 | 0.58-2.09 | 0.76 | 0.0001 | 89 | >1 |
Premenopausal | Postmenopausal | ||||||||||
n = 3 | n = 3 | ||||||||||
Test of association | Test of heterogeneity | Test of association | Test of heterogeneity | ||||||||
OR | 95% CI | pA | pB | I2 | OR | 95% CI | pA | pB | I2 | pi | |
Homozygous | 1.69 | 0.20-14.35 | 0.63 | 0.03 | 73 | 1.79 | 0.89-3.59 | 0.10 | 0.30 | 17 | >1 |
Recessive | 1.74 | 0.25-12.11 | 0.58 | 0.04 | 69 | 1.62 | 0.85-3.11 | 0.14 | 0.35 | 6 | >1 |
Dominant | 0.90 | 0.42-1.94 | 0.79 | 0.009 | 79 | 1.43 | 0.96-2.11 | 0.08 | 0.23 | 33 | >1 |
Codominant | 0.99 | 0.46-2.14 | 0.98 | 0.0007 | 86 | 1.36 | 1.01-1.83 | 0.04 | 0.11 | 54 | >1 |
OR: odds ratio; CI: confidence interval; pA and pB: p values for association and heterogeneity with significance set at < 0.05 and < 0.10, respectively; pi: p values for tests of interaction after Bonferroni correction; Values in bold are significant in the tests of association. |
Premenopausal | Postmenopausal | ||||||||||
n = 2 | n = 2 | ||||||||||
OR | 95% CI | pA | pB | I2 | OR | 95% CI | pA | pB | I2 | pi | |
Homozygous | 3.70 | 0.73-18.63 | 0.11 | 0.67 | 0 | 4.77 | 0.93-24.41 | 0.06 | 0.96 | 0 | >1 |
Recessive | 3.63 | 0.72-18.30 | 0.12 | 0.74 | 0 | 4.59 | 0.90-23.44 | 0.07 | 0.99 | 0 | >1 |
Dominant | 1.16 | 0.70-1.90 | 0.56 | 0.13 | 56 | 1.61 | 0.84-3.11 | 0.15 | 0.30 | 6 | >1 |
Codominant | 1.26 | 0.81-1.95 | 0.30 | 0.23 | 31 | 1.83 | 1.04-3.24 | 0.04 | 0.23 | 30 | >1 |
OR: odds ratio; CI: confidence interval; pA and pB: p values for association and heterogeneity with significance set at < 0.05 and < 0.10, respectively; pi: p values for tests of interaction after Bonferroni correction; Values in bold are significant in the tests of association. |
First author | Year | Indian Region | Case | Control | Case | Control | Case | Control | maf | HWE | |
M1 (1151/1527)* | TT | TC | CC | ||||||||
1 | Singh N | 2007 | North | 94 | 94 | 35 | 53 | 17 | 15 | 0.26 | 0.07 |
2 | Singh V | 2007 | North | 70 | 93 | 60 | 80 | 12 | 43 | 0.38 | 0.001 |
3 | Chacko | 2005 | South | 62 | 90 | 39 | 18 | 11 | 4 | 0.12 | 0.02 |
4 | Syamala s | 2010 | South | 93 | 243 | 87 | 106 | 39 | 18 | 0.19 | 0.22 |
5 | Syamala f | 2010 | South | 55 | 243 | 56 | 106 | 29 | 18 | 0.19 | 0.22 |
6 | Naushad | 2011 | South | 168 | 140 | 125 | 98 | 49 | 15 | 0.25 | 0.69 |
7 | Kiruthiga | 2011 | South | 30 | 36 | 17 | 12 | 3 | 2 | 0.16 | 0.45 |
M2 (1039/1022)* | AA | AG | GG | ||||||||
1 | Singh N | 2007 | North | 98 | 119 | 36 | 41 | 11 | 2 | 0.14 | 0.46 |
2 | Singh V | 2007 | North | 104 | 119 | 38 | 68 | 0 | 9 | 0.22 | 0.86 |
3 | Surekha | 2009 | South | 78 | 122 | 169 | 127 | 1 | 0 | 0.26 | <0.001 |
4 | Chacko | 2005 | South | 67 | 90 | 34 | 19 | 11 | 3 | 0.11 | 0.13 |
5 | Naushad | 2011 | South | 221 | 144 | 108 | 99 | 13 | 10 | 0.24 | 0.16 |
6 | Kiruthiga | 2011 | South | 34 | 48 | 4 | 2 | 12 | 0 | 0.02 | 0.89 |
M4 (631/611)* | CC | CA | AA | ||||||||
1 | Singh N | 2007 | North | 129 | 145 | 12 | 15 | 6 | 2 | 0.06 | 0.04 |
2 | Singh V | 2007 | North | 93 | 43 | 43 | 52 | 6 | 2 | 0.14 | 0.24 |
3 | Naushad | 2011 | South | 274 | 240 | 68 | 13 | 0 | 0 | 0.03 | 0.03 |
* CYP1A1 polymorphism (cases / controls); s: sporadic; f: familial; maf: minor allele frequency; HWE: Hardy-Weinberg Equilibrium p values. |
Premenopausal | Postmenopausal | |||||||||||||||||
First author | Year | Indian Region | Case | Control | Case | Control | Case | Control | maf | HWE | Case | Control | Case | Control | Case | Control | maf | HWE |
M1 | TT | TC | CC | TT | TC | CC | ||||||||||||
N (cases/controls) | (279/309) | (230/251) | ||||||||||||||||
Singh N | 2007 | North | 40 | 57 | 17 | 21 | 5 | 3 | 0.17 | 0.55 | 50 | 30 | 18 | 21 | 13 | 8 | 0.31 | 0.18 |
Singh V | 2007 | North | 53 | 58 | 43 | 35 | 9 | 13 | 0.31 | 0.35 | 17 | 35 | 17 | 35 | 3 | 10 | 0.34 | 0.79 |
Chacko | 2005 | South | 29 | 44 | 58 | 58 | 25 | 0 | 0.35 | 0.14 | 33 | 46 | 54 | 54 | 25 | 12 | 0.35 | 0.51 |
M2 | AA | AG | GG | AA | AG | GG | ||||||||||||
N (cases/controls) | (277/309) | (232/251) | ||||||||||||||||
Singh N | 2007 | North | 39 | 51 | 16 | 29 | 5 | 1 | 0.19 | 0.16 | 56 | 49 | 23 | 9 | 4 | 1 | 0.09 | 0.45 |
Singh V | 2007 | North | 80 | 70 | 25 | 41 | 0 | 5 | 0.22 | 0.74 | 24 | 49 | 13 | 27 | 0 | 4 | 0.22 | 0.91 |
Chacko | 2005 | South | 31 | 45 | 57 | 58 | 23 | 9 | 0.34 | 0.10 | 36 | 45 | 57 | 58 | 23 | 9 | 0.34 | 0.10 |
M4 | CC | CA | AA | CC | CA | AA | ||||||||||||
N (cases/controls) | (168/197) | (118/139) | ||||||||||||||||
Singh N | 2007 | North | 57 | 70 | 3 | 10 | 2 | 1 | 0.07 | 0.37 | 66 | 54 | 9 | 4 | 6 | 1 | 0.05 | 0.02 |
Singh V | 2007 | North | 67 | 83 | 34 | 32 | 4 | 1 | 0.15 | 0.27 | 36 | 59 | 9 | 20 | 2 | 1 | 0.14 | 0.63 |
maf: minor allele frequency; HWE: Hardy-Weinberg Equilibrium p values. |