Investigating epithelial‐mesenchymal heterogeneity of tumors and circulating tumor cells with transcriptomic analysis and biophysical modeling

Feasibility of mass cytometry proteomic characterisation of circulating tumour cells in head and neck squamous cell carcinoma for deep phenotyping

BACKGROUND

Circulating tumour cells (CTCs) are a potential cancer biomarker, but current methods of CTC analysis at single-cell resolution are limited. Here, we describe high-dimensional single-cell mass cytometry proteomic analysis of CTCs in HNSCC.

METHODS

Parsortix microfluidic-enriched CTCs from 14 treatment-naïve HNSCC patients were analysed by mass cytometry analysis using 41 antibodies. Immune cell lineage, epithelial-mesenchymal transition (EMT), stemness, proliferation and immune checkpoint expression was assessed alongside phosphorylation status of multiple signalling proteins. Patient-matched tumour gene expression and CTC EMT profiles were compared. Standard bulk CTC RNAseq was performed as a baseline comparator to assess mass cytometry data.

RESULTS

CTCs were detected in 13/14 patients with CTC counts of 2-24 CTCs/ml blood. Unsupervised clustering separated CTCs into epithelial, early EMT and advanced EMT groups that differed in signalling pathway activation state. Patient-specific CTC cluster patterns separated into immune checkpoint low and high groups. Patient tumour and CTC EMT profiles differed. Mass cytometry outperformed bulk RNAseq to detect CTCs and characterise cell phenotype.

DISCUSSION

We demonstrate mass cytometry allows high-plex proteomic characterisation of CTCs at single-cell resolution and identify common CTC sub-groups with potential for novel biomarker development and immune checkpoint inhibitor treatment stratification.

Clinical Relevancy of Circulating Tumor Cells in Breast Cancer: Epithelial or Mesenchymal Characteristics, Single Cells or Clusters?

Metastatic breast cancer (MBC) is typically an incurable disease with high mortality rates; thus, early identification of metastatic features and disease recurrence through precise biomarkers is crucial. Circulating tumor cells (CTCs) consisting of heterogeneous subpopulations with different morphology and genetic, epigenetic, and gene expression profiles represent promising candidate biomarkers for metastatic potential. The experimentally verified role of epithelial-to-mesenchymal transition in cancer dissemination has not been clearly described in BC patients, but the stemness features of CTCs strongly contributes to metastatic potency. Single CTCs have been shown to be protected in the bloodstream against recognition by the immune system through impaired interactions with T lymphocytes and NK cells, while associations of heterotypic CTC clusters with platelets, leucocytes, neutrophils, tumor-associated macrophages, and fibroblasts improve their tumorigenic behavior. In addition to single CTC and CTC cluster characteristics, we reviewed CTC evaluation methods and clinical studies in early and metastatic BCs. The variable CTC tests were developed based on specific principles and strategies. However, CTC count and the presence of CTC clusters were shown to be most clinically relevant in existing clinical trials. Despite the known progress in CTC research and sampling of BC patients, implementation of CTCs and CTC clusters in routine diagnostic and treatment strategies still requires improvement in detection sensitivity and precise molecular characterizations, focused predominantly on the role of CTC clusters for their higher metastatic potency.

Transcriptomic-Based Quantification of the Epithelial-Hybrid-Mesenchymal Spectrum across Biological Contexts

Epithelial-mesenchymal plasticity (EMP) underlies embryonic development, wound healing, and cancer metastasis and fibrosis. Cancer cells exhibiting EMP often have more aggressive behavior, characterized by drug resistance, and tumor-initiating and immuno-evasive traits. Thus, the EMP status of cancer cells can be a critical indicator of patient prognosis. Here, we compare three distinct transcriptomic-based metrics-each derived using a different gene list and algorithm-that quantify the EMP spectrum. Our results for over 80 cancer-related RNA-seq datasets reveal a high degree of concordance among these metrics in quantifying the extent of EMP. Moreover, each metric, despite being trained on cancer expression profiles, recapitulates the expected changes in EMP scores for non-cancer contexts such as lung fibrosis and cellular reprogramming into induced pluripotent stem cells. Thus, we offer a scoring platform to quantify the extent of EMP in vitro and in vivo for diverse biological applications including cancer.

Single-Cell RNA-Seq Analysis Reveals the Acquisition of Cancer Stem Cell Traits and Increase of Cell-Cell Signaling during EMT Progression

Intermediate cell states (ICSs) during the epithelial-mesenchymal transition (EMT) are emerging as a driving force of cancer invasion and metastasis. ICSs typically exhibit hybrid epithelial/mesenchymal characteristics as well as cancer stem cell (CSC) traits including proliferation and drug resistance. Here, we analyze several single-cell RNA-seq (scRNA-seq) datasets to investigate the relation between several axes of cancer progression including EMT, CSC traits, and cell-cell signaling. To accomplish this task, we integrate computational methods for clustering and trajectory inference with analysis of EMT gene signatures, CSC markers, and cell-cell signaling pathways, and highlight conserved and specific processes across the datasets. Our analysis reveals that "standard" measures of pluripotency often used in developmental contexts do not necessarily correlate with EMT progression and expression of CSC-related markers. Conversely, an EMT circuit energy that quantifies the co-expression of epithelial and mesenchymal genes consistently increases along EMT trajectories across different cancer types and anatomical locations. Moreover, despite the high context specificity of signal transduction across different cell types, cells undergoing EMT always increased their potential to send and receive signals from other cells.

Cluster size distribution of cells disseminating from a primary tumor

The first stage of the metastatic cascade often involves motile cells emerging from a primary tumor either as single cells or as clusters. These cells enter the circulation, transit to other parts of the body and finally are responsible for growth of secondary tumors in distant organs. The mode of dissemination is believed to depend on the EMT nature (epithelial, hybrid or mesenchymal) of the cells. Here, we calculate the cluster size distribution of these migrating cells, using a mechanistic computational model, in presence of different degree of EMT-ness of the cells; EMT is treated as given rise to changes in their active motile forces (μ) and cell-medium surface tension (Γ). We find that, for (μ > μ_min, Γ > 1), when the cells are hybrid in nature, the mean cluster size, N¯∼Γ2.0/μ2.8, where μ_min increases with increase in Γ. For Γ ≤ 0, N¯=1, the cells behave as completely mesenchymal. In presence of spectrum of hybrid states with different degree of EMT-ness (motility) in primary tumor, the cells which are relatively more mesenchymal (higher μ) in nature, form larger clusters, whereas the smaller clusters are relatively more epithelial (lower μ). Moreover, the heterogeneity in μ is comparatively higher for smaller clusters with respect to that for larger clusters. We also observe that more extended cell shapes promote the formation of smaller clusters. Overall, this study establishes a framework which connects the nature and size of migrating clusters disseminating from a primary tumor with the phenotypic composition of the tumor, and can lead to the better understanding of metastasis.

In the process of metastasis, tumor cells disseminate from the primary tumor either as single cells or multicellular clusters. These clusters are potential contributor to the initiation of secondary tumor in distant organs. Our computational model captures the size distribution of migrating clusters depending on the adhesion and motility of the cells (which determine the degree of their EMT nature). Furthermore, we investigate the effect of heterogeneity of cell types in the primary tumor on the resultant heterogeneity of cell types in clusters of different sizes. We believe that the understanding the formation and nature of these clusters, dangerous actors in the deadly aspect of cancer progression, will be useful for improving prognostic methods and eventually better treatments.

Systems-level network modeling deciphers the master regulators of phenotypic plasticity and heterogeneity in melanoma

Phenotypic (i.e. non-genetic) heterogeneity in melanoma drives dedifferentiation, recalcitrance to targeted therapy and immunotherapy, and consequent tumor relapse and metastasis. Various markers or regulators associated with distinct phenotypes in melanoma have been identified, but, how does a network of interactions among these regulators give rise to multiple "attractor" states and phenotypic switching remains elusive. Here, we inferred a network of transcription factors (TFs) that act as master regulators for gene signatures of diverse cell-states in melanoma. Dynamical simulations of this network predicted how this network can settle into different "attractors" (TF expression patterns), suggesting that TF network dynamics drives the emergence of phenotypic heterogeneity. These simulations can recapitulate major phenotypes observed in melanoma and explain de-differentiation trajectory observed upon BRAF inhibition. Our systems-level modeling framework offers a platform to understand trajectories of phenotypic transitions in the landscape of a regulatory TF network and identify novel therapeutic strategies targeting melanoma plasticity.

Metastasis Prevention: Focus on Metastatic Circulating Tumor Cells

Metastasis is the main cause of cancer death. Metastatic foci are derived from tumor cells that detach from the primary tumor and then enter the circulation. Circulating tumor cells (CTCs) are generally associated with a high probability of distant metastasis and a negative prognosis. Most CTCs die in the bloodstream, and only a few cells form metastases. Such metastatic CTCs have a stem-like and hybrid epithelial-mesenchymal phenotype, can avoid immune surveillance, and show increased therapy resistance. Targeting metastatic CTCs and their progenitors in primary tumors and their descendants, particularly disseminated tumor cells, represents an attractive strategy for metastasis prevention. However, current therapeutic strategies mainly target the primary tumor and only indirectly affect metastasis-initiating cells. Here, we consider potential methods for preventing metastasis based on targeting molecular and cellular features of metastatic CTCs, including CTC clusters. Also, we emphasize current knowledge gaps in CTC biology that should be addressed to develop highly effective therapeutics and strategies for metastasis suppression.

Phenotypic Heterogeneity of Triple-Negative Breast Cancer Mediated by Epithelial-Mesenchymal Plasticity

Triple-negative breast cancer (TNBC) is a subtype of breast carcinoma known for its unusually aggressive behavior and poor clinical outcome. Besides the lack of molecular targets for therapy and profound intratumoral heterogeneity, the relatively quick overt metastatic spread remains a major obstacle in effective clinical management. The metastatic colonization of distant sites by primary tumor cells is affected by the microenvironment, epigenetic state of particular subclones, and numerous other factors. One of the most prominent processes contributing to the intratumoral heterogeneity is an epithelial-mesenchymal transition (EMT), an evolutionarily conserved developmental program frequently hijacked by tumor cells, strengthening their motile and invasive features. In response to various intrinsic and extrinsic stimuli, malignant cells can revert the EMT state through the mesenchymal-epithelial transition (MET), a process that is believed to be critical for the establishment of macrometastasis at secondary sites. Notably, cancer cells rarely undergo complete EMT and rather exist in a continuum of E/M intermediate states, preserving high levels of plasticity, as demonstrated in primary tumors and, ultimately, in circulating tumor cells, representing a simplified element of the metastatic cascade. In this review, we focus on cellular drivers underlying EMT/MET phenotypic plasticity and its detrimental consequences in the context of TNBC cancer.