Genetic Variants in Nicotine Addiction and Alcohol Metabolism Genes, Oral Cancer Risk and the Propensity to Smoke and Drink Alcohol: A Replication Study in India

Abstract

Background

Genetic variants in nicotinic acetylcholine receptor and alcohol metabolism genes have been associated with propensity to smoke tobacco and drink alcohol, respectively, and also implicated in genetic susceptibility to head and neck cancer. In addition to smoking and alcohol, tobacco chewing is an important oral cancer risk factor in India. It is not known if these genetic variants influence propensity or oral cancer susceptibility in the context of this distinct etiology.

Methods

We examined 639 oral and pharyngeal cancer cases and 791 controls from two case-control studies conducted in India. We investigated six variants known to influence nicotine addiction or alcohol metabolism, including rs16969968 (CHRNA5), rs578776 (CHRNA3), rs1229984 (ADH1B), rs698 (ADH1C), rs1573496 (ADH7), and rs4767364 (ALDH2).

Results

The CHRN variants were associated with the number of chewing events per day, including in those who chewed tobacco but never smoked (P =  0.003, P =  0.01 for rs16969968 and rs578776 respectively). Presence of the variant allele contributed to approximately 13% difference in chewing frequency compared to non-carriers. While no association was observed between rs16969968 and oral cancer risk (OR =  1.01, 95% CI =  0.83– 1.22), rs578776 was modestly associated with a 16% decreased risk of oral cancer (OR =  0.84, 95% CI =  0.72– 0.98). There was little evidence for association between polymorphisms in genes encoding alcohol metabolism and oral cancer in this population.

Conclusion

The association between rs16969968 and number of chewing events implies that the effect on smoking propensity conferred by this gene variant extends to the use of smokeless tobacco.

About the Authors:

Devasena Anantharaman

Affiliation: Genetic Epidemiology Group, International Agency for Research on Cancer, Lyon, France

Amélie Chabrier

Affiliation: Genetic Cancer Susceptibility Group, International Agency for Research on Cancer, Lyon, France

Valérie Gaborieau

Affiliation: Genetic Epidemiology Group, International Agency for Research on Cancer, Lyon, France

Silvia Franceschi

Affiliation: Infections and Cancer Epidemiology Group, International Agency for Research on Cancer, Lyon, France

Rolando Herrero

Affiliation: Prevention and Implementation Group, International Agency for Research on Cancer, Lyon, France

Thangarajan Rajkumar

Affiliation: Department of Molecular Oncology, Cancer Institute, Chennai, India

Tanuja Samant

Affiliation: Mahimkar Lab, Advanced Center for Treatment Research and Education in Cancer, Tata Memorial Center, Navi Mumbai, India

Manoj B. Mahimkar

Affiliation: Mahimkar Lab, Advanced Center for Treatment Research and Education in Cancer, Tata Memorial Center, Navi Mumbai, India

Paul Brennan

Affiliation: Genetic Epidemiology Group, International Agency for Research on Cancer, Lyon, France

James D. McKay

* E-mail: [email protected]

Affiliation: Genetic Cancer Susceptibility Group, International Agency for Research on Cancer, Lyon, France

Introduction

Cancers of the oral cavity and pharynx contribute to nearly 400,000 new cases each year worldwide, more than half of which occur in India. Each year over 200,000 die of the disease, and over a third of these deaths occur in India [1]. Tobacco use and alcohol consumption are the key established risk factors for oral cancer, with the use of smokeless tobacco being particularly important in the Indian population [2]. Exposure to human papillomavirus is becoming increasingly important to cancers of the oropharynx [3], [4].

Genome-wide association studies (GWAS) have successfully identified disease susceptibility loci to various complex diseases [5]. Lung cancer GWAS and nicotine addiction studies have identified the 15q25 locus harbouring the nicotinic acetylcholine receptor (CHRN) gene cluster [6]–[8]. These genes code for receptors expressed in neuronal and other epithelial cells that bind to nicotine and nicotine derivatives [9], [10]. Two CHRN receptor subunit variants, rs16969968 and rs578776 have been consistently associated with lung cancer risk and smoking behavior in several populations [11]–[15]. Homozygous carriers of the rs1669968 rare allele have been reported to smoke approximately 1.2 cigarettes more per day [13]. Further, this variant has also been associated with increased risk of Upper Aero-Digestive Tract (UADT) cancer. UADT cancer GWAS and candidate gene association studies have identified genetic variants in the 4q (rs1229984, rs698, rs1573496) and 12q (rs4767364) loci containing genes involved in alcohol metabolism [16], [17]. The balance between alcohol dehydrogenase and aldehyde dehydrogenase activities has been suggested to regulate blood acetaldehyde concentrations that determine the unpleasant symptoms associated with alcohol consumption, thus impacting the ability to consume alcohol and potentially, cancer risk [18], [19]. The alcohol and aldehyde dehydrogenase genes (ADH and ALDH, respectively) have been associated with head and neck cancer risk [16], [17], [20]. Although alcohol is an important and established risk factor for oral cancer in India [21]–[24], there is a paucity of data on the association of ADH and ALDH2 variants in this population.

In this study we aimed to (i) determine whether the association between CHRN variants and propensity to smoke extends to non-smoking forms of tobacco use (ii) examine if CHRN genetic variants influence susceptibility to oral cancer risk in India, (iii) clarify the association between the potential causal variants in the 4q23 (ADH) and 12q (ALDH2) locus and the risk of oral cancer in India.

Materials and Methods

Ethics statement

Written informed consent was obtained from all participants. The IARC multi-center study was approved centrally by the Cancer Institute Ethical Committee at Chennai, India and the overall study was approved by the IARC Ethical Review Committee at the International Agency for Research on Cancer (IARC), Lyon, France. The Mumbai study was approved by the Hospital Ethics Committee of the Tata Memorial Hospital & Cancer Research Institute at Mumbai, India.

Study description

The present analysis utilized cases and controls from two previously conducted studies in India: the International multicenter oral cancer study (IARC study), and the Mumbai study. The IARC study was coordinated by the International Agency for Research on Cancer (IARC) across nine countries during 1996 to 1999. The study recruited cases and controls from three centers across India, including Bangalore, Chennai and Trivandrum. Cases were histologically or pathologically confirmed primary cancers of the oral cavity and oropharynx. Hospital based controls, one for each case frequency matched for age (5 year period), sex and center were recruited. The overall study participation rate was 93%, details have been described elsewhere [25]–[27]. The Mumbai study was conducted by the Cancer Research Institute, Advanced Center for Treatment, Research and Education in Cancer at the Tata Memorial Center, during 2001 to 2004. Cases were histologically confirmed diagnosis of oral cavity cancers and appropriate controls, frequency matched for age (5 year intervals), gender and tobacco habits were recruited from the Tata Memorial hospital blood bank and nearby dental clinics [28]. Lifestyle and demographic data were collected upon personal interview in both studies.

Laboratory methods

Genomic DNA was extracted from blood samples using the QIAamp 96 DNA Blood Kit (Qiagen, Hilden, Germany), or with Puregene chemistry (Gentra Systems, Minneapolis, MN, U.S.A.) on an Autopure instrument (Gentra Systems) in the IARC study. In the Mumbai study, DNA was extracted following standard protocol by Sambrook et al [29]. DNA concentrations were measured using PicoGreen double-stranded DNA quantification kits (Molecular Probes, Leiden, The Netherlands). All polymorphisms were analyzed centrally at IARC using the 5’ exonuclease assay. PCR primers and TaqMan probes were synthesized by Applied Biosystems (Foster City, CA). To ensure quality control, DNA samples from case patients and control subjects were randomly distributed on each PCR plate, and all genotyping was conducted by personnel blinded to the case-control status. Ten percent of the study subjects (i.e., both case patients and control subjects), were randomly selected for re-genotyping of each polymorphism to examine the reliability of the genotyping assays, all assays had concordance of >99%. The final analyses included samples that passed stringent quality control (QC) criteria. Subjects were excluded if they failed two or more of the six genotyping assays performed. Further, only cases with >95% genotyping success rate were included. The original IARC study recruited 582 cases and 582 controls from India. Of these, 935 subjects were genotyped, of which 872 passed the described QC criteria and were included in the present analysis. In the Mumbai study, of the 576 subjects genotyped, 226 cases and 332 controls cleared the described QC and were included in this analysis. A total of 639 cases and 791 controls were included in this analysis following exclusions described above.

Statistical methods

Data on alcohol consumption, tobacco use, chewing and smoking were collected in both the participating studies. The definition of drinker, smoker and chewer were consistent across both studies and included individuals who responded confirming ever use of alcohol, ever smoking any tobacco product (cigarette, cigar, pipe, hooka, cigarette and/or bidi), ever use of any smokeless tobacco product with or without betel quid. Several forms of smokeless tobacco are consumed in India, while khaini, mawa, mishri, zarda, gudakhu and gutkha are all formulations containing tobacco with different flavoring agents, pan masala is an areca nut product that does not contain tobacco, to which tobacco may be added externally. Duration and frequency of exposure to tobacco and alcohol were calculated in each study as described previously [25]–[28]. Briefly, smoking duration was calculated as the difference between age at start and age at stop of smoking, taking into account overlapping periods of various types of smoking (except in the Mumbai study where only cigarette and bidi use was recorded, and very few overlapping exposures were observed). Smoking frequency was calculated as number of cigarette equivalents smoked per day, with 1 cigar =  4 cigarettes, 1 pipe =  5 cigarettes and 1 bidi =  0.5 cigarettes. Chewing included use of smokeless tobacco products, including betel quid, areca nut with &/or without tobacco and pan masala in the IARC study, and tobacco with lime, pan masala, khaini, mawa, mishri, zarda, gudakhu, snuff and gutkha in the Mumbai study. Chewing frequency was calculated as the cumulative number of chewing events per day, taking into account overlapping periods of exposure. For the calculation of lifetime alcohol consumption, the percentage of ethanol was estimated as 40% for spirits (whisky, gin, rum, brandy, arrack and country liquor), 3% for beer and 5% for toddy (except in the Mumbai study). Smoking, chewing and alcohol related oral cancer risk was calculated using logistic regression models adjusted for age, sex, level of education and center. The corresponding attributable fractions were calculated based on the formula: p (ec)x(OR-1)/OR; where p(ec) is the proportion of exposed among cases [30], [31]. Chi-square tests were performed to compare allele frequencies observed in the present study with our previously published results [16], [17]. The associations between the genetic variants and oral cancer were estimated by logistic regression adjusted for age, sex and center and tobacco use where appropriate. The association between the CHRN variants and propensity to use tobacco was assessed using linear regression models with log transformed number of cigarette equivalents smoked per day, or number of chewing events per day as the outcome variable. Results represent the fold-change in the number of cigarettes smoked per day or the number of chewing events per day for each copy of the rare allele. These analyses were adjusted for age, sex, center and case-control status. Similar analyses were performed to examine the association between ADH/ALDH2 variants and the propensity to drink alcohol, measured as number of drinks of ethanol consumed per day. All statistical analyses were performed using STATA version 11 (STATA Corp, College Station, Tx), and the reported p-values are two sided.

Results

Table 1 describes the characteristics of the study population. Of the 639 cases and 791 controls included in the present analysis, the IARC international study contributed the largest proportion of subjects (61%). Majority of the cases were male (65%) and had completed lower levels of education and were cancers of the oral cavity (96%). Tobacco smoking and chewing were associated with increased risk of oral cancer (OR =  5.23, 95% CI =  3.39– 8.07 and OR =  8.30, 95% CI =  5.78– 11.93, respectively), and the combined use of smoking and smokeless tobacco products conferred a 12-fold increased risk of oral cancer (OR =  12.40, 95% CI =  7.80– 19.70). Ever use of alcohol was associated with a 1.67-fold increased risk of oral cancer (OR =  1.67, 95% CI =  1.25– 2.22).

[Figure omitted. See PDF.]

Table 1. Description of study group.

https://doi.org/10.1371/journal.pone.0088240.t001

The distribution of all six variants examined in the present analysis was consistent with expectations of genotype distributions based on the Hardy-Weinberg principle. The allele frequencies were similar across the participating studies (data not shown). We found that the variants in alcohol dehydrogenase genes- rs1229984 (ADH1B), rs698 (ADH1C) and rs1573496 (ADH7) and the nicotine dependence gene variant rs16969968 were less common among Indians compared to previous studies in predominantly western populations. On the other hand, rs578776 was more common among Indians (table2).

[Figure omitted. See PDF.]

Table 2. Variants examined in the present study.

https://doi.org/10.1371/journal.pone.0088240.t002

We examined the association between CHRN variants and the propensity to use tobacco; the results are summarized in table 3. The carriers of the rs16969968 rare allele (A) showed some tendency to smoke more, although were not statistically significant. We also observed an increase in the propensity to chew tobacco with the rs16969968 variant (p =  0.03). Presence of each copy of the A allele corresponded to an increase by 13% more chewing events per day (fold-change =  1.13, 95% CI =  1.01– 1.25) compared to non-carriers. The association with chewing remained robust upon further adjustment for smoking frequency (fold-change =  1.13, 95% CI =  1.01– 1.26) and the association was consistent among never smokers who chewed tobacco (p =  0.003) (fold-change =  1.19, 95% CI =  1.06– 1.34). We observed a similar decrease in the propensity to chew among carriers of the rs578776 rare variant (fold-change =  0.91, 95% CI =  0.84– 0.99), the association being stronger among chewers who never smoked.

[Figure omitted. See PDF.]

Table 3. CHRN variants and propensity to smoke or chew tobacco.

https://doi.org/10.1371/journal.pone.0088240.t003

No association was observed for rs16969968 and oral cancer risk in this study overall (OR =  1.01, 95% CI =  0.83– 1.22), by gender (p-heterogeneity =  0.33). Although no effect modification by smoking or chewing status was observed, a statiscally non-significant oddds ratio of 1.22 (95% CI =  0.69– 2.13) was observed in chewers and smokers (figure 1). We observe little evidence for a statistical interaction between the rs16969968 variant and tobacco use in the risk of oral cancer (OR: 1.18, 95% CI: 0.68– 2.07). On the other hand, rs578776 was associated with decreased risk of oral cancer in this study (OR =  0.84, 95% CI =  0.72– 0.98). A gene-dose effect was observed for carriers of the rare allele (p-trend =  0.02). Although no heterogeneity was observed based on tobacco habits, a decreased risk was also observed among never smokers and never chewers (OR =  0.64, 95% CI =  0.43–0.96) (figure 1).

[Figure omitted. See PDF.]

Figure 1. CHRN variants and risk of oral cancer.

Panel A shows the association between rs16969968 and oral cancer; Panel B shows the association between rs578776 and oral cancer. Forest plot represents odds ratios derived from the log additive multivariate model adjusted for age, sex, level of education and center, as appropriate.

https://doi.org/10.1371/journal.pone.0088240.g001

We examined variants in genes related to alcohol metabolism since they were strongly associated with head and neck cancer in previous studies. None of the four variants (rs1229984, rs698, rs1573496, and rs4767364) examined, showed any association with oral cancer risk (table 4). We also performed restricted analyses among ever drinkers given the important effect modification by alcohol consumption. Although the lack of association remained consistent, the observed estimates followed the direction anticipated based on expectations from previous GWAS and candidate gene studies. Stratified analyses based on gender did not reflect any heterogeneity in the observed lack of associations (table 4). Further, none of the genetic variants in alcohol metabolism genes appeared to influence alcohol consumption levels (Table S1).

[Figure omitted. See PDF.]

Table 4. Association between ADH & ALDH2 variants and risk of oral cancer.

https://doi.org/10.1371/journal.pone.0088240.t004

Discussion

The CHRN variants, rs16969968 and rs578776, have been associated with nicotine addiction and lung cancer risk and these associations appear to be independent [15], [32]. Here, we demonstrate that the association between CHRN variants also extends to the propensity to chew tobacco, a relevant oral cancer risk factor in the Indian population.

Several formulations of smokeless tobacco are commercially available and commonly used in India; these were grouped together in this study. The number of chewing events per day was recorded in both studies included in this analysis; however it is important to note that this does not reflect the amount of product chewed each time and hence, constitutes a potential source of misclassification of exposure. Given that we observe a measurable association with such a crude estimate implies that the true effect of the CHRN variants and chewing dependence could be potentially higher than that observed in the present study. Previous studies have demonstrated that UADT cancer risk associated with rs16969968 is higher among women and for cancers of the oral cavity and larynx [13], [33]. In this study, the increased propensity to chew tobacco did not translate to an increased risk of oral cancer for rs16969968; although rs578776 was associated with both decreased chewing propensity and a modest decrease in oral cancer risk. Interestingly, rs578776 was also associated with decreased oral cancer risk in never tobacco users, albeit based on small numbers. This association may be influenced by misclassification of exposure or potentially through a tobacco independent effect. Larger focused studies will be required to address this issue.

The ADH genes, correlated as a consequence of linkage disequilibrium, are located in the 4q23 locus. Three alleles, rs1229984, rs698 and rs1573496 within the ADH locus have been associated with UADT cancer risk in Europeans and other Asian populations [17], [34]–[36]. While there has been some evidence to suggest ADH variants may influence oral and neck cancer risk in India, in particular rs698 [37], we observed no association between rs1229984, rs698 and oral cancer risk. Similarly the ALDH2 variant rs4767364, identified through a recent UADT GWAS [17], was not associated with oral cancer in this study. The fraction of alcohol consumers in this study was comparable to other reports from this population, but is notably lower than reported among Europeans [38], [39]. Further, the rarity of the risk alleles in the Indian population suggests that the contribution of alcohol and alcohol related variants to oral cancer susceptibility in India is likely to be limited.

To the best of our knowledge, we report for the first time that the association between 15q25 locus and propensity to smoke extends to tobacco chewing. We also provide suggestive evidence that 15q25 region containing the CHRN gene cluster may contribute to oral cancer susceptibility in India. In this study, rs16969968 variant was associated with increased propensity to chew tobacco, and rs578776 was associated with decreased chewing propensity and oral cancer risk. Larger studies will be required to clarify the extent of risk modification mediated by propensity to use tobacco. This study also indicates that variants in alcohol metabolizing genes appear to be less important to oral cancer susceptibility in India. These data reiterate the importance of disease etiology, the relative contribution of risk factors and their impact in replication of GWAS identified risk loci across populations.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Variants in alcohol metabolizing genes and alcohol consumption levels. Fold-change were derived from linear regression models adjusted for age, sex, center and case-control status. Log transformed number of alcoholic drinks per day was treated as the outcome variable for each of the respective genetic variant as the explanatory.

https://doi.org/10.1371/journal.pone.0088240.s001

(DOCX)

Acknowledgments

The authors thank all the participants who took part in this research.

Author Contributions

Conceived and designed the experiments: JMK PB DA. Performed the experiments: AC DA JMK. Analyzed the data: DA VG AC. Contributed reagents/materials/analysis tools: PB JMK SF RH TR TS MBM. Wrote the paper: DA JMK PB SF RH TR TS MBM.

Citation: Anantharaman D, Chabrier A, Gaborieau V, Franceschi S, Herrero R, Rajkumar T, et al. (2014) Genetic Variants in Nicotine Addiction and Alcohol Metabolism Genes, Oral Cancer Risk and the Propensity to Smoke and Drink Alcohol: A Replication Study in India. PLoS ONE 9(2): e88240. https://doi.org/10.1371/journal.pone.0088240

References

1. Ferlay J, Soerjomataram I, Ervik M, Dikshit R, Eser S, et al.. (2012) GLOBOCAN 2012 v1.0, Cancer Incidence and Mortality Worldwide: IARC CancerBase No. 11 [Internet]. International Agency for Research on Cancer.

2. (2012) IARC Monographs on the Evaluation of Carcinogenic Risks to Humans: A Review of Human Carcinogens: Personal Habits and Indoor Combustions. Lyon: IARC. 579 p.

3. (2012) IARC Monographs on the Evaluation of Carcinogenic Risks to Humans: A Review of Human Carcinogens: Biological Agents Lyon: IARC. 475 p.

4. Anantharaman D, Gheit T, Waterboer T, Abedi-Ardekani B, Carreira C, et al. (2013) Human papillomavirus infections and upper aero-digestive tract cancers: the ARCAGE study. J Natl Cancer Inst 105: 536–545.

5. Hindorff LA MJEBI, Wise A, Junkins HA, Hall PN, Klemm AK, and Manolio TA A Catalog of Published Genome-Wide Association Studies.

6. Genome-wide meta-analyses identify multiple loci associated with smoking behavior. Nat Genet 42: 441–447.

7. Amos CI, Wu X, Broderick P, Gorlov IP, Gu J, et al. (2008) Genome-wide association scan of tag SNPs identifies a susceptibility locus for lung cancer at 15q25.1. Nat Genet 40: 616–622.

8. Hung RJ, McKay JD, Gaborieau V, Boffetta P, Hashibe M, et al. (2008) A susceptibility locus for lung cancer maps to nicotinic acetylcholine receptor subunit genes on 15q25. Nature 452: 633–637.

9. Minna JD (2003) Nicotine exposure and bronchial epithelial cell nicotinic acetylcholine receptor expression in the pathogenesis of lung cancer. J Clin Invest 111: 31–33.

10. Wang Y, Pereira EF, Maus AD, Ostlie NS, Navaneetham D, et al. (2001) Human bronchial epithelial and endothelial cells express alpha7 nicotinic acetylcholine receptors. Mol Pharmacol 60: 1201–1209.

11. Jaworowska E, Trubicka J, Lener MR, Masojc B, Zlowocka-Perlowska E, et al. (2011) Smoking related cancers and loci at chromosomes 15q25, 5p15, 6p22.1 and 6p21.33 in the Polish population. PLoS One 6: e25057.

12. Li MD, Yoon D, Lee JY, Han BG, Niu T, et al. (2010) Associations of variants in CHRNA5/A3/B4 gene cluster with smoking behaviors in a Korean population. PLoS One 5: e12183.

13. Lips EH, Gaborieau V, McKay JD, Chabrier A, Hung RJ, et al. (2010) Association between a 15q25 gene variant, smoking quantity and tobacco-related cancers among 17 000 individuals. Int J Epidemiol 39: 563–577.

14. Saccone NL, Culverhouse RC, Schwantes-An TH, Cannon DS, Chen X, et al.. (2010) Multiple independent loci at chromosome 15q25.1 affect smoking quantity: a meta-analysis and comparison with lung cancer and COPD. PLoS Genet 6.

15. Saccone NL, Wang JC, Breslau N, Johnson EO, Hatsukami D, et al. (2009) The CHRNA5-CHRNA3-CHRNB4 nicotinic receptor subunit gene cluster affects risk for nicotine dependence in African-Americans and in European-Americans. Cancer Res 69: 6848–6856.

16. Hashibe M, McKay JD, Curado MP, Oliveira JC, Koifman S, et al. (2008) Multiple ADH genes are associated with upper aerodigestive cancers. Nat Genet 40: 707–709.

17. McKay JD, Truong T, Gaborieau V, Chabrier A, Chuang SC, et al. (2011) A genome-wide association study of upper aerodigestive tract cancers conducted within the INHANCE consortium. PLoS Genet 7: e1001333.

18. Druesne-Pecollo N, Tehard B, Mallet Y, Gerber M, Norat T, et al. (2009) Alcohol and genetic polymorphisms: effect on risk of alcohol-related cancer. Lancet Oncol 10: 173–180.

19. Seitz HK, Stickel F (2007) Molecular mechanisms of alcohol-mediated carcinogenesis. Nat Rev Cancer 7: 599–612.

20. Chang JS, Straif K, Guha N (2012) The role of alcohol dehydrogenase genes in head and neck cancers: a systematic review and meta-analysis of ADH1B and ADH1C. Mutagenesis 27: 275–286.

21. Muwonge R, Ramadas K, Sankila R, Thara S, Thomas G, et al. (2008) Role of tobacco smoking, chewing and alcohol drinking in the risk of oral cancer in Trivandrum, India: a nested case-control design using incident cancer cases. Oral Oncol 44: 446–454.

22. Cancela Mde C, Ramadas K, Fayette JM, Thomas G, Muwonge R, et al. (2009) Alcohol intake and oral cavity cancer risk among men in a prospective study in Kerala, India. Community Dent Oral Epidemiol 37: 342–349.

23. Sankaranarayanan R, Duffy SW, Nair MK, Padmakumary G, Day NE (1990) Tobacco and alcohol as risk factors in cancer of the larynx in Kerala, India. Int J Cancer 45: 879–882.

24. Subapriya R, Thangavelu A, Mathavan B, Ramachandran CR, Nagini S (2007) Assessment of risk factors for oral squamous cell carcinoma in Chidambaram, Southern India: a case-control study. Eur J Cancer Prev 16: 251–256.

25. Balaram P, Sridhar H, Rajkumar T, Vaccarella S, Herrero R, et al. (2002) Oral cancer in southern India: the influence of smoking, drinking, paan-chewing and oral hygiene. Int J Cancer 98: 440–445.

26. Herrero R, Castellsague X, Pawlita M, Lissowska J, Kee F, et al. (2003) Human papillomavirus and oral cancer: the International Agency for Research on Cancer multicenter study. J Natl Cancer Inst 95: 1772–1783.

27. Rajkumar T, Sridhar H, Balaram P, Vaccarella S, Gajalakshmi V, et al. (2003) Oral cancer in Southern India: the influence of body size, diet, infections and sexual practices. Eur J Cancer Prev 12: 135–143.

28. Anantharaman D, Chaubal PM, Kannan S, Bhisey RA, Mahimkar MB (2007) Susceptibility to oral cancer by genetic polymorphisms at CYP1A1, GSTM1 and GSTT1 loci among Indians: tobacco exposure as a risk modulator. Carcinogenesis 28: 1455–1462.

29. Sambrook J FE, Maniatis T (1989) Molecular cloning- a laboratory Manual. New York: Cold Spring Harbor Laboratory Press.

30. Miettinen OS (1974) Proportion of disease caused or prevented by a given exposure, trait or intervention. Am J Epidemiol 99: 325–332.

31. Rothman KJ GS, Lash TL (2009) Modern Epidemiology. Philadelphia: Lippincott Williams & Wilkins.

32. Timofeeva MN, McKay JD, Smith GD, Johansson M, Byrnes GB, et al. (2011) Genetic polymorphisms in 15q25 and 19q13 loci, cotinine levels, and risk of lung cancer in EPIC. Cancer Epidemiol Biomarkers Prev 20: 2250–2261.

33. Chen D, Truong T, Gaborieau V, Byrnes G, Chabrier A, et al. (2011) A sex-specific association between a 15q25 variant and upper aerodigestive tract cancers. Cancer Epidemiol Biomarkers Prev 20: 658–664.

34. Asakage T, Yokoyama A, Haneda T, Yamazaki M, Muto M, et al. (2007) Genetic polymorphisms of alcohol and aldehyde dehydrogenases, and drinking, smoking and diet in Japanese men with oral and pharyngeal squamous cell carcinoma. Carcinogenesis 28: 865–874.

35. Ding JH, Li SP, Cao HX, Wu JZ, Gao CM, et al. (2009) Polymorphisms of alcohol dehydrogenase-2 and aldehyde dehydrogenase-2 and esophageal cancer risk in Southeast Chinese males. World J Gastroenterol 15: 2395–2400.

36. Oze I, Matsuo K, Suzuki T, Kawase T, Watanabe M, et al. (2009) Impact of multiple alcohol dehydrogenase gene polymorphisms on risk of upper aerodigestive tract cancers in a Japanese population. Cancer Epidemiol Biomarkers Prev 18: 3097–3102.

37. Solomon PR, Selvam GS, Shanmugam G (2008) Polymorphism in ADH and MTHFR genes in oral squamous cell carcinoma of Indians. Oral Dis 14: 633–639.

38. Hashibe M, Boffetta P, Zaridze D, Shangina O, Szeszenia-Dabrowska N, et al. (2006) Evidence for an important role of alcohol- and aldehyde-metabolizing genes in cancers of the upper aerodigestive tract. Cancer Epidemiol Biomarkers Prev 15: 696–703.

39. Gajalakshmi V, Hung RJ, Mathew A, Varghese C, Brennan P, et al. (2003) Tobacco smoking and chewing, alcohol drinking and lung cancer risk among men in southern India. Int J Cancer 107: 441–447.