Associations of Variants in CHRNA5/A3/B4 Gene Cluster with Smoking Behaviors in a Korean Population

Abstract

Multiple genome-wide and targeted association studies reveal a significant association of variants in the CHRNA5-CHRNA3-CHRNB4 (CHRNA5/A3/B4) gene cluster on chromosome 15 with nicotine dependence. The subjects examined in most of these studies had a European origin. However, considering the distinct linkage disequilibrium patterns in European and other ethnic populations, it would be of tremendous interest to determine whether such associations could be replicated in populations of other ethnicities, such as Asians. In this study, we performed comprehensive association and interaction analyses for 32 single-nucleotide polymorphisms (SNPs) in CHRNA5/A3/B4 with smoking initiation (SI), smoking quantity (SQ), and smoking cessation (SC) in a Korean sample (N = 8,842). We found nominally significant associations of 7 SNPs with at least one smoking-related phenotype in the total sample (SI: P = 0.015∼0.023; SQ: P = 0.008∼0.028; SC: P = 0.018∼0.047) and the male sample (SI: P = 0.001∼0.023; SQ: P = 0.001∼0.046; SC: P = 0.01). A spectrum of haplotypes formed by three consecutive SNPs located between rs16969948 in CHRNA5 and rs6495316 in the intergenic region downstream from the 5′ end of CHRNB4 was associated with these three smoking-related phenotypes in both the total and the male sample. Notably, associations of these variants and haplotypes with SC appear to be much weaker than those with SI and SQ. In addition, we performed an interaction analysis of SNPs within the cluster using the generalized multifactor dimensionality reduction method and found a significant interaction of SNPs rs7163730 in LOC123688, rs6495308 in CHRNA3, and rs7166158, rs8043123, and rs11072793 in the intergenic region downstream from the 5′ end of CHRNB4 to be influencing SI in the male sample. Considering that fewer than 5% of the female participants were smokers, we did not perform any analysis on female subjects specifically. Together, our detected associations of variants in the CHRNA5/A3/B4 cluster with SI, SQ, and SC in the Korean smoker samples provide strong evidence for the contribution of this cluster to the etiology of SI, ND, and SC in this Asian population.

About the Authors:

Ming D. Li

* E-mail: [email protected] (MDL); [email protected] (TP)

Affiliation: Department of Psychiatry and Neurobehavioral Sciences, University of Virginia, Charlottesville, Virginia, United States of America

Dankyu Yoon

Affiliation: Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul, Korea

Jong-Young Lee

Affiliation: Center for Genome Science, National Institute of Health, Seoul, Korea

Bok-Ghee Han

Affiliation: Center for Genome Science, National Institute of Health, Seoul, Korea

Tianhua Niu

Affiliation: Department of Psychiatry and Neurobehavioral Sciences, University of Virginia, Charlottesville, Virginia, United States of America

Thomas J. Payne

Affiliation: Department of Otolaryngology and Communicative Sciences, University of Mississippi Medical Center, Jackson, Mississippi, United States of America

Jennie Z. Ma

Affiliation: Department of Public Health Sciences, University of Virginia, Charlottesville, Virginia, United States of America

Taesung Park

* E-mail: [email protected] (MDL); [email protected] (TP)

Affiliation: Department of Statistics, College of Natural Science, Seoul National University, Seoul, Korea

Introduction

There are about 1.3 billion smokers worldwide, and the mortality burden from tobacco use has been estimated to exceed 6 million annually [1]. In Korea, male smoking prevalence is among the highest in the world [2]. In 2000, the prevalence was estimated to be 68% among men and 3% among women [2]. The number of deaths attributable to smoking-related diseases in Korea is about 35,000 each year, and the economic loss from premature death from smoking-related diseases exceeded 3 trillion won (approximately US $2.5 billion) in the year 2000 [3].

Nicotine is the psychoactive substance in tobacco that causes addiction. Many of those who want to quit smoking do not seek treatment but are unable to quit on their own [4]. Evidence for moderate heritabilities of smoking initiation (SI), nicotine dependence (ND), and smoking cessation (SC) has led to intensive efforts to identify susceptibility loci for these complex behavioral phenotypes [5], [6], [7], [8].

The psychopharmacologic effects of nicotine are mediated primarily by functionally diverse neuronal nicotinic acetylcholine receptors (nAChRs), a family of ligand-gated ion channels widely distributed in the brain [9], [10], [11]. Because of their unique functions, genes encoding various nAChR subunits have been proposed as plausible candidates for genetic studies of ND. Several subunit genes have been investigated for associations with ND as well as other smoking-related behaviors in human subjects (for reviews, see [12], [13]).

Initially, Saccone et al. [14] reported associations of multiple single nucleotide polymorphisms (SNPs) in the CHRNA5/A3/B4 gene cluster with ND, with the smallest P value of 0.00064 for rs16969968 (D398N) in exon 5 of CHRNA5. However, the significance of these results did not survive correction for multiple testing. Shortly after the initial publication, rs1051730 (Y215Y) in exon 5 of CHRNA3 was reported to be significantly associated at a genome-wide level with smoking quantity (SQ) [15] and lung cancer (LC) [15], [16], [17]. Several other genome-wide and candidate gene-based association studies provided further evidence for the association of variants of the CHRNA5/A3/B4 cluster with various nicotine-related behaviors [14], [18], [19], [20], [21], [22], [23], [24], [25]. In contrast, other studies have failed to reveal a significant association of this gene cluster with ND or other smoking-related phenotypes [26], [27].

Considering that (1) the participants in most of these studies were primarily European Americans or of European origin, with the exception of two studies on African Americans [18], [28], and (2) there are distinct differences in linkage disequilibrium (LD) patterns across different ethnic populations [12], it is of tremendous interest to determine whether variants in this gene cluster play any role in the etiologies of smoking behavioral phenotypes in other ethnic groups. Thus, the major objective of this study was to test for such genetic effects in a large population-based Korean sample.

Results

Description of KARE sample in its relation to smoking behaviors

Of the 8,842 subjects, 4,205 were recruited from Ansung and 4,637 from Ansan, Korea. Their average ages were 55.6±8.74 (standard deviation; SD) and 49.1±7.86 years, respectively. Although 52.7% of the participants (N = 4,659) were female, only 4.93% of these were considered either former (1.34%), light (1.29%), or habitual (2.30%) smokers. In contrast, 80.62% of the male subjects were smokers, with 31.05% being former smokers, 4.78% light smokers, and 44.79% habitual smokers. For those habitual smokers, the average number of cigarettes smoked per day (CPD) was 19.51±8.74 for male and 11.93±7.28 for female smokers. A detailed description of the characteristics of all subjects is presented in Table 1.

[Figure omitted. See PDF.]

Table 1. Demographic characteristics of study subjects.

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

Individual SNP-based association analysis

Among the 36 SNPs genotyped for the 15q24-15q25.1 region in CHRNA5/A3/B4, only 32 had a minor allele frequency (MAF) of >0.01. Considering the differences in LD patterns of the region across multiple ethnic samples and to have a better understanding of the LD landscape within this region in the Korean vs. the other ethnic populations, Table S1 provides a detailed list of all genotyped SNPs for the total and male samples. However, only those 32 SNPs with MAFs >0.01 were used in the association analyses for the three smoking-related phenotypes.

Associations of individual SNPs with the three phenotypes were determined with the PLINK program [29], and the results are shown in Table 2. Altogether, we found 7 SNPs that had nominally significant associations with at least one smoking-related phenotype in either the total or the male sample. In the total sample, we found that the A allele of rs951266 in CHRNA5 was nominally significantly associated with SI (P = 0.023; odds ratio [OR]  = 1.32; 95% confidence interval [CI]: = 1.04, 1.67) and SQ (P = 0.008; OR = 1.48; 95% CI: 1.11, 1.98), and the G allele of rs11072768 in CHRNB4 was nominally significantly associated with SI (P = 0.015; OR = 1.14; 95% CI: 1.03, 1.27), SQ (P = 0.028; OR = 1.17; 95% CI: 1.02, 1.34), and SC (P = 0.018; OR = 1.16; 95% CI: 1.03, 1.31). Furthermore, the associations of SNPs rs8043123 (C), rs4887077 (T), and rs11638372 (T) with SQ and the association of rs2869550 (C) with SC reached nominal significance.

[Figure omitted. See PDF.]

Table 2. P values for SNPs significantly associated with at least one smoking-related phenotype and their corresponding odds ratios and 95% confidence intervals under the additive and dominant model.

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

In the male sample, the G allele of rs11072768 in CHRNB4 was nominally significantly associated with SI (P = 0.001; OR = 1.22; 95% CI: 1.08, 1.37), SQ (P = 0.016; OR = 1.16; 95% CI: 1.03, 1.31), and SC (P = 0.01; OR = 1.18; 95% CI: 1.04, 1.34), and rs951266 (A) in CHRNA5 and rs6495308 (T) in CHRNA3 were nominally significantly associated with SI and SQ, respectively. In addition, we found that SNPs rs4887077 (T) and rs11638372 (T) in the intergenic region downstream from the 5′ end of CHRNB4 were nominally significantly associated with SQ. Of these associations, only that of rs951266 in CHRNA5 with SQ and rs11072768 in CHRNA4 with SI remained significant after correction for multiple testing.

We did not perform the association analysis on the female sample because fewer than 5% of these subjects smoked, such that the sample was too small to derive any meaningful conclusions.

Haplotype block structure and LD analysis

The pair-wise D' values of 32 SNPs within CHRNA5/A3/B4 were determined using the Haploview program [30]. On the basis of the block definition proposed by Gabriel et al. [31], we found three discernible haplotype blocks within the cluster in the total sample (Figure 1). The first block, with a size of about 10 kb, contains four SNPs between rs169669920 and rs7163730; the second block, with a size of 46 kb, contains 12 SNPs; the third block, with a size of 23 kb, contains only two SNPs — rs11072793 and rs11072794. Considering that (1) no haplotype block was detected within a genomic region of about 108 kb between rs481134 and rs6495316, and (2) the three discernible haplotype blocks did not correspond to the defined genes within this region, we decided to perform haplotype-based association analysis using a sliding window approach with a window size of three consecutive SNPs [32], as described in the following section.

[Figure omitted. See PDF.]

Figure 1. Haploview-generated LD patterns for 32 SNPs within the CHRNA5/A3/B4 cluster in the Korean total sample.

Pair-wise LD between all SNPs was evaluated using the Haploview program [30] with the option of determining haplotype blocks according to the criteria defined by Gabriel et al. [31]. The number in each box represents the D' value for each SNP pair.

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

Haplotype-based association analysis

Using the Haplo.stats R statistics package [33], we performed haplotype-based association analysis for the three smoking-related phenotypes in both the total and male samples (Tables 3 and 4). In the total sample, we found five major haplotypes between SNPs rs16969948 in CHRNA5 and rs6495316 in the intergenic region downstream from the 5′ end of CHRNB4 that were significantly associated with SI under the additive or dominant model. Of these haplotypes, two were positively, and the remaining three were negatively, associated with SI, with P values ranging from 0.013 for the haplotype ACT, formed by rs950776-rs11072768-rs7166158, to 0.047, for the haplotype TCC, formed by rs16969948-rs481134-rs951266.

[Figure omitted. See PDF.]

Table 3. Statistics estimates and P values for major haplotypes showing significant associations at the 0.05 level with at least one smoking-related phenotype under the additive and dominant models in the total sample.

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

[Figure omitted. See PDF.]

Table 4. Statistics estimates and P values for major haplotypes showing significant associations at the 0.05 level with at least one smoking-related phenotype under the additive and dominant models in the male sample.

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

For SQ, we found seven major haplotypes between SNPs rs16969948 in CHRNA5 and rs6495316 in the intergenic region downstream from the 5′ end of CHRNB4 that showed significant associations. Of these haplotypes, two were positively, and the remaining five were negatively, associated with SI, with a P value ranging from 0.014 for the haplotype CCT, formed by rs481134-rs951266-rs514743, to 0.044 for the haplotype GGA, formed by rs4887077-rs11638372-rs6495316. Although we performed identical association analyses for SC, we found no significant haplotypes in the total sample. By comparing significant haplotypes with SI and SQ, we found four haplotypes, GAA (rs6495308-rs950776-rs11072768), ACT (rs950776-rs11072768-rs7166158), CTC (rs11072768-rs7166158-rs8043123), and GGA (rs4887077-rs11638372-rs6495316), to be significantly associated with both phenotypes under identical genetic models, indicating that they contribute to the etiologies of these two behaviors, likely through a common uncharacterized mechanism.

Similarly, we performed haplotype-based association analysis for the male sample (Table 4). Altogether, we found that ten haplotypes between SNPs rs16969948 in CHRNA5 and rs11638372 in the intergenic region downstream from the 5′ end of CHRNB4 showed significant associations with at least one smoking-related phenotype. By comparing the haplotype-based association results of the total and male samples, we found six that were comparable in the two samples, except that the association signals appear to be much stronger in the male sample. For example, the P values of associations with SI, SQ, and SC for the haplotype ACT, formed by rs950776-rs11072768-rs7166158 under the dominant model, are, respectively, 0.001, 0.014, and 0.019 in the male sample and 0.013, 0.033, and 0.06 in the total sample. Very interestingly, we found that haplotypes ACT, formed by rs950776-rs11072768-rs7166158, and CTC, formed by rs11072768-rs7166158-rs8043123, were nominally significantly associated with all three smoking-related phenotypes in the male sample, suggesting that this region is a good target for further sequencing analysis with the hope of identifying variants that contribute to our observed association of the cluster with these smoking-related phenotypes.

Interaction analysis of CHRNA5/A3/B4 with SI, SQ, and SC

To further determine a genetic contribution of the CHRNA5/A3/B4 cluster to the three smoking-related phenotypes, we performed a gene-gene interaction analysis for the total and male samples as we did previously for nAChR genes alpha 4 (CHRNA4) and beta 2 (CHRNB2) [34] and GABA-B receptor subunits 1 (GABAB1) and 2 (GABAB2) [35]. By using the generalized multifactor dimensionality reduction (GMDR) approach [36], we performed the interaction analysis for all possible combinations of two to six SNPs within CHRNA5/A3/B4 for SI and SC in both the total and the male sample. Among these gene-gene interaction analyses, a significant result emerged only in the male sample for SI, with the best interaction model consisting of SNPs rs7163730 in LOC123688, rs6495308 in CHRNA3 and rs7166158, rs8043123, and rs11072793 from the downstream region of CHRNB4; the corresponding P value was 0.011 (Table 5). Because of the limitation of the current GMDR version, which cannot handle ordinal phenotypes such as SQ, we could not carry out gene-gene interaction analysis on this particular phenotype in our samples. However, we did perform interaction analysis on SQ by employing a cumulative logit model, which revealed no significant interactions of SNPs located in this genomic region in affecting SQ.

[Figure omitted. See PDF.]

Table 5. A significant interactive model for CHRNA5/A3/B4 cluster with SI in KARE male sample on the basis of the prediction accuracy and the sign test P value.

https://doi.org/10.1371/journal.pone.0012183.t005

Discussion

In the current study, we examined genetic associations and epistatic variants in CHRNA5/A3/B4 with SI, SQ, and SC in the total and male samples. Individual SNP-based association analyses revealed that seven SNPs within this region showed significant associations with at least one smoking behavior in either or both samples. Of these polymorphisms, rs951266 in CHRNA5 showed the strongest association with SI (P = 0.006) and SQ (P = 0.001) and rs11072768 in CHRNB4 with all the three smoking-related phenotypes in the total (SI: P = 0.015; SQ: P = 0.028; SC: P = 0.018) and male (SI: P = 0.001; SQ: P = 0.016; SC: P = 0.010) samples. Furthermore, we found multiple haplotypes between rs16969948 in CHRNA5 and rs6495316 or rs11638372 in the intergenic region downstream from the 5′ end of CHRNB4 that were associated significantly with the three smoking-related phenotypes in both the total and the male sample. Finally, considering the fact that the protein products of these nAChR subunit genes must assemble in order to form functional nAChRs, we also performed gene-gene interaction analysis on all SNPs in this genomic region and found the best interactive model involved five SNPs located in LOC123688, CHRNA3, and the intergenic region at the 5′ end of CHRNB4 that influences SI in the male sample.

To assess ND, various scales have been developed, which include the Diagnostic and Statistical Manual of Mental Disorders (DSM)-IV criteria [4], Fagerström Test for ND (FTND) [37], HSI [38], the Nicotine Dependence Syndrome Scale (NDSS) [39], and the Wisconsin Inventory of Smoking Dependence Motives (WISDM-68) [40], to name a few. Among them, FTND is the most widely used instrument, primarily because of its succinctness and ease of administration [41]. The FTND score is commonly treated as a continuous variable, and an FTND ≥4 is typically defined as highly dependent [28], [42]. CPD, the most commonly used consumption measure [38], is both highly heritable [7] and predictive of ND [24], [38]. SQ is an ordinal variable based on CPD, which provides a simple, quantified index of the amount of cigarette consumption. It should be noted that SQ and FTND are correlated but different psychometric measures of smoking behavior. To date, associations of variants in the CHRNA5/A3/B4 cluster with smoking have been detected for both SQ and FTND phenotypes, but signals appear to be stronger for SQ than for FTND (Figure 2), which is also true for linkage study on ND based on a review of more than 20 independent studies [43].

[Figure omitted. See PDF.]

Figure 2. A summary of all reported SNPs that have been significantly associated with different smoking behaviors.

Of these reported studies mentioned in the figure, most were investigated on nicotine dependence except that the study authored by Keskitalo et al. [55] was on serum cotinine level and that by Schlaepfer et al. [23] was on smoking initiation.

https://doi.org/10.1371/journal.pone.0012183.g002

There are several strengths of this study. First, although significant associations of variants within the CHRNA5/A3/B4 cluster with smoking have been replicated in multiple independent samples, almost all these studies were concentrated on smoking dependence in samples of European origin [14], [18], [19], [20], [21], [22], [23], [24], [25]. Very recently, relatively weak associations of some variants of this gene cluster with ND were reported in two African American samples [18], [28]. Among the 76 SNPs investigated by Saccone et al. [28], eight showed nominally significant association with ND (P value  = 0.0147∼0.0443), with rs16969968 in CHRNA5 being the strongest, at a P value of 0.0147. Similarly, in a prior study [18], we analyzed 22 SNPs within the gene cluster and found that only rs8040868 in CHRNA3 showed a nominally significant association with ND (P = 0.017∼0.039). However, none of these associations survived correction for multiple testing in the two African American samples [18], [28]. Although there were two studies investigating potential associations of variants in the gene cluster with LC and smoking in Chinese [44] and Japanese [45] subjects, no previous reports have examined the involvement of variants of this cluster with smoking behavior in a strictly Asian sample. Therefore, this study represents the first report on a potential association of the CHRNA5/A3/B4 cluster with smoking behaviors in an Asian population. By identifying risk alleles for smoking and LC and comparing their frequencies in different ethnic populations, it is easy to see that they are quite different, which further underscores that different variants in this genomic region contribute importantly to variations in disease risks among different ethnic populations [12].

Second, most previous association studies of the CHRNA5/A3/B4 cluster were focusing on ND or related measures [14], [15], [19], [24], [25], although a few studies were focused on SI [23], [27] or SC [26], [46], [47]. Given the significant overlap of genetic underpinnings among these three smoking-related phenotypes [48], [49], it is of great interest to investigate these phenotypes systemically in one sample. This is especially true for the CHRNA5/A3/B4 cluster, as it has been associated with ND in several reports, especially in samples of European origin [14], [18], [19], [20], [21], [22], [23], [24], [25]. However, association of the variants of this cluster with SI and SC has not been well established. By analyzing the three major smoking-related phenotypes together in the Korean sample, we found significant associations of some variants in this cluster with the three smoking-related phenotypes at both the individual SNP and the haplotype levels.

Third, in this study, we not only performed association analyses of this region with the three major smoking-related phenotypes at the SNP and haplotype levels, as in most reported association studies of the region with smoking, but also performed extensive gene-gene interaction analysis exclusively in this region. This is important in that: (1) all nAChR subunits except alpha 7 must assemble under appropriate compositions in order to form functional receptors; and (2) it is still undetermined whether other uncharacterized genes of this region contribute to smoking-related phenotypes [12]. Our gene-gene interaction results in the male sample indicate clearly that some variants exist within the region between rs7163730 in LOC123688 and rs11072793 in the intergenic region downstream from the 5′ end of CHRNB4 that are contributing to SI through gene-gene interactions. This appears consistent with findings from our recent study showing an association of the region with ND in European American and African American samples [18]. This has been the case for several other genes also. For example, we recently detected a significant interaction between variants of CHRNA4 and CHRNB2 in affecting ND, such that CHRNB2 contributes significantly to the etiology of ND together with CHRNA4 through gene-gene interaction [34]. Similarly, in another study, we found that GABAB1 contributes to ND through its interaction with GABAB2, although no significant association was observed for GABAB1 with ND directly [35].

There exist several limitations to this study. First, we did not consider Bonferroni correction for all possible multiple testing, primarily because we consider this study a replication of a positive association of this cluster with ND in an independent sample. A similar approach has been adopted by other researchers [25]. If we include correction for multiple testing, most detected associations of variants within this cluster with SI, SQ and SC at both the individual SNP and haplotype levels become non-significant, with the exceptions that rs951266 in CHRNA5 remained significant for SQ and rs11072768 in CHRNB4 for SI in the male sample (see Table 3). Similarly, we found that associations of four major haplotypes with SI and two major ones with SQ remained significant after Bonferroni correction for multiple testing of major haplotypes (Table 4). However, no significant associations of any SNPs or haplotypes with SC remained after Bonferroni correction for multiple testing. Also, we did not correct for our testing of two genetic models and three phenotypes, as they are highly related, which violates the assumption of independence for Bonferroni correction. Because of the aforementioned concerns, to some extent, we consider this study explorative, and more replication in Korean or other Asian samples is greatly needed. Second, although we defined the three smoking-related phenotypes on the basis of all the related information collected from each smoker, we are not fully convinced this is the best sample for investigation of this gene cluster and smoking behaviors, as this work was not designed originally as a genetic study on smoking, and the information collected from each smoker is limited. For example, SQ was assessed by CPD instead of ND determined by the FTND or DSM-IV criteria commonly used in other studies. Also, there was no SQ information (i.e., a 5-category ordinal trait based on CPD) for light or occasional smokers, although we do know they smoked only when they were in specific social circumstances such as during a party or gathering with their friends. Further, the SC status was based on individual self-reports, which were not biochemically verified, and this might lead to potential bias in phenotyping. Nevertheless, we do not feel this bias would greatly affect our results, as the definitions used in the current study have been used in other studies of this type as well [27], [43]. This is especially true for SQ, as it produces significantly more positive findings than any other smoking measures for ND in both linkage and association analyses [43]. Third, we investigated the association of this cluster with three smoking-related phenotypes only in the Korean total and male samples. We did not perform a similar analysis of the female sample because of the small sample size attributable to the small percentage of female smokers (3%), which is the case for many other Asian countries such as China (∼4% of Chinese women aged 15 years or older are smokers [50]). Thus, we could not determine whether this gene cluster has any significant impact on smoking behaviors in Korean women. Finally, the number of informative SNPs genotyped within this region was limited. To have a better coverage for these samples, more SNPs with a high density are greatly needed for further analysis of this cluster in relation to smoking-related phenotypes.

As shown in Table S1, the allele frequencies of an array of genetic variants (e.g., rs7168796, rs16969922, and rs1979906, to name a few) located in CHRNA5/A3/B4 differ dramatically across ethnic and geographical populations (i.e., the Caucasian [CEU], African [YRI], and Asian [CHB, JPT, and Korean] samples). Therefore, this study provides an essential basis for genetic association analysis by comparing different frequencies of SNPs of this region between smoking and non-smoking populations to examine whether such associations are consistent in different ethnicities. Indeed, different genotype frequencies have been compared in smoking and non-smoking populations in different ethnic groups for several polymorphisms of CYP2A6 [51] as well as a 44-bp insertion/deletion polymorphism of the serotonin transporter (5-HTT) gene [52], and such studies may shed new light on what molecular variants have diverging effects, and what molecular variants have converging effects in populations of different ancestries. This approach is termed “cross-population contrast mapping” based on the hypothesis that important biological mechanisms underlying a disease of interest are shared by various human populations, although differences in allele frequencies at risk loci could result in different prevalences in different populations [53]. Because smoking is a significant risk factor for LC, and the association of CHRNA5/A3/B4 with smoking or LC risk had been reported in many other populations, the present work should contribute to establishing that different polymorphic CHRNA5/A3/B4 patterns associate with the physiology of smoking and LC in ethnically and geographically different populations.

To compare the association signals detected by this study with those detected in other populations from previous studies, we have summarized both individual SNP and haplotype results in Figures 2 and 3, respectively. To our knowledge, approximately 25 SNPs in this region are reported to be associated with smoking-related and LC phenotypes. Among the 7 SNPs found to be nominally associated with SI, SQ, and SC, rs6495308 in CHRNA3 was also found associated with CPD in two European American samples [24]. All the remaining 6 SNPs (1 in CHRNA5, 1 in CHRNB4, and 4 in the intergenic regions) were unique to this Korean sample. It would be of great interest to individually genotype additional SNPs in this region in the sample, especially rs1051730 and rs16969968, which showed unequivocally significant association signals with SQ or CPD in several previous studies [15], [54], [55]. Recently, Le Marchand et al. [56] used urinary biomarkers to test whether rs1051730 and rs16969968 are associated with a higher level of nicotine and exposure to tobacco-specific carcinogenic substances per cigarette dose in a Hawaii study of 583 men and women of European, Japanese, or Native Hawaiian ancestry who were long-term smokers of more than 10 cigarettes per day. Although the T of rs1051730 and A of rs16969968 alleles were less common in Japanese Americans (3%) and Native Hawaiians (∼19%) than in European Americans (34%), those investigators found that carriers of T of rs1051730 or A of rs16969968 extract a greater amount of nicotine (P = 0.004 and 0.003, respectively) and those A carriers of the rs16969968 had a higher internal dose of total 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanol (P = 0.03) per cigarette than non-carriers. In another study with a sample size of ∼17,300 subjects from five LC studies and four upper aerodigestive tract cancer studies, Lips et al. [57] revealed no association between rs16966968 and SI or SC, age at SI, or age at SC. However, when cancer cases and controls were combined (after adjustment for case/control status), the adjusted mean difference between the two homozygote genotypes was 1.2 CPD (P<0.0001). Moreover, the rs16969968 genotype was associated with LC in both former (P<0.0001) and current (P<0.0001) smokers, and a marginally significant trend was observed in never smokers (P = 0.07).

[Figure omitted. See PDF.]

Figure 3. A summary of all reported haplotypes that have been significantly associated with different smoking behaviors in the AA, EA and Korean samples.

https://doi.org/10.1371/journal.pone.0012183.g003

At the haplotype level, in the three long-term EA smoker cohorts recruited in Utah, Wisconsin, and by the NHLBI Lung Health Study, Weiss et al. [19] found that a 5-SNP haplotype HA (CCAGA) formed by rs680244-rs569207-rs16969968-rs578776-rs1051730 (denoted EA_H1), had a risk effect on high FTND score (>6) compared with low FTND score (<4) (OR = 1.50, 95% CI: 1.21, 1.86, P = 1.3×10−4), and the 5-SNP haplotype HC (CTGAG) formed by the same SNPs as in HA had a protective effect on high ND (OR = 0.66, 95% CI: 0.52, 0.85, P = 1.07×10−3). Such findings corroborated for ND assessed by Primary Dependence Motives (PDM) score [46]. Further, our previous study in the Mid-South Tobacco Family (MSTF) cohort revealed 3 additional haplotypes in EAs [18], formed by rs6495308-rs3743075-rs8040868 (EA_H2), rs3743075-rs8040868-rs6495309 (EA_H3), and rs17408276-rs16969968-rs615470 (EA_H4), that were associated with smoking-related phenotypes. In comparison with EA_H1, although EA_H4 was physically encompassed by EA_H1, and they share 1 SNP (i.e., rs16969968), EA_H2 and EA_H3, which have 2 overlapping SNPs (i.e., rs3743075 and rs8040868), were physically separate from EA_H1. Also in the MSTF cohort, we detected two AA haplotypes formed by rs1317286-rs12914385-rs2869546 (AA_H1) and rs2869546-rs6495308-rs3743075 (AA_H2), which share 1 SNP (i.e., rs2869546) (Figure 3). Comparing the detected haplotypes in AAs and EAs, AA_H1 does not physically overlap with any of the 4 EA haplotypes, and AA_H2 shares 2 SNPs (i.e., rs6495308 and rs3743075) with EA_H2, as well as 1 SNP (i.e., rs3743075) with EA_H3. Although the overall physical region spanned by the 9 overlapping Korean haplotypes identified in this study from rs16969948 to rs6495316 does encompass all 4 EA haplotypes (i.e., EA_H1, EA_H2, EA_H3, and EA_H4) as well as both AA haplotypes (i.e., AA_H1 and AA_H2), these 9 Korean haplotypes do not share any SNPs with EA_H1 EA_H3, EA_H4 or AA_H1. Nevertheless, 3 Korean haplotypes (Korean_H3, Korean_H4, and Korean_H5) do share an SNP (i.e., rs6495308) with EA_H2 and AA_H2. Overall, the Korean haplotypes provide some support for the importance of the EA and AA haplotypes previously associated with smoking-related phenotypes However, because the SNP compositions differ across EA, AA and Korean haplotypes, the haplotype association signals are not directly comparable. To get a clearer picture of cross-population haplotypic effects, it is essential to examine the effects of haplotypes with the same SNP combinations on smoking behavior across different ethnic populations given that haplotype frequencies differ so dramatically in populations of different ancestries.

The problem of multiple hypothesis testing in biomedical studies is an important yet complex issue that needs to be considered carefully [58]. For complex traits such as smoking behavior, variants in genes of multiple biological pathways are likely to be involved in producing the phenotype through their mutual interactions and interactions with environmental factors. Therefore, genetic effects of common variants for smoking behavior are often weak secondary to genetic heterogeneity, variable expressivity, and low penetrance; and exactly how to correct for multiple testing remains a debatable and profound topic. For example, for gene-based studies such as this study focusing on the CHRNA5/A3/B4 cluster, Neale and Sham [59] have emphasized the importance of replication rather than the sole interest of detecting association signals with very low P values, and the authors have challenged the need for correcting for multiple testing in this specific scenario. Conversely, from a statistical viewpoint, when a data set is tested in multiple angles, the threshold for statistical significance should be adjusted to reduce the inflated type I error. To correct for multiple testing, as suggested by Feise [60], a prudently chosen balance needs to be reached for a study's statistical significance in consideration of multiple factors such as the magnitude of the genetic effect, biological function(s) of the marker(s) of interest, the study's quality, as well as the collective supportive evidence of the genetic locus from other independent studies [60]. Specifically, in this study, we applied the Bonferroni correction procedure for the number of SNPs analyzed within each sample for a given phenotype, but we did not correct for the number of smoking-related phenotypes (i.e., SI, SQ, and SC), the number of genetic models employed (i.e., recessive, additive, and dominant), or the number of study samples (i.e., total and male sample). The rationales underlying our decisions are that the smoking-related phenotypes are highly inter-related, the male smokers constitute the predominant fraction of all smokers in the total sample, and the three genetic models cannot be treated as totally independent. Therefore, given that the variants in CHRNA5/A3/B4 have been associated with smoking behavior in various independent samples of European origin, the correction employed was a balanced decision between the application of a relatively stringent Bonferroni correction and a precaution in guarding against an over-correction that could result in unintended losses of power and valuable information [61], [62].

In sum, this is the first genetic study aimed at investigating any association of the variants of the CHRNA5/A3/B4 cluster with SI, SQ, and SC in Korean smokers. Despite the caveats regarding correction for multiple testing, our replication of prior reported significant associations of CHRNA5/A3/B4, combined with the convergent biological data implicating the functional roles of nAChR subunit genes residing in this cluster in addiction, strengthens the notion that the multiple nominally significant association signals we detected at both the SNP and haplotype levels in the Korean sample are true positives. Moreover, our study is the first to show that the nominally significant association signals have extended beyond the 5′ end of CHRNB4 to the flanking intergenic region. Further, we found that the associations of this cluster region with SI and SQ phenotypes appeared to be stronger than those with the SC phenotype in Koreans smokers. Finally, we provided evidence for nominally significant interactions among variants studied in affecting SI in the male sample. Although these findings are novel and encouraging, they do need to be confirmed in larger, independent studies of subjects with Asian ancestry.

Materials and Methods

Ethics statement

Informed written consent was obtained in advance from all participants using a form approved by all participating Institutional Review Boards. The study likewise was approved by the Institutional Review Boards of the National Institute of Health of Korea, Seoul National University, and the University of Virginia and was in accordance with the principles of the Helsinki Declaration II.

Subjects and genotyping

The relevant information for the subjects in the current study has been reported [63]. Briefly, samples from the 10,038 participants in the Korea Association Resource (KARE) Project were obtained in two recruiting areas, Ansung and Ansan, in South Korea. Participant ages ranged from 40 to 69 years. Genomic DNA was available for 10,004 participants and was genotyped with the Affymetrix Genome-Wide Human SNP Array 5.0. Genotypes were called with Bayesian Robust Linear Modeling using the Mahalanobis Distance (BRLMM) algorithm. Those samples with low call rates (N = 401), contamination (N = 11), sex inconsistencies (N = 41), cryptic relatedness (N = 608), or serious concomitant illness (N = 101) were removed. After all the samples had been filtered with this standard quality control procedure [63], 8,842 individual samples remained for the current study (4,183 males and 4,659 females).

Although we had genotyping data available for other SNPs in these samples, in this study, we focused only on the CHRNA5/A3/B4 region of chromosome 15, as this region has received much attention for its association with smoking and LC [14], [15], [16], [17], [24]. According to the genomic locations, 36 SNPs were selected. Of them, 32 met our inclusion criteria for this study with a satisfying Hardy-Weinberg equilibrium P value >10−6, MAF >0.01, and genotype call rates >95%. The 36 SNPs are described in Table 1.

Definition of smoking-related phenotypes

On the basis of the survey questionnaire, from which information on smoking status (i.e., never smoker, former smoker, occasional or light smoker, and habitual smoker) and CPD for the habitual smokers were drawn, three smoking-related phenotypes were defined: SI, SQ, and SC. Regarding SI, the first measure (called SI-1) was defined as a binary trait comparing “never smoked” and “having regular smoking experiences,” and the second measure (called SI-2) was defined as an ordinal trait with four categories: never, former, light, and habitual. Because the association results for these two SI measures appear to be similar, in this communication, we show them together under the SI phenotype. The SQ phenotype was defined as an ordinal trait with five categories (1–5) according to CPD: non-smoking, <10 CPD, 11 to 20 CPD, 21 to 30 CPD, and >31 CPD. Such assessment of SQ not only has been commonly used in the literature but also is the most productive in terms of positive findings [43]. Finally, the SC phenotype was defined as a binary trait comparing “former smoking” and “current smoking.”

Association analysis

Statistical analysis was performed using PLINK [29] and R software. For the binary phenotypes, SI-1 and SC, the association tests were performed using logistic regression analysis, with age, sex, and geographic area of recruitment as covariates under the additive and dominant models. The cumulative logit model was fit to the ordinal phenotypes SI-2 and SQ [64]. LD analysis was performed using Haploview software [30] with the option of determining haplotype blocks according to the definitions proposed by Gabriel et al. [31]. Haplotypes were reconstructed using SNPs within the genomic region containing the CHRNA5/A3/B4 cluster with sliding window sizes of 3 to 5. Haplotype-based association analysis was carried out with the haplo.stats R statistics package [33] under additive and dominant models, adjusting for age, sex, and geographic area. As we consider this study a replication of reported significant association of variants of this gene cluster with smoking behavior in an independent sample [14], [15], [16], [17], [24], we did not correct for multiple comparisons for the two genetic models and three smoking-related phenotypes, primarily because they are related, when we interpreted our findings following the concept and approach used by other researchers [25]. However, we do recognize a potential limitation of this approach, which is addressed in the Discussion section herein.

Interaction analysis of variants in the CHRNA5/A3/B4 cluster

For the interaction analysis of variants in the CHRNA5/A3/B4 region, both logistic regression and generalized multifactor dimensionality reduction (GMDR) analyses [36] were performed for binary traits (SI-1 and SC), whereas the cumulative logit model was applied to the ordinal phenotypes (SI-2 and SQ). The logistic regression model and cumulative logit model included the main effects of SNPs, interaction effects among the SNPs, as well as adjusting covariates. The log-likelihood ratio test was performed for testing the joint significance of SNP main and interaction effects. The GMDR analysis was carried out only for the binary traits, as GMDR analysis is not applicable to ordinal traits. Age, sex, and recruiting area were considered as adjusting covariates for all interaction analyses. The best SNP combination for all two- to six-locus models was selected via the prediction rate.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Positions, nucleotide variations, and allele frequencies for SNPs on chromosome 15.

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

(0.09 MB DOC)

Acknowledgments

The authors thank Hyung-Lae Kim (Center for Genome Science, National Institute of Health, South Korea), Jong-Eun Lee (DNA Link, South Korea), Nam H Cho (Department of Preventive Medicine, Ajou University, South Korea) and Chol Shin (Department of Internal Medicine, Korea University Ansan Hospital, South Korea) for their great efforts in generating and providing this valuable data. We are grateful for the invaluable contributions of clinical information and tissue samples by the participants in this study, as well as for the dedicated work of the research staff at the different clinical sites. We are grateful to the reviewers for their constructive comments and valuable suggestions, which have helped improve the quality of the paper. We also thank Dr. David L. Bronson for his excellent English language editing.

Author Contributions

Conceived and designed the experiments: MDL DY JYL BGH TP. Performed the experiments: MDL DY JYL BGH TP. Analyzed the data: DY TN. Contributed reagents/materials/analysis tools: MDL DY JZM TP. Wrote the paper: MDL DY TN TJP JZM TP.

Citation: Li MD, Yoon D, Lee J-Y, Han B-G, Niu T, Payne TJ, et al. (2010) Associations of Variants in CHRNA5/A3/B4 Gene Cluster with Smoking Behaviors in a Korean Population. PLoS ONE5(8): e12183. https://doi.org/10.1371/journal.pone.0012183

References

1. Wipfli H, Samet JM (2009) Global economic and health benefits of tobacco control: part 1. Clin Pharmacol Ther 86: 263–271.H. WipfliJM Samet2009Global economic and health benefits of tobacco control: part 1.Clin Pharmacol Ther86263271

2. Jee SH, Samet JM, Ohrr H, Kim JH, Kim IS (2004) Smoking and cancer risk in Korean men and women. Cancer Causes Control 15: 341–348.SH JeeJM SametH. OhrrJH KimIS Kim2004Smoking and cancer risk in Korean men and women.Cancer Causes Control15341348

3. KASH (1991) Surveys on the smoking habits in Korea, 1980, 1985, and 1990. pp. 6–7.KASH1991Surveys on the smoking habits in Korea, 1980, 1985, and 199067

4. APA (1994) American Psychiatric Association. Washington DC: Diagnostic and Statistical Manual of Mental Disorders: American Psychiatric Association. APA1994American Psychiatric Association.Washington DCDiagnostic and Statistical Manual of Mental Disorders: American Psychiatric Association

5. Broms U, Silventoinen K, Madden PA, Heath AC, Kaprio J (2006) Genetic architecture of smoking behavior: a study of Finnish adult twins. Twin Res Hum Genet 9: 64–72.U. BromsK. SilventoinenPA MaddenAC HeathJ. Kaprio2006Genetic architecture of smoking behavior: a study of Finnish adult twins.Twin Res Hum Genet96472

6. Xian H, Scherrer JF, Madden PA, Lyons MJ, Tsuang M, et al. (2003) The heritability of failed smoking cessation and nicotine withdrawal in twins who smoked and attempted to quit. Nicotine Tob Res 5: 245–254.H. XianJF ScherrerPA MaddenMJ LyonsM. Tsuang2003The heritability of failed smoking cessation and nicotine withdrawal in twins who smoked and attempted to quit.Nicotine Tob Res5245254

7. Lessov CN, Martin NG, Statham DJ, Todorov AA, Slutske WS, et al. (2004) Defining nicotine dependence for genetic research: evidence from Australian twins. Psychol Med 34: 865–879.CN LessovNG MartinDJ StathamAA TodorovWS Slutske2004Defining nicotine dependence for genetic research: evidence from Australian twins.Psychol Med34865879

8. Li MD, Cheng R, Ma JZ, Swan GE (2003) A meta-analysis of estimated genetic and environmental effects on smoking behavior in male and female adult twins. Addiction 98: 23–31.MD LiR. ChengJZ MaGE Swan2003A meta-analysis of estimated genetic and environmental effects on smoking behavior in male and female adult twins.Addiction982331

9. Gaimarri A, Moretti M, Riganti L, Zanardi A, Clementi F, et al. (2007) Regulation of neuronal nicotinic receptor traffic and expression. Brain Res Rev 55: 134–143.A. GaimarriM. MorettiL. RigantiA. ZanardiF. Clementi2007Regulation of neuronal nicotinic receptor traffic and expression.Brain Res Rev55134143

10. Picciotto MR, Caldarone BJ, King SL, Zachariou V (2000) Nicotinic receptors in the brain. Links between molecular biology and behavior. Neuropsychopharmacology 22: 451–465.MR PicciottoBJ CaldaroneSL KingV. Zachariou2000Nicotinic receptors in the brain. Links between molecular biology and behavior.Neuropsychopharmacology22451465

11. Watkins SS, Koob GF, Markou A (2000) Neural mechanisms underlying nicotine addiction: acute positive reinforcement and withdrawal. Nicotine Tob Res 2: 19–37.SS WatkinsGF KoobA. Markou2000Neural mechanisms underlying nicotine addiction: acute positive reinforcement and withdrawal.Nicotine Tob Res21937

12. Li MD, Burmeister M (2009) New insights into the genetics of addiction. Nat Rev Genet 10: 225–231.MD LiM. Burmeister2009New insights into the genetics of addiction.Nat Rev Genet10225231

13. Lessov-Schlaggar CN, Pergadia ML, Khroyan TV, Swan GE (2008) Genetics of nicotine dependence and pharmacotherapy. Biochem Pharmacol 75: 178–195.CN Lessov-SchlaggarML PergadiaTV KhroyanGE Swan2008Genetics of nicotine dependence and pharmacotherapy.Biochem Pharmacol75178195

14. Saccone SF, Hinrichs AL, Saccone NL, Chase GA, Konvicka K, et al. (2007) Cholinergic nicotinic receptor genes implicated in a nicotine dependence association study targeting 348 candidate genes with 3713 SNPs. Hum Mol Genet 16: 36–49.SF SacconeAL HinrichsNL SacconeGA ChaseK. Konvicka2007Cholinergic nicotinic receptor genes implicated in a nicotine dependence association study targeting 348 candidate genes with 3713 SNPs.Hum Mol Genet163649

15. Thorgeirsson TE, Geller F, Sulem P, Rafnar T, Wiste A, et al. (2008) A variant associated with nicotine dependence, lung cancer and peripheral arterial disease. Nature 452: 638–642.TE ThorgeirssonF. GellerP. SulemT. RafnarA. Wiste2008A variant associated with nicotine dependence, lung cancer and peripheral arterial disease.Nature452638642

16. 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.RJ HungJD McKayV. GaborieauP. BoffettaM. Hashibe2008A susceptibility locus for lung cancer maps to nicotinic acetylcholine receptor subunit genes on 15q25.Nature452633637

17. 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.CI AmosX. WuP. BroderickIP GorlovJ. Gu2008Genome-wide association scan of tag SNPs identifies a susceptibility locus for lung cancer at 15q25.1.Nat Genet40616622

18. Li MD, Xu Q, Lou XY, Payne TJ, Niu T, et al. (2010) Association and interaction analysis of variants in CHRNA5/CHRNA3/CHRNB4 gene cluster with nicotine dependence in African and European Americans. Am J Med Genet B Neuropsychiatr Genet 153B: 745–756.MD LiQ. XuXY LouTJ PayneT. Niu2010Association and interaction analysis of variants in CHRNA5/CHRNA3/CHRNB4 gene cluster with nicotine dependence in African and European Americans.Am J Med Genet B Neuropsychiatr Genet153B745756

19. Weiss RB, Baker TB, Cannon DS, von Niederhausern A, Dunn DM, et al. (2008) A candidate gene approach identifies the CHRNA5-A3-B4 region as a risk factor for age-dependent nicotine addiction. PLoS Genet 4: e1000125.RB WeissTB BakerDS CannonA. von NiederhausernDM Dunn2008A candidate gene approach identifies the CHRNA5-A3-B4 region as a risk factor for age-dependent nicotine addiction.PLoS Genet4e1000125

20. Sherva R, Wilhelmsen K, Pomerleau CS, Chasse SA, Rice JP, et al. (2008) Association of a single nucleotide polymorphism in neuronal acetylcholine receptor subunit alpha 5 (CHRNA5) with smoking status and with ‘pleasurable buzz’ during early experimentation with smoking. Addiction 103: 1544–1552.R. ShervaK. WilhelmsenCS PomerleauSA ChasseJP Rice2008Association of a single nucleotide polymorphism in neuronal acetylcholine receptor subunit alpha 5 (CHRNA5) with smoking status and with ‘pleasurable buzz’ during early experimentation with smoking.Addiction10315441552

21. Chen X, Chen J, Williamson VS, An SS, Hettema JM, et al. (2009) Variants in nicotinic acetylcholine receptors alpha5 and alpha3 increase risks to nicotine dependence. Am J Med Genet B Neuropsychiatr Genet 150B: 926–933.X. ChenJ. ChenVS WilliamsonSS AnJM Hettema2009Variants in nicotinic acetylcholine receptors alpha5 and alpha3 increase risks to nicotine dependence.Am J Med Genet B Neuropsychiatr Genet150B926933

22. Greenbaum L, Kanyas K, Karni O, Merbl Y, Olender T, et al. (2006) Why do young women smoke? I. Direct and interactive effects of environment, psychological characteristics and nicotinic cholinergic receptor genes. Mol Psychiatry 11: 312–322, 223.L. GreenbaumK. KanyasO. KarniY. MerblT. Olender2006Why do young women smoke? I. Direct and interactive effects of environment, psychological characteristics and nicotinic cholinergic receptor genes.Mol Psychiatry11312322, 223

23. Schlaepfer IR, Hoft NR, Collins AC, Corley RP, Hewitt JK, et al. (2008) The CHRNA5/A3/B4 gene cluster variability as an important determinant of early alcohol and tobacco initiation in young adults. Biol Psychiatry 63: 1039–1046.IR SchlaepferNR HoftAC CollinsRP CorleyJK Hewitt2008The CHRNA5/A3/B4 gene cluster variability as an important determinant of early alcohol and tobacco initiation in young adults.Biol Psychiatry6310391046

24. Berrettini W, Yuan X, Tozzi F, Song K, Francks C, et al. (2008) Alpha-5/alpha-3 nicotinic receptor subunit alleles increase risk for heavy smoking. Mol Psychiatry 13: 368–373.W. BerrettiniX. YuanF. TozziK. SongC. Francks2008Alpha-5/alpha-3 nicotinic receptor subunit alleles increase risk for heavy smoking.Mol Psychiatry13368373

25. Bierut LJ, Stitzel JA, Wang JC, Hinrichs AL, Grucza RA, et al. (2008) Variants in Nicotinic Receptors and Risk for Nicotine Dependence. Am J Psychiatry 165: 1163–1171.LJ BierutJA StitzelJC WangAL HinrichsRA Grucza2008Variants in Nicotinic Receptors and Risk for Nicotine Dependence.Am J Psychiatry16511631171

26. Uhl GR, Liu QR, Drgon T, Johnson C, Walther D, et al. (2008) Molecular genetics of successful smoking cessation: convergent genome-wide association study results. Arch Gen Psychiatry 65: 683–693.GR UhlQR LiuT. DrgonC. JohnsonD. Walther2008Molecular genetics of successful smoking cessation: convergent genome-wide association study results.Arch Gen Psychiatry65683693

27. Vink JM, Smit AB, de Geus EJ, Sullivan P, Willemsen G, et al. (2009) Genome-wide association study of smoking initiation and current smoking. Am J Hum Genet 84: 367–379.JM VinkAB SmitEJ de GeusP. SullivanG. Willemsen2009Genome-wide association study of smoking initiation and current smoking.Am J Hum Genet84367379

28. 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.NL SacconeJC WangN. BreslauEO JohnsonD. Hatsukami2009The CHRNA5-CHRNA3-CHRNB4 nicotinic receptor subunit gene cluster affects risk for nicotine dependence in African-Americans and in European-Americans.Cancer Res6968486856

29. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.S. PurcellB. NealeK. Todd-BrownL. ThomasMA Ferreira2007PLINK: a tool set for whole-genome association and population-based linkage analyses.Am J Hum Genet81559575

30. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21: 263–265.JC BarrettB. FryJ. MallerMJ Daly2005Haploview: analysis and visualization of LD and haplotype maps.Bioinformatics21263265

31. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, et al. (2002) The structure of haplotype blocks in the human genome. Science 296: 2225–2229.SB GabrielSF SchaffnerH. NguyenJM MooreJ. Roy2002The structure of haplotype blocks in the human genome.Science29622252229

32. Lin S, Chakravarti A, Cutler DJ (2004) Exhaustive allelic transmission disequilibrium tests as a new approach to genome-wide association studies. Nat Genet 36: 1181–1188.S. LinA. ChakravartiDJ Cutler2004Exhaustive allelic transmission disequilibrium tests as a new approach to genome-wide association studies.Nat Genet3611811188

33. Schaid DJ, Rowland CM, Tines DE, Jacobson RM, Poland GA (2002) Score tests for association between traits and haplotypes when linkage phase is ambiguous. Am J Hum Genet 70: 425–434.DJ SchaidCM RowlandDE TinesRM JacobsonGA Poland2002Score tests for association between traits and haplotypes when linkage phase is ambiguous.Am J Hum Genet70425434

34. Li MD, Lou XY, Chen G, Ma JZ, Elston RC (2008) Gene-gene interactions among CHRNA4, CHRNB2, BDNF, and NTRK2 in nicotine dependence. Biol Psychiatry 64: 951–957.MD LiXY LouG. ChenJZ MaRC Elston2008Gene-gene interactions among CHRNA4, CHRNB2, BDNF, and NTRK2 in nicotine dependence.Biol Psychiatry64951957

35. Li MD, Mangold JE, Seneviratne C, Chen GB, Ma JZ, et al. (2009) Association and interaction analyses of GABBR1 and GABBR2 with nicotine dependence in European- and African-American populations. PLoS One 4: e7055.MD LiJE MangoldC. SeneviratneGB ChenJZ Ma2009Association and interaction analyses of GABBR1 and GABBR2 with nicotine dependence in European- and African-American populations.PLoS One4e7055

36. Lou XY, Chen GB, Yan L, Ma JZ, Zhu J, et al. (2007) A Generalized Combinatorial Approach for Detecting Gene-by-Gene and Gene-by-Environment Interactions with Application to Nicotine Dependence. Am J Hum Genet 80: 1125–1137.XY LouGB ChenL. YanJZ MaJ. Zhu2007A Generalized Combinatorial Approach for Detecting Gene-by-Gene and Gene-by-Environment Interactions with Application to Nicotine Dependence.Am J Hum Genet8011251137

37. Heatherton TF, Kozlowski LT, Frecker RC, Fagerstrom KO (1991) The Fagerstrom Test for Nicotine Dependence: a revision of the Fagerstrom Tolerance Questionnaire. Br J Addict 86: 1119–1127.TF HeathertonLT KozlowskiRC FreckerKO Fagerstrom1991The Fagerstrom Test for Nicotine Dependence: a revision of the Fagerstrom Tolerance Questionnaire.Br J Addict8611191127

38. Heatherton TF, Kozlowski LT, Frecker RC, Rickert W, Robinson J (1989) Measuring the heaviness of smoking: using self-reported time to the first cigarette of the day and number of cigarettes smoked per day. Br J Addict 84: 791–799.TF HeathertonLT KozlowskiRC FreckerW. RickertJ. Robinson1989Measuring the heaviness of smoking: using self-reported time to the first cigarette of the day and number of cigarettes smoked per day.Br J Addict84791799

39. Shiffman S, Waters A, Hickcox M (2004) The nicotine dependence syndrome scale: a multidimensional measure of nicotine dependence. Nicotine Tob Res 6: 327–348.S. ShiffmanA. WatersM. Hickcox2004The nicotine dependence syndrome scale: a multidimensional measure of nicotine dependence.Nicotine Tob Res6327348

40. Piper ME, Piasecki TM, Federman EB, Bolt DM, Smith SS, et al. (2004) A multiple motives approach to tobacco dependence: the Wisconsin Inventory of Smoking Dependence Motives (WISDM-68). J Consult Clin Psychol 72: 139–154.ME PiperTM PiaseckiEB FedermanDM BoltSS Smith2004A multiple motives approach to tobacco dependence: the Wisconsin Inventory of Smoking Dependence Motives (WISDM-68).J Consult Clin Psychol72139154

41. Diaz FJ, Jane M, Salto E, Pardell H, Salleras L, et al. (2005) A brief measure of high nicotine dependence for busy clinicians and large epidemiological surveys. Aust N Z J Psychiatry 39: 161–168.FJ DiazM. JaneE. SaltoH. PardellL. Salleras2005A brief measure of high nicotine dependence for busy clinicians and large epidemiological surveys.Aust N Z J Psychiatry39161168

42. Stevens VL, Bierut LJ, Talbot JT, Wang JC, Sun J, et al. (2008) Nicotinic receptor gene variants influence susceptibility to heavy smoking. Cancer Epidemiol Biomarkers Prev 17: 3517–3525.VL StevensLJ BierutJT TalbotJC WangJ. Sun2008Nicotinic receptor gene variants influence susceptibility to heavy smoking.Cancer Epidemiol Biomarkers Prev1735173525

43. Li MD (2008) Identifying susceptibility loci for nicotine dependence: 2008 update based on recent genome-wide linkage analyses. Hum Genet 123: 119–131.MD Li2008Identifying susceptibility loci for nicotine dependence: 2008 update based on recent genome-wide linkage analyses.Hum Genet123119131

44. Wu C, Hu Z, Yu D, Huang L, Jin G, et al. (2009) Genetic variants on chromosome 15q25 associated with lung cancer risk in Chinese populations. Cancer Res 69: 5065–5072.C. WuZ. HuD. YuL. HuangG. Jin2009Genetic variants on chromosome 15q25 associated with lung cancer risk in Chinese populations.Cancer Res6950655072

45. Shiraishi K, Kohno T, Kunitoh H, Watanabe S, Goto K, et al. (2009) Contribution of nicotine acetylcholine receptor polymorphisms to lung cancer risk in a smoking-independent manner in the Japanese. Carcinogenesis 30: 65–70.K. ShiraishiT. KohnoH. KunitohS. WatanabeK. Goto2009Contribution of nicotine acetylcholine receptor polymorphisms to lung cancer risk in a smoking-independent manner in the Japanese.Carcinogenesis306570

46. Baker TB, Weiss RB, Bolt D, von Niederhausern A, Fiore MC, et al. (2009) Human neuronal acetylcholine receptor A5-A3-B4 haplotypes are associated with multiple nicotine dependence phenotypes. Nicotine Tob Res 11: 785–796.TB BakerRB WeissD. BoltA. von NiederhausernMC Fiore2009Human neuronal acetylcholine receptor A5-A3-B4 haplotypes are associated with multiple nicotine dependence phenotypes.Nicotine Tob Res11785796

47. Breitling LP, Dahmen N, Mittelstrass K, Illig T, Rujescu D, et al. (2009) Smoking cessation and variations in nicotinic acetylcholine receptor subunits alpha-5, alpha-3, and beta-4 genes. Biol Psychiatry 65: 691–695.LP BreitlingN. DahmenK. MittelstrassT. IlligD. Rujescu2009Smoking cessation and variations in nicotinic acetylcholine receptor subunits alpha-5, alpha-3, and beta-4 genes.Biol Psychiatry65691695

48. Kendler KS, Neale MC, Sullivan P, Corey LA, Gardner CO, et al. (1999) A population-based twin study in women of smoking initiation and nicotine dependence. Psychol Med 29: 299–308.KS KendlerMC NealeP. SullivanLA CoreyCO Gardner1999A population-based twin study in women of smoking initiation and nicotine dependence.Psychol Med29299308

49. Maes HH, Sullivan PF, Bulik CM, Neale MC, Prescott CA, et al. (2004) A twin study of genetic and environmental influences on tobacco initiation, regular tobacco use and nicotine dependence. Psychol Med 34: 1251–1261.HH MaesPF SullivanCM BulikMC NealeCA Prescott2004A twin study of genetic and environmental influences on tobacco initiation, regular tobacco use and nicotine dependence.Psychol Med3412511261

50. Zhang H, Cai B (2003) The impact of tobacco on lung health in China. Respirology 8: 17–21.H. ZhangB. Cai2003The impact of tobacco on lung health in China.Respirology81721

51. Carter B, Long T, Cinciripini P (2004) A meta-analytic review of the CYP2A6 genotype and smoking behavior. Nicotine Tob Res 6: 221–227.B. CarterT. LongP. Cinciripini2004A meta-analytic review of the CYP2A6 genotype and smoking behavior.Nicotine Tob Res6221227

52. Lerman C, Shields PG, Audrain J, Main D, Cobb B, et al. (1998) The role of the serotonin transporter gene in cigarette smoking. Cancer Epidemiol Biomarkers Prev 7: 253–255.C. LermanPG ShieldsJ. AudrainD. MainB. Cobb1998The role of the serotonin transporter gene in cigarette smoking.Cancer Epidemiol Biomarkers Prev7253255

53. Saccone NL, Saccone SF, Goate AM, Grucza RA, Hinrichs AL, et al. (2008) In search of causal variants: refining disease association signals using cross-population contrasts. BMC Genet 9: 58.NL SacconeSF SacconeAM GoateRA GruczaAL Hinrichs2008In search of causal variants: refining disease association signals using cross-population contrasts.BMC Genet958

54. Caporaso N, Gu F, Chatterjee N, Sheng-Chih J, Yu K, et al. (2009) Genome-wide and candidate gene association study of cigarette smoking behaviors. PLoS One 4: e4653.N. CaporasoF. GuN. ChatterjeeJ. Sheng-ChihK. Yu2009Genome-wide and candidate gene association study of cigarette smoking behaviors.PLoS One4e4653

55. Keskitalo K, Broms U, Heliovaara M, Ripatti S, Surakka I, et al. (2009) Association of serum cotinine level with a cluster of three nicotinic acetylcholine receptor genes (CHRNA3/CHRNA5/CHRNB4) on chromosome 15. Hum Mol Genet 18: 4007–4012.K. KeskitaloU. BromsM. HeliovaaraS. RipattiI. Surakka2009Association of serum cotinine level with a cluster of three nicotinic acetylcholine receptor genes (CHRNA3/CHRNA5/CHRNB4) on chromosome 15.Hum Mol Genet1840074012

56. Le Marchand L, Derby KS, Murphy SE, Hecht SS, Hatsukami D, et al. (2008) Smokers with the CHRNA lung cancer-associated variants are exposed to higher levels of nicotine equivalents and a carcinogenic tobacco-specific nitrosamine. Cancer Res 68: 9137–9140.L. Le MarchandKS DerbySE MurphySS HechtD. Hatsukami2008Smokers with the CHRNA lung cancer-associated variants are exposed to higher levels of nicotine equivalents and a carcinogenic tobacco-specific nitrosamine.Cancer Res6891379140

57. 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.EH LipsV. GaborieauJD McKayA. ChabrierRJ Hung2010Association between a 15q25 gene variant, smoking quantity and tobacco-related cancers among 17 000 individuals.Int J Epidemiol39563577

58. Bender R, Lange S (2001) Adjusting for multiple testing–when and how? J Clin Epidemiol 54: 343–349.R. BenderS. Lange2001Adjusting for multiple testing–when and how?J Clin Epidemiol54343349

59. Neale BM, Sham PC (2004) The future of association studies: gene-based analysis and replication. Am J Hum Genet 75: 353–362.BM NealePC Sham2004The future of association studies: gene-based analysis and replication.Am J Hum Genet75353362

60. Ge Y, Sealfon SC, Speed TP (2008) Some Step-down Procedures Controlling the False Discovery Rate under Dependence. Stat Sin 18: 881–904.Y. GeSC SealfonTP Speed2008Some Step-down Procedures Controlling the False Discovery Rate under Dependence.Stat Sin18881904

61. Cardon LR, Bell JI (2001) Association study designs for complex diseases. Nat Rev Genet 2: 91–99.LR CardonJI Bell2001Association study designs for complex diseases.Nat Rev Genet29199

62. Zuo L, Luo X, Krystal JH, Cramer J, Charney DS, et al. (2009) The efficacies of clozapine and haloperidol in refractory schizophrenia are related to DTNBP1 variation. Pharmacogenet Genomics 19: 437–446.L. ZuoX. LuoJH KrystalJ. CramerDS Charney2009The efficacies of clozapine and haloperidol in refractory schizophrenia are related to DTNBP1 variation.Pharmacogenet Genomics19437446

63. Cho YS, Go MJ, Kim YJ, Heo JY, Oh JH, et al. (2009) A large-scale genome-wide association study of Asian populations uncovers genetic factors influencing eight quantitative traits. Nat Genet 41: 527–534.YS ChoMJ GoYJ KimJY HeoJH Oh2009A large-scale genome-wide association study of Asian populations uncovers genetic factors influencing eight quantitative traits.Nat Genet41527534

64. Bender R, Grouven U (1997) Ordinal logistic regression in medical research. J R Coll Physicians Lond 31: 546–551.R. BenderU. Grouven1997Ordinal logistic regression in medical research.J R Coll Physicians Lond31546551