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.

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.

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.

Smoking, DNA Methylation, and Lung Function: a Mendelian Randomization Analysis to Investigate Causal Pathways.

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

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.

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.

Risk factors associated with SARS-CoV-2 infection in a multiethnic cohort of United Kingdom healthcare workers (UK-REACH): A cross-sectional analysis.

BACKGROUND: Healthcare workers (HCWs), particularly those from ethnic minority groups, have been shown to be at disproportionately higher risk of infection with Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) compared to the general population. However, there is insufficient evidence on how demographic and occupational factors influence infection risk among ethnic minority HCWs. METHODS AND FINDINGS: We conducted a cross-sectional analysis using data from the baseline questionnaire of the United Kingdom Research study into Ethnicity and Coronavirus Disease 2019 (COVID-19) Outcomes in Healthcare workers (UK-REACH) cohort study, administered between December 2020 and March 2021. We used logistic regression to examine associations of demographic, household, and occupational risk factors with SARS-CoV-2 infection (defined by polymerase chain reaction (PCR), serology, or suspected COVID-19) in a diverse group of HCWs. The primary exposure of interest was self-reported ethnicity. Among 10,772 HCWs who worked during the first UK national lockdown in March 2020, the median age was 45 (interquartile range [IQR] 35 to 54), 75.1% were female and 29.6% were from ethnic minority groups. A total of 2,496 (23.2%) reported previous SARS-CoV-2 infection. The fully adjusted model contained the following dependent variables: demographic factors (age, sex, ethnicity, migration status, deprivation, religiosity), household factors (living with key workers, shared spaces in accommodation, number of people in household), health factors (presence/absence of diabetes or immunosuppression, smoking history, shielding status, SARS-CoV-2 vaccination status), the extent of social mixing outside of the household, and occupational factors (job role, the area in which a participant worked, use of public transport to work, exposure to confirmed suspected COVID-19 patients, personal protective equipment [PPE] access, aerosol generating procedure exposure, night shift pattern, and the UK region of workplace). After adjustment, demographic and household factors associated with increased odds of infection included younger age, living with other key workers, and higher religiosity. Important occupational risk factors associated with increased odds of infection included attending to a higher number of COVID-19 positive patients (aOR 2.59, 95% CI 2.11 to 3.18 for ≥21 patients per week versus none), working in a nursing or midwifery role (1.30, 1.11 to 1.53, compared to doctors), reporting a lack of access to PPE (1.29, 1.17 to 1.43), and working in an ambulance (2.00, 1.56 to 2.58) or hospital inpatient setting (1.55, 1.38 to 1.75). Those who worked in intensive care units were less likely to have been infected (0.76, 0.64 to 0.92) than those who did not. Black HCWs were more likely to have been infected than their White colleagues, an effect which attenuated after adjustment for other known risk factors. This study is limited by self-selection bias and the cross sectional nature of the study means we cannot infer the direction of causality. CONCLUSIONS: We identified key sociodemographic and occupational risk factors associated with SARS-CoV-2 infection among UK HCWs, and have determined factors that might contribute to a disproportionate odds of infection in HCWs from Black ethnic groups. These findings demonstrate the importance of social and occupational factors in driving ethnic disparities in COVID-19 outcomes, and should inform policies, including targeted vaccination strategies and risk assessments aimed at protecting HCWs in future waves of the COVID-19 pandemic. TRIAL REGISTRATION: The study was prospectively registered at ISRCTN (reference number: ISRCTN11811602).

Genotypic and Phenotypic Spectrum of Foveal Hypoplasia: A Multicenter Study.

PURPOSE: To characterize the genotypic and phenotypic spectrum of foveal hypoplasia (FH). DESIGN: Multicenter, observational study. PARTICIPANTS: A total of 907 patients with a confirmed molecular diagnosis of albinism, PAX6, SLC38A8, FRMD7, AHR, or achromatopsia from 12 centers in 9 countries (n = 523) or extracted from publicly available datasets from previously reported literature (n = 384). METHODS: Individuals with a confirmed molecular diagnosis and availability of foveal OCT scans were identified from 12 centers or from the literature between January 2011 and March 2021. A genetic diagnosis was confirmed by sequence analysis. Grading of FH was derived from OCT scans. MAIN OUTCOME MEASURES: Grade of FH, presence or absence of photoreceptor specialization (PRS+ vs. PRS-), molecular diagnosis, and visual acuity (VA). RESULTS: The most common genetic etiology for typical FH in our cohort was albinism (67.5%), followed by PAX6 (21.8%), SLC38A8 (6.8%), and FRMD7 (3.5%) variants. AHR variants were rare (0.4%). Atypical FH was seen in 67.4% of achromatopsia cases. Atypical FH in achromatopsia had significantly worse VA than typical FH (P < 0.0001). There was a significant difference in the spectrum of FH grades based on the molecular diagnosis (chi-square = 60.4, P < 0.0001). All SLC38A8 cases were PRS- (P = 0.003), whereas all FRMD7 cases were PRS+ (P < 0.0001). Analysis of albinism subtypes revealed a significant difference in the grade of FH (chi-square = 31.4, P < 0.0001) and VA (P = 0.0003) between oculocutaneous albinism (OCA) compared with ocular albinism (OA) and Hermansky-Pudlak syndrome (HPS). Ocular albinism and HPS demonstrated higher grades of FH and worse VA than OCA. There was a significant difference (P < 0.0001) in VA between FRMD7 variants compared with other diagnoses associated with FH. CONCLUSIONS: We characterized the phenotypic and genotypic spectrum of FH. Atypical FH is associated with a worse prognosis than all other forms of FH. In typical FH, our data suggest that arrested retinal development occurs earlier in SLC38A8, OA, HPS, and AHR variants and later in FRMD7 variants. The defined time period of foveal developmental arrest for OCA and PAX6 variants seems to demonstrate more variability. Our findings provide mechanistic insight into disorders associated with FH and have significant prognostic and diagnostic value.

Genetic and clinical characteristics of treatment-resistant depression using primary care records in two UK cohorts.

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.

How to design and implement a university-based COVID-19 testing programme? An evaluation of a novel RT-LAMP COVID-19 testing programme in a UK university.

BACKGROUND: Little is known about how asymptomatic testing as a method to control transmission of COVID-19 can be implemented, and the prevalence of asymptomatic infection within university populations. The objective of this study was to investigate how to effectively set-up and implement a COVID-19 testing programme using novel reverse transcriptase loop-mediated isothermal amplification (RT-LAMP) technology and to quantify the scale of asymptomatic infection on a university campus. METHODS: An observational study to describe the set-up and implementation of a novel COVID-19 testing programme on a UK university campus between September and December 2020. RT-LAMP testing was used to identify asymptomatic cases. RESULTS: A total of 1,673 tests were performed using RT-LAMP during the study period, of which 9 were positive for COVID-19, giving an overall positivity rate of 0.54%, equivalent to a rate in the tested population of 538 cases per 100,000 over the duration of testing. All positive tests were found to be positive on RT-PCR testing, giving a false positive rate of 0%. CONCLUSIONS: This study shows that it is possible to rapidly setup a universal university testing programme for COVID-19 in collaboration with local healthcare providers using RT-LAMP testing. Positive results were comparable to those in the local population, though with a different peak of infection. Further research to inform the design of the testing programme includes focus groups of those who underwent testing and further interrogation of the demographics of those opting to be tested to identify potential access problems or inequalities.

Access to personal protective equipment in healthcare workers during the COVID-19 pandemic in the United Kingdom: results from a nationwide cohort study (UK-REACH).

BACKGROUND: Healthcare workers (HCWs) are at high risk of SARS-CoV-2 infection. Effective use of personal protective equipment (PPE) reduces this risk. We sought to determine the prevalence and predictors of self-reported access to appropriate PPE (aPPE) for HCWs in the UK during the COVID-19 pandemic. METHODS: We conducted cross sectional analyses using data from a nationwide questionnaire-based cohort study administered between December 2020-February 2021. The outcome was a binary measure of self-reported aPPE (access all of the time vs access most of the time or less frequently) at two timepoints: the first national lockdown in the UK in March 2020 (primary analysis) and at the time of questionnaire response (secondary analysis). RESULTS: Ten thousand five hundred eight HCWs were included in the primary analysis, and 12,252 in the secondary analysis. 35.2% of HCWs reported aPPE at all times in the primary analysis; 83.9% reported aPPE at all times in the secondary analysis. In the primary analysis, after adjustment (for age, sex, ethnicity, migration status, occupation, aerosol generating procedure exposure, work sector and region, working hours, night shift frequency and trust in employing organisation), older HCWs and those working in Intensive Care Units were more likely to report aPPE at all times. Asian HCWs (aOR:0.77, 95%CI 0.67-0.89 [vs White]), those in allied health professional and dental roles (vs those in medical roles), and those who saw a higher number of COVID-19 patients compared to those who saw none (≥ 21 patients/week 0.74, 0.61-0.90) were less likely to report aPPE at all times. Those who trusted their employing organisation to deal with concerns about unsafe clinical practice, compared to those who did not, were twice as likely to report aPPE at all times. Significant predictors were largely unchanged in the secondary analysis. CONCLUSIONS: Only a third of HCWs in the UK reported aPPE at all times during the first lockdown and that aPPE had improved later in the pandemic. We also identified key determinants of aPPE during the first UK lockdown, which have mostly persisted since lockdown was eased. These findings have important implications for the safe delivery of healthcare during the pandemic.

Phenotypic and functional translation of IL33 genetics in asthma.

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 FEV1, FEV1/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.

A systematic analysis of protein-altering exonic variants in chronic obstructive pulmonary disease.

Genome-wide association studies (GWASs) have identified regions associated with chronic obstructive pulmonary disease (COPD). GWASs of other diseases have shown an approximately 10-fold overrepresentation of nonsynonymous variants, despite limited exonic coverage on genotyping arrays. We hypothesized that a large-scale analysis of coding variants could discover novel genetic associations with COPD, including rare variants with large effect sizes. We performed a meta-analysis of exome arrays from 218,399 controls and 33,851 moderate-to-severe COPD cases. All exome-wide significant associations were present in regions previously identified by GWAS. We did not identify any novel rare coding variants with large effect sizes. Within GWAS regions on chromosomes 5q, 6p, and 15q, four coding variants were conditionally significant (P < 0.00015) when adjusting for lead GWAS single-nucleotide polymorphisms A common gasdermin B (GSDMB) splice variant (rs11078928) previously associated with a decreased risk for asthma was nominally associated with a decreased risk for COPD [minor allele frequency (MAF) = 0.46, P = 1.8e-4]. Two stop variants in coiled-coil α-helical rod protein 1 (CCHCR1), a gene involved in regulating cell proliferation, were associated with COPD (both P < 0.0001). The SERPINA1 Z allele was associated with a random-effects odds ratio of 1.43 for COPD (95% confidence interval = 1.17-1.74), though with marked heterogeneity across studies. Overall, COPD-associated exonic variants were identified in genes involved in DNA methylation, cell-matrix interactions, cell proliferation, and cell death. In conclusion, we performed the largest exome array meta-analysis of COPD to date and identified potential functional coding variants. Future studies are needed to identify rarer variants and further define the role of coding variants in COPD pathogenesis.

Genome-Wide Gene-by-Smoking Interaction Study of Chronic Obstructive Pulmonary Disease.

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.

Resistome analyses of sputum from COPD and healthy subjects reveals bacterial load-related prevalence of target genes.

BACKGROUND: Antibiotic resistance is a major global threat. We hypothesised that the chronic obstructive pulmonary disease (COPD) airway is a reservoir of antimicrobial resistance genes (ARGs) that associate with microbiome-specific COPD subgroups. OBJECTIVE: To determine the resistance gene profiles in respiratory samples from COPD patients and healthy volunteers. METHODS: Quantitative PCR targeting 279 specific ARGs was used to profile the resistomes in sputum from subjects with COPD at stable, exacerbation and recovery visits (n=55; COPD-BEAT study), healthy controls with (n=7) or without (n=22) exposure to antibiotics in the preceding 12 months (EXCEED study) and in bronchial brush samples from COPD (n=8) and healthy controls (n=7) (EvA study). RESULTS: ARG mean (SEM) prevalence was greater in stable COPD samples (35.2 (1.6)) than in healthy controls (27.6 (1.7); p=0.004) and correlated with total bacterial abundance (r2=0.23; p<0.001). Prevalence of ARG positive signals in individuals was not related to COPD symptoms, lung function or their changes at exacerbation. In the COPD subgroups designated High γProteobacteria and High Firmicutes, ARG prevalence was not different at stable state but significantly declined from stable through exacerbation to recovery in the former (p=0.011) without changes in total bacterial abundance. The ARG patterns were similar in COPD versus health, COPD microbiome-subgroups and between sputum and bronchoscopic samples independent of antibiotic exposure in the last 12 months. CONCLUSIONS: ARGs are highly prevalent in sputum, broadly in proportion to bacterial abundance in both healthy and COPD subjects. Thus, COPD appears to be an ARG reservoir due to high levels of bacterial colonisation.

Identifying potential causal effects of age at menarche: a Mendelian randomization phenome-wide association study.

BACKGROUND: Age at menarche has been associated with various health outcomes. We aimed to identify potential causal effects of age at menarche on health-related traits in a hypothesis-free manner. METHODS: We conducted a Mendelian randomization phenome-wide association study (MR-pheWAS) of age at menarche with 17,893 health-related traits in UK Biobank (n = 181,318) using PHESANT. The exposure of interest was the genetic risk score for age at menarche. We conducted a second MR-pheWAS after excluding SNPs associated with BMI from the genetic risk score, to examine whether results might be due to the genetic overlap between age at menarche and BMI. We followed up a subset of health-related traits to investigate MR assumptions and seek replication in independent study populations. RESULTS: Of the 17,893 tests performed in our MR-pheWAS, we identified 619 associations with the genetic risk score for age at menarche at a 5% false discovery rate threshold, of which 295 were below a Bonferroni-corrected P value threshold. These included potential effects of younger age at menarche on lower lung function, higher heel bone-mineral density, greater burden of psychosocial/mental health problems, younger age at first birth, higher risk of childhood sexual abuse, poorer cardiometabolic health, and lower physical activity. After exclusion of variants associated with BMI, the genetic risk score for age at menarche was related to 37 traits at a 5% false discovery rate, of which 29 were below a Bonferroni-corrected P value threshold. We attempted to replicate findings for bone-mineral density, lung function, neuroticism, and childhood sexual abuse using 5 independent cohorts/consortia. While estimates for lung function, higher bone-mineral density, neuroticism, and childhood sexual abuse in replication cohorts were consistent with UK Biobank estimates, confidence intervals were wide and often included the null. CONCLUSIONS: The genetic risk score for age at menarche was related to a broad range of health-related traits. Follow-up analyses indicated imprecise evidence of an effect of younger age at menarche on greater bone-mineral density, lower lung function, higher neuroticism score, and greater risk of childhood sexual abuse in the smaller replication samples available; hence, these findings need further exploration when larger independent samples become available.

The sputum microbiome is distinct between COPD and health, independent of smoking history.

BACKGROUND: Airway bacterial dysbiosis is a feature of chronic obstructive pulmonary disease (COPD). However, there is limited comparative data of the lung microbiome between healthy smokers, non-smokers and COPD. METHODS: We compared the 16S rRNA gene-based sputum microbiome generated from pair-ended Illumina sequencing of 124 healthy subjects (28 smokers and 96 non-smokers with normal lung function), with single stable samples from 218 COPD subjects collected from three UK clinical centres as part of the COPDMAP consortium. RESULTS: In healthy subjects Firmicutes, Bacteroidetes and Actinobacteria were the major phyla constituting 88% of the total reads, and Streptococcus, Veillonella, Prevotella, Actinomyces and Rothia were the dominant genera. Haemophilus formed only 3% of the healthy microbiome. In contrast, Proteobacteria was the most dominant phylum accounting for 50% of the microbiome in COPD subjects, with Haemophilus and Moraxella at genus level contributing 25 and 3% respectively. There were no differences in the microbiome profile within healthy and COPD subgroups when stratified based on smoking history. Principal coordinate analysis on operational taxonomic units showed two distinct clusters, representative of healthy and COPD subjects (PERMANOVA, p = 0·001). CONCLUSION: The healthy and COPD sputum microbiomes are distinct and independent of smoking history. Our results underline the important role for Gammaproteobacteria in COPD.

Blood eosinophil count and airway epithelial transcriptome relationships in COPD versus asthma.

BACKGROUND: Whether the clinical or pathophysiologic significance of the "treatable trait" high blood eosinophil count in COPD is the same as for asthma remains controversial. We sought to determine the relationship between the blood eosinophil count, clinical characteristics and gene expression from bronchial brushings in COPD and asthma. METHODS: Subjects were recruited into a COPD (emphysema versus airway disease [EvA]) or asthma cohort (Unbiased BIOmarkers in PREDiction of respiratory disease outcomes, U-BIOPRED). We determined gene expression using RNAseq in EvA (n = 283) and Affymetrix microarrays in U-BIOPRED (n = 85). We ran linear regression analysis of the bronchial brushings transcriptional signal versus blood eosinophil counts as well as differential expression using a blood eosinophil > 200 cells/µL as a cut-off. The false discovery rate was controlled at 1% (with continuous values) and 5% (with dichotomized values). RESULTS: There were no differences in age, gender, lung function, exercise capacity and quantitative computed tomography between eosinophilic versus noneosinophilic COPD cases. Total serum IgE was increased in eosinophilic asthma and COPD. In EvA, there were 12 genes with a statistically significant positive association with the linear blood eosinophil count, whereas in U-BIOPRED, 1197 genes showed significant associations (266 positive and 931 negative). The transcriptome showed little overlap between genes and pathways associated with blood eosinophil counts in asthma versus COPD. Only CST1 was common to eosinophilic asthma and COPD and was replicated in independent cohorts. CONCLUSION: Despite shared "treatable traits" between asthma and COPD, the molecular mechanisms underlying these clinical entities are predominately different.

Epigenome-wide association study of lung function level and its change.

Previous reports link differential DNA methylation (DNAme) to environmental exposures that are associated with lung function. Direct evidence on lung function DNAme is, however, limited. We undertook an agnostic epigenome-wide association study (EWAS) on pre-bronchodilation lung function and its change in adults.In a discovery-replication EWAS design, DNAme in blood and spirometry were measured twice, 6-15 years apart, in the same participants of three adult population-based discovery cohorts (n=2043). Associated DNAme markers (p<5×10-7) were tested in seven replication cohorts (adult: n=3327; childhood: n=420). Technical bias-adjusted residuals of a regression of the normalised absolute ß-values on control probe-derived principle components were regressed on level and change of forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and their ratio (FEV1/FVC) in the covariate-adjusted discovery EWAS. Inverse-variance-weighted meta-analyses were performed on results from discovery and replication samples in all participants and never-smokers.EWAS signals were enriched for smoking-related DNAme. We replicated 57 lung function DNAme markers in adult, but not childhood samples, all previously associated with smoking. Markers not previously associated with smoking failed replication. cg05575921 (AHRR (aryl hydrocarbon receptor repressor)) showed the statistically most significant association with cross-sectional lung function (FEV1/FVC: pdiscovery=3.96×10-21 and pcombined=7.22×10-50). A score combining 10 DNAme markers previously reported to mediate the effect of smoking on lung function was associated with lung function (FEV1/FVC: p=2.65×10-20).Our results reveal that lung function-associated methylation signals in adults are predominantly smoking related, and possibly of clinical utility in identifying poor lung function and accelerated decline. Larger studies with more repeat time-points are needed to identify lung function DNAme in never-smokers and in children.