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.

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.

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.

The striking increase in COVID-19 severity in older adults provides a clear example of immunesenescence, the age-related remodelling of the immune system. To better characterise the association between convalescent immunesenescence and acute disease severity, we determined the immune phenotype of COVID-19 survivors and non-infected controls.

The study of families with rare inherited forms of hypo- and hyper-tension has been one of the most successful strategies to probe the molecular pathophysiology of blood pressure control and has revealed dysregulation of distal nephron Na+ reabsorption to be a common mechanism. FHHt (familial hyperkalaemic hypertension; also known as Gordon's syndrome) is a salt-dependent form of hypertension caused by mutations in the regulators of the thiazide-sensitive Na+-Cl- co-transporter NCC [also known as SLC12A3 (solute carrier family 12 member 3)] and is effectively treated by thiazide diuretics and/or dietary salt restriction. Variation in at least four genes can cause FHHt, including WNK1 [With No lysine (=K) 1] and WNK4, KLHL3 (kelch-like family member 3), and CUL3 (cullin 3). In the present study we have identified novel disease-causing variants in CUL3 and KLHL3 segregating in 63% of the pedigrees with previously unexplained FHHt, confirming the importance of these recently described FHHt genes. We have demonstrated conclusively, in two unrelated affected individuals, that rare intronic variants in CUL3 cause the skipping of exon 9 as has been proposed previously. KLHL3 variants all occur in kelch-repeat domains and so probably disrupt WNK complex binding. We have found no evidence of any plausible disease-causing variants within SLC4A8 (an alternative thiazide-sensitive sodium transporter) in this population. The results of the present study support the existing evidence that the CUL3 and KLHL3 gene products are physiologically important regulators of thiazide-sensitive distal nephron NaCl reabsorption, and hence potentially interesting novel anti-hypertensive drug targets. As a third of our non-WNK FHHt families do not have plausible CUL3 or KLHL3 variants, there are probably additional, as yet undiscovered, regulators of the thiazide-sensitive pathways.

Genome-Wide Association Study (GWAS) meta-analyses have identified a strong association signal for lung function, which maps to a region on 4q24 containing two oppositely transcribed genes: glutathione S-transferase, C-terminal domain containing (GSTCD) and integrator complex subunit 12 (INTS12). Both genes were found to be expressed in a range of human airway cell types. The promoter regions and transcription start sites were determined in mRNA from human lung and a novel splice variant was identified for each gene. We obtained the following evidence for GSTCD and INTS12 co-regulation and expression: (i) correlated mRNA expression was observed both via Q-PCR and in a lung expression quantitative trait loci (eQTL) study, (ii) induction of both GSTCD and INTS12 mRNA expression in human airway smooth muscle cells was seen in response to TGFβ1, (iii) a lung eQTL study revealed that both GSTCD and INTS12 mRNA levels positively correlate with percent predicted FEV1, and (iv) FEV1 GWAS associated SNPs in 4q24 were found to act as an eQTL for INTS12 in a number of tissues. In fixed sections of human lung tissue, GSTCD protein expression was ubiquitous, whereas INTS12 expression was predominantly in epithelial cells and pneumocytes. During human fetal lung development, GSTCD protein expression was observed to be highest at the earlier pseudoglandular stage (10-12 weeks) compared with the later canalicular stage (17-19 weeks), whereas INTS12 expression levels did not alter throughout these stages. Knowledge of the transcriptional and translational regulation and expression of GSTCD and INTS12 provides important insights into the potential role of these genes in determining lung function. Future work is warranted to fully define the functions of INTS12 and GSTCD.

Chronic obstructive pulmonary disease (COPD) is a leading cause of global morbidity and mortality and, whilst smoking remains the single most important risk factor, COPD risk is heritable. Of 26 independent genomic regions showing association with lung function in genome-wide association studies, eleven have been reported to show association with airflow obstruction. Although the main risk factor for COPD is smoking, some individuals are observed to have a high forced expired volume in 1 second (FEV1) despite many years of heavy smoking. We hypothesised that these "resistant smokers" may harbour variants which protect against lung function decline caused by smoking and provide insight into the genetic determinants of lung health. We undertook whole exome re-sequencing of 100 heavy smokers who had healthy lung function given their age, sex, height and smoking history and applied three complementary approaches to explore the genetic architecture of smoking resistance. Firstly, we identified novel functional variants in the "resistant smokers" and looked for enrichment of these novel variants within biological pathways. Secondly, we undertook association testing of all exonic variants individually with two independent control sets. Thirdly, we undertook gene-based association testing of all exonic variants. Our strongest signal of association with smoking resistance for a non-synonymous SNP was for rs10859974 (P = 2.34 × 10(-4)) in CCDC38, a gene which has previously been reported to show association with FEV1/FVC, and we demonstrate moderate expression of CCDC38 in bronchial epithelial cells. We identified an enrichment of novel putatively functional variants in genes related to cilia structure and function in resistant smokers. Ciliary function abnormalities are known to be associated with both smoking and reduced mucociliary clearance in patients with COPD. We suggest that genetic influences on the development or function of cilia in the bronchial epithelium may affect growth of cilia or the extent of damage caused by tobacco smoke.

In this trans-ethnic multi-omic study, we reinterpret the genetic architecture of blood pressure to identify genes, tissues, phenomes and medication contexts of blood pressure homeostasis. We discovered 208 novel common blood pressure SNPs and 53 rare variants in genome-wide association studies of systolic, diastolic and pulse pressure in up to 776,078 participants from the Million Veteran Program (MVP) and collaborating studies, with analysis of the blood pressure clinical phenome in MVP. Our transcriptome-wide association study detected 4,043 blood pressure associations with genetically predicted gene expression of 840 genes in 45 tissues, and mouse renal single-cell RNA sequencing identified upregulated blood pressure genes in kidney tubule cells.

Eosinophils are associated with airway inflammation in respiratory disease. Eosinophil production and survival is controlled partly by interleukin-5: anti-interleukin-5 agents reduce asthma and response correlates with baseline eosinophil counts. However, whether raised eosinophils are causally related to chronic obstructive pulmonary disease (COPD) and other respiratory phenotypes is not well understood.