RESUMEN
BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is a dismal disease with a survival rate of less than 7%, mainly due to the hepatic metastatic spread. Despite the importance of understanding PDAC metastases, central questions remain concerning their biology and chemosensitivity. Moreover, the transcriptomic divergence between primary tumor (PT) and hepatic metastases (HM) has been poorly studied and without a clear dissection of the confounding tumoral-surrounding tissue. METHODS: Here, to unravel key biological features not biased by the surrounding tissue, we implemented a blind source separation based on independent component analysis, ProDenICA, on a treatment-naïve cohort of PDAC paired samples and a cohort of 305 resectable patients. In addition, a time-lapse experiment was performed to assess the gemcitabine chemosensitivity profile between the PT and HM. RESULTS: We identified HM's specific transcriptomic characteristics related to the upregulation of cell cycle checkpoint, mitochondria activity, and extracellular matrix reorganization, which could be associated with metastatic niche adaptation mechanisms. Furthermore, squamous lineage emerged as a key feature linked with a downregulation in the epithelial-to-mesenchymal program that can stratifies PDAC HM independent of the classical/basal-like spectrum. Remarkably, we also demonstrated that gemcitabine response is influenced by the squamous profile, being the HM more refractory to the treatment than the PT. CONCLUSIONS: These results pointed out divergent HM aspects compared to PT and allowed their stratification through the squamous lineage. Moreover, we unravel a clinical actionable squamous signature that predicts the gemcitabine response.
