Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2012; 7(9): e46557.
Published online 2012 September 27. doi:  10.1371/journal.pone.0046557
PMCID: PMC3459914

Associations between Variation in CHRNA5-CHRNA3-CHRNB4, Body Mass Index and Blood Pressure in the Northern Finland Birth Cohort 1966

Dana C. Crawford, Editor



The CHRNA5-CHRNA3-CHRNB4 gene cluster on 15q25 has consistently been associated with smoking quantity, nicotine dependence and lung cancer. Recent research also points towards its involvement in cardiovascular homeostasis, but studies in large human samples are lacking, especially on the role of the gene cluster in blood pressure regulation.

Methodology/Principal Findings

We studied the associations between 18 single nucleotide polymorphisms (SNPs) in CHRNA5-CHRNA3-CHRNB4 and systolic blood pressure (SBP), diastolic blood pressure (DBP), and body mass index (BMI) in 5402 young adults from the Northern Finland Birth Cohort 1966. We observed some evidence for associations between two SNPs and SBP and between six SNPs and BMI; the evidence for associations with DBP was weaker. The associations with the three phenotypes were driven by different loci with low linkage disequilibrium with each other. The associations appeared more pronounced in smokers, such that the smoking-increasing alleles would predict lower SBP and BMI. Each additional copy of the rs1948 G-allele and the rs950776 A-allele reduced SBP on average by −1.21 (95% CI −2.01, −0.40) mmHg in smokers. The variants associated with BMI included rs2036534, rs6495309, rs1996371, rs6495314, rs4887077 and rs11638372 and had an average effect size of −0.38 (−0.68, −0.08) kg/m2 per an additional copy of the risk allele in smokers. Formal assessments of interactions provided weaker support for these findings, especially after adjustment for multiple testing.


Variation at 15q25 appears to interact with smoking status in influencing SBP and BMI. The genetic loci associated with SBP were in low linkage disequilibrium with those associated with BMI suggesting that the gene cluster might regulate SBP through biological mechanisms that partly differ from those regulating BMI. Further studies in larger samples are needed for more precise evaluation of the possible interactions, and to understand the mechanisms behind.


The CHRNA5-CHRNA3-CHRNB4 gene cluster at 15q25 has consistently been associated with smoking quantity and nicotine dependence [1][15], as well as with lung cancer and chronic obstructive pulmonary disease [10], [13], [16][19]. The genes in the region encode the nicotinic acetylcholine receptor (nAChR) subunits α3, α5 and β4, of which α3 and β4, and possibly also α5 are expressed in autonomic ganglia [20]. Their presence in these sites suggests also a possible involvement in the regulation of autonomic function and thus, cardiovascular homeostasis.

Animal studies exploring the role of nAChRs in cardiovascular risk factors exist but they have mainly focused on nAChR subunits α4β2 and α7. These studies have demonstrated the involvement of these subunits in central cardiovascular actions, such as the regulation of blood pressure (BP) and heart rate [21][23]. In contrast, studies on the CHRNA5-CHRNA3-CHRNB4 gene cluster and cardiovascular homeostasis in humans are scarce. Rana et al. observed in a twin study of 370 individuals that polymorphisms in the CHRNA5-CHRNA3-CHRNB4 gene cluster were associated with systolic blood pressure (SBP), circulating levels of catestatin (an endogenous antagonist of nAChR) and epinephrine. These results suggest that variation in CHRNA5-CHRNA3-CHRNB4 gene cluster, already implicated in smoking behavior, is also involved in the autonomic regulation of BP. The associations between variation in CHRNA5-CHRNA3-CHRNB4 and SBP persisted after adjustment for smoking, suggesting an effect independent of smoking. [24] Another study on humans found that polymorphisms in CHRNA3 increased susceptibility to develop peripheral arterial disease among smokers [13]. In addition, a large meta-analysis of over 24000 subjects demonstrated that CHRNA3-rs1051730 interacts with smoking status in influencing body mass index (BMI). The allele that predisposes to heavier smoking was associated with a lower BMI in those who smoked but not in non-smokers. [25]

The potential involvement of these subunits in cardiovascular homeostasis is intriguing and could at least partly explain associations between smoking and BP observed in previous epidemiological studies. Smoking has been associated with BP in several studies with smokers having on average a lower BP than non-smokers [26]. This association is often explained by the appetite-suppressing effect of smoking [27] or by smoking serving as a behavioral alternative to eating resulting in decreased food intake [28]. Both mechanisms would lead to a leaner body in smokers, and thus a lower BP. However, an alternative explanation could be a shared genetic background between smoking, BMI and BP, as suggested previously [13], [24].

Replication of published results on potential pleiotropic effects of the CHRNA5-CHRNA3-CHRNB4 gene cluster on BMI and BP in large human samples is required. It is also of interest whether the effect is purely pleiotropic or if smoking acts as a mediator or a moderator of the associations between the gene cluster and BMI/BP, as suggested by some studies [13], [25]. Evidence for associations between variation in this gene cluster and BMI/BP will also be of interest from the point of view of identifying additional, novel susceptibility loci for BMI or BP, since identification of genes involved especially in the regulation of BP has proven challenging compared to other complex traits [29].

In the present study we have tested the effects of genetic variation within CHRNA5-CHRNA3-CHRNB4 on BMI and BP and whether these effects are modified by smoking status. Our study adds further evidence to the previously reported association between variation at 15q25 and BMI in smokers. Moreover, ours is the first study, conducted in a large population-based birth cohort of over 5400 individuals, to suggest an interactive effect on BP as well.


Study population

The study sample consisted of members from the Northern Finland Birth Cohort 1966 (NFBC1966). The cohort was initiated in 1965 by inviting all pregnant mothers from the two northernmost provinces of Finland with their expected date of delivery occurring in 1966 to participate. The cohort covered over 96% of all births in the target area in 1966 (N = 12055 mothers with 12058 live-born children). [30] The children born to the cohort have been followed-up since the 24th week of gestation until adulthood, with the latest data collection being conducted at the age of 31 years. At that point, questionnaires on health, lifestyle and occupation were mailed to all cohort members alive with known address (N = 11541). The questionnaire was returned by 75% (N = 8767) of the individuals. Those living in the original target area or in the capital (Helsinki) area (N = 8463) were also invited to a clinical examination, to which 71% (N = 6033) participated. [31] Blood samples were drawn, and DNA was extracted successfully for 5753 participants. The present sample includes those with genome-wide genotype data with an informed written consent to use their data (N = 5402) (Figure 1). The Ethics Committees of the University of Oulu and Northern Ostrobothnia Hospital District approved the study.

Figure 1
Flow chart of the NFBC1966 data collections.

Phenotypic and genetic data

Phenotypic data were obtained from the questionnaire and the clinical examination at 31 years. Smoking was assessed from questions measuring current smoking behavior (during the past year) and smoking quantity. Based on those we formed a three-category variable as follows: 1) non-smokers (including never and former smokers), 2) light smokers (1–10 cigarettes/day), and 3) heavy smokers (>10 cigarettes/day). BMI (kg/m2) was calculated from measured height and weight. SBP and DBP were measured twice with a mercury sphygmomanometer with a cuff size of 14 cm×40 cm in a sitting position from the right arm after 15 minutes of rest by trained nurses using a standardized procedure and ongoing quality control [32]. We used the average of the two measurements. For those who on the same 31-year questionnaire reported as being on medication for hypertension (N = 95), we added 15 mmHg to the recorded SBP values and 10 mmHg to the recorded DBP values in accordance with earlier genome-wide association studies on BP [33].

Genome-wide genotyping was performed with Illumina HumanCNV370DUO Analysis Beadchip platform at the Broad Institute, USA. Detailed genotyping and sample quality control of the first set of data (N = 4936 individuals) have been reported previously. [34] Afterwards more samples were genotyped with the same method as before (new data including also the previously genotyped set of individuals, N = 5550). In the quality control phase, 148 samples were excluded because of contamination, gender mismatch between genotypic and phenotypic gender, duplicate genotyping, genotyping call rate <95%, relatedness (IBS pairwise sharing >0.20), consent withdrawal and low heterozygozity. After quality control 5402 individuals remained available for the analysis.

Principal components (PC) to control for potential population stratification in the subsequent association analyses were calculated from the 22 directly genotyped autosomes such that all SNPs had MAF >1%, P-value for Hardy-Weinberg Equilibrium >0.005, call rate >99.5% and no two SNPs were in LD of r2>0.2 with any other selected SNP. The data were further thinned such that 1 in 15 SNPs was selected for principal component analysis. Eigenvalue decomposition was performed using the eigen function in the statistical software R [35]. It has previously been shown that the principal components calculated from these data correspond well to the geographical background of the cohort members [34].

Statistical Analysis

For genetic data, we selected all available directly genotyped SNPs from the genome-wide data to cover the CHRNA5-CHRNA3-CHRNB4 gene cluster and its surrounding area as in our previous publication [15]. We included only common SNPs (MAF>5%) according to the established practice in many genome-wide association studies based on concerns about low statistical power to detect associations of less common variants. These selections resulted in 18 eligible SNPs for the analyses (Table 1). The association analyses between the 18 SNPs and BMI/BP were run in Plink v1.06 [36], [37] assuming an additive effect for the genotype and adjusting the analyses for gender, BMI (for the analysis with SBP/DBP as the outcome) and three first PCs to control for potential population stratification. The distributions of the residuals of BMI, SBP and DBP did not seriously deviate from normality (data not shown); thus, the outcome variables were analyzed on their original scales assuming they follow Gaussian distributions.

Table 1
Selected SNPs and their characteristics.

We performed the analyses separately in non-smokers and smokers, and in all subjects together, including a term for gene-environment (SNP×smoking) interaction, too. To account for multiple testing, we performed parametric bootstrap for gene-environment interactions [38] and applied the maxT multiple testing procedure, in which we compared each observed test statistic against the maximum of all bootstrapped test statistics within each simulation. This controls for the family-wise error rate, but also takes adequately into account the dependence structure between test statistics, i.e. the correlation structure between SNPs is preserved. Thus, a less stringent adjustment for multiple testing is obtained in comparison to the conventionally applied Bonferroni correction, which assumes all tests are mutually independent. The number of simulations was 10000, and within each simulation we took N = 5402 resamples, i.e. the original sample size. Simulations were performed in R [35]. The results are presented as estimated regression coefficients with their 95% confidence intervals (CI), and we report both the unadjusted and adjusted P-values for the interaction terms.

Additionally, we performed a series of supplemental analyses. We examined dose-response effects by separating smokers into light and heavy smokers, and differences between never, former and current smokers were tested by separating the non-smoker group into never and former smokers. Also, we tested for possible gender×SNP interactions. The analyses were performed in Plink v1.06 [36], [37] and were adjusted for the same covariates as before. Parametric bootstrap was performed as previously, and both unadjusted and adjusted P-values are reported.

Finally, we performed haplotype analyses to increase power and to refine the location of the susceptibility locus, because multiple SNPs may “tag” an untyped variant more effectively than a single typed variant [39]. We visualized the haplotype structure of the region with Haploview v3 [40] and identified two haplotype blocks by visual inspection. Possible haplotypes in the two blocks for each individual were phased with Plink v1.06 [36], [37] (using the Expectation-Maximization algorithm-based phasing), and the haplotype frequencies and consequent association analyses were estimated using the same software. We performed sliding window analyses [41] within the blocks, i.e. tested for associations between haplotypes of differing lengths and the outcome of interest, and we used only haplotypes with estimated frequency >1% in the sample. These analyses were run assuming an additive effect for the haplotype, i.e. 0, 1 or 2 copies of the haplotype, and adjusting for the same covariates as before. The analyses were run separately in non-smokers and smokers, and no interaction analyses were performed. To account for multiple testing, we performed maxT permutation of residuals [42], [43] in Plink with 10000 permutations. We report both the unadjusted and adjusted P-values.


Distributions of the study variables

Characteristics of the study sample are presented in Table 2. Among participants with both genotypic and phenotypic data available, 58% were non-smokers and 42% smokers, of whom about half were light and half heavy smokers. Men were more likely to be heavy smokers than women (31% vs. 12%), and they had on average higher SBP and DBP compared to women. There were no obvious differences in the distributions of BMI in different smoking groups in either gender. We also checked the distributions of BMI in prolonged smokers [15] by combining data on smokers from both the 14-year and the 31-year postal questionnaires, but the average BMI of persistent smokers did not differ from that of other smokers (data not shown). SBP and DBP were slightly lower in smokers compared to non-smokers in both men and women.

Table 2
Sample characteristics of the NFBC1966 study sample.

Analyses on non-smokers and smokers

Table 3 shows the estimated associations between the 18 variants in the gene cluster and SBP in non-smokers and smokers. Similar results for DBP and BMI are presented in Tables 4 and and5.5. The results are further illustrated in Figure 2, from which both the correlation structure between the SNPs and the estimated recombination rates can also be seen. There was moderate evidence for associations between two SNPs (rs1948 and rs950776) and SBP, and six SNPs (rs2036534, rs6495309, rs1996371, rs6495314, rs4887077 and rs11638372) and BMI; evidence for associations with DBP was weaker. All the observed associations were in such a direction that the smoking-increasing allele (as reported in [15]) predicted lowered SBP, DBP or BMI, and the observed associations were generally stronger among smokers. Each additional copy of the rs1948 G-allele and rs950776 A-allele were associated with a reduction of on average −1.21 (−2.01, −0.40) mmHg in SBP in smokers. The variants associated with BMI had an average effect size of −0.38 (−0.68, −0.08) kg/m2 per an additional copy of the risk allele in smokers. The results from interaction models appeared to provide some support of the interactive effects observed in the stratified analyses; however, after adjustment for multiple testing the evidence for interactions was quite weak (see unadjusted and adjusted P-values for interaction in Tables 3, ,44 and and55).

Figure 2
Associations between the CHRNA5-CHRNA3-CHRNB4 gene cluster on the three outcomes of interest.
Table 3
Estimated associations between variants in the 15q25 region and SBP by smoking status (non-smokers, smokers) in the NFBC1966.
Table 4
Estimated associations between variants in the 15q25 region and DBP by smoking status (non-smokers, smokers) in the NFBC1966.
Table 5
Estimated associations between variants in the 15q25 region and BMI by smoking status (non-smokers, smokers) in the NFBC1966.

Because the gene cluster has previously been associated with smoking quantity, we also considered the possibility that the associations found could be due to residual confounding by smoking quantity. Therefore we reran the analyses on smokers adjusted by a continuous measure of smoking quantity (cigarettes/day), but the results remained practically unchanged (data not shown).

Analyses on light and heavy smokers

For SBP, the observed effects were slightly stronger in heavy smokers than in light smokers; e.g. for rs1948 the G-allele was associated with 0.97 mmHg lower SBP (95% CI −2.20, 0.13) in light smokers whereas in heavy smokers the corresponding figures were −1.50 (−2.64, −0.37) (Table S1). For DBP, some suggestive associations that were previously seen in smokers became more evident in heavy smokers; e.g. for rs8038920 each additional copy of the G-allele predicted 0.73 mmHg lower DBP (−1.47, 0.01) in smokers while in heavy smokers the predicted reduction was −1.40 (−2.44, −0.35) mmHg/one copy of the G-allele (Table S2). In the analysis of BMI, for some of the SNPs (rs8034191, rs1051730 and rs6495309) the effects seemed to be strongest in light smokers (Table S3). The formal assessments of interaction did not give adequate support to these findings, though, especially after adjustment for multiple testing.

Analyses on never, former and current smokers

For the SNPs on which we had observed evidence for associations in current smokers (rs1948 and rs950776 with SBP; rs2036534, rs6495309, rs1996371, rs6495314, rs4887077 and rs11638372 with BMI), the associations were attenuated in former smokers (Tables S4,S5,S6). For one of the SNPs (rs12900519), however, the effect on SBP seemed to be stronger in former (−1.58 (−3.04, −0.13) mmHg/each additional copy of A-allele) than in current smokers (−0.49 (−1.57, 0.59) mmHg/each additional copy of the A-allele) (Table S4).

Gender interactions

Since the distributions of SBP, DBP and BMI differed by gender (Table 2), we studied if there is heterogeneity by gender also on the associations between the 18 SNPs and the three outcomes of interest. For the two SNPs (rs1948 and rs950776) on which we found evidence for association with SBP when the genders were analyzed together (Table 3), the effect seemed superficially to be stronger among females (Table S7). Also for an additional SNP (rs12900519) we observed a stronger effect in females. For DBP, six SNPs predicted lower DBP in females but not in males, and one SNP (rs6495309) was associated with opposing effects among males and females (Table S8). Instead for BMI, many of the associations observed when the genders were analyzed together (Table 5) seemed to be driven by stronger associations in males (Table S9). However, for all of these findings by gender, the support from interaction analyses was weak.

Multimarker analyses

The linkage disequilibrium (LD) plot of the region is displayed in Figure 3. The estimated haplotype frequencies and regression coefficients with their 95% CIs are available in Tables S10 (results for haplotypes in block1) and S11 (results for haplotypes in block2). The results are presented only for those haplotypes with an estimated frequency >1% and with an unadjusted P-value<0.05. Such results were observed only in smokers. In the single marker analysis we observed evidence for associations for SNPs from the first block and SBP, SNPs from the second block and DBP and SNPs from both blocks and BMI (Tables 35).5). The haplotype analyses were in line with this, since there was no strong evidence for haplotypes from block 2 to be associated with SBP, nor for haplotypes from block 1 with DBP, whereas there was evidence for associations between haplotypes from both blocks and BMI (Tables S10, S11).

Figure 3
The linkage disequilibrium (LD) structure of the chromosome 15q25 region.

In more detail, for SBP, and after adjustment for multiple testing, there was moderate evidence for a SBP-lowering haplotype consisting of the G-allele of rs1948 and A-allele of rs950776; the same alleles of the two SNPs that were associated with lower SBP in the single marker analyses. Instead, SBP-increasing haplotypes included the A and G alleles of the above-mentioned SNPs, respectively. The pattern for the alleles of the other SNPs included in SBP-associated haplotypes was not as clear-cut as it was for rs1948 and rs950776; the same alleles in different haplotypes could be associated either with increased or decreased SBP (Table S10).

The estimated effect sizes on SBP were bigger for many SBP-increasing haplotypes than for any of the SBP-lowering haplotypes or single SNPs in the single marker analysis (Table S10). An increase of about 6 mmHg in SBP per each additional copy of the haplotype was observed for six haplotypes of varying lengths. All of these haplotypes were estimated to be rare in the population with an allele frequency of about 1%.

Regarding BMI and the haplotypes from the first block, there was no strong evidence for associations after adjustment for multiple testing (Table S10). In the second block (Table S11), we observed some evidence for associations of the whole block-spanning haplotype GAGCGGAA or haplotypes of differing lengths but with the same alleles as in the long haplotype and lower BMI, except that in some cases the two last alleles were AA instead of GG. Again, bigger effect sizes were observed in multimarker analyses than in single marker analyses, but these effects were arising from haplotypes with low estimated population frequencies (Table S11).


Main findings

The present study suggests that in addition to the previously presented hypotheses for the association between increased smoking and lower BP, there might be shared a genetic architecture behind as well. In the single marker analyses, weak evidence was found for associations between two polymorphisms from the smoking-associated CHRNA5-CHRNA3-CHRNB4 gene cluster and lower BP. We found weak evidence for associations between smoking-increasing alleles and lower BMI as well, although the distributions of BMI did not vary across different smoking strata. If same alleles predisposed to smoking and lower BMI, we would have expected to observe on average lower BMI in smokers, which was not the case in our data. However, we believe that the association between smoking and BMI is highly confounded by socio-behavioral factors, for which we did not control when comparing the distributions of BMI across different smoking strata, because addressing the associations between smoking and BMI/BP in detail was beyond the scope of the present study.

Different SNPs were associated with BP and BMI, and even with SBP and DBP. Mainly markers from CHRNB4 were associated with SBP, pointing towards β4 subunits being mostly responsible for these effects, whereas markers from CHRNA5 and CHRNB4 were associated with BMI, emphasizing the role of α5 and β4 subunits in affecting BMI. However, further studies are needed to locate the real causal variants and the receptor subunits responsible for the observed associations. For all of the six SNPs on which we found evidence for any association with BMI, the r2 with rs1948, the most strongly associated SNP with SBP, was <0.20 indicating that the SNPs associated with BMI and SBP seem to be nearly independent of each other. The haplotype association analyses also supported this finding since mainly haplotypes from block 1 were associated with SBP, haplotypes from block 2 with DBP, and haplotypes from both blocks with BMI. The haplotype analyses also provided further support for the two SNPs found to be associated with SBP in the single marker analyses; the haplotype GA consisting of rs1948 and rs950776 predicted lower SBP even after correction for multiple testing. For BMI, the findings from the haplotype analyses were less consistent with the findings from the single marker analyses.

The results from the additional, rather explorative analyses suggested that the effects of the SNPs on SBP would be stronger in heavy smokers than in light smokers. The associations between the SNPs and SBP/BMI were attenuated in former smokers, suggesting that the associations would weaken over time after cessation. In addition, there might be heterogeneity by gender in the associations. However, additional stratifications reduced the sample sizes resulting in very wide confidence intervals, and moreover, the interaction analyses provided no adequate support for these findings. Studies with larger sample sizes are needed to find further support for these suggestive results.

We also tried to address the fact the possibility that the associations seen are due to residual confounding by smoking quantity. By adjusting the analyses in smokers by a continuous measure of cigarettes/day did not change the results obtained previously, suggesting that the associations seen are not due to residual confounding. However, in all of the analyses the measure of smoking was self-reported, possibly leading to underestimation of the true effect. For example, the availability and the use of cotinine levels might control better for the actual quantity of smoked cigarettes.

Comparison with previous studies on the gene cluster

Of the SNPs that we found evidence for an association with BP, rs1948 in CHRNB4 has previously been associated with heavy smoking in adulthood [15], nicotine dependence [44] and with early age of initiation of tobacco and alcohol use [45]. Another SNP for which we observed evidence for an association with lower BMI, rs1996371 and especially its allele G, has previously been associated with heavy smoking and nicotine dependence [7], [12], [15]. However, none of the SNPs that appeared to be associated with SBP/BMI in the present study have previously reached genome-wide significance and thus been identified as potential susceptibility loci for BMI or BP according to the NHGRI catalog of published genome-wide association studies [46] or based on our literature search.

Only two SNPs from the CHRNA5-CHRNA3-CHRNB4 gene cluster have previously been reported to be associated with BMI or SBP. SNP rs1051730 in CHRNA3 reported to be associated with BMI [25] and which is also the strongest signal previously associated with smoking quantity, was not very strongly associated either with BP or BMI in the present study. However, the effects were in the direction that the smoking-associated A-allele was predictive of decreased SBP and BMI in smokers, which is in line with the previous publication reporting lower BMI in smokers carrying the A-allele [25]. The other SNP, rs3743075 in CHRNA3, reported to be associated with lower SBP in twins [24] was not included in our directly genotyped data.

Possible mechanisms behind the associations

Co-morbidity between addictive behaviors of smoking, alcohol use and eating disorders has been reported [47], [48] suggesting a common underlying neurobiological mechanism. In this respect, it would not be surprising if same variants were associated with both smoking and BMI. In addition, nAChRs are plausible candidates for this common mechanism because they have been associated with different types of addictive behaviors: nicotine dependence [1], [2], [4], [6], [7], [9], [13], alcohol use [49] and overweight and obesity [50]. Further, evidence has been documented on an association between the CHRNA5-CHRNA3-CHRNB4 gene cluster and activity of the dorsal anterior cingulate-ventral striatum/extended amygdala circuit, which is an important circuit in modulating reward to drugs such as alcohol or smoking as well as food [51]. The associations of variation within genes encoding for nAChR subunits and BMI in the present study add evidence to this common mechanism hypothesis but future research is needed to better understand the biology behind.

The mechanisms by which the nAChRs regulate BP are also unclear. The genetic variants were still found to be associated with SBP even after adjustment for BMI suggesting a direct effect (i.e. not mediated via BMI) of these loci on SBP. In line with this idea, the genetic variants associated with SBP were different from those associated with BMI in our sample and they seemed to be nearly independent based on the low r2 between the variants associated with SBP and BMI. This suggests that the CHRNA5-CHRNA3-CHRNB4 gene cluster regulates SBP through biological mechanisms that partly differ from those regulating BMI. One plausible mechanism is that the nicotinic receptors modulate SBP by influencing the autonomic function. In this regard, Rana et al. showed that polymorphisms in the region were associated with catestatin, a peptide regulating cardiac function and blood pressure. The catestatin peptide fragment of the endogenous catecholamine secretory vesicle protein chromogranin A (CHGA) acts as an antagonist of the nicotinic cationic channels and regulates catecholamine release. This led the authors to suggest that nAChRs may regulate BP via mediating release of catestatin. Indeed, the same variants associated with catestatin were also associated with SBP, and cathecolamine levels. [24] We were not able to test this or other biologically plausible pathways due to the lack of appropriate biological markers measured in the cohort.

In the present study the associations between variation in the gene cluster and BMI/BP appeared to be more pronounced in smokers, suggesting that smoking would be a modifier of these associations. However, after adjustment for multiple testing, there was no sufficient evidence for interactive effects. The shortage of evidence for interaction might be due to limited power of our study to detect gene-environment interactions. Previous studies have reported associations between polymorphisms in the gene cluster and BMI manifest only in smokers [13], [25] or in alcohol consumers [49]. For BP, an interaction with smoking has not been reported before. Rana et al. showed associations between polymorphisms in the gene cluster and SBP adjusted for smoking, but they did not test for interactive effects in their study [24]. Although it seems that smoking is a favorable factor when carrying certain alleles as they are associated with lower BP and BMI especially in smokers, we cannot consider smoking beneficial in general because it is still an important risk factor for CVD. The identification of variants that lower BP and BMI in smokers, however, adds to the knowledge of pleiotropic effects of the variants in the gene cluster. Hopefully this could also be exploited in future development of medication by imitating the mechanisms through which nicotinic acetylcholine receptors influence BP and BMI in smokers.

Strengths and limitations

The analyses were performed in a population-based birth cohort including individuals of Northern European ancestry with a fairly homogeneous background. The study participants were of relatively young age (31 years) at the time of BMI and BP assessment, and as BP tends to increase with age, the effects could be stronger in older cohorts. In addition, age-varying associations are also possible [52]. Otherwise, previous GWA studies including NFBC have shown that the cohort is well comparable and the results are generalizable with other populations from similar ancestral background (e.g. [53], [54]).

However, replication is still required in populations of European ancestry, because although the sample size was moderately large compared to some of the previous studies, e.g. the 370 twins in the study for the gene cluster and SBP [24], the statistical evidence was not particularly strong but suggestive and especially regarding gene-environment interactions, the power was quite limited. The identification of gene-environment interactions is challenging, as was shown by Figueiredo et al. [55] in their recent study in which they compared three different methods for finding gene-environment interactions in an association study setting. Thomas et al. [56] further discuss the issues related to the search of gene-environment interactions, highlighting well-designed studies with careful measurement and efficient analysis of both genetic and environmental factors.

The replication of our suggestive results, and addressing further gene-environment, such as gender differences, or gene-gene interactions, would therefore require substantially larger numbers of individuals. Cohorts with certain biological markers will also be valuable in the search of specific pathways leading from nAChRs to cardiovascular risk factors.

Summary and conclusions

Our results suggest that variation in CHRNA5-CHRNA3-CHRNB4 is associated with cardiovascular risk factors and that smoking possibly acts as a modifier of these associations. Large collaborative efforts are needed to find further evidence for the associations, to test for their existence in other ancestral populations and age groups, and to understand the mechanisms by which the nAChRs affect BMI and BP.

Supporting Information

Table S1

Estimated associations between variants in the 15q25 region and SBP according to smoking status (non-smokers, light and heavy smokers) in the NFBC1966.


Table S2

Estimated associations between variants in the 15q25 region and DBP according to smoking status (non-smokers, light and heavy smokers) in the NFBC1966.


Table S3

Estimated associations between variants in the 15q25 region and BMI according to smoking status (non-smokers, light and heavy smokers) in the NFBC1966.


Table S4

Estimated associations between variants in the 15q25 region and SBP according to smoking status (never, former and current smokers) in NFBC1966.


Table S5

Estimated associations between variants in the 15q25 region and DBP according to smoking status (never, former and current smokers) in the NFBC1966.


Table S6

Estimated associations between variants in the 15q25 region and BMI according to smoking status (never, former and current smokers) in the NFBC1966.


Table S7

Estimated associations between variants in the 15q25 region and SBP according to gender in the NFBC1966.


Table S8

Estimated associations between variants in the 15q25 region and DBP according to gender in NFBC1966.


Table S9

Estimated associations between variants in the 15q25 region and BMI according to gender in the NFBC1966.


Table S10

Estimated haplotype frequencies and estimates from association analyses (Block 1) for SBP and BMI in smokers in the NFBC1966.


Table S11

Estimated haplotype frequencies and estimates from association analyses (Block 2) for DBP and BMI in smokers in the NFBC1966.



We thank the late Professor Paula Rantakallio (launch of NFBC1966 and initial data collection), Ms Sarianna Vaara (data collection), Ms Tuula Ylitalo (administration), Mr Markku Koiranen (data management), Ms Outi Tornwall and Ms Minttu Jussila (DNA biobanking).

Funding Statement

This work was supported by the Academy of Finland [project grants 104781, 120315, 129418, Center of Excellence in Complex Disease Genetics and Public Health Challenges Research Program (SALVE)], University Hospital Oulu, Biocenter, University of Oulu, Finland (75617), the European Commission [EURO-BLCS, Framework 5 award QLG1-CT-2000-01643], The National Heart, Lung and Blood Institute [5R01HL087679-02] through the SNP Typing for Association with Multiple Phenotypes from Existing Epidemiologic Data (STAMPEED) program [1RL1MH083268-01], The National Institute of Health/The National Institute of Mental Health [5R01MH63706:02], European Network of Genomic and Genetic Epidemiology (ENGAGE) project and grant agreement [HEALTH-F4-2007-201413], and the Medical Research Council, UK [G0500539, G0600705, PrevMetSyn/Public Health Challenges Research Program (SALVE)]. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Bierut LJ, Madden PA, Breslau N, Johnson EO, Hatsukami D, et al. (2007) Novel genes identified in a high-density genome wide association study for nicotine dependence. Hum Mol Genet 16 (1) 24–35 [PMC free article] [PubMed]
2. 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 (9) 1163–1171 [PMC free article] [PubMed]
3. 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 (4) 368–373 [PMC free article] [PubMed]
4. 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 (7) 926–933 [PMC free article] [PubMed]
5. Liu JZ, Tozzi F, Waterworth DM, Pillai SG, Muglia P, et al. (2010) Meta-analysis and imputation refines the association of 15q25 with smoking quantity. Nat Genet 42 (5) 436–440 [PubMed]
6. 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 (1) 36–49 [PMC free article] [PubMed]
7. 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 (17) 6848–6856 [PMC free article] [PubMed]
8. Saccone NL, Saccone SF, Hinrichs AL, Stitzel JA, Duan W, et al. (2009) 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 150B (4) 453–466 [PMC free article] [PubMed]
9. Saccone NL, Schwantes-An TH, Wang JC, Grucza RA, Breslau N, et al. (2010) Multiple cholinergic nicotinic receptor genes affect nicotine dependence risk in african and european americans. Genes Brain Behav 9 (7) 741–750 [PMC free article] [PubMed]
10. Saccone NL, Culverhouse RC, Schwantes-An TH, Cannon DS, Chen X, et al. (2010) Multiple independent loci at chromosome 15q25.1 affect smoking quantity: A meta-analysis and comparison with lung cancer and COPD. PLoS Genet 6 (8) e1001053. [PMC free article] [PubMed]
11. 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 (9) 1544–1552 [PMC free article] [PubMed]
12. 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 (12) 3517–3525 [PMC free article] [PubMed]
13. 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 (7187) 638–642 [PubMed]
14. Tobacco and Genetics Consortium (2010) Genome-wide meta-analyses identify multiple loci associated with smoking behavior. Nat Genet 42 (5) 441–447 [PMC free article] [PubMed]
15. Ducci F, Kaakinen M, Pouta A, Hartikainen AL, Veijola J, et al. (2011) TTC12-ANKK1-DRD2 and CHRNA5-CHRNA3-CHRNB4 influence different pathways leading to smoking behavior from adolescence to mid-adulthood. Biol Psychiatry 69 (7) 650–660 [PMC free article] [PubMed]
16. 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 (5) 616–622 [PMC free article] [PubMed]
17. 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 (7187) 633–637 [PubMed]
18. Liu P, Vikis HG, Wang D, Lu Y, Wang Y, et al. (2008) Familial aggregation of common sequence variants on 15q24-25.1 in lung cancer. J Natl Cancer Inst 100 (18) 1326–1330 [PubMed]
19. Pillai SG, Ge D, Zhu G, Kong X, Shianna KV, et al. (2009) A genome-wide association study in chronic obstructive pulmonary disease (COPD): Identification of two major susceptibility loci. PLoS Genet 5 (3) e1000421. [PMC free article] [PubMed]
20. Lindstrom JM (2003) Nicotinic acetylcholine receptors of muscles and nerves: Comparison of their structures, functional roles, and vulnerability to pathology. Ann N Y Acad Sci 998: 41–52 [PubMed]
21. Cheng PW, Lu PJ, Chen SR, Ho WY, Cheng WH, et al. (2011) Central nicotinic acetylcholine receptor involved in ca(2+) -calmodulin-endothelial nitric oxide synthase pathway modulated hypotensive effects. Br J Pharmacol 163 (6) 1203–1213 [PMC free article] [PubMed]
22. Li YF, LaCroix C, Freeling J (2009) Specific subtypes of nicotinic cholinergic receptors involved in sympathetic and parasympathetic cardiovascular responses. Neurosci Lett 462 (1) 20–23 [PMC free article] [PubMed]
23. Moore C, Wang Y, Ramage AG (2011) Nicotine's central cardiovascular actions: Receptor subtypes involved and their possible physiological role in anaesthetized rats. Eur J Pharmacol 668 (1–2) 177–183 [PubMed]
24. Rana BK, Wessel J, Mahboubi V, Rao F, Haeller J, et al. (2009) Natural variation within the neuronal nicotinic acetylcholine receptor cluster on human chromosome 15q24: Influence on heritable autonomic traits in twin pairs. J Pharmacol Exp Ther 331 (2) 419–428 [PubMed]
25. Freathy RM, Kazeem GR, Morris RW, Johnson PC, Paternoster L, et al. (2011) Genetic variation at CHRNA5-CHRNA3-CHRNB4 interacts with smoking status to influence body mass index. Int J Epidemiol 40 (6) 1617–1628 [PMC free article] [PubMed]
26. Leone A (2011) Does smoking act as a friend or enemy of blood pressure? let release pandora's box. Cardiol Res Pract 2011: 264894. [PMC free article] [PubMed]
27. Mineur YS, Abizaid A, Rao Y, Salas R, DiLeone RJ, et al. (2011) Nicotine decreases food intake through activation of POMC neurons. Science 332 (6035) 1330–1332 [PMC free article] [PubMed]
28. Audrain-McGovern J, Benowitz NL (2011) Cigarette smoking, nicotine, and body weight. Clin Pharmacol Ther 90 (1) 164–168 [PMC free article] [PubMed]
29. Newton-Cheh C, Johnson T, Gateva V, Tobin MD, Bochud M, et al. (2009) Genome-wide association study identifies eight loci associated with blood pressure. Nat Genet 41 (6) 666–676 [PMC free article] [PubMed]
30. Rantakallio P (1969) Groups at risk in low birth weight infants and perinatal mortality. Acta Paediatr Scand 193 (suppl) 1–71 [PubMed]
31. Järvelin MR, Sovio U, King V, Lauren L, Xu B, et al. (2004) Early life factors and blood pressure at age 31 years in the 1966 northern finland birth cohort. Hypertension 44 (6) 838–846 [PubMed]
32. Vartiainen E, Jousilahti P, Alfthan G, Sundvall J, Pietinen P, et al. (2000) Cardiovascular risk factor changes in finland, 1972–1997. Int J Epidemiol 29 (1) 49–56 [PubMed]
33. Tobin MD, Sheehan NA, Scurrah KJ, Burton PR (2005) Adjusting for treatment effects in studies of quantitative traits: Antihypertensive therapy and systolic blood pressure. Stat Med 24 (19) 2911–2935 [PubMed]
34. Sabatti C, Service SK, Hartikainen AL, Pouta A, Ripatti S, et al. (2009) Genome-wide association analysis of metabolic traits in a birth cohort from a founder population. Nat Genet 41 (1) 35–46 [PMC free article] [PubMed]
35. R Development Core Team. (2009) R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing.
36. 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 (3) 559–575 [PubMed]
37. Purcell S. PLINK 1.06. Accessed July 2012.
38. Buzkova P, Lumley T, Rice K (2011) Permutation and parametric bootstrap tests for gene-gene and gene-environment interactions. Ann Hum Genet 75 (1) 36–45 [PMC free article] [PubMed]
39. Liu N, Zhang K, Zhao H (2008) Haplotype-association analysis. Adv Genet 60: 335–405 [PubMed]
40. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics 21 (2) 263–265 [PubMed]
41. Zaykin DV, Westfall PH, Young SS, Karnoub MA, Wagner MJ, et al. (2002) Testing association of statistically inferred haplotypes with discrete and continuous traits in samples of unrelated individuals. Hum Hered 53 (2) 79–91 [PubMed]
42. Anderson MJ, Robinson J (2001) Permutation tests for linear models. Aust N Z J Stat 43 (1) 75–88
43. Anderson MJ (2001) Permutation tests for univariate and multivariate analysis of variance and regression. Can J Fish Aquat Sci 58: 626–639
44. Broms U, Wedenoja J, Largeau MR, Korhonen T, Pitkaniemi J, et al. (2012) Analysis of detailed phenotype profiles reveals CHRNA5-CHRNA3-CHRNB4 gene cluster association with several nicotine dependence traits. Nicotine Tob Res 14 (6) 720–733 [PMC free article] [PubMed]
45. 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 (11) 1039–1046 [PMC free article] [PubMed]
46. Hindorff LA, MacArthur J, Wise A, Junkins HA, Hall PN, et al. . A catalog of published genome-wide association studies. 2012(07/02).
47. Grucza RA, Bierut LJ (2006) Co-occurring risk factors for alcohol dependence and habitual smoking: Update on findings from the collaborative study on the genetics of alcoholism. Alcohol Res Health 29 (3) 172–178 [PubMed]
48. Petry NM, Barry D, Pietrzak RH, Wagner JA (2008) Overweight and obesity are associated with psychiatric disorders: Results from the national epidemiologic survey on alcohol and related conditions. Psychosom Med 70 (3) 288–297 [PubMed]
49. Landgren S, Engel JA, Andersson ME, Gonzalez-Quintela A, Campos J, et al. (2009) Association of nAChR gene haplotypes with heavy alcohol use and body mass. Brain Res 1305 Suppl: S72–9 [PubMed]
50. Kim J (2008) Association of CHRNA2 polymorphisms with overweight/obesity and clinical characteristics in a korean population. Clin Chem Lab Med 46 (8) 1085–1089 [PubMed]
51. Hong LE, Hodgkinson CA, Yang Y, Sampath H, Ross TJ, et al. (2010) A genetically modulated, intrinsic cingulate circuit supports human nicotine addiction. Proc Natl Acad Sci U S A 107 (30) 13509–13514 [PubMed]
52. Lasky-Su J, Lyon HN, Emilsson V, Heid IM, Molony C, et al. (2008) On the replication of genetic associations: Timing can be everything!. Am J Hum Genet 82 (4) 849–858 [PubMed]
53. International Consortium for Blood Pressure Genome-Wide Association Studies (2011) Ehret GB, Munroe PB, Rice KM, Bochud M, et al. (2011) Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk. Nature 478 (7367) 103–109 [PMC free article] [PubMed]
54. Speliotes EK, Willer CJ, Berndt SI, Monda KL, Thorleifsson G, et al. (2010) Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index. Nat Genet 42 (11) 937–948 [PMC free article] [PubMed]
55. Figueiredo JC, Lewinger JP, Song C, Campbell PT, Conti DV, et al. (2011) Genotype-environment interactions in microsatellite stable/microsatellite instability-low colorectal cancer: Results from a genome-wide association study. Cancer Epidemiol Biomarkers Prev 20 (5) 758–766 [PMC free article] [PubMed]
56. Thomas DC, Lewinger JP, Murcray CE, Gauderman WJ (2012) Invited commentary: GE-whiz! ratcheting gene-environment studies up to the whole genome and the whole exposome. Am J Epidemiol 175 (3) 203–7; discussion 208–9 [PMC free article] [PubMed]

Articles from PLoS ONE are provided here courtesy of Public Library of Science