Genome-wide association and large-scale follow up identifies 16 new loci influencing lung function.

Pulmonary function measures reflect respiratory health and are used in the diagnosis of chronic obstructive pulmonary disease. We tested genome-wide association with forced expiratory volume in 1 second and the ratio of forced expiratory volume in 1 second to forced vital capacity in 48,201 individuals of European ancestry with follow up of the top associations in up to an additional 46,411 individuals. We identified new regions showing association (combined P < 5 × 10(-8)) with pulmonary function in or near MFAP2, TGFB2, HDAC4, RARB, MECOM (also known as EVI1), SPATA9, ARMC2, NCR3, ZKSCAN3, CDC123, C10orf11, LRP1, CCDC38, MMP15, CFDP1 and KCNE2. Identification of these 16 new loci may provide insight into the molecular mechanisms regulating pulmonary function and into molecular targets for future therapy to alleviate reduced lung function.

Introduction

Pulmonary function, reliably measurable by spirometry, is a heritable trait reflecting the physiological state of the airways and lungs[1]. Pulmonary function measures are important predictors of population morbidity and mortality[2-4], and are used in the diagnosis of chronic obstructive pulmonary disease (COPD), which ranks among the leading causes of death in developed and developing countries[5,6]. A reduced ratio of forced expiratory volume in 1 second (FEV1) to forced vital capacity (FVC) is used to define airway obstruction, and a reduced FEV1 is used to grade the severity of airway obstruction[7].

Recently, two large genome-wide association studies (GWAS), each comprising discovery sets of more than 20,000 individuals of European ancestry identified novel loci for lung function[8,9]. Recognizing the need for larger datasets to increase the power to detect loci of individually modest effect size, we conducted a meta-analysis of 23 lung function GWAS comprising a total of 48,201 individuals of European ancestry (Stage 1) and followed-up potentially novel loci in 17 further studies comprising up to 46,411 individuals (Stage 2). We identified 16 additional novel loci for lung function, and provided evidence corroborating association of loci previously associated with lung function[8-11]. Our findings implicate a number of different mechanisms underlying regulation of lung function and highlight loci shared with complex traits and diseases, including height, lung cancer, and myocardial infarction.

Results

Genome-wide analysis (stage 1)

Meta-analyses for cross-sectional lung function measures were undertaken for approximately 2.5 million genotyped or imputed SNPs across 23 studies with a combined sample size of 48,201 adult individuals of European ancestry. Characteristics of the cohort participants and the genotyping are shown in Supplementary Tables 1A and 1B. FEV1 and FEV1/FVC were adjusted for ancestry principal components, age, age2, sex, and height as covariates. Association testing of the inverse-normal transformed residuals for FEV1 and FEV1/FVC assumed an additive genetic model and was stratified by ever-smoking (versus never-smoking) status. Meta-analyses of the smoking strata within study, and the study-specific results, were undertaken using inverse variance weighting (the inverse of the standard error squared was used as the weight). We applied genomic control twice at study level (to each smoking stratum separately and to the study level pooled estimates) and also at meta-analysis level to avoid inflation of test statistics due to cryptic population structure or relatedness (see Supplementary Table 1A for study level estimates). Our application of genomic control at the three stages is likely to be overly conservative because it has recently been shown that in large meta-analyses, test statistics are expected to be elevated under polygenic inheritance even when there is no population structure[12]. Test statistic inflation (λGC) prior to applying genomic control at meta-analysis level was 1.12 for FEV1 and 1.09 for FEV1/FVC. Genomic inflation estimates increase with sample size, as has been shown for other traits[13-15]; standardised estimates to a sample of 1000 individuals (λGC_1000) were 1.002 for FEV1 and 1.002 for FEV1/FVC. Plots of meta-analysis P-values for FEV1 and FEV1/FVC against a uniform distribution of P values expected under the null hypothesis showed deviations which were attenuated, but persisted, after removal of SNPs in loci reported previously, consistent with additional loci being associated with lung function (Supplementary Figure 1A).

Follow-up analysis (stage 2)

Twenty-nine new loci showing evidence of association with lung function (P<3×10−6) in Stage 1 were followed up in Stage 2 by utilizing in silico data from seven studies, and by undertaking additional genotyping in 10 studies for the 10 highest-ranked SNPs (Figure 1). Full details of the SNP selection are given in the Online Methods. Inverse variance weighting meta-analysis was performed across Stages 1 and 2, and two-sided p-values were obtained for the pooled estimates. Sixteen new loci reached genome-wide significance (P<5×10−8) and showed consistent direction of effects in both stages, comprising 12 new loci for FEV1/FVC, 3 new loci for FEV1, and one new locus reaching genome-wide significance for both traits (Figure 2, Table 1). To assess the heterogeneity across studies included in Stage 1 and Stage 2, Chi-square tests were undertaken for all 16 SNPs; and none of them was statistically significant after applying a Bonferroni correction for 16 tests. The sentinel SNPs at these loci were in or near the genes MFAP2 (1p36.13), TGFB2/LYPLAL1 (1q41), HDAC4/FLJ43879 (2q37.3), RARB (3p24.2), MECOM (EVI1) (3q26.2), SPATA9/RHOBTB3 (5q15), ARMC2 (6q21), NCR3/AIF1 (6p21.33), ZKSCAN3 (6p22.1), CDC123 (10p13), C10orf11 (10q22.3), LRP1 (12q13.3), CCDC38 (12q22), MMP15 (16q13), CFDP1 (16q23.1) and KCNE2/C21orf82 (21q22.11) (Supplementary Figures 1B and 1C). The strongest signals in AGER (rs2070600)[8,9] and two of the novel signals (rs6903823 in ZKSCAN3 and rs2857595, upstream of NCR3) lie within a~3.8Mb interval at 6p21.32-22.1 which is characterised by long-range linkage disequilibrium. Nevertheless, the leading SNPs in these regions which are within the major histocompatibility complex (MHC) were statistically independent (Supplementary Note).

Figure 1

Study design. A total of 34 SNPs showing novel evidence of association (P<3×10−6) with FEV1 and/or FEV1/FVC in a meta-analysis of the Stage 1 studies were followed up in Stage 2. Studies with a combined total of 24,737 individuals undertook genotyping and association testing of the top 10 SNPs. Seven studies (*) with a combined total of 11,275 individuals had genome-wide association data and provided results for up to 34 SNPs. #GS:SFHS undertook genotyping on a 32-SNP multiplex genotyping platform and so included the 32 top ranking SNPs (including proxies and both SNPs from regions which showed association with both FEV1 and FEV1/FVC). This assay failed for one SNP (rs3769124) which was subsequently replaced with the 33rd SNP (rs4762767). SNP rs2284746 was excluded due to poor clustering. Although rs3743563 was chosen as proxy for rs12447804 which had N effective <80% in the Stage 1 meta-analysis, BHS2 were unable to genotype rs3743563 and so undertook genotyping for rs12447804 instead. See Table 1 for definitions of abbreviations.

Figure 2

Manhattan plots of association results for (a) FEV1/FVC and (b) FEV1. Manhattan plots ordered by chromosome position. SNPs for which −log10P>5 are indicated in red. Novel regions which reached genome-wide significance after Stage 1 + Stage 2 are labelled.

Table 1

Loci associated with lung function. Shown are FEV1 and FEV1/FVC results for the leading SNPs, ordered by chromosome and position, for each independent locus associated (P <5×10−8) with FEV1 or FEV1/FVC in a joint analysis of up to 94,612 individuals of European ancestry from the SpiroMeta-CHARGE GWAS (Stage 1) and follow-up (Stage 2). Two-sided p-values are given for Stage 1, Stage 2 and the joint meta-analysis of all stages. P-values reaching genome-wide significance (P<5×10−8) in the joint meta-analysis of all stages are indicated in bold. SNPs reaching independent replication in Stage 2 (P=0.05/34 = 1.47×10−3) are indicated with their Stage 2 p-value in bold. The sample sizes (N) shown are the effective sample sizes. Effective sample size within each study is the product of sample size and imputation quality metric. Joint meta-analysis includes data from Stage 1 and Stage 2. Beta values reflect effect-size estimates on an inverse-normal transformed scale after adjustments for age, age2, sex, height and ancestry principal components. The estimated proportion of the variance explained by each SNP can be found in Supplementary Table 6.

Chr.	SNP_ID(NCBI36 position), function	Codedallelefreq.	Measure	Stage 1				Stage 2				Joint meta-analysis ofall stages
				Beta (s.e.m.)	P	Codedallele freq.	N	Beta (s.e.m.)	P	Codedallele freq.	N	Beta (s.e.m.)	P
1	rs2284746(17179262), MFAP2(intron)	G	FEV1/FVC	−0.042 (0.007)	2.47E-09	0.516	45944	−0.038 (0.007)	2.64E-07	0.522	35371	−0.04 (0.005)	7.50E-16
			FEV1	0.008 (0.007)	2.78E-01			0.006 (0.007)	3.70E-01			0.007 (0.005)	1.48E-01
1	rs993925(216926691), TGFB2(downstream)	T	FEV1/FVC	0.04 (0.007)	2.54E-07	0.308	42402	0.023 (0.01)	1.76E-02	0.348	21414	0.034 (0.006)	1.16E-08
			FEV1	0.025 (0.007)	1.51E-03			0.003 (0.007)	7.29E-01			0.014 (0.005)	8.71E-03
2	rs12477314(239542085), HDAC4(downstream)	T	FEV1/FVC	0.052 (0.008)	4.48E-09	0.202	45585	0.031 (0.008)	8.41E-05	0.206	45821	0.041 (0.006)	1.68E-12
			FEV1	0.032 (0.008)	2.77E-04			0.025 (0.007)	1.82E-04			0.028 (0.005)	1.02E-07
3	rs1529672(25495586), RARB(intron)	C	FEV1/FVC	−0.06 (0.009)	7.75E-10	0.829	40624	−0.038 (0.009)	1.16E-05	0.831	45466	−0.048 (0.006)	3.97E-14
			FEV1	−0.037 (0.009)	1.78E-04			−0.011 (0.007)	9.33E-02			−0.02 (0.006)	2.16E-04
3	rs1344555(170782913), MECOM(intron)	T	FEV1/FVC	−0.019 (0.008)	2.61E-02	0.205	46067	−0.017 (0.012)	1.55E-01	0.209	21313	−0.018 (0.007)	6.65E-03
			FEV1	−0.042 (0.008)	1.91E-06			−0.025 (0.009)	6.44E-03			−0.034 (0.006)	2.65E-08
5	rs153916(95062456), SPATA9(upstream)	T	FEV1/FVC	−0.033 (0.007)	2.06E-06	0.552	47530	−0.025 (0.009)	6.67E-03	0.535	21647	−0.031 (0.005)	2.12E-08
			FEV1	−0.001 (0.007)	8.91E-01			0.004 (0.007)	6.22E-01			0.001 (0.005)	8.20E-01
6	rs6903823(28430275), ZKSCAN3(intron)/ZNF323(intron)	G	FEV1/FVC	−0.027 (0.008)	2.28E-03	0.209	47057	−0.013 (0.011)	2.34E-01	0.246	21489	−0.021 (0.007)	1.19E-03
			FEV1	−0.046 (0.008)	2.00E-07			−0.029 (0.008)	4.75E-04			−0.037 (0.006)	2.18E-10
6	rs2857595(31676448), NCR3(upstream)	G	FEV1/FVC	0.049 (0.009)	7.86E-08	0.809	45540	0.028 (0.008)	5.36E-04	0.796	46107	0.037 (0.006)	2.28E-10
			FEV1	0.04 (0.009)	1.46E-05			0.017 (0.007)	9.41E-03			0.025 (0.005)	1.30E-06
6	rs2798641(109374743), ARMC2(intron)	T	FEV1/FVC	−0.047 (0.009)	2.81E-07	0.183	46369	−0.03 (0.012)	1.57E-02	0.179	21173	−0.041 (0.007)	8.35E-09
			FEV1	−0.046 (0.009)	5.39E-07			−0.009 (0.01)	3.35E-01			−0.03 (0.006)	4.69E-06
10	rs7068966(12317998), CDC123(intron)	T	FEV1/FVC	0.045 (0.007)	1.28E-10	0.519	47085	0.023 (0.006)	3.86E-04	0.518	46067	0.033 (0.005)	6.13E-13
			FEV1	0.04 (0.007)	1.19E-08			0.022 (0.005)	3.56E-05			0.029 (0.004)	2.82E-12
10	rs11001819(77985230), C10orf11(intron)	G	FEV1/FVC	−0.019 (0.007)	6.50E-03	0.522	45546	−0.006 (0.006)	3.17E-01	0.506	45932	−0.012 (0.005)	7.58E-03
			FEV1	−0.041 (0.007)	1.42E-08			−0.022 (0.005)	3.10E-05			−0.029 (0.004)	2.98E-12
12	rs11172113(55813550), LRP1(intron)	T	FEV1/FVC	−0.035 (0.007)	1.36E-06	0.607	45387	−0.026 (0.01)	5.83E-03	0.59	20509	−0.032 (0.006)	1.24E-08
			FEV1	−0.021 (0.007)	3.55E-03			−0.003 (0.007)	6.94E-01			−0.013 (0.005)	1.19E-02
12	rs1036429(94795559), CCDC38(intron)	T	FEV1/FVC	0.049 (0.008)	1.24E-08	0.2	47814	0.028 (0.008)	3.35E-04	0.214	46311	0.038 (0.006)	2.30E-11
			FEV1	0.01 (0.008)	2.67E-01			0.004 (0.006)	5.38E-01			0.006 (0.005)	2.26E-01
16	rs12447804(56632783), MMP15(intron)	T	FEV1/FVC	−0.053 (0.009)	7.12E-08	0.208	35123	−0.021 (0.01)	4.20E-02	0.222	24398	−0.038 (0.007)	3.59E-08
			FEV1	−0.017 (0.009)	8.02E-02			0.004 (0.007)	5.71E-01			−0.004 (0.006)	4.73E-01
16	rs2865531(73947817), CFDP1(intron)	T	FEV1/FVC	0.039 (0.007)	2.30E-08	0.418	47594	0.024 (0.006)	1.94E-04	0.409	46304	0.031 (0.005)	1.77E-11
			FEV1	0.024 (0.007)	6.30E-04			0.011 (0.005)	3.89E-02			0.016 (0.004)	1.09E-04
21	rs9978142(34574109), KCNE2(upstream)	T	FEV1/FVC	−0.048 (0.009)	8.23E-07	0.156	44577	−0.031 (0.013)	1.75E-02	0.149	20944	−0.043 (0.008)	2.65E-08
			FEV1	−0.012 (0.009)	2.47E-01			−0.015 (0.01)	1.35E-01			−0.013 (0.007)	5.57E-02

Gene expression

We investigated mRNA expression of the nearest gene for each of the 16 novel loci in human lung tissue and in a range of human primary cells including lung, brain, airway smooth muscle cells and bronchial epithelial cells. Transcripts were detected for all selected genes in lung tissue except CCDC38 and transcripts for most genes were also detected in airway smooth muscle cells and in bronchial epithelial cells (Table 2). As we were unable to detect expression of CCDC38 in any tissue, we also examined expression of SNPRF, which is the adjacent gene (Table 2), and found expression in all four cell types. TGFB2, MFAP2, EVI1 and MMP15 were expressed in one or more lung cell types but not in peripheral blood mononuclear cells providing evidence that these genes may exhibit tissue-specific expression.

Table 2

Expression profiling of candidate genes in the lung and periphery. (+) indicates the gene is expressed in the cell type used and (−) indicates the gene expression at mRNA level is not detected following 40 cycles of PCR. PCR profiling of gene transcripts in the human lung demonstrates expression of all candidates except CCDC38, for which two sets of primers were designed and tested under different optimization conditions. None of these assays detected expression of CCDC38 in the cell types analysed. The SNRPF gene neighbouring CCDC38 and harbouring SNPs with strong LD with CCDC38’s sentinel SNP was assayed instead. All PCR products were sequence verified. Glyceraldehyde-3-phosphate dehydrogenase gene (GAPDH) was used as a positive control for the cDNA and was expressed in all tissues. HASM: Human airway smooth muscle. HBEC: Human bronchial epithelial cells. PBMC: peripheral blood mononuclear cells.

Sentinel SNP(relationshipto gene)	Chr.	Gene	Putative Function of encodedprotein	Tissue
				Lung	HASM	HBEC	PBMC
rs993925(intron)	1	TGFB2	Cytokine with roles in pro-fibroticcytokine modulating epithelial repairmechanisms and extracellular matrixhomeostasis including collagendeposition[40]	+	+	−	−
rs2284746(intron)	1	MFAP2	Major antigen of elastin-associatedmicrofibrils17 and a candidate forinvolvement in the etiology of inheritedconnective tissue diseases.	+	+	+	−
rs12477314(downstream)	2	HDAC4	Deacetylase of histone surrounding DNAthus influencing transcription factoraccess to the DNA possibly repressinggene transcription.	+	+	+	+
rs1344555(intron)	3	EVI1	Zinc finger transcription factor, encodedas part of MDS1/EVI1 complex locus(MECOM).	+	+	+	−
rs1529672(intron)	3	RARB	Nuclear retinoic acid receptor responsiveto retinoic acid, a vitamin A derivativeand which also controls cell proliferationand differentiation.	+	+	+	+
rs153916(intron)	5	SPATA9	Initially identified as a mediator ofspermatogenesis, other family membersmay have a role in pancreaticdevelopment and β-cell proliferation[41]	+	+	+	+
rs2798641(intron)	6	ARMC2	Function unknown although other familymembers have been identified as havingroles in cell signalling, proteindegradation and cytoskeleton functions[42]	+	+	+	+
rs2857595(upstream)	6	NCR3	Required for efficient cytotoxicityresponses by natural killer cells againstnormal cells and tumours[43]	+	−	−	+
rs6903823(intron)	6	ZKSCAN3	Transcription factor involved in cellgrowth/cell cycle/signal transduction	+	+	+	+
rs7068966(intron)	10	CDC123	Homologue in yeast shown to be acritical control protein modulatingEukaryotic initiation factor 2 in times ofcell stress	+	+	+	+
rs11001819(intron)	10	C10orf11	Function unknown	+	+	+	+
rs11172113(intron)	12	LRP1	Potentially diverse roles including cellsignalling and migration[44]	+	+	+	+
rs1036429(intron)	12	CCDC38	Function unknown although other familymembers involved in a diverse array offunctions skeletal and motor function[45]	−	−	−	−
rs1036429(r2=0.96 withrs4762633 inSNRPF)	12	SNRPF	Small nuclear ribonucleoprotein F	+	+	+	+
rs12447804(intron)	16	MMP15	Member of a large protease family withdiverse functional roles via proteaseactivity and specificity including; tissueremodelling, wound healing,angiogenesis, and tumor invasion.	+	+	+	−
rs2865531(intron)	16	CFDP1	Craniofacial Development Protein 1	+	+	+	+
rs9978142(upstream)	21	KCNE2	KCNQ1-KCNE2 K+ channels maymodulate transepithelial anion secretionin Calu3 airway epithelial cells[46].	+	−	−	+
Referencegene	12	GAPDH		+	+	+	+

We assessed whether SNPs in these new regions, or their proxies (r2>0.6), were associated with gene expression using a database of expression-associated SNPs in lymphoblastoid cell lines[16]. Four loci showed regional (cis) effects on expression (P<1×10−7, Supplementary Note). A proxy for our sentinel SNP in CFDP1, rs2865531, coincided with the peak of the expression signal for CFDP1 and the strongest proxy for rs6903823 in ZKSCAN3 coincided with the peak of expression for ZSCAN12.

Plausible pathways for lung function involving new loci

The putative function of the genes within, or closest to, the association peaks identify a range of plausible mechanisms for impacting lung function. The most statistically significant new signal for FEV1/FVC (P=7.5×10−16) was in the gene encoding MFAP2, an antigen of elastin-associated microfibrils[17], although correlated SNPs in the region potentially implicate other genes that could plausibly influence lung function, such as CROCC, which encodes rootletin, a component of cilia[18]. Our second strongest new signal, also for FEV1/FVC, was in the gene encoding the retinoic acid receptor beta (RARB). Rarb-null knockout mice exhibit premature alveolar septation[19]. The third most statistically significant new signal for FEV1/FVC, and the most statistically significant new signal for FEV1, was in CDC123. This was the only novel region to show genome-wide association with both traits. CDC123 encodes a homologue of a yeast cell division cycle protein which plays a critical role in modulating Eukaryotic initiation factor 2 in times of cell stress[20]. The fourth signal for FEV1/FVC is downstream of HDAC4 which encodes a histone deacetylase; reductions in the expression of other histone deacetylases (specifically HDAC2, HDAC5 and HDAC8) have been noted in COPD[21]. The regions we observed in the MHC are much more difficult to localize with multiple genes being tagged by the top SNP, including non-synonymous SNPs in ZKSCAN3, PGBD1, ZSCAN12, ZNF323, TCF19, LTA, C6orf15 and GPANK1 (also known as BAT4) (Supplementary Table 2). At 6p21.33, the strongest association with lung function was observed for rs2857595, which is in linkage disequilibrium (LD, r2=0.47) with a non-synonymous SNP in LTA (encoding lymphotoxin alpha) and with a SNP in the upstream promoter region of TNFA (encoding tumour necrosis factor alpha, r2=0.86), both of which are plausible candidates[22,23]. Our top SNP in MMP15 is in strong LD (r2=1) with a non-synonymous SNP (rs3743563, which has an association with FEV1/FVC at P=1.8×10−7) within the same gene. Plausible mechanisms implicated by the other novel signals of association with lung function reported here include TGF-beta signalling; TGFB2 expression is upregulated in bronchial epithelial cells in asthma[24]. The putative function of key genes (as defined by LD with the leading SNP) in each of the 16 loci, and relevant findings from animal models, are summarised in Table 2 and detailed in Supplementary Table 2.

Associations with lung function in children

Alleles representing 11 of the 16 novel loci showed directionally consistent effects on lung function in 6,281 children (7 to 9 years of age) (Supplementary Table 3A) suggesting that genetic determination of lung function in adults may in part act via effects on lung development, or alternatively, that some genetic determinants of lung growth and lung function decline are shared.

Association of lung function loci with other traits

Although we stratified for ever-smoking versus never-smoking, we did not adjust for the amount smoked. In order to investigate the possibility that the associations at any of our 16 novel regions were driven by an effect of the SNP on smoking behaviour, we evaluated in silico data for associations with smoking amount from the Ox-GSK consortium[25] for the leading SNPs in these 16 regions. None of these 16 SNPs showed statistically significant association with the number of cigarettes smoked per day (Supplementary Table 3B).

In addition, in our Stage 1 and Stage 2 datasets combined, we assessed whether the estimated effect sizes of the variants on lung function phenotypes differed substantially between ever-smokers and never-smokers (Supplementary Table 4) across the 16 loci. For the most strongly associated trait at each locus, we tested the SNP interaction with ever-smoking (versus never-smoking). None of the 16 novel loci showed a significant interaction (Bonferroni corrected threshold for 16 independent SNPs P=0.003125). These analyses suggest that the genetic effects we have identified underlie lung function variability irrespective of smoking exposure.

Our lung function associations were adjusted for height, but there are some overlaps between loci associated with height and those associated with lung function. Therefore, we evaluated in silico data for height associations of our novel regions in the GIANT consortium[14] dataset. The G allele of rs2284746 (MFAP2, intron), which was associated with decreased FEV1/FVC was associated with increased height (Supplementary Table 3C).

Given reported associations between lung cancer and either COPD or lung function decline, we also assessed in silico data for sentinel or proxy SNPs in these 16 regions for associations with lung cancer in the International Lung Cancer Consortium (ILCCO) GWAS meta-analysis[26]. Alleles associated with reduced lung function were associated with risk of lung cancer at the strongest available proxy SNP for rs2857595 (upstream of NCR3) at 6p21.33 (rs3099844, r2=0.67), and the strongest proxy SNP for rs6903823 (SNP in intron of ZKSCAN3 and ZNF323) at 6p22.1 (rs209181, r2=0.69) (lung cancer associations, P=2.2×10−7 and P=3.4×10−5, respectively, Supplementary Table 3D). No significant associations with lung cancer were seen at the other new loci (proxy SNPs were available for 15 of the 16 loci, Bonferroni corrected P<0.0033).

In addition to the effects on height, smoking and lung cancer described above, we examined the literature for evidence for associations with other traits for each of the 16 new loci (detailed in Supplementary Table 2). Genome-wide significant associations (P<5×10−8) have been reported in KCNE2 with myocardial infarction[27], and at 6p21.33 near NCR3/AIF1, with neonatal lupus[28] and with systemic lupus erythematosus[29]. Other significant complex disease associations have also been noted in the regions of CDC123 (type 2 diabetes[30]), CFDP1 (type 1 diabetes[31]) and MECOM (blood pressure[32,33]), but with weaker LD (r2<0.3) between the reported SNP and the sentinel SNP for lung function in the region (Supplementary Table 2).

Proportion of variance explained by loci discovered to date

Associations in 10 loci previously reported for lung function[8,9] reached genome-wide significance (P<5×10−8) in our Stage 1 data, namely loci in or near TNS1, FAM13A, GSTCD/NPNT, HHIP, HTR4, ADAM19, AGER, GPR126, PTCH1, and TSHD4 (Supplementary Table 5A). Thus, a total of 26 regions showed genome-wide significant association with lung function in our study. In aggregate, variants at these 26 regions explain approximately 3.2% of the additive polygenic variance for FEV1/FVC and 1.5% for FEV1 (see Supplementary Note). Following the approach described by Park et al.[34] we estimated that there is a total of 102 (95% CI 57-155) independent variants with similar effect sizes to the 26 variants we report. In combination these 102 variants, comprising 26 discovered variants and 76 putative undiscovered variants, collectively explain around 7.5% of the additive polygenic variance for FEV1/FVC and 3.4% for FEV1 (see Supplementary Table 6, Online Methods and Supplementary Note).

Discussion

In meta-analysis of 23 studies comprising 48,201 individuals of European ancestry and follow-up in 17 studies comprising up to 46,411 individuals, we report genome-wide significant associations with an additional 12 regions for FEV1/FVC, an additional three regions for FEV1 and one additional region associated with both FEV1 and FEV1/FVC. We also confirm genome-wide association with 10 regions previously associated with lung function, bringing to 26 the total number of loci associated with lung function in these data. Most of the new loci are in regions not previously suspected to have been involved in lung development, the control of pulmonary function or risk of developing COPD. Elucidating the mechanisms through which these regions influence lung function should lead to a more complete understanding of lung function regulation and the pathogenesis of COPD. Four of the new loci (MFAP2, ZKSCAN3, near NCR3 and near KCNE2) we show to be associated with lung function are also associated with other complex traits and diseases (P<5×10−8 for the other trait at a SNP with r2>0.3 with the top lung function SNP in the region). Understanding the intermediates underlying these pleiotropic effects could also reveal crucial insights into the pathophysiology of lung disease. One potential explanation is that these loci underlie control of the mechanisms regulating the development and resolution of inflammation and subsequent tissue remodelling in a range of tissues.

The effect sizes of the variants in the 26 loci associated with lung function explain a modest proportion of the additive genetic variance in FEV1/FVC and in FEV1, even after accounting for putative undetected variants with a similar distribution of effect sizes[34]. Our findings are consistent with those from other common complex traits, where it is thought that many as yet unidentified common and rare sequence variants, and potentially structural variants could explain the remaining heritability[35]. That our study more than doubled the number of loci known to be associated with lung function underlines the utility of large sample sizes to achieve the power to detect common variants associated with complex traits. Nevertheless, it is likely that additional variants with similar effect sizes remain undiscovered[14]. In addition, our study was not designed to detect rare variants or structural variants associated with lung function. Identification of rare variants associated with lung function could be helpful in narrowing the scope of ongoing functional work to those genes most likely to be causally related to the association signals we detected.

Our study focused on cross-sectional measures of lung function. Adult lung function at a particular time point is influenced by the peak lung function achieved by 25-35 years of age, as well as the rate of decline of lung function after that peak[36]. The 26 loci now confirmed to be associated with lung function could affect either pre- or post-natal lung development and growth or decline in lung function during adulthood, or both. We showed consistent directions of estimated effects on lung function between adults and children at 7-9 years of age for SNPs at 11 of the 16 new loci, and eight of 10 previously reported loci (Supplementary Table 3A). The results we show for lung function in children provide some indication that these loci affect lung function development, although studies in larger populations of children would provide greater clarity for SNPs in the new loci. Further investigations will be required in large populations with longitudinal data to delineate the influence of these variants on the rates of development of, and decline in, lung function and on the risk of developing COPD.

Of the sentinel SNPs at the 16 new loci associated with lung function, only rs2284746 (MFAP2) was associated with height in the GIANT consortium[14] dataset. The G allele of rs2284746 was associated with both increased height and reduced lung function. A similar relationship between lung function and height was previously reported for the G allele of rs3817928 in GPR126[8,14], which is associated with decreased height, but with increased FEV1/FVC. A further three of the 180 loci found to be associated with height[14] showed association (for 180 loci, Bonferroni corrected threshold P=2.8×10−4) with either FEV1 (CLIC4 and BMP6) or FEV1/FVC (PIP4K2B) (Supplementary Table 3E). In each case, the allele associated with an increase in height was associated with a decrease in lung function. This is not the case for the association of rs1032296 near HHIP, which has shown consistent directions of effects on lung function and height[14,11]. However, the strongest SNP associated with height in the HHIP region lies within an intron of HHIP but shows no association with FEV1 or FEV1/FVC. Furthermore, while height is an important predictor of FEV1, this is not true for its ratio to FVC[37]. These observations argue against the associations with lung function at these loci being simply due to incomplete adjustment for height.

We stratified by ever- and never-smoker status in our analyses and in our investigation of amount smoked in the Ox-GSK consortium[25] none of the sentinel SNPs in the 16 new regions showed association with the number of cigarettes smoked per day. Additionally, none of these regions was associated with ever-smoking in the Ox-GSK consortium data (Supplementary Table 3B). Thus the SNP associations with lung function we observed are unlikely to have arisen simply as a consequence of inadequate adjustment for smoking.

We did not observe any interactions with ever-smoking for any of the sentinel SNPs in the 16 new regions that exceeded a Bonferroni-corrected significance level (for 16 SNPs). Thus, the effects on lung function of the novel variants we identified are apparent in both ever-smokers and in never-smokers, and the effects of smoking and of these genetic variants may be independent and additive.

In other common complex diseases, follow up studies that incorporate common genetic risk variants into models to predict disease have not been shown to add substantially to existing risk models, particularly when such models already include family history[38,39]. The same may also prove to be true for the 26 genetic variants described in this paper, as the effect size of any individual variant is small, but further work is required in this area. The major utility of our findings will be in the knowledge they provide about previously unknown pathways underlying lung function. Elucidating the mechanisms that these genes are involved in will lead to improved understanding of the regulation of lung function and potentially to new therapeutic targets for COPD.