Age-associated differences in the human lung extracellular matrix.

Extracellular matrix (ECM) remodeling has been associated with chronic lung diseases. However, information about specific age-associated differences in lung ECM is currently limited. In this study, we aimed to identify and localize age-associated ECM differences in human lungs using comprehensive transcriptomic, proteomic, and immunohistochemical analyses. Our previously identified age-associated gene expression signature of the lung was re-analyzed limiting it to an aging signature based on 270 control patients (37-80 years) and focused on the Matrisome core geneset using geneset enrichment analysis. To validate the age-associated transcriptomic differences on protein level, we compared the age-associated ECM genes (false discovery rate, FDR < 0.05) with a profile of age-associated proteins identified from a lung tissue proteomics dataset from nine control patients (49-76 years) (FDR < 0.05). Extensive immunohistochemical analysis was used to localize and semi-quantify the age-associated ECM differences in lung tissues from 62 control patients (18-82 years). Comparative analysis of transcriptomic and proteomic data identified seven ECM proteins with higher expression with age at both gene and protein levels: COL1A1, COL6A1, COL6A2, COL14A1, FBLN2, LTBP4, and LUM. With immunohistochemistry, we demonstrated higher protein levels with age for COL6A2 in whole tissue, parenchyma, airway wall, and blood vessel, for COL14A1 and LUM in bronchial epithelium, and COL1A1 in lung parenchyma. Our study revealed that higher age is associated with lung ECM remodeling, with specific differences occurring in defined regions within the lung. These differences may affect lung structure and physiology with aging and as such may increase susceptibility to developing chronic lung diseases.NEW & NOTEWORTHY We identified seven age-associated extracellular matrix (ECM) proteins, i.e., COL1A1, COL6A1, COL6A2 COL14A1, FBLN2, LTBP4, and LUM with higher transcript and protein levels in human lung tissue with age. Extensive immunohistochemical analysis revealed significant age-associated differences for COL6A2 in whole tissue, parenchyma, airway wall, and vessel, for COL14A1 and LUM in bronchial epithelium, and COL1A1 in parenchyma. Our findings lay a new foundation for the investigation of ECM differences in age-associated chronic lung diseases.

Genetic Associations and Architecture of Asthma-COPD Overlap.

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.

Prioritization of candidate causal genes for asthma in susceptibility loci derived from UK Biobank.

To identify candidate causal genes of asthma, we performed a genome-wide association study (GWAS) in UK Biobank on a broad asthma definition (n = 56,167 asthma cases and 352,255 controls). We then carried out functional mapping through transcriptome-wide association studies (TWAS) and Mendelian randomization in lung (n = 1,038) and blood (n = 31,684) tissues. The GWAS reveals 72 asthma-associated loci from 116 independent significant variants (PGWAS < 5.0E-8). The most significant lung TWAS gene on 17q12-q21 is GSDMB (PTWAS = 1.42E-54). Other TWAS genes include TSLP on 5q22, RERE on 1p36, CLEC16A on 16p13, and IL4R on 16p12, which all replicated in GTEx lung (n = 515). We demonstrate that the largest fold enrichment of regulatory and functional annotations among asthma-associated variants is in the blood. We map 485 blood eQTL-regulated genes associated with asthma and 50 of them are causal by Mendelian randomization. Prioritization of druggable genes reveals known (IL4R, TSLP, IL6, TNFSF4) and potentially new therapeutic targets for asthma.

Detection of respiratory syncytial virus defective genomes in nasal secretions is associated with distinct clinical outcomes.

Respiratory syncytial virus (RSV) causes respiratory illness in children, immunosuppressed individuals and the elderly. However, the viral factors influencing the clinical outcome of RSV infections remain poorly defined. Defective viral genomes (DVGs) can suppress virus replication by competing for viral proteins and by stimulating antiviral immunity. We studied the association between detection of DVGs of the copy-back type and disease severity in three RSV A-confirmed cohorts. In hospitalized children, detection of DVGs in respiratory samples at or around the time of admission associated strongly with more severe disease, higher viral load and a stronger pro-inflammatory response. Interestingly, in experimentally infected adults, the presence of DVGs in respiratory secretions differentially associated with RSV disease severity depending on when DVGs were detected. Detection of DVGs early after infection associated with low viral loads and mild disease, whereas detection of DVGs late after infection, especially if DVGs were present for prolonged periods, associated with high viral loads and severe disease. Taken together, we demonstrate that the kinetics of DVG accumulation and duration could predict clinical outcome of RSV A infection in humans, and thus could be used as a prognostic tool to identify patients at risk of worse clinical disease.

Global Landscape of Clostridioides Difficile Phylogeography, Antibiotic Susceptibility, and Toxin Polymorphisms by Post-Hoc Whole-Genome Sequencing from the MODIFY I/II Studies.

INTRODUCTION: Clostridioides (Clostridium) difficile infection, the leading cause of healthcare-associated diarrhea, represents a significant burden on global healthcare systems. Despite being a global issue, information on C. difficile from a global perspective is lacking. The aim of this study is to model the global phylogeography of clinical C. difficile. METHODS: Using samples collected from the MODIFY I and II studies (NCT01241552, NCT01513239), we performed whole-genome sequencing of 1501 clinical isolates including 37 novel sequence types (STs), representing the largest worldwide collection to date. RESULTS: Our data showed ribotypes, multi-locus sequence typing clades, and whole-genome phylogeny were in good accordance. The clinical C. difficile genome was found to be more conserved than previously reported (61% core genes), and modest recombination rates of 1.4-5.0 were observed across clades. We observed a significant continent distribution preference among five C. difficile clades (Benjamini-Hochberg corrected Fisher's exact test P < 0.01); moreover, weak association between geographic and genetic distance among ribotypes suggested sources beyond healthcare-related transmission. Markedly different trends of antibiotic susceptibility among lineages and regions were identified, and three novel mutations (in pyridoxamine 5'-phosphate oxidase family protein: Tyr130Ser, Tyr130Cys, and a promoter SNP) associated with metronidazole-reduced susceptibility were discovered on a nim-related gene and its promotor by genome-wide association study. Toxin gene polymorphisms were shown to vary within and between prevalent ribotypes, and novel severe mutations were found on the tcdC toxin regulator protein. CONCLUSION: Our systematic characterization of a global set of clinical trial C. difficile isolates from infected individuals demonstrated the complexity of the genetic makeup of this pathogenic organism. The geographic variability of clades, variability in toxin genes, and mutations associated with antibiotic susceptibility indicate a highly complex interaction of C. difficile between host and environment. This dataset will provide a useful resource for validation of findings and future research of C. difficile.

Multi-omics highlights ABO plasma protein as a causal risk factor for COVID-19.

SARS-CoV-2 is responsible for the coronavirus disease 2019 (COVID-19) and the current health crisis. Despite intensive research efforts, the genes and pathways that contribute to COVID-19 remain poorly understood. We, therefore, used an integrative genomics (IG) approach to identify candidate genes responsible for COVID-19 and its severity. We used Bayesian colocalization (COLOC) and summary-based Mendelian randomization to combine gene expression quantitative trait loci (eQTLs) from the Lung eQTL (n = 1,038) and eQTLGen (n = 31,784) studies with published COVID-19 genome-wide association study (GWAS) data from the COVID-19 Host Genetics Initiative. Additionally, we used COLOC to integrate plasma protein quantitative trait loci (pQTL) from the INTERVAL study (n = 3,301) with COVID-19 loci. Finally, we determined any causal associations between plasma proteins and COVID-19 using multi-variable two-sample Mendelian randomization (MR). The expression of 18 genes in lung and/or blood co-localized with COVID-19 loci. Of these, 12 genes were in suggestive loci (PGWAS < 5 × 10-05). LZTFL1, SLC6A20, ABO, IL10RB and IFNAR2 and OAS1 had been previously associated with a heightened risk of COVID-19 (PGWAS < 5 × 10-08). We identified a causal association between OAS1 and COVID-19 GWAS. Plasma ABO protein, which is associated with blood type in humans, demonstrated a significant causal relationship with COVID-19 in the MR analysis; increased plasma levels were associated with an increased risk of COVID-19 and, in particular, severe COVID-19. In summary, our study identified genes associated with COVID-19 that may be prioritized for future investigations. Importantly, this is the first study to demonstrate a causal association between plasma ABO protein and COVID-19.

Aryl hydrocarbon receptor deficiency causes the development of chronic obstructive pulmonary disease through the integration of multiple pathogenic mechanisms.

Emphysema, a component of chronic obstructive pulmonary disease (COPD), is characterized by irreversible alveolar destruction that results in a progressive decline in lung function. This alveolar destruction is caused by cigarette smoke, the most important risk factor for COPD. Only 15%-20% of smokers develop COPD, suggesting that unknown factors contribute to disease pathogenesis. We postulate that the aryl hydrocarbon receptor (AHR), a receptor/transcription factor highly expressed in the lungs, may be a new susceptibility factor whose expression protects against COPD. Here, we report that Ahr-deficient mice chronically exposed to cigarette smoke develop airspace enlargement concomitant with a decline in lung function. Chronic cigarette smoke exposure also increased cleaved caspase-3, lowered SOD2 expression, and altered MMP9 and TIMP-1 levels in Ahr-deficient mice. We also show that people with COPD have reduced expression of pulmonary and systemic AHR, with systemic AHR mRNA levels positively correlating with lung function. Systemic AHR was also lower in never-smokers with COPD. Thus, AHR expression protects against the development of COPD by controlling interrelated mechanisms involved in the pathogenesis of this disease. This study identifies the AHR as a new, central player in the homeostatic maintenance of lung health, providing a foundation for the AHR as a novel therapeutic target and/or predictive biomarker in chronic lung disease.

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.

Rhinovirus and asthma: Challenges and opportunities.

Human rhinoviruses (RVs) are the primary aetiological agent of the common cold. Generally, the associated infection is mild and self-limiting, but may also be associated with bronchiolitis in infants, pneumonia in the immunocompromised and exacerbation in patients with pulmonary conditions such as asthma or chronic obstructive pulmonary disease. Viral infection accounts for as many as two thirds of asthma exacerbations in children and more than half in adults. Allergy and asthma are major risk factors for more frequent and severe RV-related illnesses. The prevalence of RV-induced wheezing will likely continue to increase given that asthma affects a significant proportion of the population, with allergic asthma accounting for the majority. Several new respiratory viruses and their subgroups have been discovered, with various degrees of relevance. This review will focus on RV infection in the context of the epidemiologic evidence, genetic variability, pathobiology, clinical studies in the context of asthma, differences with other viruses including COVID-19 and current treatment interventions.

Gene expression network analysis provides potential targets against SARS-CoV-2.

Cell entry of SARS-CoV-2, the novel coronavirus causing COVID-19, is facilitated by host cell angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2). We aimed to identify and characterize genes that are co-expressed with ACE2 and TMPRSS2, and to further explore their biological functions and potential as druggable targets. Using the gene expression profiles of 1,038 lung tissue samples, we performed a weighted gene correlation network analysis (WGCNA) to identify modules of co-expressed genes. We explored the biology of co-expressed genes using bioinformatics databases, and identified known drug-gene interactions. ACE2 was in a module of 681 co-expressed genes; 10 genes with moderate-high correlation with ACE2 (r > 0.3, FDR < 0.05) had known interactions with existing drug compounds. TMPRSS2 was in a module of 1,086 co-expressed genes; 31 of these genes were enriched in the gene ontology biologic process 'receptor-mediated endocytosis', and 52 TMPRSS2-correlated genes had known interactions with drug compounds. Dozens of genes are co-expressed with ACE2 and TMPRSS2, many of which have plausible links to COVID-19 pathophysiology. Many of the co-expressed genes are potentially targetable with existing drugs, which may accelerate the development of COVID-19 therapeutics.

Neutrophilic inflammation in the respiratory mucosa predisposes to RSV infection.

The variable outcome of viral exposure is only partially explained by known factors. We administered respiratory syncytial virus (RSV) to 58 volunteers, of whom 57% became infected. Mucosal neutrophil activation before exposure was highly predictive of symptomatic RSV disease. This was associated with a rapid, presymptomatic decline in mucosal interleukin-17A (IL-17A) and other mediators. Conversely, those who resisted infection showed presymptomatic activation of IL-17- and tumor necrosis factor-related pathways. Vulnerability to infection was not associated with baseline microbiome but was reproduced in mice by preinfection chemokine-driven airway recruitment of neutrophils, which caused enhanced disease mediated by pulmonary CD8+ T cell infiltration. Thus, mucosal neutrophilic inflammation at the time of RSV exposure enhances susceptibility, revealing dynamic, time-dependent local immune responses before symptom onset and explaining the as-yet unpredictable outcomes of pathogen exposure.

Genetic regulation of gene expression of MIF family members in lung tissue.

Macrophage migration inhibitory factor (MIF) is a cytokine found to be associated with chronic obstructive pulmonary disease (COPD). However, there is no consensus on how MIF levels differ in COPD compared to control conditions and there are no reports on MIF expression in lung tissue. Here we studied gene expression of members of the MIF family MIF, D-Dopachrome Tautomerase (DDT) and DDT-like (DDTL) in a lung tissue dataset with 1087 subjects and identified single nucleotide polymorphisms (SNPs) regulating their gene expression. We found higher MIF and DDT expression in COPD patients compared to non-COPD subjects and found 71 SNPs significantly influencing gene expression of MIF and DDTL. Furthermore, the platform used to measure MIF (microarray or RNAseq) was found to influence the splice variants detected and subsequently the direction of the SNP effects on MIF expression. Among the SNPs found to regulate MIF expression, the major LD block identified was linked to rs5844572, a SNP previously found to be associated with lower diffusion capacity in COPD. This suggests that MIF may be contributing to the pathogenesis of COPD, as SNPs that influence MIF expression are also associated with symptoms of COPD. Our study shows that MIF levels are affected not only by disease but also by genetic diversity (i.e. SNPs). Since none of our significant eSNPs for MIF or DDTL have been described in GWAS for COPD or lung function, MIF expression in COPD patients is more likely a consequence of disease-related factors rather than a cause of the disease.

Genetic Association Reveals Protection against Recurrence of Clostridium difficile Infection with Bezlotoxumab Treatment.

Bezlotoxumab is a human monoclonal antibody against Clostridium difficile toxin B, indicated to prevent recurrence of C. difficile infection (rCDI) in high-risk adults receiving antibacterial treatment for CDI. An exploratory genome-wide association study investigated whether human genetic variation influences bezlotoxumab response. DNA from 704 participants who achieved initial clinical cure in the phase 3 MODIFY I/II trials was genotyped. Single nucleotide polymorphisms (SNPs) and human leukocyte antigen (HLA) imputation were performed using IMPUTE2 and HIBAG, respectively. A joint test of genotype and genotype-by-treatment interaction in a logistic regression model was used to screen genetic variants associated with response to bezlotoxumab. The SNP rs2516513 and the HLA alleles HLA-DRB1*07:01 and HLA-DQA1*02:01, located in the extended major histocompatibility complex on chromosome 6, were associated with the reduction of rCDI in bezlotoxumab-treated participants. Carriage of a minor allele (homozygous or heterozygous) at any of the identified loci was related to a larger difference in the proportion of participants experiencing rCDI versus placebo; the effect was most prominent in the subgroup at high baseline risk for rCDI. Genotypes associated with an improved bezlotoxumab response showed no association with rCDI in the placebo cohort. These data suggest that a host-driven, immunological mechanism may impact bezlotoxumab response. Trial registration numbers are as follows: NCT01241552 (MODIFY I) and NCT01513239 (MODIFY II).IMPORTANCEClostridium difficile infection is associated with significant clinical morbidity and mortality; antibacterial treatments are effective, but recurrence of C. difficile infection is common. In this genome-wide association study, we explored whether host genetic variability affected treatment responses to bezlotoxumab, a human monoclonal antibody that binds C. difficile toxin B and is indicated for the prevention of recurrent C. difficile infection. Using data from the MODIFY I/II phase 3 clinical trials, we identified three genetic variants associated with reduced rates of C. difficile infection recurrence in bezlotoxumab-treated participants. The effects were most pronounced in participants at high risk of C. difficile infection recurrence. All three variants are located in the extended major histocompatibility complex on chromosome 6, suggesting the involvement of a host-driven immunological mechanism in the prevention of C. difficile infection recurrence.

Phenotypic and functional translation of IL1RL1 locus polymorphisms in lung tissue and asthmatic airway epithelium.

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.

Letermovir Resistance Analysis in a Clinical Trial of Cytomegalovirus Prophylaxis for Hematopoietic Stem Cell Transplant Recipients.

BACKGROUND: Letermovir (LET), a cytomegalovirus (CMV) deoxyribonucleic acid (DNA) terminase inhibitor, was recently approved for prophylaxis of CMV infection in adult CMV-seropositive recipients of allogeneic hematopoietic stem cell transplantation. Cytomegalovirus genotyping was performed to identify LET-resistance-associated variants (RAVs) among subjects in a Phase 3 trial. METHODS: The CMV UL56 and UL89 genes, encoding subunits of CMV DNA terminase, were sequenced from plasma collected from subjects with clinically significant CMV infection (CS-CMVi). Novel variants were evaluated by recombinant phenotyping to assess their potential to confer resistance to LET. RESULTS: Genotyping was successful for 50 of 79 LET subjects with CS-CMVi. Resistance-associated variants (encoding pUL56 V236M and C325W) were detected independently in subjects 1 and 3 who experienced CS-CMVi while receiving LET prophylaxis, and 2 other variants (encoding pUL56 E237G and R369T) were detected >3 weeks after subjects 2 and 3, respectively, had discontinued LET prophylaxis and received preemptive therapy with ganciclovir. CONCLUSIONS: The detected incidence of CMV resistance among subjects who received LET as prophylaxis in this Phase 3 trial was low. The LET RAVs that were detected mapped to the CMV UL56 gene at positions associated with reduced susceptibility to LET based on resistance selections in cell culture.

Genome-Wide Association Study of Susceptibility to Idiopathic Pulmonary Fibrosis.

Rationale: Idiopathic pulmonary fibrosis (IPF) is a complex lung disease characterized by scarring of the lung that is believed to result from an atypical response to injury of the epithelium. Genome-wide association studies have reported signals of association implicating multiple pathways including host defense, telomere maintenance, signaling, and cell-cell adhesion.Objectives: To improve our understanding of factors that increase IPF susceptibility by identifying previously unreported genetic associations.Methods: We conducted genome-wide analyses across three independent studies and meta-analyzed these results to generate the largest genome-wide association study of IPF to date (2,668 IPF cases and 8,591 controls). We performed replication in two independent studies (1,456 IPF cases and 11,874 controls) and functional analyses (including statistical fine-mapping, investigations into gene expression, and testing for enrichment of IPF susceptibility signals in regulatory regions) to determine putatively causal genes. Polygenic risk scores were used to assess the collective effect of variants not reported as associated with IPF.Measurements and Main Results: We identified and replicated three new genome-wide significant (P < 5 × 10-8) signals of association with IPF susceptibility (associated with altered gene expression of KIF15, MAD1L1, and DEPTOR) and confirmed associations at 11 previously reported loci. Polygenic risk score analyses showed that the combined effect of many thousands of as yet unreported IPF susceptibility variants contribute to IPF susceptibility.Conclusions: The observation that decreased DEPTOR expression associates with increased susceptibility to IPF supports recent studies demonstrating the importance of mTOR signaling in lung fibrosis. New signals of association implicating KIF15 and MAD1L1 suggest a possible role of mitotic spindle-assembly genes in IPF susceptibility.

Genome-wide interaction study of gene-by-occupational exposures on respiratory symptoms.

Respiratory symptoms are important indicators of respiratory diseases. Both genetic and environmental factors contribute to respiratory symptoms development but less is known about gene-environment interactions. We aimed to assess interactions between single nucleotide polymorphisms (SNPs) and occupational exposures on respiratory symptoms cough, dyspnea and phlegm. As identification cohort LifeLines I (n = 7976 subjects) was used. Job-specific exposure was estimated using the ALOHA + job exposure matrix. SNP-by-occupational exposure interactions on respiratory symptoms were tested using logistic regression adjusted for gender, age, and current smoking. SNP-by-exposure interactions with a p-value <10-4 were tested for replication in two independent cohorts: LifeLines II (n = 5260) and the Vlagtwedde-Vlaardingen cohort (n = 1529). The interaction estimates of the replication cohorts were meta-analyzed using PLINK. Replication was achieved when the meta-analysis p-value was <0.05 and the interaction effect had the same direction as in the identification cohort. Additionally, we assessed whether replicated SNPs associated with gene expression by analyzing if they were cis-acting expression quantitative trait loci (eQTL) in lung tissue. In the replication meta-analysis, sixteen out of 477 identified SNP-by-occupational exposure interactions had a p-value <0.05 and 9 of these interactions had the same direction as in the identification cohort. Several identified loci were plausible candidates for respiratory symptoms, such as TMPRSS9, SERPINH1, TOX3, and ARHGAP18. Three replicated SNPs were cis-eQTLs for FCER1A, CHN1, and TIMM13 in lung tissue. Taken together, this genome-wide SNP-by-occupational exposure interaction study in relation to cough, dyspnea, and phlegm identified several suggestive susceptibility genes. Further research should determine if these genes are true susceptibility loci for respiratory symptoms in relation to occupational exposures.

Moderate-to-severe asthma in individuals of European ancestry: a genome-wide association study.

BACKGROUND: Few genetic studies that focus on moderate-to-severe asthma exist. We aimed to identity novel genetic variants associated with moderate-to-severe asthma, see whether previously identified genetic variants for all types of asthma contribute to moderate-to-severe asthma, and provide novel mechanistic insights using expression analyses in patients with asthma. METHODS: In this genome-wide association study, we used a two-stage case-control design. In stage 1, we genotyped patient-level data from two UK cohorts (the Genetics of Asthma Severity and Phenotypes [GASP] initiative and the Unbiased BIOmarkers in PREDiction of respiratory disease outcomes [U-BIOPRED] project) and used data from the UK Biobank to collect patient-level genomic data for cases and controls of European ancestry in a 1:5 ratio. Cases were defined as having moderate-to-severe asthma if they were taking appropriate medication or had been diagnosed by a doctor. Controls were defined as not having asthma, rhinitis, eczema, allergy, emphysema, or chronic bronchitis as diagnosed by a doctor. For stage 2, an independent cohort of cases and controls (1:5) was selected from the UK Biobank only, with no overlap with stage 1 samples. In stage 1 we undertook a genome-wide association study of moderate-to-severe asthma, and in stage 2 we followed up independent variants that reached the significance threshold of p less than 1 × 10-6 in stage 1. We set genome-wide significance at p less than 5 × 10-8. For novel signals, we investigated their effect on all types of asthma (mild, moderate, and severe). For all signals meeting genome-wide significance, we investigated their effect on gene expression in patients with asthma and controls. FINDINGS: We included 5135 cases and 25 675 controls for stage 1, and 5414 cases and 21 471 controls for stage 2. We identified 24 genome-wide significant signals of association with moderate-to-severe asthma, including several signals in innate or adaptive immune-response genes. Three novel signals were identified: rs10905284 in GATA3 (coded allele A, odds ratio [OR] 0·90, 95% CI 0·88-0·93; p=1·76 × 10-10), rs11603634 in the MUC5AC region (coded allele G, OR 1·09, 1·06-1·12; p=2·32 × 10-8), and rs560026225 near KIAA1109 (coded allele GATT, OR 1·12, 1·08-1·16; p=3·06 × 10-9). The MUC5AC signal was not associated with asthma when analyses included mild asthma. The rs11603634 G allele was associated with increased expression of MUC5AC mRNA in bronchial epithelial brush samples via proxy SNP rs11602802; (p=2·50 × 10-5) and MUC5AC mRNA was increased in bronchial epithelial samples from patients with severe asthma (in two independent analyses, p=0·039 and p=0·022). INTERPRETATION: We found substantial shared genetic architecture between mild and moderate-to-severe asthma. We also report for the first time genetic variants associated with the risk of developing moderate-to-severe asthma that regulate mucin production. Finally, we identify candidate causal genes in these loci and provide increased insight into this difficult to treat population. FUNDING: Asthma UK, AirPROM, U-BIOPRED, UK Medical Research Council, and Rosetrees Trust.

The Overlap of Lung Tissue Transcriptome of Smoke Exposed Mice with Human Smoking and COPD.

Genome-wide mRNA profiling in lung tissue from human and animal models can provide novel insights into the pathogenesis of chronic obstructive pulmonary disease (COPD). While 6 months of smoke exposure are widely used, shorter durations were also reported. The overlap of short term and long-term smoke exposure in mice is currently not well understood, and their representation of the human condition is uncertain. Lung tissue gene expression profiles of six murine smoking experiments (n = 48) were obtained from the Gene Expression Omnibus (GEO) and analyzed to identify the murine smoking signature. The "human smoking" gene signature containing 386 genes was previously published in the lung eQTL study (n = 1,111). A signature of mild COPD containing 7 genes was also identified in the same study. The lung tissue gene signature of "severe COPD" (n = 70) contained 4,071 genes and was previously published. We detected 3,723 differentially expressed genes in the 6 month-exposure mice datasets (FDR <0.1). Of those, 184 genes (representing 48% of human smoking) and 1,003 (representing 27% of human COPD) were shared with the human smoking-related genes and the COPD severity-related genes, respectively. There was 4-fold over-representation of human and murine smoking-related genes (P = 6.7 × 10-26) and a 1.4 fold in the severe COPD -related genes (P = 2.3 × 10-12). There was no significant enrichment of the mice and human smoking-related genes in mild COPD signature. These data suggest that murine smoke models are strongly representative of molecular processes of human smoking but less of COPD.