Abstract
To elucidate the dynamic evolution of cancer cell characteristics within the tumor microenvironment (TME), we developed an integrative approach combining single-cell tracking, cell fate simulation, and 3D TME modeling. We began our investigation by analyzing the spatiotemporal behavior of individual cancer cells in cultured pancreatic (MiaPaCa2) and cervical (HeLa) cancer cell lines, with a focus on the α2-6 sialic acid (α2-6Sia) modification on glycans, which is associated with cell stemness. Our findings revealed that MiaPaCa2 cells exhibited significantly higher levels of α2-6Sia modification, correlating with enhanced reproductive capabilities, whereas HeLa cells showed less prevalence of this modification. To accommodate the in vivo variability of α2-6Sia levels, we employed a cell fate simulation algorithm that digitally generates cell populations based on our observed data while varying the level of sialylation, thereby simulating cell growth patterns. Subsequently, we performed a 3D TME simulation with these deduced cell populations, considering the microenvironment that could impact cancer cell growth. Immune cell landscape information derived from 193 cervical and 172 pancreatic cancer cases was used to estimate the degree of the positive or negative impact. Our analysis suggests that the deduced cells generated based on the characteristics of MiaPaCa2 cells are less influenced by the immune cell landscape within the TME compared to those of HeLa cells, highlighting that the fate of cancer cells is shaped by both the surrounding immune landscape and the intrinsic characteristics of the cancer cells.