The exostosin glycosyltransferase 1/STAT3 axis is a driver of breast cancer aggressiveness.

The epithelial-mesenchymal transition (EMT) program is crucial for transforming carcinoma cells into a partially mesenchymal state, enhancing their chemoresistance, migration, and metastasis. This shift in cell state is tightly regulated by cellular mechanisms that are not yet fully characterized. One intriguing EMT aspect is the rewiring of the proteoglycan landscape, particularly the induction of heparan sulfate proteoglycan (HSPG) biosynthesis. This proteoglycan functions as a co-receptor that accelerates cancer-associated signaling pathways through its negatively-charged residues. However, the precise mechanisms through which EMT governs HSPG biosynthesis and its role in cancer cell plasticity remain elusive. Here, we identified exostosin glycosyltransferase 1 (EXT1), a central enzyme in HSPG biosynthesis, to be selectively upregulated in aggressive tumor subtypes and cancer cell lines, and to function as a key player in breast cancer aggressiveness. Notably, ectopic expression of EXT1 in epithelial cells is sufficient to induce HSPG levels and the expression of known mesenchymal markers, subsequently enhancing EMT features, including cell migration, invasion, and tumor formation. Additionally, EXT1 loss in MDA-MB-231 cells inhibits their aggressiveness-associated traits such as migration, chemoresistance, tumor formation, and metastasis. Our findings reveal that EXT1, through its role in HSPG biosynthesis, governs signal transducer and activator of transcription 3 (STAT3) signaling, a known regulator of cancer cell aggressiveness. Collectively, we present the EXT1/HSPG/STAT3 axis as a central regulator of cancer cell plasticity that directly links proteoglycan synthesis to oncogenic signaling pathways.

Perspectives on polarity - exploring biological asymmetry across scales.

In this Perspective, Journal of Cell Science invited researchers working on cell and tissue polarity to share their thoughts on unique, emerging or open questions relating to their field. The goal of this article is to feature 'voices' from scientists around the world and at various career stages, to bring attention to innovative and thought-provoking topics of interest to the cell biology community. These voices discuss intriguing questions that consider polarity across scales, evolution, development and disease. What can yeast and protists tell us about the evolution of cell and tissue polarity in animals? How are cell fate and development influenced by emerging dynamics in cell polarity? What can we learn from atypical and extreme polarity systems? How can we arrive at a more unified biophysical understanding of polarity? Taken together, these pieces demonstrate the broad relevance of the fascinating phenomenon of cell polarization to diverse fundamental biological questions.

p53 amyloid pathology is correlated with higher cancer grade irrespective of the mutant or wild-type form.

p53 (also known as TP53) mutation and amyloid formation are long associated with cancer pathogenesis; however, the direct demonstration of the link between p53 amyloid load and cancer progression is lacking. Using multi-disciplinary techniques and 59 tissues (53 oral and stomach cancer tumor tissue samples from Indian individuals with cancer and six non-cancer oral and stomach tissue samples), we showed that p53 amyloid load and cancer grades are highly correlated. Furthermore, next-generation sequencing (NGS) data suggest that not only mutant p53 (e.g. single-nucleotide variants, deletions, and insertions) but wild-type p53 also formed amyloids either in the nucleus (50%) and/or in the cytoplasm in most cancer tissues. Interestingly, in all these cancer tissues, p53 displays a loss of DNA-binding and transcriptional activities, suggesting that the level of amyloid load correlates with the degree of loss and an increase in cancer grades. The p53 amyloids also sequester higher amounts of the related p63 and p73 (also known as TP63 and TP73, respectively) protein in higher-grade tumor tissues. The data suggest p53 misfolding and/or aggregation, and subsequent amyloid formation, lead to loss of the tumor-suppressive function and the gain of oncogenic function, aggravation of which might determine the cancer grade.

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.

Spatial heterogeneity in tumor adhesion qualifies collective cell invasion.

Collective cell invasion (CCI), a canon of most invasive solid tumors, is an emergent property of the interactions between cancer cells and their surrounding extracellular matrix (ECM). However, tumor populations invariably consist of cells expressing variable levels of adhesive proteins that mediate such interactions, disallowing an intuitive understanding of how tumor invasiveness at a multicellular scale is influenced by spatial heterogeneity of cell-cell and cell-ECM adhesion. Here, we have used a Cellular Potts model-based multiscale computational framework that is constructed on the histopathological principles of glandular cancers. In earlier efforts on homogenous cancer cell populations, this framework revealed the relative ranges of interactions, including cell-cell and cell-ECM adhesion that drove collective, dispersed, and mixed multimodal invasion. Here, we constitute a tumor core of two separate cell subsets showing distinct intra- and inter-subset cell-cell or cell-ECM adhesion strengths. These two subsets of cells are arranged to varying extents of spatial intermingling, which we call the heterogeneity index (HI). We observe that low and high inter-subset cell adhesion favors invasion of high-HI and low-HI intermingled populations with distinct intra-subset cell-cell adhesion strengths, respectively. In addition, for explored values of cell-ECM adhesion strengths, populations with high HI values collectively invade better than those with lower HI values. We then asked how spatial invasion is regulated by progressively intermingled cellular subsets that are epithelial, i.e., showed high cell-cell but poor cell-ECM adhesion, and mesenchymal, i.e., with reversed adhesion strengths to the former. Here too, inter-subset adhesion plays an important role in contextualizing the proportionate relationship between HI and invasion. An exception to this relationship is seen for cases of heterogeneous cell-ECM adhesion where sub-maximal HI patterns with higher outer localization of cells with stronger ECM adhesion collectively invade better than their relatively higher-HI counterparts. Our simulations also reveal how adhesion heterogeneity qualifies collective invasion, when either cell-cell or cell-ECM adhesion type is varied but results in an invasive dispersion when both adhesion types are simultaneously altered.

Oncogenic gain of function due to p53 amyloids occurs through aberrant alteration of cell cycle and proliferation.

Transcription factor p53 (also known as TP53) has been shown to aggregate into cytoplasmic and nuclear inclusions, compromising its native tumor suppressive functions. Recently, p53 has been shown to form amyloids, which play a role in conferring cancerous properties to cells, leading to tumorigenesis. However, the exact pathways involved in p53 amyloid-mediated cellular transformations are unknown. Here, using an in cellulo model of full-length p53 amyloid formation, we demonstrate the mechanism of loss of p53 tumor-suppressive function with concomitant oncogenic gain of functions. Global gene expression profiling of cells suggests that p53 amyloid formation dysregulates genes associated with the cell cycle, proliferation, apoptosis and senescence along with major signaling pathways. This is further supported by a proteome analysis, showing a significant alteration in levels of p53 target proteins and enhanced metabolism, which enables the survival of cells. Our data indicate that specifically targeting the key molecules in pathways affected by p53 amyloid formation, such as cyclin-dependent kinase-1, leads to loss of the oncogenic phenotype and induces apoptosis of cells. Overall, our work establishes the mechanism of the transformation of cells due to p53 amyloids leading to cancer pathogenesis. This article has an associated First Person interview with the first author of the paper.

Intrinsically Disordered Proteins: Critical Components of the Wetware.

Despite the wealth of knowledge gained about intrinsically disordered proteins (IDPs) since their discovery, there are several aspects that remain unexplored and, hence, poorly understood. A living cell is a complex adaptive system that can be described as a wetware─a metaphor used to describe the cell as a computer comprising both hardware and software and attuned to logic gates─capable of "making" decisions. In this focused Review, we discuss how IDPs, as critical components of the wetware, influence cell-fate decisions by wiring protein interaction networks to keep them minimally frustrated. Because IDPs lie between order and chaos, we explore the possibility that they can be modeled as attractors. Further, we discuss how the conformational dynamics of IDPs manifests itself as conformational noise, which can potentially amplify transcriptional noise to stochastically switch cellular phenotypes. Finally, we explore the potential role of IDPs in prebiotic evolution, in forming proteinaceous membrane-less organelles, in the origin of multicellularity, and in protein conformation-based transgenerational inheritance of acquired characteristics. Together, these ideas provide a new conceptual framework to discern how IDPs may perform critical biological functions despite their lack of structure.

Identification of EMT signaling cross-talk and gene regulatory networks by single-cell RNA sequencing.

The epithelial-to-mesenchymal transition (EMT) plays a critical role during normal development and in cancer progression. EMT is induced by various signaling pathways, including TGF-ß, BMP, Wnt-ß-catenin, NOTCH, Shh, and receptor tyrosine kinases. In this study, we performed single-cell RNA sequencing on MCF10A cells undergoing EMT by TGF-ß1 stimulation. Our comprehensive analysis revealed that cells progress through EMT at different paces. Using pseudotime clustering reconstruction of gene-expression profiles during EMT, we found sequential and parallel activation of EMT signaling pathways. We also observed various transitional cellular states during EMT. We identified regulatory signaling nodes that drive EMT with the expression of important microRNAs and transcription factors. Using a random circuit perturbation methodology, we demonstrate that the NOTCH signaling pathway acts as a key driver of TGF-ß-induced EMT. Furthermore, we demonstrate that the gene signatures of pseudotime clusters corresponding to the intermediate hybrid EMT state are associated with poor patient outcome. Overall, this study provides insight into context-specific drivers of cancer progression and highlights the complexities of the EMT process.

Emerging perspectives on growth factor metabolic relationships in the ovarian cancer ascites environment.

The ascites ecosystem in ovarian cancer is inhabited by complex cell types and is bathed in an environment rich in cytokines, chemokines, and growth factors that directly and indirectly impact metabolism of cancer cells and tumor associated cells. This milieu of malignant ascites, provides a 'rich' environment for the disease to thrive, contributing to every aspect of advanced ovarian cancer, a devastating gynecological cancer with a significant gap in targeted therapeutics. In this perspective we focus our discussions on the 'acellular' constituents of this liquid malignant tumor microenvironment, and how they influence metabolic pathways. Growth factors, chemokines and cytokines are known modulators of metabolism and have been shown to impact nutrient uptake and metabolic flexibility of tumors, yet few studies have explored how their enrichment in malignant ascites of ovarian cancer patients contributes to the metabolic requirements of ascites-resident cells. We focus here on TGF-ßs, VEGF and ILs, which are frequently elevated in ovarian cancer ascites and have all been described to have direct or indirect effects on metabolism, often through gene regulation of metabolic enzymes. We summarize what is known, describe gaps in knowledge, and provide examples from other tumor types to infer potential unexplored roles and mechanisms for ovarian cancer. The distribution and variation in acellular ascites components between patients poses both a challenge and opportunity to further understand how the ascites may contribute to disease heterogeneity. The review also highlights opportunities for studies on ascites-derived factors in regulating the ascites metabolic environment that could act as a unique signature in aiding clinical decisions in the future.

Effects of microRNA-mediated negative feedback on gene expression noise.

MicroRNAs (miRNAs) are small noncoding RNAs that regulate gene expression post-transcriptionally in eukaryotes by binding with target mRNAs and preventing translation. miRNA-mediated feedback motifs are ubiquitous in various genetic networks that control cellular decision making. A key question is how such a feedback mechanism may affect gene expression noise. To answer this, we have developed a mathematical model to study the effects of a miRNA-dependent negative-feedback loop on mean expression and noise in target mRNAs. Combining analytics and simulations, we show the existence of an expression threshold demarcating repressed and expressed regimes in agreement with earlier studies. The steady-state mRNA distributions are bimodal near the threshold, where copy numbers of mRNAs and miRNAs exhibit enhanced anticorrelated fluctuations. Moreover, variation of negative-feedback strength shifts the threshold locations and modulates the noise profiles. Notably, the miRNA-mRNA binding affinity and feedback strength collectively shape the bimodality. We also compare our model with a direct auto-repression motif, where a gene produces its own repressor. Auto-repression fails to produce bimodal mRNA distributions as found in miRNA-based indirect repression, suggesting the crucial role of miRNAs in creating phenotypic diversity. Together, we demonstrate how miRNA-dependent negative feedback modifies the expression threshold and leads to a broader parameter regime of bimodality compared to the no-feedback case.

Unique sphingolipid signature identifies luminal and triple-negative breast cancer subtypes.

Breast cancer (luminal and triple-negative breast cancer [TNBC]) is the most common cancer among women in India and worldwide. Altered sphingolipid levels have emerged as a common phenomenon during cancer progression. However, these alterations are yet to be translated into robust diagnostic and prognostic markers for cancer. Here, we present the quantified sphingolipids of tumor and adjacent-normal tissues from patients of luminal (n = 70) and TNBC (n = 42) subtype from an Indian cohort using targeted liquid chromatography mass spectrometry. We recorded unique sphingolipid profiles that distinguished luminal and TNBC tumors in comparison to adjacent normal tissue by six-sphingolipid signatures. Moreover, systematic comparison of the profiles of luminal and TNBC tumors provided a unique five-sphingolipid signature distinguishing the two subtypes. We further identified key sphingolipids that can stratify grade II and grade III tumors of luminal and TNBC subtype as well as their lymphovascular invasion status. Therefore, we provide the right evidence to develop these candidate sphingolipids as widely acceptable marker/s capable of diagnosing luminal vs TNBC subtype of breast cancer, and predicting the disease severity by identifying the tumor grade.

Theoretical and computational tools to model multistable gene regulatory networks.

The last decade has witnessed a surge of theoretical and computational models to describe the dynamics of complex gene regulatory networks, and how these interactions can give rise to multistable and heterogeneous cell populations. As the use of theoretical modeling to describe genetic and biochemical circuits becomes more widespread, theoreticians with mathematical and physical backgrounds routinely apply concepts from statistical physics, non-linear dynamics, and network theory to biological systems. This review aims at providing a clear overview of the most important methodologies applied in the field while highlighting current and future challenges. It also includes hands-on tutorials to solve and simulate some of the archetypical biological system models used in the field. Furthermore, we provide concrete examples from the existing literature for theoreticians that wish to explore this fast-developing field. Whenever possible, we highlight the similarities and differences between biochemical and regulatory networks and 'classical' systems typically studied in non-equilibrium statistical and quantum mechanics.

Coupled mutual inhibition and mutual activation motifs as tools for cell-fate control.

Multistability is central to biological systems as it plays a crucial role in adaptation, evolvability, and differentiation. The presence of positive feedback loops can enable multistability. The simplest of such feedback loops are a) a mutual inhibition loop (MI), b) a mutual activation loop (MA), and c) self-activation, all three of them known to give rise to bistability. However, the characteristic differences in the bistability exhibited by these motifs are relatively less understood. Here, we use dynamical simulations across a large ensemble of parameter sets and initial conditions to study the bistability characteristics of these motifs. Furthermore, we investigate the utility of these motifs for achieving coordinated expression through cyclic and parallel coupling amongst them. Our analysis revealed that MI-based architectures offer discrete and robust control over gene expression, multistability, and coordinated expression among multiple genes, as compared to MA-based architectures. We then devised a combination of MI and MA architectures to improve coordination and multistability. Such designs help improve our understanding of the control structures involved in robust cell-fate decisions and provide a way to achieve controlled decision-making in synthetic systems.

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.

Calcium signaling induces a partial EMT.

Epithelial plasticity, or epithelial-to-mesenchymal transition (EMT), is a well-recognized form of cellular plasticity, which endows tumor cells with invasive properties and alters their sensitivity to various agents, thus representing a major challenge to cancer therapy. It is increasingly accepted that carcinoma cells exist along a continuum of hybrid epithelial-mesenchymal (E-M) states and that cells exhibiting such partial EMT (P-EMT) states have greater metastatic competence than those characterized by either extreme (E or M). We described recently a P-EMT program operating in vivo by which carcinoma cells lose their epithelial state through post-translational programs. Here, we investigate the underlying mechanisms and report that prolonged calcium signaling induces a P-EMT characterized by the internalization of membrane-associated E-cadherin (ECAD) and other epithelial proteins as well as an increase in cellular migration and invasion. Signaling through Gαq-associated G-protein-coupled receptors (GPCRs) recapitulates these effects, which operate through the downstream activation of calmodulin-Camk2b signaling. These results implicate calcium signaling as a trigger for the acquisition of hybrid/partial epithelial-mesenchymal states in carcinoma cells.

Lhx2 in germ cells suppresses endothelial cell migration in the developing ovary.

LIM-homeobox genes play multiple roles in developmental processes, but their roles in gonad development are not completely understood. Herein, we report that Lhx2, Ils2, Lmx1a, and Lmx1b are expressed in a sexually dimorphic manner in mouse, rat, and human gonads during sex determination. Amongst these, Lhx2 has female biased expression in the developing gonads of species with environmental and genetic modes of sex determination. Single-cell RNAseq analysis revealed that Lhx2 is exclusively expressed in the germ cells of the developing mouse ovaries. To elucidate the roles of Lhx2 in the germ cells, we analyzed the phenotypes of Lhx2 knockout XX gonads. While the gonads developed appropriately in Lhx2 knockout mice and the somatic cells were correctly specified in the developing ovaries, transcriptome analysis revealed enrichment of genes in the angiogenesis pathway. There was an elevated expression of several pro-angiogenic factors in the Lhx2 knockout ovaries. The elevated expression of pro-angiogenic factors was associated with an increase in numbers of endothelial cells in the Lhx2-/- ovaries at E13.5. Gonad recombination assays revealed that the increased numbers of endothelial cells in the XX gonads in absence of Lhx2 was due to ectopic migration of endothelial cells in a cell non-autonomous manner. We also found that, there was increased expression of several endothelial cell-enriched male-biased genes in Lhx2 knockout ovaries. Also, in absence of Lhx2, the migrated endothelial cells formed an angiogenic network similar to that of the wild type testis, although the coelomic blood vessel did not form. Together, our results suggest that Lhx2 in the germ cells is required to suppress vascularization in the developing ovary. These results suggest a need to explore the roles of germ cells in the control of vascularization in developing gonads. Preprint version of the article is available on BioRxiv at https://doi.org/10.1101/2022.03.07.483280.

Robustness in phenotypic plasticity and heterogeneity patterns enabled by EMT networks.

Epithelial-mesenchymal plasticity (EMP) is a key arm of cancer metastasis and is observed across many contexts. Cells undergoing EMP can reversibly switch between three classes of phenotypes: epithelial (E), mesenchymal (M), and hybrid E/M. While a large number of multistable regulatory networks have been identified to be driving EMP in various contexts, the exact mechanisms and design principles that enable robustness in driving EMP across contexts are not yet fully understood. Here, we investigated dynamic and structural robustness in EMP networks with regard to phenotypic heterogeneity and plasticity. We use two different approaches to simulate these networks: a computationally inexpensive, parameter-independent continuous state space Boolean model, and an ODE-based parameter-agnostic framework (RACIPE), both of which yielded similar phenotypic distributions. While the latter approach is useful for measurements of plasticity, the former model enabled us to extensively investigate robustness in phenotypic heterogeneity. Using perturbations to network topology and by varying network parameters, we show that multistable EMP networks are structurally and dynamically more robust compared with their randomized counterparts, thereby highlighting their topological hallmarks. These features of robustness are governed by a balance of positive and negative feedback loops embedded in these networks. Using a combination of the number of negative and positive feedback loops weighted by their lengths, we identified a metric that can explain the structural and dynamical robustness of these networks. This metric enabled us to compare networks across multiple sizes, and the network principles thus obtained can be used to identify fragilities in large networks without simulating their dynamics. Our analysis highlights a network topology-based approach to quantify robustness in the phenotypic heterogeneity and plasticity emergent from EMP networks.

Twist1 induces chromosomal instability (CIN) in colorectal cancer cells.

Twist1 is a basic helix-loop-helix transcription factor, essential during early development in mammals. While Twist1 induces epithelial-to-mesenchymal transition (EMT), here we show that Twist1 overexpression enhances nuclear and mitotic aberrations. This is accompanied by an increase in whole chromosomal copy number gains and losses, underscoring the role of Twist1 in inducing chromosomal instability (CIN) in colorectal cancer cells. Array comparative genomic hybridization (array CGH) analysis further shows sub-chromosomal deletions, consistent with an increased frequency of DNA double strand breaks (DSBs). Remarkably, Twist1 overexpression downmodulates key cell cycle checkpoint factors-Bub1, BubR1, Mad1 and Mad2-that regulate CIN. Mathematical simulations using the RACIPE tool show a negative correlation of Twist1 with E-cadherin and BubR1. Data analyses of gene expression profiles of patient samples from The Cancer Genome Atlas (TCGA) reveal a positive correlation between Twist1 and mesenchymal genes across cancers, whereas the correlation of TWIST1 with CIN and DSB genes is cancer subtype-specific. Taken together, these studies highlight the mechanistic involvement of Twist1 in the deregulation of factors that maintain genome stability during EMT in colorectal cancer cells. Twist1 overexpression enhances genome instability in the context of EMT that further contributes to cellular heterogeneity. In addition, these studies imply that Twist1 downmodulates nuclear lamins that further alter spatiotemporal organization of the cancer genome and epigenome. Notwithstanding their genetic background, colorectal cancer cells nevertheless maintain their overall ploidy, while the downstream effects of Twist1 enhance CIN and DNA damage enriching for sub-populations of aggressive cancer cells.

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 modeling, and highlights a set of questions that remain unanswered about the dynamics of hybrid E/M states.

Roadmap on plasticity and epigenetics in cancer.

The role of plasticity and epigenetics in shaping cancer evolution and response to therapy has taken center stage with recent technological advances including single cell sequencing. This roadmap article is focused on state-of-the-art mathematical and experimental approaches to interrogate plasticity in cancer, and addresses the following themes and questions: is there a formal overarching framework that encompasses both non-genetic plasticity and mutation-driven somatic evolution? How do we measure and model the role of the microenvironment in influencing/controlling non-genetic plasticity? How can we experimentally study non-genetic plasticity? Which mathematical techniques are required or best suited? What are the clinical and practical applications and implications of these concepts?