Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Psychiatr Genet. Author manuscript; available in PMC 2011 October 1.
Published in final edited form as:
PMCID: PMC3107398

Genetic association of bipolar disorder with the β3 nicotinic receptor subunit gene

Sarah M. Hartz,a Peng Lin,a Howard J. Edenberg,j Xiaoling Xuei,j Nanette Rochberg,a Scott Saccone,a Wade Berrettini,g Elliot Nelson,a John Nurnberger,k Laura J. Bierut,a John P. Rice,a and NIMH Genetics Initiative Bipolar Disorder Consortiuma,b,c,d,e,f,g,h,i,j,k,l



Owing to the clinical relationship between bipolar disorder and nicotine dependence, we investigated two research questions: (i) are genetic associations with nicotine dependence different in individuals with bipolar disorder as compared with individuals without bipolar disorder, and (ii) do loci earlier associated with nicotine dependence have pleiotropic effects on these two diseases.


Our study consisted of 916 cases with bipolar disorder and 1028 controls. On the basis of known associations with nicotine dependence, we genotyped eight single-nucleotide polymorphisms (SNPs) on chromosome 8 (three bins) in the regions of CHRNB3 and CHRNA6, and six SNPs on chromosome 15 (three bins) in the regions of CHRNA5 and CHRNA3.


To determine whether the genetic associations with nicotine dependence are different in bipolar disorder than in the general population, we compared allele frequencies of candidate SNPs between individuals with nicotine dependence only and individuals with both nicotine dependence and bipolar disorder. There were no statistical differences between these frequencies, indicating that genetic association with nicotine dependence is similar in individuals with bipolar disorder as in the general population. In the investigation of pleiotropic effects of these SNPs on bipolar disorder, two highly correlated synonymous SNPs in CHRNB3, rs4952 and rs4953, were significantly associated with bipolar disorder (odds ratio 1.7, 95% confidence interval: 1.2–2.4, P = 0.001). This association remained significant both after adjusting for a smoking covariate and analyzing the association in nonsmokers only.


Our results suggest that (i) bipolar disorder does not modify the association between nicotine dependence and nicotinic receptor subunit genes, and (ii) variants in CHRNB3/CHRNA6 are independently associated with bipolar disorder. Psychiatr Genet 00:000–000.

Keywords: analyses, CHRNA3, CHRNA5, CHRNA6, genetic association, nicotine, tobacco use disorder


There is significant comorbidity between nicotine dependence and bipolar disorder. Not only are individuals with bipolar disorder prone to smoking (Gonzalez-Pinto et al., 1998), but lifetime smoking in bipolar disorder is associated with earlier age of onset of symptoms, greater severity of symptoms, decreased functioning, history of suicide attempts, and comorbid anxiety and substance use disorders (Ostacher et al., 2006). Additional studies described lower quality of life in bipolar disorder associated with nicotine dependence (Gutierrez-Rojas et al., 2008), and poor treatment outcome in bipolar mania for smokers as compared with nonsmokers (Berk et al., 2008).

Evidence for genetic factors contributing to the risk of smoking behaviors and nicotine dependence is seen in the clustering of heavy smoking and nicotine dependence in families and the similarity of smoking behaviors in identical twins (Bierut et al., 1998; Li, 2006). The most biologically compelling associations with nicotine dependence have been found in the α5 nicotinic receptor subunit gene CHRNA5 on chromosome 15 (Bierut et al., 2007; Saccone et al., 2007). This region also encompasses the nicotinic receptor subunit genes CHRNA3 and CHRNB4 in addition to the genes IREB2 and PSMA4. This association has been replicated in many other independent studies (Berrettini et al., 2008; Spitz et al., 2008; Thorgeirsson et al., 2008; Weiss et al., 2008; Saccone et al., 2009). A second region on chromosome 8 including the nicotinic subunit receptor genes CHRNB3 and CHRNA6 has also been associated with nicotine dependence (Bierut et al., 2007; Saccone et al., 2007). This association was recently confirmed in a large metaanalysis (Thorgeirsson et al., 2010).

The molecular and genetic basis for bipolar disorder is just beginning to emerge. Family and twin studies have found strong evidence of heritability (Smoller and Finn, 2003; Kieseppa et al., 2004), but despite significant effort, only two consistent signals were replicated:α-1 subunit of the L-type voltage-gated calcium channel, and Ankyrin 3 (The Welcome Trust Case Control Consortium, 2007; Baum et al., 2008; Ferreira et al., 2008; Sklar et al., 2008; Smith et al., 2009, in preparation).

In this study, we investigated the relationship between nicotine dependence and bipolar disorder, focusing on the nicotinic receptor subunit genes associated with nicotine dependence: CHRNA5/CHRNA3/CHRNB4 (chromosome 15q25), and CHRNB3/CHRNA6 (chromosome 8p11) (Bierut et al., 2007; Saccone et al., 2007, 2008). First, we evaluated whether the association between these genes and nicotine dependence is the same in bipolar disorder as it is in the general population. Then, we investigated the possibility of pleiotropy by evaluating associations between bipolar disorder and variants in CHRNA5/ CHRNA3/CHRNB4 on chromosome 15, and CHRNB3/ CHRNA6 on chromosome 8.

Materials and methods


Three different study populations were used. A sample of 916 European Americans (EA) with bipolar disorder and 1028 EA controls were recruited for the primary analysis; this is subsequently referred to as the primary data set. The remaining two samples were recruited for replication: 595 EA patients with bipolar disorder and 394 EA controls, and 389 African American (AA) participants with bipolar disorder and 671 AA controls. All bipolar patients were recruited at 11 data collection sites as the fifth wave of data collection by the Bipolar Consortium. Participants were interviewed with the Diagnostic Interview for Genetic Studies (Nurnberger et al., 1994), medical records were obtained, and final bipolar diagnosis was based on the APA Diagnostic and Statistical Manual, Fourth edition (American Psychiatric Association, American Psychiatric Association, Task Force on DSM-IV, 2000). Control participants were collected separately through a volunteer panel and web-based psychiatric interview (Sanders et al., 2008). The replication sample (both the EA sample and AA sample) was collected and genotyped as part of the Genetic Association Information Network initiative. These data and corresponding documentation are available at dbGaP ( The specific population and study characteristics are described in detail elsewhere (Dick et al., 2003; Smith et al., 2009, in preparation).

Smoking covariates were collected in all participants. Nicotine dependence was defined by a commonly used definition of Fagerstrom Test of Nicotine Dependence (FTND) score of four or more (Heatherton et al., 1991), and smoker is defined as having smoked ≥100 cigarettes lifetime. Demographics and smoking variables for the primary data set are given in Table 1.

Table 1
Demographic characteristics of the primary data set

The mean age of participants with bipolar disorder is significantly lower than controls and there are more female participants with bipolar disorder than controls. For this reason, we adjusted for age and sex in our analyses. The participants with bipolar disorder had a higher proportion of smoking (58% of participants with bipolar disorder versus 51% of participants without bipolar disorder) and nicotine dependence (69% of smokers with bipolar disorder are nicotine dependent versus 52% of smokers without bipolar disorder).

Single-nucleotide polymorphism selection

The criteria for selection of candidate genes for the primary analysis were based on earlier studies of the genetics of nicotine dependence and related traits. Eight single-nucleotide polymorphisms (SNPs) were chosen in the regions of CHRNB3 and CHRNA6 on chromosome 8 and six SNPs were chosen in the region of CHRNA5 and CHRNA3 on chromosome 15 based on the top results from the nicotine dependence genome-wide association study (GWAS) of Bierut et al. (2007) and the focused follow-up study of Saccone et al. (2007). One hundred and seventy-three SNPs were selected as stratification SNPs. A companion study was designed to evaluate the association of bipolar disorder with several linkage peaks on chromosome 6, as a confirmatory study based on promising results published by other investigators (Schulze et al., 2004). It used 1139 candidate SNPs in an 8Mb region of chromosome 6; results are described elsewhere (Knight et al., 2010). Initial genotyping was carried out by Illumina using GoldenGate Assay and BeadArray technology (

Genotyping of the replication sample was performed using the Sequenom MassArray system (Sequenom, San Diego, California, USA). The alleles of the SNPs were discriminated by mass spectrometry. Sixteen candidate SNPs were chosen based on being in perfect linkage disequilibrium (LD) (r2=1) with our original findings (rs4952, and rs4953) in the Hapmap CEU population (Frazer et al., 2007).

Quality control

Both SNPs and participants were evaluated for accuracy of data. Sixty-seven participants were removed from the primary data when genetic data did not match with the reported data: six participants were not phenotypically European, 41 had incomplete diagnoses, and 20 had phenotypic sex that did not match genotypic sex. All SNPs were required to have call rates greater than 99%, have minor allele frequencies greater than 1%, and be consistent with the Hardy–Weinberg equilibrium (P>0.001).

Population stratification analysis

To avoid the detection of artificial associations, we used EIGENSTRAT software (Price et al., 2006) in the primary data to identify subpopulations that are genetically similar. Genotype data for 173 high-performance SNPs were analyzed across all participants. The initial results yielded evidence for three genetic clusters of individuals. As one of the principal components used to construct the clusters was statistically different between participants with bipolar disorder and controls [P=0.005 from a Wilcoxon signed-rank test, computed using SAS software (SAS, 2004)], the association analyses used this principal component as a covariate. To further identify population admixture, the population was anchored with the Hapmap samples of EAs, Yorubans, and Asians. Three individuals (all participants with bipolar disorder) were identified as outliers from the group of EAs. They were therefore dropped from the analysis.

Linkage disequilibrium

We estimated D′ and r2 correlation for all pairs of the 14 SNPs on the same chromosome using an EM algorithm as implemented in the computer program Haploview (version 4.0 (Barrett et al., 2005). The D0 between any two of the six or eight SNPs on the same chromosome is greater than 0.92. The r2 values are more variable, ranging from 0.15 to 1.0.

Using the program Tagger ( (De Bakker et al., 2005), SNPs were grouped into LD bins; every bin contains at least one tag SNP where the minimum r2 between each member of the bin is 0.8. Using an r2 threshold of 0.8, six tag SNPs captured most of the variance of the 14 SNPs. The group of association signals from each bin can then be viewed as a single signal.

Correcting for multiple testing

Although we examined a total of 14 SNPs in this study, the SNPs sharing a chromosome (both the eight SNPs on chromosome 8 and the six SNPs on chromosome 15) were highly correlated. Specifically, there were six LD bins formed by the 14 SNPs. The six bins (three on chromosome 8 and three on chromosome 15) were further correlated with the bins on the same chromosome. As the positive correlation results in fewer than six independent statistical tests, we conservatively set the α cut-off at 0.01 for our initial analysis. Statistical significance for the replication sample was set at α=0.05.

Genetic association analysis

To evaluate whether genetic associations for nicotine dependence are different in bipolar disorder as compared with the general population, a x2 test was used to compare allele frequencies between participants with nicotine dependence and bipolar disorder and participants with nicotine dependence, but no bipolar disorder. To confirm genetic association between the SNPs and nicotine dependence within bipolar disorder, we evaluated the subset of the primary data with the diagnosis of bipolar disorder and adequate smoking exposure (lifetime cigarettes ≥100). Nicotine dependence was the dependent variable. A logistic regression model was fit correcting for sex, age, and a population stratification principal component. The SNPs were coded 0, 1, and 2, indicating the number of minor alleles.

To evaluate the association between SNPs and bipolar disorder in the primary data, we used a logistic regression model correcting for sex, age, and a population stratification principal component. The SNPs were coded 0, 1, and 2, indicating the number of minor alleles. As we were looking exclusively at genes that had shown association with nicotine dependence, the interaction with nicotine dependence was of interest. For this reason, we ran three separate analyses. First, the analysis was run with the baseline model (bipolar disorder association with SNP adjusted for age, sex and population stratification), then with a three-level covariate for nicotine dependence (0, never smoker; 1, smoker with FTND <4; 2, FTND ≥4), and finally the baseline model using the subsample of nonsmokers only. Logistic regression parameters [including odds ratios (OR) and confidence intervals (CI)] were estimated using PLINK (version 1.05 (Purcell et al., 2007).

Imputation and confirmatory analyses

The goal of the confirmatory analysis was to evaluate the association between bipolar disorder and rs4952/rs4953 in an independent population. The replication sample was originally genotyped using the Affymetrix 6.0 chip (Affymetrix, Santa Clara, California, USA), which does not include rs4952 or rs4953. Using the program SNAP (Johnson et al., 2008) (, we found 16 SNPs in the genome that had an r2 of 1 with rs4952 and rs4953. One of these SNPs, rs10104038, is on the Affymetrix 6.0 platform, but it did not pass quality control and hence was not present in either the EA or the AA replication data sets.

The program IMPUTE (Marchini et al., 2007) was used to impute rs4952, rs4953, and the 16 related SNPs (r2 of 1 with rs4952 and rs4953) in both the replication data sets. EAs were imputed separately from AAs. Initially, we evaluated the imputation in EA based on the original data set. Owing to low minor allele frequency (MAF) (4%), the quality of the imputation was poor as measured by the imputation quality score (Lin et al., 2008, in preparation), the computed imputation accuracy adjusted for chance. However, imputation was carried out in AAs in the SNPs with higher allele frequency. As the association was replicated in the AA population in these imputed SNPs, we genotyped the SNPs in the remainder of the sample and tested both the imputation and association. Power calculations were performed by using Quanto (, assuming a 1 : 1 case–control design and a log-additive model.


To determine whether the association with nicotine dependence was the same in bipolar disorder as in general populations, we ran two analyses. First, we statistically compared allele frequencies in participants with nicotine dependence and bipolar disorder to allele frequencies in participants with nicotine dependence, but no bipolar disorder, and there were no statistically significant differences (Table 2). Second, we evaluated the association between nicotine dependence and the individual SNPs in a subset of the n=570 participants with both bipolar disorder and adequate exposure to nicotine (smoked >100 cigarettes lifetime). Owing to the lack of power, no statistically significant associations were found. However, the association with rs16969968, a nonsynonymous SNP in CHRNA5, had an estimated OR of 1.3 (P=0.09) for the additive model and an estimated OR of 1.54 (P=0.02) for the recessive model; ORs similar to published results (Saccone et al., 2007, 2008). In addition, the estimated ORs of the remaining SNPs were in the same direction as published results (Saccone et al., 2008). As these SNPs were chosen for their association with nicotine dependence, the combination of consistent associations with nicotine dependence and lack of differences in allele frequencies for nicotine dependence participants with and without bipolar disorder indicate the risk for nicotine dependence imparted by this genetic region is the same in bipolar disorder as it is in the general population.

Table 2
Comparison of allele frequencies in participants with nicotine dependence with and without bipolar disorder

To evaluate whether the candidate SNPs are associated with bipolar disorder, three separate analyses were performed (Table 3). In the first analysis, individual SNPs were tested for association with bipolar disorder using the covariates age, sex, and a population stratification factor. An association was seen between bipolar disorder and rs4952 and rs4953. To verify that the result was not as a result of confounding by smoking status, a second analysis added a three-level smoking covariate, which retained statistical significance and did not appreciably change the OR for rs4952 and rs4953. As there may be a component of nicotine dependence or cigarette smoking not accounted for by this three-level smoking covariate, and thus driving the observed association with rs4952 and rs4953, we stratified the data by smoking status and looked only at the nonsmokers. Interestingly, the effects of rs4952 and rs4953 remained strong and statistically significant both by adding the smoking covariate and by evaluating only the nonsmokers. The MAF of both rs4952 and rs4953 in the combined sample of cases and controls is 5.5%.

Table 3
Association analyses between bipolar disorder and candidate single-nucleotide polymorphisms

To validate these results, in the two confirmatory data sets, we imputed rs4952, rs4953, and the 16 SNPs that are in full LD (r2=1) with these SNPs (based on the Hapmap CEU population). The overlapping controls were used to evaluate the imputation results in EAs. Owing to the low allele frequency in the population (3% in the n=636 controls with known genotypes), the imputation quality score (Lin et al., 2008, in preparation), a measure of imputation quality adjusted for allele frequency, was tested and indicated poor quality of imputation of that region in EAs. Owing to varying allele frequencies in AAs for these 16 SNPs, the imputation in AAs seemed to be successful for the nine SNPs with minor allele frequencies greater than 5%. Four SNPs (rs13419479, rs10104038, rs16891521, and rs16891530) forming one LD bin in AAs were significantly associated with bipolar disorder (P=0.007, OR=1.8). The MAF for these SNPs in AAs is 6% overall, 5% in controls, and 9% in cases. The remaining five SNPs, forming a second LD bin, were not significantly associated with bipolar disorder (P=0.15, OR=1.3, MAF overall 8%, MAF bipolar 9%, and MAF controls 7%).

To validate the imputation, we genotyped rs4952, rs4953 and all 16 SNPs in both replication samples (EA and AA). This group of SNPs defines one LD bin in EA, with a MAF of 4%, and three LD bins in AA, where one bin has a MAF of 8% (the bin tagged by rs4952/rs4953 has a 0% MAF in AA). The OR between bipolar disorder and the bin tagged by rs10104038 in AA was similar to our original sample at 1.3 (CI: 1.1–2.1, P=0.03). However, the OR between bipolar disorder and this bin in EA was 0.9 (CI: 0.76–2.0, P=0.46) likely because of low power from the combination of small sample size and low MAF. The sample size necessary to detect an OR of 1.3 with a population MAFof 4% is 5700 participants (equal number of cases and controls).


Despite strong evidence of the heritability of bipolar disorder (Smoller and Finn, 2003; Kieseppa et al., 2004), GWASs have yielded limited results (The Welcome Trust Case Control Consortium, 2007; Baum et al., 2008; Ferreira et al., 2008; Sklar et al., 2008). As participants with bipolar disorder are much more likely to smoke than controls (Gonzalez-Pinto et al., 1998) and robust genetic associations between nicotine dependence and nicotinic receptor subunit genes have been found, we evaluated whether the genetic mechanism for nicotine dependence is similar for participants with bipolar disorder as in the controls, and whether nicotinic receptor subunit genes are associated with bipolar disorder.

Our first aim was to investigate the genetic associations of nicotine dependence in participants with bipolar disorder as compared with the general population. Our hypothesis was that the strong genetic associations in the general population would be replicated in a sample of participants with bipolar disorder. Studies have shown a strong association of rs16969968 with nicotine dependence (Bierut et al., 2007, 2008; Saccone et al., 2007, 2008), with evidence for a recessive model (Saccone et al., 2007, 2008). In addition, a functional study showed decreased response to a nicotine agonist in the risk allele of rs16969968 (Bierut et al., 2008). Despite the small sample size, we replicated the recessive finding in rs16969968, and found consistent directions of the ORs for the other SNPs. In addition, participants with nicotine dependence had statistically similar allele frequencies both with and without bipolar disorder. This leads us to believe that the functional association between nicotine dependence and the α5, α3, α6, and β3 nicotinic receptor subunits is similar for participants with bipolar disorder as it is in the general population.

Our second aim was to evaluate the association of nicotinic receptor subunit genes with bipolar disorder. We were interested in understanding whether there are associations between bipolar disorder and the candidate SNPs and if those associations are confounded by smoking or nicotine dependence. The strongest association we found was with the bin tagged by rs4953, a synonymous SNP in the coding region of the β3 nicotinic receptor subunit gene (Fig. 1). Although these SNPs have been significantly associated with nicotine dependence in the past (Saccone et al., 2007), the association with bipolar disorder seems to be independent of nicotine dependence. The association remains statistically significant both after controlling for nicotine dependence using a smoking covariate, and by performing the association test on a subset of the data containing only nonsmokers. This association was replicated with rs10104038, a SNP in perfect LD with rs4953 in EA located approximately 20 kb upstream of CHRNB3, using a replication sample of AA participants (the second replication sample of EA participants lacked the power to confirm the findings). This indicates that there is likely to be an association of the β3 nicotinic receptor subunit gene with bipolar disorder independent of smoking or nicotine dependence, and independent of any environmental contribution of smoking to bipolar disorder. The causal variant is unclear because of the size of the LD bin in both EA and AA, but these data suggest that there is a causal variant in LD with these SNPs.

Fig. 1
Plot of r2 values with rs13439479 (y-axis), with respect to chromosome 8 position (x-axis) of 10 single-nucleotide polymorphisms (SNPs) based on (i) the HapMap CEU population (European ancestry) and (ii) the HapMap YRI population (African ancestry). The ...

As both nicotine dependence and bipolar disorder are typically diagnosed in young adulthood, and some of the participants were as young as 17 years upon data collection, the younger participants may later develop these disorders, and are thus ‘incorrectly classified’ as controls at this time. This bias, however, would reduce power by contaminating the ‘control’ population. Therefore, despite this potential bias in the data, the conclusions are likely valid.

One study has specifically evaluated the association between bipolar disorder and nicotinic receptors. Shi et al. (2007) tested the association between CHRNB3 and bipolar disorder, and did not find a positive association. A careful inspection of their study reveals that the study did not test any SNPs correlated with our findings, and, subsequently would not be expected to find these associations.

Several bipolar disorder GWAS and meta-analyses have been conducted (The Welcome Trust Case Control Consortium, 2007; Ferreira et al., 2008; Sklar et al., 2008; Scott et al., 2009; Smith et al., 2009, in preparation). One study was genotyped on the Affymetrix 6.0 chip, which includes a single SNP in LD with rs4953, but it did not pass quality control (Smith et al., 2009, in preparation). Unfortunately, the remainder of the studies were genotyped on the Affymetrix 500K chip or the Illumina HumanHap550 Beadchip, neither of which include SNPs in LD with rs4953 (The Welcome Trust Case Control Consortium, 2007; Ferreira et al., 2008; Sklar et al., 2008; Scott et al., 2009). Furthermore, we showed in this study that imputation in EAs is inadequate in this region because of low allele frequency and unavailability of correlated SNPs in the region. Therefore, the fact that this region has not been highlighted in earlier studies does not discount our present findings.

The β3 nicotinic subunit receptor has particular functional relevance from the perspective of bipolar disorder. β3 is expressed in almost all dopamine neurons and is coassembled with α6 in dopamine striatal terminals (Gotti et al., 2006). The dopaminergic neurons are often targets of treatment of bipolar mania with antipsychotics. In addition, experiments using knockout mice indicate that loss of the β3 subunit indirectly causes defects in nicotinic acetylcholine receptor assembly, degradation, and/or trafficking (Gotti et al., 2005). The biological relationship between the β3 nicotinic subunit receptor and bipolar disorder adds to the strength of these findings.

The association of a nicotinic receptor subunit gene with bipolar disorder has marked clinical implications. It suggests that there may be a joint susceptibility to both nicotine dependence and bipolar disorder, explaining the rate of concordance of the two disorders. As schizophrenia has also increased rates of nicotine dependence (Lasser et al., 2000), is highly heritable (Sullivan et al., 2003), and is more common in families with bipolar disorder (Cardno et al., 2002), nicotinic receptors may play a role in susceptibility to other mental illnesses including schizophrenia. Although we have access to the Genetic Association Information Network schizophrenia dataset (Suarez et al., 2006), we were unable to directly test the hypothesis that this variant contributes to schizophrenia because these SNPs were not genotyped in that dataset (and imputation of these SNPs has been unreliable, as showed earlier).

The cholinergic system is also a target for treatments in attention-deficit hyperactivity disorder, depression, and dementia. Therefore, the association of nicotinic receptor subunit genes with bipolar disorder may not only have implications for diagnosis and treatment of bipolar disorder, but may also impact our understanding of the interrelationships between many psychiatric disorders.


The authors thank the participants of the study, without whom this study would not be possible. Data and biomaterials were collected as part of 10 projects that participated in the National Institute of Mental Health (NIMH) Bipolar Disorder Genetics Initiative. From 1999 to 2007, the principal investigators and coinvestigators were: Indiana University, Indianapolis, IN, R01 MH59545, John Nurnberger, MD, PhD, Marvin J. Miller, MD, Elizabeth S. Bowman, MD, N. Leela Rau, MD, P. Ryan Moe, MD, Nalini Samavedy, MD, Rif El-Mallakh, MD (University of Louisville), Husseini Manji, MD (Johnson and Johnson), Debra A. Glitz, MD (Wayne State University), Eric T. Meyer, PhD, MS (Oxford University, UK) Carrie Smiley, RN, Tatiana Foroud, PhD, Leah Flury, MS, Danielle M. Dick, PhD (Virginia Commonwealth University), Howard Edenberg, PhD; Washington University, St Louis, MO, R01 MH059534, John Rice, PhD, Theodore Reich, MD, Allison Goate, PhD, Laura Bierut, MD, K02 DA21237; Johns Hopkins University, Baltimore, MD, R01 MH59533, Melvin McInnis MD, J. Raymond DePaulo, Jr, MD, Dean F. MacKinnon, MD, Francis M. Mondimore, MD, James B. Potash, MD, Peter P. Zandi, PhD, Dimitrios Avramopoulos, and Jennifer Payne; University of Pennsylvania, PA, R01 MH59553, Wade Berrettini MD, PhD; University of California at San Francisco, CA, R01 MH60068, William Byerley MD, and Sophia Vinogradov MD; University of Iowa, IA, R01 MH059548, William Coryell MD, and Raymond Crowe MD; University of Chicago, IL, R01 MH59535, Elliot Gershon, MD, Judith Badner PhD, Francis McMahon MD, Chunyu Liu PhD, Alan Sanders MD, Maria Caserta, Steven Dinwiddie MD, Tu Nguyen, Donna Harakal; University of California at San Diego, CA, R01MH59567, John Kelsoe,MD, Rebecca McKinney, BA; Rush University, IL, R01 MH059556, William Scheftner MD, Howard M. Kravitz, DO, MPH, Diana Marta, BS, Annette Vaughn-Brown, MSN, RN, and Laurie Bederow, MA; NIMH Intramural Research Program, Bethesda, MD, 1Z01MH002810-01, Francis J. McMahon, MD, Layla Kassem, PsyD, Sevilla Detera-Wadleigh, PhD, Lisa Austin, PhD, Dennis L. Murphy, MD; Howard University, William B. Lawson MD, PhD, Evarista Nwulia MD and Maria Hipolito MD.

This study was supported by the NIH grants P50CA89392 from the National Cancer Institute and 5K02DA021237 National Institute of Drug Abuse.


Financial disclosure: Drs L.J. Bierut, A.M. Goate, A.J. Hinrichs, J.P. Rice, S.F. Saccone, and J.C. Wang are listed as inventors on a patent (US 20070258898) held by Perlegen Sciences, Inc., covering the use of certain SNPs in determining the diagnosis, prognosis, and treatment of addiction. Dr Bierut has acted as a consultant for Pfizer, Inc., in 2008.


  • American Psychiatric Association, American Psychiatric Association, Task Force on DSM-IV. Diagnostic and statistical manual of mental disorders: DSM-IV-TR. 4. Washington, DC: American Psychiatric Association; 2000. p. 943.
  • Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–265. [PubMed]
  • Baum AE, Akula N, Cabanero M, Cardona I, Corona W, Klemens B, et al. A genome-wide association study implicates diacylglycerol kinase eta (DGKH) and several other genes in the etiology of bipolar disorder. Mol Psychiatry. 2008;13:197–207. [PMC free article] [PubMed]
  • Berk M, Ng F, Wang WV, Tohen M, Lubman DI, Vieta E, Dodd S. Going up in smoke: tobacco smoking is associated with worse treatment outcomes in mania. J Affect Disord. 2008;110:126–134. [PubMed]
  • Berrettini W, Yuan X, Tozzi F, Song K, Francks C, Chilcoat H, et al. Alpha- 5/alpha-3 nicotinic receptor subunit alleles increase risk for heavy smoking. Mol Psychiatry. 2008;13:368–373. [PMC free article] [PubMed]
  • Bierut LJ, Dinwiddie SH, Begleiter H, Crowe RR, Hesselbrock V, Nurnberger JI, Jr, et al. Familial transmission of substance dependence: alcohol, marijuana, cocaine, and habitual smoking: a report from the collaborative study on the Genetics of Alcoholism. Arch Gen Psychiatry. 1998;55:982–988. [PubMed]
  • Bierut LJ, Madden PA, Breslau N, Johnson EO, Hatsukami D, Pomerleau OF, et al. Novel genes identified in a high-density genome wide association study for nicotine dependence. Hum Mol Genet. 2007;16:24–35. [PMC free article] [PubMed]
  • Bierut LJ, Stitzel JA, Wang JC, Hinrichs AL, Grucza RA, Xuei X, et al. Variants in nicotinic receptors and risk for nicotine dependence. Am J Psychiatry. 2008;165:1163–1171. [PMC free article] [PubMed]
  • Cardno AG, Rijsdijk FV, Sham PC, Murray RM, McGuffin P. A twin study of genetic relationships between psychotic symptoms. Am J Psychiatry. 2002;159:539–545. [PubMed]
  • De Bakker PI, Yelensky R, Pe’er I, Gabriel SB, Daly MJ, Altshuler D. Efficiency and power in genetic association studies. Nat Genet. 2005;37:1217–1223. [PubMed]
  • Dick DM, Foroud T, Flury L, Bowman ES, Miller MJ, Rau NL, et al. Genomewide linkage analyses of bipolar disorder: a new sample of 250 pedigrees from the National Institute of Mental Health Genetics Initiative. Am J Hum Genet. 2003;73:107–114. [PubMed]
  • Ferreira MA, O’Donovan MC, Meng YA, Jones IR, Ruderfer DM, Jones L, et al. Collaborative genome-wide association analysis supports a role for ANK3 and CACNA1C in bipolar disorder. Nat Genet. 2008;40:1056–1058. [PMC free article] [PubMed]
  • Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, Gibbs RA, et al. A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007;449:851–861. [PMC free article] [PubMed]
  • Gonzalez-Pinto A, Gutierrez M, Ezcurra J, Aizpuru F, Mosquera F, Lopez P, de Leon J. Tobacco smoking and bipolar disorder. J Clin Psychiatry. 1998;59:225–228. [PubMed]
  • Gotti C, Moretti M, Clementi F, Riganti L, McIntosh JM, Collins AC, et al. Expression of nigrostriatal alpha 6-containing nicotinic acetylcholine receptors is selectively reduced, but not eliminated, by beta 3 subunit gene deletion. Mol Pharmacol. 2005;67:2007–2015. [PubMed]
  • Gotti C, Zoli M, Clementi F. Brain nicotinic acetylcholine receptors: native subtypes and their relevance. Trends Pharmacol Sci. 2006;27:482–491. [PubMed]
  • Gutierrez-Rojas L, Gurpegui M, Ayuso-Mateos JL, Gutierrez-Ariza JA, Ruiz-Veguilla M, Jurado D. Quality of life in bipolar disorder patients: a comparison with a general population sample. Bipolar Disord. 2008;10:625–634. [PubMed]
  • Heatherton TF, Kozlowski LT, Frecker RC, Fagerstrom KO. The Fagerstrom test for nicotine dependence: a revision of the Fagerstrom tolerance questionnaire. Br J Addict. 1991;86:1119–1127. [PubMed]
  • Johnson AD, Handsaker RE, Pulit SL, Nizzari MM, O’Donnell CJ, de Bakker Paul IW. SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics. 2008;24:2938–2939. [PMC free article] [PubMed]
  • Kieseppa T, Partonen T, Haukka J, Kaprio J, Lonnqvist J. High concordance of bipolar I disorder in a nationwide sample of twins. Am J Psychiatry. 2004;161:1814–1821. [PubMed]
  • Knight JO, Rochberg N, Saccone SF, Nurnberger J, Jr, Rice J. An investigation of candidate regions for association with bipolar disorder. Am J Med Gen Neuropsychiatr Gen. 2010;153B:1292–1297. [PMC free article] [PubMed]
  • Lasser K, Boyd JW, Woolhandler S, Himmelstein DU, McCormick D, Bor DH. Smoking and mental illness: a population-based prevalence study. JAMA. 2000;284:2606–2610. [PubMed]
  • Li MD. The genetics of nicotine dependence. Curr Psychiatry Rep. 2006;8:158–164. [PubMed]
  • Marchini J, Howie B, Myers S, McVean G, Donnelly P. A new multipoint method for genome-wide association studies by imputation of genotypes. Nat Genet. 2007;39:906–913. [PubMed]
  • Nurnberger JI, Jr, Blehar MC, Kaufmann CA, York-Cooler C, Simpson SG, Harkavy-Friedman J, et al. Diagnostic interview for genetic studies. Rationale, unique features, and training. NIMH Genetics Initiative. Arch Gen Psychiatry. 1994;51:849–859. [PubMed]
  • Ostacher MJ, Nierenberg AA, Perlis RH, Eidelman P, Borrelli DJ, Tran TB, et al. The relationship between smoking and suicidal behavior, comorbidity, and course of illness in bipolar disorder. J Clin Psychiatry. 2006;67:1907–1911. [PubMed]
  • 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:904–909. [PubMed]
  • Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–575. [PubMed]
  • Saccone SF, Hinrichs AL, Saccone NL, Chase GA, Konvicka K, Madden PA, 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:36–49. [PMC free article] [PubMed]
  • Saccone NL, Saccone SF, Hinrichs AL, Stitzel JA, Duan W, Pergadia M, et al. Multiple distinct risk loci for nicotine dependence identified by dense coverage of the complete family of nicotinic receptor subunit (CHRN) genes. Am J Med Gen Neuropsychiatr Gen. 2008;150B:453–466. [PMC free article] [PubMed]
  • Saccone N, Wang JC, Breslau N, Johnson EO, Hatsukami D, Saccone S, et al. The CHRNA5-CHRNA3-CHRNB4 nicotinic receptor subunit gene cluster affects risk for nicotine dependence in African-Americans and in European-Americans. Cancer Res. 2009;69:6848–6856. [PMC free article] [PubMed]
  • Sanders AR, Duan J, Levinson DF, Shi J, He D, Hou C, et al. No significant association of 14 candidate genes with schizophrenia in a large European ancestry sample: implications for psychiatric genetics. Am J Psychiatry. 2008;165:497–506. [PubMed]
  • SAS. SAS Release 9.1.3. Cary, NC: SAS Institute; 2004.
  • Schulze TG, Buervenich S, Badner JA, Steele CJ, Detera-Wadleigh SD, Dick D, et al. Loci on chromosomes 6q and 6p interact to increase susceptibility to bipolar affective disorder in the national institute of mental health genetics initiative pedigrees. Biol Psychiatry. 2004;56:18–23. [PubMed]
  • Scott LJ, Muglia P, Kong XQ, Guan W, Flickinger M, Upmanyu R, et al. Genome-wide association and meta-analysis of bipolar disorder in individuals of European ancestry. Proc Natl Acad Sci U S A. 2009;106:7501–7506. [PubMed]
  • Shi J, Hattori E, Zou H, Badner JA, Christian SL, Gershon ES, Liu C. No evidence for association between 19 cholinergic genes and bipolar disorder. Am J Med Genet B Neuropsychiatr Genet. 2007;144B:715–723. [PMC free article] [PubMed]
  • Sklar P, Smoller JW, Fan J, Ferreira MA, Perlis RH, Chambert K, et al. Whole-genome association study of bipolar disorder. Mol Psychiatry. 2008;13:558–569. [PubMed]
  • Smoller JW, Finn CT. Family, twin, and adoption studies of bipolar disorder. Am J Med Genet C Semin Med Genet. 2003;123C:48–58. [PubMed]
  • 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:1552–1556. [PMC free article] [PubMed]
  • Suarez BK, Duan J, Sanders AR, Hinrichs AL, Jin CH, Hou C, et al. Genomewide linkage scan of 409 European-ancestry and African American families with schizophrenia: suggestive evidence of linkage at 8p23.3-p21.2 and 11p13.1-q14.1 in the combined sample. Am J Hum Genet. 2006;78:315–333. [PubMed]
  • Sullivan PF, Kendler KS, Neale MC. Schizophrenia as a complex trait: evidence from a meta-analysis of twin studies. Arch Gen Psychiatry. 2003;60:1187–1192. [PubMed]
  • The Welcome Trust Case Control Consortium. Genome-wide association study of 14 000 cases of seven common diseases and 3000 shared controls. Nature. 2007;447:661–678. [PMC free article] [PubMed]
  • Thorgeirsson TE, Geller F, Sulem P, Rafnar T, Wiste A, Magnusson KP, et al. A variant associated with nicotine dependence, lung cancer and peripheral arterial disease. Nature. 2008;452:638–642. [PubMed]
  • Thorgeirsson TE, Gudbjartsson DF, Surakka I, Vink JM, Amin N, Geller F, et al. Sequence variants at CHRNB3-CHRNA6 and CYP2A6 affect smoking behavior. Nat Genet. 2010;42:448–453. [PMC free article] [PubMed]
  • Weiss RB, Baker TB, Cannon DS, von Niederhausern A, Dunn DM, Matsunami N, et al. A candidate gene approach identifies the CHRNA5-A3-B4 region as a risk factor for age-dependent nicotine addiction. PLoS Genet. 2008;4:e1000125. [PMC free article] [PubMed]