Citations

Abstract RATIONALE. Significant sex-differences exist in asthma across the lifespan. Asthma prevalence and severity are greater in pre-pubescent boys compared to girls but greater in women after puberty. In women, asthma may improve after menopause when sex hormones wane. To address potential mechanisms underlying these changes, we compared the airway transcriptomics signatures of post-menopausal vs. premenopausal women. METHODS. We used bulk RNA seq data from bronchial epithelial cells from three different cohorts (SARP III, GEO-GSE85567, GEO-GSE201955) to perform a comprehensive analysis of differential gene expression and regulation. Data were combined into a meta-analysis, to compare pre-menopausal (n=114 cases with asthma, 48 controls) to post-menopausal women (n=52 cases with asthma, 13 controls). Differential gene expression between pre- and postmenopausal women was analyzed separately for asthma (Figure 1A) and control groups (not shown). Significant differentially expressed genes (DEGs) were identified by contrasting effect sizes (log fold-changes) between asthma and control groups using a Z-score interaction test. Accordingly, a scatter plot (Figure 1B) was created to represent whether differential gene expression between pre- and post-menopausal conditions is modified by disease status (asthma vs. control). The scatter plot also highlights that these findings are related to asthma per. se. and not to normal aging and menopausal transition. Pathway enrichment analysis was conducted using the differential expression data, followed by differential regulation analysis of the identified pathways. Statistical significance was assessed using a false discovery rate (FDR) threshold of 5%. RESULTS. Out of 15,468 genes that passed quality control, 19 genes were differentially expressed in bronchial epithelial cells from pre- and post-menopausal women with asthma (as compared to controls) (Figure 1A). Of these, 5 genes were down-regulated (ITLN1, POSTN, FCRL3, CCL22, PCSK6), and 14 genes were up-regulated (SPP1, S100A9, A2ML1, TNFRSF10C, CSF3R, OLFM1, C1QL2, NLRP3, HAL, PI3, RNASE1, TMEM178B, TNNI3, RNF175) when comparing pre- and post-menopausal women with asthma (Figure 1A). Of these 19 DEGs, 14 were differentially expressed in pre- as compared to post-menopausal women with asthma, but not in controls (Figure 1B). Comparing pre- to post-menopausal women with asthma, biologic pathways related to the innate and adaptive immune system activation and mitochondrial function (Electron transport Chain oxidative phosphorylation and ATP production) were both differentially regulated and expressed. CONCLUSIONS. This analysis highlights several DEGs and pathways that characterize asthma pathogenesis when comparing pre- and post-menopausal women, and emphasize the role of female sex hormones on the immune system response and mitochondrial function in asthma pathogenesis.

Abstract Introduction: The airway inflammation of asthma is characterized by immune phenotypes, often classified as Type 2 inflammation high or low. Social determinants of health (SDoH) can influence asthma risk and severity, but whether or not SDoH influence immune phenotypes is unclear. The social deprivation index (SDI) is a composite measure of seven neighborhood level characteristics: poverty level, less than 12 years of education, single parent household, living in rented housing unit, living in overcrowded housing unit, households without a car, and non-employed adults under 65 years of age. The scale of SDI ranges from 1 to 100, correlating with increasing deprivation. We hypothesized that social deprivation influences immune phenotype of asthma, and that deprivation effects vary by sex and race. Methods: Using the American Community Survey summary file spanning 2015 to 2019, we linked SDI to participant zip codes for 680 Black and White participants in the Severe Asthma Research Program (SARP), a NHLBI-multicenter network phenotyping asthma across the US from 2001-2024 (table 1). We compared SDI among the cohort by sex, race and asthma severity, and correlated SDI to superoxide dismutase (SOD) and T2-high biomarkers, such as eosinophils, IgE and Fraction of Exhaled Nitric Oxide (FENO). Results: SDI is similar across healthy and asthmatic participants (SDI: median [interquartile range]: healthy: 40 [56]; asthmatics: 51 [53], P = 0.53). When evaluated by race, SDI is significantly higher in Black asthmatics, independent of sex (SDI: Black male: 76 [4], White male: 35 [45]; P<0.0001; Black females: 77 [35], White females: 40 [45]; P<0.00010). Within SARP asthmatics, SDI shows a positive association with the T2 inflammatory markers, log IgE (R=0.13 [90%CI: 0.1, 0.2)] and log FENO (R=0.15 [90%CI: 0.1, 0.2]), while it correlates negatively with log of absolute neutrophil count (R= -0.12 [90%CI: -0.19,-0.1]). The antioxidant enzyme SOD increases with SDI (R=0.11, [90% CI: 0.04, 0.17]). Notably, within Black asthmatic males only, SDI is associated with Type 2 biomarkers, showing a positive correlation with log FENO (R=0.24 [90%CI: 0.1,0.4]) and negative correlation with log IgE levels (R= -0.25, [90% CI: -0.4, -0.1]). No association was found between SDI and eosinophils or with SOD in Black males. Conclusions: Neighborhood level social deprivation is associated with a unique asthma phenotype in Black males. Black males living at high social deprivation have a lower IgE and higher FENO, with no change in SOD activity.

Abstract Rationale: The club cell secretory protein (CC16), encoded by the SCGB1A1 gene, has an anti-inflammatory role in airways diseases, including asthma (Li, AJRCCM 2023; Voraphani, AJRCCM 2023). We hypothesize that additional pathway genes regulated or co-expressed with CC16 have cumulative effects on asthma severity outcomes by modulating airway inflammation and maintaining lung homeostasis. Methods: As published previously, four genes were downregulated in CC16 knockout experiments and correlated with CC16 levels: BPIFA1, SFTPD, LTF, and LYZ (Iannuzo, Front. Immunol 2023). We developed a four-gene biomarker score in 94 patients with asthma from the NHLBI-sponsored Severe Asthma Research Program (SARP3) who underwent research bronchoscopy for lower airway epithelial brushings, RNA isolation, and bulk RNA sequencing. We log-transformed the raw values of each transcript; obtained standardized residuals from regression models that included sex, age, and batch effects; and summed the standardized residuals of all four proteins to estimate a standardized sum score (Zhai, JACI Allergy 2024). Logistic regression models tested the four-gene score, three-gene score without LTF, and individual transcripts for associations with exacerbations during 12-month follow-up and zero-inflated negative binomial regression for number of exacerbations. Results: 30.85% percent of 94 patients experienced at least one exacerbation at the end of follow-up. The CC16-associated four-gene biomarker score showed protective effects on asthma exacerbation risk with a 54% percent reduction in exacerbation risk for every one-fold higher score (OR=0.46, 95% CI: 0.27-0.78, p =0.0040, Figure 1). Transcript expression of BPIFA1 (OR=0.87, 95% CI:0.66-1.16, p=0.337), SFTPD (OR=0.33, 95% CI:0.04-2.52, p=0.283), LYZ (OR=0.62, 95%CI:0.28-1.39, p=0.248), and LTF (OR=0.36, 95% CI:0.18-0.74, p=0.0053) all indicated protective directionality with respect to exacerbation risk. However, only LTF demonstrated a statistically significant association. When we excluded LTF, the remaining three-gene score showed a trend toward association with risk (OR=0.64, 95% CI:0.40-1.03, p=0.068) and was significantly associated with number of exacerbations (RR=0.49, 95% CI:0.26-0.95, p=0.034). Compared to LTF alone (OR=0.48, PPV=0.50, NPV=0.71, Sensitivity=0.14, AUC=0.65), the combined four-gene score (OR=0.46, PPV=0.60, NPV= 0.73, Sensitivity=0.21, AUC=0.68) showed improved predictive performance. Conclusions: We identified a CC16-associated transcriptomic predictive score for asthma exacerbations, with LTF showing the strongest effect. Although the other three genes showed individually weak and non-significant associations with exacerbations, the three gene score was significant for number of exacerbations.We demonstrate the potential of this transcriptomic score as a predictive biomarker for asthma severity-related outcomes. Further research investigating these genes and their interaction with other CC16-associated immune pathways could provide insight into therapeutic targets.

Abstract Asthma is a heterogeneous disease with variable presentation and characteristics. There is a critical need to identify underlying molecular endotypes of asthma. We performed the largest transcriptomic analysis of 808 bronchial epithelial cell samples across 11 independent cohorts, including 3 cohorts from the Severe Asthma Research Program. Using seven datasets (218 patients with asthma, 148 healthy control subjects) as discovery cohorts, we identified 505 differentially expressed genes, which we validated in the remaining four datasets. Unsupervised clustering using the 505 differentially expressed genes identified four reproducible clusters of patients with asthma across all datasets, corresponding to healthy control subjects, patients with mild/moderate asthma, and patients with severe asthma with significant differences in several clinical markers of severity, including pulmonary function, Type 2 inflammation, fractional exhaled nitric oxide, and maximum bronchodilator reversibility. Importantly, we found the same clusters in pediatric patients using nasal lavage fluid cells, demonstrating the gene signature and clusters are not confounded by age and are conserved in both lower and upper airways. The four asthma clusters may represent a unifying framework for understanding the molecular heterogeneity of asthma. Further study could potentially enable a precision medicine approach of matching therapies with patients with asthma most likely to benefit.

AbstractIntroductionPrevious bronchoalveolar lavage fluid (BALF) proteomic analysis has evaluated limited numbers of subjects for only a few proteins of interest, which may differ between asthma and normal controls. Our objective was to examine a more comprehensive inflammatory biomarker panel in quantitative proteomic analysis for a large asthma cohort to identify molecular phenotypes distinguishing severe from nonsevere asthma.MethodsBronchoalveolar lavage fluid from 48 severe and 77 nonsevere adult asthma subjects were assessed for 75 inflammatory proteins, normalized to BALF total protein concentration. Validation of BALF differences was sought through equivalent protein analysis of autologous sputum. Subjects' data, stratified by asthma severity, were analysed by standard statistical tests, principal component analysis and 5 machine learning algorithms.ResultsThe severe group had lower lung function and greater health care utilization. Significantly increased BALF proteins for severe asthma compared to nonsevere asthma were fibroblast growth factor 2 (FGF2), TGFα, IL1Ra, IL2, IL4, CCL8, CCL13 and CXCL7 and significantly decreased were platelet‐derived growth factor a‐a dimer (PDGFaa), vascular endothelial growth factor (VEGF), interleukin 5 (IL5), CCL17, CCL22, CXCL9 and CXCL10. Four protein differences were replicated in sputum. FGF2, PDGFaa and CXCL7 were independently identified by 5 machine learning algorithms as the most important variables for discriminating severe and nonsevere asthma. Increased and decreased proteins identified for the severe cluster showed significant protein–protein interactions for chemokine and cytokine signalling, growth factor activity, and eosinophil and neutrophil chemotaxis differing between subjects with severe and nonsevere asthma.ConclusionThese inflammatory protein results confirm altered airway remodelling and cytokine/chemokine activity recruiting leukocytes into the airways of severe compared to nonsevere asthma as important processes even in stable status.

AbstractBackgroundHeterozygote carriers of potentially pathogenic variants in the cystic fibrosis transmembrane conductance regulator (CFTR) gene have increased asthma risk. However, the frequency and impact of CFTR variation among individuals with asthma is unknown.ObjectiveTo determine whether potentially pathogenic CFTR variants associate with disease severity and whether individuals with two potentially pathogenic variants exist in a severe asthma‐enriched cohort.MethodsWe analyzed sequencing data spanning a 190.5Kb region of CFTR in participants from the Severe Asthma Research Program (SARP1‐3). Potentially pathogenic, rare CFTR variants (frequency < 0.05) were classified as CF‐causing or of varying clinical consequences (VVCC) (CFTR2. org). Regression‐based models tested for association between CFTR genotypes (0‐2 potentially pathogenic variants) and severity outcomes.ResultsOf 1401 participants, 9.5% (134) had one potentially pathogenic variant, occurring more frequently in non‐Hispanic white (NHW, 10.1% [84 of 831]) compared to African American individuals (AA, 5.2% [22 of 426]). We found ≥2 potentially pathogenic CFTR variants in 1.4% (19); 0.5% (4) of NHW and 2.8% (12) of AA. Potentially pathogenic CFTR variant genotypes (≥1 or ≥2 variants) were not cumulatively associated with lung function or exacerbations. In NHW, we found three F508del compound heterozygotes with F508del and a VVCC (two 5 T; TG12[c.1210‐11 T > G] and one Arg1070Trp) and a homozygote for the VVCC, 5 T; TG12.ConclusionsWe found potentially pathogenic CFTR variants within a severe asthma‐enriched cohort, including three compound heterozygote genotypes variably associated with CF in NHW individuals. These findings provide the rationale for CFTR sequencing and phenotyping of CF‐related traits in individuals with severe asthma.

Background: Heterozygote carriers of potentially pathogenic variants in the cystic fibrosis transmembrane conductance regulator ( CFTR ) gene have increased asthma risk. However, the frequency and impact of CFTR variation among individuals with asthma is unknown. Objective: To determine whether potentially pathogenic CFTR variants associate with disease severity and whether individuals with two potentially pathogenic variants exist in a severe asthma-enriched cohort . Methods : We analyzed sequencing data spanning a 190.5Kb region of CFTR in participants from the Severe Asthma Research Program (SARP1-3). Potentially pathogenic, rare CFTR variants (frequency<0.05) were classified as CF-causing or of varying clinical consequences (VVCC) (CFTR2.org). Regression-based models tested for association between CFTR genotypes (0-2 potentially pathogenic variants) and severity outcomes. Results: Of 1401 participants, 9.5% (134) had one potentially pathogenic variant, occurring more frequently in non-Hispanic white (NHW, 10.1% [84 of 831]) compared to African American individuals (AA, 5.2% [22 of 426]). We found ≥2 potentially pathogenic CFTR variants in 1.4% (19); 0.5% (4) of NHW and 2.8% (12) of AA. Potentially pathogenic CFTR variant genotypes (≥1 or ≥2 variants) were not cumulatively associated with lung function or exacerbations. In NHW, we found three F508del compound heterozygotes with F508del and a VVCC (two 5T;TG12[c.1210-11T>G] and one Arg1070Trp) and a homozygote for the VVCC, 5T;TG12. Conclusions: We found potentially pathogenic CFTR variants within a severe asthma-enriched cohort , including three compound heterozygote genotypes variably associated with CF in NHW individuals. These findings provide the rationale for CFTR sequencing and phenotyping of CF-related traits in individuals with severe asthma.