Single-cell analysis of an engineered organoid-based model of pancreatic cancer identifies hypoxia as a contributing factor in the determination of transcriptional subtypes

Abstract

Pancreatic ductal adenocarcinoma (PDAC) is a high-mortality cancer characterized by its aggressive, treatment-resistant phenotype and a complex tumour microenvironment (TME) featuring significant hypoxia. Bulk transcriptomic analysis has identified the classical and basal-like transcriptional subtypes which have prognostic value in PDAC; however, it remains unclear how microenvironmental heterogeneity contributes to the expression of these transcriptional signatures. Here, we used single cell transcriptome analysis of the organoid TRACER platform to explore the effect of oxygen and other microenvironmental gradients on PDAC organoid cells. We found that the microenvironmental gradients present in TRACER significantly impact the distribution of organoid transcriptional phenotypes and the enrichment of gene sets linked to cancer progression and treatment resistance. More significantly, we found that microenvironmental gradients drive changes in the expression of the classical and basal-like transcriptional subtype gene signatures. This effect is likely dominated by the oxygen gradients in TRACER, as hypoxia alone induced decreases in the expression of classical marker GATA6 at both the gene and protein level in PDAC cells. This work suggests that hypoxia contributes to determining transcriptional subtypes in PDAC and broadly underscores the importance of considering microenvironmental gradients in organoid-based transcriptomic studies of PDAC.