Genome-wide analysis in over 1 million individuals of European ancestry yields improved polygenic risk scores for blood pressure traits.

Hypertension affects more than one billion people worldwide. Here we identify 113 novel loci, reporting a total of 2,103 independent genetic signals (P < 5 × 10-8) from the largest single-stage blood pressure (BP) genome-wide association study to date (n = 1,028,980 European individuals). These associations explain more than 60% of single nucleotide polymorphism-based BP heritability. Comparing top versus bottom deciles of polygenic risk scores (PRSs) reveals clinically meaningful differences in BP (16.9 mmHg systolic BP, 95% CI, 15.5-18.2 mmHg, P = 2.22 × 10-126) and more than a sevenfold higher odds of hypertension risk (odds ratio, 7.33; 95% CI, 5.54-9.70; P = 4.13 × 10-44) in an independent dataset. Adding PRS into hypertension-prediction models increased the area under the receiver operating characteristic curve (AUROC) from 0.791 (95% CI, 0.781-0.801) to 0.826 (95% CI, 0.817-0.836, ∆AUROC, 0.035, P = 1.98 × 10-34). We compare the 2,103 loci results in non-European ancestries and show significant PRS associations in a large African-American sample. Secondary analyses implicate 500 genes previously unreported for BP. Our study highlights the role of increasingly large genomic studies for precision health research.

Air pollution, lung function and mortality: survival and mediation analyses in UK Biobank.

Background: Air pollution is associated with lower lung function, and both are associated with premature mortality and cardiovascular disease (CVD). Evidence remains scarce on the potential mediating effect of impaired lung function on the association between air pollution and mortality or CVD. Methods: We used data from UK Biobank (n∼200 000 individuals) with 8-year follow-up to mortality and incident CVD. Exposures to particulate matter <10 µm (PM10), particulate matter <2.5 µm (PM2.5) and nitrogen dioxide (NO2) were assessed by land-use regression modelling. Lung function (forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and the FEV1/FVC ratio) was measured between 2006 and 2010 and transformed to Global Lung Function Initiative (GLI) z-scores. Adjusted Cox proportional hazards and causal proportional hazards mediation analysis models were fitted, stratified by smoking status. Results: Lower FEV1 and FVC were associated with all-cause and CVD mortality, and incident CVD, with larger estimates in ever- than never-smokers (all-cause mortality hazard ratio per FEV1 GLI z-score decrease 1.29 (95% CI 1.24-1.34) for ever-smokers and 1.16 (95% CI 1.12-1.21) for never-smokers). Long-term exposure to PM2.5 or NO2 was associated with incident CVD, with similar effect sizes for ever- and never-smokers. Mediated proportions of the air pollution-all-cause mortality estimates driven by FEV1 were 18% (95% CI 2-33%) for PM2.5 and 27% (95% CI 3-51%) for NO2. Corresponding mediated proportions for incident CVD were 9% (95% CI 4-13%) for PM2.5 and 16% (95% CI 6-25%) for NO2. Conclusions: Lung function may mediate a modest proportion of associations between air pollution and mortality and CVD outcomes. Results likely reflect the extent of either shared mechanisms or direct effects relating to lower lung function caused by air pollution.

Association study of human leukocyte antigen variants and idiopathic pulmonary fibrosis.

Introduction: Idiopathic pulmonary fibrosis (IPF) is a chronic interstitial pneumonia marked by progressive lung fibrosis and a poor prognosis. Recent studies have highlighted the potential role of infection in the pathogenesis of IPF, and a prior association of the HLA-DQB1 gene with idiopathic fibrotic interstitial pneumonia (including IPF) has been reported. Owing to the important role that the human leukocyte antigen (HLA) region plays in the immune response, here we evaluated if HLA genetic variation was associated specifically with IPF risk. Methods: We performed a meta-analysis of associations of the HLA region with IPF risk in individuals of European ancestry from seven independent case-control studies of IPF (comprising 5159 cases and 27 459 controls, including a prior study of fibrotic interstitial pneumonia). Single nucleotide polymorphisms, classical HLA alleles and amino acids were analysed and signals meeting a region-wide association threshold of p<4.5×10-4 and a posterior probability of replication >90% were considered significant. We sought to replicate the previously reported HLA-DQB1 association in the subset of studies independent of the original report. Results: The meta-analysis of all seven studies identified four significant independent single nucleotide polymorphisms associated with IPF risk. However, none met the posterior probability for replication criterion. The HLA-DQB1 association was not replicated in the independent IPF studies. Conclusion: Variation in the HLA region was not consistently associated with risk in studies of IPF. However, this does not preclude the possibility that other genomic regions linked to the immune response may be involved in the aetiology of IPF.

Asthma exacerbations and eosinophilia in the UK Biobank: a genome-wide association study.

Background: Asthma exacerbations reflect disease severity, affect morbidity and mortality, and may lead to declining lung function. Inflammatory endotypes (e.g. T2-high (eosinophilic)) may play a key role in asthma exacerbations. We aimed to assess whether genetic susceptibility underlies asthma exacerbation risk and additionally tested for an interaction between genetic variants and eosinophilia on exacerbation risk. Methods: UK Biobank data were used to perform a genome-wide association study of individuals with asthma and at least one exacerbation compared to individuals with asthma and no history of exacerbations. Individuals with asthma were identified using self-reported data, hospitalisation data and general practitioner records. Exacerbations were identified as either asthma-related hospitalisation, general practitioner record of asthma exacerbation or an oral corticosteroid burst prescription. A logistic regression model adjusted for age, sex, smoking status and genetic ancestry via principal components was used to assess the association between genetic variants and asthma exacerbations. We sought replication for suggestive associations (p<5×10-6) in the GERA cohort. Results: In the UK Biobank, we identified 11 604 cases and 37 890 controls. While no variants reached genome-wide significance (p<5×10-8) in the primary analysis, 116 signals were suggestively significant (p<5×10-6). In GERA, two single nucleotide polymorphisms (rs34643691 and rs149721630) replicated (p<0.05), representing signals near the NTRK3 and ABCA13 genes. Conclusions: Our study has identified reproducible associations with asthma exacerbations in the UK Biobank and GERA cohorts. Confirmation of these findings in different asthma subphenotypes in diverse ancestries and functional investigation will be required to understand their mechanisms of action and potentially inform therapeutic development.

Effect of Atypical Sleep EEG Patterns on Weaning From Prolonged Mechanical Ventilation.

BACKGROUND: Approximately one-third of acute ICU patients display atypical sleep patterns that cannot be interpreted by using standard EEG criteria for sleep. Atypical sleep patterns have been associated with poor weaning outcomes in acute ICUs. RESEARCH QUESTION: Do patients being weaned from prolonged mechanical ventilation experience atypical sleep EEG patterns, and are these patterns linked with weaning outcomes? STUDY DESIGN AND METHODS: EEG power spectral analysis during wakefulness and overnight polysomnogram were performed on alert, nondelirious patients at a long-term acute care facility. RESULTS: Forty-four patients had been ventilated for a median duration of 38 days at the time of the polysomnogram study. Eleven patients (25%) exhibited atypical sleep EEG. During wakefulness, relative EEG power spectral analysis revealed higher relative delta power in patients with atypical sleep than in patients with usual sleep (53% vs 41%; P < .001) and a higher slow-to-fast power ratio during wakefulness: 4.39 vs 2.17 (P < .001). Patients with atypical sleep displayed more subsyndromal delirium (36% vs 6%; P = .027) and less rapid eye movement sleep (4% vs 11% total sleep time; P < .02). Weaning failure was more common in the atypical sleep group than in the usual sleep group: 91% vs 45% (P = .013). INTERPRETATION: This study provides the first evidence that patients in a long-term acute care facility being weaned from prolonged ventilation exhibit atypical sleep EEG patterns that are associated with weaning failure. Patients with atypical sleep EEG patterns had higher rates of subsyndromal delirium and slowing of the wakeful EEG, suggesting that these two findings represent a biological signal for brain dysfunction.

PCSK9 genetic variants and risk of vascular and non-vascular diseases in Chinese and UK populations.

AIMS: Lowering low-density lipoprotein cholesterol (LDL-C) through PCSK9 inhibition represents a new therapeutic approach to preventing and treating cardiovascular disease (CVD). Phenome-wide analyses of PCSK9 genetic variants in large biobanks can help to identify unexpected effects of PCSK9 inhibition. METHODS AND RESULTS: In the prospective China Kadoorie Biobank, we constructed a genetic score using three variants at the PCSK9 locus associated with directly measured LDL-C [PCSK9 genetic score (PCSK9-GS)]. Logistic regression gave estimated odds ratios (ORs) for PCSK9-GS associations with CVD and non-CVD outcomes, scaled to 1 SD lower LDL-C. PCSK9-GS was associated with lower risks of carotid plaque [n = 8340 cases; OR = 0.61 (95% confidence interval: 0.45-0.83); P = 0.0015], major occlusive vascular events [n = 15 752; 0.80 (0.67-0.95); P = 0.011], and ischaemic stroke [n = 11 467; 0.80 (0.66-0.98); P = 0.029]. However, PCSK9-GS was also associated with higher risk of hospitalization with chronic obstructive pulmonary disease [COPD: n = 6836; 1.38 (1.08-1.76); P = 0.0089] and with even higher risk of fatal exacerbations amongst individuals with pre-existing COPD [n = 730; 3.61 (1.71-7.60); P = 7.3 × 10-4]. We also replicated associations for a PCSK9 variant, reported in UK Biobank, with increased risks of acute upper respiratory tract infection (URTI) [pooled OR after meta-analysis of 1.87 (1.38-2.54); P = 5.4 × 10-5] and self-reported asthma [pooled OR of 1.17 (1.04-1.30); P = 0.0071]. There was no association of a polygenic LDL-C score with COPD hospitalization, COPD exacerbation, or URTI. CONCLUSION: The LDL-C-lowering PCSK9 genetic variants are associated with lower risk of subclinical and clinical atherosclerotic vascular disease but higher risks of respiratory diseases. Pharmacovigilance studies may be required to monitor patients treated with therapeutic PCSK9 inhibitors for exacerbations of respiratory diseases or respiratory tract infections. LAY SUMMARY: Genetic analyses of over 100 000 participants of the China Kadoorie Biobank, mimicking the effect of new drugs intended to reduce cholesterol by targeting the PCSK9 protein, have identified potential severe effects of lower PCSK9 activity in patients with existing respiratory disease.PCSK9 genetic variants that are associated with lower cholesterol and reduced rates of cardiovascular disease are also associated with increased risk of a range of respiratory diseases, including asthma, upper respiratory tract infections, and hospitalization with chronic obstructive pulmonary disease (COPD).These genetic variants are not associated with whether or not individuals have COPD; instead, they are specifically associated with an increase in the chance of those who already have COPD being hospitalized and even dying, suggesting that careful monitoring of such patients should be considered during development of and treatment with anti-PCSK9 medication.

Genome-wide SNP-sex interaction analysis of susceptibility to idiopathic pulmonary fibrosis.

Background: Idiopathic pulmonary fibrosis (IPF) is a chronic lung condition that is more prevalent in males than females. The reasons for this are not fully understood, with differing environmental exposures due to historically sex-biased occupations, or diagnostic bias, being possible explanations. To date, over 20 independent genetic variants have been identified to be associated with IPF susceptibility, but these have been discovered when combining males and females. Our aim was to test for the presence of sex-specific associations with IPF susceptibility and assess whether there is a need to consider sex-specific effects when evaluating genetic risk in clinical prediction models for IPF. Methods: We performed genome-wide single nucleotide polymorphism (SNP)-by-sex interaction studies of IPF risk in six independent IPF case-control studies and combined them using inverse-variance weighted fixed effect meta-analysis. In total, 4,561 cases (1,280 females and 2,281 males) and 23,500 controls (8,360 females and 14,528 males) of European genetic ancestry were analysed. We used polygenic risk scores (PRS) to assess differences in genetic risk prediction between males and females. Findings: Three independent genetic association signals were identified. All showed a consistent direction of effect across all individual IPF studies and an opposite direction of effect in IPF susceptibility between females and males. None had been previously identified in IPF susceptibility genome-wide association studies (GWAS). The predictive accuracy of the PRSs were similar between males and females, regardless of whether using combined or sex-specific GWAS results. Interpretation: We prioritised three genetic variants whose effect on IPF risk may be modified by sex, however these require further study. We found no evidence that the predictive accuracy of common SNP-based PRSs varies significantly between males and females.

Genome-wide association study of preserved ratio impaired spirometry (PRISm).

BACKGROUND: Preserved ratio impaired spirometry (PRISm) is defined as a forced expiratory volume in 1 s (FEV1) <80% predicted and FEV1/forced vital capacity ≥0.70. PRISm is associated with respiratory symptoms and comorbidities. Our objective was to discover novel genetic signals for PRISm and see if they provide insight into the pathogenesis of PRISm and associated comorbidities. METHODS: We undertook a genome-wide association study (GWAS) of PRISm in UK Biobank participants (Stage 1), and selected single nucleotide polymorphisms (SNPs) reaching genome-wide significance for replication in 13 cohorts (Stage 2). A combined meta-analysis of Stage 1 and Stage 2 was done to determine top SNPs. We used cross-trait linkage disequilibrium score regression to estimate genome-wide genetic correlation between PRISm and pulmonary and extrapulmonary traits. Phenome-wide association studies of top SNPs were performed. RESULTS: 22 signals reached significance in the joint meta-analysis, including four signals novel for lung function. A strong genome-wide genetic correlation (rg) between PRISm and spirometric COPD (rg=0.62, p<0.001) was observed, and genetic correlation with type 2 diabetes (rg=0.12, p=0.007). Phenome-wide association studies showed that 18 of 22 signals were associated with diabetic traits and seven with blood pressure traits. CONCLUSION: This is the first GWAS to successfully identify SNPs associated with PRISm. Four of the signals, rs7652391 (nearest gene MECOM), rs9431040 (HLX), rs62018863 (TMEM114) and rs185937162 (HLA-B), have not been described in association with lung function before, demonstrating the utility of using different lung function phenotypes in GWAS. Genetic factors associated with PRISm are strongly correlated with risk of both other lung diseases and extrapulmonary comorbidity.

Genome-wide association study of thyroid-stimulating hormone highlights new genes, pathways and associations with thyroid disease.

Thyroid hormones play a critical role in regulation of multiple physiological functions and thyroid dysfunction is associated with substantial morbidity. Here, we use electronic health records to undertake a genome-wide association study of thyroid-stimulating hormone (TSH) levels, with a total sample size of 247,107. We identify 158 novel genetic associations, more than doubling the number of known associations with TSH, and implicate 112 putative causal genes, of which 76 are not previously implicated. A polygenic score for TSH is associated with TSH levels in African, South Asian, East Asian, Middle Eastern and admixed American ancestries, and associated with hypothyroidism and other thyroid disease in South Asians. In Europeans, the TSH polygenic score is associated with thyroid disease, including thyroid cancer and age-of-onset of hypothyroidism and hyperthyroidism. We develop pathway-specific genetic risk scores for TSH levels and use these in phenome-wide association studies to identify potential consequences of pathway perturbation. Together, these findings demonstrate the potential utility of genetic associations to inform future therapeutics and risk prediction for thyroid diseases.

Association study of human leukocyte antigen (HLA) variants and idiopathic pulmonary fibrosis.

Introduction: Idiopathic pulmonary fibrosis (IPF) is a chronic interstitial pneumonia marked by progressive lung fibrosis and a poor prognosis. Recent studies have highlighted the potential role of infection in the pathogenesis of IPF and a prior association of the HLA-DQB1 gene with idiopathic fibrotic interstitial pneumonia (including IPF) has been reported. Due to the important role that the Human Leukocyte Antigen (HLA) region plays in the immune response, here we evaluated if HLA genetic variation was associated specifically with IPF risk. Methods: We performed a meta-analysis of associations of the HLA region with IPF risk in individuals of European ancestry from seven independent case-control studies of IPF (comprising a total of 5,159 cases and 27,459 controls, including the prior study of fibrotic interstitial pneumonia). Single nucleotide polymorphisms, classical HLA alleles and amino acids were analysed and signals meeting a region-wide association threshold p<4.5×10-4 and a posterior probability of replication >90% were considered significant. We sought to replicate the previously reported HLA-DQB1 association in the subset of studies independent of the original report. Results: The meta-analysis of all seven studies identified four significant independent single nucleotide polymorphisms associated with IPF risk. However, none met the posterior probability for replication criterion. The HLA-DQB1 association was not replicated in the independent IPF studies. Conclusion: Variation in the HLA region was not consistently associated with risk in studies of IPF. However, this does not preclude the possibility that other genomic regions linked to the immune response may be involved in the aetiology of IPF.

Identification and characterisation of a rare MTTP variant underlying hereditary non-alcoholic fatty liver disease.

Background & Aims: Non-alcoholic fatty liver disease (NAFLD) is a complex trait with an estimated prevalence of 25% globally. 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. Methods: Exome 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 with 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 and restored to wild-type using CRISPR-Cas9. The phenotype was assessed using imaging, targeted RNA analysis, and molecular expression arrays. Results: We identified a rare causal variant c.1691T>C p.I564T (rs745447480) in MTTP, encoding microsomal triglyceride transfer protein (MTP), associated with progressive NAFLD, unrelated to metabolic syndrome and without characteristic features of abetalipoproteinaemia. HLCs derived from a homozygote donor had significantly lower MTP activity and lower lipoprotein ApoB secretion than wild-type cells, while having similar levels of MTP mRNA and protein. Cytoplasmic triglyceride accumulation in HLCs triggered endoplasmic reticulum stress, secretion of pro-inflammatory mediators, and production of reactive oxygen species. Conclusions: We have identified and characterised a rare causal variant in MTTP, and homozygosity for MTTP p.I564T is associated with progressive NAFLD without any other manifestations of abetalipoproteinaemia. Our findings provide insights into mechanisms driving progressive NAFLD. Impact and Implications: A rare genetic variant in the gene MTTP has been identified as responsible for the development of severe non-alcoholic fatty liver disease in a four-generation family with no typical disease risk factors. A cell line culture created harbouring this variant gene was characterised to understand how this genetic variation leads to a defect in liver cells, which results in accumulation of fat and processes that promote disease. This is now a useful model for studying the disease pathways and to discover new ways to treat common types of fatty liver disease.

A Polygenic Risk Score for Idiopathic Pulmonary Fibrosis and Interstitial Lung Abnormalities.

Rationale: In addition to rare genetic variants and the MUC5B locus, common genetic variants contribute to idiopathic pulmonary fibrosis (IPF) risk. The predictive power of common variants outside the MUC5B locus for IPF and interstitial lung abnormalities (ILAs) is unknown. Objectives: We tested the predictive value of IPF polygenic risk scores (PRSs) with and without the MUC5B region on IPF, ILA, and ILA progression. Methods: We developed PRSs that included (PRS-M5B) and excluded (PRS-NO-M5B) the MUC5B region (500-kb window around rs35705950-T) using an IPF genome-wide association study. We assessed PRS associations with area under the receiver operating characteristic curve (AUC) metrics for IPF, ILA, and ILA progression. Measurements and Main Results: We included 14,650 participants (1,970 IPF; 1,068 ILA) from six multi-ancestry population-based and case-control cohorts. In cases excluded from genome-wide association study, the PRS-M5B (odds ratio [OR] per SD of the score, 3.1; P = 7.1 × 10-95) and PRS-NO-M5B (OR per SD, 2.8; P = 2.5 × 10-87) were associated with IPF. Participants in the top PRS-NO-M5B quintile had ∼sevenfold odds for IPF compared with those in the first quintile. A clinical model predicted IPF (AUC, 0.61); rs35705950-T and PRS-NO-M5B demonstrated higher AUCs (0.73 and 0.7, respectively), and adding both genetic predictors to a clinical model yielded the highest performance (AUC, 0.81). The PRS-NO-M5B was associated with ILA (OR, 1.25) and ILA progression (OR, 1.16) in European ancestry participants. Conclusions: A common genetic variant risk score complements the MUC5B variant to identify individuals at high risk of interstitial lung abnormalities and pulmonary fibrosis.

Genome-wide association study of chronic sputum production implicates loci involved in mucus production and infection.

BACKGROUND: Chronic sputum production impacts on quality of life and is a feature of many respiratory diseases. Identification of the genetic variants associated with chronic sputum production in a disease agnostic sample could improve understanding of its causes and identify new molecular targets for treatment. METHODS: We conducted a genome-wide association study (GWAS) of chronic sputum production in UK Biobank. Signals meeting genome-wide significance (p<5×10-8) were investigated in additional independent studies, were fine-mapped and putative causal genes identified by gene expression analysis. GWASs of respiratory traits were interrogated to identify whether the signals were driven by existing respiratory disease among the cases and variants were further investigated for wider pleiotropic effects using phenome-wide association studies (PheWASs). RESULTS: From a GWAS of 9714 cases and 48 471 controls, we identified six novel genome-wide significant signals for chronic sputum production including signals in the human leukocyte antigen (HLA) locus, chromosome 11 mucin locus (containing MUC2, MUC5AC and MUC5B) and FUT2 locus. The four common variant associations were supported by independent studies with a combined sample size of up to 2203 cases and 17 627 controls. The mucin locus signal had previously been reported for association with moderate-to-severe asthma. The HLA signal was fine-mapped to an amino acid change of threonine to arginine (frequency 36.8%) in HLA-DRB1 (HLA-DRB1*03:147). The signal near FUT2 was associated with expression of several genes including FUT2, for which the direction of effect was tissue dependent. Our PheWAS identified a wide range of associations including blood cell traits, liver biomarkers, infections, gastrointestinal and thyroid-associated diseases, and respiratory disease. CONCLUSIONS: Novel signals at the FUT2 and mucin loci suggest that mucin fucosylation may be a driver of chronic sputum production even in the absence of diagnosed respiratory disease and provide genetic support for this pathway as a target for therapeutic intervention.

Identifying COPD subtypes using multi-trait genetics.

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.

DeepPheWAS: an R package for phenotype generation and association analysis for phenome-wide association studies.

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.

PCSK6 and Survival in Idiopathic Pulmonary Fibrosis.

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.

Mendelian randomisation of eosinophils and other cell types in relation to lung function and disease.

RATIONALE: 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. OBJECTIVES: We investigated causality between eosinophils and: lung function, acute exacerbations of COPD, asthma-COPD overlap (ACO), moderate-to-severe asthma and respiratory infections. METHODS: We performed Mendelian randomisation (MR) using 151 variants from genome-wide association studies of blood eosinophils in UK Biobank/INTERVAL, and respiratory traits in UK Biobank/SpiroMeta, using methods relying on different assumptions for validity. We performed multivariable analyses using eight cell types where there was possible evidence of causation by eosinophils. MEASUREMENTS AND MAIN RESULTS: Causal estimates derived from individual variants were highly heterogeneous, which may arise from pleiotropy. The average effect of raising eosinophils was to increase risk of ACO (weighted median OR per SD eosinophils, 1.44 (95%CI 1.19 to 1.74)), and moderate-severe asthma (weighted median OR 1.50 (95%CI 1.23 to 1.83)), and to reduce forced expiratory volume in 1 s (FEV1)/forced vital capacity (FVC) and FEV1 (weighted median estimator, SD FEV1/FVC: -0.054 (95% CI -0.078 to -0.029), effect only prominent in individuals with asthma). CONCLUSIONS: Broad consistency across MR methods may suggest causation by eosinophils (although of uncertain magnitude), yet heterogeneity necessitates caution: other important mechanisms may be responsible for the impairment of respiratory health by these eosinophil-raising variants. These results could suggest that anti-IL5 agents (designed to lower eosinophils) may be valuable in treating other respiratory conditions, including people with overlapping features of asthma and COPD.