Erzurum, S C

Abstract Introduction Insulin resistance is prevalent in pulmonary arterial hypertension (PAH) and is associated with worse health outcomes. The triglyceride-glucose (TyG) index as a measure of insulin resistance is superior to HOMA-IR in predicting metabolic syndrome. This study utilizes the Redefining Pulmonary Hypertension through Pulmonary Vascular Disease Phenomics (PVDOMICS) data registry to further examine the relationship between TyG index and reported quality of life (QOL) in patients with PAH. Methods 327 patients with a diagnosis of PAH in PVDOMICS were stratified into insulin resistant (IR) (TyG index >8.55) and insulin sensitive (IS) (TyG < 8.55) groups where TyG = (ln[triglyceride*glucose/2]. QOL measures assessed were the Minnesota Living With Heart Failure (MLHF) Physical and Emotional domain surveys, Short Form (36) (SF-36) Physical and Mental Component Health Surveys (PCS and MCS, respectively), and emPHasis-10 (E10). Reported measures between groups were compared using least squares means adjusted for age, sex, fat mass index (FMI), PAH etiology, 6-minute walk distance (6MWD) and NT-pro-brain natriuretic peptide (NTproBNP). Kaplan-Meier curves and Cox proportional hazards models adjusted for these same variables were performed to assess difference in survival between groups. Results 155 patients with PAH (47%) were classified as IR. Patients with IR were more likely to be female (79% vs 69%, p = 0.04), more likely to be non-Hispanic white, and less likely to be non-Hispanic Black (p = 0.045). There were no differences in age or etiology of PAH between groups. Patients with IR had higher body mass index (30.9 vs 27.3 kg/m2, p < 0.001) with higher FMI (16.2 vs 13.2 kg/m2, p < 0.001). Patients with IR were more likely to have diabetes mellitus (24% vs 11%, p = 0.004). Compared to their IS counterparts and after adjusting for possible confounders, patients with IR had higher (worse QOL) MLHF physical domain (p = 0.04), higher (worse QOL) emPhasis-10 scores (p = 0.03), and lower (worse QOL) SF36 physical component scores (p = 0.03) (Figure). IR and IS groups had no significant differences in the MLHF emotional domain and SF36 mental component scores. After multivariable adjustment, there was no difference in survival between groups (HR 0.91, 95% CI 0.56-1.46). Conclusion Patients with PAH and IR had worse QOL compared to their IS counterparts but no difference in survival in this large prospective cohort. This abstract is funded by: This work was supported by the National Institutes of Health grants K23HL 141584, R01 HL173533

Abstract Introduction Interstitial lung disease-associated pulmonary hypertension (ILD-PH) is a clinically challenging entity that lies on a pathophysiologic spectrum between pulmonary arterial hypertension (PAH) and interstitial lung disease (ILD). We compared the metabolomes of patients with ILD-PH to those with PAH and ILD to characterize the PAH-ILD axis and determine whether ILD-PH patients can be phenotyped based on their metabolomic similarity to PAH. Methods Using untargeted metabolomics profiles from patients with PAH, ILD, and ILD-PH in the large, U.S.-based, NIH-funded, multi-center PVDOMICS cohort, we trained an anchor-based classifier on PAH and ILD metabolomes using supervised PLS-DA with a 75/25 train/test split. Primary and secondary PH or risk-group classifications in PVDOMICS were based on expert adjudication of clinical data including right heart catheterization hemodynamics, pulmonary function testing, chest imaging, and laboratory results. We then applied the trained model to each ILD-PH patient to assign a PAH similarity score. We classified patients as having a PAH-like or ILD-like metabolomic phenotype using an optimal PAH similarity score cutoff of 0.5 then compared PAH metabolomic similarity to expert-adjudicated primary and secondary WSPH classifications in PVDOMICS. We implemented FELLA, a network-based metabolomics enrichment method, to holistically characterize the biologic underpinnings of PAH, ILD-PH, and ILD metabolomes. Results We analyzed 244 patients with PAH, 54 with ILD, and 47 with ILD-PH. The anchor-based classifier demonstrated high train-test performance with a c-statistic of 0.81 and 78% agreement with expert-adjudicated diagnoses. Among the ILD-PH patients, 23 were classified as PAH-like and 24 as ILD-like (74% agreement with expert adjudication, Cohen’s kappa 0.49, Hedge’s g estimate of 0.86). FELLA enrichment revealed core glycolytic and bioenergetic pathways were enriched in both pure PAH and ILD-PH. Comparison of diffusion Z-scores between groups demonstrated subtle but biologically coherent differences in several reaction-level nodes, including lactate, NADPH oxidoreductase, and phosphofructokinase, which were enriched in ILD-PH. Conclusion Patients with ILD-PH demonstrate clinically relevant metabolomic phenotypes. Examination of the network biology underlying ILD-PH may point to potential therapeutic targets in this vulnerable subgroup with lagging clinical outcomes. This abstract is funded by: Pulmonary Fibrosis Foundation PFF Scholars Program

Abstract Rationale The pathophysiology and clinical impact of liver injury in pulmonary hypertension (PH) is poorly understood. Liver injury in PH is often attributed to hepatic congestion, but some studies have failed to find strong associations with hemodynamics, suggesting alternative mechanisms. We sought to define the prevalence of liver injury and patterns of liver injury across the spectrum of PH and to determine whether patterns of liver injury were associated with clinical characteristics, hemodynamics and outcomes. Methods We performed an analysis of the PVDOMICS study, a multicenter, prospective, longitudinal cohort study that enrolled incident and prevalent patients with PH and disease comparators. The prevalence of liver injury, defined as aspartate or alanine aminotransferase (AST or ALT) or total bilirubin above the upper limit of normal, and different patterns of liver injury [hepatocellular (elevated AST or ALT with normal bilirubin), cholestatic (elevated bilirubin) and no liver injury] was defined among participants with known or suspected PH. Associations between patterns of liver injury with demographics, hemodynamics, right ventricular function, exercise capacity and quality of life (QoL) were assessed. We used Cox proportional hazards models adjusted for age and sex to assess the association of liver injury patterns with transplant-free survival. Results Liver injury was present in 27.8% (288/1036) of participants (19.0% cholestatic and 8.8% hepatocellular). Patients with liver injury were less likely to be female and have prevalent disease with differences in PH classification (Table 1). Patients with hepatocellular liver injury had a higher body mass index. A cholestatic liver injury pattern was most consistently associated with worse functional class, biventricular filling pressures, pulmonary vascular resistance, right ventricular function, natriuretic peptide elevation, exercise capacity, and REVEAL risk scores. Total bilirubin levels and right atrial pressure as a marker of systemic congestion were weakly correlated (r2=0.02, p < 0.0001). Compared to no liver injury, cholestatic liver injury was associated with an increased risk of death in age and sex adjusted models (HR 1.6, 95% CI 1.3-2.1, p < 0.001) but hepatocellular liver injury was not (p > 0.05). Conclusions Liver injury is common among patients with known or suspected PH. In contrast to a hepatocellular injury pattern, a cholestatic liver injury pattern was most consistently associated with worse hemodynamics, exercise capacity, QoL and survival. However, the weak relationship between right atrial pressure and bilirubin suggests that systemic congestion alone is unlikely to completely explain cholestatic liver injury. Further study is needed to uncouple the pathophysiological basis of cholestasis in PH. This abstract is funded by: NHLBI

Abstract Rationale Prior studies demonstrate pulmonary artery (PA) enlargement on computed tomography (CT) in pulmonary hypertension (PH). Some patients exhibit markedly enlarged or “aneurysmal” PA dilation, but its clinical relevance remains unclear. We aimed to characterize clinical characteristics, hemodynamics, and survival probability of patients with aneurysmal PA size. Methods Demographics, clinical characteristics, and CT imaging were obtained for individuals with group 1 or 2 PH and matched controls in the Pulmonary Vascular Disease Phenomics cohort. PA diameter was measured using Osirix software. PA diameter was categorized based on previously reported values as normal (<27 mm [females], <29 mm [males]), enlarged (27-39 mm [females], 29-39 mm [males]), or aneurysmal (≥40 mm). Outcomes included mean PA pressure (mPAP) on right heart catheterization (RHC) and time to death or heart/lung transplant. Associations with mortality and mPAP were evaluated using multivariable Cox and linear regression models, adjusting for age, sex, and BMI. Descriptive statistics are presented as median [IQR]. Results 690 subjects met the inclusion criteria for survival analysis (364 WSPH group 1 or matched comparator, 241 group 2 or comparator, 85 healthy control), with 594 undergoing RHC. PA diameter was normal in 163 participants, enlarged in 414, and aneurysmal in 113. Aneurysmal PA, compared with normal or enlarged PA, was associated with younger age (53 ± 14 vs 55 ± 16 and 59 ± 15 years, respectively; p < 0.001), non-White race (18.9% vs 8% and 17.2%; p = 0.034), pro-BNP (1950 ± 5708 vs 531 ±1651 and 110 ± 2258 pg/mL; p < 0.001), PH duration (6.78 ± 6.86 vs 2.24 ± 4.39 and 3.82 ± 5.47 years; p < 0.001), group 1 disease (77% vs 37% and 61%), and congenital heart disease (CHD)-PAH (10.6% vs 0.0% and 4.8%, p < 0.001). Aneurysmal PA subjects had higher mPAP (46 [36-57] vs 19 [16- 25] and 35 [28-47] mmHg, p < 0.001), pulmonary vascular resistance (PVR) (5.5 [3.6-9.4] vs 1.7 [1.0-2.3] and 4.1 [2.4-7.1] Woods units, p < 0.001), cardiac index (CI; 2.73 [2.34-3.29] vs 2.55 [2.13-2.98] and 2.52 [2.12-3.10] L/min/m2, p < 0.02) and likelihood of death or transplant (p < 0.001; Figure 1). Conclusion “Aneurysmal” PA was more common in younger, non-White subjects with longer WSPH group 1 disease duration, CHD, and higher pro-BNP levels. Aneurysmal PA was associated with elevated mPAP, PVR, and CI, and reduced likelihood of survival. These findings highlight the need to further elucidate the pathophysiology and prognostic implications of aneurysmal PA remodeling in PH. This abstract is funded by: NHLBI and the Pulmonary Hypertension Association

Abstract Background We recently reported a strong association between CT markers of pulmonary artery (PA) dilation and mean PA pressure (mPAP) in subjects with World Symposium on Pulmonary Hypertension (WSPH) group 1 and 2 pulmonary hypertension (PH). The relationship between PA size and mPAP in subjects with WSPH group 3 PH, in which parenchymal changes may influence PA dilation, remains uncertain. We sought to investigate this relationship in subjects with and without PH who have significant parenchymal pulmonary disease or hypoxia. Methods We measured PA diameter and PA:aorta (PA:A) ratio on CT for all subjects with WSPH group 3 PH, matched disease comparators, and healthy controls in the Pulmonary Vascular Disease Phenomics (PVDOMICS) cohort. PH and comparator subjects were categorized by group 3 etiology into obstructive (including chronic obstructive pulmonary disease [COPD] and cystic fibrosis), interstitial (including idiopathic pulmonary fibrosis, unspecified interstitial lung disease, combined pulmonary fibrosis and emphysema, and chronic hypersensitivity pneumonitis), and non-parenchymal (including obesity hypoventilation syndrome, obstructive sleep apnea, and extrinsic restrictive disease) subgroups. Primary outcome was mPAP on right heart catheterization (RHC). Association with mPAP was evaluated using multivariable linear regression adjusted for age, sex, and BMI. Results 339 subjects met inclusion criteria (101 WSPH group 3 PH, 153 disease comparators, 85 healthy controls). 83 subjects had obstructive disease, 125 had interstitial disease, and 46 had non-parenchymal disease. PA diameter and PA:A were associated with mPAP in all subjects and in each subgroup (p < 0.001; see figure 1). The association was notably stronger in the non-parenchymal group (ρ = 0.707 for diameter, 0.750 for PA:A) compared to obstructive (ρ = 0.566 for diameter, 0.529 for PA:A) and interstitial (ρ = 0.432 for diameter, 0.373 for PA:A) groups. Change in PA size between 25th and 75th percentile mPAP was lower in the interstitial group (3.88 mm, ratio 0.12) than in the obstructive (4.88 mm, ratio 0.17) and non-parenchymal (4.80 mm, ratio 0.19) groups (see figure 1). Conclusions PA size on CT is associated with mPAP in subjects with chronic lung disease with and without PH. This association is strongest in subjects without parenchymal lung disease but with features of hypercarbia or hypoxia alone. These data suggest that features of parenchymal lung disease in certain patients may limit expected PA dilation in response to elevated PA pressures. Further study is needed regarding the implications of these findings for early identification of PH in patients with chronic lung disease. This abstract is funded by: The NHLBI and Pulmonary Hypertension Association

Abstract Rationale Pulmonary hypertension (PH) is a heterogeneous condition classified into 5 WSPH groups but lacking detailed molecular distinction. The PVDOMICS cohort includes all PH groups alongside healthy and disease comparators, is deeply clinically phenotyped and has RNA sequencing data. Pseudogenes are conserved genome elements that due to sequence-defects cannot be translated into functional proteins. Some pseudogenes never-the-less produce RNA-transcripts which can regulate the expression of other, protein coding genes. Therefore, our study was focused on identifying differentially expressed pseudogene-transcripts within the dataset from the PVDOMICS study and gaining insights into their potential function. Methods Peripheral venous whole blood was collected from PVDOMICS participants [n = 525 PH patients; n = 312 controls without PH who were comparators (n = 256) or healthy (n = 56)]. RNA was isolated, RNA sequencing was performed. After quality control, differential transcript expression between PH and control groups was analyzed using the limma-voom pipeline, rSEA with adjustment for age, sex, race, smoking, WSPH group, latent factors and false discovery rate. We focused on a) transcripts expressed at increased levels (n = 218): the first 200 DEG (differentially-expressed-genes)-transcripts based on log2 difference plus increased DEG within the first 200 most significant p-values; and b) transcripts expressed at decreased levels (n = 292): last 200 DEG based on log2 difference plus decreased DEG within the first 200 most significant p-values. Results Six pseudogene transcripts were among the increased and 35 pseudogene transcripts were among the decreased DEG in the PH group compared to controls. Some of the differentially expressed pseudogenes (DEP) belonged to gene families, like the ribosomal protein family (7 DEP), mitochondria expressed transcripts (3 DEP), zinc finger proteins (2 DEP). Associations of pseudogenes with biological function or disease were extracted via GeneCards(R) and information from GWAS analysis. The DUSP8P5 pseudogene was among the DEP and is known for its association with left heart disease and hypertension; 2 additional DEP have known association with hypertension. Further known associations were with inhibitory or inflammatory immune processes (8 DEP), bone density (6 DEP), metabolism and diabetes (6 DEP), respiratory function or disease (4 DEP), erythrocyte biology (3 DEP), cancer (5 DEP), kidney disease (1 DEP), and demographic parameters like eye or hair color (2 DEP). Ten DEP had no known association with function or disease. Conclusions We identified pseudogene RNA-transcripts associated with PH across WSPH groups. The local vs. systemic effects of these pseudogene transcripts on PH remain to be determined and regulators of pseudogene-transcript expression (e.g. environmental exposures) identified. This abstract is funded by: 1 X01 HL148611 (TopMed), PVDOMICS NHLBI, PVDOMICS Pulmonary Hypertension Association, HT9425-25-1-0533 (GG), 1R21OH012955-01 (GG)

Abstract Background Prostaglandin E2 (PGE2) exerts paradoxical effects in asthma, acting as both a protective mediator and a contributor to pathology depending on context. Its role is shaped by type 2 (T2)-status, eosinophilia, and exacerbation frequency. In epithelial cell cultures from asthmatic and healthy individuals, interleukin (IL)-13 suppresses PGE2 synthesis genes (PTGS2 and PGES) while upregulating HPGD, which encodes for the PGE2 catabolic enzyme 15-prostaglandin dehydrogenase (15-PGDH) (PMID:16751012), suggesting diminished PGE2 availability during T2 inflammation. Similarly, in aspirin-exacerbated respiratory disease (AERD), impaired cyclooxygenase disrupts PGE2 biosynthesis, fueling unchecked leukotriene production and eosinophilic inflammation (PMID: 30516547), underscoring PGE2’s protective role in limiting disease progression. Conversely, in nasal polyps, tuft cells are activated by IL-13 to produce PGE2 (PMID: 35608904), highlighting PGE2’s potential pathogenicity. Objective To assess associations between PGE2-pathway gene expression in bronchial brushings and asthma severity, and to determine how nasal polyp history modifies these relationships. Methods We analyzed bronchial epithelial brushings from 109 adults with asthma (mean age 42.4±12.1, 64% female, 52% severe asthma, 18% nasal polyps) using normalized and batch-corrected RNA-seq data from the baseline visit in the Severe Asthma Research Program (SARP)3 (PMID: 36796454). Age-adjusted and nasal polyp-adjusted relationships were calculated between ratios of PGE2 synthesis and metabolism genes and clinical outcomes (ASSESS score, FEV1 %predicted). The composite PGE2 score (PMID: 35608904), derived from nasal polyp analysis, was calculated from the average scaled log-normalized expression of six PGE2-responsive genes. All data were analyzed using a linear regression mixed model with a random effect for site. Results Ratios of PGE2 synthesis genes (PTGS1/2, PTGES2) to HPGD were significantly associated with greater disease severity (change in ASSESS Score) (Fig.1A, solid lines) and lower lung function (change in FEV1 %predicted) (Fig.1B, solid lines). These relationships persisted after adjusting for nasal polyps (Fig.1, dashed lines). The composite PGE2 score was also related to disease severity (Fig.1A, black solid line) but attenuated when nasal polyps were accounted for (Fig.1A, black dashed lines). Interestingly, lower lung function remained associated with the PGE2 score even after nasal polyp adjustment (Fig.1B, dashed black lines). Conclusions Our analysis demonstrates that the balance between PGE2 synthesis and degradation in the bronchial epithelium is closely linked to asthma severity and lung function, independent of nasal polyp status. Future work will integrate transcriptomic and clinical data longitudinally to establish whether PGE2 could serve as a biomarker for asthma trajectory. This abstract is funded by: HL115118 (LCD), U10HL109152 (S.E.W.), U10HL109172 (to E.I. and B.D.L), U10HL109168 (N.N.J.), U10HL109250 (S.C.E), U10HL109164 (W.C.M), U10HL109257 (M.C), and U10HL109146 (to P.G.W. and S.A.C.)

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

Abstract Rationale The activity of fucosyltransferase 2 (FUT2) determines if mucin glycan structures in epithelial mucus include histo-blood group antigens (HBGAs). People with active FUT2 express HBGAs in epithelial mucus secretions and are terms “secretors”, while non-secretors lack HBGAs in their secretions because they have a FUT2 genotype that inactivates the glycosyl transferase. Prior research suggests that secretors who are blood type O (type O secretors) are susceptible to asthma exacerbations1 and that FUT2 knockout mice (functional non-secretors) are protected from mucus plugging in allergen challenge models2,3. We set out here to determine if secretor status influences mucus plug prevalence in asthma and to confirm if type O secretors are susceptible to asthma exacerbations. Methods We determined secretor status in asthma participants from the Severe Asthma Research Program-3 (SARP-3) cohort (N = 507) in two ways, a specific fucose lectin binding assay in sputum and FUT2 genotyping. Blood group type was determined using ABO genotyping. We analyzed mucus plug segment scores, asthma exacerbations rates and other clinical and biological features in patients stratified by secretor status and blood type (non-secretors, type O secretors, and non-type O secretors). Results Secretor status determined by lectin binding in sputum was strongly concordant with genotypic analysis. Within the SARP-3 cohort, 26.2% were non-secretors, 32.3% type were O secretors, and 41.4% were non-type O secretors. Mucus plug segment scores, annualized asthma exacerbation rates, and other measures of asthma severity were similar in non-secretors, type O secretors, and non-type O secretors (Table). Although type O secretors had increased values for fractional exhaled nitric oxide (FeNO) and higher expression of type 2 gene network genes in their sputum cells, measures of sputum eosinophil inflammation were not higher in this subgroup (Table). Conclusion Secretor status does not influence the prevalence of mucus plugging in asthma patients, and we did not confirm a prior report that type O secretors have an increased risk of asthma exacerbation. Although type O secretors have modestly increased FeNO levels and increased expression of type 2 immune genes in their sputum, it is not clear if this is a biologically important signal. 1. Innes et al. AJRCCM; 2011. 2. Saku et al. JACI; 2019. 3. Hara et al. AJRCMB; 2025 This abstract is funded by: The study was supported by the National Institutes of Health grants R01 HL080414 (JVF), P01 HL107202 (JVF), R01 HL164787 (JVF), U10 HL109172 (EI), U10 HL109152 (SEW), U10 HL109257 (MC), U10 HL109250 (SCE), U10 HL109146 (JVF), U10 HL109164 (ERB), U01 HL146002 (BDL), U10 HL109086 (DTM) and Sandler Asthma Basic Research Center at UCSF (JVF). The following companies provided financial support for study activities at the Coordinating and Clinical Centers beyond the third year of patient follow-up: AstraZeneca, Boehringer-Ingelheim, Genentech, GlaxoSmithKline, Sanofi–Genzyme Regeneron, and Teva. These companies had no role in study design or data analysis, and the only restriction on the funds was that they be used to support the SARP initiative.

Abstract Rationale We hypothesize that corticosteroid resistance causes uncontrolled inflammation and severe asthma. The genetic and molecular mechanisms underlying steroid responsiveness in asthma remain poorly understood (PMID28493293). This study aims to characterize transcriptomic profiles distinguishing steroid responders from steroid resistant asthma. Methods Longitudinal bulk RNA-seq of SARP3 sputum samples was analyzed before and after an “induced phenotype” with 40 mg triamcinolone acetonide (TCA) injection (PMID 27967215). Participants were classified as responders or non-responders based on whether their FEV₁ change after TCA exceeded the median. Paired differential gene expression analyses were performed separately for each group using limma and voom, with participant as a factor and adjustment for sputum cell count differentials. Z-scores were calculated to compare gene expression changes after TCA injection between responders and non-responders, and pathway enrichment was assessed using the fgsea package. Results Baseline demographics and clinical characteristics of 696 SARP3 asthma participants receiving TCA are shown in Table 1. The median [IQR] FEV1 change post-TCA was 1.81% [-2.16, 7.17]. Responders had higher baseline FEV1%, FEV1/FVC ratio, blood eosinophil counts, and sputum eosinophil percentages, but lower serum CRP levels. Steroid responsiveness was associated with asthma severity (Table 1). Among 154 participants with paired sputum RNA-seq data, 70 (45%) exhibited >1.81% improvement in FEV₁ and were classified as responders. Following TCA treatment, non-responders showed no significant changes in gene expression. In responders, 18 genes were upregulated and 203 downregulated. Downregulated genes were associated with Th2 immunity (IL13, IL1RL1, CCL17, SIGLEC8), IgE signaling (FCER1A, FCER2), epithelial barrier function and mucus production (MUC2, MUC7, ALOX15, TFF3), tissue remodeling and fibrosis (MMP10, TIMP3), ciliary function (DNAH5, CFAP43, FOXJ1), and airway hyperreactivity (HRH4, RAMP1). Key genes notably downregulated in responders compared to non-responders included IL13, CCL5, PTGDR2, FCER2, and CLEC10A. Pathway analysis revealed enhanced IL-10-mediated anti-inflammatory signaling, suppression of pro-inflammatory pathways (NF-κB, IL-1, Th1/Th2, asthma, arachidonic acid metabolism, chemokine, NOTCH4), and regulation of mitochondrial oxidative phosphorylation. Conclusions Steroid resistance impairs glucocorticoid-mediated suppression of pro-inflammatory pathways, including NF-κB, IL-1, and Th1/Th2 signaling. This persistent airway inflammation limits improvements in lung function and contributes to worse clinical outcomes, including increased airflow obstruction, and symptom burden. Glucocorticoid resistance in asthma may be driven by receptor dysfunction, airway inflammation, epigenetic changes, and environmental exposures resulting in persistent airway inflammation and remodeling. Elucidating the molecular mechanisms underlying steroid resistance is critical for developing personalized targeted therapies in asthma. This abstract is funded by: NIH-NHLBI - R01 HL161674 (PI: Joe Zein)