A mechanistic model captures the emergence and implications of non-genetic heterogeneity and reversible drug resistance in ER+ 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.

Self-consistent signal transduction analysis for modeling context-specific signaling cascades and perturbations

Biological signal transduction networks are central to information processing and regulation of gene expression across all domains of life. Dysregulation is known to cause a wide array of diseases, including cancers. Here I introduce self-consistent signal transduction analysis, which utilizes genome-scale -omics data (specifically transcriptomics and/or proteomics) in order to predict the flow of information through these networks in an individualized manner. I apply the method to the study of endocrine therapy in breast cancer patients, and show that drugs that inhibit estrogen receptor α elicit a wide array of antitumoral effects, and that their most clinically-impactful ones are through the modulation of proliferative signals that control the genes GREB1, HK1, AKT1, MAPK1, AKT2, and NQO1. This method offers researchers a valuable tool in understanding how and why dysregulation occurs, and how perturbations to the network (such as targeted therapies) effect the network itself, and ultimately patient outcomes.

A resource-based mechanistic framework for castration-resistant prostate cancer (CRPC)

Cancer therapy often leads to the selective elimination of drug-sensitive cells from the tumour. This can favour the growth of cells resistant to the therapeutic agent, ultimately causing a tumour relapse. Castration-resistant prostate cancer (CRPC) is a well-characterised instance of this phenomenon. In CRPC, after systemic androgen deprivation therapy (ADT), a subset of drug-resistant cancer cells autonomously produce testosterone, thus enabling tumour regrowth. A previous theoretical study has shown that such a tumour relapse can be delayed by inhibiting the growth of drug-resistant cells using biotic competition from drug-sensitive cells. In this context, the centrality of resource dynamics to intra-tumour competition in the CRPC system indicates clear scope for the construction of theoretical models that can explicitly incorporate the underlying mechanisms of tumour ecology. In the current study, we use a modified logistic framework to model cell-cell interactions in terms of the production and consumption of resources. Our results show that steady state composition of CRPC can be understood as a composite function of the availability and utilisation efficiency of two resources-oxygen and testosterone. In particular, we show that the effect of changing resource availability or use efficiency is conditioned by their general abundance regimes. Testosterone typically functions in trace amounts and thus affects steady state behaviour of the CRPC system differently from oxygen, which is usually available at higher levels. Our data thus indicate that explicit consideration of resource dynamics can produce novel and useful mechanistic understanding of CRPC. Furthermore, such a modelling approach also incorporates variables into the system's description that can be directly measured in a clinical context. This is therefore a promising avenue of research in cancer ecology that could lead to therapeutic approaches that are more clearly rooted in the biology of CRPC.

Semi-mechanistic modeling of resistance development to β-lactam and β-lactamase-inhibitor combinations

The use of β-lactam (BL) and β-lactamase inhibitor (BLI) combinations, such as piperacillin-tazobactam (PIP-TAZ) is an effective strategy to combat infections by extended-spectrum β-lactamase-producing bacteria. However, in Gram-negative bacteria, resistance (both mutational and adaptive) to BL-BLI combination can still develop through multiple mechanisms. These mechanisms may include increased β-lactamase activity, reduced drug influx, and increased drug efflux. Understanding the relative contribution of these mechanisms during resistance development helps identify the most impactful mechanism to target in designing a treatment to counter BL-BLI resistance. This study used semi-mechanistic mathematical modeling in combination with antibiotic sensitivity assays to assess the potential impact of different resistance mechanisms during the development of PIP-TAZ resistance in a Klebsiella pneumoniae isolate expressing CTX-M-15 and SHV-1 β-lactamases. The mathematical models were used to evaluate the potential impact of several cellular changes as a sole mediator of PIP-TAZ resistance. Our semi-mechanistic model identified 2 out of the 13 inspected mechanisms as key resistance mechanisms that may independently support the observed magnitude of PIP-TAZ resistance, namely porin loss and efflux pump up-regulation. Simulation using the resulting models also suggested the possible adjustment of PIP-TAZ dose outside its commonly used 8:1 dosing ratio. The current study demonstrated how theory-based mechanistic models informed by experimental data can be used to support hypothesis generation regarding potential resistance mechanisms, which may guide subsequent experimental studies.

Proneural-mesenchymal antagonism dominates the patterns of phenotypic heterogeneity in glioblastoma

The aggressive nature of glioblastoma (GBM) - one of the deadliest forms of brain tumors - is majorly attributed to underlying phenotypic heterogeneity. Early attempts to classify this heterogeneity at a transcriptomic level in TCGA GBM cohort proposed the existence of four distinct molecular subtypes: Proneural, Neural, Classical, and Mesenchymal. Further, a single-cell RNA sequencing (scRNA-seq) analysis of primary tumors also reported similar four subtypes mimicking neurodevelopmental lineages. However, it remains unclear whether these four subtypes identified via bulk and single-cell transcriptomics are mutually exclusive or not. Here, we perform pairwise correlations among individual genes and gene signatures corresponding to these proposed subtypes and show that the subtypes are not distinctly mutually antagonistic in either TCGA or scRNA-seq data. We observed that the proneural (or neural progenitor-like)-mesenchymal axis is the most prominent antagonistic pair, with the other two subtypes lying on this spectrum. These results are reinforced through a meta-analysis of over 100 single-cell and bulk transcriptomic datasets as well as in terms of functional association with metabolic switching, cell cycle, and immune evasion pathways. Finally, this proneural-mesenchymal antagonistic trend percolates to the association of relevant transcription factors with patient survival. These results suggest rethinking GBM phenotypic characterization for more effective therapeutic targeting efforts.

Transcriptional state dynamics lead to heterogeneity and adaptive tumor evolution in urothelial bladder carcinoma

Intra-tumor heterogeneity contributes to treatment failure and poor survival in urothelial bladder carcinoma (UBC). Analyzing transcriptome from a UBC cohort, we report that intra-tumor transcriptomic heterogeneity indicates co-existence of tumor cells in epithelial and mesenchymal-like transcriptional states and bi-directional transition between them occurs within and between tumor subclones. We model spontaneous and reversible transition between these partially heritable states in cell lines and characterize their population dynamics. SMAD3, KLF4 and PPARG emerge as key regulatory markers of the transcriptional dynamics. Nutrient limitation, as in the core of large tumors, and radiation treatment perturb the dynamics, initially selecting for a transiently resistant phenotype and then reconstituting heterogeneity and growth potential, driving adaptive evolution. Dominance of transcriptional states with low PPARG expression indicates an aggressive phenotype in UBC patients. We propose that phenotypic plasticity and dynamic, non-genetic intra-tumor heterogeneity modulate both the trajectory of disease progression and adaptive treatment response in UBC.

Cell population growth kinetics in the presence of stochastic heterogeneity of cell phenotype

Recent studies at individual cell resolution have revealed phenotypic heterogeneity in nominally clonal tumor cell populations. The heterogeneity affects cell growth behaviors, which can result in departure from the idealized uniform exponential growth of the cell population. Here we measured the stochastic time courses of growth of an ensemble of populations of HL60 leukemia cells in cultures, starting with distinct initial cell numbers to capture a departure from the uniform exponential growth model for the initial growth ("take-off"). Despite being derived from the same cell clone, we observed significant variations in the early growth patterns of individual cultures with statistically significant differences in growth dynamics, which could be explained by the presence of inter-converting subpopulations with different growth rates, and which could last for many generations. Based on the hypothesis of existence of multiple subpopulations, we developed a branching process model that was consistent with the experimental observations.

A systems-level analysis of the mutually antagonistic roles of RKIP and BACH1 in dynamics of cancer cell plasticity

Epithelial-mesenchymal transition (EMT) is an important axis of phenotypic plasticity-a hallmark of cancer metastasis. Raf kinase-B inhibitor protein (RKIP) and BTB and CNC homology 1 (BACH1) are reported to influence EMT. In breast cancer, they act antagonistically, but the exact nature of their roles in mediating EMT and associated other axes of plasticity remains unclear. Here, analysing transcriptomic data, we reveal their antagonistic trends in a pan-cancer manner in terms of association with EMT, metabolic reprogramming and immune evasion via PD-L1. Next, we developed and simulated a mechanism-based gene regulatory network that captures how RKIP and BACH1 engage in feedback loops with drivers of EMT and stemness. We found that RKIP and BACH1 belong to two antagonistic 'teams' of players-while BACH1 belonged to the one driving pro-EMT, stem-like and therapy-resistant cell states, RKIP belonged to the one enabling pro-epithelial, less stem-like and therapy-sensitive phenotypes. Finally, we observed that low RKIP levels and upregulated BACH1 levels associated with worse clinical outcomes in many cancer types. Together, our systems-level analysis indicates that the emergent dynamics of underlying regulatory network enable the antagonistic patterns of RKIP and BACH1 with various axes of cancer cell plasticity, and with patient survival data.

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.

A Nexus between Genetic and Non-Genetic Mechanisms Guides KRAS Inhibitor Resistance in Lung Cancer

Several studies in the last few years have determined that, in contrast to the prevailing dogma that drug resistance is simply due to Darwinian evolution-the selection of mutant clones in response to drug treatment-non-genetic changes can also lead to drug resistance whereby tolerant, reversible phenotypes are eventually relinquished by resistant, irreversible phenotypes. Here, using KRAS as a paradigm, we illustrate how this nexus between genetic and non-genetic mechanisms enables cancer cells to evade the harmful effects of drug treatment. We discuss how the conformational dynamics of the KRAS molecule, that includes intrinsically disordered regions, is influenced by the binding of the targeted therapies contributing to conformational noise and how this noise impacts the interaction of KRAS with partner proteins to rewire the protein interaction network. Thus, in response to drug treatment, reversible drug-tolerant phenotypes emerge via non-genetic mechanisms that eventually enable the emergence of irreversible resistant clones via genetic mutations. Furthermore, we also discuss the recent data demonstrating how combination therapy can help alleviate KRAS drug resistance in lung cancer, and how new treatment strategies based on evolutionary principles may help minimize or even preclude the emergence of drug resistance.

Toward a systems-level probing of tumor clonality

Cancer has been described as a genetic disease that clonally evolves in the face of selective pressures imposed by cell-intrinsic and extrinsic factors. Although classical models based on genetic data predominantly propose Darwinian mechanisms of cancer evolution, recent single-cell profiling of cancers has described unprecedented heterogeneity in tumors providing support for alternative models of branched and neutral evolution through both genetic and non-genetic mechanisms. Emerging evidence points to a complex interplay between genetic, non-genetic, and extrinsic environmental factors in shaping the evolution of tumors. In this perspective, we briefly discuss the role of cell-intrinsic and extrinsic factors that shape clonal behaviors during tumor progression, metastasis, and drug resistance. Taking examples of pre-malignant states associated with hematological malignancies and esophageal cancer, we discuss recent paradigms of tumor evolution and prospective approaches to further enhance our understanding of this spatiotemporally regulated process.

The role of tumor microenvironment in drug resistance: emerging technologies to unravel breast cancer heterogeneity

Breast cancer is a highly heterogeneous disease, at both inter- and intra-tumor levels, and this heterogeneity is a crucial determinant of malignant progression and response to treatments. In addition to genetic diversity and plasticity of cancer cells, the tumor microenvironment contributes to tumor heterogeneity shaping the physical and biological surroundings of the tumor. The activity of certain types of immune, endothelial or mesenchymal cells in the microenvironment can change the effectiveness of cancer therapies via a plethora of different mechanisms. Therefore, deciphering the interactions between the distinct cell types, their spatial organization and their specific contribution to tumor growth and drug sensitivity is still a major challenge. Dissecting intra-tumor heterogeneity is currently an urgent need to better define breast cancer biology and to develop therapeutic strategies targeting the microenvironment as helpful tools for combined and personalized treatment. In this review, we analyze the mechanisms by which the tumor microenvironment affects the characteristics of tumor heterogeneity that ultimately result in drug resistance, and we outline state of the art preclinical models and emerging technologies that will be instrumental in unraveling the impact of the tumor microenvironment on resistance to therapies.

Unraveling non-genetic heterogeneity in cancer with dynamical models and computational tools

Individual cells within an otherwise genetically homogenous population constantly undergo fluctuations in their molecular state, giving rise to non-genetic heterogeneity. Such diversity is being increasingly implicated in cancer therapy resistance and metastasis. Identifying the origins of non-genetic heterogeneity is therefore crucial for making clinical breakthroughs. We discuss with examples how dynamical models and computational tools have provided critical multiscale insights into the nature and consequences of non-genetic heterogeneity in cancer. We demonstrate how mechanistic modeling has been pivotal in establishing key concepts underlying non-genetic diversity at various biological scales, from population dynamics to gene regulatory networks. We discuss advances in single-cell longitudinal profiling techniques to reveal patterns of non-genetic heterogeneity, highlighting the ongoing efforts and challenges in statistical frameworks to robustly interpret such multimodal datasets. Moving forward, we stress the need for data-driven statistical and mechanistically motivated dynamical frameworks to come together to develop predictive cancer models and inform therapeutic strategies.

Emergent dynamics of underlying regulatory network links EMT and androgen receptor-dependent resistance in prostate cancer

Advanced prostate cancer patients initially respond to hormone therapy, be it in the form of androgen deprivation therapy or second-generation hormone therapies, such as abiraterone acetate or enzalutamide. However, most men with prostate cancer eventually develop hormone therapy resistance. This resistance can arise through multiple mechanisms, such as through genetic mutations, epigenetic mechanisms, or through non-genetic pathways, such as lineage plasticity along epithelial-mesenchymal or neuroendocrine-like axes. These mechanisms of hormone therapy resistance often co-exist within a single patient's tumor and can overlap within a single cell. There exists a growing need to better understand how phenotypic heterogeneity and plasticity results from emergent dynamics of the regulatory networks governing androgen independence. Here, we investigated the dynamics of a regulatory network connecting the drivers of androgen receptor (AR) splice variant-mediated androgen independence and those of epithelial-mesenchymal transition. Model simulations for this network revealed four possible phenotypes: epithelial-sensitive (ES), epithelial-resistant (ER), mesenchymal-resistant (MR) and mesenchymal-sensitive (MS), with the latter phenotype occurring rarely. We observed that well-coordinated "teams" of regulators working antagonistically within the network enable these phenotypes. These model predictions are supported by multiple transcriptomic datasets both at single-cell and bulk levels, including in vitro EMT induction models and clinical samples. Further, our simulations reveal spontaneous stochastic switching between the ES and MR states. Addition of the immune checkpoint molecule, PD-L1, to the network was able to capture the interactions between AR, PD-L1, and the mesenchymal marker SNAIL, which was also confirmed through quantitative experiments. This systems-level understanding of the driver of androgen independence and EMT could aid in understanding non-genetic transitions and progression of such cancers and help in identifying novel therapeutic strategies or targets.

Cell Population Growth Kinetics in the Presence of Stochastic Heterogeneity of Cell Phenotype

Recent studies at individual cell resolution have revealed phenotypic heterogeneity in nominally clonal tumor cell populations. The heterogeneity affects cell growth behaviors, which can result in departure from the idealized exponential growth. Here we measured the stochastic time courses of growth of an ensemble of populations of HL60 leukemia cells in cultures, starting with distinct initial cell numbers to capture the departure from the exponential growth model in the initial growth phase. Despite being derived from the same cell clone, we observed significant variations in the early growth patterns of individual cultures with statistically significant differences in growth kinetics and the presence of subpopulations with different growth rates that endured for many generations. Based on the hypothesis of existence of multiple inter-converting subpopulations, we developed a branching process model that captures the experimental observations.

"In medio stat virtus": Insights into hybrid E/M phenotype attitudes

Epithelial-mesenchymal plasticity (EMP) refers to the ability of cells to dynamically interconvert between epithelial (E) and mesenchymal (M) phenotypes, thus generating an array of hybrid E/M intermediates with mixed E and M features. Recent findings have demonstrated how these hybrid E/M rather than fully M cells play key roles in most of physiological and pathological processes involving EMT. To this regard, the onset of hybrid E/M state coincides with the highest stemness gene expression and is involved in differentiation of either normal and cancer stem cells. Moreover, hybrid E/M cells are responsible for wound healing and create a favorable immunosuppressive environment for tissue regeneration. Nevertheless, hybrid state is responsible of metastatic process and of the increasing of survival, apoptosis and therapy resistance in cancer cells. The present review aims to describe the main features and the emerging concepts regulating EMP and the formation of E/M hybrid intermediates by describing differences and similarities between cancer and normal hybrid stem cells. In particular, the comprehension of hybrid E/M cells biology will surely advance our understanding of their features and how they could be exploited to improve tissue regeneration and repair.

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.

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.

Tumor Hybrid Cells: Nature and Biological Significance

Metastasis is the leading cause of cancer death and can be realized through the phenomenon of tumor cell fusion. The fusion of tumor cells with other tumor or normal cells leads to the appearance of tumor hybrid cells (THCs) exhibiting novel properties such as increased proliferation and migration, drug resistance, decreased apoptosis rate, and avoiding immune surveillance. Experimental studies showed the association of THCs with a high frequency of cancer metastasis; however, the underlying mechanisms remain unclear. Many other questions also remain to be answered: the role of genetic alterations in tumor cell fusion, the molecular landscape of cells after fusion, the lifetime and fate of different THCs, and the specific markers of THCs, and their correlation with various cancers and clinicopathological parameters. In this review, we discuss the factors and potential mechanisms involved in the occurrence of THCs, the types of THCs, and their role in cancer drug resistance and metastasis, as well as potential therapeutic approaches for the prevention, and targeting of tumor cell fusion. In conclusion, we emphasize the current knowledge gaps in the biology of THCs that should be addressed to develop highly effective therapeutics and strategies for metastasis suppression.

Intrinsically disordered proteins: Ensembles at the limits of Anfinsen's dogma

Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure. Hence, they are often misconceived to present a challenge to Anfinsen's dogma. However, IDPs exist as ensembles that sample a quasi-continuum of rapidly interconverting conformations and, as such, may represent proteins at the extreme limit of the Anfinsen postulate. IDPs play important biological roles and are key components of the cellular protein interaction network (PIN). Many IDPs can interconvert between disordered and ordered states as they bind to appropriate partners. Conformational dynamics of IDPs contribute to conformational noise in the cell. Thus, the dysregulation of IDPs contributes to increased noise and "promiscuous" interactions. This leads to PIN rewiring to output an appropriate response underscoring the critical role of IDPs in cellular decision making. Nonetheless, IDPs are not easily tractable experimentally. Furthermore, in the absence of a reference conformation, discerning the energy landscape representation of the weakly funneled IDPs in terms of reaction coordinates is challenging. To understand conformational dynamics in real time and decipher how IDPs recognize multiple binding partners with high specificity, several sophisticated knowledge-based and physics-based in silico sampling techniques have been developed. Here, using specific examples, we highlight recent advances in energy landscape visualization and molecular dynamics simulations to discern conformational dynamics and discuss how the conformational preferences of IDPs modulate their function, especially in phenotypic switching. Finally, we discuss recent progress in identifying small molecules targeting IDPs underscoring the potential therapeutic value of IDPs. Understanding structure and function of IDPs can not only provide new insight on cellular decision making but may also help to refine and extend Anfinsen's structure/function paradigm.

Quantifying the Patterns of Metabolic Plasticity and Heterogeneity along the Epithelial–Hybrid–Mesenchymal Spectrum in Cancer

Cancer metastasis is the leading cause of cancer-related mortality and the process of the epithelial-to-mesenchymal transition (EMT) is crucial for cancer metastasis. Both partial and complete EMT have been reported to influence the metabolic plasticity of cancer cells in terms of switching among the oxidative phosphorylation, fatty acid oxidation and glycolysis pathways. However, a comprehensive analysis of these major metabolic pathways and their associations with EMT across different cancers is lacking. Here, we analyse more than 180 cancer cell datasets and show the diverse associations of these metabolic pathways with the EMT status of cancer cells. Our bulk data analysis shows that EMT generally positively correlates with glycolysis but negatively with oxidative phosphorylation and fatty acid metabolism. These correlations are also consistent at the level of their molecular master regulators, namely AMPK and HIF1α. Yet, these associations are shown to not be universal. The analysis of single-cell data for EMT induction shows dynamic changes along the different axes of metabolic pathways, consistent with general trends seen in bulk samples. Further, assessing the association of EMT and metabolic activity with patient survival shows that a higher extent of EMT and glycolysis predicts a worse prognosis in many cancers. Together, our results reveal the underlying patterns of metabolic plasticity and heterogeneity as cancer cells traverse through the epithelial–hybrid–mesenchymal spectrum of 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?

NRF2-dependent Epigenetic Regulation can Promote the Hybrid Epithelial/Mesenchymal Phenotype

The epithelial-mesenchymal transition (EMT) is a cellular process critical for wound healing, cancer metastasis and embryonic development. Recent efforts have identified the role of hybrid epithelial/mesenchymal states, having both epithelial and mesehncymal traits, in enabling cancer metastasis and resistance to various therapies. Also, previous work has suggested that NRF2 can act as phenotypic stability factor to help stablize such hybrid states. Here, we incorporate a phenomenological epigenetic feedback effect into our previous computational model for EMT signaling. We show that this type of feedback can stabilize the hybrid state as compared to the fully mesenchymal phenotype if NRF2 can influence SNAIL at an epigenetic level, as this link makes transitions out of hybrid state more difficult. However, epigenetic regulation on other NRF2-related links do not significantly change the EMT dynamics. Finally, we considered possible cell division effects in our epigenetic regulation model, and our results indicate that the degree of epigenetic inheritance does not appear to be a critical factor for the hybrid E/M state stabilizing behavior of NRF2.

Immunosuppressive Traits of the Hybrid Epithelial/Mesenchymal Phenotype

Recent preclinical and clinical data suggests enhanced metastatic fitness of hybrid epithelial/mesenchymal (E/M) phenotypes, but mechanistic details regarding their survival strategies during metastasis remain unclear. Here, we investigate immune-evasive strategies of hybrid E/M states. We construct and simulate the dynamics of a minimalistic regulatory network encompassing the known associations among regulators of EMT (epithelial-mesenchymal transition) and PD-L1, an established immune-suppressor. Our simulations for the network consisting of SLUG, ZEB1, miR-200, CDH1 and PD-L1, integrated with single-cell and bulk RNA-seq data analysis, elucidate that hybrid E/M cells can have high levels of PD-L1, similar to those seen in cells with a full EMT phenotype, thus obviating the need for cancer cells to undergo a full EMT to be immune-evasive. Specifically, in breast cancer, we show the co-existence of hybrid E/M phenotypes, enhanced resistance to anti-estrogen therapy and increased PD-L1 levels. Our results underscore how the emergent dynamics of interconnected regulatory networks can coordinate different axes of cellular fitness during metastasis.

Emerging roles of epithelial-mesenchymal plasticity in invasion-metastasis cascade and therapy resistance

Strong association of cancer incidence and its progression with mortality highlights the need to decipher the cellular and molecular mechanisms that drive tumor cells to rapidly progress to metastatic disease and therapy resistance. Epithelial-mesenchymal plasticity (EMP) emerged as a key regulator of metastatic outgrowth. It allows neoplastic epithelial cells to delaminate from their neighbors either individually or collectively, traverse the extracellular matrix (ECM) barrier, enter into the circulation, and establish distal metastases. Plasticity between epithelial and mesenchymal states and the existence of hybrid epithelial/mesenchymal (E/M) phenotypes are increasingly being reported in different tumor contexts. Small subset of cancer cells with stemness called cancer stem cells (CSCs) exhibit plasticity, possess high tumorigenic potential, and contribute to high degree of tumoral heterogeneity. EMP characterized by the presence of dynamic intermediate states is reported to be influenced by (epi)genomic reprograming, growth factor signaling, inflammation, and low oxygen generated by tumor stromal microenvironment. EMP alters the genotypic and phenotypic characteristics of tumor cells/CSCs, disrupts tissue homeostasis, induces the reprogramming of angiogenic and immune recognition functions, and renders tumor cells to survive hostile microenvironments and resist therapy. The present review summarizes the roles of EMP in tumor invasion and metastasis and provides an update on therapeutic strategies to target the metastatic and refractory cancers.

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.

Protein conformational dynamics and phenotypic switching

Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure but exist as conformational ensembles. Because of their structural plasticity, they can interact with multiple partners. The protein interactions between IDPs and their partners form scale-free protein interaction networks (PINs) that facilitate information flow in the cell. Because of their plasticity, IDPs typically occupy hub positions in cellular PINs. Furthermore, their conformational dynamics and propensity for post-translational modifications contribute to "conformational" noise which is distinct from the well-recognized transcriptional noise. Therefore, upregulation of IDPs in response to a specific input, such as stress, contributes to increased noise and, hence, an increase in stochastic, "promiscuous" interactions. These interactions lead to activation of latent pathways or can induce "rewiring" of the PIN to yield an optimal output underscoring the critical role of IDPs in regulating information flow. We have used PAGE4, a highly intrinsically disordered stress-response protein as a paradigm. Employing a variety of experimental and computational techniques, we have elucidated the role of PAGE4 in phenotypic switching of prostate cancer cells at a systems level. These cumulative studies over the past decade provide a conceptual framework to better understand how IDP conformational dynamics and conformational noise might facilitate cellular decision-making.

Impact of Epithelial-Mesenchymal Transition on the Immune Landscape in Breast Cancer

The impact of epithelial-mesenchymal transition (EMT) signature on the immune infiltrate present in the breast cancer tumor microenvironment (TME) is still poorly understood. Since there is mounting interest in the use of immunotherapy for the treatment of subsets of breast cancer patients, it is of major importance to understand the fundamental tumor characteristics which dictate the inter-tumor heterogeneity in immune landscapes. We aimed to assess the impact of EMT-related markers on the nature and magnitude of the inflammatory infiltrate present in breast cancer TME and their association with the clinicopathological parameters. Tissue microarrays were constructed from 144 formalin-fixed paraffin-embedded invasive breast cancer tumor samples. The protein expression patterns of Snail, Twist, ZEB1, N-cadherin, Vimentin, GRHL2, E-cadherin, and EpCAM were examined by immunohistochemistry (IHC). The inflammatory infiltrate in the TME was assessed semi-quantitatively on hematoxylin and eosin (H&E)-stained whole sections and was characterized using IHC. The inflammatory infiltrate was more intense in poorly differentiated carcinomas and triple-negative carcinomas in which the expression of E-cadherin and GRHL2 was reduced, while EpCAM was overexpressed. Most EMT-related markers correlated with plasma cell infiltration of the TME. Taken together, our findings reveal that the EMT signature might impact the immune response in the TME.

The role of memory in non-genetic inheritance and its impact on cancer treatment resistance

Intra-tumour heterogeneity is a leading cause of treatment failure and disease progression in cancer. While genetic mutations have long been accepted as a primary mechanism of generating this heterogeneity, the role of phenotypic plasticity is becoming increasingly apparent as a driver of intra-tumour heterogeneity. Consequently, understanding the role of this plasticity in treatment resistance and failure is a key component of improving cancer therapy. We develop a mathematical model of stochastic phenotype switching that tracks the evolution of drug-sensitive and drug-tolerant subpopulations to clarify the role of phenotype switching on population growth rates and tumour persistence. By including cytotoxic therapy in the model, we show that, depending on the strategy of the drug-tolerant subpopulation, stochastic phenotype switching can lead to either transient or permanent drug resistance. We study the role of phenotypic heterogeneity in a drug-resistant, genetically homogeneous population of non-small cell lung cancer cells to derive a rational treatment schedule that drives population extinction and avoids competitive release of the drug-tolerant sub-population. This model-informed therapeutic schedule results in increased treatment efficacy when compared against periodic therapy, and, most importantly, sustained tumour decay without the development of resistance.

Lineage Plasticity in Cancer: The Tale of a Skin-Walker

Lineage plasticity, the switching of cells from one lineage to another, has been recognized as a cardinal property essential for embryonic development, tissue repair and homeostasis. However, such a highly regulated process goes awry when cancer cells exploit this inherent ability to their advantage, resulting in tumorigenesis, relapse, metastasis and therapy resistance. In this review, we summarize our current understanding on the role of lineage plasticity in tumor progression and therapeutic resistance in multiple cancers. Lineage plasticity can be triggered by treatment itself and is reported across various solid as well as liquid tumors. Here, we focus on the importance of lineage switching in tumor progression and therapeutic resistance of solid tumors such as the prostate, lung, hepatocellular and colorectal carcinoma and the myeloid and lymphoid lineage switch observed in leukemias. Besides this, we also discuss the role of epithelial-mesenchymal transition (EMT) in facilitating the lineage switch in biphasic cancers such as aggressive carcinosarcomas. We also discuss the mechanisms involved, current therapeutic approaches and challenges that lie ahead in taming the scourge of lineage plasticity in cancer.