OBJECTIVE: Subphenotypes of asthma may be determined by age onset and atopic status. We sought to characterize early or late onset atopic asthma with fungal or non-fungal sensitization (AAFS or AANFS) and non-atopic asthma (NAA) in children and adults in the Severe Asthma Research Program (SARP). SARP is an ongoing project involving well-phenotyped patients with mild to severe asthma. METHODS: Phenotypic comparisons were performed using Kruskal-Wallis or chi-square test. Genetic association analyses were performed using logistic or linear regression. RESULTS: Airway hyper-responsiveness, total serum IgE levels, and T2 biomarkers showed an increasing trend from NAA to AANFS and then to AAFS. Children and adults with early onset asthma had greater % of AAFS than adults with late onset asthma (46% and 40% vs. 32%; P < 0.00001). In children, AAFS and AANFS had lower % predicted FEV1 (86% and 91% vs. 97%) and greater % of patients with severe asthma than NAA (61% and 59% vs. 43%). In adults with early or late onset asthma, NAA had greater % of patients with severe asthma than AANFS and AAFS (61% vs. 40% and 37% or 56% vs. 44% and 49%). The G allele of rs2872507 in GSDMB had higher frequency in AAFS than AANFS and NAA (0.63 vs. 0.55 and 0.55), and associated with earlier age onset and asthma severity. CONCLUSIONS: Early or late onset AAFS, AANFS, and NAA have shared and distinct phenotypic characteristics in children and adults. AAFS is a complex disorder involving genetic susceptibility and environmental factors.

OBJECTIVE: Genome-wide association studies (GWASs) have identified single nucleotide polymorphisms (SNPs) in chr11p15.5 region associated with asthma and idiopathic interstitial pneumonias (IIPs). We sought to identify functional genes for asthma by combining SNPs and mRNA expression in bronchial epithelial cells (BEC) in the Severe Asthma Research Program (SARP). METHODS: Correlation analyses of mRNA expression of six candidate genes (AP2A2, MUC6, MUC2, MUC5AC, MUC5B, and TOLLIP) and asthma phenotypes were performed in the longitudinal cohort (n = 156) with RNAseq in BEC, and replicated in the cross-sectional cohort (n = 155). eQTL (n = 114) and genetic association analysis of asthma severity (426 severe vs. 531 non-severe asthma) were performed, and compared with previously published GWASs of IIPs and asthma. RESULTS: Higher expression of AP2A2 and MUC5AC and lower expression of MUC5B in BEC were correlated with asthma, asthma exacerbations, and T2 biomarkers (P < 0.01). SNPs associated with asthma and IIPs in previous GWASs were eQTL SNPs for MUC5AC, MUC5B, or TOLLIP, however, they were not in strong linkage disequilibrium. The risk alleles for asthma or protective alleles for IIPs were associated with higher expression of MUC5AC and lower expression of MUC5B. rs11603634, rs12788104, and rs28415845 associated with moderate-to-severe asthma or adult onset asthma in previous GWASs were not associated with asthma severity (P > 0.8). CONCLUSIONS: SNPs associated with asthma in chr11p15.5 region are not associated with asthma severity neither with IIPs. Higher expression of MUC5AC and lower expression of MUC5B are risk for asthma but protective for IIPs.

Background

The Chr17q12-21.2 region is the strongest and most consistently associated region with asthma susceptibility. The functional genes or single nucleotide polymorphisms (SNPs) are not obvious due to linkage disequilibrium.

Objectives

We sought to comprehensively investigate whole-genome sequence and RNA sequence from human bronchial epithelial cells to dissect functional genes/SNPs for asthma severity in the Severe Asthma Research Program.

Methods

Expression quantitative trait loci analysis (n = 114), correlation analysis (n = 156) of gene expression and asthma phenotypes, and pathway analysis were performed in bronchial epithelial cells and replicated. Genetic association for asthma severity (426 severe vs 531 nonsevere asthma) and longitudinal asthma exacerbations (n = 273) was performed.

Results

Multiple SNPs in gasdermin B (GSDMB) associated with asthma severity (odds ratio, >1.25) and longitudinal asthma exacerbations (P < .05). Expression quantitative trait loci analyses identified multiple SNPs associated with expression levels of post-GPI attachment to proteins 3, GSDMB, or gasdermin A (3.1 × 10^-9 -4). Higher expression levels of GSDMB correlated with asthma and greater number of exacerbations (P < .05). Expression levels of GSDMB correlated with genes involved in IFN signaling, MHC class I antigen presentation, and immune system pathways (false-discovery rate-adjusted P < .05). rs1031458 and rs3902920 in GSDMB colocalized with IFN regulatory factor binding sites and associated with GSDMB expression, asthma severity, and asthma exacerbations (P < .05).

Conclusions

By using a unique set of gene expression data from lung cells obtained using bronchoscopy from comprehensively characterized subjects with asthma, we show that SNPs in GSDMB associated with asthma severity, exacerbations, and GSDMB expression levels. Furthermore, its expression levels correlated with asthma exacerbations and antiviral pathways. Thus, GSDMB is a functional gene for both asthma susceptibility and severity.

Combination of IL-6 (non-Type 2 asthma) and FeNO or blood eosinophil count (Type 2 asthma) identified asthma endotypes related to asthma severity, exacerbations, and responsiveness to corticosteroids and potential for response to anti-Type 2 and anti-IL-6 treatment.

Rationale: CC16 is a protein mainly produced by nonciliated bronchial epithelial cells (BECs) that participates in host defense. Reduced CC16 protein concentrations in BAL and serum are associated with asthma susceptibility. Objectives: Few studies have investigated the relationship between CC16 and asthma progression, and none has focused on BECs. In this study, we sought to determine if CC16 mRNA expression levels in BECs are associated with asthma severity. Methods: Association analyses between CC16 mRNA expression levels in BECs (242 asthmatics and 69 control subjects) and asthma-related phenotypes in Severe Asthma Research Program were performed using a generalized linear model. Measurements and Main Results: Low CC16 mRNA expression levels in BECs were significantly associated with asthma susceptibility and asthma severity, high systemic corticosteroids use, high retrospective and prospective asthma exacerbations, and low pulmonary function. Low CC16 mRNA expression levels were significantly associated with high T2 inflammation biomarkers (fractional exhaled nitric oxide and sputum eosinophils). CC16 mRNA expression levels were negatively correlated with expression levels of Th2 genes (IL1RL1, POSTN, SERPINB2, CLCA1, NOS2, and MUC5AC) and positively correlated with expression levels of Th1 and inflammation genes (IL12A and MUC5B). A combination of two nontraditional T2 biomarkers (CC16 and IL-6) revealed four asthma endotypes with different characteristics of T2 inflammation, obesity, and asthma severity. Conclusions: Our findings indicate that low CC16 mRNA expression levels in BECs are associated with asthma susceptibility, severity, and exacerbations, partially through immunomodulation of T2 inflammation. CC16 is a potential nontraditional T2 biomarker for asthma development and progression.

Background

Recent meta-analyses of genome-wide association studies in general populations of European descent have identified 28 loci for lung function.

Objective

We sought to identify novel lung function loci specifically for asthma and to confirm lung function loci identified in general populations.

Methods

Genome-wide association studies of lung function (percent predicted FEV1 [ppFEV1], percent predicted forced vital capacity, and FEV1/forced vital capacity ratio) were performed in 4 white populations of European descent (n = 1544), followed by meta-analyses.

Results

Seven of 28 previously identified lung function loci (HHIP, FAM13A, THSD4, GSTCD, NOTCH4-AGER, RARB, and ZNF323) identified in general populations were confirmed at single nucleotide polymorphism (SNP) levels (P < .05). Four of 32 loci (IL12A, IL12RB1, STAT4, and IRF2) associated with ppFEV1 (P < 10(-4)) belong to the TH1 or IL-12 cytokine family pathway. By using a linear additive model, these 4 TH1 pathway SNPs cumulatively explained 2.9% to 7.8% of the variance in ppFEV1 values in 4 populations (P = 3 × 10(-11)). Genetic scores of these 4 SNPs were associated with ppFEV1 values (P = 2 × 10(-7)) and the American Thoracic Society severe asthma classification (P = .005) in the Severe Asthma Research Program population. TH2 pathway genes (IL13, TSLP, IL33, and IL1RL1) conferring asthma susceptibility were not associated with lung function.

Conclusion

Genes involved in airway structure/remodeling are associated with lung function in both general populations and asthmatic subjects. TH1 pathway genes involved in anti-virus/bacterial infection and inflammation modify lung function in asthmatic subjects. Genes associated with lung function that might affect asthma severity are distinct from those genes associated with asthma susceptibility.