PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of jnciLink to Publisher's site
 
J Natl Cancer Inst. Nov 5, 2008; 100(21): 1488–1491.
Published online Nov 5, 2008. doi:  10.1093/jnci/djn380
PMCID: PMC2610344
NIHMSID: NIHMS82231
Intermediacy and Gene–Environment Interaction: The Example of CHRNA5-A3 Region, Smoking, Nicotine Dependence, and Lung Cancer
Sholom Wacholder,corresponding author Nilanjan Chatterjee, and Neil Caporaso
Affiliation of authors: Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, MD
corresponding authorCorresponding author.
Correspondence to: Sholom Wacholder, Division of Cancer Epidemiology and Genetics, National Cancer Institute, Bethesda, MD 20892 (e-mail: WacholdS/at/mail.NIH.gov).
During the “candidate gene era,” many studies with a small number of cases investigated a small number of variants in a limited number of genes chosen on the basis of partial knowledge of their function. These studies established only a few genetic risk factors for cancer (1) and other diseases. Remarkably, recent genome-wide association (GWA) studies have produced strong evidence that variation in a number of genomic regions affects risk of disease (2). Successful GWA studies most often used an agnostic approach to interrogating as many of the more common genetic variants as feasible (up to a million now) and large numbers of cases. They adopted a de facto criterion for genome-wide statistical significance of P values below 10−7 (3), thus guaranteeing a low chance of reporting a false positive even for quite low prior probabilities of association and weak effects (4), and they used statistical methods, such as principal components analysis (5), which protects not only against population stratification but also against poor control selection (6).
GWA studies are the culmination of advances in laboratory methodologies for detecting genetic variation, genomics, statistics, and informatics, yet they leave important scientific gaps (7). They can establish the presence of association but, by themselves, can identify neither the causal variants nor their function. Currently, they do not have power to capture the effects of the same variants that are most difficult to identify in family studies: variants with very small effects and low minor allele frequencies.
GWA studies of lung cancer are an instructive model for the study of gene the variant should not affect risk of lung cancer in those variants with small effects and low minor allele frequencies and behavioral determinants of disease because of the dominant role of smoking in the etiology of this disease. Like GWA studies of cancers of the breast and prostate, the first three GWA studies of lung cancer (810) identify several genetic variants strongly associated with disease. But a GWA study of lung cancer without smoking information cannot distinguish among three possibilities for a genetic variant found to be strongly associated with disease but about which nothing was previously known: 1) the variant increases risk of lung cancer solely through effects on smoking behavior; 2) the variant induces carcinogenesis through a molecular mechanism relevant to lung cancer without affecting smoking behavior; or 3) the variant affects both carcinogenesis and smoking behavior. Indeed, the essence of the controversy over the role of the chromosome 15q24/25.1 or CHRNA5-A3 region in lung cancer risk is the role of smoking: Does the association of 15q24/25.1 with nicotine dependence (11,12) arise because smoking is an intermediate factor, sometimes called a mediator or endophenotype, in the causal path between the genetic variant and disease? Does 15q24/25.1 itself affect lung carcinogenesis directly? Are both possibilities correct?
With the knowledge about the possible role of 15q24/25.1 on smoking behavior and the smoking data (11,12), the authors of the GWA studies (810) offer differing etiologic interpretations of the finding that 15q24/25 was associated with lung cancer. Now, following the spirit of the advice of Chanock and Hunter (13), Spitz et al. (14) use new data to argue for the third choice mentioned above, that variant rs1051370 in the 15q24/25.1 region affects both smoking behavior and risk of lung cancer. Their analysis focuses on direct effects of the variant on measured smoking phenotypes in approximately 3600 case and control subjects “from the same source as their GWA population.” They find strong evidence for a small effect on smoking intensity, as measured by cigarettes smoked per day, and provide further support for an effect on nicotine dependence (9,10,15), as measured by the Fagerstrom Test of Nicotine Dependence, but only weak or no supporting evidence for an association with other smoking phenotypes: age at initiation of smoking, duration of smoking, and sustained smoking cessation. The lack of association between the locus and bladder and renal cancers in a combined case set (14) and head and neck cancers (8), without adjustment for smoking, is consistent with a specific effect on lung carcinogenesis, but provides no further support for a more general effect on cancers related to tobacco. Overall, the receptor's ability to bind nicotine and, possibly, downstream carcinogens (8,16) may account in part for the seeming inconsistency of the results for various endpoints in Spitz (14) and earlier studies (810). For example, some laboratory and population data suggest both central (an impact of smoking mediated at least in part by an effect on nicotine dependency) and peripheral (an effect on lung carcinogenesis), mechanisms as evidenced by expression of the nicotinic acetylcholine receptors in bronchial cells (17) and that the receptors are ligands for tobacco-specific carcinogens (18).
If the only effect of a variant on disease is through nicotine dependence, there should be no effect of the variant on lung cancer in those without exposure to smoking. In never smokers, Spitz et al. (14) saw no evidence of increased risk in carriers of the variant associated with lung cancer in smokers. Unmeasured risk factors for smoking and for lung cancer, however, can distort the stratified analysis (19). Measurement errors in identifying the functional variants and in reports on smoking are perhaps more important sources of bias than confounding (20). For instance, carriers and noncarriers of a variant may metabolize tobacco differently or behave differently, so that the carcinogenic dose from smoking varies by carrier status due perhaps to depth of inhalation, even if the reports of history of cigarettes smoked per day were perfectly accurate. Also, a GWA study cannot determine the specific variant or variants in the region of the tagging SNP associated with smoking and whether the same SNP or SNPs, or one or more additional proximal variant are associated directly with disease, especially in a region of high LD like 15q24/25.
Despite impressive work on 15q24/25.1 and lung cancer (810,14), we are not close to understanding the precise mechanisms underlying the genotypic association. The lack of statistically significant associations of the region with other cancers (8,14) and with lung cancer in nonsmokers (14) does not establish lack of association, just as a statistically significant result does not prove association. Even a convincing demonstration that there is no effect of a variant on lung cancer in nonsmokers does not imply that the variant acts on disease solely by its effect on addiction to smoking; we expect the noncausal effect of a consequence of smoking, like tobacco-stained teeth, on lung cancer would also to disappear on adjustment for smoking. Variants in a gene encoding a metabolic activity that affects carcinogenicity of tobacco components also will have no effect in nonsmokers, under the assumption that there is no effect from passive smoking or exposure to other carcinogens that bind to the receptors. Or a metabolic variant might affect behavior by providing smokers with pleasant or unpleasant feedback—like ALDH2 variants for alcohol (21)—rather than be related to addiction.
Table 1 lists several other examples where there are similar fundamental questions about the role of genes and behavioral or endogenous factors in cancer etiology. In each example, intermediacy, where factor A greatly influences factor B, itself a cause of outcome C, but A does not affect C except through B, is plausible. Although intermediacy is a special case (22) of interaction, which addresses how risk factors act together to cause disease, it is not always considered in theoretical discussions (23,24).
Table 1
Table 1
Plausible intermediate factors in causal pathways from genetic variation to disease
Statistical approaches (2528) may help distinguish among potential mechanisms. They can help articulate the assumptions and define precisely what can be estimated from standard multivariate analyses as well as alternative methods to address these questions. Both epidemiological (19,22,28,29) and clinical trial (30) literature discuss the statistical problem of whether and how much of the causal effect of an exposure or treatment on an outcome is mediated through an intermediate variable. Causal modeling frameworks, including theory of counterfactual outcomes (31,32) and directed acyclic graphs (19,22,33), can examine the complex interplay among an exposure, a possible intermediate and an outcome. These models can also be useful to understand the precise assumptions needed and the pitfalls, for example, unmeasured confounders, in the standard stratified or multivariate regression methods to measure the “direct” effect of an exposure on an outcome that is not mediated through an intermediate. We anticipate that in future such causal modeling framework will provide more insight into the interrelationship between 15q, smoking, and risk of lung cancer and other examples in Table 1.
The agnostic approach underpinning the design and analysis of GWA studies is integral to their success. Now we face the challenge of unraveling the meaning of the associations we have discovered and leveraging their findings. For the chromosome 15–smoking–lung cancer association, inter-disciplinary teams will require high-quality information on environmental causes of disease; careful use of perhaps unfamiliar statistical methods; and more gritty, perhaps hypothesis-based investigation of molecular and behavioral mechanisms.
1. Dong LM, Potter JD, White E, Ulrich CM, Cardon LR, Peters U. Genetic susceptibility to cancer: the role of polymorphisms in candidate genes. JAMA. 2008;299(20):2423–2436. [PMC free article] [PubMed]
2. Manolio TA, Brooks LD, Collins FS. A HapMap harvest of insights into the genetics of common disease. J Clin Invest. 2008;118(5):1590–1605. [PMC free article] [PubMed]
3. Chanock SJ, Manolio T, Boehnke M, et al. Replicating genotype-phenotype associations. Nature. 2007;447(7145):655–660. [PubMed]
4. Wacholder S, Chanock S, Garcia-Closas M, El GL, Rothman N. Assessing the probability that a positive report is false: an approach for molecular epidemiology studies. J Natl Cancer Inst. 2004;96(6):434–442. [PubMed]
5. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006;38(8):904–909. [PubMed]
6. Yu K, Wang Z, Li Q, et al. Population substructure and control selection in genome-wide association studies. PLoS ONE. 2008;3(7)) e2551. [PMC free article] [PubMed]
7. Hemminki K, Forsti A, Lorenzo BJ. New cancer susceptibility loci: population and familial risks. Int J Cancer. 2008;123(7):1726–1729. [PubMed]
8. Hung RJ, McKay JD, Gaborieau V, et al. A susceptibility locus for lung cancer maps to nicotinic acetylcholine receptor subunit genes on 15q25. Nature. 2008;452(7187):633–637. [PubMed]
9. Amos CI, Wu X, Broderick P, et al. Genome-wide association scan of tag SNPs identifies a susceptibility locus for lung cancer at 15q25.1. Nat Genet. 2008;40(5):616–622. [PMC free article] [PubMed]
10. Thorgeirsson TE, Geller F, Sulem P, et al. A variant associated with nicotine dependence, lung cancer and peripheral arterial disease. Nature. 2008;452(7187):638–642. [PubMed]
11. Bierut LJ, Madden PA, Breslau N, et al. Novel genes identified in a high-density genome wide association study for nicotine dependence. Hum Mol Genet. 2007;16(1):24–35. [PMC free article] [PubMed]
12. Saccone SF, Hinrichs AL, Saccone NL, et al. Cholinergic nicotinic receptor genes implicated in a nicotine dependence association study targeting 348 candidate genes with 3713 SNPs. Hum Mol Genet. 2007;16(1):36–49. [PMC free article] [PubMed]
13. Chanock SJ, Hunter DJ. Genomics: when the smoke clears. Nature. 2008;452(7187):537–538. [PubMed]
14. Spitz MR, Amos CI, Dong Q, Lin J, Wu X. The CHRNA5-A3 region on chromosome 15q24-25.1 is a risk factor both for nicotine dependence and for lung cancer. J Natl Cancer Inst. 2008;100(21):1552–1556. [PMC free article] [PubMed]
15. Berrettini W, Yuan X, Tozzi F, et al. Alpha-5/alpha-3 nicotinic receptor subunit alleles increase risk for heavy smoking. Mol Psychiatry. 2008;13(4):368–373. [PMC free article] [PubMed]
16. Schuller HM. Cell type specific, receptor-mediated modulation of growth kinetics in human lung cancer cell lines by nicotine and tobacco-related nitrosamines. Biochem Pharmacol. 1989;38(20):3439–3442. [PubMed]
17. Wang Y, Pereira EF, Maus AD, et al. Human bronchial epithelial and endothelial cells express alpha7 nicotinic acetylcholine receptors. Mol Pharmacol. 2001;60(6):1201–1209. [PubMed]
18. Schuller HM, Orloff M. Tobacco-specific carcinogenic nitrosamines. Ligands for nicotinic acetylcholine receptors in human lung cancer cells. Biochem Pharmacol. 1998;55(9):1377–1384. [PubMed]
19. Cole SR, Hernan MA. Fallibility in estimating direct effects. Int J Epidemiol. 2002;31(1):163–165. [PubMed]
20. Blakely T. Commentary: estimating direct and indirect effects-fallible in theory, but in the real world? Int J Epidemiol. 2002;31(1):166–167. [PubMed]
21. Yokoyama A, Kato H, Yokoyama T, et al. Genetic polymorphisms of alcohol and aldehyde dehydrogenases and glutathione S-transferase M1 and drinking, smoking, and diet in Japanese men with esophageal squamous cell carcinoma. Carcinogenesis. 2002;23(11):1851–1859. [PubMed]
22. Weinberg CR. Can DAGs clarify effect modification? Epidemiology. 2007;18(5):569–572. [PMC free article] [PubMed]
23. Ottman R. Gene-environment interaction: definitions and study designs. Prev Med. 1996;25(6):764–770. [PMC free article] [PubMed]
24. Khoury MJ, Beatty TH, Cohen BH. Study of Genetic Factors in Disease. Oxford: Oxford University Press; 1993.
25. Taylor JM, Wang Y, Thiebaut R. Counterfactual links to the proportion of treatment effect explained by a surrogate marker. Biometrics. 2005;61(4):1102–1111. [PubMed]
26. Ditlevsen S, Christensen U, Lynch J, Damsgaard MT, Keiding N. The mediation proportion: a structural equation approach for estimating the proportion of exposure effect on outcome explained by an intermediate variable. Epidemiology. 2005;16(1):114–120. [PubMed]
27. Robins JM, Hernan MA, Brumback B. Marginal structural models and causal inference in epidemiology. Epidemiology. 2000;11(5):550–560. [PubMed]
28. Petersen ML, Sinisi SE, van der Laan MJ. Estimation of direct causal effects. Epidemiology. 2006;17(3):276–284. [PubMed]
29. Schatzkin A, Freedman LS, Dorgan J, McShane LM, Schiffman MH, Dawsey SM. Surrogate end points in cancer research: a critique. Cancer Epidemiol Biomarkers Prev. 1996;5(12):947–953. [PubMed]
30. Prentice RL. Surrogate endpoints in clinical trials: definition and operational criteria. Stat Med. 1989;8(4):431–440. [PubMed]
31. Rosenbaum PR. From association to causation in observational studies: the role of tests of strongly ignorable treatment assignment. J Am Stat Assoc. 1984;79((385)):41–48.
32. Kaufman JS, Cooper RS. Seeking causal explanations in social epidemiology. Am J Epidemiol. 1999;150(2):113–120. [PubMed]
33. Greenland S, Pearl J, Robins JM. Causal diagrams for epidemiologic research. Epidemiology. 1999;10(1):37–48. [PubMed]
34. Grucza RA, Wang JC, Stitzel JA, et al. A risk allele for nicotine dependence in CHRNA5 is a protective allele for cocaine dependence. Biol Psychiatry. 2008 PMID: 18519132. [PMC free article] [PubMed]
35. Popkin MK. Exacerbation of recurrent depression as a result of treatment with varenicline. Am J Psychiatry. 2008;165(6):774. [PubMed]
36. Ioannidis JP, Contopoulos-Ioannidis DG, Rosenberg PS, et al. Effects of CCR5-delta32 and CCR2-64I alleles on disease progression of perinatally HIV-1-infected children: an international meta-analysis. AIDS. 2003;17(11):1631–1638. [PubMed]
37. Feigelson HS, Cox DG, Cann HM, et al. Haplotype analysis of the HSD17B1 gene and risk of breast cancer: a comprehensive approach to multicenter analyses of prospective cohort studies. Cancer Res. 2006;66(4):2468–2475. [PubMed]
38. Wang JC, Grucza R, Cruchaga C, et al. Genetic variation in the CHRNA5 gene affects mRNA levels and is associated with risk for alcohol dependence. Mol Psychiatry. 2008 doi:10.1038/mp.2008.42. [PubMed]
39. Pal P, Xi H, Sun G, et al. Tagging SNPs in the kallikrein genes 3 and 2 on 19q13 and their associations with prostate cancer in men of European origin. Hum Genet. 2007;122(3–4):251–259. [PubMed]
40. Thomas G, Jacobs KB, Yeager M, et al. Multiple loci identified in a genome-wide association study of prostate cancer. Nat Genet. 2008;40(3):310–315. [PubMed]
41. Eeles RA, Kote-Jarai Z, Giles GG, et al. Multiple newly identified loci associated with prostate cancer susceptibility. Nat Genet. 2008;40(3):316–321. [PubMed]
42. Ahn J, Berndt S, Wacholder S, Kraft P, Kibel AS, Yeager M. Variation in KLK genes, prostate-specific antigen and risk of prostate cancer. Nat Genet. 2008;40((9)):1032–1034. [PMC free article] [PubMed]
43. Eeles R, Giles G, Neal D, Muir K, Easton DF. Reply to “Variation in KLK genes, prostate-specific antigen and risk of prostate cancer” Nat Genet. 2008;40((9)):1035–1036. [PMC free article] [PubMed]
Articles from JNCI Journal of the National Cancer Institute are provided here courtesy of
Oxford University Press