Introduction {#Sec1}
============

Chronic obstructive pulmonary disease (COPD) is the third cause of death in the US after cancers and cardiovascular diseases^[@CR1]^ and is among the leading causes of hospitalization in industrialized countries^[@CR2],[@CR3]^. It was recently estimated that the absolute number of COPD cases in developed countries will increase by more than 150% from 2010 to 2030^[@CR4]^, yet there is no curative therapy for COPD^[@CR1],[@CR5]^. A lack of understanding of the molecular mechanisms in the pathogenesis of COPD has hampered efforts to develop new biomarkers and effective therapies.

COPD has multiple risk factors and complex etiology, where cigarette smoking and genetic susceptibility are among the main risk factors^[@CR1],[@CR6]^. Tobacco smoking accounted for about 5.1 million deaths globally in 2004, and it is observed recent increases in smoking prevalence in developing countries^[@CR7]^. But only 20--25% of smokers develop clinically significant airflow obstruction^[@CR8]^. Smoking behavior is partially genetically determined, and at least genome-wide significant associated loci were identified in European ancestry subjects^[@CR9],[@CR10]^, where the strongest and most consistent association reported is at the 15q25.1 locus (CHRNA3-CHRNA5-CHRNB4 gene cluster). Genes associated with smoking quantity were enriched for cholinergic receptors, sensory perception of smell, and Retinoid Binding^[@CR11]^.

In parallel, there is strong evidence that genetic variants (both common and rare) contribution to COPD and pulmonary function. Candidate gene and genome-wide association studies (GWAS) have identified genetic variants associated with COPD^[@CR12]--[@CR16]^. The latest GWAS was performed by the International COPD Genetics Consortium (ICGC)^[@CR17]^ identified 22 COPD susceptibility loci at genome-wide significance. Exome sequencing study found rare variants on CHRNA3, CHRNA5, and CHRNB4 genes were associated with COPD^[@CR18]^. Lung function (e.g. FEV~1~, FVC and FEV~1~/FVC ratio) are traits closely related to COPD, in fact, lung function parameters are often used in defining COPD cases^[@CR6]^. Three key parameters are commonly used in characterizing lung function: expired volume in 1st second (FEV~1~), forced vital capacity (FVC) and FEV~1~/FVC. FVC is the volume of air that can forcibly be blown out after full inspiration^[@CR19]^, and is the most basic maneuver in spirometry tests. FEV~1~ is the volume of air exhaled in the first second during forced exhalation after maximal inspiration^[@CR19]^. FEV~1~/FVC is the ratio of these two values, where the ratio ≥ 80% is considered normal^[@CR19]^. To date, 97 independent genetic loci were identified influencing lung function traits (FEV~1~, FVC and FEV~1~/FVC)^[@CR6]^.

Genetic pleiotropy is the phenomenon where a DNA variant influences multiple traits^[@CR20]^. Proposed by our group and others, available GWAS summary data can be used to detect genetic pleiotropy and pinpoint the variants, gene and pathways underlying the shared etiology^[@CR21]--[@CR23]^. Given smoking is a strong cause of COPD and lung function decline, genetic variants associated with smoking behavior or nicotine dependence (ND) are likely also in association with lung function (LF), which can be viewed as genetic pleiotropy. The genetic pleiotropy could have two possible mechanisms: (1) a genetic locus causes ND and then in turn causes LF decline (Fig. [1A](#Fig2){ref-type="fig"}); and (2) a genetic locus causes ND and LF decline independently (Fig. [1B](#Fig2){ref-type="fig"}). These two mechanisms can be distinguished by testing the association between genetic locus and smoking adjusted LF (ie, LF~adj~). If mechanism (1) is true, the ND locus will show no association with LF~adj~ (Fig. [1C](#Fig2){ref-type="fig"}).Figure 1Models of genetic pleiotropy. Figure 2QQplot showed SNPs associated with nicotine dependence are enriched for small p values in lung function GWAS. The y axis of the plot of the top row (plot A, B, C and D) denotes the observed p value of smoking adjusted FEV~1~ in UKBB study; the y axis of the plot of the middle row (plot E, F, G and H) denotes the observed p value of smoking adjusted FVC in UKBB study; the y axis of the plot of the bottom row (plot I, J, K and L) denotes the observed p value of smoking adjusted FEV~1~/FVC Ratio in UKBB study. In the 1st column (plot A, E and I), we investigated on SNPs associated with cigarettes per day (CPD) at ≤1e-2 (red dots), ≤1e-3 (blue dots), ≤1e-4 (green dots), and ≤1e-5 (black dots); in the 2nd column (plot B, F and J), we investigated on SNPs associated with ever smoking (EVERSMK) at ≤1e-2 (red dots), ≤1e-3 (blue dots), ≤1e-4 (green dots), and ≤1e-5 (black dots); in the 3rd column (plot C, G and K), we investigated on SNPs associated with former smoking (FORMER) at ≤1e-2 (red dots), ≤1e-3 (blue dots), ≤1e-4 (green dots), and ≤1e-5 (black dots); in the 4th column (plot D, H and L), we investigated on SNPs associated with log of smoking onset (ONSET) at ≤1e-2 (red dots), ≤1e-3 (blue dots), ≤1e-4 (green dots), and ≤1e-5 (black dots).

A recent GWAS meta-analysis employed 48,943 UK Biobank individuals in discovery phase and 95,375 individuals in follow-up phase, and increased the independent signals for lung function from 54 to 97^[@CR6]^. Importantly, the study carefully adjusted for smoking behavior, and among smokers, the pack-year was also adjustment^[@CR6]^. The GWAS results of smoking-adjusted lung function offer a unique opportunity to decipher the genetic pleiotropy between nicotine dependence and function as well as the etiological mechanisms.

Results {#Sec2}
=======

Genetic predisposition of nicotine dependence is associated with COPD risk {#Sec3}
--------------------------------------------------------------------------

We leverage the individual level phenotype (COPD affected status and smoking behavior) and genotype data of COPDgene cohort (Materials and Methods) to test association between genetic predisposition of nicotine dependence (quantitatively summarized as polygenic score, or PGS) and COPD risk. On each COPDgene study subject, we computed the PGS for ND CPD trait (i.e., PGS~ND-CDP~). The PGS~ND-CDP~ compiled at 1e-3 and 1e-4 p-value thresholds were associated with that COPD affected status at p-value 0.025 and 1.03e-6, respectively (Table [1](#Tab1){ref-type="table"}), indicating significant pleiotropy of ND and COPD. The positive association coefficient indicate genetic predisposition of higher cigarettes per day leads to higher COPD risk (Table [1](#Tab1){ref-type="table"}). PGS of former smoking (ie, PGS~ND-FORMER~) were also associated with COPD, and the negative coefficient suggested the genetic predisposition of smoking cessation reduced COPD risk. Next, we adjusted the smoking covariates and found the PGS~ND-CDP~ and PGS~ND-FORMER~ were still associated with COPD (Tables [2](#Tab2){ref-type="table"} and [3](#Tab3){ref-type="table"}). Importantly, the association p value and coefficient were little changed before and after adjustment, indicating an independent pleiotropy model (Fig. [1B](#Fig2){ref-type="fig"}), where the effect of PGS~ND-CDP~ and PGS~ND-FORMER~ on COPD risk was not mediated by smoking behavior of the study participants.Table 1Polygenic score of nicotine dependence is associated with COPD case/control status**PGS Training Trait1e-3 p-value threshold in PGS construction1e-4 p-value threshold in PGS construction**^**\#**^ **SNPsβp**.**value**^**\#**^ **SNPsβp**.**value**Cigarettes Per Day8910.0860.0254951160.1901.03E-06Ever Smoked11100.0820.0308471420.0440.257046Former Smoker858−0.1382.37E-04114−0.0690.078856LogOnset7970.0280.4331876−0.0230.430121Polygenic score (PGS) formulation is established based on TAG GWAS summary data using 1e-3 and 1e-4 p-value threshold on nicotine dependence (ND) traits. The ND PGS was computed on COPDgene dataset (sample size = 4903), and tested for association with COPD case-control status. ^\#^SNPs, the number of SNPs used in PGS computation, β, association coefficient (ie, log Odds Ratio). Table 2Polygenic score of nicotine dependence is associated with COPD case/control status adjusted for smoking status**PGS Training Trait1e-3 p-value threshold in PGS construction1e-4 p-value threshold in PGS construction**^**\#**^ **SNPsβp**.**value**^**\#**^ **SNPsβp**.**value**Cigarettes Per Day8910.0830.0314821160.1852.07E-06Ever Smoked11100.0900.0175211420.0480.216163Former Smoker858−0.1431.49E-04114−0.0660.091961LogOnset7970.0290.41949776−0.0230.440345Polygenic score (PGS) formulation is established based on TAG GWAS summary data using 1e-3 and 1e-4 p value threshold on nicotine dependence (ND) traits. The ND PGS was computed on COPDgene dataset (sample size = 4903), and tested for association with COPD case-control status, adjusted for smoking status (current smoker vs. non-smoker). ^\#^SNPs, the number of SNPs used in PGS computation, β, association coefficient (ie, log Odds Ratio). Table 3Polygenic score of nicotine dependence is associated with COPD case/control status adjusted for smoking duration**PGS Training Trait1e-3 p-value threshold in PGS construction1e-4 p-value threshold in PGS construction**^**\#**^ **SNPsβp**.**value**^**\#**^ **SNPsβp**.**value**Cigarettes Per Day8910.0850.03761160.1971.82E-06Ever Smoked11100.0720.0704471420.0570.163214Former Smoker858−0.1329.46E-04114−0.0670.106736LogOnset7970.0220.55550576−0.0190.544609Polygenic score (PGS) formulation is established based on TAG GWAS summary data using 1e-3 and 1e-4 p value threshold on nicotine dependence (ND) traits. The ND PGS was computed on COPDgene dataset (sample size = 4903), and tested for association with COPD case-control status, adjusted for smoking duration. ^\#^SNPs, the number of SNPs used in PGS computation, β, association coefficient (ie, log Odds Ratio).

Nicotine dependence genome-wide significant loci and their association with lung function (LF) {#Sec4}
----------------------------------------------------------------------------------------------

Firstly, we retrieved six independent ND GWAS loci of genome-wide significance (p \< 5e-8) in European ancestry^[@CR24]^ from the NHGRI-EBI catalog^[@CR25]^, and examined the association between these loci and smoking adjusted lung function (LF~adj~) in UKBB meta-analysis cohort (Table [4](#Tab4){ref-type="table"}). In UKBB study, the smoking status was adjusted with the regression model in the form of indicator variables denoting current smoker, and denoting former smokers; as well as quantitative variable denoting pack-year (in smokers). At nominal p value (≤0.05), two ND loci showed association with at least one LF~adj~ traits (FEV~1~, FVC, and FEV~1~/FVC Ratio). The lead SNP, rs1051730, of the 15q25.1 locus was strongly associated with all three LF~adj~ trait, and its association with smoking adjusted FEV~1~ reached genome-wide significant (p = 5e-8), further, the 19q13.2 locus, known in association with COPD^[@CR26]^, is also associated with smoking adjusted FEV~1~ (Table [4](#Tab4){ref-type="table"}), suggesting the 15q25.1 and 19q13.2 loci could influence ND and LF independently (Fig. [1B](#Fig2){ref-type="fig"}). In contrast, four genome-wide significant ND loci (10q23.32, 8p11.21, 11p14.1, 9q34.2) showed no association with LF~adj~ traits (Table [4](#Tab4){ref-type="table"} and Fig [S1](#MOESM1){ref-type="media"}), suggesting causal pleiotropy model (Fig. [1 A](#Fig2){ref-type="fig"}), where the genetic variants influence LF, mediated by smoking behavior.Table 4Nicotine dependence loci and in lung function.**Locus\*Lead SNPChrPositionReported GenesND PAssociation with LFFEV** ~**1**~**FVCRATIO**15q25.1rs10517301578894339CHRNA33.0E-73**YYY**10q23.32rs13296501093348120LOC1001889476.0E-10NNN8p11.21rs6474412842550498CHRNB3, CHRNA61.0E-08NNN19q13.2rs37338291941310571EGLN2, CYP2A6, RAB4D1.0E-08**Y**N**Y**11p14.1rs62651127679916BDNF2.0E-08NNN9q34.2rs30253439136478355DBH4.0E-08NNN^\*^Genome-wide significant loci of nicotine dependence in European ancestry, 15q25.1^[@CR32]^, 10q23.32^[@CR32]^, 19q13.2^[@CR10],[@CR32]^, 11p14.1^[@CR32]^, 9q34.2^[@CR32]^, 8p11.21^[@CR10]^ were summarized by NHGRI-EBI GWAS catalog. The reported lead SNP was checked for association with smoking adjusted lung function in UKBB GWAS meta-analysis^[@CR6]^. ND, nicotine dependence; LF, lung function.

Identification of SNPs associated with lung function and nicotine dependence at suggestive threshold {#Sec5}
----------------------------------------------------------------------------------------------------

There is a general consensus among GWAS studies that a p-value less than 5e-8 corresponds to genome-wide significance^[@CR24]^. NHGRI-EBI GWAS catalog focuses on highly significant loci, but true association loci may only show suggestive p values given the limited power in GWAS. We leverage the summary data of with UKBB smoking adjusted LF (LF~adj~) and TAG ND (Materials and Methods). Applied various p thresholds (ranging from 1e-8 to 1e-2), we identified SNPs in associated with LF~adj~ or ND (Table [S1](#MOESM1){ref-type="media"}). The UKBB LF~adj~ GWAS yield many SNPs highly significantly associated with smoking adjusted FEV~1~, FVC and FEV~1~/FVC Ratio. On the other hand, the TAG ND study captured strong genetic signals for CPD (cigarettes per day) and moderate signal for former smoking (Table [S1](#MOESM1){ref-type="media"}). But no SNP associated ever smoking or log-onset at p \< 1e-7 level, indicating these two traits either are not controlled by genetic factors or the TAG study did not have sufficient power (Table [S1](#MOESM1){ref-type="media"}).

Among the \~2 million SNPs tested in both LF~adj~ and ND GWAS studies, the proportion of associated SNPs at any given p-value threshold was greater than alpha level (Table [5](#Tab5){ref-type="table"}). For example, at p-value ≤ 1e-3 threshold, 11,716 SNPs and 2,791 SNPs were associated with FEV~1~ and CPD, respectively, and 71 SNPs were associated with both trait, which is substantially greater than random chance (enrichment fold = 4.61). In fact, we observed excess overlapping of FEV~1~ and CPD GWAS SNPs were across multiple p-value thresholds (Table [5](#Tab5){ref-type="table"} and Fig. [2](#Fig1){ref-type="fig"}). Further, at p value ≤ 1e-3 threshold, 139 SNPs were associated with both FEV~1~ and former smoking (Table [6](#Tab6){ref-type="table"}), 11.8 folds of enrichment than random chance. While, GWAS signals of ever smoking and smoking onset showed little or very modest overlap with LF~adj~ GWAS (Tables [S2](#MOESM1){ref-type="media"} and [S3](#MOESM1){ref-type="media"}). The SNPs associated with at least one of the nicotine dependence and lung function traits (cigarettes per day, former smoking, FEV~1~, FVC, FEV~1~/FVC) at p value ≤ 1e-3 level were listed in Table [S4](#MOESM2){ref-type="media"}; while SNPs associated with all of the five nicotine dependence and lung function traits were detailed in Table [S4](#MOESM2){ref-type="media"}.Table 5Overlap of GWAS signals between nicotine dependence (cigarettes per day) and lung function traits.**Lung Function TraitsGWAS P value thresholdEnrichment FoldN of overlap SNPsN of overlap SNPs of divergent risk direction% SNPs of divergent risk directionp-value of divergent risk direction**FEV~1~1.0E-0626126261006.32E-081.0E-0512726261006.32E-081.0E-0454.942421001.35E-121.0E-034.61716591.554.34E-131.0E-021.6782645154.601.29E-21.0E-011.0329,06715,40452.999.00E-24FVC1.0E-06---00------1.0E-05---00------1.0E-0419.113131003.91E-041.0E-033.50494795.921.17E-111.0E-021.8186949156.502.40E-041.0E-011.0729,86515,84053.044.47E-25FEV~1~/FVC Ratio1.0E-06---00------1.0E-0563.9014141002.08E-041.0E-0414.6014141002.08E-041.0E-032.2839391009.80E-121.0E-021.3166940760.845.23E-081.0E-011.0328,46814,47450.846.70E-03^\*^Consistent direction allele, the SNPs where the same allele association with higher cigarettes per day and higher lung function trait (e.g. FEV~1~). Table 6Overlap of GWAS signals between nicotine dependence (former smoking) and lung function traits.**Lung Function TraitsGWAS P value thresholdEnrichment FoldN. overlap SNPsN of overlap SNPs of divergent risk direction% SNPs of divergent risk directionp-value of divergent risk direction**FEV~1~1.0E-06---00------1.0E-0522.94001.53E-011.0E-0417.120003.59E-061.0E-0311.81932110.886.69E-301.0E-022.29120229024.135.74E-741.0E-011.07303281349044.484.54E-81FVC1.0E-06---00------1.0E-05---00------1.0E-048.248001.13E-021.0E-039.681442114.583.94E-181.0E-022.14109933630.576.23E-381.0E-011.06300501345344.772.26E-72FEV~1~/FVC Ratio1.0E-06---00------1.0E-0523.54001.53E-011.0E-0417.323004.77E-071.0E-039.3516831.793.49E-441.0E-021.98106634332.185.01E-311.0E-011.04289801392348.047.26E-11^\*^Consistent direction allele, the SNPs where the same allele association with higher odds of former smoking (ie, smoking cessation) and higher lung function trait (e.g FEV~1~).

Association direction of shared LF~adj~ and ND SNPs {#Sec6}
---------------------------------------------------

The shared LF~adj~/ND SNPs can be stratified into categories according whether the allele associated with higher LF~adj~ trait is also associated with higher ND trait or ND event odds (ie, consistent direction SNP) or vice versa (ie, divergent direction SNP). We stratified the shared LF~adj~/ND SNPs by directions, and found all three LF~adj~ traits primarily have divergent direction with CDP (Table [5](#Tab5){ref-type="table"}). At 1e-3 threshold, 71 SNPs are shared by FEV~1~ and CPD, where 65 SNPs (91.55%) with divergent directions (binomial p value = 4.34e-13), indicating alleles lead to higher smoking quantity (cigarettes per day) tends to associated with lower LF~adj~ traits. In contrast, the LF~adj~ traits primarily have consistent direction with former smoking (Table [6](#Tab6){ref-type="table"}). At 1e-3 threshold, 193 SNPs were shared by FEV~1~ and former smoking, and only 21 SNPs (10.88%) with divergent directions but 172 SNPs are of a consistent direction (binomial p-value = 4.34e-13), indicating alleles lead to smoking cessation tend to associated with higher LF~adj~ traits.

Further, we stratified the SNPs into the consistent direction SNP and divergent SNP subcategories, regardless of GWAS test p-values criteria (Tables [S5](#MOESM1){ref-type="media"}--[S8](#MOESM1){ref-type="media"}). We found the overlap of CPD and LF~adj~ GWAS signals primarily occur among the SNPs of divergent allele direction, and observed little overlap among the SNPs of consistent allele direction. For example, at a GWAS p-value \< 1e-3, SNPs of divergent allele direction for FEV~1~ and CPD showed substantial overlap (6.94 fold enrichment), in contrast, the overlap for SNPs of consistent allele direction were not different from random chance.

Gene ontology influenced by the shared SNPs of LF~adj~ and ND {#Sec7}
-------------------------------------------------------------

Shared LF~adj~/ND SNPs may undergo complicated pathways and lead to disease risk. Identifying the genes influenced by such shared SNPs would be the key step elucidating the SNPs' function and pathogenic pathways. To investigate potential functional impacts of LF~adj~/ND shared SNPs, we used a "relax" p-value threshold criterion of 1e-2 for GWAS. Among the 826 shared FEV~1~ - CPD SNPs, 375 are characterized by having consistent direction and 451 are characterized by divergent direction. Based on SNP position and dbSNP annotation^[@CR27]^, we identified genes located with on or near the share GWAS SNPs of LF~adj~ and ND, and then carried out the Gene Ontology (GO) cellular processes enrichment analysis using METACORE suite (Fig. [3](#Fig3){ref-type="fig"}). The top 10 enriched processes (p-value ≤ 1e-6) were regulation of synapse vesicles as well as regulation of respiratory systems, highly relevant to lung function and nicotine dependence.Figure 3Genes harbor or close to share SNPs between cigarettes per day and smoking adjusted FEV~1~ are enriched for certain Gene Ontology cellular processes.

Co-localization of 15q25.1 locus underlying LF~adj~ and ND {#Sec8}
----------------------------------------------------------

SNPs on 15q25.1 locus are associated with both LF~adj~ and ND at genome-wide significant level (Table [4](#Tab4){ref-type="table"}). Also, the CHRNA3-CHRNA5-CHRNB4 gene cluster in this locus is functional relevant to both traits. However, given that neighboring SNPs were often in tight linkage disequilibrium (LD), the overlap of GWAS signals do not guarantee that the disease risks of the two traits are caused by the same variant. Recently developed methods allow more advanced integration of GWAS summary data to co-localize GWAS signals^[@CR28]^. The co-localization methods^[@CR28]^ evaluated 5 hypotheses (Materials and Methods), where we were particularly interested in hypothesis 4 (H4: the two phenotypic traits were caused by the same SNP in the locus), and H4 posterior probability over 0.75 was considered as supporting evidence to the corresponding hypothesis^[@CR29],[@CR30]^. All the three LF~adj~ trait (FEV~1~, FVC and FEV~1~/FVC Ratio) were co-localized at 15q25.1 (Table [7](#Tab7){ref-type="table"}), indicating they are controlled by the same genetic variant. The cigarettes per day and former smoking traits were also co-localized with LF~adj~ (Table [8](#Tab8){ref-type="table"}). For example, between CPD and FEV~1~, the coloc H~4~ posterior probability was 0.908, and between former smoking and FEV~1~, the coloc H~4~ posterior probability was 0.942. The clear co-localization of LF~adj~ and ND at 15q25.1 suggested same genetic variant influence both nicotine dependence and lung function, through an independent pleiotropy model.Table 7Co-localization of lung function GWAS signal at 15q25.1 locus\*.**Trait 1Trait 2N SNPsPP**.**H0PP**.**H1PP**.**H2PP**.**H3PP**.**H4**FEV~1~FVC130400.0001490.06530.0002660.1154610.818824FEV~1~FEV~1~/FVC Ratio130414.85E-060.0021219.77E-050.0417580.956019FVCFEV~1~/FVC Ratio130400.0040640.0072360.0818520.1449940.761854^\*^chr15:77,801,394-79,801,394 (HG19). Five COLOC hypothesis (H0\~H4) were evaluated (Materials and Methods). Table 8Co-localization of lung function and nicotine dependence GWAS signal at 15q25.1 locus^\*^.**Trait 1Trait 2N SNPsPP**.**H0PP**.**H1PP**.**H2PP**.**H3PP**.**H4**FEV~1~Cigarettes Per Day17648.56E-321.46E-290.0005370.0910250.908438Ever Smoked17620.0053810.9209540.0004010.0686630.004601Former Smoker17589.76E-050.0167080.0005020.0849860.897706LogOnset17600.0056010.958570.0001720.0294060.006252FVCCigarettes Per Day17642.52E-291.24E-290.1578980.0770690.765033Ever Smoked17620.6221320.3066750.0463870.0228640.001942Former Smoker17580.0484710.0238890.2491480.1222380.556253LogOnset17600.6485520.3196650.01990.0098060.002077FEV~1~/FVC RatioCigarettes Per Day17641.17E-308.20E-300.0073710.0505540.942076Ever Smoked17620.1160320.8106570.0086520.0604390.00422Former Smoker17580.0020550.014360.0105650.0729120.900107LogOnset17600.1207590.8436420.0037050.025880.006014^\*^chr15:77,801,394-79,801,394 (HG19). Five COLOC hypothesis (H0\~H4) were evaluated (Materials and Methods).

Discussion {#Sec9}
==========

In this report, we leverage the latest large GWAS data sets and investigated the genetic pleotropic effect between nicotine dependence and respiratory outcomes (ie, lung function and COPD). It is known smoking is a major cause of reduced lung function and COPD^[@CR8]^. Smoking behavior is at least partially controlled by genetic factors^[@CR7],[@CR9]^. To date, the pleiotropic effects of genetic risk of smoking behavior and lung function/COPD have not been systematically explored.

Employing several analytical approaches, this paper addressed two questions, (1) whether the genetic predisposition of nicotine dependence influence COPD risk and lung function; and (2) the genetic pleiotropy follow causal or independent model (Fig. [1](#Fig2){ref-type="fig"}). On COPDgene cohort, we found the polygenic score of nicotine dependence (ie, PGS~ND~) was associated with COPD case/control status, demonstrating the genetic pleiotropy of the two conditions. We also investigated the association between PGS~ND~ and COPD while adjusting smoking behavior. Interestingly, the crude and adjusted results were very similar, indicating a mainly independent pleiotropy model. That is the shared genetic factors directly modify COPD risk, not mediated by influencing the individual's smoking behavior.

We zoomed in the known ND loci of genome-wide significance (Table [4](#Tab4){ref-type="table"}), and found both causal and independent pleiotropy models may exist in certain loci. Two genome-wide significant ND loci (15q25.1 and 19q13.2) were associated with LF~adj~ (smoking adjusted lung function), supporting independent pleiotropy model. While, four ND loci (10q23.32, 8p11.21, 11p14.1 and 9q34.2) were not associated with smoking adjusted lung function. At various p value threshold (1e-6 to 1e-1), we found the LF~adj~ traits share association SNP with cigarettes per day and former smoking substantially more than random chance, indicating a large number of genetic variants contribute to the genetic pleiotropy. Importantly, the lung function and cigarettes per day mainly share SNPs of divergent direction, meaning genetic predisposition of higher smoking dosage leads to lower lung function. In contrast, the lung function and former mainly share SNPs of consistent direction, meaning genetic predisposition of smoking cessation leads to higher lung function.

In summary, we used empirical data of largest cohorts to date and showed the genetic pleiotropy between nicotine dependence and COPD or lung function. The pleiotropic effect exist even COPD status or lung function is adjusted for smoking behavior. Further, we found the pleiotropic effect is attributable not only to the genome-wide significant loci, but also loci associated to ND and COPD/LF at suggestive p value (e.g. 1e-3), suggesting a large number of variants influence both ND and respiratory outcome, and among which many variants functions through independent genetic pleiotropy model.

Materials and Methods {#Sec10}
=====================

Genome-wide meta-analyses on nicotine dependence (ND) {#Sec11}
-----------------------------------------------------

The Tobacco and Genetics (TAG) Consortium conducted GWAS meta-analyses across 16 studies^[@CR10]^. We examined four carefully harmonized smoking phenotypes: smoking initiation (ever versus never been a regular smoker), age of smoking initiation, smoking quantity (number of cigarettes smoked per day, CPD) and smoking cessation (former versus current smokers) among people of European ancestry. Smoking cessation contrasted former versus current smokers, where current smokers reported at interview that they presently smoked and former smokers had quit smoking at least 1 year before interview. Smokers who had quit smoking for less than 1 year at interview were excluded from the analysis to minimize misclassification. Genotype imputation resulted in a common set of \~2.5 million SNPs that entered the GWAS and meta-analysis, where summary statistics were used in this report.

COPDGene dataset {#Sec12}
----------------

Individual level genotype and phenotype data of COPDgene study were retrieved from dbGap (accession: phs000179.v1.p1). COPD case/control status, indicator variable of current smoking, indicator variable of former smoking, quantitative variable of smoking duration, and SNP genotype data were available on 4,903 individuals. We conducted genotype imputation using HRC reference^[@CR31]^, and in total, 19,932,879 SNPs of high quality score entered the current study. In polygenic score analysis, the smoking status variable (indicator and quantitative variables) were adjusted within the regression model.

UK Biobank Lung Function (LF) GWAS {#Sec13}
----------------------------------

Recently we reported a large GWAS study on lung function using the UK biobank samples^[@CR6]^. Genome-wide association analyses of forced expired volume in 1 s (FEV~1~), forced vital capacity (FVC) and FEV~1~/FVC were undertaken in 48,943 individuals from the UK BiLEVE study^[@CR7]^ who were selected from the extremes of the lung function distribution in UK Biobank (total n = 502,682). Association tests were conducted on 27,624,732 variants, where linear regression of age, age^[@CR2]^, sex, height, the first ten principal components of genetic ancestry and pack-years of smoking (in smokers), and summary statistics were used in this report.

Polygenic Score {#Sec14}
---------------

We analyzed GWAS summary statistics data from the TAG nicotine dependence study, COPDgene dataset case control status, and COPDgene dataset smoking related covariate: smoking status (binary variable) and smoking duration (continuous variable). We computed the nicotine dependence polygenic score (PGS) on each COPDgene study subject in following steps: (1) identify shared SNPs in TAG GWAS summary data and COPDgene imputed genotype; (2) align alleles strands to the 1000 G panel (hg19), and adjust β coefficients TAG nicotine dependence GWAS accordingly; (3) filter TAG GWAS data by p value threshold (1e-3 and 1e-4 as shown in Table [5](#Tab5){ref-type="table"}) and prune the SNPs by linkage disequilibrium (LD) based on 1000 G EUR reference; lastly (4) for each ND traits, and for each p value threshold, we computed a PGS for every subject as a linear combination of the imputed doses of the selected coefficients. Then we tested the association between nicotine dependent PGS and COPD case/control status using a logistic regression model with or without adjusting for smoking covariates.

Identification of SNPs associated with both LF and ND {#Sec15}
-----------------------------------------------------

The effect size attributable to individual genetic variants for a given complex disorder is typically modest, suggesting that individual genetic variants may only explain a very small amount of the genetic risk and heritability of complex disorders^[@CR21],[@CR31]^. Therefore, genetic contributions to complex conditions such as LF and ND are likely derived from a large number of genetic causal variants, each contributing a small genetic risk. We surveyed a number of p value thresholds, and identified SNPs that are associated with LF~adj~ and ND in order to more comprehensively capture SNPs with modest effect sizes. For a shared SNP, we term it "consistent allele direction" if a specific allele that is associated with increased LF~adj~ traits and that specific SNP allele is also associated with higher ND trait value or event odds. We term "divergent allele direction" for SNP where a specific allele associated with increased lung function, and lower ND trait value or ND event odds.

Pathway Enrichment Analyses {#Sec16}
---------------------------

To further characterize the regulatory nature, enrichment analysis of the genes influenced by shared SNPs were performed using the METACORE integrated software suite (<http://thomsonreuters.com/metacore/>).

Co-localization of lung function and nicotine dependence GWAS top SNPs {#Sec17}
----------------------------------------------------------------------

LF~adj~ and ND GWAS results were used in co-localization analysis, which is performed using COLOC version 2.3--6 in R^[@CR28]^. Our analysis focuses on the 15q25.1 locus. Default priors of the software were used. In total, 5 hypotheses were evaluated. H0: No association with either trait 1 or trait 2; H1: Association with trait 1, not with trait 2; H2: Association with trait 2, not with trait 1; H3: Association with trait 1 and trait 2, two independent SNPs; H4: Association with trait 1 and trait 2, one shared SNP. Genes that demonstrated a high posterior probability of hypothesis 4 (PP.H4 \> 75%) indicate the disease risk and placenta gene expression were controlled by the same genetic variant; and genes that demonstrated a high posterior probability of hypothesis 3 (PP.H3 \> 75%) indicate the disease risk and placenta gene expression were controlled by distinct genetic variant at the locus.

Electronic supplementary material
=================================

 {#Sec18}

Supplementary Figures and Tables Table S4

**Electronic supplementary material**

**Supplementary information** accompanies this paper at 10.1038/s41598-017-16964-4.

**Publisher\'s note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

We thank Drs. Alison Goate and Manav Kapoor for the valuable discussions and suggestions. This work is partially supported by NIH 1R41DA042464--01, NIH 1U01HD079068-01, National Natural Science Foundation of China (Grant No. 21477087, 91643201) and by the Ministry of Science and Technology of China (Grant No. 2016YFC0206507).

Z.J., P.S., C.H., D.N.A. and H.K. are responsible for designing and conducting the study. Z.J., P.S., Y.N., D.N.A., and H.K. wrote the manuscript. All authors read and approved the final manuscript.

Competing Interests {#FPar1}
===================

The authors declare that they have no competing interests.
