Intrinsic Differences in Spatiotemporal Organization and Stromal Cell Interactions Between Isogenic Lung Cancer Cells of Epithelial and Mesenchymal Phenotypes Revealed by High-Dimensional Single-Cell Analysis of Heterotypic 3D Spheroid Models.

The lack of inadequate preclinical models remains a limitation for cancer drug development and is a primary contributor to anti-cancer drug failures in clinical trials. Heterotypic multicellular spheroids are three-dimensional (3D) spherical structures generated by self-assembly from aggregates of two or more cell types. Compared to traditional monolayer cell culture models, the organization of cells into a 3D tissue-like structure favors relevant physiological conditions with chemical and physical gradients as well as cell-cell and cell-extracellular matrix (ECM) interactions that recapitulate many of the hallmarks of cancer . Epidermal growth factor receptor (EGFR) mutations are prevalent in non-small cell lung cancer (NSCLC), yet various mechanisms of acquired resistance, including epithelial-to-mesenchymal transition (EMT), limit the clinical benefit of EGFR tyrosine kinase inhibitors (EGFRi). Improved preclinical models that incorporate the complexity induced by epithelial-to-mesenchymal plasticity (EMP) are urgently needed to advance new therapeutics for clinical NSCLC management. This study was designed to provide a thorough characterization of multicellular spheroids of isogenic cancer cells of various phenotypes and demonstrate for the applicability of the presented spheroid model to evaluate the impact of cancer cell phenotype in drug screening experiments through high-dimensional and spatially resolved imaging mass cytometry (IMC) analyses. First, we developed and characterized 3D homotypic and heterotypic spheroid models comprising EGFRi-sensitive or EGFRi-resistant NSCLC cells. We observed that the degree of EMT correlated with the spheroid generation efficiency in monocultures. In-depth characterization of the multicellular heterotypic spheroids using immunohistochemistry and high-dimensional single-cell analyses by IMC revealed intrinsic differences between epithelial and mesenchymal-like cancer cells with respect to self-sorting, spatiotemporal organization, and stromal cell interactions when co-cultured with fibroblasts. While the carcinoma cells harboring an epithelial phenotype self-organized into a barrier sheet surrounding the fibroblasts, mesenchymal-like carcinoma cells localized to the central hypoxic and collagen-rich areas of the compact heterotypic spheroids. Further, deep-learning-based single-cell segmentation of IMC images and application of dimensionality reduction algorithms allowed a detailed visualization and multiparametric analysis of marker expression across the different cell subsets. We observed a high level of heterogeneity in the expression of EMT markers in both the carcinoma cell populations and the fibroblasts. Our study supports further application of these models in pre-clinical drug testing combined with complementary high-dimensional single-cell analyses, which in turn can advance our understanding of the impact of cancer-stroma interactions and epithelial phenotypic plasticity on innate and acquired therapy resistance in NSCLC.