Analytical challenges in omics research on asthma and allergy: A National Institute of Allergy and Infectious Diseases workshop.

Studies of asthma and allergy are generating increasing volumes of omics data for analysis and interpretation. The National Institute of Allergy and Infectious Diseases (NIAID) assembled a workshop comprising investigators studying asthma and allergic diseases using omics approaches, omics investigators from outside the field, and NIAID medical and scientific officers to discuss the following areas in asthma and allergy research: genomics, epigenomics, transcriptomics, microbiomics, metabolomics, proteomics, lipidomics, integrative omics, systems biology, and causal inference. Current states of the art, present challenges, novel and emerging strategies, and priorities for progress were presented and discussed for each area. This workshop report summarizes the major points and conclusions from this NIAID workshop. As a group, the investigators underscored the imperatives for rigorous analytic frameworks, integration of different omics data types, cross-disciplinary interaction, strategies for overcoming current limitations, and the overarching goal to improve scientific understanding and care of asthma and allergic diseases.

REALGAR: a web app of integrated respiratory omics data.

MOTIVATION: In the post genome-wide association study (GWAS) era, omics techniques have characterized information beyond genomic variants to include cell and tissue type-specific gene transcription, transcription factor binding sites, expression quantitative trait loci (eQTL) and many other biological layers. Analysis of omics data and its integration has in turn improved the functional interpretation of disease-associated genetic variants. Over 170 000 transcriptomic and epigenomic datasets corresponding to studies of various cell and tissue types under specific disease, treatment and exposure conditions are available in the Gene Expression Omnibus resource. Although these datasets are valuable to guide the design of experimental validation studies to understand the function of disease-associated genetic loci, in their raw form, they are not helpful to experimental researchers who lack adequate computational resources or experience analyzing omics data. We sought to create an integrated re-source of tissue-specific results from omics studies that is guided by disease-specific knowledge to facilitate the design of experiments that can provide biologically meaningful insights into genetic associations. RESULTS: We designed the Reducing Associations by Linking Genes and omics Results web app to provide multi-layered omics information based on results from GWAS, transcriptomic, epigenomic and eQTL studies for gene-centric analysis and visualization. With a focus on asthma datasets, the integrated omics results it contains facilitate the formulation of hypotheses related to airways disease-associated genes and can be addressed with experimental validation studies. AVAILABILITY AND IMPLEMENTATION: The REALGAR web app is available at: http://realgar.org/. The source code is available at: https://github.com/HimesGroup/realgar. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Multiomics analysis identifies BIRC3 as a novel glucocorticoid response-associated gene.

BACKGROUND: Inhaled corticosteroid (ICS) response among patients with asthma is influenced by genetics, but biologically actionable insights based on associations have not been found. Various glucocorticoid response omics data sets are available to interrogate their biological effects. OBJECTIVE: We sought to identify functionally relevant ICS-response genetic associations by integrating complementary multiomics data sets. METHODS: Variants with P values less than 10-4 from a previous ICS-response genome-wide association study were reranked on the basis of integrative scores determined from (1) glucocorticoid receptor- and (2) RNA polymerase II-binding regions inferred from ChIP-Seq data for 3 airway cell types, (3) glucocorticoid response element motifs, (4) differentially expressed genes in response to glucocorticoid exposure according to 20 transcriptomic data sets, and (5) expression quantitative trait loci from GTEx. Candidate variants were tested for association with ICS response and asthma in 6 independent studies. RESULTS: Four variants had significant (q value < 0.05) multiomics integrative scores. These variants were in a locus consisting of 52 variants in high linkage disequilibrium (r2 ≥ 0.8) near glucocorticoid receptor-binding sites by the gene BIRC3. Variants were also BIRC3 expression quantitative trait loci in lung, and 2 were within/near putative glucocorticoid response element motifs. BIRC3 had increased RNA polymerase II occupancy and gene expression, with glucocorticoid exposure in 2 ChIP-Seq and 13 transcriptomic data sets. Some BIRC3 variants in the 52-variant locus were associated (P < .05) with ICS response in 3 independent studies and others with asthma in 1 study. CONCLUSIONS: BIRC3 should be prioritized for further functional studies of ICS response.

Genetic Variation and Sensory Perception of a Pediatric Formulation of Ibuprofen: Can a Medicine Taste Too Good for Some?

There is wide variation in how individuals perceive the chemosensory attributes of liquid formulations of ibuprofen, encompassing both adults and children. To understand personal variation in the taste and chemesthesis properties of this medicine, and how to measure it, our first scientific strategy centered on utilizing trained adult panelists, due to the complex and time-consuming psychophysical tasks needed at this initial stage. We conducted a double-blind cohort study in which panelists underwent whole-genome-wide genotyping and psychophysically evaluated an over-the-counter pediatric medicine containing ibuprofen. Associations between sensory phenotypes and genetic variation near/within irritant and taste receptor genes were determined. Panelists who experienced the urge to cough or throat sensations found the medicine less palatable and sweet, and more irritating. Perceptions varied with genetic ancestry; panelists of African genetic ancestry had fewer chemesthetic sensations, rating the medicine sweeter, less irritating, and more palatable than did those of European genetic ancestry. We discovered a novel association between TRPA1 rs11988795 and tingling sensations, independent of ancestry. We also determined for the first time that just tasting the medicine allowed predictions of perceptions after swallowing, simplifying future psychophysical studies on diverse populations of different age groups needed to understand genetic, cultural-dietary, and epigenetic factors that influence individual perceptions of palatability and, in turn, adherence and the risk of accidental ingestion.

Inhibition of ABCC1 Decreases cAMP Egress and Promotes Human Airway Smooth Muscle Cell Relaxation.

In most living cells, the second-messenger roles for adenosine 3',5'-cyclic monophosphate (cAMP) are short-lived, confined to the intracellular space, and tightly controlled by the binary switch-like actions of Gαs (stimulatory G protein)-activated adenylyl cyclase (cAMP production) and cAMP-specific PDE (cAMP breakdown). Here, by using human airway smooth muscle (HASM) cells in culture as a model, we report that activation of the cell-surface ß2AR (ß2-adrenoceptor), a Gs-coupled GPCR (G protein-coupled receptor), evokes cAMP egress to the extracellular space. Increased extracellular cAMP levels ([cAMP]e) are long-lived in culture and are induced by receptor-dependent and receptor-independent mechanisms in such a way as to define a universal response class of increased intracellular cAMP levels ([cAMP]i). We find that HASM cells express multiple ATP-binding cassette (ABC) membrane transporters, with ABCC1 (ABC subfamily member C 1) being the most highly enriched transcript mapped to MRPs (multidrug resistance-associated proteins). We show that pharmacological inhibition or downregulation of ABCC1 with siRNA markedly reduces ß2AR-evoked cAMP release from HASM cells. Furthermore, inhibition of ABCC1 activity or expression decreases basal tone and increases ß-agonist-induced HASM cellular relaxation. These findings identify a previously unrecognized role for ABCC1 in the homeostatic regulation of [cAMP]i in HASM that may be conserved traits of the Gs-GPCRs (Gs-coupled family of GPCRs). Hence, the general features of this activation mechanism may uncover new disease-modifying targets in the treatment of airflow obstruction in asthma. Surprisingly, we find that serum cAMP levels are elevated in a small cohort of patients with asthma as compared with control subjects, which warrants further investigation.

Skin fibroblast metabolomic profiling reveals that lipid dysfunction predicts the severity of Friedreich's ataxia.

Friedreich's ataxia (FRDA) is an autosomal recessive neurodegenerative disorder caused by a triplet guanine-adenine-adenine (GAA) repeat expansion in intron 1 of the FXN gene, which leads to decreased levels of the frataxin protein. Frataxin is involved in the formation of iron-sulfur (Fe-S) cluster prosthetic groups for various metabolic enzymes. To provide a better understanding of the metabolic status of patients with FRDA, here we used patient-derived fibroblast cells as a surrogate tissue for metabolic and lipidomic profiling by liquid chromatography-high resolution mass spectrometry. We found elevated HMG-CoA and ß-hydroxybutyrate-CoA levels, implying dysregulated fatty acid oxidation, which was further demonstrated by elevated acyl-carnitine levels. Lipidomic profiling identified dysregulated levels of several lipid classes in FRDA fibroblast cells when compared with non-FRDA fibroblast cells. For example, levels of several ceramides were significantly increased in FRDA fibroblast cells; these results positively correlated with the GAA repeat length and negatively correlated with the frataxin protein levels. Furthermore, stable isotope tracing experiments indicated increased ceramide synthesis, especially for long-chain fatty acid-ceramides, in FRDA fibroblast cells compared with ceramide synthesis in healthy control fibroblast cells. In addition, PUFA-containing triglycerides and phosphatidylglycerols were enriched in FRDA fibroblast cells and negatively correlated with frataxin levels, suggesting lipid remodeling as a result of FXN deficiency. Altogether, we demonstrate patient-derived fibroblast cells exhibited dysregulated metabolic capabilities, and their lipid dysfunction predicted the severity of FRDA, making them a useful surrogate to study the metabolic status in FRDA.

Integrative genomic analysis in African American children with asthma finds three novel loci associated with lung function.

Bronchodilator (BD) drugs are commonly prescribed for treatment and management of obstructive lung function present with diseases such as asthma. Administration of BD medication can partially or fully restore lung function as measured by pulmonary function tests. The genetics of baseline lung function measures taken before BD medication have been extensively studied, and the genetics of the BD response itself have received some attention. However, few studies have focused on the genetics of post-BD lung function. To address this gap, we analyzed lung function phenotypes in 1103 subjects from the Study of African Americans, Asthma, Genes, and Environment, a pediatric asthma case-control cohort, using an integrative genomic analysis approach that combined genotype, locus-specific genetic ancestry, and functional annotation information. We integrated genome-wide association study (GWAS) results with an admixture mapping scan of three pulmonary function tests (forced expiratory volume in 1 s [FEV1 ], forced vital capacity [FVC], and FEV1 /FVC) taken before and after albuterol BD administration on the same subjects, yielding six traits. We identified 18 GWAS loci, and five additional loci from admixture mapping, spanning several known and novel lung function candidate genes. Most loci identified via admixture mapping exhibited wide variation in minor allele frequency across genotyped global populations. Functional fine-mapping revealed an enrichment of epigenetic annotations from peripheral blood mononuclear cells, fetal lung tissue, and lung fibroblasts. Our results point to three novel potential genetic drivers of pre- and post-BD lung function: ADAMTS1, RAD54B, and EGLN3.

Epinephrine evokes shortening of human airway smooth muscle cells following ß2 adrenergic receptor desensitization.

Epinephrine (EPI), an endogenous catecholamine involved in the body's fight-or-flight responses to stress, activates α1-adrenergic receptors (α1ARs) expressed on various organs to evoke a wide range of physiological functions, including vasoconstriction. In the smooth muscle of human bronchi, however, the functional role of EPI on α1ARs remains controversial. Classically, evidence suggests that EPI promotes bronchodilation by stimulating ß2-adrenergic receptors (ß2ARs). Conventionally, the selective ß2AR agonism of EPI was thought to be, in part, due to a predominance of ß2ARs and/or a sparse, or lack of α1AR activity in human airway smooth muscle (HASM) cells. Surprisingly, we find that HASM cells express a high abundance of ADRA1B (the α1AR subtype B) and identify a spontaneous "switch-like" activation of α1ARs that evokes intracellular calcium, myosin light chain phosphorylation, and HASM cell shortening. The switch-like responses, and related EPI-induced biochemical and mechanical signals, emerged upon pharmacological inhibition of ß2ARs and/or under experimental conditions that induce ß2AR tachyphylaxis. EPI-induced procontractile effects were abrogated by an α1AR antagonist, doxazosin mesylate (DM). These data collectively uncover a previously unrecognized feed-forward mechanism driving bronchospasm via two distinct classes of G protein-coupled receptors (GPCRs) and provide a basis for reexamining α1AR inhibition for the management of stress/exercise-induced asthma and/or ß2-agonist insensitivity in patients with difficult-to-control, disease subtypes.

CEBPD modulates the airway smooth muscle transcriptomic response to glucocorticoids.

BACKGROUND: CCAAT/Enhancer Binding Protein D (CEBPD), a pleiotropic glucocorticoid-responsive transcription factor, modulates inflammatory responses. Of relevance to asthma, expression of CEBPD in airway smooth muscle (ASM) increases with glucocorticoid exposure. We sought to characterize CEBPD-mediated transcriptomic responses to glucocorticoid exposure in ASM by measuring changes observed after knockdown of CEBPD and its impact on asthma-related ASM function. METHODS: Primary ASM cells derived from four donors were transfected with CEBPD or non-targeting (NT) siRNA and exposed to vehicle control, budesonide (100 nM, 18 h), TNFα (10 ng/ml, 18 h), or both budesonide and TNFα. Subsequently, RNA-Seq was used to measure gene expression levels, and pairwise differential expression results were obtained for exposures versus vehicle and knockdown versus control conditions. Weighted gene co-expression analysis was performed to identify groups of genes with similar expression patterns across the various experimental conditions (i.e., CEBPD knockdown status, exposures). RESULTS: CEBPD knockdown altered expression of 3037 genes under at least one exposure (q-value < 0.05). Co-expression analysis identified sets of 197, 152 and 290 genes that were correlated with CEBPD knockdown status, TNFα exposure status, and both, respectively. JAK-STAT signaling pathway genes, including IL6R and SOCS3, were among those influenced by both TNFα and CEBPD knockdown. Immunoblot assays revealed that budesonide-induced IL-6R protein expression and augmented IL-6-induced STAT3 phosphorylation levels were attenuated by CEBPD knockdown in ASM. CONCLUSIONS: CEBPD modulates glucocorticoid responses in ASM, in part via modulation of IL-6 receptor signaling.

Genome-wide CRISPR screen identifies suppressors of endoplasmic reticulum stress-induced apoptosis.

Sensing misfolded proteins in the endoplasmic reticulum (ER), cells initiate the ER stress response and, when overwhelmed, undergo apoptosis. However, little is known about how cells prevent excessive ER stress response and cell death to restore homeostasis. Here, we report the identification and characterization of cellular suppressors of ER stress-induced apoptosis. Using a genome-wide CRISPR library, we screen for genes whose inactivation further increases ER stress-induced up-regulation of C/EBP homologous protein 10 (CHOP)-the transcription factor central to ER stress-associated apoptosis. Among the top validated hits are two interacting components of the polycomb repressive complex (L3MBTL2 [L(3)Mbt-Like 2] and MGA [MAX gene associated]), and microRNA-124-3 (miR-124-3). CRISPR knockout of these genes increases CHOP expression and sensitizes cells to apoptosis induced by multiple ER stressors, while overexpression confers the opposite effects. L3MBTL2 associates with the CHOP promoter in unstressed cells to repress CHOP induction but dissociates from the promoter in the presence of ER stress, whereas miR-124-3 directly targets the IRE1 branch of the ER stress pathway. Our study reveals distinct mechanisms that suppress ER stress-induced apoptosis and may lead to a better understanding of diseases whose pathogenesis is linked to overactive ER stress response.

Shared genetic and experimental links between obesity-related traits and asthma subtypes in UK Biobank.

BACKGROUND: Clinical and epidemiologic studies have shown that obesity is associated with asthma and that these associations differ by asthma subtype. Little is known about the shared genetic components between obesity and asthma. OBJECTIVE: We sought to identify shared genetic associations between obesity-related traits and asthma subtypes in adults. METHODS: A cross-trait genome-wide association study (GWAS) was performed using 457,822 subjects of European ancestry from the UK Biobank. Experimental evidence to support the role of genes significantly associated with both obesity-related traits and asthma through a GWAS was sought by using results from obese versus lean mouse RNA sequencing and RT-PCR experiments. RESULTS: We found a substantial positive genetic correlation between body mass index and later-onset asthma defined by asthma age of onset at 16 years or greater (Rg = 0.25, P = 9.56 × 10-22). Mendelian randomization analysis provided strong evidence in support of body mass index causally increasing asthma risk. Cross-trait meta-analysis identified 34 shared loci among 3 obesity-related traits and 2 asthma subtypes. GWAS functional analyses identified potential causal relationships between the shared loci and Genotype-Tissue Expression (GTEx) quantitative trait loci and shared immune- and cell differentiation-related pathways between obesity and asthma. Finally, RNA sequencing data from lungs of obese versus control mice found that 2 genes (acyl-coenzyme A oxidase-like [ACOXL] and myosin light chain 6 [MYL6]) from the cross-trait meta-analysis were differentially expressed, and these findings were validated by using RT-PCR in an independent set of mice. CONCLUSIONS: Our work identified shared genetic components between obesity-related traits and specific asthma subtypes, reinforcing the hypothesis that obesity causally increases the risk of asthma and identifying molecular pathways that might underlie both obesity and asthma.

Gene Expression Profiles in Children With Suspected Sepsis.

STUDY OBJECTIVE: Sepsis recognition is a clinical challenge in children. We aim to determine whether peripheral blood gene expression profiles are associated with pathogen type and sepsis severity in children with suspected sepsis. METHODS: This was a prospective pilot observational study in a tertiary pediatric emergency department with a convenience sample of children enrolled. Participants were older than 56 days and younger than 18 years, had suspected sepsis, and had not received broad-spectrum antibiotics in the previous 4 hours. Primary outcome was source pathogen, defined as confirmed bacterial source from sterile body fluid or confirmed viral source. Secondary outcome was sepsis severity, defined as maximum therapy required for shock reversal in the first 3 hospital days. We drew peripheral blood for ribonucleic acid isolation at the sepsis protocol activation, obtained gene expression measures with the GeneChip Human Gene 2.0 ST Array, and conducted differential expression analysis. RESULTS: We collected ribonucleic acid samples from a convenience sample of 122 children with suspected sepsis and 12 healthy controls. We compared the 66 children (54%) with confirmed bacterial or viral infection and found 558 differentially expressed genes, many related to interferon signaling or viral immunity. We did not find statistically significant gene expression differences in patients according to sepsis severity. CONCLUSION: The study demonstrates feasibility of evaluating gene expression profiling data in children evaluated for sepsis in the pediatric emergency department setting. Our results suggest that gene expression profiling may facilitate identification of source pathogen in children with suspected sepsis, which could ultimately lead to improved tailoring of sepsis treatment and antimicrobial stewardship.

Airway Smooth Muscle-Specific Transcriptomic Signatures of Glucocorticoid Exposure.

Glucocorticoids, commonly used asthma controller medications, decrease symptoms in most patients, but some remain symptomatic despite high-dose treatment. The physiological basis underlying the glucocorticoid response, especially in asthma patients with severe, refractory disease, is not fully understood. We sought to identify differences between the transcriptomic response of airway smooth muscle (ASM) cells derived from donors with fatal asthma and donors without asthma to glucocorticoid exposure and to compare ASM-specific changes with those observed in other cell types. In cells derived from nine donors with fatal asthma and eight donors without asthma, RNA sequencing was used to measure ASM transcriptome changes after exposure to budesonide (100 nM 24 h) or control vehicle (DMSO). Differential expression results were obtained for this dataset, as well as 13 publicly available glucocorticoid-response transcriptomic datasets corresponding to seven cell types. Specific genes were differentially expressed in response to glucocorticoid exposure (7,835 and 6,957 in ASM cells derived from donors with fatal asthma and donors without asthma, respectively; adjusted P value < 0.05). Transcriptomic changes in response to glucocorticoid exposure were similar in ASM derived from donors with fatal asthma and donors without asthma, with enriched ontological pathways that included cytokine- and chemokine-related categories. A comparison of glucocorticoid-induced changes in the nonasthma ASM transcriptome with those observed in six other cell types showed that ASM has a distinct glucocorticoid-response signature that is also present in ASM cells from donors with fatal asthma.

Whole blood RNA sequencing reveals a unique transcriptomic profile in patients with ARDS following hematopoietic stem cell transplantation.

BACKGROUND: The acute respiratory distress syndrome (ARDS) is characterized by the acute onset of hypoxemia and bilateral lung infiltrates in response to an inciting event, and is associated with high morbidity and mortality. Patients undergoing allogeneic hematopoietic stem cell transplantation (HSCT) are at increased risk for ARDS. We hypothesized that HSCT patients with ARDS would have a unique transcriptomic profile identifiable in peripheral blood compared to those that did not undergo HSCT. METHODS: We isolated RNA from banked peripheral blood samples from a biorepository of critically ill ICU patients. RNA-Seq was performed on 11 patients with ARDS (5 that had undergone HSCT and 6 that had not) and 12 patients with sepsis without ARDS (5 that that had undergone HCST and 7 that had not). RESULTS: We identified 687 differentially expressed genes between ARDS and ARDS-HSCT (adjusted p-value < 0.01), including IFI44L, OAS3, LY6E, and SPATS2L that had increased expression in ARDS vs. ARDS-HSCT; these genes were not differentially expressed in sepsis vs sepsis-HSCT. Gene ontology enrichment analysis revealed that many differentially expressed genes were related to response to type I interferon. CONCLUSIONS: Our findings reveal significant differences in whole blood transcriptomic profiles of patients with non-HSCT ARDS compared to ARDS-HSCT patients and point toward different immune responses underlying ARDS and ARDS-HSCT that contribute to lung injury.

Current Status and Future Opportunities in Lung Precision Medicine Research with a Focus on Biomarkers. An American Thoracic Society/National Heart, Lung, and Blood Institute Research Statement.

BACKGROUND: Thousands of biomarker tests are either available or under development for lung diseases. In many cases, adoption of these tests into clinical practice is outpacing the generation and evaluation of sufficient data to determine clinical utility and ability to improve health outcomes. There is a need for a systematically organized report that provides guidance on how to understand and evaluate use of biomarker tests for lung diseases. METHODS: We assembled a diverse group of clinicians and researchers from the American Thoracic Society and leaders from the National Heart, Lung, and Blood Institute with expertise in various aspects of precision medicine to review the current status of biomarker tests in lung diseases. Experts summarized existing biomarker tests that are available for lung cancer, pulmonary arterial hypertension, idiopathic pulmonary fibrosis, asthma, chronic obstructive pulmonary disease, sepsis, acute respiratory distress syndrome, cystic fibrosis, and other rare lung diseases. The group identified knowledge gaps that future research studies can address to efficiently translate biomarker tests into clinical practice, assess their cost-effectiveness, and ensure they apply to diverse, real-life populations. RESULTS: We found that the status of biomarker tests in lung diseases is highly variable depending on the disease. Nevertheless, biomarker tests in lung diseases show great promise in improving clinical care. To efficiently translate biomarkers into tests used widely in clinical practice, researchers need to address specific clinical unmet needs, secure support for biomarker discovery efforts, conduct analytical and clinical validation studies, ensure tests have clinical utility, and facilitate appropriate adoption into routine clinical practice. CONCLUSIONS: Although progress has been made toward implementation of precision medicine for lung diseases in clinical practice in certain settings, additional studies focused on addressing specific unmet clinical needs are required to evaluate the clinical utility of biomarkers; ensure their generalizability to diverse, real-life populations; and determine their cost-effectiveness.

Whole-Genome Sequencing of Pharmacogenetic Drug Response in Racially Diverse Children with Asthma.

RATIONALE: Albuterol, a bronchodilator medication, is the first-line therapy for asthma worldwide. There are significant racial/ethnic differences in albuterol drug response. OBJECTIVES: To identify genetic variants important for bronchodilator drug response (BDR) in racially diverse children. METHODS: We performed the first whole-genome sequencing pharmacogenetics study from 1,441 children with asthma from the tails of the BDR distribution to identify genetic association with BDR. MEASUREMENTS AND MAIN RESULTS: We identified population-specific and shared genetic variants associated with BDR, including genome-wide significant (P < 3.53 × 10-7) and suggestive (P < 7.06 × 10-6) loci near genes previously associated with lung capacity (DNAH5), immunity (NFKB1 and PLCB1), and ß-adrenergic signaling (ADAMTS3 and COX18). Functional analyses of the BDR-associated SNP in NFKB1 revealed potential regulatory function in bronchial smooth muscle cells. The SNP is also an expression quantitative trait locus for a neighboring gene, SLC39A8. The lack of other asthma study populations with BDR and whole-genome sequencing data on minority children makes it impossible to perform replication of our rare variant associations. Minority underrepresentation also poses significant challenges to identify age-matched and population-matched cohorts of sufficient sample size for replication of our common variant findings. CONCLUSIONS: The lack of minority data, despite a collaboration of eight universities and 13 individual laboratories, highlights the urgent need for a dedicated national effort to prioritize diversity in research. Our study expands the understanding of pharmacogenetic analyses in racially/ethnically diverse populations and advances the foundation for precision medicine in at-risk and understudied minority populations.

A functional splice variant associated with decreased asthma risk abolishes the ability of gasdermin B to induce epithelial cell pyroptosis.

BACKGROUND: Genetic variants in the chromosomal region 17q21 are consistently associated with asthma. However, mechanistic studies have not yet linked any of the associated variants to a function that could influence asthma, and as a result, the identity of the asthma gene(s) remains elusive. OBJECTIVES: We sought to identify and characterize functional variants in the 17q21 locus. METHODS: We used the Exome Aggregation Consortium browser to identify coding (amino acid-changing) variants in the 17q21 locus. We obtained asthma association measures for these variants in both the Genetic Epidemiology Research in Adult Health and Aging (GERA) cohort (16,274 cases and 38,269 matched controls) and the EVE Consortium study (5,303 asthma cases and 12,560 individuals). Gene expression and protein localization were determined by quantitative RT-PCR and fluorescence immunostaining, respectively. Molecular and cellular studies were performed to determine the functional effects of coding variants. RESULTS: Two coding variants (rs2305480 and rs11078928) of the gasdermin B (GSDMB) gene in the 17q21 locus were associated with lower asthma risk in both GERA (odds ratio, 0.92; P = 1.01 × 10-6) and EVE (odds ratio, 0.85; joint PEVE = 1.31 × 10-13). In GERA, rs11078928 had a minor allele frequency (MAF) of 0.45 in unaffected (nonasthmatic) controls and 0.43 in asthma cases. For European Americans in EVE, the MAF of rs2305480 was 0.45 for controls and 0.39 for cases; for all EVE subjects, the MAF was 0.32 for controls and 0.27 for cases. GSDMB is highly expressed in differentiated airway epithelial cells, including the ciliated cells. We found that, when the GSDMB protein is cleaved by inflammatory caspase-1 to release its N-terminal fragment, potent pyroptotic cell death is induced. The splice variant rs11078928 deletes the entire exon 6, which encodes 13 amino acids in the critical N-terminus, and abolishes the pyroptotic activity of the GSDMB protein. CONCLUSIONS: Our study identified a functional asthma variant in the GSDMB gene of the 17q21 locus and implicates GSDMB-mediated epithelial cell pyroptosis in pathogenesis.

TGF-ß1 Evokes Human Airway Smooth Muscle Cell Shortening and Hyperresponsiveness via Smad3.

Transforming growth factor ß1 (TGF-ß1), a cytokine whose levels are elevated in the airways of patients with asthma, perpetuates airway inflammation and modulates airway structural cell remodeling. However, the role of TGF-ß1 in excessive airway narrowing in asthma, or airway hyperresponsiveness (AHR), remains unclear. In this study, we set out to investigate the direct effects of TGF-ß1 on human airway smooth muscle (HASM) cell shortening and hyperresponsiveness. The dynamics of AHR and single-cell excitation-contraction coupling were measured in human precision-cut lung slices and in isolated HASM cells using supravital microscopy and magnetic twisting cytometry, respectively. In human precision-cut lung slices, overnight treatment with TGF-ß1 significantly augmented basal and carbachol-induced bronchoconstriction. In isolated HASM cells, TGF-ß1 increased basal and methacholine-induced cytoskeletal stiffness in a dose- and time-dependent manner. TGF-ß1-induced single-cell contraction was corroborated by concomitant increases in myosin light chain and myosin phosphatase target subunit 1 phosphorylation levels, which were attenuated by small interfering RNA-mediated knockdown of Smad3 and pharmacological inhibition of Rho kinase. Strikingly, these physiological effects of TGF-ß1 occurred through a RhoA-independent mechanism, with little effect on HASM cell [Ca2+]i levels. Together, our data suggest that TGF-ß1 enhances HASM excitation-contraction coupling pathways to induce HASM cell shortening and hyperresponsiveness. These findings reveal a potential link between airway injury-repair responses and bronchial hyperreactivity in asthma, and define TGF-ß1 signaling as a potential target to reduce AHR in asthma.

PDE8 Is Expressed in Human Airway Smooth Muscle and Selectively Regulates cAMP Signaling by ß2-Adrenergic Receptors and Adenylyl Cyclase 6.

Two cAMP signaling compartments centered on adenylyl cyclase (AC) exist in human airway smooth muscle (HASM) cells, one containing ß2-adrenergic receptor AC6 and another containing E prostanoid receptor AC2. We hypothesized that different PDE isozymes selectively regulate cAMP signaling in each compartment. According to RNA-sequencing data, 18 of 24 PDE genes were expressed in primary HASM cells derived from age- and sex-matched donors with and without asthma. PDE8A was the third most abundant of the cAMP-degrading PDE genes, after PDE4A and PDE1A. Knockdown of PDE8A using shRNA evoked twofold greater cAMP responses to 1 µM forskolin in the presence of 3-isobutyl-1-methylxanthine. Overexpression of AC2 did not alter this response, but overexpression of AC6 increased cAMP responses an additional 80%. We examined cAMP dynamics in live HASM cells using a fluorescence sensor. PF-04957325, a PDE8-selective inhibitor, increased basal cAMP concentrations by itself, indicating a significant basal level of cAMP synthesis. In the presence of an AC inhibitor to reduce basal signaling, PF-04957325 accelerated cAMP production and increased the inhibition of cell proliferation induced by isoproterenol, but it had no effect on cAMP concentrations or cell proliferation regulated by prostaglandin E2. Lipid raft fractionation of HASM cells revealed PDE8A immunoreactivity in buoyant fractions containing caveolin-1 and AC5/6 immunoreactivity. Thus, PDE8 is expressed in lipid rafts of HASM cells, where it specifically regulates ß2-adrenergic receptor AC6 signaling without effects on signaling by the E prostanoid receptors 2/4-AC2 complex. In airway diseases such as asthma and chronic obstructive pulmonary disease, PDE8 may represent a novel therapeutic target to modulate HASM responsiveness and airway remodeling.