Integrative multi-omics and single-cell analysis identifies EGFR pathway activation and metabolic reprogramming as...

Although fibroblast growth factor receptor (FGFR) inhibitors (FGFRi) have demonstrated clinical promise, the inevitable emergence of acquired resistance remains a critical bottleneck, severely compromising their long-term clinical efficacy. The pan-cancer molecular landscape and heterogeneous mechanisms driving this resistance, ranging from genetic alterations to dynamic network rewiring, remain poorly understood.

We integrated large-scale pharmacogenomic profiling of the FGFR inhibitor AZD4547 from the GDSC2 and PRISM databases with single-cell RNA sequencing to dissect the multi-omics landscape of FGFRi resistance across 312 cell lines from 8 cancer types. This multi-omics framework was further extended by machine learning modeling and systematic synthetic lethality screening to uncover actionable therapeutic targets. In vitro viability assays and western blot analysis were subsequently conducted to experimentally evaluate the predicted FGFR-EGFR synthetic lethality.

Our dual-database analysis unveiled a multi-dimensional atlas of FGFRi resistance. We identified cancer-specific genomic drivers, such as ELF4 amplification in glioblastoma, alongside key transcriptomic markers including UCP2 and FSCN1, highlighting a shift towards metabolic reprogramming and epithelial-mesenchymal transition (EMT). Single-cell analysis unveiled that resistance is linked to the heterogeneous enrichment of baseline subpopulations characterized by distinct metaprograms, including cell-cycle dysregulation. Furthermore, a random forest model built on a LASSO-derived transcriptomic signature was constructed, demonstrating promising predictive capability for AZD4547 sensitivity (mean test-set AUC = 0.73, 95% CI [0.63, 0.80]); the signature generalized well to erdafitinib but showed limited transferability to some other FGFR inhibitors (e.g. pemigatinib, BGJ398). Most notably, our synthetic lethal screening revealed a convergent reliance on compensatory RTK signaling (specifically EGFR pathway enrichment) and downstream MAPK/PI3K cascades in resistant phenotypes, providing converging computational evidence for EGFR pathway activation as an adaptive bypass mechanism. This predicted synthetic lethality was experimentally supported in two FGFR-dependent cell line models (RT112 and CCLP1), in which combined FGFR-EGFR inhibition produced marked synergistic antiproliferative effects.

This study establishes a comprehensive multi-omics atlas of resistance to the FGFR inhibitor AZD4547, delineating convergent mechanisms of metabolic reprogramming and EGFR-mediated bypass signaling. Our findings characterize the resistance as a dynamic network rewiring and nominate rational combination strategies to overcome this therapeutic bottleneck. While FGFR-EGFR co-inhibition is experimentally supported, metabolic co-targeting remains a computationally derived, hypothesis-generating strategy.