Background {#Sec1}
==========

Chronic obstructive pulmonary disease (COPD) is a major cause of morbidity and mortality in the world and encompasses emphysema, chronic bronchitis, and small airways disease \[[@CR1], [@CR2]\].  The diagnosis of COPD is largely based on the presence of airflow obstruction, measured by the spirometric assessment (post-bronchodilator) of the ratio between forced expiratory volume in one second and forced vital capacity (FEV~1~/FVC).  The Global initiative for chronic Obstructive Lung Disease (GOLD) recommends to use a fixed cut-off for defining airflow obstruction, namely an FEV~1~/FVC ratio below 70% \[[@CR3]\], whereas the American Thoracic Society/European Respiratory Society (ATS/ERS) guidelines recommend to define airflow obstruction as FEV~1~/FVC below the lower limit of normal (LLN) \[[@CR4]\].  The LLN is a reference value based on sex, age, height and ethnicity and is calculated as the lower fifth percentile of a healthy reference population \[[@CR5]\].  There is a considerable controversy about which definition should be used in research and clinical practice, since both may lead to misclassifications \[[@CR5]--[@CR8]\].  This has important implications, since misclassifications may lead to inappropriate medication and therapies \[[@CR9], [@CR10]\]. 

It is generally accepted that both genetic susceptibility and environmental factors contribute to airflow obstruction.  Genetic variants associated with airflow obstruction have been identified by several genome-wide association studies (GWAS), but different definitions of airflow obstruction and populations were used.  \[[@CR11]--[@CR16]\] As an illustration, the case-control study including only smokers with \> 2.5 pack-years by Pillai et al.  used the fixed ratio (FEV~1~/FVC \< 70%) to define airflow obstruction, while the population based study including both ever- and never-smokers by Wilk et al.  used the lower limit of normal (LLN) \[[@CR15], [@CR16]\].  Only few regions were identified in both studies, namely the *CHRNA5/3* and *HHIP* regions.  We therefore aimed to assess the genetic overlap between the two definitions of airflow obstruction in the same individuals.  We stratified by smoking status to assess the overlap between the two airflow obstruction definitions in never- and ever-smokers separately.  We used the Lifelines Cohort Study as discovery sample and the Vlagtwedde-Vlaardingen study to replicate our observations.  In addition, genetic loci associated with both airflow obstruction definitions could indicate robust genetic associations with airflow obstruction, which could potentially be novel loci.  We therefore, as a secondary aim, validated the top overlapping single-nucleotide polymorphisms (SNPs) between the two airflow obstruction discovery analyses in an independent SNP validation sample and assessed if they were acting as expression quantitative trait loci (eQTLs) in a lung tissue sample. 

Materials and methods {#Sec2}
=====================

Study populations {#Sec3}
-----------------

To study the overlap between the two airflow obstruction definitions, all subjects with available genotypic data were included from the Dutch LifeLines Cohort Study (discovery sample) and the Vlagtwedde-Vlaardingen study (replication sample) \[[@CR17]--[@CR19]\].  In addition, subjects from the Vlagtwedde-Vlaardingen study and the three independent cohorts of the Rotterdam Study (RS I to III) were selected to validate the top overlapping SNPs from LifeLines (SNP validation sample), thereby increasing the SNP validation sample size \[[@CR20]\].  All subjects provided written informed consent and the studies were approved by local medical ethics committees.  Smoking status was based on self-reported smoking history and pack-years smoked.  In the stratified analyses never-smokers having smoked 0 pack-years and ever-smokers having smoked \> 5 pack-years were included, thereby excluding subjects with \> 0 and ≤ 5 pack-years.  Subjects were defined as having airflow obstruction based on having a pre-bronchodilator FEV~1~/FVC ratio (%) \< 70% or \< LLN (based on Global Lung Initiative 2012 (GLI-2012)) \[[@CR21]\].  All subjects completed pulmonary function testing according to ATS or ERS criteria \[[@CR22]\].  Additional details are provided in Additional file [1](#MOESM1){ref-type="media"}. 

Genotyping {#Sec4}
----------

The IlluminaCytoSNP-12 arrays were used to genotype blood samples in LifeLines and the Vlagtwedde-Vlaardingen study.  SNPs with a genotype call-rate ≥ 95%, minor allele frequency ≥ 1% and Hardy-Weinberg *p*-value ≥10^− 4^ were included.  Non-Caucasian samples and first-degree relatives were excluded based on self-reporting, outlier (Identity By State) and principal component analysis.  After quality control, 227,981 genotyped SNPs were included in the discovery analyses (LifeLines) and 242,926 genotyped SNPs were included in the replication analyses (Vlagtwedde-Vlaardingen).  Only genotyped SNPs were included in the analyses to prevent introducing bias, since it is known that imputation can reduce the effect size estimation, especially if healthy controls are used as reference \[[@CR23]\].  Blood samples in the Rotterdam study were genotyped with the 610 K and 660 K Illumina arrays and similar QC criteria as in the other cohorts were applied. 

Statistical analysis {#Sec5}
--------------------

Four separate GWASes were performed assessing the genetic associations between the two definitions of airflow obstruction, stratified by smoking status, for both the LifeLines (discovery) and Vlagtwedde-Vlaardingen (replication) studies.  Logistic regression (additive genetic model) was performed using PLINK (v1.07) \[[@CR24]\].  The "FEV~1~/FVC \< 70%" model was adjusted for sex, age and height.  The "FEV~1~/FVC \< LLN" model was not adjusted for these variables, since they are included in the LLN calculation.  In ever-smokers, the models were additionally adjusted for pack-years and current-smoking.  We used different *p*-value thresholds to assess the number of overlapping SNPs between the two definitions.  In addition, to quantify the agreement of the results between the two definitions and between never-and ever-smokers, we calculated the Pearson correlation coefficient between the *p*-values and between the ORs. 

Power simulations {#Sec6}
-----------------

Or study has a relative small sample size (*n* = 5070) and therefore relative low power.  We assessed the effect of low power on the overlap between the two definitions by increasing our never-smoking discovery sample (LifeLines) 2 (*n* = 10,140) and 4 (*n* = 20,280) times.  In addition, to assess if our results were spurious, we used our never-smoking discovery sample and randomly allocated 10 times the airflow obstruction cases but keeping the same distribution as in our original dataset (FEV~1~/FVC \< 70%: *n* = 548, FEV~1~/FVC \< LLN: *n* = 401, overlapping cases: *n* = 371 (64%)).  For both simulation studies, we repeated the GWAS analyses on both airflow obstruction definitions in the created datasets and compared the number of overlapping SNPs. 

Validation of overlapping SNPs {#Sec7}
------------------------------

Only the top overlapping SNPs between the two airflow obstruction definitions in the discovery sample (LifeLines) were evaluated in the SNP validation sample, the Vlagtwedde-Vlaardingen study and RS I to III.  A fixed-effects meta-analysis of the effect estimates weighted by the inverse of the standard errors from all four validation cohorts was performed using METAL (v2011) \[[@CR25]\].  We considered replication if the meta-analysis *p*-value was below the Bonferroni corrected p-value defined as 0.05/number of overlapping SNPs and, in addition, had the same direction of effect in all cohorts.  In addition, SNP\*ever-smoking interactions were estimated and we assessed if the overlapping SNPs were associated with gene expression levels in lung tissue within a 4 Mb window around the SNP (2 Mb on either side of the SNP), using data from the lung eQTL consortium \[[@CR26]\].  In total, 1087 subjects were included in the linear regression model, adjusted for disease status, age, sex, smoking, and cohort specific principal components.  SNPs with a p-value below the Bonferroni corrected threshold (*p* = 0.05/number of probesets) were considered significant eQTLs.  See Additional file [1](#MOESM1){ref-type="media"} and GEO accession numbers GSE23546 and GPL10379 for additional information. 

Results {#Sec8}
=======

Population characteristics {#Sec9}
--------------------------

The LifeLines cohort (discovery sample) included 5070 never-smokers and 4855 ever-smokers with complete data on all covariates (see Table [1](#Tab1){ref-type="table"}).  Of the never-smokers in LifeLines, 96% had a concordant airflow obstruction classification for the two definitions: 89% did not have airflow obstruction and 7% did have airflow obstruction.  The remaining 4% had a discordant classification (see Additional file [1](#MOESM1){ref-type="media"}: Table S1A).  Figure [1](#Fig1){ref-type="fig"}a shows that of all never-smoking subjects with airflow obstruction based on at least one airflow obstruction definition (*n* = 578), 36% had a discordant airflow obstruction classification.  Of the ever-smokers, 93% was classified concordantly: 77% did not have airflow obstruction and 17% did have airflow obstruction.  The remaining 7% had a discordant classification (see Additional file [1](#MOESM1){ref-type="media"}: Table S1B).  Of all ever-smoking subjects with airflow obstruction based on at least one definition (*n* = 1138), 30% had a discordant airflow obstruction classification (see Fig.  [1](#Fig1){ref-type="fig"}a).  Subjects with an FEV~1~/FVC \< 70% and \> LLN were aged between 41 and 85, and subjects with an FEV~1~/FVC \> 70% and \< LLN were aged between 22 and 43.  These and other characteristics of the airflow obstruction groups separately for never- and ever-smokers in LifeLines are shown in Additional file [1](#MOESM1){ref-type="media"}: Table S2.Table 1Characteristics of never- and ever-smokers included in the current studyNever-smokersEver-smokersLifeLinesVla-Vla^b^RS IRS IIRS IIILifeLinesVla-Vla^b^RS IRS IIRS IIIN with no missing data507043240837974748558236405831088Males, N (%)1942 (38)103 (23)81 (20)124 (33)308 (41)2312 (48)590 (72)375 (59)343 (59)531 (49)Age (yrs), median (min-max)46 (18--89)54 (36--79)78 (72--94)71 (65--98)62 (51--93)49 (22--85)53 (35--79)79 (72--95)71 (65--93)62 (52--93)Height (cm), mean (SD)174 (9)166 (9)163 (9)166 (9)171 (9)175 (9)173 (8)168 (9)171 (9)172 (9)Current-smokers, N (%)----------2171 (45)478 (58)99 (16)108 (19)298 (27)Pack-years (yrs), mean (SD)----------17 (11)27 (21)29 (21)29 (21)27 (20)Pulmonary function, mean (SD) FEV~1~%predicted (%), mean (SD)^a^98 (13)91 (13)101 (18)102 (16)101 (16)94 (14)85 (15)98 (23)96 (21)99 (19) FEV~1~/FVC (%), mean (SD)^b^78 (7)76 (7)77 (7)79 (6)78 (6)75 (8)72 (10)74 (9)74 (9)75 (8) FEV~1~/FVC \<70%, N (%)548 (11)77 (18)51 (13)27 (7)75 (10)1107 (23)287 (35)163 (26)151 (26)235 (22) FEV~1~/FVC \< LLN, N (%)401 (8)49 (11)9 (2)9 (2)25 (3)833 (17)212 (26)63 (10)62 (11)122 (11) Both FEV~1~/FVC \<70% and \< LLN, N (%)371 (7)49 (11)9 (2)9 (2)25 (3)802 (17)207 (25)63 (10)62 (11)122 (11)Moderate/severe COPD, N (%)^*c*^106 (2)26 (6)21 (5)9 (2)16 (2)310 (6)145 (18)90 (14)90 (15)104 (10)Discovery sample = LifeLines cohort study, replication samples = Vlagtwedde-Vlaardingen (Vla-Vla) study, and SNP validation sample = Vla-Vla and RS I to III*FEV*~*1*~ forced expiratory volume in one second, *FVC* forced vital capacity^a^FEV~1~%predicted is based on the reference equation by GLI-2012 \[[@CR21]\]^b^FEV~1~%IVC (inspiratory vital capacity) for the Vlagtwedde-Vlaardingen study (Vla-Vla)^c^COPD GOLD stage 2 and up (FEV~1~/FVC \<70% and FEV~1~%*p* \< 80% based on pre-bronchodilator measurements)Fig.  1Venn diagrams showing the overlap between the two definitions of airflow obstruction for the number of subjects classified as having airflow obstruction (**a**) and the number of identified SNPs with *p* \< 10^− 4^ (**b**) in LifeLines (discovery sample)

The Vlagtwedde-Vlaardingen study (replication sample) included 432 never-smokers and 823 ever-smokers (see Table [1](#Tab1){ref-type="table"}).  Of the Vlagtwedde-Vlaardingen study, 94% and 90% of the never- and ever-smokers were classified concordantly based on the two definitions (see Additional file [1](#MOESM1){ref-type="media"}: Table S1 C-D). 

The SNP validation sample used for the SNP validation meta-analysis included 1966 never-smokers and 3134 ever-smokers from the Vlagtwedde-Vlaardingen study and RS I to III (see Table [1](#Tab1){ref-type="table"}). 

GWAS results {#Sec10}
------------

There was minimal population stratification in all analyses of LifeLines, indicated by the genomic inflation factor lambda (λ: 1.0002--1.0217, see Additional file [1](#MOESM1){ref-type="media"}: Figure S1).  The results based on a *p* \< 10^− 4^ of all four analyses in LifeLines are given in See (Additional file [1](#MOESM1){ref-type="media"}: Tables S3-S6), including the Manhattan plots (see Additional file [1](#MOESM1){ref-type="media"}: Figures S2 and S3).  For comparison, the effect estimates of both airflow obstruction definition are given in these tables.  Summary statistics (*p* values, betas, and standard errors for all SNPs that were tested) of the GWAS result of both the Lifelines Cohort Study and the Vlagtwedde-Vlaardingen study are provided in Additional file [2](#MOESM2){ref-type="media"}. 

Overlap between the results {#Sec11}
---------------------------

We used several *p*-value thresholds to assess the overlap between the GWAS results of both airflow obstruction definitions separately in never- and ever-smokers of LifeLines (see Table [2](#Tab2){ref-type="table"}).  A threshold of 0.05 resulted in the observation that 26% and 29% of the SNPs were overlapping between the two airflow obstruction definitions in never- and ever-smokers, respectively.  Three percent of the SNPs were overlapping between never- and ever-smokers for both definitions.  A smaller *p*-value threshold resulted in a lower percentage of overlap e.g. a threshold of p \< 10^− 4^ resulted in 4% and 6% overlapping SNPs between the two airflow obstruction definitions in never- and ever-smokers, respectively (see Fig. [ 1](#Fig1){ref-type="fig"}b), and zero overlap between never- and ever-smokers using the same definition of airflow obstruction.  Similar observations were made in the replication sample the Vlagtwedde-Vlaardingen study (see Table [2](#Tab2){ref-type="table"}), since at *p* \< 0.05 the overlap between the definitions was 24% and 25% in never- and ever-smokers, respectively, and 2% of the SNPs were overlapping between never- and ever-smokers for both definitions. Table 2Table showing the number of SNPs with a *p*-value below the mentioned threshold for both FEV~1~/FVC \< 70% and \< LLN analysis and the overlapThresholdNever-smokersEver-smokers (\> 5 py)Overlap never- and ever-smokers70%LLNOverlap70%LLNOverlap70%LLN*Discovery analysis (LifeLines)* \< 0.0511,37711,4754755 (26%)12,11412,3665445 (29%)625 (2.7%)616 (2.7%) \<0.0122322197673 (18%)25222493824 (20%)23 (0.5%)21 (0.4%) \< 10^−3^22223354 (13%)24626676 (17%)1 (0.2%)0 \< 10^− 4^31272 (4%)21293 (6%)00 \< 10^−5^440321 (25%)00 \< 10^− 6^10000000 \<Bonferroni\*00000000*Replication analysis (Vlagtwedde-Vlaardingen study)* \<0.0510,70210,2954026 (24%)12,59212,5714976 (25%)487 (2.1%)488 (2.2%) \<0.0118571807544 (17%)25672609754 (17%)17 (0.4%)12 (0.3%) \< 10^−3^16417438 (13%)22326241 (9%)00 \< 10^−4^18196 (19%)1221000 \< 10^−5^732 (25%)00000 \< 10^−6^10000000 \<Bonferroni\*00000000\**P* \< 2.19 × 10^−7^

The correlations between the SNP-specific *p*-values and ORs from the two airflow obstruction definitions were 0.48 (*p*-value) and 0.78 (OR) in never-smokers, and 0.51 (*p*-value) and 0.81 (OR) in ever-smokers (see Fig.  [2](#Fig2){ref-type="fig"}).  Between never- and ever-smokers the correlations of the SNP-specific p-values were 0.0008 for FEV~1~/FVC \< 70% and 0.002 for FEV~1~/FVC \< LLN, and for the OR the correlation was − 0.02 for both definitions.  Similar observations were made in the replication sample, the Vlagtwedde-Vlaardingen study.  The correlations between the two definitions were 0.45 (*p*-value) and 0.76 (OR) in never-smokers, and 0.41 (*p*-value) and 0.74 (OR) in ever-smokers.  Between never- and ever-smokers the correlations were − 0.001 (*p*-value) and 0.015 (OR) for FEV~1~/FVC \< 70% and − 0.003 (*p*-value) and 0.004 (OR) for FEV~1~/FVC \< LLN.Fig.  2Pearson correlation between the p-values (**a**/**c**) or OR (B/D) of FEV~1~/FVC \< 70% and \< LLN analyses separately for never- (**a**/**b**) and ever-smokers (**c**/**d**) in LifeLines (discovery sample)

Power simulations {#Sec12}
-----------------

We found that the percentage of overlap increased when we expanded our never-smoking identification sample 2 and 4 times (see Additional file [1](#MOESM1){ref-type="media"}: Table S7).  The overlap between SNPs with *p* \< 10^− 4^ was 3.6% in the original dataset, 16.5% in the 2x dataset and 26.6% in the 4x dataset.  In addition, when we randomly allocated cases 10 times in our identification sample, we found the percentage of overlap between the two definitions at *p* \< 10^− 4^ varied between 0 to 16%, compared to 4% in the original dataset (see Additional file [1](#MOESM1){ref-type="media"}: Table S8). 

Validation of overlapping SNPs {#Sec13}
------------------------------

In never-smokers of LifeLines, two SNPs were overlapping between the FEV~1~/FVC \< 70% and \< LLN definitions at a threshold of p \< 10^− 4^ (see Table [3](#Tab3){ref-type="table"}).  The first SNP (rs7519348) is located in an intron of the gene nuclear transcription factor Y subunit C (*NFYC*)*,* and the second SNP (rs6913003) is located in an intron of fatty acid binding protein 7 (*FABP7*, see Additional file [1](#MOESM1){ref-type="media"}: Figure S4 and S5 for LocusZoom plots).  The minor alleles of both SNPs were associated with a higher risk of airflow obstruction and had comparable odds ratios in both analyses.  The SNP in *NFYC* (rs7519348) was significantly associated with FEV~1~/FVC \< LLN in the SNP validation meta-analysis (*p* = 0.034), but did not pass the multiple testing correction (0.05/2 = 0.025), and was not significantly associated with FEV~1~/FVC \< 70% in the SNP validation meta-analysis (*p* = 0.07).  The SNP in *FABP7* (rs6913003) was not significantly associated with FEV~1~/FVC \< 70% or \< LLN in the SNP validation meta-analyses (*p* = 0.08 in both), although the direction of effect was the same in all independent cohorts.  Both SNPs did not reach genome-wide significance according to the Bonferroni-corrected threshold (*p* \< 2.19 × 10^− 7^) in the discovery analysis (LifeLines) or meta-analysis of both the discovery and SNP validation samples (see Table [3](#Tab3){ref-type="table"} and Additional file [1](#MOESM1){ref-type="media"}: Table S9).  Yet, the odds ratios were comparable between all analyses (see Additional file 1: Table S10 and Figure S6).  These two overlapping SNPs were not associated with airflow obstruction in ever-smokers and these associations were significantly different between ever- and never-smokers as shown in the interaction analysis (see Additional file [1](#MOESM1){ref-type="media"}: Table S11).Table 3Results of the overlapping SNPs identified in both genome-wide association studies on FEV~1~/FVC \< 70% and FEV~1~/FVC \< LLN in never- and ever-smokersDiscovery analysisSNP validation meta-analysisDirection of effect in the independent cohorts^a^SNPChrA1MAFGeneTestORSEPORSEPI^2^*Never-smokers(n = 5070)(n = 1966)* rs75193481A33%*NFYC* (intronic)\<70%1.360.074.92 × 10^−6^1.210.110.070.0++++ 0\<LLN1.370.072.27 × 10^−5^1.400.160.030.0+++++ rs69130036T4%*FABP7* (intronic)\< 70%1.900.138.99 × 10^−7^1.480.220.080.0+++++\<LLN1.830.141.83 × 10^−5^1.720.310.080.0+++++*Ever-smokers(n = 4855)(n = 3134)* rs131180834A45%*HHIP* (342 kb 5′)\<70%1.230.054.36 × 10^−5^0.970.070.6017.0+ 0 − − +\<LLN1.260.052.27 × 10^−5^1.050.080.5746.7++ 0 − + rs707421010G18%*ST8SIA6* (62 kb 5′)\<70%1.350.063.08 × 10^−6^1.050.080.5764.7+ − +++\<LLN1.330.073.23 × 10^−5^0.990.100.9564.5+ − + − + rs493039011G24%*C11orf80* (intronic)\<70%0.760.069.40 × 10^−6^1.020.070.740.0− + 0 -- 0\<LLN0.730.076.68 × 10^−6^1.010.090.9110.1−++ − −SNPs were selected based on having a p-value \< 10^−4^ in both the discovery analyses on the fixed ratio of 70% and LLN.  The logistic regression model of FEV~1~/FVC \< 70% was adjusted for sex, age and height, the LLN model was not adjusted.  Ever-smoking models were additionally adjusted for pack-years and current-smoking.  Discovery sample = LifeLines cohort study, and SNP validation sample = Vlagtwedde-Vlaardingen and RS I to III.  A1 = minor allele (effect allele), *MAF* = minor allele frequency, OR = Odds Ratio, SE = standard error and P = *p*-value, I^2^ = heterogeneity measure^a^ Order: LifeLines, Vlagtwedde-Vlaardingen, and Rotterdam Study I to III.  + represents an OR \> 1, − represents an OR \< 1, and 0 represents is an OR between 0.95 and 1.05 (no effect)

In ever-smokers of LifeLines, three SNPs were overlapping between the two analyses in at *p* \< 10^− 4^ (see Table [3](#Tab3){ref-type="table"}).  The first SNP (rs13118083) is annotated to hedgehog interacting protein (*HHIP*, 342 kb away), but is located within the long non-coding RNA *LOC105377462* according to the SNP database by NCBI (<https://www.ncbi.nlm.nih.gov/SNP/>).  The second SNP (rs7074210) is located approximately 62 kb from ST8 Alpha-N-Acetyl-Neuraminide Alpha-2,8-Sialyltransferase 6 (*ST8SIA6*), and the last SNP (rs4930390) is annotated to Chromosome 11 Open Reading Frame 80 (*C11orf80*).  The minor alleles of the first 2 SNPs were associated with a higher risk of airflow obstruction and the minor allele of rs4930390 with a lower risk.  The effect was significantly different between never- and ever-smokers for SNP rs4930390 according to both definitions and for rs7074210 in the FEV~1~/FVC \< 70% analyses (see Additional file [1](#MOESM1){ref-type="media"}: Table S11).  The three SNPs were not replicated in the SNP validation sample (see Table [3](#Tab3){ref-type="table"} and Additional file [1](#MOESM1){ref-type="media"}: Tables S9-S10). 

Gene expression in lung tissue {#Sec14}
------------------------------

The minor allele (G) of rs7519348 (overlapping SNP in never-smokers) was associated with higher gene expression of *NFYC* Antisense RNA 1 (*NFYC-AS1*) in lung tissue (Fig.  [3](#Fig3){ref-type="fig"}).  Summary statistics of the eQTL analysis for all overlapping SNPs at *p* \< 10^− 4^ are provided in Additional file [2](#MOESM2){ref-type="media"}.Fig.  3Results of eQTL analysis in lung tissue for rs7519348, an overlapping SNP in never-smokers.  The unadjusted mean log2 microarray intensity and 95% CI are plotted, obtained from a meta-analysis of three cohorts included in the lung eQTL dataset

Discussion {#Sec15}
==========

We investigated the genetic overlap between GWASes using two airflow obstruction definitions in the same population (FEV~1~/FVC \< 70 or \< LLN).  We expected a reasonable overlap in associated SNPs between the two definitions, since 96% of the never-smokers and 93% of the ever-smokers were classified the same way in the discovery sample LifeLines.  Surprisingly, only a very small proportion (4% and 6%) of SNPs was overlapping at *p* \< 10^− 4^ (see Fig. [ 1](#Fig1){ref-type="fig"}).  Even with different significance thresholds the overlap was limited (26% and 29% at *p* \< 0.05) (see Table [2](#Tab2){ref-type="table"}).  The same observation was made in the replication sample, the Vlagtwedde-Vlaardingen study.  In this cohort, 94% and 90% of the never- and ever-smokers, respectively, were classified concordantly, but at p \< 0.05 only 24% or 25% of the SNPs were overlapping.  In addition, the effect estimates for the two airflow obstruction definitions correlated strongly in both cohorts but the *p*-values showed more variation and correlated only moderately resulting in different top-hits depending on the obstruction definition (see Fig. [ 2](#Fig2){ref-type="fig"}).  Thus, the chosen strategy and definition of airflow obstruction had a substantial influence on the GWAS results.  This implies that in a discovery-replication design with a predetermined selection p-value, different genetic variants would be followed-up depending on the definition used.  In addition, there was no correlation between the p-values nor between the ORs of never- and ever-smokers in both cohorts.  None of the selected SNPs overlapped between never- and ever-smokers at p \< 10^− 4^, and at *p* \< 0.05 the overlap was only 3% in LifeLines (discovery sample) and 2% in Vlagtwedde-Vlaardingen (replication sample, see Table [2](#Tab2){ref-type="table"}).  The current study therefore also highlights the importance of stratifying the analysis according to smoking status. 

The difference between results from the two definitions might be explained by the fact that obstructive airway diseases are heterogeneous diseases with multiple phenotypes, symptoms and comorbidities.  It might thus be beneficial for future GWA studies to focus more on specific COPD subtypes rather than on a broad definition of airflow obstruction or COPD that can be caused by multiple underlying physiologic and genetic mechanisms.  In previous GWA studies, in mainly smokers, on classical COPD phenotypes like emphysema and chronic bronchitis, the well-known general COPD genes (*HHIP*, *CHRNA* and *FAM13A*) were consistently identified \[[@CR27]--[@CR32]\].  Perhaps, to identify specific genetic pathways underlying specific COPD phenotypes we should not study the classical COPD phenotypes, but rather clinical COPD subtypes based on symptoms, comorbidities or pathology. 

The *CHRNA5/3* and *HHIP* regions were overlapping between six previous GWA studies on airflow obstruction, using different airflow obstruction definitions and populations \[[@CR11]--[@CR16]\].  In the current study, two of the identified SNPs in ever-smokers were located in the *CHRNA5* and *HHIP* regions as well, pointing towards a robust genetic association of these regions with airflow obstruction and COPD (see Additional file [1](#MOESM1){ref-type="media"}: Table S6).  Likewise, most of previously identified regions associated with airflow obstruction or COPD were nominal significant (*p* \< 0.05) in the current study (see Additional file [1](#MOESM1){ref-type="media"}: Table S12).  Out of the 22 loci identified by the study of Hobbs et al, SNPs in 18 loci were associated with at least one of the airflow definitions at a nominal significance (10 SNPs in never-smokers and 12 SNPs in ever smokers) \[[@CR14]\].  In never-smokers, 6 of the 10 SNPs were significantly associated with both definitions and in ever-smokers 7 of the 12 SNPs were significantly associated with both definitions.  Some SNPs were significant in both never- and ever-smokers (e.g. *HHIP, PID1* and *THSD4*), while others were either only significant in never-smokers (e.g. *FAM13A, DSP* and *RIN3*) or in ever-smokers (e.g. *CHRNA5, TET2* and *ADGRG6*).  In addition, many of the loci previously associated with lung function outcomes (FEV~1~, FVC, and FEV~1~/FVC) were also nominal significant (p \< 0.05) in the current study (see Additional file [1](#MOESM1){ref-type="media"}: Table S13).  Specifically, of the loci reported by Wain et al.,  23 out of 28 loci for FEV~1~, 10 out of 17 loci for FVC and 38 out of 51 loci for FEV~1~/FVC were associated with at least one of the airflow definitions at a nominal significance \[[@CR33]\].  Lastly, we also checked if the top overlapping SNPs were associated with lung function outcomes in our previous GWA studies on FEV~1~, FEV~1~/FVC and FEF~25--75~ \[[@CR34], [@CR35]\].  A SNP annotated to *HHIP* was associated with FEV~1~/FVC and FEF~25--75~ in both never- and ever-smokers (results were replicated) and the *CHRNA5/3* region was only associated with FEV~1~/FVC in ever-smokers.  The *NFYC* and *FABP7* regions were associated with FEV~1~/FVC (*p* = 4.40 × 10^− 4^ and *p* = 1.87 × 10^− 4^) in never-smokers, and the *FABP7* SNP was also associated with FEF~25--75~ levels (*p* = 0.026).  Interestingly, the *NFYC* region was also overlapping between the current study and the study by Pillai et al.  We identified multiple SNPs annotated to *NFYC,* whereas Pillai et al.  identified a SNP (rs3767943) in the gene *KCNQ4*, which is located on the right side (3′) of *NFYC* \[[@CR15]\].  The *NFYC* region might therefore be an interesting region to further study the underlying mechanisms of its association with airflow obstruction. 

A SNP in the intron of *NFYC* and a SNP in *FABP7* were the two overlapping SNPs between the airflow obstruction definitions at *p* \< 10^− 4^ in never-smokers and showed the same direction of effect in the five independent cohorts.  The minor allele of the SNP in *NFYC* (rs7519348) was associated with a higher risk of airflow obstruction.  This gene is a highly conserved transcription factor that is predicted by GeneGlobe to bind promoter regions of 218 genes (see Additional file [1](#MOESM1){ref-type="media"}: Table S14) including genes previously associated with lung related outcomes, like *ADORA2B, AKAP9, CD163, ELMOD2, HLA-DPB1, ITPR2, KLF10* and *SERPINA6* \[[@CR27], [@CR36]--[@CR42]\].  In more detail, *HLA-DPB1* is a known COPD gene related to disease severity, *SERPINA6* was associated with emphysema, a deletion in *ADORA2B* was shown to be associated with a decrease in lung fibrosis and pulmonary hypertension, and *ELMOD2* is a candidate gene for familial idiopathic pulmonary fibrosis \[[@CR27], [@CR36], [@CR39], [@CR40]\].  The identified SNP was not associated with expression levels of *NFYC* in lung tissue, but was an eQTL for a probeset annotated to *NFYC-AS1.*  The function of this specific antisense-RNA, which are generally thought to have a regulatory role, is still unknown. 

The minor allele of the SNP in *FABP7* (rs6913003) was also associated with a higher risk of airflow obstruction in never-smokers.  This SNP was not associated with the expression of *FABP7* or other genes in lung tissue. * FABP7* is an intracellular lipid-binding protein, involved in long-chain fatty acids transport and cell proliferation \[[@CR43]\].  It may be involved in abnormal pulmonary development, since lower expression of *FABP7* was found in patients with congenital cystic adenomatoid malformation \[[@CR44]\].  In addition, higher expression of *FABP7* was seen in clear cell renal cell carcinoma and the authors suggested that the gene activates the *ERK* and *STAT3* signalling pathways \[[@CR45]\]. * STAT3* was implicated to play a role in pulmonary inflammation and thus *FABP7* might indirectly be involved in airflow obstruction \[[@CR46]\]. 

We were aware of the risk for spurious findings due to the low power of our study and thus we validated our top overlapping SNPs in 4 independent validation cohorts.  We furthermore investigated the effect of low power on the overlap between the two definitions by increasing our dataset 2 and 4 times.  We found that the percentage of overlap increases when the sample size increases, but still the number of SNPs that do not overlap remains high, i.e. 73.4% when the sample size increased 4-fold.  So even when the study power is greatly increased, different SNPs will be found depending on the airflow obstruction definition tested.  We also performed a simulation study by 10 times randomly allocating airflow obstruction cases and based on this simulation, we have to conclude that the differences and overlap we found could be chance findings, but that is why we validated the overlapping SNPs in 4 independent validation cohorts. 

We only assessed a modest number of SNPs (*n* = 227,981 SNPs) compared to previous large GWAS studies (*n* \> 1 million SNPs), since we only included genotyped SNPs to prevent any bias by imputation.  The disadvantage of this approach is that we may have a lower genomic coverage.  Another limitation of the current study is the use of pre-bronchodilator measurements to define airflow obstruction, which preferably should be based on post-bronchodilator measurements.  Especially subjects with asthma could be misclassified as having airflow obstruction, but the results of the overlapping SNPs did not change in a sensitivity analysis excluding asthmatics or adjusting for asthma (see Additional file [1](#MOESM1){ref-type="media"}: Table S15).  Moreover, only a low number of never-smoking subjects had an FEV~1~/FVC \< LLN in the three Rotterdam Study cohorts, but nevertheless results were replicated in these never-smokers.  Finally, the "FEV~1~/FVC \< 70%" model was adjusted for sex, age and height, but the "FEV~1~/FVC \< LLN" model was not adjusted for these variables, since they are included in the LLN calculation.  If we do however adjust the "FEV~1~/FVC \< LLN" model for these variables, the results do not change.  The top SNPs are the same and the correlation between *p*-values for the LLN models adjusted and not adjusted is 0.98.  In addition, the reported correlation in never-smokers between the two definitions was 0.48 for p-values and 0.78 for OR.  If we use the LLN adjusted model the correlation is 0.48 and 0.79, respectively.  This confirms that we used appropriate models to assess the genetic overlap between the two airflow definitions. 

Conclusions {#Sec16}
===========

The definition of airflow obstruction and the population under study may be important determinants of which SNPs are associated with airflow obstruction, and thus on which variants are selected for replication.  It is therefore important to use the same definition of airflow obstruction in future studies, especially in consortia.  In addition, future studies should focus more on specific COPD subtypes and subgroups (e.g. based on smoking status), since there was no overlap in results between never- and ever-smokers, pointing towards possible different underlying mechanisms.  Finally, our results suggest that the genes *FABP7* and *NFYC*(-*AS1*) could play a role in the pathogenesis of airflow obstruction in never-smokers. 

Additional files
================

 {#Sec17}

Additional file 1:Supplementary methods, tables and figures. ( DOCX 1790 kb) Additional file 2:GWAS summary statistics. ( XLSX 139809 kb)

COPD

:   Chronic Obstructive Pulmonary Disease

eQTL

:   Expression quantitative trait locus

*FABP7*

:   Fatty acid binding protein 7

FEV~1~/FVC

:   Ratio of forced expiratory volume in one second to forced vital capacity

GWAS

:   Genome-wide association study

MAF

:   Minor allele frequency

*NFYC*

:   Nuclear transcription factor Y subunit C

SNP

:   Single-nucleotide polymorphism

We thank Rob Bieringa, Joost Keers, René Oostergo, Rosalie Visser, and Jan Schouten for their work related to data collection and validation in the LifeLines cohort study and the Vlagtwedde-Vlaardingen study.  The authors would like to thank the staff at the Respiratory Health Network Tissue Bank of the FRQS for their valuable assistance with the lung eQTL dataset at Laval University.  We are grateful to the study participants and all staff involved in the LifeLines cohort study, Vlagtwedde-Vlaardingen Study, Rotterdam Study and lung eQTL database.  We would like to thank Anis Abuseiris, Karol Estrada, Dr. Tobias A. Knoch, and Rob de Graaf as well as their institutions Biophysical Genomics, Erasmus MC Rotterdam, The Netherlands, and especially the national German MediGRID and Services\@MediGRID part of the German D-Grid, both funded by the German Bundesministerium fuer Forschung und Technology under grants \#01 AK 803 A-H and \# 01 IG 07015 G for access to their grid resources. 

Funding {#FPar1}
=======

This study is sponsored by grant number 4.1.13.007 of Lung Foundation (Longfonds), the Netherlands.  DAvdP, KdJ and NA are supported by grant number 4.1.13.007 of Longfonds.  LL is a Postdoctoral Fellow of the Research Foundation - Flanders (FWO).  The LifeLines Biobank initiative has been made possible by funds from FES (Fonds Economische Structuurversterking), SNN (Samenwerkingsverband Noord Nederland) and REP (Ruimtelijke Economisch Programma).  The Vlagtwedde-Vlaardingen cohort study was supported by the Ministry of Health and Environmental Hygiene of the Netherlands and the Netherlands Asthma Fund (grant 187) and the genetic data of the cohort were funded by the Netherlands Asthma Fund (grant no.  3.2.02.51), the Stichting Astma Bestrijding, BBMRI-NL (Complementiation project) and the European Respiratory Society COPD research award 2011 to H.M. Boezen.  The Rotterdam Study was supported by the Erasmus MC and Erasmus University Rotterdam; the Netherlands Organisation for Scientific Research (NWO); the Netherlands Organisation for Health Research and Development (ZonMw); the Research Institute for Diseases in the Elderly (RIDE); the Netherlands Genomics Initiative (NGI); the Ministry of Education, Culture and Science; the Ministry of Health, Welfare and Sports; the European Commission (DG XII); and the Municipality of Rotterdam.  Genotyping and gene expression for the lung eQTL study was funded by Merck & Co. Inc. The lung eQTL study at Laval University was supported by the Fondation de l'Institut universitaire de cardiologie et de pneumologie de Québec, the Respiratory Health Network of the FRQS, the Canadian Institutes of Health Research (MOP - 123369).  Y.B. holds a Canada Research Chair in Genomics of Heart and Lung Diseases.  The sponsors of this study played no role in the design of the study, data collection, analysis, interpretation or in the writing and submission of the manuscript. 

Availability of data and materials {#FPar2}
==================================

Summary statistics of all analyses are available in Additional file [2](#MOESM2){ref-type="media"}.  The eQTL lung tissue datasets analysed during the current study are available at the Gene Expression Omnibus (GEO) repository, GEO accession numbers GSE23546 and GPL10379.  LifeLines data is available (at costs) to all scientists.  Scientists can apply for access to Lifelines data and samples by submitting a research proposal to the LifeLines biobank ([www.lifelines.net](http://www.lifelines.net)). 

DAvdP participated in the study design, analysis and interpretation of the data, and drafting of the manuscript, tables and figures.  HMB, DSP, CCvD and CMVD obtained funding.  KdJ, HMB, DSP, JMV, CCvD, CMVD, IN and NA determined the study design, participated in the analysis and interpretation of data, and critically supervised writing of the manuscript.  LL and GGB participated in collecting and analysing the Rotterdam Study data.  AF, YB, CAB, DSP, KH, and PDP participated in setting up the lung tissue database and analyses.  All authors approved the final version of the manuscript. 

Ethics approval and consent to participate {#FPar3}
==========================================

Written informed consent was provided by all subjects of all included cohorts.  The Lifelines Cohort Study and the Vlagtwedde-Vlaardingen study were approved by the Medical Ethics Committee of the University Medical Center Groningen, Groningen, the Netherlands.  The Rotterdam Study was approved by the medical ethics committee of Erasmus University.  The eQTL lung tissue database was approved by the ethics committees of the Institut universitaire de cardiologie et de pneumologie de Québec (Laval) and the UBC-Providence Health Care Research Institute Ethics Board (UBC).  The study protocol was consistent with the Research Code of the University Medical Center Groningen and Dutch national ethical and professional guidelines. 

Consent for publication {#FPar4}
=======================

Not applicable. 

Competing interests {#FPar5}
===================

LL reports personal fees from Boehringer Ingelheim GmbH, non-financial support from Novartis, grants from AstraZeneca, grants and non-financial support from European Respiratory Society and Belgian Respiratory Society, outside the submitted work.  The University of Groningen has received money for Professor Postma (DSP) regarding a grant for research from Astra Zeneca, Chiesi, Genentec, GSK and Roche.  Fees for consultancies were given to the University of Groningen by Astra Zeneca, Boehringer Ingelheim, Chiesi, GSK, Takeda and TEVA.  All other authors declare no competing interests. 

Publisher's Note {#FPar6}
================

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.