Morphological State Transition Dynamics in EGF-Induced Epithelial to Mesenchymal Transition

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.

Data- and theory-driven approaches for understanding paths of epithelial-mesenchymal transition

Reversible transitions between epithelial and mesenchymal cell states are a crucial form of epithelial plasticity for development and disease progression. Recent experimental data and mechanistic models showed multiple intermediate epithelial-mesenchymal transition (EMT) states as well as trajectories of EMT underpinned by complex gene regulatory networks. In this review, we summarize recent progress in quantifying EMT and characterizing EMT paths with computational methods and quantitative experiments including omics-level measurements. We provide perspectives on how these studies can help relating fundamental cell biology to physiological and pathological outcomes of EMT.

Single-cell tracking as a tool for studying EMT-phenotypes

This article demonstrates that label-free single-cell video tracking is a useful approach for in vitro studies of Epithelial-Mesenchymal Transition (EMT). EMT is a highly heterogeneous process, involved in wound healing, embryogenesis and cancer. The process promotes metastasis, and increased understanding can aid development of novel therapeutic strategies. The role of EMT-associated biomarkers depends on biological context, making it challenging to compare and interpret data from different studies. We demonstrate single-cell video tracking for comprehensive phenotype analysis. In this study we performed single-cell video tracking on 72-h long recordings. We quantified several behaviours at a single-cell level during induced EMT in MDA-MB-468 cells. This revealed notable variations in migration speed, with different dose-response patterns and varying distributions of speed. By registering cell morphologies during the recording, we determined preferred paths of morphological transitions. We also found a clear association between migration speed and cell morphology. We found elevated rates of cell death, diminished proliferation, and an increase in mitotic failures followed by re-fusion of sister-cells. The method allows tracking of phenotypes in cell lineages, which can be particularly useful in epigenetic studies. Sister-cells were found to have significant similarities in their speeds and morphologies, illustrating the heritability of these traits.

Simultaneous Proximity DNAzyme-Activated Duplexed Protein-Specific Glycosylation Imaging on Cell Surface via Bioorthogonal Chemistry

Due to the scarcity of strategies to evaluate the multiple subtype monosaccharides in one specific protein simultaneously within a single assay, understanding the glycosylation mechanisms and revealing their roles in disease development become extremely challenging. Herein, a strategy of proximity DNAzyme-activated fluorescence imaging of multiplex saccharides in a protein on the cell surface via bio-orthogonal chemistry is reported. The multichannel proximity DNAzyme-activated fluorescence recovery enabled the highly selective and effective imaging analysis of multiplexed protein-specific glycosylation in situ and has been demonstrated. This strategy is successfully applied to visualize the sialylation and fucosylation in four specific proteins on different cell lines and evaluate the variations of protein-specific glycosylation in response to the alterations of the cellular physiological status. More importantly, the quantitative tracking of the terminal sialyation and fucosylation changes at the single-protein level is realized by assigning the target protein as the native reference, which has the potential to be a versatile platform for glycobiology research and clinical diagnosis.

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.

Dynamics of Epithelial-Mesenchymal Plasticity: What Have Single-Cell Investigations Elucidated So Far?

Epithelial-mesenchymal plasticity (EMP) is a key driver of cancer metastasis and therapeutic resistance, through which cancer cells can reversibly and dynamically alter their molecular and functional traits along the epithelial-mesenchymal spectrum. While cells in the epithelial phenotype are usually tightly adherent, less metastatic, and drug-sensitive, those in the hybrid epithelial/mesenchymal and/or mesenchymal state are more invasive, migratory, drug-resistant, and immune-evasive. Single-cell studies have emerged as a powerful tool in gaining new insights into the dynamics of EMP across various cancer types. Here, we review many recent studies that employ single-cell analysis techniques to better understand the dynamics of EMP in cancer both in vitro and in vivo. These single-cell studies have underlined the plurality of trajectories cells can traverse during EMP and the consequent heterogeneity of hybrid epithelial/mesenchymal phenotypes seen at both preclinical and clinical levels. They also demonstrate how diverse EMP trajectories may exhibit hysteretic behavior and how the rate of such cell-state transitions depends on the genetic/epigenetic background of recipient cells, as well as the dose and/or duration of EMP-inducing growth factors. Finally, we discuss the relationship between EMP and patient survival across many cancer types. We also present a next set of questions related to EMP that could benefit much from single-cell observations and pave the way to better tackle phenotypic switching and heterogeneity in clinic.

Nanosilver-functionalized polysaccharides as a platform for wound dressing

Polysaccharides that are naturally sourced have enormous promise as wound dressings, due to their wider availability and reasonable cost and good biocompatibility. Furthermore, nanosilver extensively applied in wound treatment is attributed to its broad spectrum of antimicrobial effects and lesser drug resistance. Consequently, wound dressings in corporating nanosilver have attracted wide-scale interest in wound healing, and nanosilver-functionalized polysaccharide-based wound dressings present an affordable option for healing of chronic wounds. This review encompasses preparation methods, classification, and antibacterial performances of nanosilver wound dressings. The prospective research arenas of nanosilver-based wound polysaccharide dressings are also elaborated. The review attempts to include a summary of the most recent advancements in silver nanotechnology as well as guidance for the investigation of nanosilver-functionalized polysaccharide-based wound dressings.

Metabolomic analysis of circulating tumor cells derived liver metastasis of colorectal cancer

Metabolic reprogramming is one of the essential features of tumor that may dramatically contribute to metastasis and collapse. The metabolic profiling is investigated on the patient derived tissue and cancer cell line derived mouse metastasis xenograft. As well-recognized "seeds" for remote metastasis of tumor, role of circulating tumor cells (CTCs) in the study of metabolic reprogramming feature of tumor is yet to be elucidated. More specifically, whether there is difference of metabolic features of liver metastasis in colorectal cancer (CRC) derived from either CTCs or cancer cell line is still unknown. In this study, comprehensive untargeted metabolomics was performed using high performance liquid chromatography-mass spectrometry (HPLC-MS) in liver metastasis tissues from CT26 cells and CTCs derived mouse models. We identified 288 differential metabolites associated with the pathways such as one carbon pool by folate, folate biosynthesis and histidine metabolism through bioinformation analysis. Multiple gene expression was upregulated in the CTCs derived liver metastasis, specifically some specific enzymes. These results indicated that the metabolite phenotype and corresponding gene expression in the CTCs derived liver metastasis tissues was different from the parental CT26 cells, displaying a specific up-regulation of mRNAs involved in the above metabolism-related pathways. The metabolic profile of CTCs was characterized on the liver metastatic process in colorectal cancer. The invasion ability and chemo drug tolerance of the CTCs derived tumor and metastasis was found to be overwhelming higher than cell line derived counterpart. Identification of the differential metabolites will lead to a better understanding of the hallmarks of the cancer progression and metastasis, which may suggest potential attractive target for treating metastatic CRC.

Reconstructing data-driven governing equations for cell phenotypic transitions: integration of data science and systems biology

Cells with the same genome can exist in different phenotypes and can change between distinct phenotypes when subject to specific stimuli and microenvironments. Some examples include cell differentiation during development, reprogramming for induced pluripotent stem cells and transdifferentiation, cancer metastasis and fibrosis progression. The regulation and dynamics of cell phenotypic conversion is a fundamental problem in biology, and has a long history of being studied within the formalism of dynamical systems. A main challenge for mechanism-driven modeling studies is acquiring sufficient amount of quantitative information for constraining model parameters. Advances in quantitative experimental approaches, especially high throughput single-cell techniques, have accelerated the emergence of a new direction for reconstructing the governing dynamical equations of a cellular system from quantitative single-cell data, beyond the dominant statistical approaches. Here I review a selected number of recent studies using live- and fixed-cell data and provide my perspective on future development.

Reconstructing data-driven governing equations for cell phenotypic transitions: integration of data science and systems biology

Cells with the same genome can exist in different phenotypes and can change between distinct phenotypes when subject to specific stimuli and microenvironments. Some examples include cell differentiation during development, reprogramming for induced pluripotent stem cells and transdifferentiation, cancer metastasis and fibrosis progression. The regulation and dynamics of cell phenotypic conversion is a fundamental problem in biology, and has a long history of being studied within the formalism of dynamical systems. A main challenge for mechanism-driven modeling studies is acquiring sufficient amount of quantitative information for constraining model parameters. Advances in quantitative experimental approaches, especially high throughput single-cell techniques, have accelerated the emergence of a new direction for reconstructing the governing dynamical equations of a cellular system from quantitative single-cell data, beyond the dominant statistical approaches. Here I review a selected number of recent studies using live- and fixed-cell data and provide my perspective on future development.

Novel cost-efficient protein-membrane basedsystem for cells co-cultivation and modeling theintercellular communication

In vitro systems serve as compact and manipulate models to investigate interactions between different cell types. A homogeneous population of cells predictably and uniformly responds to external factors. In a heterogeneous cell population, the effect of external growth factors is perceived in the context of intercellular interactions. Indirect cell cocultivation allows one to observe the paracrine effects of cells and separately analyze cell populations. The article describes an application of custom-made cell co-cultivation systems based on protein membranes separated from the bottom of the vessel by the 3d printed holder or kept afloat by a magnetic field. Using the proposed co-cultivation system, we analyzed the interaction of A549 cells and fibroblasts, in the presence and absence of growth factors. During co-cultivation of cells, the expression of genes of the activation for epithelial and mesenchymal transitions decreases. The article proposes the application of a newly available system for the co-cultivation of different cell types. This article is protected by copyright. All rights reserved.

Interconnected high-dimensional landscapes of epithelial-mesenchymal plasticity and stemness in cancer

Establishing macrometastases at distant organs is a highly challenging process for cancer cells, with extremely high attrition rates. A very small percentage of disseminated cells have the ability to dynamically adapt to their changing micro-environments through reversibly switching to another phenotype, aiding metastasis. Such plasticity can be exhibited along one or more axes-epithelial-mesenchymal plasticity (EMP) and cancer stem cells (CSCs) being the two most studied, and often tacitly assumed to be synonymous. Here, we review the emerging concepts related to EMP and CSCs across multiple cancers. Both processes are multi-dimensional in nature; for instance, EMP can be defined on morphological, molecular and functional changes, which may or may not be synchronized. Similarly, self-renewal, multi-lineage potential, and resistance to anoikis and/or therapy may not all occur simultaneously in CSCs. Thus, understanding the complexity in defining EMP and CSCs is essential if we are to understand their contribution to cancer metastasis. This will require a more comprehensive understanding of the non-linearity of these processes. These processes are dynamic, reversible, and semi-independent in nature; cells traverse the inter-connected high-dimensional EMP and CSC landscapes in diverse paths, each of which may exhibit a distinct EMP-CSC coupling. Our proposed model offers a potential unifying framework for elucidating the coupled decision-making along these dimensions and highlights a key set of open questions to be answered.

Biophysical and biochemical attributes of hybrid epithelial/mesenchymal phenotypes

The Epithelial-Mesenchymal Transition (EMT) is a biological phenomenon associated with explicit phenotypic and molecular changes in cellular traits. Unlike the earlier-held popular belief of it being a binary process, EMT is now thought of as a landscape including diverse hybrid E/M phenotypes manifested by varying degrees of the transition. These hybrid cells can co-express both epithelial and mesenchymal markers and/or functional traits, and can possess the property of collective cell migration, enhanced tumor-initiating ability, and immune/targeted therapy-evasive features, all of which are often associated with worse patient outcomes. These characteristics of the hybrid E/M cells have led to a surge in studies that map their biophysical and biochemical hallmarks that can be helpful in exploiting their therapeutic vulnerabilities. This review discusses recent advances made in investigating hybrid E/M phenotype(s) from diverse biophysical and biochemical aspects by integrating live cell-imaging, cellular morphology quantification and mathematical modelling, and highlights a set of questions that remain unanswered about the dynamics of hybrid E/M states.

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.

Immunofluorescence Image Feature Analysis and Phenotype Scoring Pipeline for Distinguishing Epithelial-Mesenchymal Transition

Epithelial–mesenchymal transition (EMT) is an essential biological process, also implicated in pathological settings such as cancer metastasis, in which epithelial cells transdifferentiate into mesenchymal cells. We devised an image analysis pipeline to distinguish between tissues comprised of epithelial and mesenchymal cells, based on extracted features from immunofluorescence images of differing biochemical markers. Mammary epithelial cells were cultured with 0 (control), 2, 4, or 10 ng/mL TGF-β1, a well-established EMT-inducer. Cells were fixed, stained, and imaged for E-cadherin, actin, fibronectin, and nuclei via immunofluorescence microscopy. Feature selection was performed on different combinations of individual cell markers using a Bag-of-Features extraction. Control and high-dose images comprised the training data set, and the intermediate dose images comprised the testing data set. A feature distance analysis was performed to quantify differences between the treatment groups. The pipeline was successful in distinguishing between control (epithelial) and the high-dose (mesenchymal) groups, as well as demonstrating progress along the EMT process in the intermediate dose groups. Validation using quantitative PCR (qPCR) demonstrated that biomarker expression measurements were well-correlated with the feature distance analysis. Overall, we identified image pipeline characteristics for feature extraction and quantification of immunofluorescence images to distinguish progression of EMT.

Machine learning reveals mesenchymal breast carcinoma cell adaptation in response to matrix stiffness

Epithelial-mesenchymal transition (EMT) and its reverse process, mesenchymal-epithelial transition (MET), are believed to play key roles in facilitating the metastatic cascade. Metastatic lesions often exhibit a similar epithelial-like state to that of the primary tumour, in particular, by forming carcinoma cell clusters via E-cadherin-mediated junctional complexes. However, the factors enabling mesenchymal-like micrometastatic cells to resume growth and reacquire an epithelial phenotype in the target organ microenvironment remain elusive. In this study, we developed a workflow using image-based cell profiling and machine learning to examine morphological, contextual and molecular states of individual breast carcinoma cells (MDA-MB-231). MDA-MB-231 heterogeneous response to the host organ microenvironment was modelled by substrates with controllable stiffness varying from 0.2kPa (soft tissues) to 64kPa (bone tissues). We identified 3 distinct morphological cell types (morphs) varying from compact round-shaped to flattened irregular-shaped cells with lamellipodia, predominantly populating 2-kPa and >16kPa substrates, respectively. These observations were accompanied by significant changes in E-cadherin and vimentin expression. Furthermore, we demonstrate that the bone-mimicking substrate (64kPa) induced multicellular cluster formation accompanied by E-cadherin cell surface localisation. MDA-MB-231 cells responded to different substrate stiffness by morphological adaptation, changes in proliferation rate and cytoskeleton markers, and cluster formation on bone-mimicking substrate. Our results suggest that the stiffest microenvironment can induce MET.

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.

Measuring and Modelling the Epithelial- Mesenchymal Hybrid State in Cancer: Clinical Implications

The epithelial-mesenchymal (E/M) hybrid state has emerged as an important mediator of elements of cancer progression, facilitated by epithelial mesenchymal plasticity (EMP). We review here evidence for the presence, prognostic significance, and therapeutic potential of the E/M hybrid state in carcinoma. We further assess modelling predictions and validation studies to demonstrate stabilised E/M hybrid states along the spectrum of EMP, as well as computational approaches for characterising and quantifying EMP phenotypes, with particular attention to the emerging realm of single-cell approaches through RNA sequencing and protein-based techniques.

A Theoretical Approach to Coupling the Epithelial-Mesenchymal Transition (EMT) to Extracellular Matrix (ECM) Stiffness via LOXL2

Simple Summary

Epithelial-mesenchymal transition (EMT) is a key process in cancer progression through which cells weaken their cell-cell adhesion and gain mobility and invasive traits. Besides chemical signaling, recent studies have established the connection of EMT to mechanical microenvironment, such as the stiffness of extracellular matrix (ECM). LOXL2 is representative of a family of enzymes that promotes fiber cross-linking in ECM. With increased cross-linking comes increased stiffness, which induces EMT that can, in turn, elevate LOXL2 levels. As such, a positive feedback loop among EMT, LOXL2, and ECM stiffness can be formed. We built a mathematical model on a core biochemical reaction network featuring this feedback loop, and showed how strongly it drives EMT. We also illustrated mechanistically how cross-linking connects with stiffness, using a mechanical model of collagen (a major component of ECM). Using this theoretical framework, we demonstrated the heterogeneity of LOXL2/stiffness and its implications on migrating cancer cells that could seed metastasis, the growth of secondary malignant tumors. This framework can inspire experimental studies of more fine-grained mechanotransduction and biomechanical heterogeneity in cancers.

Abstract

The epithelial-mesenchymal transition (EMT) plays a critical role in cancer progression, being responsible in many cases for the onset of the metastatic cascade and being integral in the ability of cells to resist drug treatment. Most studies of EMT focus on its induction via chemical signals such as TGF-β or Notch ligands, but it has become increasingly clear that biomechanical features of the microenvironment such as extracellular matrix (ECM) stiffness can be equally important. Here, we introduce a coupled feedback loop connecting stiffness to the EMT transcription factor ZEB1, which acts via increasing the secretion of LOXL2 that leads to increased cross-linking of collagen fibers in the ECM. This increased cross-linking can effectively increase ECM stiffness and increase ZEB1 levels, thus setting a positive feedback loop between ZEB1 and ECM stiffness. To investigate the impact of this non-cell-autonomous effect, we introduce a computational approach capable of connecting LOXL2 concentration to increased stiffness and thereby to higher ZEB1 levels. Our results indicate that this positive feedback loop, once activated, can effectively lock the cells in a mesenchymal state. The spatial-temporal heterogeneity of the LOXL2 concentration and thus the mechanical stiffness also has direct implications for migrating cells that attempt to escape the primary tumor.

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.

Hybrid E/M Phenotype(s) and Stemness: A Mechanistic Connection Embedded in Network Topology

Metastasis remains an unsolved clinical challenge. Two crucial features of metastasizing cancer cells are (a) their ability to dynamically move along the epithelial-hybrid-mesenchymal spectrum and (b) their tumor initiation potential or stemness. With increasing functional characterization of hybrid epithelial/mesenchymal (E/M) phenotypes along the spectrum, recent in vitro and in vivo studies have suggested an increasing association of hybrid E/M phenotypes with stemness. However, the mechanistic underpinnings enabling this association remain unclear. Here, we develop a mechanism-based mathematical modeling framework that interrogates the emergent nonlinear dynamics of the coupled network modules regulating E/M plasticity (miR-200/ZEB) and stemness (LIN28/let-7). Simulating the dynamics of this coupled network across a large ensemble of parameter sets, we observe that hybrid E/M phenotype(s) are more likely to acquire stemness relative to "pure" epithelial or mesenchymal states. We also integrate multiple "phenotypic stability factors" (PSFs) that have been shown to stabilize hybrid E/M phenotypes both in silico and in vitro-such as OVOL1/2, GRHL2, and NRF2-with this network, and demonstrate that the enrichment of hybrid E/M phenotype(s) with stemness is largely conserved in the presence of these PSFs. Thus, our results offer mechanistic insights into recent experimental observations of hybrid E/M phenotype(s) that are essential for tumor initiation and highlight how this feature is embedded in the underlying topology of interconnected EMT (Epithelial-Mesenchymal Transition) and stemness networks.

Hypoxia, partial EMT and collective migration: Emerging culprits in metastasis

Epithelial-mesenchymal transition (EMT) is a cellular biological process involved in migration of primary cancer cells to secondary sites facilitating metastasis. Besides, EMT also confers properties such as stemness, drug resistance and immune evasion which can aid a successful colonization at the distant site. EMT is not a binary process; recent evidence suggests that cells in partial EMT or hybrid E/M phenotype(s) can have enhanced stemness and drug resistance as compared to those undergoing a complete EMT. Moreover, partial EMT enables collective migration of cells as clusters of circulating tumor cells or emboli, further endorsing that cells in hybrid E/M phenotypes may be the 'fittest' for metastasis. Here, we review mechanisms and implications of hybrid E/M phenotypes, including their reported association with hypoxia. Hypoxia-driven activation of HIF-1α can drive EMT. In addition, cyclic hypoxia, as compared to acute or chronic hypoxia, shows the highest levels of active HIF-1α and can augment cancer aggressiveness to a greater extent, including enriching for a partial EMT phenotype. We also discuss how metastasis is influenced by hypoxia, partial EMT and collective cell migration, and call for a better understanding of interconnections among these mechanisms. We discuss the known regulators of hypoxia, hybrid EMT and collective cell migration and highlight the gaps which needs to be filled for connecting these three axes which will increase our understanding of dynamics of metastasis and help control it more effectively.

Epigenetic feedback and stochastic partitioning during cell division can drive resistance to EMT

Epithelial-mesenchymal transition (EMT) and its reverse process mesenchymal-epithelial transition (MET) are central to metastatic aggressiveness and therapy resistance in solid tumors. While molecular determinants of both processes have been extensively characterized, the heterogeneity in the response of tumor cells to EMT and MET inducers has come into focus recently, and has been implicated in the failure of anti-cancer therapies. Recent experimental studies have shown that some cells can undergo an irreversible EMT depending on the duration of exposure to EMT-inducing signals. While the irreversibility of MET, or equivalently, resistance to EMT, has not been studied in as much detail, evidence supporting such behavior is slowly emerging. Here, we identify two possible mechanisms that can underlie resistance of cells to undergo EMT: epigenetic feedback in ZEB1/GRHL2 feedback loop and stochastic partitioning of biomolecules during cell division. Identifying the ZEB1/GRHL2 axis as a key determinant of epithelial-mesenchymal plasticity across many cancer types, we use mechanistic mathematical models to show how GRHL2 can be involved in both the abovementioned processes, thus driving an irreversible MET. Our study highlights how an isogenic population may contain subpopulation with varying degrees of susceptibility or resistance to EMT, and proposes a next set of questions for detailed experimental studies characterizing the irreversibility of MET/resistance to EMT.

DEBay: A computational tool for deconvolution of quantitative PCR data for estimation of cell type-specific gene expression in a mixed population

The expression of a gene is commonly estimated by quantitative PCR (qPCR) using RNA isolated from a large number of pooled cells. Such pooled samples often have subpopulations of cells with different levels of expression of the target gene. Estimation of gene expression from an ensemble of cells obscures the pattern of expression in different subpopulations. Physical separation of various subpopulations is a demanding task. We have developed a computational tool, Deconvolution of Ensemble through Bayes-approach (DEBay), to estimate cell type-specific gene expression from qPCR data of a mixed population. DEBay estimates Normalized Gene Expression Coefficient (NGEC), which is a relative measure of the expression of the target gene in each cell type in a population. NGEC has a direct algebraic correspondence with the normalized fold change in gene expression measured by qPCR. DEBay can deconvolute both time-dependent and -independent gene expression profiles. It uses the Bayesian method of model selection and parameter estimation. We have evaluated DEBay using synthetic and real experimental data. DEBay is implemented in Python. A GUI of DEBay and its source code are available for download at SourceForge (https://sourceforge.net/projects/debay).

Spheroid-Based Approach to Assess the Tissue Relevance of Analysis of Dispersed-Settled Tissue Cells by Cytometry of the Reaction Rate Constant

Cytometry of Reaction Rate Constant (CRRC) uses time-lapse fluorescence microscopy to measure a rate constant of a catalytic reaction in individual cells and, thus, facilitate accurate size determination for cell subpopulations with distinct efficiencies of this reaction. Reliable CRRC requires uniform exposure of cells to the reaction substrate followed by their uniform imaging, which in turn, requires that a tissue sample be disintegrated into a suspension of dispersed cells, and these cells settle on the support surface before being analyzed by CRRC. We call such cells "dispersed-settled" to distinguish them from cells cultured as a monolayer. Studies of the dispersed-settled cells can be tissue-relevant only if the cells maintain their 3D tissue state during the multi-hour CRRC procedure. Here, we propose an approach for assessing tissue relevance of the CRRC-based analysis of the dispersed-settled cells. Our approach utilizes cultured multicellular spheroids as a 3D cell model and cultured cell monolayers as a 2D cell model. The CRRC results of the dispersed-settled cells derived from spheroids are compared to those of spheroids and monolayers in order to find if the dispersed-settled cells are representative of the spheroids. To demonstrate its practical use, we applied this approach to a cellular reaction of multidrug resistance (MDR) transport, which was followed by extrusion of a fluorescent substrate from the cells. The approach proved to be reliable and revealed long-term maintenance of MDR transport in the dispersed-settled cells obtained from cultured ovarian cancer spheroids. Accordingly, CRRC can be used to determine accurately the size of a cell subpopulation with an elevated level of MDR transport in tumor samples, which makes CRRC a suitable method for the development of MDR-based predictors of chemoresistance. The proposed spheroid-based approach for validation of CRRC is applicable to other types of cellular reactions and, thus, will be an indispensable tool for transforming CRRC from an experimental technique into a practical analytical tool.

Effect of Kelulut Honey on the Cellular Dynamics of TGFβ-Induced Epithelial to Mesenchymal Transition in Primary Human Keratinocytes

Over-induction of epithelial to mesenchymal transition (EMT) by tumor growth factor beta (TGFβ) in keratinocytes is a key feature in keloid scar. The present work seeks to investigate the effect of Kelulut honey (KH) on TGFβ-induced EMT in human primary keratinocytes. Image analysis of the real time observation of TGFβ-induced keratinocytes revealed a faster wound closure and individual migration velocity compared to the untreated control. TGFβ-induced keratinocytes also have reduced circularity and display a classic EMT protein expression. Treatment of 0.0015% (v/v) KH reverses these effects. In untreated keratinocytes, KH resulted in slower initial wound closure and individual migration velocity, which sped up later on, resulting in greater wound closure at the final time point. KH treatment also led to greater directional migration compared to the control. KH treatment caused reduced circularity in keratinocytes but displayed a partial EMT protein expression. Taken together, the findings suggest the therapeutic potential of KH in preventing keloid scar by attenuating TGFβ-induced EMT.

Comparative Study of Transcriptomics-Based Scoring Metrics for the Epithelial-Hybrid-Mesenchymal Spectrum

The Epithelial-mesenchymal transition (EMT) is a cellular process implicated in embryonic development, wound healing, and pathological conditions such as cancer metastasis and fibrosis. Cancer cells undergoing EMT exhibit enhanced aggressive behavior characterized by drug resistance, tumor-initiation potential, and the ability to evade the immune system. Recent in silico, in vitro, and in vivo evidence indicates that EMT is not an all-or-none process; instead, cells can stably acquire one or more hybrid epithelial/mesenchymal (E/M) phenotypes which often can be more aggressive than purely E or M cell populations. Thus, the EMT status of cancer cells can prove to be a critical estimate of patient prognosis. Recent attempts have employed different transcriptomics signatures to quantify EMT status in cell lines and patient tumors. However, a comprehensive comparison of these methods, including their accuracy in identifying cells in the hybrid E/M phenotype(s), is lacking. Here, we compare three distinct metrics that score EMT on a continuum, based on the transcriptomics signature of individual samples. Our results demonstrate that these methods exhibit good concordance among themselves in quantifying the extent of EMT in a given sample. Moreover, scoring EMT using any of the three methods discerned that cells can undergo varying extents of EMT across tumor types. Separately, our analysis also identified tumor types with maximum variability in terms of EMT and associated an enrichment of hybrid E/M signatures in these samples. Moreover, we also found that the multinomial logistic regression (MLR)-based metric was capable of distinguishing between "pure" individual hybrid E/M vs. mixtures of E and M cells. Our results, thus, suggest that while any of the three methods can indicate a generic trend in the EMT status of a given cell, the MLR method has two additional advantages: (a) it uses a small number of predictors to calculate the EMT score and (b) it can predict from the transcriptomic signature of a population whether it is comprised of "pure" hybrid E/M cells at the single-cell level or is instead an ensemble of E and M cell subpopulations.

Comparative Study of Transcriptomics-Based Scoring Metrics for the Epithelial-Hybrid-Mesenchymal Spectrum

The Epithelial-mesenchymal transition (EMT) is a cellular process implicated in embryonic development, wound healing, and pathological conditions such as cancer metastasis and fibrosis. Cancer cells undergoing EMT exhibit enhanced aggressive behavior characterized by drug resistance, tumor-initiation potential, and the ability to evade the immune system. Recent in silico, in vitro, and in vivo evidence indicates that EMT is not an all-or-none process; instead, cells can stably acquire one or more hybrid epithelial/mesenchymal (E/M) phenotypes which often can be more aggressive than purely E or M cell populations. Thus, the EMT status of cancer cells can prove to be a critical estimate of patient prognosis. Recent attempts have employed different transcriptomics signatures to quantify EMT status in cell lines and patient tumors. However, a comprehensive comparison of these methods, including their accuracy in identifying cells in the hybrid E/M phenotype(s), is lacking. Here, we compare three distinct metrics that score EMT on a continuum, based on the transcriptomics signature of individual samples. Our results demonstrate that these methods exhibit good concordance among themselves in quantifying the extent of EMT in a given sample. Moreover, scoring EMT using any of the three methods discerned that cells can undergo varying extents of EMT across tumor types. Separately, our analysis also identified tumor types with maximum variability in terms of EMT and associated an enrichment of hybrid E/M signatures in these samples. Moreover, we also found that the multinomial logistic regression (MLR)-based metric was capable of distinguishing between "pure" individual hybrid E/M vs. mixtures of E and M cells. Our results, thus, suggest that while any of the three methods can indicate a generic trend in the EMT status of a given cell, the MLR method has two additional advantages: (a) it uses a small number of predictors to calculate the EMT score and (b) it can predict from the transcriptomic signature of a population whether it is comprised of "pure" hybrid E/M cells at the single-cell level or is instead an ensemble of E and M cell subpopulations.

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.