AUTHOR=Lotsberg Maria L. , Røsland Gro V. , Rayford Austin J. , Dyrstad Sissel E. , Ekanger Camilla T. , Lu Ning , Frantz Kirstine , Stuhr Linda E. B. , Ditzel Henrik J. , Thiery Jean Paul , Akslen Lars A. , Lorens James B. , Engelsen Agnete S. T. 

TITLE=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

JOURNAL=Frontiers in Oncology

VOLUME=Volume 12 - 2022

YEAR=2022

URL=https://www.frontiersin.org/journals/oncology/articles/10.3389/fonc.2022.818437

DOI=10.3389/fonc.2022.818437

ISSN=2234-943X

ABSTRACT=The lack of adequate preclinical models remains a limitation for cancer drug development and is a primary contributor to anti-cancer drug failures in clinical trials. 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 proof-of-principle 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 homotypic and 3D 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 imaging mass cytometry (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.