Aberrant mucin expression and keratinization distinguishing severe from mild asthma revealed by interpretable machine learning.

Type 2 (T2) immune cells dominate the airways of patients with mild-moderate asthma (MMA) with a more complex type 1 (T1)-T2 mixed immune response evident in treatment-refractory severe asthma (SA). We hypothesized that comparing the transcriptomes of the airway epithelium of patients with SA and MMA would reveal molecular signatures associated with more severe disease in the context of a complex immune response. Using our interpretable machine learning tool, SLIDE, meaningful latent factors (context-specific gene co-expression networks) were revealed that distinguished SA from MMA. Unexpectedly, an aberrant high expression of normally host-protective, membrane-tethered, and IFN-inducible mucins, MUC1 and MUC4, was identified in SA. Gene networks in the significant latent factors discriminating SA from MMA corresponded to enrichment of a keratinization program in SA airways. Keratinization was marked by increased expression of the stress keratin KRT16, signifying squamous metaplasia suggesting adaptive reprogramming of the airway epithelium in response to chronic stress. These mucins and KRT16 were inversely associated with lung function in 2 separate asthma cohorts. Imaging of endobronchial biopsies revealed significantly higher KRT16 protein expression in SA compared with MMA that strongly correlated with MUC1 protein expression. Our study identifies dysregulated host-protective and maladaptive repair responses in SA distinguishing from MMA.