Search

Summary 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. Availability and implementation The DeepPheWAS R package is freely available under GNU general public licence v3.0 from at https://github.com/Richard-Packer/DeepPheWAS. Supplementary information Supplementary data are available at Bioinformatics online.

Chronic Obstructive Pulmonary Disease (COPD) has a simple physiological diagnostic criterion but a wide range of clinical characteristics. The mechanisms underlying this variability in COPD phenotypes are unclear. To investigate the potential contribution of genetic variants to phenotypic heterogeneity, we examined the association of genome-wide associated lung function, COPD, and asthma variants with other phenotypes using phenome-wide association results derived in the UK Biobank. Our clustering analysis of the variants-phenotypes association matrix identified three clusters of genetic variants with different effects on white blood cell counts, height, and body mass index (BMI). To assess the potential clinical and molecular effects of these groups of variants, we investigated the association between cluster-specific genetic risk scores and phenotypes in the COPDGene cohort. We observed differences in steroid use, BMI, lymphocyte counts, chronic bronchitis, and differential gene and protein expression across the three genetic risk scores. Our results suggest that multi-phenotype analysis of obstructive lung disease-related risk variants may identify genetically driven phenotypic patterns in COPD.

Treatment-resistant depression (TRD) is a major contributor to the disability caused by major depressive disorder (MDD). Primary care electronic health records provide an easily accessible approach to investigate TRD clinical and genetic characteristics. MDD defined from primary care records in UK Biobank (UKB) and EXCEED studies was compared with other measures of depression and tested for association with MDD polygenic risk score (PRS). Using prescribing records, TRD was defined from at least two switches between antidepressant drugs, each prescribed for at least 6 weeks. Clinical-demographic characteristics, SNP-based heritability (h2SNP) and genetic overlap with psychiatric and non-psychiatric traits were compared in TRD and non-TRD MDD cases. In 230,096 and 8926 UKB and EXCEED participants with primary care data, respectively, the prevalence of MDD was 8.7% and 14.2%, of which 13.2% and 13.5% was TRD, respectively. In both cohorts, MDD defined from primary care records was strongly associated with MDD PRS, and in UKB it showed overlap of 71–88% with other MDD definitions. In UKB, TRD vs healthy controls and non-TRD vs healthy controls h2SNP was comparable (0.25 [SE = 0.04] and 0.19 [SE = 0.02], respectively). TRD vs non-TRD was positively associated with the PRS of attention deficit hyperactivity disorder, with lower socio-economic status, obesity, higher neuroticism and other unfavourable clinical characteristics. This study demonstrated that MDD and TRD can be reliably defined using primary care records and provides the first large scale population assessment of the genetic, clinical and demographic characteristics of TRD.

Background: Asthma is a complex disease with multiple phenotypes that may differ in disease pathobiology and treatment response. IL33 single nucleotide polymorphisms (SNPs) have been reproducibly associated with asthma. IL33 levels are elevated in sputum and bronchial biopsies of patients with asthma. The functional consequences of IL33 asthma SNPs remain unknown. Objective: This study sought to determine whether IL33 SNPs associate with asthma-related phenotypes and with IL33 expression in lung or bronchial epithelium. This study investigated the effect of increased IL33 expression on human bronchial epithelial cell (HBEC) function. Methods: Association between IL33 SNPs (Chr9: 5,815,786-6,657,983) and asthma phenotypes (Lifelines/DAG [Dutch Asthma GWAS]/GASP [Genetics of Asthma Severity & Phenotypes] cohorts) and between SNPs and expression (lung tissue, bronchial brushes, HBECs) was done using regression modeling. Lentiviral overexpression was used to study IL33 effects on HBECs. Results: We found that 161 SNPs spanning the IL33 region associated with 1 or more asthma phenotypes after correction for multiple testing. We report a main independent signal tagged by rs992969 associating with blood eosinophil levels, asthma, and eosinophilic asthma. A second, independent signal tagged by rs4008366 presented modest association with eosinophilic asthma. Neither signal associated with FEV 1, FEV 1/forced vital capacity, atopy, and age of asthma onset. The 2 IL33 signals are expression quantitative loci in bronchial brushes and cultured HBECs, but not in lung tissue. IL33 overexpression in vitro resulted in reduced viability and reactive oxygen species–capturing of HBECs, without influencing epithelial cell count, metabolic activity, or barrier function. Conclusions: We identify IL33 as an epithelial susceptibility gene for eosinophilia and asthma, provide mechanistic insight, and implicate targeting of the IL33 pathway specifically in eosinophilic asthma.

Lung function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry GWAS meta-analysis of lung function to date, comprising 580,869 participants, 1020 independent association signals identified 559 genes supported by ≥2 criteria from a systematic variant-to-gene mapping framework. These genes were enriched in 29 pathways. Individual variants showed heterogeneity across ancestries, age and smoking groups, and collectively as a genetic risk score (GRS) showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies (PheWAS) for selected associated variants, and trait and pathway-specific GRS to infer possible consequences of intervening in pathways underlying lung function. We highlight new putative causal variants, genes, proteins and pathways, including those targeted by existing drugs. These findings bring us closer to understanding the mechanisms underlying lung function and COPD, and should inform functional genomics experiments and potentially future COPD therapies.

Background: Some people have characteristics of both asthma and COPD (asthma-COPD overlap), and evidence suggests they experience worse outcomes than those with either condition alone.Research Question: What is the genetic architecture of asthma-COPD overlap, and do the determinants of risk for asthma-COPD overlap differ from those for COPD or asthma?Study Design and Methods: We conducted a genome-wide association study in 8,068 asthma-COPD overlap case subjects and 40,360 control subjects without asthma or COPD of European ancestry in UK Biobank (stage 1). We followed up promising signals (P < 5 × 10–6) that remained associated in analyses comparing (1) asthma-COPD overlap vs asthma-only control subjects, and (2) asthma-COPD overlap vs COPD-only control subjects. These variants were analyzed in 12 independent cohorts (stage 2).Results: We selected 31 independent variants for further investigation in stage 2, and discovered eight novel signals (P < 5 × 10–8) for asthma-COPD overlap (meta-analysis of stage 1 and 2 studies). These signals suggest a spectrum of shared genetic influences, some predominantly influencing asthma (FAM105A, GLB1, PHB, TSLP), others predominantly influencing fixed airflow obstruction (IL17RD, C5orf56, HLA-DQB1). One intergenic signal on chromosome 5 had not been previously associated with asthma, COPD, or lung function. Subgroup analyses suggested that associations at these eight signals were not driven by smoking or age at asthma diagnosis, and in phenome-wide scans, eosinophil counts, atopy, and asthma traits were prominent.Interpretation: We identified eight signals for asthma-COPD overlap, which may represent loci that predispose to type 2 inflammation, and serious long-term consequences of asthma.

Risk of chronic obstructive pulmonary disease (COPD) is determined by both cigarette smoking and genetic susceptibility, but little is known about gene-by-smoking interactions. We performed a genome-wide association analysis of 179,689 controls and 21,077 COPD cases from UK Biobank subjects of European ancestry recruited from 2006 to 2010, considering genetic main effects and gene-by-smoking interaction effects simultaneously (2-degrees-of-freedom (df) test) as well as interaction effects alone (1-df interaction test). We sought to replicate significant results in COPDGene (United States, 2008–2010) and SpiroMeta Consortium (multiple countries, 1947–2015) data. We considered 2 smoking variables: 1) ever/never and 2) current/noncurrent. In the 1-df test, we identified 1 genome-wide significant locus on 15q25.1 (cholinergic receptor nicotinic β4 subunit, or CHRNB4) for ever- and current smoking and identified PI*Z allele (rs28929474) of serpin family A member 1 (SERPINA1) for ever-smoking and 3q26.2 (MDS1 and EVI1 complex locus, or MECOM) for current smoking in an analysis of previously reported COPD loci. In the 2-df test, most of the significant signals were also significant for genetic marginal effects, aside from 16q22.1 (sphingomyelin phosphodiesterase 3, or SMPD3) and 19q13.2 (Egl-9 family hypoxia inducible factor 2, or EGLN2). The significant effects at 15q25.1 and 19q13.2 loci, both previously described in prior genome-wide association studies of COPD or smoking, were replicated in COPDGene and SpiroMeta. We identified interaction effects at previously reported COPD loci; however, we failed to identify novel susceptibility loci.

Background: 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. Methods: We constructed a polygenic risk score using a genome-wide association study of lung function (FEV1 and FEV1/forced vital capacity [FVC]) from the UK Biobank and SpiroMeta. We tested this polygenic risk score in nine cohorts of multiple ethnicities for an association with moderate-to-severe COPD (defined as FEV1/FVC 1

Genome-wide association studies have identified numerous common genetic variants associated with spirometric measures of pulmonary function, including forced expiratory volume in one second (FEV1), forced vital capacity, and their ratio. However, variants with lower minor allele frequencies are less explored. We conducted a large-scale gene-smoking interaction meta-analysis on exonic rare and low-frequency variants involving 44,429 individuals of European ancestry in the discovery stage and sought replication in the UK BiLEVE study with 45,133 European ancestry samples and UK Biobank study with 59,478 samples. We leveraged data on cigarette smoking, the major environmental risk factor for reduced lung function, by testing gene-by-smoking interaction effects only and simultaneously testing the genetic main effects and interaction effects. The most statistically significant signal that replicated was a previously reported low-frequency signal in GPR126, distinct from common variant associations in this gene. Although only nominal replication was obtained for a top rare variant signal rs142935352 in one of the two studies, interaction and joint tests for current smoking and PDE3B were significantly associated with FEV1. This study investigates the utility of assessing gene-by-smoking interactions and underscores their effects on potential pulmonary function.

Background and aimsNon-alcoholic fatty liver disease (NAFLD) is a complex trait that has a global prevalence estimated as 25%. We aimed to identify the genetic variant underlying a four-generation family with progressive NAFLD leading to cirrhosis, decompensation and development of hepatocellular carcinoma in the absence of common risk factors such as obesity and type 2 diabetes.MethodsExome sequencing and genome comparisons were used to identify the likely causal variant. We extensively characterised the clinical phenotype and post-prandial metabolic responses of family members with the identified novel variant in comparison to healthy non-carriers and wild type patients with NAFLD. Variant-expressing hepatocyte-like cells (HLCs) were derived from human induced pluripotent stem cells generated from homozygous donor skin fibroblasts. The phenotype was assessed using imaging, targeted RNA analysis and molecular expression arrays.ResultsWe identified a rare causal variant in MTTP, c.1691T>C p.I564T (rs745447480) encoding microsomal triglyceride transfer protein (MTP) associated with progressive non-alcoholic fatty liver disease, unrelated to metabolic syndrome. Although other described mutations in MTTP cause abetalipoproteinemia, neither homozygotes nor heterozygotes exhibited characteristic manifestations of this severe disease. HLCs derived from a homozygote donor had lower lipoprotein ApoB secretion, compared to wild type cells. Cytoplasmic triglyceride accumulation in HLCs triggered endoplasmic reticulum stress, secretion of pro-inflammatory mediators and production of reactive oxygen species.ConclusionWe have identified and characterized a rare causal variant in MTTP and homozygosity for MTTP p.I564T is associated with progressive NAFLD without any other manifestations of abetalipoproteinemia.

Impaired lung function is predictive of mortality and is a key component in the diagnosis of chronic obstructive pulmonary disease. Lung function has a strong genetic component but is also affected by environmental factors such as increased exposure to air pollution. How genetic factors and air pollution interact to affect lung function is however less understood.We conducted a genome-wide gene-air pollution interaction analysis of spirometry measures with three measures of air pollution at home address: particulate matter (PM2.5& PM10) and nitrogen dioxide (NO2), in approximately 300,000 unrelated European individuals from UK Biobank. We explored air pollution interactions with previously identified lung function signals and determined their combined interaction effect using a polygenic risk score (PRS).We identified seven genome-wide interaction signals (P< 5 × 10−8), and a further ten suggestive interaction signals (P< 5 × 10−7). We found statistical evidence of interaction with PM2.5for previous lung function signal, rs10841302, nearAEBP2, suggesting increased susceptibility of FEV1/FVC to PM2.5, as copies of the G allele increased (interaction beta: −0.073 percentage points, 95%CI: −0.105, −0.041). There was no observed interaction between air pollutants and the weighted genetic risk score.We carried out the largest genome-wide gene-air pollution interaction study of lung function and identified effects of clinically relevant size and significance. We observed up to 440ml lower lung function for certain genotypes associated with mean levels of outdoor air pollution at baseline, which is approximately equivalent to nine years of normal loss of lung function.

Reduced lung function predicts mortality and is key to the diagnosis of chronic obstructive pulmonary disease (COPD). In a genome-wide association study in 400,102 individuals of European ancestry, we define 279 lung function signals, 139 of which are new. In combination, these variants strongly predict COPD in independent patient populations. Furthermore, the combined effect of these variants showed generalizability across smokers and never-smokers, and across ancestral groups. We highlight biological pathways, known and potential drug targets for COPD and, in phenome-wide association studies, autoimmune-related and other pleiotropic effects of lung function associated variants. This new genetic evidence has potential to improve future preventive and therapeutic strategies for COPD., Editorial summary: A genome-wide association study in over 400,000 individuals identifies 139 new signals for lung function. These variants can predict chronic obstructive pulmonary disease in independent, trans-ethnic cohorts.

Variability in SARS-CoV-2 susceptibility and COVID-19 disease severity between individuals is partly due to genetic factors. Here, we identify 4 genomic loci with suggestive associations for SARS-CoV-2 susceptibility and 19 for COVID-19 disease severity. Four of these 23 loci likely have an ethnicity-specific component. Genome-wide association study (GWAS) signals in 11 loci colocalize with expression quantitative trait loci (eQTLs) associated with the expression of 20 genes in 62 tissues/cell types (range: 1:43 tissues/gene), including lung, brain, heart, muscle, and skin as well as the digestive system and immune system. We perform genetic fine mapping to compute 99% credible SNP sets, which identify 10 GWAS loci that have eight or fewer SNPs in the credible set, including three loci with one single likely causal SNP. Our study suggests that the diverse symptoms and disease severity of COVID-19 observed between individuals is associated with variants across the genome, affecting gene expression levels in a wide variety of tissue types., Graphical abstract, D’Antonio et al. characterize associations between GWAS signals for COVID-19 disease and eQTLs in 69 human tissues to identify causal variants and their underlying molecular mechanisms. They show that diverse symptoms and disease severity of COVID-19 are associated with variants affecting gene expression in a wide variety of tissues.

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P-8) interactions with sex on lung function, and 21 showed suggestive interactions (P-6). The strongest signal, from rs7697189 (chr4:145436894) on forced expiratory volume in 1 second (FEV1) (P=3.15x10-15), was replicated (P=0.016) in SpiroMeta. The C allele increased FEV1 more in males (untransformed FEV1 β=0.028 [SE 0.0022] litres) than females (β=0.009 [SE 0.0014] litres), and this effect was not accounted for by differential effects on height, smoking or pubertal age. rs7697189 resides upstream of the hedgehog-interacting protein (HHIP) gene and was previously associated with lung function and HHIP lung expression. We found HHIP expression was significantly different between the sexes (P=6.90x10-6), but we could not detect sex differential effects of rs7697189 on expression. Conclusions: We identified a novel genotype-by-sex interaction at a putative enhancer region upstream of the HHIP gene. Establishing the mechanism by which HHIP SNPs have different effects on lung function in males and females will be important for our understanding of lung health and diseases in both sexes.

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. Tweetable abstract @ERSpublications click to tweet Heterozygosity for the SERPINA1 Z allele is associated with higher lung function in nonsmokers (but this advantage is abolished by smoking) as well as a variety of nonrespiratory diseases and traits, including liver-, heart- and bone-related conditions https://bit.ly/37XC0Yg Introduction Homozygosity for the SERPINA1 Z allele (rs28929474(T)) is the commonest cause of severe α1-antitrypsin deficiency (AATD) and is a well-established genetic risk factor for lung diseases such as chronic obstructive pulmonary disease (COPD). However, the health consequences of heterozygosity for the Z allele are not as well-understood [1]. Given that approximately one in 30 Europeans is heterozygous for the Z allele, the phenotypic consequences of carriage of this allele could have important public health implications. Some previous studies have sought to characterise the effect of Z allele heterozygosity on nonrespiratory traits, particularly liver diseases [2–6]. However, these have often been carried out in small sample sizes and/or clinical subgroups. The recent development of phenome-wide association studies (PheWASs) and the availability of the well-phenotyped UK Biobank population-based cohort provides a platform for systematic investigation of the effects of heterozygosity for the Z allele on nonrespiratory traits. PheWASs test the association between genetic variants and a large number of phenotypic traits, including diseases and their subtypes, and potential intermediate phenotypes [7]. This differs from genome-wide association studies, which test a large number of variants across the genome for association with only one trait. The effect of Z allele heterozygosity on lung function traits and lung disease has been the subject of many studies. Recent COPD case–control and family-based studies have shown reduced lung function and increased risk of COPD in heterozygous current and former smokers [8–11]. However, a population-based study demonstrated no significant reductions in lung function in heterozygous smokers, despite having greater numbers of heterozygous smokers compared to previous studies [12]. It also showed enhanced lung function in heterozygous individuals overall, partially explained by strong association of the Z allele with increased height. This discrepancy may be due, in part, to the fact that identifying and recruiting study participants based on their health status (as in case–control studies [8, 9]) or based on the health status of a family member [10, 11] can lead to causal estimates that are subject to ascertainment bias [13], whereas population-based studies [12, 14] overcome these biases. The effects of the Z allele on lung function in relation to smoking status therefore remain uncertain. Finally, despite evidence for sex-differential effects of Z allele homozygosity on lung function [15, 16], we are not aware of any studies that have compared the effect of Z allele heterozygosity on lung function in males versus females. We therefore systematically evaluated the effects of Z allele heterozygosity in the UK biobank population, which in total includes >18 000 Z allele heterozygotes. We aimed: 1) to undertake the most extensive PheWAS to date, including blood biomarkers, for Z allele heterozygosity and homozygosity to identify effects beyond the respiratory system; and 2) to fully define the effects of Z allele heterozygosity on lung function measures (forced expiratory volume in 1 s (FEV1) and forced vital capacity (FVC)) in smokers and nonsmokers and in males and females.

Some individuals have characteristics of both asthma and COPD (asthma-COPD overlap, ACO), and evidence suggests they experience worse outcomes than those with either condition alone. Improved knowledge of the genetic architecture would contribute to understanding whether determinants of risk in this group differ from those in COPD or asthma.We conducted a genome-wide association study in 8,068 cases and 40,360 controls of European ancestry from UK Biobank (stage 1). After excluding variants only associated with asthma or COPD we selected 31 variants for further investigation in 12 additional cohorts (stage 2), and discovered eight novel signals for ACO in a meta-analysis of stage 1 and 2 studies.Our signals include an intergenic signal on chromosome 5 not previously associated with asthma, COPD or lung function, and suggest a spectrum of shared and distinct genetic influences in asthma, COPD and ACO. A number of signals may represent loci that predispose to serious long-term consequences in people with asthma.

Treatment-resistant depression (TRD) is a major contributor to the disability caused by major depressive disorder (MDD). Using primary care electronic health records from UK Biobank and EXCEED studies, we defined MDD and TRD, providing an easily accessible approach to investigate their clinical and genetic characteristics.MDD defined from primary care records was compared with other measures of depression and validated using the MDD polygenic risk score (PRS). Using prescribing records, TRD was defined from at least two switches between antidepressant drugs, each prescribed for at least six weeks. Clinical-demographic characteristics, SNP-heritability and genetic overlap with psychiatric and non-psychiatric traits were compared in TRD and non-TRD MDD cases.In 230,096 and 8,926 UKB and EXCEED participants with primary care data, respectively, the prevalence of MDD was 8.7% and 14.2%, of which 13.2% and 13.5% was TRD (2,430 and 159 cases), respectively. In both cohorts, MDD defined from primary care records was strongly associated with MDD PRS, and in UKB it showed overlap of 72%-88% with other MDD definitions. In UKB, TRD and non-TRD heritability was comparable (h2SNP = 0.25 [SE=0.04] and 0.19 [SE=0.02], respectively). TRD was positively associated with the polygenic risk score (PRS) of attention deficit hyperactivity disorder and negatively associated with the PRS of intelligence compared to non-TRD. It was more strongly associated with unfavourable clinical-demographic variables than non-TRD.This study demonstrated that MDD and TRD can be reliably defined using primary care records and provides the first large scale population assessment of the genetic, clinical and demographic characteristics of TRD.

BackgroundEosinophils are granulocytes associated with airway inflammation in respiratory disease. Eosinophil production and survival is controlled by interleukin-5: anti-interleukin-5 agents reduce asthma and COPD exacerbation frequency, and response correlates with baseline eosinophil counts. However, causal relationships between eosinophils and other respiratory phenotypes are less studied.MethodsWe investigated causality between eosinophils and: lung function, acute exacerbations of COPD (AECOPD), asthma-COPD overlap (ACO), moderate-to-severe asthma, and respiratory infections. We performed Mendelian randomization (MR) using 151 genetic variants from genome-wide association studies of blood eosinophil counts in UK Biobank/INTERVAL, and respiratory data from UK Biobank, using MR methods relying on different assumptions for validity. Multivariable MR using eight blood cell type exposures was performed for outcomes showing evidence of causation by eosinophils.FindingsThere was evidence that higher eosinophils reduce FEV1/FVC and FEV1 (weighted median estimator, SD change FEV1/FVC per SD eosinophils: −0.054 [95%CI −0.078,−0.029]. There was also evidence that eosinophils cause ACO (weighted median OR 1.44 [95%CI 1.19,1.74]), and asthma (weighted median OR 1.50 [95%CI 1.23,1.83]). Multivariable MR for FEV1/FVC, FEV1, ACO and asthma suggested that eosinophils were the cell type with the most important effect. Causal estimates of individual variants were heterogeneous, which may arise from pleiotropy.InterpretationWe found evidence that eosinophils reduce lung function, and increase ACO and asthma risk, on average over the set of genetic variants studied. Eosinophils appear to be causal determinants of fixed airflow obstruction among individuals with features of both asthma and COPD.FundingWellcome, BHF, MRC, BBSRC CASE studentship with GSK, GSK/BLF.

Homozygosity for theSERPINA1Z allele causes alpha-1 antitrypsin deficiency, a rare condition that can cause lung and liver disease. However, the effects of Z allele heterozygosity on non-respiratory phenotypes, and on lung function in the general population, remain unclear. We conducted the largest population-based study to date to determine Z allele effects on >2,400 phenotypes using the UK Biobank study (N>303,353). We detected strong associations between heterozygosity and non-respiratory phenotypes including increased height, 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 lung function in non-smokers, but smoking appears to abolish this protective effect. Individuals heterozygous for the Z allele may therefore have altered risk of smoking-induced lung disease and other, non-respiratory conditions.

Background: Lung function is highly heritable and differs between the sexes throughout life. However, little is known about sex-differential genetic effects on lung function. We aimed to conduct the first genome-wide genotype-by-sex interaction study on lung function to identify genetic effects that differ between males and females. Methods: We tested for interactions between 7,745,864 variants and sex on spirometry-based measures of lung function in UK Biobank (N=303,612), and sought replication in 75,696 independent individuals from the SpiroMeta consortium. Results: Five independent single-nucleotide polymorphisms (SNPs) showed genome-wide significant (P

The IL1RL1 (ST2) gene locus is robustly associated with asthma; however, the contribution of single nucleotide polymorphisms (SNPs) in this locus to specific asthma subtypes and the functional mechanisms underlying these associations remain to be defined. We tested for association between IL1RL1 region SNPs and characteristics of asthma as defined by clinical and immunological measures and addressed functional effects of these genetic variants in lung tissue and airway epithelium. Utilizing 4 independent cohorts (Lifelines, Dutch Asthma GWAS [DAG], Genetics of Asthma Severity and Phenotypes [GASP], and Manchester Asthma and Allergy Study [MAAS]) and resequencing data, we identified 3 key signals associated with asthma features. Investigations in lung tissue and primary bronchial epithelial cells identified context-dependent relationships between the signals and IL1RL1 mRNA and soluble protein expression. This was also observed for asthma-associated IL1RL1 nonsynonymous coding TIR domain SNPs. Bronchial epithelial cell cultures from asthma patients, exposed to exacerbation-relevant stimulations, revealed modulatory effects for all 4 signals on IL1RL1 mRNA and/ or protein expression, suggesting SNP-environment interactions. The IL1RL1 TIR signaling domain haplotype affected IL-33-driven NF-κB signaling, while not interfering with TLR signaling. In summary, we identify that IL1RL1 genetic signals potentially contribute to severe and eosinophilic phenotypes in asthma, as well as provide initial mechanistic insight, including genetic regulation of IL1RL1 isoform expression and receptor signaling.

Chronic obstructive pulmonary disease (COPD), diagnosed by reduced lung function, is a leading cause of morbidity and mortality. We performed whole genome sequence (WGS) analysis of lung function and COPD in a multi-ethnic sample of 11,497 participants from population- and family-based studies, and 8499 individuals from COPD-enriched studies in the NHLBI Trans-Omics for Precision Medicine (TOPMed) Program. We identify at genome-wide significance 10 known GWAS loci and 22 distinct, previously unreported loci, including two common variant signals from stratified analysis of African Americans. Four novel common variants within the regions of PIAS1, RGN (two variants) and FTO show evidence of replication in the UK Biobank (European ancestry n ~ 320,000), while colocalization analyses leveraging multi-omic data from GTEx and TOPMed identify potential molecular mechanisms underlying four of the 22 novel loci. Our study demonstrates the value of performing WGS analyses and multi-omic follow-up in cohorts of diverse ancestry., Chronic obstructive pulmonary disease is a leading cause of morbidity and mortality. Here, the authors analyse whole genome sequence data and find new loci associated with lung function and COPD.

This cohort study assesses the interaction of cigarette smoking and polygenic risk score in association with reduced lung function among participants in the UK Biobank cohort., Key Points Question Does cigarette smoking interact with genetic risk on the percent of forced vital capacity exhaled in the first second (FEV1/FVC)? Findings In this UK Biobank cohort study of 319 730 UK citizens, FEV1/FVC was associated with polygenic risk score-by-smoking interactions, and smoking was detrimental across all categories of estimated genetic risk, although it was worse for those with the highest estimated genetic risks. For every reported 20 pack-years of smoking, individuals in the top decile compared with the bottom decile of genetic risk showed nearly twice the reduction in FEV1/FVC. Meaning These findings suggest that elucidating mechanisms for the interaction between smoking and genetic risk could yield greater insight into the chronic obstructive pulmonary disease pathogenesis., Importance The risk of airflow limitation and chronic obstructive pulmonary disease (COPD) is influenced by combinations of cigarette smoking and genetic susceptibility, yet it remains unclear whether gene-by-smoking interactions are associated with quantitative measures of lung function. Objective To assess the interaction of cigarette smoking and polygenic risk score in association with reduced lung function. Design, Setting, and Participants This UK Biobank cohort study included UK citizens of European ancestry aged 40 to 69 years with genetic and spirometry data passing quality control metrics. Data was analyzed from July 2020 to March 2021. Exposures PRS of combined forced expiratory volume in 1 second (FEV1) and percent of forced vital capacity exhaled in the first second (FEV1/FVC), self-reported pack-years of smoking, ever- vs never-smoking status, and current- vs former- or never-smoking status. Main Outcomes and Measures FEV1/FVC was the primary outcome. Models were used to test for interactions with models, including the main effects of PRS, different smoking variables, and their cross-product terms. The association between pack-years of smoking and FEV1/FVC were compared for those in the highest vs lowest decile of estimated genetic risk for low lung function. Results We included 319 730 individuals, of whom 24 915 (8%) had moderate-to-severe COPD cases, and 44.4% were men. Participants had a mean (SD) age 56.5 of (8.02) years. The PRS and pack-years were significantly associated with lower FEV1/FVC (PRS: β, −0.03; 95% CI, −0.031 to −0.03; pack-years: β, −0.0064; 95% CI, −0.0064 to −0.0063) and the interaction term (β, −0.0028; 95% CI, −0.0029 to −0.0026). A stepwise increment in estimated effect sizes for these interaction terms was observed per 10 pack-years of smoking exposure. The interaction of PRS with 11 to 20, 31 to 40, and more than 50 pack-years categories were β (interaction) −0.0038 (95% CI, −0.0046 to −0.0031); −0.013 (95% CI, −0.014 to −0.012); and −0.017 (95% CI, −0.019 to −0.016), respectively. There was evidence of significant interaction between PRS with ever- or never- smoking status (β, interaction; −0.0064; 95% CI, −0.0068 to −0.0060) and current or not-current smoking (β, interaction; −0.0091; 95% CI, −0.0097 to −0.0084). For any given level of pack-years of smoking exposure, FEV1/FVC was significantly lower for individuals in the tenth decile (ie, highest risk) than the first decile (ie, lowest risk) of genetic risk. For every 20 pack-years of smoking, those in the tenth decile compared with the first decile of genetic risk showed nearly a 2-fold reduction in FEV1/FVC. Conclusions and Relevance COPD is characterized by diminished lung function, and our analyses suggest there is substantial interaction between genome-wide PRS and smoking exposures. While smoking was associated with decreased lung function across all genetic risk categories, the associations were strongest in individuals with higher estimated genetic risk.

Background: Globally, respiratory infections contribute to significant morbidity and mortality. However, genetic determinants of respiratory infections are understudied and remain poorly understood. Methods: We conducted a genome-wide association study in 19,459 hospitalised respiratory infection cases and 101,438 controls from UK Biobank. We followed-up well-imputed top signals from the UK Biobank discovery analysis in 50,912 respiratory infection cases and 150,442 controls from 11 cohorts. We aggregated effect estimates across studies using inverse variance-weighted meta-analyses. Additionally, we investigated the function of the top signals in order to gain understanding of the underlying biological mechanisms. Results: In the discovery analysis, we report 56 signals at P×10-6, one of which was genome-wide significant (P×10-8). The genome-wide significant signal was in an intron of PBX3, a gene that encodes pre-B-cell leukaemia transcription factor 3, a homeodomain-containing transcription factor. Further, the genome-wide significant signal was found to colocalise with gene-specific expression quantitative trait loci (eQTLs) affecting expression of PBX3 in lung tissue, where the respiratory infection risk alleles were associated with decreased PBX3 expression in lung tissue, highlighting a possible biological mechanism. Of the 56 signals, 40 were well-imputed in UK Biobank and were investigated in the 11 follow-up cohorts. None of the 40 signals replicated, with effect estimates attenuated. Conclusions: Our discovery analysis implicated PBX3 as a candidate causal gene and suggests a possible role of transcription factor binding activity in respiratory infection susceptibility. However, the PBX3 signal, and the other well-imputed signals, did not replicate when aggregating effect estimates across 11 independent cohorts. Significant phenotypic heterogeneity and differences in study ascertainment may have contributed to this lack of statistical replication. Overall, our study highlighted putative associations and possible biological mechanisms that may provide insight into respiratory infection susceptibility.

A trend towards earlier menarche in women has been associated with childhood factors (e.g. obesity) and hypothesised environmental exposures (e.g. endocrine disruptors present in household products). Observational evidence has shown detrimental effects of early menarche on various health outcomes including adult lung function, but these might represent spurious associations due to confounding. To address this we used Mendelian randomization where genetic variants are used as proxies for age at menarche, since genetic associations are not affected by classical confounding. We estimated the effects of age at menarche on forced vital capacity (FVC), a proxy for restrictive lung impairment, and ratio of forced expiratory volume in one second to FVC (FEV1/FVC), a measure of airway obstruction, in both adulthood and adolescence. We derived SNP-age at menarche association estimates for 122 variants from a published genome-wide meta-analysis (N = 182,416), with SNP-lung function estimates obtained by meta-analysing three studies of adult women (N = 46,944) and two of adolescent girls (N = 3025). We investigated the impact of departures from the assumption of no pleiotropy through sensitivity analyses. In adult women, in line with previous evidence, we found an effect on restrictive lung impairment with a 24.8 mL increase in FVC per year increase in age at menarche (95% CI 1.8–47.9; p = 0.035); evidence was stronger after excluding potential pleiotropic variants (43.6 mL; 17.2–69.9; p = 0.001). In adolescent girls we found an opposite effect (−56.5 mL; −108.3 to −4.7; p = 0.033), suggesting that the detrimental effect in adulthood may be preceded by a short-term post-pubertal benefit. Our secondary analyses showing results in the same direction in men and boys, in whom age at menarche SNPs have also shown association with sexual development, suggest a role for pubertal timing in general rather than menarche specifically. We found no effect on airway obstruction (FEV1/FVC). Electronic supplementary material The online version of this article (doi:10.1007/s10654-017-0272-9) contains supplementary material, which is available to authorized users.

Chronic obstructive pulmonary disease (COPD) is a leading cause of mortality worldwide(1). We performed a genetic association study in 15,256 cases and 47,936 controls, with replication of select top results (P

BACKGROUND: Genome-wide association studies have identified numerous genetic regions that influence cross-sectional lung function. Longitudinal decline in lung function also includes a heritable component but the genetic determinants have yet to be defined.OBJECTIVES: We aimed to determine whether regions associated with cross-sectional lung function were also associated with longitudinal decline and to seek novel variants which influence decline.METHODS: We analysed genome-wide data from 4167 individuals from the Busselton Health Study cohort, who had undergone spirometry (12 695 observations across eight time points). A mixed model was fitted and weighted risk scores were calculated for the joint effect of 26 known regions on baseline and longitudinal changes in FEV1 and FEV1/FVC. Potential additional regions of interest were identified and followed up in two independent cohorts.RESULTS: The 26 regions previously associated with cross-sectional lung function jointly showed a strong effect on baseline lung function (p=4.44×10-16 for FEV1/FVC) but no effect on longitudinal decline (p=0.160 for FEV1/FVC). This was replicated in an independent cohort. 39 additional regions of interest (48 variants) were identified; these associations were not replicated in two further cohorts.CONCLUSIONS: Previously identified genetic variants jointly have a strong effect on cross-sectional lung function in adults but little or no effect on the rate of decline of lung function. It is possible that they influence COPD risk through lung development. Although no genetic variants have yet been associated with lung function decline at stringent genome-wide significance, longitudinal change in lung function is heritable suggesting that there is scope for future discoveries.

Chronic obstructive pulmonary disease (COPD) is characterized by reduced lung function and is the third leading cause of death globally. Through genome-wide association discovery in 48,943 individuals, selected from extremes of the lung function distribution in UK Biobank, and follow-up in 95,375 individuals, we increased the yield of independent signals for lung function from 54 to 97. A genetic risk score was associated with COPD susceptibility (odds ratio per 1 s.d. of the risk score (similar to 6 alleles) (95% confidence interval) = 1.24 (1.20-1.27), P = 5.05 x 10(-49)), and we observed a 3.7-fold difference in COPD risk between individuals in the highest and lowest genetic risk score deciles in UK Biobank. The 97 signals show enrichment in genes for development, elastic fibers and epigenetic regulation pathways. We highlight targets for drugs and compounds in development for COPD and asthma (genes in the inositol phosphate metabolism pathway and CHRM3) and describe targets for potential drug repositioning from other clinical indications.

Adult lung function is highly heritable and 279 genetic loci were recently reported as associated with spirometry-based measures of lung function. Though lung development and function differ between males and females throughout life, there has been no genome-wide study to identify genetic variants with differential effects on lung function in males and females. Here, we present the first genome-wide genotype-by-sex interaction study on four lung function traits in 303,612 participants from the UK Biobank. We detected five SNPs showing genome-wide significant (P−8) interactions with sex on lung function, as well as 21 suggestively significant interactions (P−6). The strongest sex interaction signal came from rs7697189 at 4:145436894 on forced expiratory volume in 1 second (FEV1) (P = 3.15 × 10−15), and was replicated (P = 0.016) in 75,696 individuals in the SpiroMeta consortium. Sex-stratified analyses demonstrated that the minor (C) allele of rs7697189 increased lung function to a greater extent in males than females (untransformed FEV1β = 0.028 [SE 0.0022] litres in males vs β = 0.009 [SE 0.0014] litres in females), and this effect was not accounted for by differential effects on height, smoking or age at puberty. This SNP resides upstream of the gene encoding hedgehog-interacting protein (HHIP) and has previously been reported for association with lung function andHHIPexpression in lung tissue. In our analyses, whileHHIPexpression in lung tissue was significantly different between the sexes with females having higher expression (most significant probeset P=6.90 × 10−6) after adjusting for age and smoking, rs7697189 did not demonstrate sex differential effects on expression. Establishing the mechanism by whichHHIPSNPs have different effects on lung function in males and females will be important for our understanding of lung health and diseases, such as chronic obstructive pulmonary disease (COPD), in both sexes.

The UK Biobank (UKB) has recently released genotypes on 152,328 individuals together with extensive phenotypic and lifestyle information. We present a new phasing method, SHAPEIT3, that can handle such biobank-scale data sets and results in switch error rates as low as ∼0.3%. The method exhibits O(NlogN) scaling with sample size N, enabling fast and accurate phasing of even larger cohorts.

EXCEED is a longitudinal population-based cohort which facilitates investigation of genetic, environmental and lifestyle-related determinants of a broad range of diseases and of multiple morbidity through data collected at baseline and via electronic healthcare record linkage. Recruitment has taken place in Leicester, Leicestershire and Rutland since 2013 and is ongoing, with 9,840 participants aged 30-69 to date. The population of Leicester is diverse and additional recruitment from the local South Asian community is ongoing. Participants have consented to follow-up for up to 25 years through electronic health records and additional bespoke data collection is planned. Data available includes baseline demographics, anthropometry, spirometry, lifestyle factors (smoking and alcohol use) and longitudinal health information from primary care records, with additional linkage to other EHR datasets planned. Patients have consented to be contacted for recall-by-genotype and recall-by-phenotype sub-studies, providing an important resource for precision medicine research. We welcome requests for collaboration and data access by contacting the study management team via exceed@le.ac.uk.