Abstract Rationale Type-1 inflammation in asthma leads to significant airway level changes independent of Type-2 inflammation and may underlie the severity seen in combined Type1/2 high disease. As IFNs can cause airway epithelial injury and impaired repair, we applied a novel interpretable machine learning algorithm, SLIDE, to analyze the transcriptome of airway epithelial brushings of mild-moderate (MMA) and severe asthma (SA) participants in the Immune Mechanisms in Severe Asthma (IMSA) cohort to identify latent/hidden factors of disease that associate with disease severity. Methods The expression of genes identified by SLIDE as part of a co-expression program was assessed in the full IMSA cohort by severity and in an independent cohort of the NHLBI Severe Asthma Research Program (SARP). We also examined associations of these genes with clinical asthma outcomes in both cohorts. Results SLIDE identified aberrant increased expression of membrane-tethered mucin genes, MUC1 and MUC4. MUC1 appeared in each of 5 significant latent factors associated with greater severity and MUC4 in 2/5. Standard DGE analysis also identified the gene CEACAM5 as a predictor of severity. Pathway analysis of latent factor genes revealed airway keratinization which may contribute to airway remodeling and worsening obstruction over time. Expression of MUC1, MUC4, and CEACAM5 were increased in SA compared to HC in both the full IMSA and separate SARP cohorts (p 0.001 for all genes in both cohorts). Expression in MMA was consistently lower than in SA in both cohorts, reaching significance in IMSA (p 0.01) but not consistently in SARP (MUC1 p = 0.04, MUC4 p = 0.27, CEACAM5 p = 0.27). Protein analysis of endobronchial biopsies by immunofluorescence staining showed a trend towards increased MUC1 protein content with evidence of secretion. MUC1, MUC4 and CEACAM5 were negatively correlated with FEV1 % predicted in both cohorts (IMSA: MUC1 r=-0.47, p 0.01, MUC4 r=-0.47, p 0.01, CEACAM5 r=-0.62, p 0.01; SARP: MUC1 r=-0.19, p = 0.05, MUC4 r=-0.29, p 0.01, CEACAM5 r=-0.51, p 0.01). All three genes positively correlated with prospective exacerbations in both cohorts (IMSA: MUC1 r = 0.345, p = 0.01, MUC4 r = 0.216, p = 0.08, CEACAM5 r = 0.368, p 0.01; SARP: MUC1 r = 0.22, p = 0.02, MUC4 r = 0.32, p 0.01, CEACAM5 r = 0.38, p 0.01). Conclusions The expression of MUC1, MUC4 and CEACAM5 increased with asthma severity in the airway epithelium and correlated with both prospective exacerbations and obstruction. As pathway analysis of genes in the significant latent factors were associated with airway keratinization aberrant expression of these mucins (with known signaling functions) may contribute to maladaptive repair and airway remodeling that drives obstruction over time, contributing to worse asthma outcomes. This abstract is funded by: NIH
Gauthier et al. (Fri,) studied this question.