Understanding the genetic basis of airflow obstruction and smoking behaviour is key to determining the pathophysiology of chronic obstructive pulmonary disease (COPD). We used UK Biobank data to study the genetic causes of smoking behaviour and lung health.

Chronic obstructive pulmonary disease (COPD) is the leading cause of respiratory mortality worldwide. Genetic risk loci provide new insights into disease pathogenesis. We performed a genome-wide association study in 35,735 cases and 222,076 controls from the UK Biobank and additional studies from the International COPD Genetics Consortium. We identified 82 loci associated with P < 5 × 10-8; 47 of these were previously described in association with either COPD or population-based measures of lung function. Of the remaining 35 new loci, 13 were associated with lung function in 79,055 individuals from the SpiroMeta consortium. Using gene expression and regulation data, we identified functional enrichment of COPD risk loci in lung tissue, smooth muscle, and several lung cell types. We found 14 COPD loci shared with either asthma or pulmonary fibrosis. COPD genetic risk loci clustered into groups based on associations with quantitative imaging features and comorbidities. Our analyses provide further support for the genetic susceptibility and heterogeneity of COPD.

Background: Over 90 regions of the genome have been associated with lung function to date, many of which have also been implicated in chronic obstructive pulmonary disease. Methods: We carried out meta-analyses of exome array data and three lung function measures: forced expiratory volume in one second (FEV 1), forced vital capacity (FVC) and the ratio of FEV 1 to FVC (FEV 1/FVC). These analyses by the SpiroMeta and CHARGE consortia included 60,749 individuals of European ancestry from 23 studies, and 7,721 individuals of African Ancestry from 5 studies in the discovery stage, with follow-up in up to 111,556 independent individuals. Results: We identified significant (P<2·8x10 -7) associations with six SNPs: a nonsynonymous variant in RPAP1, which is predicted to be damaging, three intronic SNPs ( SEC24C, CASC17 and UQCC1) and two intergenic SNPs near to LY86 and FGF10. Expression quantitative trait loci analyses found evidence for regulation of gene expression at three signals and implicated several genes, including TYRO3 and PLAU. Conclusions: Further interrogation of these loci could provide greater understanding of the determinants of lung function and pulmonary disease.

Genome-wide association studies (GWAS) identified the chromosome 15q25.1 locus as a leading susceptibility region for lung cancer. However, the pathogenic pathways, through which susceptibility SNPs within chromosome 15q25.1 affects lung cancer risk, have not been explored. We analyzed three cohorts with GWAS data consisting 42,901 individuals and lung expression quantitative trait loci (eQTL) data on 409 individuals to identify and validate the underlying pathways and to investigate the combined effect of genes from the identified susceptibility pathways. The KEGG neuroactive ligand receptor interaction pathway, two Reactome pathways, and 22 Gene Ontology terms were identified and replicated to be significantly associated with lung cancer risk, with P values less than 0.05 and FDR less than 0.1. Functional annotation of eQTL analysis results showed that the neuroactive ligand receptor interaction pathway and gated channel activity were involved in lung cancer risk. These pathways provide important insights for the etiology of lung cancer.

Chronic obstructive pulmonary disease (COPD) independently associates with an increased risk of coronary artery disease (CAD), but it has not been fully investigated whether this co-morbidity involves shared pathophysiological mechanisms. To identify potential common pathways across the two diseases, we tested all recently published single nucleotide polymorphisms (SNPs) associated with human lung function (spirometry) for association with carotid intima-media thickness (cIMT) in 3,378 subjects with multiple CAD risk factors, and for association with CAD in a case-control study of 5,775 CAD cases and 7,265 controls. SNPs rs2865531, located in the CFDP1 gene, and rs9978142, located in the KCNE2 gene, were significantly associated with CAD. In addition, SNP rs9978142 and SNP rs3995090 located in the HTR4 gene, were associated with average and maximal cIMT measures. Genetic risk scores combining the most robustly spirometry-associated SNPs from the literature were modestly associated with CAD, (odds ratio (OR) (95% confidence interval (CI95) = 1.06 (1.03, 1.09); P-value = 1.5 × 10(-4), per allele). In conclusion, our study suggests that some genetic loci implicated in determining human lung function also influence cIMT and susceptibility to CAD. The present results should help elucidate the molecular underpinnings of the co-morbidity observed across COPD and CAD.

Deep sequence-based imputation can enhance the discovery power of genome-wide association studies by assessing previously unexplored variation across the common- and low-frequency spectra. We applied a hybrid whole-genome sequencing (WGS) and deep imputation approach to examine the broader allelic architecture of 12 anthropometric traits associated with height, body mass, and fat distribution in up to 267,616 individuals. We report 106 genome-wide significant signals that have not been previously identified, including 9 low-frequency variants pointing to functional candidates. Of the 106 signals, 6 are in genomic regions that have not been implicated with related traits before, 28 are independent signals at previously reported regions, and 72 represent previously reported signals for a different anthropometric trait. 71% of signals reside within genes and fine mapping resolves 23 signals to one or two likely causal variants. We confirm genetic overlap between human monogenic and polygenic anthropometric traits and find signal enrichment in cis expression QTLs in relevant tissues. Our results highlight the potential of WGS strategies to enhance biologically relevant discoveries across the frequency spectrum.

Thiazide diuretics are among the most widely used treatments for hypertension, but thiazide-induced hyponatremia (TIH), a clinically significant adverse effect, is poorly understood. Here, we have studied the phenotypic and genetic characteristics of patients hospitalized with TIH. In a cohort of 109 TIH patients, those with severe TIH displayed an extended phenotype of intravascular volume expansion, increased free water reabsorption, urinary prostaglandin E2 excretion, and reduced excretion of serum chloride, magnesium, zinc, and antidiuretic hormone. GWAS in a separate cohort of 48 TIH patients and 2,922 controls from the 1958 British birth cohort identified an additional 14 regions associated with TIH. We identified a suggestive association with a variant in SLCO2A1, which encodes a prostaglandin transporter in the distal nephron. Resequencing of SLCO2A1 revealed a nonsynonymous variant, rs34550074 (p.A396T), and association with this SNP was replicated in a second cohort of TIH cases. TIH patients with the p.A396T variant demonstrated increased urinary excretion of prostaglandin E2 and metabolites. Moreover, the SLCO2A1 phospho-mimic p.A396E showed loss of transporter function in vitro. These findings indicate that the phenotype of TIH involves a more extensive metabolic derangement than previously recognized. We propose one mechanism underlying TIH development in a subgroup of patients in which SLCO2A1 regulation is altered.

Genome-wide association studies have identified numerous genetic loci for spirometic measures of pulmonary function, forced expiratory volume in one second (FEV(1)), and its ratio to forced vital capacity (FEV(1)/FVC). Given that cigarette smoking adversely affects pulmonary function, we conducted genome-wide joint meta-analyses (JMA) of single nucleotide polymorphism (SNP) and SNP-by-smoking (ever-smoking or pack-years) associations on FEV(1) and FEV(1)/FVC across 19 studies (total N = 50,047). We identified three novel loci not previously associated with pulmonary function. SNPs in or near DNER (smallest P(JMA = )5.00×10(-11)), HLA-DQB1 and HLA-DQA2 (smallest P(JMA = )4.35×10(-9)), and KCNJ2 and SOX9 (smallest P(JMA = )1.28×10(-8)) were associated with FEV(1)/FVC or FEV(1) in meta-analysis models including SNP main effects, smoking main effects, and SNP-by-smoking (ever-smoking or pack-years) interaction. The HLA region has been widely implicated for autoimmune and lung phenotypes, unlike the other novel loci, which have not been widely implicated. We evaluated DNER, KCNJ2, and SOX9 and found them to be expressed in human lung tissue. DNER and SOX9 further showed evidence of differential expression in human airway epithelium in smokers compared to non-smokers. Our findings demonstrated that joint testing of SNP and SNP-by-environment interaction identified novel loci associated with complex traits that are missed when considering only the genetic main effects.

Several infrequent genetic polymorphisms in the SERPINA1 gene are known to substantially reduce concentration of alpha1-antitrypsin (AAT) in the blood. Since low AAT serum levels fail to protect pulmonary tissue from enzymatic degradation, these polymorphisms also increase the risk for early onset chronic obstructive pulmonary disease (COPD). The role of more common SERPINA1 single nucleotide polymorphisms (SNPs) in respiratory health remains poorly understood. We present here an agnostic investigation of genetic determinants of circulating AAT levels in a general population sample by performing a genome-wide association study (GWAS) in 1392 individuals of the SAPALDIA cohort. Five common SNPs, defined by showing minor allele frequencies (MAFs) >5%, reached genome-wide significance, all located in the SERPINA gene cluster at 14q32.13. The top-ranking genotyped SNP rs4905179 was associated with an estimated effect of β = -0.068 g/L per minor allele (P = 1.20*10(-12)). But denser SERPINA1 locus genotyping in 5569 participants with subsequent stepwise conditional analysis, as well as exon-sequencing in a subsample (N = 410), suggested that AAT serum level is causally determined at this locus by rare (MAF<1%) and low-frequent (MAF 1-5%) variants only, in particular by the well-documented protein inhibitor S and Z (PI S, PI Z) variants. Replication of the association of rs4905179 with AAT serum levels in the Copenhagen City Heart Study (N = 8273) was successful (P<0.0001), as was the replication of its synthetic nature (the effect disappeared after adjusting for PI S and Z, P = 0.57). Extending the analysis to lung function revealed a more complex situation. Only in individuals with severely compromised pulmonary health (N = 397), associations of common SNPs at this locus with lung function were driven by rarer PI S or Z variants. Overall, our meta-analysis of lung function in ever-smokers does not support a functional role of common SNPs in the SERPINA gene cluster in the general population.

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10(-7) study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r(2) = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10(-7) at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies.

Genetic factors influence chronic obstructive pulmonary disease (COPD) risk, but the individual variants that have been identified have small effects. We hypothesised that a polygenic risk score using additional variants would predict COPD and associated phenotypes.

Whether smoking-associated DNA methylation has a causal effect on lung function has not been thoroughly evaluated. We first investigated the causal effects of 474 smoking-associated CpGs on forced expiratory volume in 1 s (FEV1) in UK Biobank (n = 321,047) by using two-sample Mendelian randomization (MR) and then replicated this investigation in the SpiroMeta Consortium (n = 79,055). Second, we used two-step MR to investigate whether DNA methylation mediates the effect of smoking on FEV1. Lastly, we evaluated the presence of horizontal pleiotropy and assessed whether there is any evidence for shared causal genetic variants between lung function, DNA methylation, and gene expression by using a multiple-trait colocalization ("moloc") framework. We found evidence of a possible causal effect for DNA methylation on FEV1 at 18 CpGs (p < 1.2 × 10-4). Replication analysis supported a causal effect at three CpGs (cg21201401 [LIME1 and ZGPAT], cg19758448 [PGAP3], and cg12616487 [EML3 and AHNAK] [p < 0.0028]). DNA methylation did not clearly mediate the effect of smoking on FEV1, although DNA methylation at some sites might influence lung function via effects on smoking. By using "moloc", we found evidence of shared causal variants between lung function, gene expression, and DNA methylation. These findings highlight potential therapeutic targets for improving lung function and possibly smoking cessation, although larger, tissue-specific datasets are required to confirm these results.

Considerable clinical heterogeneity in idiopathic pulmonary fibrosis (IPF) suggests the existence of multiple disease endotypes. Identifying these endotypes would improve our understanding of the pathogenesis of IPF and could allow for a biomarker-driven personalised medicine approach. We aimed to identify clinically distinct groups of patients with IPF that could represent distinct disease endotypes.

DeepPheWAS is an R package for phenome-wide association studies that creates clinically curated composite phenotypes and integrates quantitative phenotypes from primary care data, longitudinal trajectories of quantitative measures, disease progression and drug response phenotypes. Tools are provided for efficient analysis of association with any genetic input, under any genetic model, with optional sex-stratified analysis, and for developing novel phenotypes.

Rationale: Idiopathic pulmonary fibrosis (IPF) is a devastating disease characterized by limited treatment options and high mortality. A better understanding of the molecular drivers of IPF progression is needed. Objectives: To identify and validate molecular determinants of IPF survival. Methods: A staged genome-wide association study was performed using paired genomic and survival data. Stage I cases were drawn from centers across the United States and Europe and stage II cases from Vanderbilt University. Cox proportional hazards regression was used to identify gene variants associated with differential transplantation-free survival (TFS). Stage I variants with nominal significance (P < 5 × 10-5) were advanced for stage II testing and meta-analyzed to identify those reaching genome-wide significance (P < 5 × 10-8). Downstream analyses were performed for genes and proteins associated with variants reaching genome-wide significance. Measurements and Main Results: After quality controls, 1,481 stage I cases and 397 stage II cases were included in the analysis. After filtering, 9,075,629 variants were tested in stage I, with 158 meeting advancement criteria. Four variants associated with TFS with consistent effect direction were identified in stage II, including one in an intron of PCSK6 (proprotein convertase subtilisin/kexin type 6) reaching genome-wide significance (hazard ratio, 4.11 [95% confidence interval, 2.54-6.67]; P = 9.45 × 10-9). PCSK6 protein was highly expressed in IPF lung parenchyma. PCSK6 lung staining intensity, peripheral blood gene expression, and plasma concentration were associated with reduced TFS. Conclusions: We identified four novel variants associated with IPF survival, including one in PCSK6 that reached genome-wide significance. Downstream analyses suggested that PCSK6 protein plays a potentially important role in IPF progression.

Many adults hospitalised with COVID-19 have persistent symptoms such as fatigue, breathlessness and brain fog that limit day-to-day activities. These symptoms can last over 2 years. Whilst there is limited controlled studies on interventions that can support those with ongoing symptoms, there has been some promise in rehabilitation interventions in improving function and symptoms either using face-to-face or digital methods, but evidence remains limited and these studies often lack a control group.

Although several genetic associations with scleroderma (SSc) are defined, very little is known on genetic susceptibility to SSc-associated interstitial lung disease (SSc-ILD). A number of common polymorphisms have been associated with SSc-ILD, but most have not been replicated in separate populations. Four SNPs in IRF5, and one in each of STAT4, CD226 and IRAK1, selected as having been previously the most consistently associated with SSc-ILD, were genotyped in 612 SSc patients, of European descent, of whom 394 had ILD. The control population (n = 503) comprised individuals of European descent from the 1000 Genomes Project. After Bonferroni correction, two of the IRF5 SNPs, rs2004640 (OR (95% CI)1.30 (1.10-1.54), pcorr = 0.015) and rs10488631 (OR 1.48 (1.14-1.92), pcorr = 0.022), and the STAT4 SNP rs7574865 (OR 1.43 (1.18-1.73), pcorr = 0.0015) were significantly associated with SSc compared with controls. However, none of the SNPs were significantly different between patients with SSc-ILD and controls. Two SNPs in IRF5, rs10488631 (OR 1.72 (1.24-2.39), pcorr = 0.0098), and rs2004640 (OR 1.39 (1.11-1.75), pcorr = 0.03), showed a significant difference in allele frequency between controls and patients without ILD, as did STAT4 rs7574865 (OR 1.86 (1.45-2.38), pcorr = 6.6 × 10-6). A significant difference between SSc with and without ILD was only observed for STAT4 rs7574865, being less frequent in patients with ILD (OR 0.66 (0.51-0.85), pcorr = 0.0084). In conclusion, IRF5 rs2004640 and rs10488631, and STAT4 rs7574865 were significantly associated with SSc as a whole. Only STAT4 rs7574865 showed a significant difference in allele frequency in SSc-ILD, with the T allele being protective against ILD.Key points• We confirm the associations of the IRF5 SNPs rs2004640 and rs10488631, and the STAT4 SNP rs7574865, with SSc as a whole.• None of the tested SNPs were risk factors for SSc-ILD specifically.• The STAT4 rs7574865 T allele was protective against the development of lung fibrosis in SSc patients.• Further work is required to understand the genetic basis of lung fibrosis in association with scleroderma.

Specific common and rare single nucleotide variants (SNVs) increase the likelihood of developing sporadic idiopathic pulmonary fibrosis (IPF). We performed target-enriched sequencing on three loci previously identified by a genome-wide association study to gain a deeper understanding of the full spectrum of IPF genetic risk and performed a two-stage case-control association study.

Homozygosity for the SERPINA1 Z allele causes α1-antitrypsin deficiency, a rare condition that can cause lung and liver disease. However, the effects of Z allele heterozygosity on nonrespiratory phenotypes, and on lung function in the general population, remain unclear. We conducted a large, population-based study to determine Z allele effects on >2400 phenotypes in the UK Biobank (N=303 353). Z allele heterozygosity was strongly associated with increased height (β=1.02 cm, p=3.91×10-68), and with other nonrespiratory phenotypes including increased risk of gall bladder disease, reduced risk of heart disease and lower blood pressure, reduced risk of osteoarthritis and reduced bone mineral density, increased risk of headache and enlarged prostate, as well as with blood biomarkers of liver function. Heterozygosity was associated with higher height-adjusted forced expiratory volume in 1 s (FEV1) (β=19.36 mL, p=9.21×10-4) and FEV1/forced vital capacity (β=0.0031, p=1.22×10-5) in nonsmokers, whereas in smokers, this protective effect was abolished. Furthermore, we show for the first time that sex modifies the association of the Z allele on lung function. We conclude that Z allele heterozygosity and homozygosity exhibit opposing effects on lung function in the UK population, and that these associations are modified by smoking and sex. In exploratory analyses, heterozygosity for the Z allele also showed pleiotropic associations with nonrespiratory health-related traits and disease risk.

Gastro-oesophageal reflux disease (GORD) is associated with idiopathic pulmonary fibrosis (IPF) in observational studies. It is not known if this association arises because GORD causes IPF or because IPF causes GORD, or because of confounding by factors, such as smoking, associated with both GORD and IPF. We used bidirectional Mendelian randomisation (MR), where genetic variants are used as instrumental variables to address issues of confounding and reverse causation, to examine how, if at all, GORD and IPF are causally related.