Four SNPs in the CHRNA3/5 Alpha-Neuronal Nicotinic Acetylcholine Receptor Subunit Locus Are Associated with COPD Risk Based on Meta-Analyses

Abstract

Background

Several single nucleotide polymorphisms (SNPs) in an α-neuronal nicotinic acetylcholine receptor subunit (CHRNA3/5) were identified to be associated with chronic obstructive pulmonary disease (COPD) in a study based on a Norwegian population. However, results from subsequent studies have been controversial, particularly in studies recruiting Asians. In the present study, we conducted a comprehensive search and meta-analyses to identify susceptibility SNPs for COPD in the CHRNA3/5 locus.

Methods

A comprehensive literature search was conducted to find studies that have reported an association between SNPs in the CHRNA3/5 locus and COPD risk. Pooled odds ratios (ORs) with 95% confidence intervals (CIs) for each SNP were calculated with the major allele or genotype as the reference group. The influence of individual studies on pooled measures was assessed, in addition to publication bias.

Results

A total of 12 articles with 14 eligible studies were included in this analysis. Association between 4 SNPs in the CHRNA3/5 locus and COPD was evaluated and included rs1051730, rs8034191, rs6495309, and rs16969968. Significant associations between the 4 SNPs and COPD were identified under allele (rs1051730: OR = 1.14, 95%CI = 1.10–1.18; rs8034191: OR = 1.29, 95%CI = 1.18–1.41; rs6495309: OR = 1.26, 95%CI = 1.09–1.45; rs16969968: OR = 1.27, 95%CI = 1.17–1.39) and genotype models. Subgroup analysis conducted for rs1051730 showed a significant association between this SNP and COPD risk in non-Asians (OR = 1.14, 95%CI = 1.10–1.18), but not Asians (OR = 1.23, 95%CI = 0.91–1.67). Rs1051730 and rs6495309 were also significantly associated with COPD after adjusting for multiple variables, including age and smoking status.

Conclusion

Our results indicate that 4 SNPs in the CHRNA3/5 locus are associated with COPD risk. Rs1051730 was particularly associated with COPD in non-Asians, but its role in Asians still needs to be verified. Additional studies will be necessary to assess the effect of rs6495309 on COPD. Although rs1051730 and rs6495309 were shown to be independent risk factors for COPD, validation studies should be performed.

About the Authors:

Kai Cui

Affiliation: School of Public Health, Liaoning Medical University, Jinzhou, Liaoning, P.R. China

Xiaoyan Ge

Affiliation: School of Public Health, Liaoning Medical University, Jinzhou, Liaoning, P.R. China

Honglin Ma

* E-mail: [email protected]

Affiliation: School of Public Health, Liaoning Medical University, Jinzhou, Liaoning, P.R. China

Introduction

Tobacco smoking is a major risk factor for development of chronic obstructive pulmonary disease (COPD) [1]. Nicotine is the major reinforcing component of tobacco smoking and acts through neuronal nicotinic acetylcholine receptors (nAChRs) [2]. Studies have shown that genetic variants in the α-nAChR 3/5 subunit (CHRNA3/5) locus can influence nicotine dependence, smoking behavior, and lung cancer [3]–[6]. Pillai et al identified 2 susceptibility single nucleotide polymorphisms (SNPs) in the CHRNA3/5 locus, rs1051730 and rs8034191, as significant locations that may influence COPD in a recent genome-wide association (GWA) study [7]. Several subsequent studies have assessed the association between polymorphisms in this locus and COPD risk in different ethnicities and identified several other SNPs in the CHRNA3/5 locus as being associated with COPD risk. However, several studies have yielded contradictory results, particularly regarding Asian ethnic populations, possibly due to limited sample sizes and variations in study design. While most association studies recruit Caucasians, several SNPs, such as rs1051730 and rs8034191, are extremely rare in Asians [8]–[10]. Thus, meta-analyses are necessary to quantitatively synthesize the results. In the present study, we conducted a comprehensive search and meta-analyses to confirm the effect of CHRNA3/5 variations on COPD risk.

Materials and Methods

Search strategy

A comprehensive search was performed in several databases, including Pubmed, ISI Web of Science, CNKI (China National Knowledge Infrastructure), CBM (China Biology Medical Literature database), and Wangfang Data for potentially relevant studies published in English or Chinese. An upper date limit of Dec 1, 2013 was applied, and a lower date limit was not specified. A combination of the following Medical Subject Headings (MeSH) terms and key words were used: “pulmonary disease, chronic obstructive”, “lung disease, obstructive”, “COPD”, “chronic obstructive pulmonary disease”, “receptors, nicotinic”, “nicotinic receptor subunit”, “nAChR”, “CHRNA”, and “15q25”. Bibliographies of relevant articles were manually reviewed to identify additional studies that were not captured in the key word search.

Study eligibility

Studies that assessed the association between any SNP in the CHRNA3/5 locus and COPD risk and reported genotype count, allele frequency, or odds ratios (ORs) with 95% confidence intervals (CIs) were included in this study. Case-control, cohort, or population-based studies were eligible, although studies without a control group or based on family or sibling pairs were excluded. The study with the largest sample size was included when 2 or more studies contained the same data set. All potentially available articles were independently identified by 2 reviewers to determine study eligibility. Consensus with a third reviewer was used to resolve disagreements regarding eligibility.

Data extraction and quality assessment

Extracted data included the following: name of first author, publication year, country in which the study was conducted, ethnicity of the study population, type of study design, source of control subjects, mean age, male percentage, smoking status, genotyping methods, definition of COPD, matched methods, total number of cases and controls, genotype count, allele frequency, and ORs with 95%CIs adjusted for multiple variables.

The Newcastle-Ottawa quality assessment scale was used to evaluate the quality of each study. Assessment for case-control studies was based on 3 major components, including selection of cases and controls, comparability of cases and controls, and exposure. Prospective cohort studies were evaluated on another scale consisting of 3 components, including selection of exposed and non-exposed populations, comparability of cohorts, and outcome. The assessment was conducted independently by 2 reviewers, and disagreements were resolved by consensus with a third reviewer.

Statistical analysis

SNPs reported by more than 3 independent studies underwent further meta-analysis. Fisher's exact test was used to analyze deviations from Hardy-Weinberg equilibrium (HWE) for each genotype distribution in the control group. The inverse variance-weighted mean was used to calculate pooled ORs with 95%CIs for each SNP to represent association with risk of COPD under allele, homozygote, and heterozygote genotype models. The major allele or genotype was used as the reference group. The I² statistic was used to evaluate heterogeneity among studies, which describes the proportion of total variation attributable to between-study heterogeneity. I²>50% was considered significant statistical heterogeneity, and pooled ORs were calculated using the DerSimonian and Laird random effect model. For I²<50%, the fixed effect model was used. If needed, subgroup analyses were conducted to explore the source of heterogeneity. One-way sensitivity analyses in which a single study was omitted were employed to assess robustness of the pooled estimators in response to removal of an individual study. When the point estimate of the omitted analysis was outside the 95%CI of the pooled measures, the individual study was suspected of excessive influence. Begg's test was conducted to estimate publication bias and displayed as a funnel plot. All statistical analyses were performed with STATA version 9.2 (Stata Corporation, College Station, TX, USA). All reported probabilities were two-tailed, with P<0.05 considered statistically significant.

Results

Characteristics of included studies

A study identification flow diagram is shown in Figure 1. A total of 12 articles [7], [9], [11]–[20] consisting of 14 eligible studies satisfied the inclusion criteria. Only 1 article was published in Chinese [14]. Four SNPs were studied, including 9 studies on rs1051730 [11]–[13], [15]–[18], [20], 5 on rs8034191 [7], [13], [18], [20], 4 on rs6495309 [9], [14], [17], and 3 on rs16969968 [11], [18]–[19]. A total of 6 studies assessed associations between SNPs in the CHRNA3/5 locus and COPD risk in Asian populations [9], [14]–[15], [17]–[18]. A total of 10 studies were case-control, while all others were population-based. Most of the studies used the Global Initiative for Chronic Obstructive Lung Disease (GOLD) or American Thoracic Society (ATS) to diagnose or define COPD patients. General characteristics, quality assessments for each study, and allele and genotype distributions of each SNP are shown in Tables S1 and S2.

[Figure omitted. See PDF.]

Figure 1. Study identification flow diagram.

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

Quantitative synthesis

As shown in Table S2, none of the studies deviated from HWE for each genotype distribution in the control group. For rs1051730, 9 studies were included, and the pooled OR with 95%CI was calculated with the C allele as the reference group. The I² was 14.0%, and a fixed effect model was used. Allele distribution was significantly different between the groups, with an OR of 1.14 for the T vs C allele (95%CI = 1.10–1.18, P = 0.000) (Figure 2). Considering the genetic diversity between different ethnic populations, subgroup analysis was employed according to the ethnicity of the study population (Asian and non-Asian). In non-Asians, the T allele was associated with a significantly increased risk of COPD (OR = 1.14, 95%CI = 1.10–1.18, P = 0.000). However, a statistically significant difference was not observed in Asians (OR = 1.23, 95%CI = 0.91–1.67, P = 0.18).

[Figure omitted. See PDF.]

Figure 2. Forest plot for rs1051730 T vs C.

The fixed effect model was used in both overall and subgroup analyses. Each box represents the OR point estimate, and its area is proportional to the weight of the study. Diamond and broken lines represent the overall estimate, with 95%CI represented by the width. The solid vertical line is set at the null value (OR = 1.0).

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

For rs8034191, 5 studies were included. The I² was 0.0%, and the fixed effect model was used to calculate association of the C allele with COPD risk compared with the T allele. The C allele was significantly associated with the risk of COPD (OR = 1.29, 95%CI = 1.18–1.41, P = 0.000) (Figure 3).

[Figure omitted. See PDF.]

Figure 3. Forest plot for rs8034191 C vs T allele.

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

For rs6495309, 4 studies were included, and all were based on an Asian population. Marginal heterogeneity was observed (I² = 50.4%), and a random effect model was used to calculate association of the C allele with COPD risk compared with the T allele. The C allele was associated with a significantly increased risk of COPD (OR = 1.26, 95%CI = 1.09–1.45, P = 0.001) (Figure 4).

[Figure omitted. See PDF.]

Figure 4. Forest plot for rs6495309 C vs T allele.

https://doi.org/10.1371/journal.pone.0102324.g004

For rs16969968, 3 studies were included. The I² was 0.0%, and the OR for the A allele regarding risk of COPD compared with the G allele was 1.27 (95%CI = 1.17–1.39, P = 0.000) (Figure 5).

[Figure omitted. See PDF.]

Figure 5. Forest plot for rs16969968 C vs T allele.

https://doi.org/10.1371/journal.pone.0102324.g005

ORs with 95%CIs for heterozygote and homozygote genotypes were also calculated (Table 1). Under both heterozygote or homozygote genotype models of each SNP, results were similar to those of the allele models, except for the rs6495309 heterozygote genotype. The CT genotype increased the risk of COPD, but the difference was not significant (OR = 1.12, 95%CI = 0.85–1.49, P = 0.424).

[Figure omitted. See PDF.]

Table 1. ORs with 95%CIs for heterozygote and homozygote genotypes of each SNP.

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

To adjust for other confounding factors, such as smoking, age, and gender, we conducted meta-analyses for eligible studies that reported ORs with 95%CIs adjusted for multiple variables. For rs1051730 and rs6495309, 3 studies were included in the meta-analyses. Pooled ORs with 95%CIs were 1.11 (1.01–1.20) and 1.31 (1.08–1.54), respectively (Table 2).

[Figure omitted. See PDF.]

Table 2. Pooled ORs with 95%CIs for multiple variable adjustment.

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

Sensitivity analyses and publication bias

No strong influence of any single study was detected after it was omitted, except for the rs1051730 heterozygote genotype (Figure 6). After the 2012 study by Kaur-Knudsen et al was omitted, the pooled OR was outside the 95%CI of the pooled measures, suggesting that this study had an excessive influence on the results. The pooled OR increased after omission of this individual study (OR = 1.24, 95%CI = 1.13–1.38, P = 0.000). Considering the small number of studies included in the meta-analyses, Begg's test was conducted to assess potential publication bias when more than 4 studies were included. No obvious evidence of publication bias was detected for any of the allele, homozygote, and heterozygote models in each SNP, except for rs6495309 (Table 3).

[Figure omitted. See PDF.]

Figure 6. Sensitivity analysis for the rs1051730 heterozygote genotype (CT vs CC).

https://doi.org/10.1371/journal.pone.0102324.g006

[Figure omitted. See PDF.]

Table 3. Publication bias in meta-analyses for each inheritance model.

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

Discussion

In 2011, Zhang et al conducted a meta-analysis to assess the effect of rs1051730 in the CHRNA3/5 locus on COPD [21] based on data from published studies and 2 GWA studies (Bergen and COPCETI cohort [7], [22]). The results indicated that rs1051730 was a susceptibility variant for COPD. However, studies included in the Zhang et al meta-analyses were all based on Caucasians. Several subsequent association studies recruiting Asians have been conducted, but results have been controversial and inconclusive. It may be because polymorphism frequencies, such as rs1051073 and rs8034191, are lower in Asians than in Caucasians [8]–[10]. This implies differences in genetic susceptibility between the 2 ethnic populations and suggests that there might be other variants in the CHRNA3/5 locus that could influence COPD. In this study, we conducted a comprehensive search for studies about genetic variants in the CHRNA3/5 locus based on both Caucasian and Asian populations to assess the association of these polymorphisms with COPD risk. COPD likely results from a cumulative effect of environmental and genetic factors and gene-by-environment interactions [23]. The subjects in recent studies are former or current smokers, and as such it is difficult to identify whether variants in the CHRNA3/5 locus have a direct influence on COPD independent of smoking behavior or nicotine dependence.

While searching for eligible studies, we found 17 SNPs in the CHRNA3/5 locus that were reported to be associated with COPD, including rs1051730, rs8034191, rs6495309, rs16969968, rs578776, rs12910984, rs660652, rs1317286, rs569207, rs12914008, rs931794, rs938682, rs17486278, rs7180002, rs951266, rs684513, and rs588765. However, only the former 4 SNPs were studied in at least 5 individual studies. For each SNP included in this study, we compared allele and genotype distribution between cases and controls. Our results showed that the T allele in rs1051730 was associated with COPD risk, which was similar for both homozygote and heterozygote genotypes. This conclusion was consistent with the studies by Pillai et al and Zhang et al. In the present study, we found the polymorphism frequency was far lower in Asians. The minor allele frequency (MAF) in the non-Asian subgroup was greater than 0.48. The largest MAF was 0.031 from Guo's study in an Asian subgroup. We further evaluated the association between rs1051730 and COPD in different ethnic subgroups and showed that the T allele increased COPD risk, particularly in non-Asians. However, based on the results, rs1051730 does not appear to be associated with COPD in Asians. Because only 3 of the studies included an Asian subgroup and the MAF was lower than in non-Asians, additional studies are still needed to confirm the association between this polymorphism and COPD risk in Asians.

Rs8034191 is another susceptibility SNP identified in the study of Pillai et al. Our results showed that the rs8034191 C allele was associated with COPD risk under any of the inheritance models, including allele, homozygote, and heterozygote genotypes. Given that there was only one study based on an Asian population, subgroup analysis for different ethnicities was not conducted. Similar to rs1051730, the MAF of this SNP was 0.031 from Zhou's study based on Chinese subjects, while the lowest MAF in the studies with non-Asian participants was 0.51 (study of Kim et al), which was based on Caucasians. We suspect that there is a stronger association between rs8034191 and COPD risk in non-Asians than Asians, which should be confirmed through additional studies.

Based on the results of rs6495309 and rs16969968, an association between SNPs and COPD risk was detected. The rs6495309 C allele and rs16969968 A allele were considered risk alleles that increased the risk of COPD. The studies included in meta-analyses for rs6495309 were mainly based on Chinese or Korean populations, and none of the studies reported an association between rs6495309 and COPD in Caucasians or non-Asians, including the GWA studies, indicating that this SNP is likely an important susceptibility variant for COPD in Asians.

We tried to detect whether these SNPs in the CHRNA3/5 locus were associated with COPD risk by direct or indirect influence. We pooled the ORs with 95%CIs adjusted for multiple variables, including smoking status and pack-years of smoking. Eligible data were found in the rs1051730 heterozygote genotype and rs6495309 homozygote genotype. Although the number of included studies was small, the results still suggest that these variants are associated with COPD risk, independent of smoking behavior or nicotine dependence. This conclusion needs to be further confirmed with additional studies based on never smokers.

In sensitivity analyses, the study by Kaur-Knudsen et al had a strong influence on pooled measures under an inheritance model of rs1051730, with the pooled OR increasing after omitting this individual study. However, the pooled OR was still statistically significant at P = 0.05 and the results after omitting the study still confirmed an association between the polymorphism and COPD. Publication bias was detected under some inheritance models of rs6495309. In this meta-analysis, only 4 studies on the association of rs6495309 with COPD were included, and the bias could be caused by the small number of studies. According to our results, rs6495309 is associated with COPD, but future studies will need to confirm this conclusion regarding publication bias.

Based on our study, 4 SNPs in the CHRNA3/5 locus were associated with COPD risk. Among the 4 SNPs, the rs1051730 T allele was considered a risk allele, particularly in non-Asians. However, the effect of both rs1051730 in Asians and rs6495309 on COPD risk needs to be confirmed by further studies. Although rs1051730 and rs6495309 were found to be independent risk factors for COPD, this effect will need to be confirmed.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

General characteristics and quality assessment of each study.

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

(DOC)

Table S2.

Distribution of the SNP alleles and genotypes.

https://doi.org/10.1371/journal.pone.0102324.s002

(DOC)

Checklist S1.

PRISMA meta-analysis checklist.

https://doi.org/10.1371/journal.pone.0102324.s003

(DOCX)

Author Contributions

Conceived and designed the experiments: KC XG. Performed the experiments: KC XG HM. Analyzed the data: KC XG HM. Contributed reagents/materials/analysis tools: KC XG HM. Wrote the paper: KC XG HM.

Citation: Cui K, Ge X, Ma H (2014) Four SNPs in the CHRNA3/5 Alpha-Neuronal Nicotinic Acetylcholine Receptor Subunit Locus Are Associated with COPD Risk Based on Meta-Analyses. PLoS ONE 9(7): e102324. https://doi.org/10.1371/journal.pone.0102324

References

1. Hogg JC, Timens W (2009) The pathology of chronic obstructive pulmonary disease. Annu Rev Pathol 4: 435–459.

2. Picciotto MR, Kenny PJ (2013) Molecular mechanisms underlying behaviors related to nicotine addiction. Cold Spring Harb Perspect Med 3: a012112.

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

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

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

6. Amos CI, Wu X, Broderick P, Gorlov IP, Gu J, et al. (2008) Genome-wide association scan of tag SNPs identifies a susceptibility locus for lung cancer at 15q25.1. Nat Genet 40: 616–622.

7. 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: e1000421.

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

9. Lee JY, Yoo SS, Kang HG, Jin G, Bae EY, et al. (2012) A functional polymorphism in the CHRNA3 gene and risk of chronic obstructive pulmonary disease in a Korean population. J Korean Med Sci 27: 1536–1540.

10. Niu X, Chen Z, Shen S, Liu Y, Zhou D, et al. (2010) Association of the CHRNA3 locus with lung cancer risk and prognosis in Chinese Han population. J Thorac Oncol 5: 658–666.

11. Young RP, Hophins RJ, Hay BA, Epton MJ, Black PN, et al. (2008) Lung cancer gene associated with COPD: triple whammy or possible confounding effect? European Respiratory Journal 32: 1158–1164.

12. Kaur-Knudsen D, Bojesen SE, Tybjrg-Hansen A, Nordestgaard BG (2011) Nicotinic Acetylcholine Receptor Polymorphism, Smoking Behavior, and Tobacco-Related Cancer and Lung and Cardiovascular Diseases: A Cohort Study. Journal of Clinical Oncology 29: 2875–2882.

13. Kim DK, Hersh CP, Washko GR, Hokanson JE, Lynch DA, et al.. (2011) Epidemiology, radiology, and genetics of nicotine dependence in COPD. Respiratory Research 12.

14. Du M, Zhang W (2012) CHRNA3 gene promoter polymorphism with susceptibility to chronic obstructive pulmonary disease in smokers. International Journal of Respiration 32: 641–647.

15. Guo Y, Gong Y, Pan CM, Qian YR, Shi GC, et al.. (2012) Association of genetic polymorphisms with chronic obstructive pulmonary disease in the Chinese Han population: a case-control study. Bmc Medical Genomics 5.

16. Kaur-Knudsen D, Nordestgaard BG, Bojesen SE (2012) CHRNA3 genotype, nicotine dependence, lung function and disease in the general population. European Respiratory Journal 40: 1538–1544.

17. Yang L, Qiu F, Lu X, Huang D, Ma G, et al. (2012) Functional polymorphisms of CHRNA3 predict risks of chronic obstructive pulmonary disease and lung cancer in Chinese. PLoS One 7: e46071.

18. Zhou HX, Yang J, Li DX, Xiao J, Wang B, et al. (2012) Association of IREB2 and CHRNA3/5 polymorphisms with COPD and COPD-related phenotypes in a Chinese Han population. Journal of Human Genetics 57: 738–746.

19. Gabrielsen ME, Romundstad P, Langhammer A, Krokan HE, Skorpen F (2013) Association between a 15q25 gene variant, nicotine-related habits, lung cancer and COPD among 56 307 individuals from the HUNT study in Norway. European Journal of Human Genetics 21: 1293–1299.

20. Hoesein FAAM, Wauters E, Janssens W, Groen HJM, Smolonska J, et al.. (2013) Variants in the 15q24/25 Locus Associate with Lung Function Decline in Active Smokers. PLoS One 8.

21. Zhang J, Summah H, Zhu YG, Qu JM (2011) Nicotinic acetylcholine receptor variants associated with susceptibility to chronic obstructive pulmonary disease: a meta-analysis. Respir Res 12: 158.

22. Lambrechts D, Buysschaert I, Zanen P, Coolen J, Lays N, et al. (2010) The 15q24/25 susceptibility variant for lung cancer and chronic obstructive pulmonary disease is associated with emphysema. Am J Respir Crit Care Med 181: 486–493.

23. Rabe KF, Hurd S, Anzueto A, Barnes PJ, Buist SA, et al. (2007) Global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease: GOLD executive summary. Am J Respir Crit Care Med 176: 532–555.