Published OnlineFirst August 25, 2009; DOI: 10.1158/0008-5472.CAN-09-0786 Published Online First on August 25, 2009 as 10.1158/0008-5472.CAN-09-0786
Epidemiology
The CHRNA5-CHRNA3-CHRNB4 Nicotinic Receptor Subunit Gene Cluster Affects Risk for Nicotine Dependence in African-Americans and in European-Americans
Nancy L. Saccone,1 Jen C. Wang,2 Naomi Breslau,5 Eric O. Johnson,6 Dorothy Hatsukami,7 Scott F. Saccone,2 Richard A. Grucza,2 Lingwei Sun,2 Weimin Duan,1 John Budde,2 Robert C. Culverhouse,3 Louis Fox,2 Anthony L. Hinrichs,2 Joseph Henry Steinbach,4 Meng Wu,1 John P. Rice,1,2 Alison M. Goate,1,2 and Laura J. Bierut2
Departments of 1Genetics, 2Psychiatry, 3Medicine, and 4Anesthesiology, Washington University, St.Louis, Missouri; 5Department of Epidemiology, Michigan State University, East Lansing, Michigan; 6Research Triangle Institute International, Research Triangle Park, North Carolina; and 7Department of Psychiatry, University of Minnesota, Minneapolis, Minnesota
Abstract Introduction Genetic association studies have shown the importance of Cigarette smoking is a major public health problem and is the variants in the CHRNA5-CHRNA3-CHRNB4 cholinergic nicotin- single largest cause of preventable death in the world, contributing ic receptor subunit gene cluster on chromosome 15q24-25.1 to >5 million deaths a year (1).In the United States, smoking causes for the risk of nicotine dependence, smoking, and lung cancer f443,000 deaths each year including an estimated 82% of lung in populations of European descent. We have carried out a cancer deaths (2). detailed study of this region using dense genotyping in both Cigarette smoking is common in European-Americans and European-Americans and African-Americans. We genotyped African-Americans, with a similar prevalence of current smoking 75 known single nucleotide polymorphisms (SNPs) and one among adults: 21.4% [95% confidence interval (CI), 20.4–22.4] in sequencing-discovered SNP in an African-American sample European-Americans and 19.8% (19.2–21.4) in African-Americans (N = 710) and in a European-American sample (N = 2,062). (3).In past-month smokers, the rates of current nicotine Cases were nicotine-dependent and controls were nondepen- dependence, defined by the Fagerstro¨m Test for Nicotine dent smokers. The nonsynonymous CHRNA5 SNP rs16969968 Dependence (FTND; ref.4), were also similar although slightly is the most significant SNP associated with nicotine depen- lower in African-Americans: 57.3% (55.1–59.5) versus 61.7% (61.0– dence in the full sample of 2,772 subjects [P = 4.49 10 8; 62.4) in European-Americans (5). Despite these similarities, odds ratio (OR), 1.42; 95% confidence interval (CI), 1.25–1.61] African-Americans have a higher incidence of lung cancer than as well as in African-Americans only (P = 0.015; OR, 2.04; European-Americans (82.7 per 100,000 versus 64.3 per 100,000) and 7 1.15–3.62) and in European-Americans only (P = 4.14 10 ; higher lung cancer mortality (64.1 per 100,000 versus 54.1 per OR, 1.40; 1.23–1.59). Other SNPs that have been shown to 100,000; ref.6).Many issues could contribute to these differences affect the mRNA levels of CHRNA5 in European-Americans are in lung cancer, including discrepancies in quit rates and disparities associated with nicotine dependence in African-Americans in medical care (7, 8).Genetic factors such as differences in allele but not in European-Americans. The CHRNA3 SNP rs578776, frequencies and linkage disequilibrium (LD) might also play a role. which has a low correlation with rs16969968, is associated We previously carried out a high-density whole-genome and with nicotine dependence in European-Americans but not in candidate gene association study of nicotine dependence in African-Americans. Less common SNPs (frequency V 5%) are European-descent samples (9, 10).Among the top association also associated with nicotine dependence. In summary, signals were chromosome 15q24-25.1 variants, including a multiple variants in this gene cluster contribute to nicotine biologically compelling finding at rs16969968, a nonsynonymous dependence risk, and some are also associated with functional single nucleotide polymorphism (SNP) in the a5 cholinergic effects on CHRNA5. The nonsynonymous SNP rs16969968, a nicotinic receptor subunit gene (CHRNA5; ref.10).The risk allele known risk variant in populations of European-descent, is also at rs16969968 causes an amino acid change (D398N, aspartate to significantly associated with risk in African-Americans. asparagine) in the a5 subunit and functional changes in the Additional SNPs contribute to risk in distinct ways in these receptor in vitro (11).Independent replication of this association two populations. [Cancer Res 2009;69(17):6848–56] with nicotine dependence or smoking has now been confirmed in multiple populations of European descent, using rs16969968 or SNPs highly correlated (r2 z 0.8) with rs16969968 (11–19). Fine- mapping in our original sample continued to highlight rs16969968 and other highly correlated SNPs, which extend across CHRNA5- CHRNA3-CHRNB4 and into neighboring genes IREB2, PSMA4, and LOC123688 (20). Note: Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org/). Also in this region, a locus tagged by rs578776 and rs3743078 is Requests for reprints: Nancy L.Saccone, Department of Genetics, Washington significantly associated with nicotine dependence and smoking in University School of Medicine, Box 8232, St.Louis, MO 63110.Phone: 314-747-3263; European-descent samples (10, 12, 14, 16).This locus has been Fax: 314-747-2489; E-mail: [email protected]. I2009 American Association for Cancer Research. described as protective because the minor allele is more frequent doi:10.1158/0008-5472.CAN-09-0786 in controls than in cases.There is evidence that this locus is
Cancer Res 2009; 69: (17). September 1, 2009 6848 www.aacrjournals.org
Downloaded from cancerres.aacrjournals.org on October 1, 2021. © 2009 American Association for Cancer Research. Published OnlineFirst August 25, 2009; DOI: 10.1158/0008-5472.CAN-09-0786
CHRNA5-CHRNA3-CHRNB4 and Nicotine Dependence in African-Americans statistically distinct from rs16969968: in two independent European- were dropped.All remaining DNA samples had call rates above 90% and ancestry samples, the correlation between rs578776 and rs16969968 were retained; 99.5% of DNA samples had call rates z95%.Twenty samples is low (r2 =0.2),withD¶ = 1.0 indicating repulsion phase (the minor were included in duplicate.SNPs with more than one nonconcordant allele at one locus co-occurs with the major allele at the other), and duplicate sample genotype were dropped; the 33 SNPs reported here were 100% concordant.SNPs were also required to pass a call rate threshold of both SNPs are significant in joint analyses (14, 20). 98%.Self-reported race was verified using EIGENSTRAT (27).An additional Genome-wide studies of lung cancer report strong association 9 2 42 SNPs were genotyped using Sequenom MassArray iPLEX technology. with rs16969968 and its r correlates (16, 21–23).The locus tagged These SNPs passed a call rate threshold of 95%.Map positions were by rs578776 is also associated with lung cancer.Thus, there has obtained from the National Center for Biotechnology Information Human been a convergence of genetic findings for nicotine dependence Reference Build 36.2 and dbSNP Build 129. and for lung cancer. These 75 SNPs were selected to include the same SNPs used to cover this The association studies cited above were carried out in region in our previous publications (9, 10, 20).Forty-nine of the genotyped European-descent samples.Studying this region in individuals of SNPs lie between 76642961 bp and 76722642 bp, encompassing CHRNA5- ¶ ¶ African descent provides an opportunity to confirm its role in CHRNA3-CHRNB4 and the 2 kb flanking either side (5 of CHRNA5 and 5 of an important and understudied population.Also, the contrasting CHRNB4), matching the dbSNP definition of in or near a gene.These SNPs provide good coverage of the 72 SNPs in this region that are common genetic architecture in Africans and Europeans can be leveraged [minor allele frequency (MAF) z0.05] in either the HapMap CEU or YRI to narrow down the potential functional source of the disease samples.Of the 61 SNPs that have MAF z5% in CEU, 54 (93%) are tagged by associations (24).In the European ancestry population, many at least one genotyped SNP at an r2 z 0.8 according to CEU LD data; all 61 highly correlated SNPs could often detect the same association are tagged at an r2 z 0.5. Of the 65 SNPs in this region that have MAF z5% signal, and determining which variant(s) are ‘‘causal’’ is not in YRI, 46 (71%) are tagged at an r2 z 0.8 in the YRI LD data; 60 (92%) are straightforward.LD and allele frequencies can differ between tagged at an r2 z 0.5 (HapMap release 23a). Twenty-six additional SNPs European and African populations.For instance, in HapMap, flanking the cluster or in IREB2, LOC123688, and PSMA4 were included 2 rs16969968 is nonpolymorphic in the YRI (Yoruba in Ibadan, because of the high r extending from CHRNA5-CHRNA3-CHRNB4. Nigeria) sample; for rs3743078, the minor allele in YRI (‘‘C’’, 33%) is Genotype data were then combined for this new sample and the the major allele (77%) in CEU (CEPH, Utah residents with ancestry previously genotyped sample of 1,608 European-Americans.Allele corre- spondence was checked.All association analyses are reported using the from northern and western Europe; ref.25). major allele in the European-American sample as the reference allele. Here, we present an association study of nicotine dependence DNA resequencing in the CHRNA5 and CHRNA3 intergenic region. across the CHRNA5-CHRNA3-CHRNB4 region in a new African- Human CHRNA3 and CHRNA5 are physically linked and partially overlap American sample and an expanded European-American sample. tail-to-tail (28).We targeted a 26-kb region spanning exon 5 of CHRNA5 to We also analyze novel variants discovered from DNA resequencing exon 4 of CHRNA3 and sequenced 10 cases and 12 controls of European of European-Americans in a 26 kb region between CHRNA5 and descent.The goal was to identify novel variants, genotype them in the larger CHRNA3.Our goal is to determine whether any variants sample, and test them for association with nicotine dependence. consistently affect nicotine dependence risk in both European- Sequencing was performed at the Genome Center at Washington University 10 Americans and African-Americans.Additionally, we are interested using standard protocols. We used Sequencher (Gene Codes) to assemble in evidence of variants that affect risk in the African-American and scan sequence traces for variation from the reference sequence.Novel variants were genotyped using Sequenom MassArray iPLEX technology. sample alone. Linkage disequilibrium. LD was calculated in European-Americans and African-Americans using PLINK (29) and Haploview (30).LD ( r2) plots were 11 Materials and Methods generated using WGAviewer (31). Genetic association analyses. Association analysis of case-control Study design and sample. All individuals (N = 2,772) were recruited by status was performed on the entire sample using logistic regression with the Collaborative Genetic Study of Nicotine Dependence (COGEND), a gender and race as covariates.This analysis highlights SNPs that have a multi-site project from the United States.All cases and controls reported similar effect in both European-Americans and African-Americans.The full smoking z100 cigarettes lifetime, the threshold classically used to define a genetic model is smoker (26).In the recruitment of new subjects, cases were nicotine- z dependent, defined by current FTND 4, and controls were required to P have a lifetime FTND of 0 or 1.The African-American genetic sample in ln ¼ a þ b g þ b s þ b G; ðAÞ 1 P 1 2 3 COGEND consists of 461 cases and 249 controls.The European-American genetic sample of 1,063 cases and 999 controls consists of 454 new subjects (266 cases and 188 controls) together with 1,608 original COGEND where P is the probability of being a case, g is gender (0 = male, 1 = female), European-American individuals previously analyzed in refs.[9, 10; 797 s is race (0 = European-American, 1 = African-American), and G is genotype, cases (FTND z4) and 811 controls (FTND = 0)]. coded as the number of copies of the minor allele in the European- The study complies with the Code of Ethics of the World Medical American sample, corresponding to a log additive (multiplicative) model. 2 Association and obtained informed consent from participants and approval The term b3G is tested for significance by the standard likelihood m with from the appropriate institutional review boards. 1 df.We then analyzed African-Americans and European-Americans Genotyping and quality control. Thirty-three SNPs covering the separately with gender as a covariate.Association tests were carried out CHRNA5-CHRNA3-CHRNB4 region were genotyped in the additional using PLINK (29) and SAS. subjects by the Genome Center at Washington University.8 These SNPs Cross-population heterogeneity analysis. Heterogeneity of the odds were genotyped using Illumina GoldenGate assay as part of a set of 1,536 ratios (OR) in European-Americans and African-Americans was tested SNPs that included SNPs for checking population stratification and SNPs to by adding the interaction term b4Gs to Eq.A (24).This test is of interest follow-up several genomic regions.Cleaning was performed using all 1,536 for SNPs that are associated with disease in at least one population.If SNPs.Four DNA samples had poor call rates ( V63% of the 1,536 SNPs) and
9 www.sequenom.com 10 http://genome.wustl.edu/activity/med_seq/protocols.cgi 8 http://genome.wustl.edu/services/microarray.cgi 11 http://people.genome.duke.edu/~dg48/WGAViewer/ www.aacrjournals.org 6849 Cancer Res 2009; 69: (17). September 1, 2009
Downloaded from cancerres.aacrjournals.org on October 1, 2021. © 2009 American Association for Cancer Research. Published OnlineFirst August 25, 2009; DOI: 10.1158/0008-5472.CAN-09-0786
Cancer Research
Table 1. Selected genotyped SNPs and allele frequencies, sorted by position
SNP Chr Position (bp) Gene Reference Minor allele Reference allele frequency allele* European-American African-American European-American African-American sample sample sample sample
c rs17483548 15 76517368 IREB2 G A A 0.662 0.941 c rs17405217 15 76518204 IREB2 C T T 0.662 0.944 c rs2656052 15 76527987 IREB2 T G G 0.653 0.565 c rs17484235 15 76548469 IREB2 C G G 0.662 0.944 c rs7164594 15 76590112 LOC123688 C T T 0.785 0.584 c rs8034191 15 76593078 LOC123688 T C C 0.648 0.840 c rs10519203 15 76601101 LOC123688 A G G 0.643 0.697 rs2036534 15 76614003 LOC123688 T C C 0.789 0.778 c rs12916483 15 76619452 PSMA4 G A A 0.575 0.835 c rs3813570 15 76619887 PSMA4 A G G 0.789 0.739 c rs11858230 15 76622607 PSMA4 G A A 0.574 0.771 c rs2036527 15 76638670 C T T 0.645 0.778 c rs667282 15 76650527 CHRNA5 T C C 0.779 0.704 c rs588765 15 76652480 CHRNA5 C T T 0.576 0.705 c rs6495306 15 76652948 CHRNA5 A G G 0.576 0.702 c rs17486278 15 76654537 CHRNA5 A C C 0.648 0.710 rs601079 15 76656634 CHRNA5 T A A 0.574 0.578 rs680244 15 76658343 CHRNA5 G A A 0.573 0.576 c rs637137 15 76661031 CHRNA5 T A A 0.779 0.701 c rs951266 15 76665596 CHRNA5 C T T 0.649 0.888 c rs555018 15 76666297 CHRNA5 T C C 0.576 0.699 c rs647041 15 76667536 CHRNA5 C T T 0.578 0.704 c rs16969968 15 76669980 CHRNA5 G A A 0.650 0.948 c rs660652 15 76674887 CHRNA3 G A A 0.632 0.717 b c rs578776 15 76675455 CHRNA3 CT C 0.726 0.451 rs6495307 15 76677376 CHRNA3 C T T 0.580 0.595 c rs1051730 15 76681394 CHRNA3 C T T 0.651 0.882 b c rs3743078 15 76681814 CHRNA3 CG C 0.774 0.416 c rs3743077 15 76681951 CHRNA3 G A A 0.579 0.847 c rs1317286 15 76683184 CHRNA3 A G G 0.645 0.753 c rs938682 15 76683602 CHRNA3 T C C 0.778 0.711 c rs11637630 15 76686774 CHRNA3 A G G 0.778 0.697 c GSC3_26 15 76693692 CHRNA3 A G G 0.938 0.986 c rs7177514 15 76694461 CHRNA3 C G G 0.777 0.698 c rs6495308 15 76694711 CHRNA3 T C C 0.777 0.702 c rs8042059 15 76694914 CHRNA3 A C C 0.777 0.702 rs8042374 15 76695087 CHRNA3 A G G 0.780 0.750 c rs4887069 15 76696125 CHRNA3 A G G 0.772 0.698 c rs8192475 15 76698285 CHRNA3 G A A 0.950 0.991 c rs6495309 15 76702300 CHRNB4 C T T 0.794 0.748 c rs12914008 15 76710560 CHRNB4 G A A 0.955 0.994 c rs17487223 15 76711042 CHRNB4 C T T 0.629 0.897 c rs1996371 15 76743861 A G G 0.587 0.929
*The reference allele is the major allele in the European-American sample. cAllele frequencies differ significantly between European-Americans and African-Americans by Fisher’s exact test (P <10 3). bMinor allele differs in European-Americans and African-Americans. the interaction term is significant (e.g., P V 0.05, unadjusted for multiple Results tests), we conclude that the genetic effect is different in the two population samples. Genetic Association Analyses A significant difference in ORs may occur because the genetic effect is Supplementary Table S1 shows the allele frequencies of the 75 present in only one of the two populations.In this case, one interpretation known SNPs, and one novel SNP from our SNP discovery project, in is that this SNP is not biologically causal, although in the associated the European-American and African-American samples; selected population it may tag the functional SNP.This interpretation is based on SNPs, including those most associated with nicotine dependence the hypothesis that the underlying biological processes leading to disease in each sample, are in Table 1.Sixty-three SNPs show significant risk are shared in common across human populations (32). allele frequency differences between European-Americans and
Cancer Res 2009; 69: (17). September 1, 2009 6850 www.aacrjournals.org
Downloaded from cancerres.aacrjournals.org on October 1, 2021. © 2009 American Association for Cancer Research. Published OnlineFirst August 25, 2009; DOI: 10.1158/0008-5472.CAN-09-0786
CHRNA5-CHRNA3-CHRNB4 and Nicotine Dependence in African-Americans