- Open Access
Associations between the orexin (hypocretin) receptor 2 gene polymorphism Val308Ile and nicotine dependence in genome-wide and subsequent association studies
Molecular Brain volume 8, Article number: 50 (2015)
Many genetic and environmental factors are involved in the etiology of nicotine dependence. Although several candidate gene variations have been reported by candidate gene studies or genome-wide association studies (GWASs) to be associated with smoking behavior and the vulnerability to nicotine dependence, such studies have been mostly conducted with subjects with European ancestry. However, genetic factors have rarely been investigated for the Japanese population as GWASs. To elucidate genetic factors involved in nicotine dependence in Japanese, the present study comprehensively explored genetic contributors to nicotine dependence by using whole-genome genotyping arrays with more than 200,000 markers in Japanese subjects.
The subjects for the GWAS and replication study were 148 and 374 patients, respectively. A two-stage GWAS was conducted using the Fagerström Test for Nicotine Dependence (FTND), Tobacco Dependence Screener (TDS), and number of cigarettes smoked per day (CPD) as indices of nicotine dependence. For the additional association analyses, patients who underwent major abdominal surgery, patients with methamphetamine dependence/psychosis, and healthy subjects with schizotypal personality trait data were recruited. Autopsy specimens with various diseases were also evaluated. After the study of associations between more than 200,000 marker single-nucleotide polymorphisms (SNPs) and the FTND, TDS, and CPD, the nonsynonymous rs2653349 SNP (located on the gene that encodes orexin [hypocretin] receptor 2) was selected as the most notable SNP associated with FTND, with a p value of 0.0005921 in the two-stage GWAS. This possible association was replicated for the remaining 374 samples. This SNP was also associated with postoperative pain, the initiation of methamphetamine use, schizotypal personality traits, and susceptibility to goiter.
Although the p value did not reach a conventional genome-wide level of significance in our two-stage GWAS, we obtained significant results in the subsequent analyses that suggest that the rs2653349 SNP (Val308Ile) could be a genetic factor that is related to nicotine dependence and possibly pain, schizotypal personality traits, and goiter in the Japanese population.
Despite the worldwide spread of cigarette smoking, many tobacco users are reportedly unaware of the harmful chemicals in tobacco products and tobacco smoke and the wide spectrum of specific illnesses that are caused by tobacco use . They also frequently do not know that smoking can cause cancer (other than lung cancer), heart disease, stroke, and many other diseases . Direct tobacco smoking is currently responsible for the death of approximately 5 million people worldwide each year, with many of these deaths occurring prematurely [3, 4]. In Japan, the proportion of deaths that are attributable to tobacco is estimated to be 17 % for ages 30 and over . Therefore, decreasing the frequency of cigarette smoking is imperative to prevent diseases caused by tobacco use. However, many smokers are not easily able to quit smoking because of the addictive property of nicotine, a major component of tobacco.
The heritability of nicotine dependence is estimated to range from 0.40 to 0.75 [5–7], and many genetic and environmental factors have been implicated in the etiology of the disorder. To date, several candidate gene variations have been reported to be associated with smoking behavior and the vulnerability to nicotine dependence. These were found mostly by investigating single-nucleotide polymorphisms (SNPs) on each gene related to specific phenotypes or genome-wide association studies (GWASs), which evaluate genetic markers in the entire genome. By taking the former candidate gene approach, several SNPs have been found to be associated with smoking behavior or nicotine dependence [8, 9]. The latter genome-wide approach has also successfully identified numerous candidate SNPs that may contribute to the etiology of nicotine dependence or smoking-related behavior [10–23]. The strongest signal of association appears to be SNPs at 15q25, represented by a synonymous rs1051730 SNP or nonsynonymous rs16969968 SNP in the region of CHRNA5–CHRNA3–CHRNB4, which encodes neuronal nicotinic acetylcholine receptor subunits [16–21]. However, these studies were conducted with subjects with mostly European or African ancestry [10–23]. Although a recent study that used Korean samples found loci that influence smoking initiation and nicotine dependence other than those identified in studies with samples of European origin , the applicability of the results to Japanese or other Asian populations is unknown.
To elucidate the genetic factors that are involved in nicotine dependence in the Japanese population, the present study comprehensively explored genetic contributors to nicotine dependence using whole-genome genotyping arrays with more than 200,000 markers in Japanese subjects.
Genome-wide association study identifies a locus associated with nicotine dependence in a Japanese population
In the first exploratory stage analysis of the two-stage GWAS, 11,968, 10,865, and 11,244 SNPs showed p values less than the threshold (p = 0.05) and were selected as candidates for the second-stage analysis for the FTND, TDS, and CPD. In the second stage analysis and LD-based SNP pruning, 273, 250, and 248 SNPs showed p values less than the threshold (p = 0.05) for the FTND, TDS, and CPD, respectively. These SNPs were considered potent candidates, and SNPs for the replication study were selected from among these SNPs. The top 50 candidate SNPs with the lowest p values for the FTND, TDS, and CPD are presented in Tables 1, 2, and 3, respectively. The top 51–100 candidate SNPs for the FTND, TDS, and CPD are presented in Additional file 1: Table S1, Additional file 2: Table S2, and Additional file 3: Table S3, respectively. Several candidate SNPs after the second-stage analysis were found to be located within or around the gene regions that were already identified as candidate loci associated with nicotine dependence or smoking behavior (Additional file 4: Table S4). Additionally, several SNPs or related genes that were selected after the second-stage analysis commonly appeared in the lists of candidates that resulted from the analysis for different phenotypic measures of nicotine dependence (Additional file 4: Table S4).
Although the p values for the top candidate SNPs were not sufficiently low when considering the conventional genome-wide level of significance, the list of the top 50 candidate SNPs (Tables 1, 2, 3) was found to include a few nonsynonymous SNPs that involve amino-acid changes and are expected to have a greater impact on the functions of gene products than other SNPs. Therefore, the best of these candidates for the FTND and CPD, rs2653349 (Rank 50) and rs726016 (Rank 1), respectively, were considered for the further confirmatory replication study (Tables 1 and 3). In the replication study, of the two SNPs, only the rs2653349 SNP showed a significant association with FTND score (χ2 = 5.094, q = 0.048; Table 4). This SNP was within the region of the HCRTR2 gene, which encodes orexin (hypocretin) receptor 2.
Imputation analysis was further conducted for SNPs around the rs2653349 SNP for the GWAS samples. As a result, several SNPs whose genotypes were imputed exhibited strong associations similarly to the rs2653349 SNP. These SNPs (rs2653342, rs2653343, rs2653344, and rs2653347) were located at a slightly upstream region of the rs2653349 SNP (Fig. 1). However, our database search with SNPinfo Web Server (see Additional file 5: Supporting Information) revealed that all of these SNPs were in absolute linkage disequilibrium with the rs2653349 SNP (D′ = r 2 = 1). Given that genotype call ratios by whole-genome genotyping or imputation for the rs2653342, rs2653343, rs2653344, rs2653347, and rs2653349 SNPs were 99.0 %, 97.7 %, 100.0 %, 98.0 %, 100.0 %, and 100.0 %, respectively, the apparently stronger associations that were observed for the rs2653342, rs2653343, and rs2653347 SNPs than for the rs2653349 SNP were merely attributable to the incomplete genotype data for these SNPs. Therefore, the contribution of the rs2653342, rs2653343, rs2653344, and rs2653347 SNPs to nicotine dependence was regarded as statistically equal to that of the rs2653349 SNP. Considering that rs2653349 is a nonsynonymous SNP and the others are located in the intronic region of the HCRTR2 gene, one could assume that this SNP is most likely the causative SNP.
Association between rs2653349 SNP and postoperative pain in patients who underwent major open abdominal surgery
To examine whether the rs2653349 SNP that was identified in our GWAS affects individual differences in the sensitivity to some drugs and pain, we investigated the association between the rs2653349 SNP and postoperative opioid requirements and pain in another cohort who underwent major open abdominal surgery (Additional file 6: Table S5). The subjects included 112 patients who underwent major open abdominal surgery under combined general and epidural anesthesia (Additional file 6: Table S5). Significant associations were not found for the frequency of analgesic administration (U = 495, p = 0.1373; Additional file 6: Table S5) or the total dose administered (U = 480, p = 0.1091; Additional file 6: Table S5). However, a significant difference in pain intensity at rest during the first 24-h postoperative period, reflected by scores on a numerical rating scale (NRS; i.e., a 5-point verbal rating scale, in which 0 = no pain, 1 = mild pain, 2 = moderate pain, 3 = severe pain, and 4 = extremely severe pain), was found between subjects with the A/G genotype and subjects with the G/G genotype in the rs2653349 SNP. Furthermore, carriers of the G/G genotype complained more about pain than carriers of the A/G genotype (U = 212, p = 0.0017; Additional file 6: Table S5), suggesting that G/G carriers are more sensitive to postoperative pain than A/G carriers.
Association between rs2653349 SNP and the initiation of drug use in patients with METH dependence
We also investigated the association between the rs2653349 SNP and two measures of the severity of METH dependence: age at first use of the drug (years) and number of drugs used, including METH. Significant associations were not found for the number of drugs used (U = 1894, p = 0.2601; Additional file 7: Table S6). However, a significant difference was found in the first use of METH between the subjects with the A/G genotype and subjects with the G/G genotype in the rs2653349 SNP. Carriers of the A/G genotype first used METH at a significantly younger age than carriers of the G/G genotype (U = 1097, p = 0.0056; Additional file 7: Table S6), suggesting that A/G carriers might harbor a high risk of initiating METH use and becoming regular users at younger ages than G/G carriers.
Association between rs2653349 SNP and schizotypal personality trait in healthy subjects
To further examine the impact of the rs2653349 SNP on individual differences in the vulnerability to develop other psychiatric diseases in healthy subjects, we examined the possible effect of the rs2653349 genotype in the HCRTR2 gene on total Schizotypal Personality Questionnaire (SPQ) scores [25–27] using a one-way analysis of covariance (ANCOVA). To control confounding factors, age, sex, and years of education were used as covariates . This analysis indicated a significant effect of genotype (F 1,306 = 5.80, p = 0.017; Additional file 8: Table S7). We then investigated the effect of genotype on the three SPQ factors: cognitive/perceptual, interpersonal, and disorganization . Significant effects of genotype were found on the cognitive/perceptual factor (F 1,306 = 4.37, p = 0.037) and disorganization factor (F 1,306 = 4.19, p = 0.042; Additional file 9: Table S8), but no effect of genotype was found on the interpersonal factor (p > 0.09). Subjects with the risk A/G genotype of the rs2653349 SNP had higher schizotypal trait scores than subjects with the G/G genotype (Additional file 9: Table S8).
Association between rs2653349 SNP and goiter in samples of autopsy specimens
We also investigated whether the rs2653349 SNP affects the susceptibility to diseases other than psychiatric diseases in another cohort by evaluating the association between this SNP and a total of 104 diseases/phenotypes using autopsy specimens (Additional file 10: Table S9). Nominally significant associations were found between goiter, aortic aneurysm, myeloma, and amyotrophic lateral sclerosis and the rs2653349 SNP. Carriers of the A/G genotype were more frequently in the disease group than in the control group (Additional file 11: Table S10). The association between this SNP and goiter was much stronger than the associations with the other three diseases (Additional file 11: Table S10). These results suggest that A/G carriers may harbor a higher risk of developing these diseases than G/G carriers.
Previous GWASs on nicotine dependence that used samples of mainly European origin identified numerous candidate SNPs [10–23]. Several of the candidate SNPs in the present study after the second-stage analysis were found to be located within or around the gene regions that have already been identified as candidate loci associated with nicotine dependence or smoking behavior (Additional file 4: Table S4). These candidate genetic factors may be common to some populations other than Japanese. Additionally, several SNPs or related genes that were selected after the second-stage analysis commonly appeared in the lists of candidates that resulted from the analysis of different phenotypic measures of nicotine dependence (Additional file 4: Table S4). These SNPs/genes might be candidates that possibly contribute to several aspects of dependence simultaneously. Among these candidate SNPs/genes were CSMD1 and IMMP2L, which were selected as candidate genes to test associations with the FTND, TDS, and CPD (Additional file 4: Table S4). CSMD1, a tumor suppressor gene , is a known candidate for nicotine dependence . The SNPs in the region of CHRNA5–CHRNA3–CHRNB4 at 15q25, which have been reported to show strong associations with CPD in studies with subjects with European ancestry [16–18], were not included in our list of the top 100 candidates (Tables 1, 2, 3, Additional file 1: Table S1, Additional file 2: Table S2, Additional file 3: Table S3). Although definitive conclusions cannot be drawn until further studies are conducted with larger sample sizes, the SNPs in this region might not greatly affect the severity of nicotine dependence in the Japanese population.
The best candidate SNP, rs2653349, is located within the HCRTR2 gene on chromosome 6 (Table 4). The HCRTR2 gene encodes orexin receptor 2. The orexin receptor has two subtypes (orexin receptor 1 [HCRTR1] and orexin receptor 2 [HCRTR2]) and two known neuropeptide ligands (orexin A, which binds both HCRTR1 and HCRTR2, and orexin B, which binds only HCRTR2). Orexin peptides, whose mRNA initially identified them as a hypothalamus-specific peptides with neuroexcitatory activity , and their receptors are distributed in various tissues, including the hypothalamus and pituitary [32–34]. They are well known to regulate feeding behavior and states of sleep/wakefulness [32, 35]. Orexin knockout mice were reported to exhibit an attenuation of morphine dependence . The hypocretin receptor 1 antagonist SB334867 and hypocretin receptor 1/2 antagonist almorexant both dose-dependently reduced nicotine self-administration , suggesting the involvement of the orexin system in the rewarding effects.
Our further study that examined whether the rs2653349 SNP affects individual differences in the susceptibility to other diseases or phenotypic traits suggested that A/G genotype carriers of this SNP are less sensitive to postoperative pain (Additional file 6: Table S5), have a high risk of initiating METH use regularly at younger ages (Additional file 7: Table S6), have higher schizotypal trait scores (Additional file 9: Table S8), and harbor a higher risk of developing goiter, aortic aneurysm, and myeloma than G/G carriers (Additional file 11: Table S10). For this SNP , several previous reports demonstrated an association with cluster headache [39–41]. Intriguingly, carriers of the G/G genotype have been reported to have a higher disease risk than the other genotypes. Given the modulation of antinociceptive effects by the orexin system in various animal models of pain , one speculation is that G/G carriers are generally more sensitive to pain because of dampened activity of the orexin system. The severer levels of nicotine dependence in A/G genotype carriers of this SNP may be attributable to increased rewarding effects related to the orexin system, possibly leading to a higher risk of initiating METH use regularly at younger ages. The higher schizotypal trait scores in A/G carriers than in G/G carriers in the present study may be explained by orexin-induced activation of the mesolimbic dopamine system, which may underlie some of its actions in schizophrenia . Previous reports indicated that nicotine dependence was more frequent in smokers with schizophrenia compared with the general population . Chronic schizophrenia patients also have lower visual analogue scale postoperative pain scores than control patients , which appears to be consistent with the trend observed in the present study. Less understood is the involvement of the orexin system in goiter, aortic aneurysm, myeloma, and amyotrophic lateral sclerosis. However, the increased frequency of cigarette smoking associated with severer levels of nicotine dependence in A/G carriers of this SNP may cause an elevation of thiocyanate synthesis , which would then enhance the risk of goiter [47–49].
Our database search estimated that the possible impact of the amino-acid substitution from valine to isoleucine by the base substitution from G to A in the rs2653349 SNP on the structure and function of human protein is predicted to be “benign.” Indeed, few studies have investigated the functional alterations caused by the rs2653349 SNP. A previous report only showed that the SNP, located in the exon 5 region, had no effect on orexin binding sites, but the corresponding residue Val308Ile might be at the dimer interface and influence the dimerization process of the receptor . Further functional studies will be needed to provide a clearer explanation for the associations observed in human studies between the SNP and phenotypic traits.
With regard to our initial GWAS that was conducted as a consecutive two-stage analysis, the lowest exploratory p value for the entire sample was > ~1 × 10−6 (Tables 1, 2, 3), which would have been deemed genome-wide nonsignificant if only a single-stage analysis was performed to calculate conventional Bonferroni- or FDR-corrected p values for the total samples to determine statistical significance. The situation would not have been improved even if only non-synonymous SNPs had been considered in the analysis as in previous studies [51, 52] because a total of 1,413 non-synonymous SNPs were included among the 225,591 SNPs that were used in our two-stage GWAS, and the conventionally corrected p value of the rs2653349 SNP would have exceeded 0.05. However, “significant” results that are obtained as conventionally corrected p values may not always represent true associations, meaning that the results may not be necessarily replicated in other studies, and vice versa. For example, data from the National Human Genome Research Institute GWAS catalog (as of January 31, 2009) showed 1,321 entries of discovered associations with p values < 10−5, but only 550 of these entries had p values < 5 × 10−8 , which is known as a conventionally corrected conservative threshold for declaring a significant association in a GWAS [54, 55]. In both cases, truly potent candidate SNPs may be included in the outcome of the studies. Furthermore, inconsistencies between genome-wide significance and truly meaningful “significance” have been exemplified by numerous studies. For example, Bierut et al.  identified SNPs in the β3 nicotinic cholinergic receptor gene as candidates that contribute to the development of nicotine dependence based on p values that were genome-wide nonsignificant. Thorgeirsson et al.  subsequently performed a genome-wide association meta-analysis using independent samples and found that SNPs in this gene region were genome-wide significantly associated with CPD (p < 5 × 10−8). Altogether, these reports suggest that conventionally corrected p values for combined samples may not be the sole criteria for identifying true associations between SNPs and specific phenotypes. Nevertheless, we did not obtain results with markedly low p values in the present GWAS. Future studies are required to draw more definitive conclusions about associations between the rs2653349 SNP and nicotine dependence and the other traits that were examined in the present study.
Although the p value did not reach a conventional genome-wide level of significance in our two-stage GWAS, we obtained results in the subsequent analyses that suggest that carriers of the A/G genotype of the nonsynonymous rs2653349 SNP are more likely to exhibit more severe levels of nicotine dependence than carriers of the G/G genotype in the Japanese population. This result suggests that this SNP may be a genetic factor that is associated with the exacerbation of nicotine dependence and may serve as a marker that predicts the severity of nicotine dependence in the Japanese population.
The participants in the present study were a total of 999 patients who were ambulatory, able to communicate orally, and 60 years of age or older. Of these 999 patients, 148 were included in the two-stage GWAS to explore the association between polymorphisms in the entire human genome and the severity of nicotine dependence. Numerous participants in this study had various smoking habits and completed a questionnaire that contained questions about lifestyle, including alcohol consumption, smoking, diet, and cancer history [56, 57]. The detailed demographic and clinical characteristics of the subjects are provided in Additional file 12: Table S11.
The subsequent study that investigated the associations between the candidate SNPs and sensitivity to pain and analgesia included 112 patient subjects who underwent major open abdominal surgery (age, 28–80 years; 60 males and 52 females), mostly gastrectomy for gastric cancer and colectomy for colorectal cancer [58, 59]. Additionally, 203 patients with methamphetamine (METH) dependence (age, 18–69 years; 165 males and 38 females) who were recruited for a collaborative study of the Japanese Genetics Initiative for Drug Abuse (JGIDA) [60, 61] were also enrolled in the present study. Healthy subjects (age, 18–66 years; 129 males and 182 females) were also recruited for the association study between candidate SNPs and schizotypal personality traits. To investigate the impact of the candidate SNPs on the susceptibility to various diseases, samples of 2,305 consecutive Japanese autopsy cases (1,266 male and 1,039 female; mean age at death, 80.6 ± 8.8 years) who were registered in the Internet database of Japanese SNPs for geriatric research  were also used. The demographic and clinical data of all of these subjects are described in the Supporting Information (Additional file 6: Table S5, Additional file 7: Table S6, Additional file 8: Table S7, Additional file 10: Table S9) and previous reports [58–65].
Smoking behavior data
The results of the questionnaire, especially questions related to smoking, were used in the analysis of the severity of nicotine dependence. Some of them included the Fagerstrӧm Test for Nicotine Dependence (FTND; a test that yields a continuous measure of nicotine dependence)  and Tobacco Dependence Screener (TDS; a screening questionnaire for tobacco/nicotine dependence according to the International Statistical Classification of Diseases and Related Health Problems, 10th revision [ICD-10], Diagnostic and Statistical Manual of Mental Disorders, 3rd edition, revised [DSM-III-R], and DSM-IV), which consists of 10 questions . The questionnaire also included questions about the number of cigarettes smoked per day (CPD) and other smoking-related traits (see Additional file 5: Supporting Information). In the present study, the FTND, TDS, and CPD were used as measures of nicotine dependence (Additional file 12: Table S11).
Whole-genome genotyping and quality control
A total of 300 DNA samples from the patients with smoking behavior data who visited Iwata City Hospital were used for genotyping. Whole-genome genotyping was performed using the Infinium assay II (Illumina, San Diego, CA, USA) with HumanCytoSNP-12 v2.0 (total markers: 301,232) BeadChips (see Additional file 5: Supporting Information).
A total of 225,602 SNP markers survived the initial filtration process as described in the Additional file 5: Supporting Information, and a total of 225,591 SNP markers were finally used for the GWAS after filtration using the Hardy-Weinberg equilibrium test (Additional file 13: Figure S1).
Genotyping for specific candidate SNPs
The TaqMan allelic discrimination assay (Life Technologies, Carlsbad, CA, USA) was adopted for a total of 112 and 203 DNA samples from patients who underwent major abdominal surgery and patients with METH dependence/psychosis, respectively. Additional file 14: Table S12 provides the genotype distribution of the candidate SNPs. The details of genomic DNA extraction and the TaqMan allelic discrimination assay are described in the Additional file 5: Supporting Information and previous reports .
Array-based genotyping procedure for data extraction of a specific SNP
To genotype the rs2653349 SNP for a total of 311 and 2,305 DNA samples from healthy subjects with schizotypal personality trait data and autopsy specimens, respectively, array-based genotyping was conducted. The genotype data for this SNP were extracted from the overall results. In the former case, genotyping was performed using Affymetrix Genome-Wide Human SNP Array 6.0 (Affymetrix, Santa Clara, CA, USA) as previously described . In the latter case, all of the samples were genotyped using Illumina Infinium HumanExome Beadchips (Illumina, CA, USA, USA) with iScan. The details are described in the Additional file 5: Supporting Information.
Genome-wide association study
A two-stage GWAS was conducted for patients with smoking behavior data to initially investigate the association between genetic variations and the severity of nicotine dependence as an exploratory study to narrow down candidate SNPs. Among the 300 subjects genotyped, approximately half lacked smoking behavior data that were related to nicotine dependence. Therefore, a total of 148 subjects were used for our two-stage GWAS (68 and 80 subjects for the first- and second-stage analyses, respectively). In our two-stage GWAS, 225,591 SNPs were selected for the analyses.
The FTND, TDS, and CPD data were used as an index of nicotine dependence for the GWAS. Prior to the analyses, the subjects were divided into two subgroups: low dependence (FTND < 4, TDS < 3, CPD ≤ 20) and high dependence (FTND ≥ 4, TDS ≥ 3, CPD > 20). The association with SNPs was explored using Cochran-Armitage trend tests. All of the statistical analyses were performed using gPLINK v. 2.050, PLINK v. 1.07 (http://pngu.mgh.harvard.edu/purcell/plink/; accessed April 19, 2014)  and Haploview v. 4.1 .
The GWAS procedure is summarized in Additional file 13: Figure S1. In the first-stage analysis with 68 samples, the SNPs that showed p < 0.05 were selected as candidate SNPs for the second-stage analysis among the 225,591 SNPs. For these SNPs, the second-stage analysis was conducted with 80 samples, and the SNPs that showed p < 0.05 for the single analysis of this stage and combined analysis of the first and second stages were considered potent candidates. The SNPs for further replication analysis were selected from these SNPs after linkage disequilibrium (LD)-based SNP pruning to remove redundant SNPs due to strong LD with each other (see Additional file 5: Supporting Information). In the confirmatory replication study to corroborate the possible association between the SNPs that were selected after the second stage and phenotypic trait, the q values of the false discovery rate were calculated to correct for multiple testing, in addition to p values based on previous reports [72, 73]. The SNPs that showed q < 0.05 in the analysis were considered statistically significant. A log quantile-quantile (QQ) p-value plot based on the results of the GWAS for the combined samples was subsequently drawn to check the pattern of the generated p-value distribution (Additional file 15: Figure S2).
Association study for candidate SNP
For all of the genotype data used in the additional analyses after the GWAS, the distributions were checked using the χ2 test, and the absence of significant deviation from the theoretical distribution expected from Hardy-Weinberg equilibrium was confirmed (Additional file 14: Table S12). To explore the association between the rs2653349 SNP and phenotypes, the Mann–Whitney U-test, the χ2 test, and Fisher’s exact test were performed. The statistical methods and software used in the analyses for the candidate SNP are provided in the Additional file 5: Supporting Information. The criterion for significance was set at p < 0.05.
Additional in silico analysis
Organization WH. WHO report on the global tobacco epidemic, 2011. 2011.
Siahpush M, McNeill A, Hammond D, Fong GT. Socioeconomic and country variations in knowledge of health risks of tobacco smoking and toxic constituents of smoke: results from the 2002 International Tobacco Control (ITC) Four Country Survey. Tob Control. 2006;15 Suppl 3:iii65–70. doi:10.1136/tc.2005.013276.
Organization WH. WHO global report: mortality attributable to tobacco. 2011.
Mathers CD, Loncar D. Projections of global mortality and burden of disease from 2002 to 2030. PLoS Med. 2006;3(11), e442. doi:10.1371/journal.pmed.0030442.
Li MD, Cheng R, Ma JZ, Swan GE. A meta-analysis of estimated genetic and environmental effects on smoking behavior in male and female adult twins. Addiction. 2003;98(1):23–31.
Vink JM, Willemsen G, Boomsma DI. Heritability of smoking initiation and nicotine dependence. Behav Genet. 2005;35(4):397–406. doi:10.1007/s10519-004-1327-8.
Lessov-Schlaggar CN, Pergadia ML, Khroyan TV, Swan GE. Genetics of nicotine dependence and pharmacotherapy. Biochem Pharmacol. 2008;75(1):178–95. doi:10.1016/j.bcp.2007.08.018.
Batra V, Patkar AA, Berrettini WH, Weinstein SP, Leone FT. The genetic determinants of smoking. Chest. 2003;123(5):1730–9.
MacLeod SL, Chowdhury P. The genetics of nicotine dependence: relationship to pancreatic cancer. World J Gastroenterol. 2006;12(46):7433–9.
Uhl GR, Liu QR, Drgon T, Johnson C, Walther D, Rose JE. Molecular genetics of nicotine dependence and abstinence: whole genome association using 520,000 SNPs. BMC Genet. 2007;8:10. doi:10.1186/1471-2156-8-10.
Uhl GR, Liu QR, Drgon T, Johnson C, Walther D, Rose JE, et al. Molecular genetics of successful smoking cessation: convergent genome-wide association study results. Arch Gen Psychiatry. 2008;65(6):683–93. doi:10.1001/archpsyc.65.6.683.
Uhl GR, Drgon T, Johnson C, Ramoni MF, Behm FM, Rose JE. Genome-wide association for smoking cessation success in a trial of precessation nicotine replacement. Mol Med. 2010;16(11–12):513–26. doi:10.2119/molmed.2010.00052.
Drgon T, Montoya I, Johnson C, Liu QR, Walther D, Hamer D, et al. Genome-wide association for nicotine dependence and smoking cessation success in NIH research volunteers. Mol Med. 2009;15(1–2):21–7. doi:10.2119/molmed.2008.00096.
Drgon T, Johnson C, Walther D, Albino AP, Rose JE, Uhl GR. Genome-wide association for smoking cessation success: participants in a trial with adjunctive denicotinized cigarettes. Mol Med. 2009;15(7–8):268–74. doi:10.2119/molmed.2009.00040.
Liu YZ, Pei YF, Guo YF, Wang L, Liu XG, Yan H, et al. Genome-wide association analyses suggested a novel mechanism for smoking behavior regulated by IL15. Mol Psychiatry. 2009;14(7):668–80. doi:10.1038/mp.2009.3.
Liu JZ, Tozzi F, Waterworth DM, Pillai SG, Muglia P, Middleton L, et al. Meta-analysis and imputation refines the association of 15q25 with smoking quantity. Nat Genet. 2010;42(5):436–40. doi:10.1038/ng.572.
Consortium TTaG. Genome-wide meta-analyses identify multiple loci associated with smoking behavior. Nat Genet. 2010;42(5):441–7. doi:10.1038/ng.571.
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(5):448–53. doi:10.1038/ng.573.
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(1):24–35. doi:10.1093/hmg/ddl441.
Saccone NL, Saccone SF, Hinrichs AL, Stitzel JA, Duan W, Pergadia ML, 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 Genet B Neuropsychiatr Genet. 2009;150B(4):453–66. doi:10.1002/ajmg.b.30828.
Saccone NL, Culverhouse RC, Schwantes-An TH, Cannon DS, Chen X, Cichon S et al. Multiple independent loci at chromosome 15q25.1 affect smoking quantity: a meta-analysis and comparison with lung cancer and COPD. PLoS Genet. 2010;6(8). doi:10.1371/journal.pgen.1001053.
Caporaso N, Gu F, Chatterjee N, Sheng-Chih J, Yu K, Yeager M, et al. Genome-wide and candidate gene association study of cigarette smoking behaviors. PLoS ONE. 2009;4(2), e4653. doi:10.1371/journal.pone.0004653.
Lind PA, Macgregor S, Vink JM, Pergadia ML, Hansell NK, de Moor MH, et al. A genomewide association study of nicotine and alcohol dependence in Australian and Dutch populations. Twin Res Hum Genet. 2010;13(1):10–29. doi:10.1375/twin.13.1.10.
Yoon D, Kim YJ, Cui WY, Van der Vaart A, Cho YS, Lee JY, et al. Large-scale genome-wide association study of Asian population reveals genetic factors in FRMD4A and other loci influencing smoking initiation and nicotine dependence. Hum Genet. 2012;131(6):1009–21. doi:10.1007/s00439-011-1102-x.
Iijima Y, Sasaki J, Bando N, Asai T, Mouri I, Tanno Y. Development of a Japanese Version of the Schizotypal Personality Questionnaire and Factor Structure of Schizotypy. Koudouryouhoukenkyu. 2010;36:29–41.
Someya T, Sasaki T, Takahashi S. Reliability and validity of schizotypal personality questionnaire. The Proceeding of the 32nd Scientific Meeting of the University Health Care in Japan. 1994:286–90.
Raine A. The SPQ: a scale for the assessment of schizotypal personality based on DSM-III-R criteria. Schizophr Bull. 1991;17(4):555–64.
Ma X, Sun J, Yao J, Wang Q, Hu X, Deng W, et al. A quantitative association study between schizotypal traits and COMT, PRODH and BDNF genes in a healthy Chinese population. Psychiatry Res. 2007;153(1):7–15. doi:10.1016/j.psychres.2007.02.003.
Raine A, Reynolds C, Lencz T, Scerbo A, Triphon N, Kim D. Cognitive-perceptual, interpersonal, and disorganized features of schizotypal personality. Schizophr Bull. 1994;20(1):191–201.
Ma C, Quesnelle KM, Sparano A, Rao S, Park MS, Cohen MA, et al. Characterization CSMD1 in a large set of primary lung, head and neck, breast and skin cancer tissues. Cancer Biol Ther. 2009;8(10):907–16.
de Lecea L, Kilduff TS, Peyron C, Gao X, Foye PE, Danielson PE, et al. The hypocretins: hypothalamus-specific peptides with neuroexcitatory activity. Proc Natl Acad Sci U S A. 1998;95(1):322–7.
Sakurai T, Amemiya A, Ishii M, Matsuzaki I, Chemelli RM, Tanaka H, et al. Orexins and orexin receptors: a family of hypothalamic neuropeptides and G protein-coupled receptors that regulate feeding behavior. Cell. 1998;92(4):573–85.
Blanco M, Lopez M, Garcia-Caballero T, Gallego R, Vazquez-Boquete A, Morel G, et al. Cellular localization of orexin receptors in human pituitary. J Clin Endocrinol Metab. 2001;86(4):1616–9.
Adeghate E. Orexins: tissue localization, functions, and its relation to insulin secretion and diabetes mellitus. Vitam Horm. 2012;89:111–33. doi:10.1016/B978-0-12-394623-2.00007-X.
Chemelli RM, Willie JT, Sinton CM, Elmquist JK, Scammell T, Lee C, et al. Narcolepsy in orexin knockout mice: molecular genetics of sleep regulation. Cell. 1999;98(4):437–51.
Georgescu D, Zachariou V, Barrot M, Mieda M, Willie JT, Eisch AJ, et al. Involvement of the lateral hypothalamic peptide orexin in morphine dependence and withdrawal. J Neurosci. 2003;23(8):3106–11.
LeSage MG, Perry JL, Kotz CM, Shelley D, Corrigall WA. Nicotine self-administration in the rat: effects of hypocretin antagonists and changes in hypocretin mRNA. Psychopharmacology (Berl). 2010;209(2):203–12. doi:10.1007/s00213-010-1792-0.
Thompson MD, Xhaard H, Sakurai T, Rainero I, Kukkonen JP. OX1 and OX2 orexin/hypocretin receptor pharmacogenetics. Front Neurosci. 2014;8:57. doi:10.3389/fnins.2014.00057.
Rainero I, Gallone S, Valfre W, Ferrero M, Angilella G, Rivoiro C, et al. A polymorphism of the hypocretin receptor 2 gene is associated with cluster headache. Neurology. 2004;63(7):1286–8.
Schurks M, Kurth T, Geissler I, Tessmann G, Diener HC, Rosskopf D. Cluster headache is associated with the G1246A polymorphism in the hypocretin receptor 2 gene. Neurology. 2006;66(12):1917–9. doi:10.1212/01.wnl.0000215852.35329.34.
Rainero I, Rubino E, Valfre W, Gallone S, De Martino P, Zampella E, et al. Association between the G1246A polymorphism of the hypocretin receptor 2 gene and cluster headache: a meta-analysis. J Headache Pain. 2007;8(3):152–6. doi:10.1007/s10194-007-0383-x.
Chiou LC, Lee HJ, Ho YC, Chen SP, Liao YY, Ma CH, et al. Orexins/hypocretins: pain regulation and cellular actions. Curr Pharm Des. 2010;16(28):3089–100.
Borgland SL, Labouebe G. Orexin/hypocretin in psychiatric disorders: present state of knowledge and future potential. Neuropsychopharmacology. 2010;35(1):353–4. doi:10.1038/npp.2009.119.
de Leon J, Diaz FJ. A meta-analysis of worldwide studies demonstrates an association between schizophrenia and tobacco smoking behaviors. Schizophr Res. 2005;76(2–3):135–57. doi:10.1016/j.schres.2005.02.010.
Kudoh A, Ishihara H, Matsuki A. Current perception thresholds and postoperative pain in schizophrenic patients. Reg Anesth Pain Med. 2000;25(5):475–9. doi:10.1053/rapm.2000.7617.
Erdogan MF. Thiocyanate overload and thyroid disease. Biofactors. 2003;19(3–4):107–11.
Brauer VF, Below H, Kramer A, Fuhrer D, Paschke R. The role of thiocyanate in the etiology of goiter in an industrial metropolitan area. Eur J Endocrinol. 2006;154(2):229–35. doi:10.1530/eje.1.02076.
Bertelsen JB, Hegedus L. Cigarette smoking and the thyroid. Thyroid. 1994;4(3):327–31.
Hegedus L, Karstrup S, Veiergang D, Jacobsen B, Skovsted L, Feldt-Rasmussen U. High frequency of goitre in cigarette smokers. Clin Endocrinol (Oxf). 1985;22(3):287–92.
Rainero I, Gallone S, Rubino E, Ponzo P, Valfre W, Binello E, et al. Haplotype analysis confirms the association between the HCRTR2 gene and cluster headache. Headache. 2008;48(7):1108–14. doi:10.1111/j.1526-4610.2008.01080.x.
Nho K, Corneveaux JJ, Kim S, Lin H, Risacher SL, Shen L, et al. Whole-exome sequencing and imaging genetics identify functional variants for rate of change in hippocampal volume in mild cognitive impairment. Mol Psychiatry. 2013;18(7):781–7. doi:10.1038/mp.2013.24.
Kobayashi D, Nishizawa D, Takasaki Y, Kasai S, Kakizawa T, Ikeda K, et al. Genome-wide association study of sensory disturbances in the inferior alveolar nerve after bilateral sagittal split ramus osteotomy. Mol Pain. 2013;9:34. doi:10.1186/1744-8069-9-34.
Ioannidis JP, Thomas G, Daly MJ. Validating, augmenting and refining genome-wide association signals. Nat Rev Genet. 2009;10(5):318–29. doi:10.1038/nrg2544.
Hirschhorn JN, Daly MJ. Genome-wide association studies for common diseases and complex traits. Nat Rev Genet. 2005;6(2):95–108.
Risch N, Merikangas K. The future of genetic studies of complex human diseases. Science. 1996;273(5281):1516–7.
Sato N, Kageyama S, Chen R, Suzuki M, Mori H, Tanioka F, et al. Association between neuropeptide Y receptor 2 polymorphism and the smoking behavior of elderly Japanese. J Hum Genet. 2010;55(11):755–60. doi:10.1038/jhg.2010.108.
Ella E, Sato N, Nishizawa D, Kageyama S, Yamada H, Kurabe N, et al. Association between dopamine beta hydroxylase rs5320 polymorphism and smoking behaviour in elderly Japanese. J Hum Genet. 2012;57(6):385–90. doi:10.1038/jhg.2012.40.
Hayashida M, Nagashima M, Satoh Y, Katoh R, Tagami M, Ide S, et al. Analgesic requirements after major abdominal surgery are associated with OPRM1 gene polymorphism genotype and haplotype. Pharmacogenomics. 2008;9(11):1605–16. doi:10.2217/146224126.96.36.1995.
Nishizawa D, Nagashima M, Katoh R, Satoh Y, Tagami M, Kasai S, et al. Association between KCNJ6 (GIRK2) gene polymorphisms and postoperative analgesic requirements after major abdominal surgery. PLoS ONE. 2009;4(9), e7060. doi:10.1371/journal.pone.0007060.
Ujike H, Harano M, Inada T, Yamada M, Komiyama T, Sekine Y, et al. Nine- or fewer repeat alleles in VNTR polymorphism of the dopamine transporter gene is a strong risk factor for prolonged methamphetamine psychosis. Pharmacogenomics J. 2003;3(4):242–7. doi:10.1038/sj.tpj.6500189.
Ujike H. Japanese Genetics Initiative for Drug Abuse (JGIDA). Nihon Shinkei Seishin Yakurigaku Zasshi. 2004;28(5):299–302.
Sawabe M, Arai T, Kasahara I, Esaki Y, Nakahara K, Hosoi T, et al. Developments of geriatric autopsy database and Internet-based database of Japanese single nucleotide polymorphisms for geriatric research (JG-SNP). Mech Ageing Dev. 2004;125(8):547–52. doi:10.1016/j.mad.2004.06.005.
Ujike H, Kishimoto M, Okahisa Y, Kodama M, Takaki M, Inada T, et al. Association Between 5HT1b Receptor Gene and Methamphetamine Dependence. Curr Neuropharmacol. 2011;9(1):163–8. doi:10.2174/157015911795017137.
Kishi T, Kitajima T, Kawashima K, Okochi T, Yamanouchi Y, Kinoshita Y, et al. Association Analysis of Nuclear Receptor Rev-erb Alpha Gene (NR1D1) and Japanese Methamphetamine Dependence. Curr Neuropharmacol. 2011;9(1):129–32. doi:10.2174/157015911795017065.
Kishi T, Kitajima T, Tsunoka T, Okumura T, Kawashima K, Okochi T, et al. Lack of association between prokineticin 2 gene and Japanese methamphetamine dependence. Curr Neuropharmacol. 2011;9(1):133–6. doi:10.2174/157015911795016994.
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(9):1119–27.
Kawakami N, Takatsuka N, Inaba S, Shimizu H. Development of a screening questionnaire for tobacco/nicotine dependence according to ICD-10, DSM-III-R, and DSM-IV. Addict Behav. 1999;24(2):155–66.
Nishizawa D, Hayashida M, Nagashima M, Koga H, Ikeda K. Genetic polymorphisms and human sensitivity to opioid analgesics. Methods Mol Biol. 2010;617:395–420. doi:10.1007/978-1-60327-323-7_29.
Hashimoto R, Ikeda M, Ohi K, Yasuda Y, Yamamori H, Fukumoto M, et al. Genome-wide association study of cognitive decline in schizophrenia. A J Psychiatry. 2013;170(6):683–4. doi:10.1176/appi.ajp.2013.12091228.
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(3):559–75. doi:S0002-9297(07)61352-4.
Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;25(2):263–5. doi:10.1093/bioinformatics/bth457.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995;57:289–300.
Storey J. The positive False Discovery Rate: a Bayesian interpretation and the q-value. Ann Statist. 2003;31(6):2013–35.
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7(4):248–9. doi:10.1038/nmeth0410-248.
Faul F, Erdfelder E, Lang AG, Buchner A. G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Methods. 2007;39(2):175–91.
Cohen J. Statistical power analysis for the behavioral sciences (revised edition). New York: Academic; 1977.
We acknowledge Mr. Michael Arends for his assistance with editing the manuscript. We are grateful to the volunteers for their participation in this study and the anesthesiologists, surgeons, and psychiatrists at related hospitals for collecting the clinical data.
This work was supported by grants from the MEXT KAKENHI (Tokyo, Japan; no. 22790518, 23390377, 24659549, 24790544, 25116532, 26293347, and 26860360), a Grant-in-Aid for Scientific Research on Innovative Areas (Comprehensive Brain Science Network) from the MEXT (to RH), the Ministry of Health, Labour and Welfare (MHLW) of Japan (Tokyo, Japan; no. H21-3jigan-ippan-011, H22-Iyaku-015, H25-Iyaku-020, H26-Kakushintekigan-ippan-060, and 14524680), Grants-in-Aid for the U.S.-Japan Cooperative Medical Science Program, National Cancer Center Research and Development Fund, Smoking Research Foundation, and Astellas Foundation for Research on Metabolic Disorders. This work was also supported by Grants-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science, and Technology of Japan and JSPS KAKENHI Grant Numbers (A-22240072, A-25242062, B-21390459, C-21590411 and C-26461480) by Grant-in-Aid for Challenging Exploratory Research (26670481) by Grants-in-Aid for Research on Intractable Diseases (Mitochondrial Disorders) from the Ministry of Health, Labour, Welfare of Japan (23–016, 23–116 and 24–005) and by grants for scientific research from the Takeda Science Foundation.
The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
DN, SK, MH, HS, and KI were responsible for the study concept and design. NaS, HY, FT, MN, RK, YS, MeT, HU, NO, TI, NI, IS, MasaoI, MitY, NK, MJW, NN, KUA, MasanI, TA, SM, MS, MNM, YY, MikY, NoS, MM, MaT, MH, HS, and KI contributed to the acquisition of clinical data and DNA samples. DN, JH, KO, RH, KT, TA, SM, MS, MNM, YY, MikY, NoS, MM, MaT, YIT, YCS, and TS were responsible for the acquisition of experimental data. DN, KO, TA, SM, MS, MNM, YY, MikY, NoS, MM, MaT, YIT, YCS, and TS performed the statistical analysis. SK assisted with data analysis and interpretation of findings. DN, KO, RH, TA, SM, MS, MNM, YY, MikY, NoS, MM, MaT, and KI drafted the manuscript. All authors critically reviewed content and approved final version for publication.
Top 51–100 candidate SNPs possibly associated with nicotine dependence (FTND score). (DOC 131 kb)
Top 51–100 candidate SNPs possibly associated with nicotine dependence (TDS score). (DOC 127 kb)
Top 51–100 candidate SNPs possibly associated with nicotine dependence (CPD). (DOC 126 kb)
SNPs or genes commonly included in the GWAS candidates for different phenotypes. (DOC 174 kb)
Supporting information (Materials and Methods). (DOC 124 kb)
Demographic and clinical data of the subjects who underwent major abdominal surgery. (DOC 74 kb)
Demographic and clinical data of patient subjects with methamphetamine dependence/psychosis. (DOC 72 kb)
Demographic variables for subjects included in the SPQ analysis (mean ± SD). (DOC 52 kb)
Impact of the HCRTR2 gene risk variant on schizotypal personality traits (mean ± SD). (DOC 54 kb)
Demographic and clinical data of patient subjects with autopsy specimens. (DOC 139 kb)
Results of analysis of associations between the rs2653349 SNP and diseases/phenotypes. (DOC 221 kb)
Demographic data of patient subjects with clinical data related to smoking behavior. (DOC 77 kb)
Schematic illustration of the two-stage GWAS. Candidate SNPs associated with nicotine dependence in Japanese were selected. (TIFF 4809 kb)
Distribution of genotypes of the rs2653349 and rs726016 SNPs and results of tests for Hardy-Weinberg equilibrium. (DOC 69 kb)
Log quantile-quantile (QQ) p-value plot as a result of the genome-wide association study for the combined samples. (a) Plot of the results from the FTND. (b) Plot of the results from the TDS. (c) Plot of the results from the CPD. (TIFF 411 kb)