Mathematical modelling of phenotypic plasticity and conversion to a stem-cell state under hypoxia

Hypoxia, or oxygen deficiency, is known to be associated with breast tumour progression, resistance to conventional therapies and poor clinical prognosis. The epithelial-mesenchymal transition (EMT) is a process that confers invasive and migratory capabilities as well as stem cell properties to carcinoma cells thus promoting metastatic progression. In this work, we examined the impact of hypoxia on EMT-associated cancer stem cell (CSC) properties, by culturing transformed human mammary epithelial cells under normoxic and hypoxic conditions, and applying in silico mathematical modelling to simulate the impact of hypoxia on the acquisition of CSC attributes and the transitions between differentiated and stem-like states. Our results indicate that both the heterogeneity and the plasticity of the transformed cell population are enhanced by exposure to hypoxia, resulting in a shift towards a more stem-like population with increased EMT features. Our findings are further reinforced by gene expression analyses demonstrating the upregulation of EMT-related genes, as well as genes associated with therapy resistance, in hypoxic cells compared to normoxic counterparts. In conclusion, we demonstrate that mathematical modelling can be used to simulate the role of hypoxia as a key contributor to the plasticity and heterogeneity of transformed human mammary epithelial cells.