Identifying Molecular Signatures of Distinct Modes of Collective Migration in Response to the Microenvironment Using Three-Dimensional Breast Cancer Models.

Collective cell migration is a key feature of transition of ductal carcinoma in situ (DCIS) to invasive ductal carcinoma (IDC) among many other cancers, yet the microenvironmental factors and underlying mechanisms that trigger collective migration remain poorly understood. Here, we investigated two microenvironmental factors, tumor-intrinsic hypoxia and tumor-secreted factors (secretome), as triggers of collective migration using three-dimensional (3D) discrete-sized microtumor models that recapitulate hallmarks of DCIS-IDC transition. Interestingly, the two factors induced two distinct modes of collective migration: directional and radial migration in the 3D microtumors generated from the same breast cancer cell line model, T47D. Without external stimulus, large (600 µm) T47D microtumors exhibited tumor-intrinsic hypoxia and directional migration, while small (150 µm), non-hypoxic microtumors exhibited radial migration only when exposed to the secretome of large microtumors. To investigate the mechanisms underlying hypoxia- and secretome-induced directional vs. radial migration modes, we performed differential gene expression analysis of hypoxia- and secretome-induced migratory microtumors compared with non-hypoxic, non-migratory small microtumors as controls. We propose unique gene signature sets related to tumor-intrinsic hypoxia, hypoxia-induced epithelial-mesenchymal transition (EMT), as well as hypoxia-induced directional migration and secretome-induced radial migration. Gene Set Enrichment Analysis (GSEA) and protein-protein interaction (PPI) network analysis revealed enrichment and potential interaction between hypoxia, EMT, and migration gene signatures for the hypoxia-induced directional migration. In contrast, hypoxia and EMT were not enriched in the secretome-induced radial migration, suggesting that complete EMT may not be required for radial migration. Survival analysis identified unique genes associated with low survival rate and poor prognosis in TCGA-breast invasive carcinoma dataset from our tumor-intrinsic hypoxia gene signature (CXCR4, FOXO3, LDH, NDRG1), hypoxia-induced EMT gene signature (EFEMP2, MGP), and directional migration gene signature (MAP3K3, PI3K3R3). NOS3 was common between hypoxia and migration gene signature. Survival analysis from secretome-induced radial migration identified ATM, KCNMA1 (hypoxia gene signature), and KLF4, IFITM1, EFNA1, TGFBR1 (migration gene signature) to be associated with poor survival rate. In conclusion, our unique 3D cultures with controlled microenvironments respond to different microenvironmental factors, tumor-intrinsic hypoxia, and secretome by adopting distinct collective migration modes and their gene expression analysis highlights the phenotypic heterogeneity and plasticity of epithelial cancer cells.