Citations

Abstract Rationale Non-selective beta-blockers (BB) are the standard of care for variceal prophylaxis in patients with cirrhosis. It is unclear whether BB negatively affect exercise capacity in patients with portopulmonary hypertension (POPH) or hepatopulmonary syndrome (HPS). We hypothesize that BB use is associated with shorter six-minute walk distance (6MWD) in HPS and POPH than in patients without pulmonary vascular complications of liver disease. Methods The Pulmonary Vascular Complications of Liver Disease 2 (PVCLD2) Study enrolled patients evaluated for liver transplantation or pulmonary hypertension at eight U.S. centers (2013-2017). We conducted a retrospective analysis of 333 patients with HPS, POPH, and portal hypertension without lung complications to evaluate whether disease state modified the association of BB on 6MWD. Multivariable linear regression adjusted for a priori covariates including age, sex, body mass index (BMI), MELD-Na, and WHO functional class. An interaction term between disease category and BB use was included. Results We included 85 patients with HPS, 48 with POPH, and 200 patients with liver disease. Compared to liver disease controls, those with HPS were younger, had higher MELD-Na scores, and more often had worse WHO functional class. POPH patients were similar to controls in age and liver disease severity but had higher WHO functional class. Sex, BMI, and left ventricular ejection fraction did not differ significantly between groups. Higher age, BMI, MELD-Na, and WHO functional class were all significantly associated with lower 6MWD while male sex was associated with higher 6MWD. Adjusting for age, sex, BMI, MELD-Na score, and WHO functional class, there was a significant interaction between BB use and disease category (p = 0.024). Patients with HPS had a 53.4 m longer 6MWD with BB use (95% CI 4.21- 102.5, p = 0.03 for interaction) compared to liver disease controls, whereas patients with POPH had a 72.4 m longer 6MWD (95% CI 7.60- 137.2, p = 0.03 for interaction) compared to liver disease controls. Conclusion The association of BB with 6MWD in advanced liver disease appears to be modified by whether a person has a pulmonary vascular complication including HPS or POPH. The mechanisms for these differences and their clinical implications require further study. This abstract is funded by: NIH T32 Grant

Abstract Introduction Sotatercept was FDA-approved for Group 1 pulmonary arterial hypertension (PAH) in spring 2024, but patients with portopulmonary hypertension (POPH) were not included in major clinical trials that led to its approval. Therefore, safety and efficacy of sotatercept is not well-understood among patients with POPH. We present a case of severe POPH treated with sotatercept. Case Report A 61 year-old female presented to pulmonary hypertension (PH) clinic for follow-up of POPH. History notable for cirrhosis in setting of Hepatis C virus and portosystemic shunts status post embolization. She was on triple PAH therapy with sildenafil, macitentan, and IV epoprostenol. She had functional class 2 dyspnea on exertion, and transthoracic echocardiogram (TTE) demonstrated mild-moderate RV enlargement and mild dysfunction, with estimated right ventricular systolic pressure (RVSP) 76 mmHg. Right heart catheterization (RHC) revealed mean pulmonary arterial pressure (MPAP) 46 mmHg, pulmonary vascular resistance (PVR) 7.9 wood units (WU), cardiac output (CO) 4.9 L/minute, and low biventricular filling pressures. She was started on sotatercept, with goal of improving hemodynamics, RV function, and ultimately qualifying for POPH Model for End-Stage Liver Disease (MELD) exception. She returned for assessment several months later and reported progressive dyspnea with decreasing oxygen saturations. Arterial blood gas (ABG) revealed decreased PaO2 (56 mmHg, compared with 64 mmHg approximately 3 months prior before sotatercept initiation). She was noted to have cutaneous telangiectasias on face, chest, upper extremities, and back (Figure 1). TTE showed decreased RV size and improved function and RVSP (62 mmHg) but large intrapulmonary shunt per agitated saline (“bubble”) study. RHC demonstrated improvement in hemodynamics but remained inadequate to meet POPH MELD exception criteria (MPAP 36 mmHg, PVR 4.6 WU). Sotatercept was discontinued. Follow-up ABG demonstrated increase in PaO2 (63 mmHg), similar to prior baseline, and TTE approximately 6 weeks after stopping sotatercept demonstrated notable decrease in intrapulmonary shunting, now described as small. Telangiectasias improved as well. Discussion The use of sotatercept in patient with POPH and/or liver disease remains an exciting area of investigation, and additional clinical studies are warranted. Potential sequelae such as increased intrapulmonary shunting may particularly affect this patient population, given their unique and complex multiorgan pathology. In current clinical practice, close monitoring is essential when using sotatercept among patients with POPH. This abstract is funded by: n/a

Abstract Rationale Rural populations in the US experience a higher lung cancer burden compared to urban residents, potentially due to differences in risk factors and access to preventative care and early detection. We sought to characterize tobacco use, smoking cessation, and lung cancer screening patterns by rurality to identify modifiable gaps in the lung cancer care continuum. Methods We analyzed data from current and former smoking adults in the Southern Community Cohort Study, a large prospective cohort recruited from community health centers across 12 southeastern states. Participants aged 65-80 years with at least 9 months of Medicare Part A/B and less than 3 months of Part C coverage between 2014 and 2020 were included. Data were linked with the National Death Index and Medicare claims for outcome ascertainment. Rurality was defined using Rural-Urban Commuting Area (RUCA) codes. We compared demographic characteristics, smoking patterns, use of smoking cessation interventions, and receipt of low-dose computed tomography (LDCT) screening between rural and urban participants over time. Results Among 10,114 adults included (56% urban; 46% rural), 54% were female, 55% identified as Black/African American, about a third reported less than a high school education, and 45% lived in areas in the highest quartile of deprivation. Urban residents were of similar age (66 vs 67 years), and more likely to currently smoke (42% vs 35%), but had fewer baseline smoking pack-years (27 vs 30 pack-years) than rural residents. Approximately 2% died of lung cancer during follow-up. Using the USPSTF 2013 lung cancer screening criteria, 27% of urban and 26% of rural residents were eligible. This increased unevenly to 38% urban and 35% rural under the USPSTF 2021 criteria. LDCT uptake was low, increasing from 0% in 2014 to 6% in 2020 among screening-eligible participants. Use of cessation interventions was similarly low. Among current smoking participants, 7% received smoking cessation counseling, and 3% had a varenicline prescription. These patterns remained relatively stable over time, with no significant rural-urban differences in LDCT and cessation intervention use. Conclusion Despite modest differences in smoking prevalence and screening eligibility, we found comparable lung cancer mortality and similarly low utilization of LDCT and cessation interventions among rural and urban residents. Our findings suggest that, among people in lower socioeconomic groups, shared socioeconomic constraints may extend beyond the rural-urban divide. Further research is needed to understand how social determinants of health may be addressed to improve lung cancer outcomes across the screening continuum. This abstract is funded by: NIH

Abstract Rationale Interstitial lung abnormalities (ILA) are abnormal findings on chest CT scans that may represent early-stages of interstitial lung disease. The presence of fibrosis (traction bronchiectasis and/or honeycombing) identifies those with ILA at an increased risk for mortality, however, determining who will progress remains a challenge, and the role of blood-based biomarkers in ILA risk-stratification is unknown. Methods Data from four prospective cohort studies were used to identify associations between plasma proteins measured with DNA aptamers (SomaScan) and ILA, fibrosis, and mortality. The Genetic Epidemiology of COPD Study (COPDGene) was used to identify protein associations, and findings were replicated by meta-analysis of the Age, Gene/Environment Susceptibility-Reykjavik study (AGES-Reykjavik), SubPopulations and InteRmediate Outcome Measures in COPD Study (SPIROMICS), and Framingham Heart Study (FHS). Protein associations with ILA and with fibrotic ILA compared to non-fibrotic ILA were performed using multivariable logistic regression; among those with ILA, Cox proportional hazards models were used to identify protein associations with mortality. Models were adjusted for age, gender, race, BMI, smoking status, pack-year smoking history, and for mortality, COPD. Benjamini-Hochberg false discovery rate p-values <0.05 were considered significant. Results 524 (10.7%) of 4,899 COPDGene participants had ILA, of whom 57 (10.9%) had fibrosis; of 8,146 AGES-Reykjavik, SPIROMICS and FHS participants included in the meta-analysis, 666 (8.2%) had ILA, of whom 198 (29.7%) had fibrosis. In COPDGene, 230 proteins were associated with ILA; the most significant of which were SFTPD (adjusted odds ratio per standard deviation [SD] of protein [OR]=2.21, p = 3.7E-41) and SFTPB (OR = 1.82, p = 7.9E-37). In the meta-analysis, 159 of the ILA-associated proteins were replicated. Of the 230 ILA-associated proteins in COPDGene, 9 were associated with fibrosis, and 5 of these 9 proteins were associated with mortality among those with ILA. The 5 proteins associated with all three outcomes were WFDC2 (ORILA=1.90, p = 2.6E-34; ORFibrosis=2.00, p = 0.0001; adjusted hazards of mortality per SD of protein [HR]=1.59, p = 3.4E-05), GDF15 (ORILA=1.67, p = 3.6E-23; ORFibrosis=1.78, p = 0.0006; HR = 1.63, p = 1.4E-06), SCGB3A1 (ORILA=1.52, p = 4.5e-14; ORFibrosis=1.80, p = 0.0002; HR = 1.25, p = 0.02), TREM1 (ORILA=1.44, p = 1.4e-14; ORFibrosis=1.62, p = 0.001; HR = 1.25, p = 0.02) ICAM5 (ORILA=1.42, p = 1.2e-11; ORFibrosis=1.79, p = 0.0001; HR = 1.31, p = 0.002). In the meta-analysis, 6 of the 9 fibrosis-associated proteins and all 5 ILA mortality-associated proteins replicated. Conclusion This represents the first plasma proteomic screen demonstrating reproducible biomarker associations not only with ILA presence, but also with measures of severity and adverse outcomes. These results provide compelling targets for further study, including their prospective evaluation as clinical biomarkers for early disease identification and prognostication. This abstract is funded by: Multiple Sources

Abstract Rationale Deficiency of alpha-1 antitrypsin (AAT), the protein encoded by SERPINA1, remains the strongest genetic risk factor for COPD. Individuals heterozygote for the type Z variant, termed the PI*MZ genotype (estimated global prevalence n = 35,000,000), experience varying respiratory morbidity beyond expected from tobacco smoking alone. AAT is a pleiotropic protein which inhibits not only neutrophil elastase but also neutrophil chemotaxis mediated by leukotriene B4 (LTB4) [by antagonizing the LTB4 receptor (LTB4R)] and interleukin-8 (IL-8). Although LTB-4 and IL-8 gene variants and protein levels have independently been implicated in COPD, their potential role in PI*MZ is unknown. Methods COPDGene and SPIROMICS are observational tobacco-smoke and COPD-enriched cohort studies. Participants were phenotyped by spirometry and CT imaging, with whole genome sequencing data available via TOPMed. We analyzed participants with PI*MM (COPDGene n = 9,732; SPIROMICS n = 2,105) and PI*MZ (COPDGene n = 300; SPIROMICS n = 83). LTB-4 (ALOX5AP, LTA4H, LTB4R) and IL-8 (CXCL8, CXCR1, CXCR2) associated gene variants (exonic and intronic) with a minor allele frequency of ≥ 0.05 were extracted. We assessed statistical interactions between MZ status and chest CT outcomes or FEV1% predicted using the SPAGxECCT package with adjustment for relatedness (genetic related matrix) and ancestry (top ten PCs). Baseline protein LTB4R (SPIROMICS only) and IL-8 (both cohorts) protein levels were measured via SomaScan, and statistical interactions between MZ status and AAT concentration were explored. Models were adjusted for age, sex, BMI, smoking status, and pack-years. Given the candidate gene discovery approach, significance was defined as nominal p-value<0.05 Results We identified the following number of variants in genes of the LTB-4 or IL-8 pathways: 87 shared between cohorts; 133 unique to COPDGene; 1 unique to SPIROMICS. LTB-4-related polymorphisms were associated (p < 0.05) with any radiographic outcomes in 22 variants (COPDGene) vs. 23 variants (SPIROMICS), with a single overlap (rs17074984 in ALOX5AP). Two LTB-4-pathway variants were associated with FEV1% predicted in SPIROMICS (both in ALOX5AP), but none in COPDGene. AAT plasma concentrations displayed a statistically significant interaction with FEV1% predicted (p = 0.046) and PRMemphysema (p = 0.013) for LTBR4, but not for IL-8. There was not an interaction by MM/MZ status and protein levels. Conclusions Predominantly LTB-4, and not IL-8, pathway-related genetic variants and protein levels may contribute to disease heterogeneity among PI*MZ individuals. The association of circulating LTB4R protein with reduced lung function, particularly among those with low AAT levels, further suggest interactions between AAT and the leukotriene pathway in PI*MZ individuals. Thus, PI*MZ individuals may uniquely derive benefit from leukotriene-targeting pharmacotherapy. This abstract is funded by: NIH NHLBI K23HL173570

Abstract Rationale Alpha-1 antitrypsin deficiency (AATD) continues to be underdiagnosed. The estimated prevalence of the most common severe deficient genotype, PI*ZZ, is ∼1:3500 or 100,000 Americans but only ∼10,000 such individuals have been recognized. Diagnostic delays persist and these delays are associated with worse disease. High-resolution CT (HRCT) imaging may contain latent phenotypic AATD signatures that are not visibly apparent to clinicians. We hypothesized a deep learning model could identify imaging patterns predictive of AATD from inspiratory HRCT scans. Methods Confirmed genotypes and HRCT scans were pooled from COPDGene, SPIROMICS, and GRADS, yielding 350 MM, 102 MZ, and 69 ZZ participants. Only scans reconstructed with smooth kernels (GE STANDARD, Philips B, Siemens B31f/B35f) and slice thickness < =1.25mm were used. The cohort was selected to have similar distributions of quantitative emphysema, sex, and age by matching five MZ and two MM subjects to each ZZ subject (<10-year age difference, same sex, and <3% difference in %LAA). Tobacco smoke-exposure was more prevalent among the MM/MZ groups given these were predominantly from the tobacco smoke-enriched COPDGene and SPIROMICS cohorts. A DenseNet121 convolutional neural network (CNN) was trained using the MONAI platform to predict genotype by outputting probabilities for MM, MZ, and ZZ. Input to the CNN was a three-channel LAA-910, LAA-950, and LAA-970 map downsampled to 2.0mm3 resolution to remove bias due to differences in scanner and image reconstruction parameters. Five-fold cross-validation was performed. Area under the curve (AUC) was calculated for classification of ZZ vs other genotypes (MM/MZ) and MZ vs MM, with statistical significance assessed at p < 0.05. Results The three groups were similar across age (mean±SD years: MM 57.9±8.8, MZ: 57.5±8.9, ZZ: 55.0±10.7), sex (% female: MM: 48.2%, MZ: 46.7%, ZZ: 47.6), and emphysema (LAA-950%: MM: 13.0±12.2, MZ: 12.4±11.8, ZZ: 13.1±11.9). There was a higher prevalence of ever smokers among MM (100%) and MZ (91.7%), compared to ZZ (47.6%). The model achieved an AUC of 0.806 (p < 0.001) for classification of ZZ vs other genotypes, demonstrating robust discriminatory ability for homozygous AATD. Classification of MZ vs other genotypes did not achieve an AUC significantly greater than 0.5. Conclusions This proof-of-concept study demonstrates that deep learning applied to routine inspiratory HRCT can accurately identify homozygous AATD. These findings support the potential role of imaging-based machine learning as a non-invasive screening tool for undiagnosed ZZ individuals. Future work will focus on cohort enrichment, integration with clinical parameters, and explainability of model-derived imaging features. This abstract is funded by: NIH NHLBI K23HL173570

Abstract Rationale Chronic lower respiratory diseases (CLRD), including chronic obstructive pulmonary disease (COPD), were the third-leading causes of death globally in 2019. Quantitative computed tomography (CT) measures of airway wall and lumen areas are associated with greater respiratory symptoms, reduced FEV1, accelerated lung function decline, and higher risk of incident COPD. Even in individuals without COPD, greater airway wall thickness—corresponding to smaller absolute wall and lumen areas—predicts incident COPD and CLRD-related hospitalizations and mortality. Aims To identify protein signatures and biological pathways underlying airway lumen and wall areas using cross-sectional plasma protein associations with CT measures. Methods SPIROMICS is a multicenter longitudinal case-control study that recruited 2,983 adults with COPD and a 20+ packyear smoking history and controls aged 40-80 years between 2010-2015. Standardized chest CT scans were labeled with airway measures. A total of 1,255 plasma proteins measured by SomaScan were included. Linear regression models adjusted for demographic, clinical, and imaging covariates. Significance was defined by a False Discovery Rate (FDR) p-value (q-value) <0.05 using the Benjamini-Hochberg method. Replication was performed in two external cohorts with inverse variance-weighted meta-analysis. Proteins meeting nominal replication underwent elastic net regression and pathway enrichment (Gene Ontology [GO] and Reactome). Results Among 2,065 participants with CT and proteomic data (mean age 63.2 years, 46.3% women, 61.6% with COPD), proteins meeting replication were significantly associated with lumen area (2 proteins; top hit a disintegrin and metalloproteinase with thrombospondin motifs 5, q = 0.016), wall area (3 proteins; top hit leptin, q = 0.011), and percent wall area (16 proteins; top hit roundabout homolog 2, q = 0.002). Elastic net selected 2 proteins for lumen area, 3 for wall area, and 12 for percent wall area. Pathway enrichment revealed FDR-significant pathways only for percent wall area via GO (50 pathways; top hit axogenesis) and Reactome (1 pathway; diseases of glycosylation). Conclusions This study identified distinct proteins, signatures, and biological pathways associated with CT airway measures of lumen area, wall area, and percent wall area. Proteins linked to airway lumen and wall area are involved in metabolic and vasoactive processes, while those for percent wall area are involved in vascular remodeling, cell adhesion, and extracellular matrix organization—highlighting structural pathophysiologic features of COPD. These findings may inform biomarker discovery and prediction of disease onset and progression. This abstract is funded by: Funding: Daniel E. Guzman, MD, MSE was funded by T32-HL144442. SPIROMICS was supported by contracts from the NIH/NHLBI (HHSN268200900013C, HHSN268200900014C, HHSN268200900015C, HHSN268200900016C, HHSN268200900017C, HHSN268200900018C, HHSN268200900019C, HHSN268200900020C, 75N92024D00012), grants from the NIH/NHLBI (U01HL137880, U24HL141762, R01HL182622, R01HL144718, and R01HL093081), and supplemented by contributions made through the Foundation for the NIH and the COPD Foundation from Amgen; AstraZeneca/MedImmune; Bayer; Bellerophon Therapeutics; Boehringer-Ingelheim Pharmaceuticals, Inc.; Bristol Myers Squibb; Chiesi Farmaceutici S.p.A.; Forest Research Institute, Inc.; Genentech; GlaxoSmithKline; Grifols Therapeutics, Inc.; Ikaria, Inc.; MGC Diagnostics; Novartis Pharmaceuticals Corporation; Nycomed GmbH; Polarean; ProterixBio; Regeneron Pharmaceuticals, Inc.; Sanofi; Sunovion; Takeda Pharmaceutical Company; Theravance Biopharma; Verona; and Mylan/Viatris.

Abstract Rationale Alpha-1 antitrypsin (AAT) deficiency, primarily caused by homozygous SERPINA1 Z mutation, increases the risk of COPD. Individuals carrying one Z and one normal M allele (MZ) have higher COPD risk than MM individuals, although many never develop disease. COPD is a polygenic condition influenced by numerous genetic factors across the genome and by environmental exposures such as smoking. We aimed to examine how polygenic factors and smoking modify risk of airflow obstruction in MZ compared with MM. Methods We included participants with the MZ and MM genotypes from the Framingham Heart Study. We first tested the association between AAT genotype (MZ coded as 1 and MM as 0) and FEV1/FVC, stratified by ever-smoking status, and further evaluated interactions between AAT genotype and smoking pack-years. Next, for each participant we calculated a polygenic risk score (PRS) for FEV1/FVC, which was derived from external genome-wide association studies in the UK Biobank and SpiroMeta Consortium. To reflect risk for airflow limitation, higher PRSs correspond to lower FEV1/FVC. We then tested the interaction between the PRS and AAT genotype, stratified by ever-smoking. Models were adjusted for age, sex, height, current smoking, pack-years (among ever-smokers), familial relatedness, and in PRS models, principal components of genetic ancestry. Results Among 4,220 participants (mean age 53.7 years, 54.2% female, 50.2% ever-smokers, mean FEV1/FVC 74.6%), 133 (3.2%) had MZ genotype. In never-smokers, MZ individuals had 1.72% higher FEV1/FVC compared with MM (P = 0.032), whereas no difference was observed among ever-smokers (β= -0.39%; P = 0.65). Cumulative smoking modified the effect of MZ genotype on FEV1/FVC (interaction β= -0.13% per pack-year; P = 0.003); this interaction persisted among ever-smokers (β= -0.11% per pack-year; P = 0.044). A one-standard-deviation increase in the FEV1/FVC PRS was associated with lower FEV1/FVC in both never- (β= -2.37%; P < 0.0001) and ever-smokers (β= -2.46%; P < 0.0001). A higher PRS was disproportionately associated with reduced FEV1/FVC in MZ compared to MM (interaction β= -1.87% per standard-deviation PRS increase, P = 0.027) among ever-smokers, but not in never-smokers (P = 0.15). Conclusions In this community-based cohort, the AAT MZ genotype was associated with higher FEV1/FVC among non-smokers compared to MM. The interaction between MZ genotype and pack-years suggested a dose-dependent modifying effect of smoking on the MZ-related COPD risk. Among ever-smokers, higher polygenic risk further amplified susceptibility to airflow obstruction in MZ compared with MM. These findings highlight the roles of both polygenic and environmental risk factors in COPD susceptibility in MZ. Replication of the results is needed. This abstract is funded by: This work is supported by an Alpha-1 Foundation Grant and a Boston University Department of Medicine Career Investment Award

Abstract Background Muscle loss in COPD results from both disuse atrophy and systemic inflammation. We sought to identify inflammatory biomarkers associated with loss of pectoralis muscle area in COPD using serial CT scans and peripheral blood proteomics. Design We analyzed five-year longitudinal data from COPDGene, an ongoing observational cohort of current and former tobacco smokers with or without airflow obstruction. We quantified pectoralis major muscle (PMM) area on chest CT scans using a 3D convolutional neural network (CNN) trained on CT volumes and annotations provided by the TotalSegmentator project. We compared participants with CT PMM data and plasma proteins by SOMAScan (1K platform) from Phase 2 and 3 (paired n = 1,191). The correlation coefficient (r) between longitudinal changes in protein and in PMM area was plotted against the false discovery rate (FDR)-adjusted p-values. Proteins significantly associated with muscle loss (blue) or gain (red) were visualized using volcano plots. We then analyzed a subset (paired n = 500) of individuals with BMI>25 in Phase 2 whose BMI declined in Phase 3. Proteomic models were adjusted for age and clinical center. Pathway analysis utilized G:Profiler. Results Of the participants with paired proteomic and PMM data, mean PMM area was 21.7 ± 17.3 cm² at Phase 2 and 21.3 ± 17.9 cm² at Phase 3, with an annual change in PMM area of -0.31 + 0.99 cm². Multiple proteins were associated with muscle loss from Phase 2 to 3 (including leptin, FABP, and growth hormone receptor) and muscle gain (including F172A, GMEB2, and RNF41) (Fig 1a). Significant pathways were related to cancer, growth factor receptor dysregulation, and PI3K/AKT signaling. When restricting to those with BMI>25 who experienced BMI decline at Phase 3, muscle loss was associated with increased leptin and PACAP, versus association of gain with increased GDF118, EF1β, and FJF6 (Fig 1b). In the entire group, leptin declined between Phase 2 to Phase 3 in n = 571 participants and increased or was unchanged in n = 620. Compared with leptin decline, unchanged or increased leptin was associated with significantly greater interval PMM loss (−2.25 ± 9.75 cm² vs + 0.98 ± 9.43 cm²; p < 0.001) (Fig 1c). Conclusion Loss of PMM over five years was associated with several plasma proteins, including FABP and growth hormone receptor. The top biomarker, leptin, is a pleiotropic adipokine that regulates energy balance and inflammation. Integrating imaging and proteomics may identify future therapeutic targets to mitigate skeletal muscle loss in COPD. This abstract is funded by: This work was supported by NHLBI grants U01 HL089897 and U01 HL089856 and by NIH contract 75N92023D00011

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