Interrogation of Phenotypic Plasticity between Epithelial and Mesenchymal States in Breast Cancer

Mutually exclusive teams-like patterns of gene regulation characterize phenotypic heterogeneity along the noradrenergic-mesenchymal axis in neuroblastoma

Neuroblastoma is the most frequent extracranial pediatric tumor and leads to 15% of all cancer-related deaths in children. Tumor relapse and therapy resistance in neuroblastoma are driven by phenotypic plasticity and heterogeneity between noradrenergic (NOR) and mesenchymal (MES) cell states. Despite the importance of this phenotypic plasticity, the dynamics and molecular patterns associated with these bidirectional cell-state transitions remain relatively poorly understood. Here, we analyze multiple RNA-seq datasets at both bulk and single-cell resolution, to understand the association between NOR- and MES-specific factors. We observed that NOR-specific and MES-specific expression patterns are largely mutually exclusive, exhibiting a "teams-like" behavior among the genes involved, reminiscent of our earlier observations in lung cancer and melanoma. This antagonism between NOR and MES phenotypes was also associated with metabolic reprogramming and with immunotherapy targets PD-L1 and GD2 as well as with experimental perturbations driving the NOR-MES and/or MES-NOR transition. Further, these "teams-like" patterns were seen only among the NOR- and MES-specific genes, but not in housekeeping genes, possibly highlighting a hallmark of network topology enabling cancer cell plasticity.

Cell-state transitions and density-dependent interactions together explain the dynamics of spontaneous epithelial-mesenchymal heterogeneity

Cancer cell populations comprise phenotypes distributed among the epithelial-mesenchymal (E-M) spectrum. However, it remains unclear which population-level processes give rise to the observed experimental distribution and dynamical changes in E-M heterogeneity, including (1) differential growth, (2) cell-state switching, and (3) population density-dependent growth or state-transition rates. Here, we analyze the necessity of these three processes in explaining the dynamics of E-M population distributions as observed in PMC42-LA and HCC38 breast cancer cells. We find that, while cell-state transition is necessary to reproduce experimental observations of dynamical changes in E-M fractions, including density-dependent growth interactions (cooperation or suppression) better explains the data. Further, our models predict that treatment of HCC38 cells with transforming growth factor β (TGF-β) signaling and Janus kinase 2/signal transducer and activator of transcription 3 (JAK2/3) inhibitors enhances the rate of mesenchymal-epithelial transition (MET) instead of lowering that of E-M transition (EMT). Overall, our study identifies the population-level processes shaping the dynamics of spontaneous E-M heterogeneity in breast cancer cells.

Increased prevalence of hybrid epithelial/mesenchymal state and enhanced phenotypic heterogeneity in basal breast cancer

Intra-tumoral phenotypic heterogeneity promotes tumor relapse and therapeutic resistance and remains an unsolved clinical challenge. Decoding the interconnections among different biological axes of plasticity is crucial to understand the molecular origins of phenotypic heterogeneity. Here, we use multi-modal transcriptomic data-bulk, single-cell, and spatial transcriptomics-from breast cancer cell lines and primary tumor samples, to identify associations between epithelial-mesenchymal transition (EMT) and luminal-basal plasticity-two key processes that enable heterogeneity. We show that luminal breast cancer strongly associates with an epithelial cell state, but basal breast cancer is associated with hybrid epithelial/mesenchymal phenotype(s) and higher phenotypic heterogeneity. Mathematical modeling of core underlying gene regulatory networks representative of the crosstalk between the luminal-basal and epithelial-mesenchymal axes elucidate mechanistic underpinnings of the observed associations from transcriptomic data. Our systems-based approach integrating multi-modal data analysis with mechanism-based modeling offers a predictive framework to characterize intra-tumor heterogeneity and identify interventions to restrict it.

Cancer cell plasticity defines response to immunotherapy in cutaneous squamous cell carcinoma

Immune checkpoint blockade (ICB) approaches have changed the therapeutic landscape for many tumor types. However, half of cutaneous squamous cell carcinoma (cSCC) patients remain unresponsive or develop resistance. Here, we show that, during cSCC progression in male mice, cancer cells acquire epithelial/mesenchymal plasticity and change their immune checkpoint (IC) ligand profile according to their features, dictating the IC pathways involved in immune evasion. Epithelial cancer cells, through the PD-1/PD-L1 pathway, and mesenchymal cancer cells, through the CTLA-4/CD80 and TIGIT/CD155 pathways, differentially block antitumor immune responses and determine the response to ICB therapies. Accordingly, the anti-PD-L1/TIGIT combination is the most effective strategy for blocking the growth of cSCCs that contain both epithelial and mesenchymal cancer cells. The expression of E-cadherin/Vimentin/CD80/CD155 proteins in cSCC, HNSCC and melanoma patient samples predicts response to anti-PD-1/PD-L1 therapy. Collectively, our findings indicate that the selection of ICB therapies should take into account the epithelial/mesenchymal features of cancer cells.

Mesenchymal stem cells derived from CHIR99021 and TGF‑β induction remained on the colicomentum and improved cardiac function of a rat model of acute myocardium infarction

Human induced pluripotent stem cells (hiPSCs) have been regarded as a potential stem cell source for cell therapy. However, the production of cells with mesenchymal potential from hiPSCs through spontaneous differentiation is time consuming and laborious. In the present study, the combined use of the GSK-3 inhibitor CHIR99021 and TGF-β was used to obtain mesenchymal stem cell (MSC)-like cells from hiPSCs. During the induction process, the transcription of epithelial-mesenchymal transition (EMT)-related genes N-cadherin and Vimentin in the transformed cells was upregulated, whereas the transcription of E-cadherin and pluripotency-related transcription factors SOX2, OCT4 and NANOG did not change significantly. This indicated that whilst cells were pluripotent, EMT was initiated by the upregulation of transcription of EMT promoting genes. Both SMAD-dependent and independent signalling pathways were significantly activated by the combined induction treatment compared with the single factor induction. The hiPSC-derived MSC-like cells (hiPSC-MSCs) expressed MSC-related markers and acquired osteogenic, chondrogenic and adipogenic differentiation potentials. After being injected into the peritoneal cavity of rats, the hiPSC-MSCs secreted angiogenic and immune-regulatory factors and remained on the colicomentum for 3 weeks. Within an 11-week period, four intraperitoneal hiPSC-MSC injections (1x107 cells/injection) into acute myocardial infarction (AMI) model rats significantly increased the left ventricular ejection fraction, left ventricular fractional shortening and angiogenesis and significantly reduced scar size and the extent of apoptosis in the infarcted area compared with that of the control PBS injection. Symptoms of hiPSC-MSC-induced immune reaction or tumour formation were not observed over the course of the experiment in the hiSPC-MSC treated rats. In conclusion, the CHIR99021 and TGF-β combined induction was a rapid and effective method to obtain MSC-like cells from hiPSCs and multiple high dose intraperitoneal injections of hiPSC-derived MSCs were safe and effective at restoring cardiac function in an AMI rat model.

Docetaxel Inhibits Epithelial-Mesenchymal Transition in Human Mammary Cells

Epithelial-mesenchymal transition (EMT) is a reversible and dynamic biological process in which epithelial cells acquire mesenchymal characteristics including enhanced stemness and migratory ability. EMT can facilitate cancer metastasis and is a known driver of cellular resistance to common chemotherapeutic drugs, such as docetaxel. Current chemotherapeutic practices such as docetaxel treatment can promote EMT and increase the chance of tumor recurrence and resistance, calling for new approaches in cancer treatment. Here we show that prolonged docetaxel treatment at a sub-IC50 concentration inhibits EMT in immortalized human mammary epithelial (HMLE) cells. Using immunofluorescence, flow cytometry, and bulk transcriptomic sequencing to assess EMT progression, we analyzed a range of cellular markers of EMT in docetaxel-treated cells and observed an upregulation of epithelial markers and downregulation of mesenchymal markers in the presence of docetaxel. This finding suggests that docetaxel may have clinical applications not only as a cytotoxic drug but also as an inhibitor of EMT-driven metastasis and multidrug resistance depending on the concentration of its use.

Stochastic epithelial-mesenchymal transitions diversify non-cancerous lung cell behaviours

Epithelial-mesenchymal plasticity (EMP) is a hallmark of cancer. By enabling cells to shift between different morphological and functional states, EMP promotes invasion, metastasis and therapy resistance. We report that near-diploid non-cancerous human epithelial lung cells spontaneously shift along the EMP spectrum without genetic changes. Strikingly, more than half of single cell-derived clones adopt a mesenchymal morphology. We independently characterise epithelial-like and mesenchymal-like clones. Mesenchymal clones lose epithelial markers, display larger cell aspect ratios and lower motility, with mostly unaltered proliferation rates. Stemness marker expression and metabolic rewiring diverge independently of phenotypes. In 3D culture, more epithelial clones become mesenchymal-like. Thus, non-cancerous epithelial cells may acquire cancer metastasis-associated features prior to genetic alterations and cancerous transformation.

Dynamical hallmarks of cancer: Phenotypic switching in melanoma and epithelial-mesenchymal plasticity

Phenotypic plasticity was recently incorporated as a hallmark of cancer. This plasticity can manifest along many interconnected axes, such as stemness and differentiation, drug-sensitive and drug-resistant states, and between epithelial and mesenchymal cell-states. Despite growing acceptance for phenotypic plasticity as a hallmark of cancer, the dynamics of this process remains poorly understood. In particular, the knowledge necessary for a predictive understanding of how individual cancer cells and populations of cells dynamically switch their phenotypes in response to the intensity and/or duration of their current and past environmental stimuli remains far from complete. Here, we present recent investigations of phenotypic plasticity from a systems-level perspective using two exemplars: epithelial-mesenchymal plasticity in carcinomas and phenotypic switching in melanoma. We highlight how an integrated computational-experimental approach has helped unravel insights into specific dynamical hallmarks of phenotypic plasticity in different cancers to address the following questions: a) how many distinct cell-states or phenotypes exist?; b) how reversible are transitions among these cell-states, and what factors control the extent of reversibility?; and c) how might cell-cell communication be able to alter rates of cell-state switching and enable diverse patterns of phenotypic heterogeneity? Understanding these dynamic features of phenotypic plasticity may be a key component in shifting the paradigm of cancer treatment from reactionary to a more predictive, proactive approach.

Multi-modal transcriptomic analysis unravels enrichment of hybrid epithelial/mesenchymal state and enhanced phenotypic heterogeneity in basal breast cancer

Intra-tumoral phenotypic heterogeneity promotes tumor relapse and therapeutic resistance and remains an unsolved clinical challenge. It manifests along multiple phenotypic axes and decoding the interconnections among these different axes is crucial to understand its molecular origins and to develop novel therapeutic strategies to control it. Here, we use multi-modal transcriptomic data analysis - bulk, single-cell and spatial transcriptomics - from breast cancer cell lines and primary tumor samples, to identify associations between epithelial-mesenchymal transition (EMT) and luminal-basal plasticity - two key processes that enable heterogeneity. We show that luminal breast cancer strongly associates with an epithelial cell state, but basal breast cancer is associated with hybrid epithelial/mesenchymal phenotype(s) and higher phenotypic heterogeneity. These patterns were inherent in methylation profiles, suggesting an epigenetic crosstalk between EMT and lineage plasticity in breast cancer. Mathematical modelling of core underlying gene regulatory networks representative of the crosstalk between the luminal-basal and epithelial-mesenchymal axes recapitulate and thus elucidate mechanistic underpinnings of the observed associations from transcriptomic data. Our systems-based approach integrating multi-modal data analysis with mechanism-based modeling offers a predictive framework to characterize intra-tumor heterogeneity and to identify possible interventions to restrict it.

The epithelial-mesenchymal plasticity landscape: principles of design and mechanisms of regulation

Epithelial-mesenchymal plasticity (EMP) enables cells to interconvert between several states across the epithelial-mesenchymal landscape, thereby acquiring hybrid epithelial/mesenchymal phenotypic features. This plasticity is crucial for embryonic development and wound healing, but also underlies the acquisition of several malignant traits during cancer progression. Recent research using systems biology and single-cell profiling methods has provided novel insights into the main forces that shape EMP, which include the microenvironment, lineage specification and cell identity, and the genome. Additionally, key roles have emerged for hysteresis (cell memory) and cellular noise, which can drive stochastic transitions between cell states. Here, we review these forces and the distinct but interwoven layers of regulatory control that stabilize EMP states or facilitate epithelial-mesenchymal transitions (EMTs) and discuss the therapeutic potential of manipulating the EMP landscape.

Statistical inference of the rates of cell proliferation and phenotypic switching in cancer

Recent evidence suggests that nongenetic (epigenetic) mechanisms play an important role at all stages of cancer evolution. In many cancers, these mechanisms have been observed to induce dynamic switching between two or more cell states, which commonly show differential responses to drug treatments. To understand how these cancers evolve over time, and how they respond to treatment, we need to understand the state-dependent rates of cell proliferation and phenotypic switching. In this work, we propose a rigorous statistical framework for estimating these parameters, using data from commonly performed cell line experiments, where phenotypes are sorted and expanded in culture. The framework explicitly models the stochastic dynamics of cell division, cell death and phenotypic switching, and it provides likelihood-based confidence intervals for the model parameters. The input data can be either the fraction of cells or the number of cells in each state at one or more time points. Through a combination of theoretical analysis and numerical simulations, we show that when cell fraction data is used, the rates of switching may be the only parameters that can be estimated accurately. On the other hand, using cell number data enables accurate estimation of the net division rate for each phenotype, and it can even enable estimation of the state-dependent rates of cell division and cell death. We conclude by applying our framework to a publicly available dataset.

Fabrication of a microfluidic device for probiotic drug's dosage screening: Precision Medicine for Breast Cancer Treatment

Breast cancer is the most common cancer in women; it has been affecting the lives of millions each year globally and microfluidic devices seem to be a promising method for the future advancements in this field. This research uses a dynamic cell culture condition in a microfluidic concentration gradient device, helping us to assess breast anticancer activities of probiotic strains against MCF-7 cells. It has been shown that MCF-7 cells could grow and proliferate for at least 24 h; however, a specific concentration of probiotic supernatant could induce more cell death signaling population after 48 h. One of our key findings was that our evaluated optimum dose (7.8 mg/L) was less than the conventional static cell culture treatment dose (12 mg/L). To determine the most effective dose over time and the percentage of apoptosis versus necrosis, flowcytometric assessment was performed. Exposing the MCF-7 cells to probiotic supernatant after 6, 24 and 48 h, confirmed that the apoptotic and necrotic cell death signaling were concentration and time dependent. We have shown a case that these types of microfluidics platforms performing dynamic cell culture could be beneficial in personalized medicine and cancer therapy.

Epigenetic memory acquired during long-term EMT induction governs the recovery to the epithelial state

Epithelial-mesenchymal transition (EMT) and its reverse mesenchymal-epithelial transition (MET) are critical during embryonic development, wound healing and cancer metastasis. While phenotypic changes during short-term EMT induction are reversible, long-term EMT induction has been often associated with irreversibility. Here, we show that phenotypic changes seen in MCF10A cells upon long-term EMT induction by TGFβ need not be irreversible, but have relatively longer time scales of reversibility than those seen in short-term induction. Next, using a phenomenological mathematical model to account for the chromatin-mediated epigenetic silencing of the miR-200 family by ZEB family, we highlight how the epigenetic memory gained during long-term EMT induction can slow the recovery to the epithelial state post-TGFβ withdrawal. Our results suggest that epigenetic modifiers can govern the extent and time scale of EMT reversibility and advise caution against labelling phenotypic changes seen in long-term EMT induction as 'irreversible'.

Network topology metrics explaining enrichment of hybrid epithelial mesenchymal phenotypes in metastasis

Epithelial to Mesenchymal Transition (EMT) and its reverse—Mesenchymal to Epithelial Transition (MET) are hallmarks of metastasis. Cancer cells use this reversible cellular programming to switch among Epithelial (E), Mesenchymal (M), and hybrid Epithelial/Mesenchymal (hybrid E/M) state(s) and seed tumors at distant sites. Hybrid E/M cells are often more aggressive and metastatic than the “pure” E and M cells. Thus, identifying mechanisms to inhibit hybrid E/M cells can be promising in curtailing metastasis. While multiple gene regulatory networks (GRNs) based mathematical models for EMT/MET have been developed recently, identifying topological signatures enriching hybrid E/M phenotypes remains to be done. Here, we investigate the dynamics of 13 different GRNs and report an interesting association between “hybridness” and the number of negative/positive feedback loops across the networks. While networks having more negative feedback loops favor hybrid phenotype(s), networks having more positive feedback loops (PFLs) or many HiLoops–specific combinations of PFLs, support terminal (E and M) phenotypes. We also establish a connection between “hybridness” and network-frustration by showing that hybrid phenotypes likely result from non-reinforcing interactions among network nodes (genes) and therefore tend to be more frustrated (less stable). Our analysis, thus, identifies network topology-based signatures that can give rise to, as well as prevent, the emergence of hybrid E/M phenotype in GRNs underlying EMP. Our results can have implications in terms of targeting specific interactions in GRNs as a potent way to restrict switching to the hybrid E/M phenotype(s) to curtail metastasis.

Stochastic population dynamics of cancer stemness and adaptive response to therapies

Intratumoral heterogeneity can exist along multiple axes: Cancer stem cells (CSCs)/non-CSCs, drug-sensitive/drug-tolerant states, and a spectrum of epithelial-hybrid-mesenchymal phenotypes. Further, these diverse cell-states can switch reversibly among one another, thereby posing a major challenge to therapeutic efficacy. Therefore, understanding the origins of phenotypic plasticity and heterogeneity remains an active area of investigation. While genomic components (mutations, chromosomal instability) driving heterogeneity have been well-studied, recent reports highlight the role of non-genetic mechanisms in enabling both phenotypic plasticity and heterogeneity. Here, we discuss various processes underlying phenotypic plasticity such as stochastic gene expression, chromatin reprogramming, asymmetric cell division and the presence of multiple stable gene expression patterns ('attractors'). These processes can facilitate a dynamically evolving cell population such that a subpopulation of (drug-tolerant) cells can survive lethal drug exposure and recapitulate population heterogeneity on drug withdrawal, leading to relapse. These drug-tolerant cells can be both pre-existing and also induced by the drug itself through cell-state reprogramming. The dynamics of cell-state transitions both in absence and presence of the drug can be quantified through mathematical models. Such a dynamical systems approach to elucidating patterns of intratumoral heterogeneity by integrating longitudinal experimental data with mathematical models can help design effective combinatorial and/or sequential therapies for better clinical outcomes.

Interleukin-3-Receptor-α in Triple-Negative Breast Cancer (TNBC): An Additional Novel Biomarker of TNBC Aggressiveness and a Therapeutic Target

Simple Summary

Molecular and histological profiling is crucial for biomarker and therapeutic target discovery, for example, in TNBC. We demonstrated that IL-3Rα expression led to the identification of a subgroup of TNBC patients displaying a poor overall survival. Moreover, we refined TNBC molecular annotation and drew a model including IL-3Rα, PD-L1, and genes related to EMT, which finely discriminates cancer aggressiveness. Finally, we first demonstrated that IL-3Rα is instrumental in granting tumour adaptation and progression by reprogramming TNBC cells to form large dysfunctional vessels and reshaping PD-L1 expression in primary tumours and metastases. Therefore, the IL-3/IL-3Rα axis may be proposed as a marker of TNBC aggressiveness, as a novel TNBC therapeutic challenge.

Abstract

Tumour molecular annotation is mandatory for biomarker discovery and personalised approaches, particularly in triple-negative breast cancer (TNBC) lacking effective treatment options. In this study, the interleukin-3 receptor α (IL-3Rα) was investigated as a prognostic biomarker and therapeutic target in TNBC. IL-3Rα expression and patients’ clinical and pathological features were retrospectively analysed in 421 TNBC patients. IL-3Rα was expressed in 69% human TNBC samples, and its expression was associated with nodal metastases (p = 0.026) and poor overall survival (hazard ratio = 1.50; 95% CI = 1.01–2.2; p = 0.04). The bioinformatics analysis on the Breast Invasive Carcinoma dataset of The Cancer Genome Atlas (TCGA) proved that IL-3Rα was highly expressed in TNBC compared with luminal breast cancers (p = 0.017, padj = 0.026). Functional studies demonstrated that IL-3Rα activation induced epithelial-to-endothelial and epithelial-to-mesenchymal transition, promoted large blood lacunae and lung metastasis formation, and increased programmed-cell death ligand-1 (PD-L1) in primary tumours and metastases. Based on the TCGA data, IL-3Rα, PD-L1, and EMT coding genes were proposed to discriminate against TNBC aggressiveness (AUC = 0.86 95% CI = 0.82–0.89). Overall, this study identified IL-3Rα as an additional novel biomarker of TNBC aggressiveness and provided the rationale to further investigate its relevance as a therapeutic target.

The Quasimesenchymal Pancreatic Ductal Epithelial Cell Line PANC-1-A Useful Model to Study Clonal Heterogeneity and EMT Subtype Shifting

Simple Summary

Malignant tumors often escape therapy due to clonal propagation of a subfraction of drug-resistant cancer cells. The underlying phenomenon of intratumoral heterogeneity is driven by epithelial–mesenchymal plasticity (EMP) involving the developmental programs, epithelial–mesenchymal transition (EMT), in which epithelial cells are converted to invasive mesenchymal cells, and the reverse process, mesenchymal–epithelial transition (MET), which allows for metastatic outgrowth at distant sites. For therapeutic targeting of EMP, a better understanding of this process is required; however, cellular models with which to study EMP in pancreatic ductal adenocarcinoma (PDAC) are scarce. Using single-cell clonal analysis, we have found the PDAC cell line, PANC-1, to consist of cells with different E/M phenotypes and functional attributes. Parental PANC-1 cultures could be induced in vitro to shift towards either a more mesenchymal or a more epithelial phenotype, and this bidirectional shift was controlled by the small GTPases RAC1 and RAC1b, together identifying PANC-1 cells as a useful model with which to study EMP.

Abstract

Intratumoral heterogeneity (ITH) is an intrinsic feature of malignant tumors that eventually allows a subfraction of resistant cancer cells to clonally evolve and cause therapy failure or relapse. ITH, cellular plasticity and tumor progression are driven by epithelial–mesenchymal transition (EMT) and the reverse process, MET. During these developmental programs, epithelial (E) cells are successively converted to invasive mesenchymal (M) cells, or back to E cells, by passing through a series of intermediate E/M states, a phenomenon termed E–M plasticity (EMP). The induction of MET has clinical potential as it can block the initial EMT stages that favor tumor cell dissemination, while its inhibition can curb metastatic outgrowth at distant sites. In pancreatic ductal adenocarcinoma (PDAC), cellular models with which to study EMP or MET induction are scarce. Here, we have generated single cell-derived clonal cultures of the quasimesenchymal PDAC-derived cell line, PANC-1, and found that these differ strongly with respect to cell morphology and EMT marker expression, allowing for their tentative classification as E, E/M or M. Interestingly, the different EMT phenotypes were found to segregate with differences in tumorigenic potential in vitro, as measured by colony forming and invasive activities, and in circadian clock function. Moreover, the individual clones the phenotypes of which remained stable upon prolonged culture also responded differently to treatment with transforming growth factor (TGF)β1 in regard to regulation of growth and individual TGFβ target genes, and to culture conditions that favour ductal-to-endocrine transdifferentiation as a more direct measure for cellular plasticity. Of note, stimulation with TGFβ1 induced a shift in parental PANC-1 cultures towards a more extreme M and invasive phenotype, while exposing the cells to a combination of the proinflammatory cytokines IFNγ, IL1β and TNFα (IIT) elicited a shift towards a more E and less invasive phenotype resembling a MET-like process. Finally, we show that the actions of TGFβ1 and IIT both converge on regulating the ratio of the small GTPase RAC1 and its splice isoform, RAC1b. Our data provide strong evidence for dynamic EMT–MET transitions and qualify this cell line as a useful model with which to study EMP.

Population Dynamics of Epithelial-Mesenchymal Heterogeneity in Cancer Cells

Phenotypic heterogeneity is a hallmark of aggressive cancer behaviour and a clinical challenge. Despite much characterisation of this heterogeneity at a multi-omics level in many cancers, we have a limited understanding of how this heterogeneity emerges spontaneously in an isogenic cell population. Some longitudinal observations of dynamics in epithelial-mesenchymal heterogeneity, a canonical example of phenotypic heterogeneity, have offered us opportunities to quantify the rates of phenotypic switching that may drive such heterogeneity. Here, we offer a mathematical modeling framework that explains the salient features of population dynamics noted in PMC42-LA cells: (a) predominance of EpCAM^high subpopulation, (b) re-establishment of parental distributions from the EpCAM^high and EpCAM^low subpopulations, and (c) enhanced heterogeneity in clonal populations established from individual cells. Our framework proposes that fluctuations or noise in content duplication and partitioning of SNAIL—an EMT-inducing transcription factor—during cell division can explain spontaneous phenotypic switching and consequent dynamic heterogeneity in PMC42-LA cells observed experimentally at both single-cell and bulk level analysis. Together, we propose that asymmetric cell division can be a potential mechanism for phenotypic heterogeneity.

Deregulation of Transcriptional Enhancers in Cancer

Simple Summary

One of the major challenges in cancer treatments is the dynamic adaptation of tumor cells to cancer therapies. In this regard, tumor cells can modify their response to environmental cues without altering their DNA sequence. This cell plasticity enables cells to undergo morphological and functional changes, for example, during the process of tumour metastasis or when acquiring resistance to cancer therapies. Central to cell plasticity, are the dynamic changes in gene expression that are controlled by a set of molecular switches called enhancers. Enhancers are DNA elements that determine when, where and to what extent genes should be switched on and off. Thus, defects in enhancer function can disrupt the gene expression program and can lead to tumour formation. Here, we review how enhancers control the activity of cancer-associated genes and how defects in these regulatory elements contribute to cell plasticity in cancer. Understanding enhancer (de)regulation can provide new strategies for modulating cell plasticity in tumour cells and can open new research avenues for cancer therapy.

Abstract

Epigenetic regulations can shape a cell’s identity by reversible modifications of the chromatin that ultimately control gene expression in response to internal and external cues. In this review, we first discuss the concept of cell plasticity in cancer, a process that is directly controlled by epigenetic mechanisms, with a particular focus on transcriptional enhancers as the cornerstone of epigenetic regulation. In the second part, we discuss mechanisms of enhancer deregulation in adult stem cells and epithelial-to-mesenchymal transition (EMT), as two paradigms of cell plasticity that are dependent on epigenetic regulation and serve as major sources of tumour heterogeneity. Finally, we review how genetic variations at enhancers and their epigenetic modifiers contribute to tumourigenesis, and we highlight examples of cancer drugs that target epigenetic modifications at enhancers.

Analysis of immune subtypes across the epithelial-mesenchymal plasticity spectrum

Epithelial-mesenchymal plasticity plays a critical role in many solid tumor types as a mediator of metastatic dissemination and treatment resistance. In addition, there is also a growing appreciation that the epithelial/mesenchymal status of a tumor plays a role in immune evasion and immune suppression. A deeper understanding of the immunological features of different tumor types has been facilitated by the availability of large gene expression datasets and the development of methods to deconvolute bulk RNA-Seq data. These resources have generated powerful new ways of characterizing tumors, including classification of immune subtypes based on differential expression of immunological genes. In the present work, we combine scoring algorithms to quantify epithelial-mesenchymal plasticity with immune subtype analysis to understand the relationship between epithelial plasticity and immune subtype across cancers. We find heterogeneity of epithelial-mesenchymal transition (EMT) status both within and between cancer types, with greater heterogeneity in the expression of EMT-related factors than of MET-related factors. We also find that specific immune subtypes have associated EMT scores and differential expression of immune checkpoint markers.

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.

A Wnt-mediated phenotype switch along the epithelial-mesenchymal axis defines resistance and invasion downstream of ionising radiation in oral squamous cell carcinoma

BACKGROUND

Dynamic transitions of tumour cells along the epithelial-mesenchymal axis are important in tumorigenesis, metastasis and therapy resistance.

METHODS

In this study, we have used cell lines, 3D spheroids and tumour samples in a variety of cell biological and transcriptome analyses to highlight the cellular and molecular dynamics of OSCC response to ionising radiation.

RESULTS

Our study demonstrates a prominent hybrid epithelial-mesenchymal state in oral squamous cell carcinoma cells and tumour samples. We have further identified a key role for levels of E-cadherin in stratifying the hybrid cells to compartments with varying levels of radiation response and radiation-induced epithelial-mesenchymal transition. The response to radiation further entailed the generation of a new cell population with low expression levels of E-cadherin, and positive for Vimentin (ECAD^Low/Neg-VIM^Pos), a phenotypic signature that showed an enhanced capacity for radiation resistance and invasion. At the molecular level, transcriptome analysis of spheroids in response to radiation showed an initial burst of misregulation within the first 30 min that further declined, although still highlighting key alterations in gene signatures. Among others, pathway analysis showed an over-representation for the Wnt signalling pathway that was further confirmed to be functionally involved in the generation of ECAD^Low/Neg-VIM^Pos population, acting upstream of radiation resistance and tumour cell invasion.

CONCLUSION

This study highlights the functional significance and complexity of tumour cell remodelling in response to ionising radiation with links to resistance and invasive capacity. An area of less focus in conventional radiotherapy, with the potential to improve treatment outcomes and relapse-free survival.

Coupled Feedback Loops Involving PAGE4, EMT and Notch Signaling Can Give Rise to Non-genetic Heterogeneity in Prostate Cancer Cells

Non-genetic heterogeneity is emerging as a crucial factor underlying therapy resistance in multiple cancers. However, the design principles of regulatory networks underlying non-genetic heterogeneity in cancer remain poorly understood. Here, we investigate the coupled dynamics of feedback loops involving (a) oscillations in androgen receptor (AR) signaling mediated through an intrinsically disordered protein PAGE4, (b) multistability in epithelial-mesenchymal transition (EMT), and c) Notch-Delta-Jagged signaling mediated cell-cell communication, each of which can generate non-genetic heterogeneity through multistability and/or oscillations. Our results show how different coupling strengths between AR and EMT signaling can lead to monostability, bistability, or oscillations in the levels of AR, as well as propagation of oscillations to EMT dynamics. These results reveal the emergent dynamics of coupled oscillatory and multi-stable systems and unravel mechanisms by which non-genetic heterogeneity in AR levels can be generated, which can act as a barrier to most existing therapies for prostate cancer patients.

Multi-Stability and Consequent Phenotypic Plasticity in AMPK-Akt Double Negative Feedback Loop in Cancer Cells

Adaptation and survival of cancer cells to various stress and growth factor conditions is crucial for successful metastasis. A double-negative feedback loop between two serine/threonine kinases AMPK (AMP-activated protein kinase) and Akt can regulate the adaptation of breast cancer cells to matrix-deprivation stress. This feedback loop can significantly generate two phenotypes or cell states: matrix detachment-triggered pAMPK^high/ pAkt^low state, and matrix (re)attachment-triggered pAkt^high/ pAMPK^low state. However, whether these two cell states can exhibit phenotypic plasticity and heterogeneity in a given cell population, i.e., whether they can co-exist and undergo spontaneous switching to generate the other subpopulation, remains unclear. Here, we develop a mechanism-based mathematical model that captures the set of experimentally reported interactions among AMPK and Akt. Our simulations suggest that the AMPK-Akt feedback loop can give rise to two co-existing phenotypes (pAkt^high/ pAMPK^low and pAMPK^high/pAkt^low) in specific parameter regimes. Next, to test the model predictions, we segregated these two subpopulations in MDA-MB-231 cells and observed that each of them was capable of switching to another in adherent conditions. Finally, the predicted trends are supported by clinical data analysis of The Cancer Genome Atlas (TCGA) breast cancer and pan-cancer cohorts that revealed negatively correlated pAMPK and pAkt protein levels. Overall, our integrated computational-experimental approach unravels that AMPK-Akt feedback loop can generate multi-stability and drive phenotypic switching and heterogeneity in a cancer cell population.

EpCAM-Mediated Cellular Plasticity Promotes Radiation Resistance and Metastasis in Breast Cancer

Substantial number of breast cancer (BC) patients undergoing radiation therapy (RT) develop local recurrence over time. During RT therapy, cells can gradually acquire resistance implying adaptive radioresistance. Here we probe the mechanisms underlying this acquired resistance by first establishing radioresistant lines using ZR-75-1 and MCF-7 BC cells through repeated exposure to sub-lethal fractionated dose of 2Gy up to 15 fractions. Radioresistance was found to be associated with increased cancer stem cells (CSCs), and elevated EpCAM expression in the cell population. A retrospective analysis of TCGA dataset indicated positive correlation of high EpCAM expression with poor response to RT. Intriguingly, elevated EpCAM expression in the radioresistant CSCs raise the bigger question of how this biomarker expression contributes during radiation treatment in BC. Thereafter, we establish EpCAM overexpressing ZR-75-1 cells (ZR-75-1^EpCAM), which conferred radioresistance, increased stemness through enhanced AKT activation and induced a hybrid epithelial/mesenchymal phenotype with enhanced contractility and invasiveness. In line with these observations, orthotopic implantation of ZR-75-1^EpCAM cells exhibited faster growth, lesser sensitivity to radiation therapy and increased lung metastasis than baseline ZR-75-1 cells in mice. In summary, this study shows that similar to radioresistant BC cells, EpCAM overexpressing cells show high degree of plasticity and heterogeneity which ultimately induces radioresistant and metastatic behavior of cancer cells, thus aggravating the disease condition.

Cancer Stem Cell Plasticity - A Deadly Deal

Intratumoral heterogeneity is a major ongoing challenge in the effective therapeutic targeting of cancer. Accumulating evidence suggests that a fraction of cells within a tumor termed Cancer Stem Cells (CSCs) are primarily responsible for this diversity resulting in therapeutic resistance and metastasis. Adding to this complexity, recent studies have shown that there can be different subpopulations of CSCs with varying biochemical and biophysical traits resulting in varied dissemination and drug-resistance potential. Moreover, cancer cells can exhibit a high level of plasticity or the ability to dynamically switch between CSC and non-CSC states or among different subsets of CSCs. In addition, CSCs also display extensive metabolic plasticity. The molecular mechanisms underlying these different interconnected axes of plasticity has been under extensive investigation and the trans-differentiation process of Epithelial to Mesenchymal transition (EMT) has been identified as a major contributing factor. Besides genetic and epigenetic factors, CSC plasticity is also shaped by non-cell-autonomous effects such as the tumor microenvironment (TME). In this review, we discuss the latest developments in decoding mechanisms and implications of CSC plasticity in tumor progression at biochemical and biophysical levels, and the latest in silico approaches being taken for characterizing cancer cell plasticity. These efforts can help improve existing therapeutic approaches by taking into consideration the contribution of cellular plasticity/heterogeneity in enabling drug resistance.

New Insights Into the Role of Phenotypic Plasticity and EMT in Driving Cancer Progression

Tumor cells demonstrate substantial plasticity in their genotypic and phenotypic characteristics. Epithelial-mesenchymal plasticity (EMP) can be characterized into dynamic intermediate states and can be orchestrated by many factors, either intercellularly via epigenetic reprograming, or extracellularly via growth factors, inflammation and/or hypoxia generated by the tumor stromal microenvironment. EMP has the capability to alter phenotype and produce heterogeneity, and thus by changing the whole cancer landscape can attenuate oncogenic signaling networks, invoke anti-apoptotic features, defend against chemotherapeutics and reprogram angiogenic and immune recognition functions. We discuss here the role of phenotypic plasticity in tumor initiation, progression and metastasis and provide an update of the modalities utilized for the molecular characterization of the EMT states and attributes of cellular behavior, including cellular metabolism, in the context of EMP. We also summarize recent findings in dynamic EMP studies that provide new insights into the phenotypic plasticity of EMP flux in cancer and propose therapeutic strategies to impede the metastatic outgrowth of phenotypically heterogeneous tumors.

Enhanced antitumour efficacy of functionalized doxorubicin plus schisandrin B co-delivery liposomes via inhibiting epithelial-mesenchymal transition

Non-small cell lung cancer (NSCLC) is a malignant cancer characterized by easy invasion, metastasis and poor prognosis, so that conventional chemotherapy cannot inhibit its invasion and metastasis. Doxorubicin (DOX), as a broad-spectrum antitumour drug, cannot be widely used in clinic because of its poor targeting, short half-life, strong toxicity and side effects. Therefore, the aim of our study is to construct a kind of PFV modified DOX plus schisandrin B liposomes to solve the above problems, and to explore its potential mechanism of inhibiting NSCLC invasion and metastasis. The antitumour efficiency of the targeting liposomes was carried out by cytotoxicity, heating ablation, wound healing, transwell, vasculogenic mimicry channels formation and metastasis-related protein tests in vitro. Pharmacodynamics were evaluated by tumour inhibition rate, HE staining and TUNEL test in vivo. The enhanced anti-metastatic mechanism of the targeting liposomes was attributed to the downregulation of vimentin, vascular endothelial growth factor, matrix metalloproteinase 9 and upregulation of E-cadherin. In conclusion, the PFV modified DOX plus schisandrin B liposomes prepared in this study provided a treatment strategy with high efficiency for NSCLC.

Cell Fate Forecasting: A Data-Assimilation Approach to Predict Epithelial-Mesenchymal Transition

Epithelial-mesenchymal transition (EMT) is a fundamental biological process that plays a central role in embryonic development, tissue regeneration, and cancer metastasis. Transforming growth factor-β (TGFβ) is a potent inducer of this cellular transition, which is composed of transitions from an epithelial state to intermediate or partial EMT state(s) to a mesenchymal state. Using computational models to predict cell state transitions in a specific experiment is inherently difficult for reasons including model parameter uncertainty and error associated with experimental observations. In this study, we demonstrate that a data-assimilation approach using an ensemble Kalman filter, which combines limited noisy observations with predictions from a computational model of TGFβ-induced EMT, can reconstruct the cell state and predict the timing of state transitions. We used our approach in proof-of-concept “synthetic” in silico experiments, in which experimental observations were produced from a known computational model with the addition of noise. We mimic parameter uncertainty in in vitro experiments by incorporating model error that shifts the TGFβ doses associated with the state transitions and reproduces experimentally observed variability in cell state by either shifting a single parameter or generating “populations” of model parameters. We performed synthetic experiments for a wide range of TGFβ doses, investigating different cell steady-state conditions, and conducted parameter studies varying properties of the data-assimilation approach including the time interval between observations and incorporating multiplicative inflation, a technique to compensate for underestimation of the model uncertainty and mitigate the influence of model error. We find that cell state can be successfully reconstructed and the future cell state predicted in synthetic experiments, even in the setting of model error, when experimental observations are performed at a sufficiently short time interval and incorporate multiplicative inflation. Our study demonstrates the feasibility and utility of a data-assimilation approach to forecasting the fate of cells undergoing EMT.

The Physics of Cellular Decision Making During Epithelial-Mesenchymal Transition

The epithelial-mesenchymal transition (EMT) is a process by which cells lose epithelial traits, such as cell-cell adhesion and apico-basal polarity, and acquire migratory and invasive traits. EMT is crucial to embryonic development and wound healing. Misregulated EMT has been implicated in processes associated with cancer aggressiveness, including metastasis. Recent experimental advances such as single-cell analysis and temporal phenotypic characterization have established that EMT is a multistable process wherein cells exhibit and switch among multiple phenotypic states. This is in contrast to the classical perception of EMT as leading to a binary choice. Mathematical modeling has been at the forefront of this transformation for the field, not only providing a conceptual framework to integrate and analyze experimental data, but also making testable predictions. In this article, we review the key features and characteristics of EMT dynamics, with a focus on the mathematical modeling approaches that have been instrumental to obtaining various useful insights.

Controversies around epithelial-mesenchymal plasticity in cancer metastasis

Experimental evidence accumulated over decades has implicated epithelial-mesenchymal plasticity (EMP), which collectively encompasses epithelial-mesenchymal transition and the reverse process of mesenchymal-epithelial transition, in tumour metastasis, cancer stem cell generation and maintenance, and therapeutic resistance. However, the dynamic nature of EMP processes, the apparent need to reverse mesenchymal changes for the development of macrometastases and the likelihood that only minor cancer cell subpopulations exhibit EMP at any one time have made such evidence difficult to accrue in the clinical setting. In this Perspectives article, we outline the existing preclinical and clinical evidence for EMP and reflect on recent controversies, including the failure of initial lineage-tracing experiments to confirm a major role for EMP in dissemination, and discuss accumulating data suggesting that epithelial features and/or a hybrid epithelial-mesenchymal phenotype are important in metastasis. We also highlight strategies to address the complexities of therapeutically targeting the EMP process that give consideration to its spatially and temporally divergent roles in metastasis, with the view that this will yield a potent and broad class of therapeutic agents.

Dynamics of Phenotypic Heterogeneity Associated with EMT and Stemness during Cancer Progression

Genetic and phenotypic heterogeneity contribute to the generation of diverse tumor cell populations, thus enhancing cancer aggressiveness and therapy resistance. Compared to genetic heterogeneity, a consequence of mutational events, phenotypic heterogeneity arises from dynamic, reversible cell state transitions in response to varying intracellular/extracellular signals. Such phenotypic plasticity enables rapid adaptive responses to various stressful conditions and can have a strong impact on cancer progression. Herein, we have reviewed relevant literature on mechanisms associated with dynamic phenotypic changes and cellular plasticity, such as epithelial-mesenchymal transition (EMT) and cancer stemness, which have been reported to facilitate cancer metastasis. We also discuss how non-cell-autonomous mechanisms such as cell-cell communication can lead to an emergent population-level response in tumors. The molecular mechanisms underlying the complexity of tumor systems are crucial for comprehending cancer progression, and may provide new avenues for designing therapeutic strategies.

Multi-Omics Characterization of the Spontaneous Mesenchymal-Epithelial Transition in the PMC42 Breast Cancer Cell Lines

Epithelial-mesenchymal plasticity (EMP), encompassing epithelial-mesenchymal transition (EMT) and mesenchymal-epithelial transition (MET), are considered critical events for cancer metastasis. We investigated chromosomal heterogeneity and chromosomal instability (CIN) profiles of two sister PMC42 breast cancer (BC) cell lines to assess the relationship between their karyotypes and EMP phenotypic plasticity. Karyotyping by GTG banding and exome sequencing were aligned with SWATH quantitative proteomics and existing RNA-sequencing data from the two PMC42 cell lines; the mesenchymal, parental PMC42-ET cell line and the spontaneously epithelially shifted PMC42-LA daughter cell line. These morphologically distinct PMC42 cell lines were also compared with five other BC cell lines (MDA-MB-231, SUM-159, T47D, MCF-7 and MDA-MB-468) for their expression of EMP and cell surface markers, and stemness and metabolic profiles. The findings suggest that the epithelially shifted cell line has a significantly altered ploidy of chromosomes 3 and 13, which is reflected in their transcriptomic and proteomic expression profiles. Loss of the TGFβR2 gene from chromosome 3 in the epithelial daughter cell line inhibits its EMT induction by TGF-β stimulus. Thus, integrative 'omics' characterization established that the PMC42 system is a relevant MET model and provides insights into the regulation of phenotypic plasticity in breast cancer.