Mapping Phenotypic Plasticity upon the Cancer Cell State Landscape Using Manifold Learning

Intercellular signaling reinforces single-cell level phenotypic transitions and facilitates robust re-equilibrium of heterogeneous cancer cell populations

Background

Cancer cells within tumors exhibit a wide range of phenotypic states driven by non-genetic mechanisms in addition to extensively studied genetic alterations. Conversions among cancer cell states can result in intratumoral heterogeneity which contributes to metastasis and development of drug resistance. However, mechanisms underlying the initiation and/or maintenance of such phenotypic plasticity are poorly understood. In particular, the role of intercellular communications in phenotypic plasticity remains elusive.

Methods

In this study, we employ a multiscale inference-based approach using single-cell RNA sequencing (scRNA-seq) data to explore how intercellular interactions influence phenotypic dynamics of cancer cells, particularly cancers undergoing epithelial-mesenchymal transition. In addition, we use mathematical models based on our data-driven findings to interrogate the roles of intercellular communications at the cell populations from the viewpoint of dynamical systems.

Results

Our inference approach reveals that signaling interactions between cancerous cells in small cell lung cancer (SCLC) result in the reinforcement of the phenotypic transition in single cells and the maintenance of population-level intratumoral heterogeneity. Additionally, we find a recurring signaling pattern across multiple types of cancer in which the mesenchymal-like subtypes utilize signals from other subtypes to support its new phenotype, further promoting the intratumoral heterogeneity. Our models show that inter-subtype communication both accelerates the development of heterogeneous tumor populations and confers robustness to their steady state phenotypic compositions.

Conclusions

Our work highlights the critical role of intercellular signaling in sustaining intratumoral heterogeneity, and our approach of computational analysis of scRNA-seq data can infer inter- and intra-cellular signaling networks in a holistic manner.

Tumor dormancy and relapse: understanding the molecular mechanisms of cancer recurrence

Cancer recurrence, driven by the phenomenon of tumor dormancy, presents a formidable challenge in oncology. Dormant cancer cells have the ability to evade detection and treatment, leading to relapse. This review emphasizes the urgent need to comprehend tumor dormancy and its implications for cancer recurrence. Despite notable advancements, significant gaps remain in our understanding of the mechanisms underlying dormancy and the lack of reliable biomarkers for predicting relapse. This review provides a comprehensive analysis of the cellular, angiogenic, and immunological aspects of dormancy. It highlights the current therapeutic strategies targeting dormant cells, particularly combination therapies and immunotherapies, which hold promise in preventing relapse. By elucidating these mechanisms and proposing innovative research methodologies, this review aims to deepen our understanding of tumor dormancy, ultimately facilitating the development of more effective strategies for preventing cancer recurrence and improving patient outcomes.

Design principles of regulatory networks underlying epithelial mesenchymal plasticity in cancer cells

Phenotypic plasticity is a hallmark of cancer and drives metastatic disease and drug resistance. The dynamics of epithelial mesenchymal plasticity is driven by complex interactions involving multiple feedback loops in underlying networks operating at multiple regulatory levels such as transcriptional and epigenetic. The past decade has witnessed a surge in systems level analysis of structural and dynamical traits of these networks. Here, we highlight the key insights elucidated from such efforts-a) multistability in gene regulatory networks and the co-existence of many hybrid phenotypes, thus enabling a landscape with multiple 'attractors', b) mutually antagonistic 'teams' of genes in these networks, shaping the rates of cell state transition in this landscape, and c) chromatin level changes that can alter the landscape, thus controlling reversibility of cell state transitions, allowing cellular memory in the context of epithelial mesenchymal plasticity in cancer cells. Such approaches, in close integration with high-throughput longitudinal data, have improved our understanding of the dynamics of cell state transitions implicated in tumor cell plasticity.

Knockdown of WISP1/DKK1 restrains phenotypic plasticity in esophageal squamous cell carcinoma by suppressing epithelial-mesenchymal transition and stemness

OBJECTIVE

Wnt-induced signaling protein 1 (WISP1) and Dickkopf-1 (DKK1) are highly expressed in esophageal squamous cell carcinoma (ESCC), but no direct connection was identified between them. Phenotypic plasticity is a hallmark of ESCC. This research intended to identify the association between WISP1 and DKK1 and their roles in the phenotypic plasticity of ESCC.

METHODS

Genes differentially expressed in esophageal carcinoma were analyzed in the GEO database, followed by analyses of GO and KEGG enrichment to screen the hub gene. WISP1 expression and DKK1 secretion was assessed in ESCC tissues and cells. The tumor xenograft and in vivo metastasis models were established by injecting ESCC cells into nude mice. Functional deficiency and rescue experiments were conducted, followed by assays for cell proliferation, migration/invasion, stemness, epithelial-mesenchymal transition (EMT), and apoptosis, as well as tumor volume, weight, proliferation, stemness, and lung metastasis. The binding relationship and co-expression of WISP1 and DKK1 were determined.

RESULTS

WISP1 and DKK1 were upregulated in ESCC cells and tissues, and WISP1 was enriched in the cell stemness and Wnt pathways. WISP1 knockdown subdued proliferation, migration/invasion, EMT activity, and stemness but enhanced apoptosis in ESCC cells. WISP1 knockdown restrained ESCC growth, proliferation, stemness, and metastasis in vivo. WISP1 bound to DKK1 in ESCC. DKK1 overexpression abolished the repressive impacts of WISP1 knockdown on the malignant behaviors of ESCC cells in vitro and of ESCC tumor in vivo.

CONCLUSION

Knockdown of WISP1/DKK1 restrains the phenotypic plasticity in esophageal squamous cell carcinoma by suppressing epithelial-mesenchymal transition and stemness.

Chromatin-based memory as a self-stabilizing influence on cell identity

Cell types are traditionally thought to be specified and stabilized by gene regulatory networks. Here, we explore how chromatin memory contributes to the specification and stabilization of cell states. Through pervasive, local, feedback loops, chromatin memory enables cell states that were initially unstable to become stable. Deeper appreciation of this self-stabilizing role for chromatin broadens our perspective of Waddington's epigenetic landscape from a static surface with islands of stability shaped by evolution, to a plasticine surface molded by experience. With implications for the evolution of cell types, stabilization of resistant states in cancer, and the widespread plasticity of complex life.

Heterogeneity-driven phenotypic plasticity and treatment response in branched-organoid models of pancreatic ductal adenocarcinoma

In patients with pancreatic ductal adenocarcinoma (PDAC), intratumoural and intertumoural heterogeneity increases chemoresistance and mortality rates. However, such morphological and phenotypic diversities are not typically captured by organoid models of PDAC. Here we show that branched organoids embedded in collagen gels can recapitulate the phenotypic landscape seen in murine and human PDAC, that the pronounced molecular and morphological intratumoural and intertumoural heterogeneity of organoids is governed by defined transcriptional programmes (notably, epithelial-to-mesenchymal plasticity), and that different organoid phenotypes represent distinct tumour-cell states with unique biological features in vivo. We also show that phenotype-specific therapeutic vulnerabilities and modes of treatment-induced phenotype reprogramming can be captured in phenotypic heterogeneity maps. Our methodology and analyses of tumour-cell heterogeneity in PDAC may guide the development of phenotype-targeted treatment strategies.

Serial single-cell RNA sequencing unveils drug resistance and metastatic traits in stage IV breast cancer

Metastasis is a complex process that remains poorly understood at the molecular levels. We profiled single-cell transcriptomic, genomic, and epigenomic changes associated with cancer cell progression, chemotherapy resistance, and metastasis from a Stage IV breast cancer patient. Pretreatment- and posttreatment-specimens from the primary tumor and distant metastases were collected for single-cell RNA sequencing and subsequent cell clustering, copy number variation (CNV) estimation, transcriptomic factor estimation, and pseudotime analyses. CNV analysis revealed that a small population of pretreatment cancer cells resisted chemotherapy and expanded. New clones including Metastatic Precursor Cells (MPCs), emerged in the posttreatment primary tumors in CNV similar to metastatic cells. MPCs exhibited expression profiles indicative of epithelial-mesenchymal transition. Comparison of MPCs with metastatic cancer cells also revealed dynamic changes in transcription factors and calcitonin pathway gene expression. These findings demonstrate the utility of single-patient clinical sample analysis for understanding tumor drug resistance, regrowth, and metastasis.

Inertial effect of cell state velocity on the quiescence-proliferation fate decision

Energy landscapes can provide intuitive depictions of population heterogeneity and dynamics. However, it is unclear whether individual cell behavior, hypothesized to be determined by initial position and noise, is faithfully recapitulated. Using the p21-/Cdk2-dependent quiescence-proliferation decision in breast cancer dormancy as a testbed, we examined single-cell dynamics on the landscape when perturbed by hypoxia, a dormancy-inducing stress. Combining trajectory-based energy landscape generation with single-cell time-lapse microscopy, we found that a combination of initial position and velocity on a p21/Cdk2 landscape, but not position alone, was required to explain the observed cell fate heterogeneity under hypoxia. This is likely due to additional cell state information such as epigenetic features and/or other species encoded in velocity but missing in instantaneous position determined by p21 and Cdk2 levels alone. Here, velocity dependence manifested as inertia: cells with higher cell cycle velocities prior to hypoxia continued progressing along the cell cycle under hypoxia, resisting the change in landscape towards cell cycle exit. Such inertial effects may markedly influence cell fate trajectories in tumors and other dynamically changing microenvironments where cell state transitions are governed by coordination across several biochemical species.

Modeling tumors as complex ecosystems

Many cancers resist therapeutic intervention. This is fundamentally related to intratumor heterogeneity: multiple cell populations, each with different phenotypic signatures, coexist within a tumor and its metastases. Like species in an ecosystem, cancer populations are intertwined in a complex network of ecological interactions. Most mathematical models of tumor ecology, however, cannot account for such phenotypic diversity or predict its consequences. Here, we propose that the generalized Lotka-Volterra model (GLV), a standard tool to describe species-rich ecological communities, provides a suitable framework to model the ecology of heterogeneous tumors. We develop a GLV model of tumor growth and discuss how its emerging properties provide a new understanding of the disease. We discuss potential extensions of the model and their application to phenotypic plasticity, cancer-immune interactions, and metastatic growth. Our work outlines a set of questions and a road map for further research in cancer ecology.

Cancer plasticity in therapy resistance: Mechanisms and novel strategies

Therapy resistance poses a significant obstacle to effective cancer treatment. Recent insights into cell plasticity as a new paradigm for understanding resistance to treatment: as cancer progresses, cancer cells experience phenotypic and molecular alterations, corporately known as cell plasticity. These alterations are caused by microenvironment factors, stochastic genetic and epigenetic changes, and/or selective pressure engendered by treatment, resulting in tumor heterogeneity and therapy resistance. Increasing evidence suggests that cancer cells display remarkable intrinsic plasticity and reversibly adapt to dynamic microenvironment conditions. Dynamic interactions between cell states and with the surrounding microenvironment form a flexible tumor ecosystem, which is able to quickly adapt to external pressure, especially treatment. Here, this review delineates the formation of cancer cell plasticity (CCP) as well as its manipulation of cancer escape from treatment. Furthermore, the intrinsic and extrinsic mechanisms driving CCP that promote the development of therapy resistance is summarized. Novel treatment strategies, e.g., inhibiting or reversing CCP is also proposed. Moreover, the review discusses the multiple lines of ongoing clinical trials globally aimed at ameliorating therapy resistance. Such advances provide directions for the development of new treatment modalities and combination therapies against CCP in the context of therapy resistance.

Chromatin as an old and new anticancer target

Recent genome-wide analyses identified chromatin modifiers as one of the most frequently mutated classes of genes across all cancers. However, chemotherapies developed for cancers involving DNA damage remain the standard of care for chromatin-deranged malignancies. In this review we address this conundrum by establishing the concept of 'chromatin damage': the non-genetic damage to protein-DNA interactions induced by certain small molecules. We highlight anthracyclines, a class of chemotherapeutic agents ubiquitously applied in oncology, as an example of overlooked chromatin-targeting agents. We discuss our current understanding of this phenomenon and explore emerging chromatin-damaging agents as a basis for further studies to maximize their impact in modern cancer treatment.

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.

Assessing RNA-Seq Workflow Methodologies Using Shannon Entropy

RNA-seq faces persistent challenges due to the ongoing, expanding array of data processing workflows, none of which have yet achieved standardization to date. It is imperative to determine which method most effectively preserves biological facts. Here, we used Shannon entropy as a tool for depicting the biological status of a system. Thus, we assessed the measurement of Shannon entropy by several RNA-seq workflow approaches, such as DESeq2 and edgeR, but also by combining nine normalization methods with log2 fold change on paired samples of TCGA RNA-seq representing datasets of 515 patients and spanning 12 different cancer types with 5-year overall survival rates ranging from 20% to 98%. Our analysis revealed that TPM, RLE, and TMM normalization, coupled with a threshold of log2 fold change ≥1, for identifying differentially expressed genes, yielded the best results. We propose that Shannon entropy can serve as an objective metric for refining the optimization of RNA-seq workflows and mRNA sequencing technologies.

Single-cell multiomics profiling reveals heterogeneous transcriptional programs and microenvironment in DSRCTs

Desmoplastic small round cell tumor (DSRCT) is a rare, aggressive sarcoma driven by the EWSR1::WT1 chimeric transcription factor. Despite this unique oncogenic driver, DSRCT displays a polyphenotypic differentiation of unknown causality. Using single-cell multi-omics on 12 samples from five patients, we find that DSRCT tumor cells cluster into consistent subpopulations with partially overlapping lineage- and metabolism-related transcriptional programs. In vitro modeling shows that high EWSR1::WT1 DNA-binding activity associates with most lineage-related states, in contrast to glycolytic and profibrotic states. Single-cell chromatin accessibility analysis suggests that EWSR1::WT1 binding site variability may drive distinct lineage-related transcriptional programs, supporting some level of cell-intrinsic plasticity. Spatial transcriptomics reveals that glycolytic and profibrotic states specifically localize within hypoxic niches at the periphery of tumor cell islets, suggesting an additional role of tumor cell-extrinsic microenvironmental cues. We finally identify a single-cell transcriptomics-derived epithelial signature associated with improved patient survival, highlighting the clinical relevance of our findings.

Phenotypic noise and plasticity in cancer evolution

Non-genetic alterations can produce changes in a cell's phenotype. In cancer, these phenomena can influence a cell's fitness by conferring access to heritable, beneficial phenotypes. Herein, we argue that current discussions of 'phenotypic plasticity' in cancer evolution ignore a salient feature of the original definition: namely, that it occurs in response to an environmental change. We suggest 'phenotypic noise' be used to distinguish non-genetic changes in phenotype that occur independently from the environment. We discuss the conceptual and methodological techniques used to identify these phenomena during cancer evolution. We propose that the distinction will guide efforts to define mechanisms of phenotype change, accelerate translational work to manipulate phenotypes through treatment, and, ultimately, improve patient outcomes.

Modeling tumors as species-rich ecological communities

Many advanced cancers resist therapeutic intervention. This process is fundamentally related to intra-tumor heterogeneity: multiple cell populations, each with different mutational and phenotypic signatures, coexist within a tumor and its metastatic nodes. Like species in an ecosystem, many cancer cell populations are intertwined in a complex network of ecological interactions. Most mathematical models of tumor ecology, however, cannot account for such phenotypic diversity nor are able to predict its consequences. Here we propose that the Generalized Lotka-Volterra model (GLV), a standard tool to describe complex, species-rich ecological communities, provides a suitable framework to describe the ecology of heterogeneous tumors. We develop a GLV model of tumor growth and discuss how its emerging properties, such as outgrowth and multistability, provide a new understanding of the disease. Additionally, we discuss potential extensions of the model and their application to three active areas of cancer research, namely phenotypic plasticity, the cancer-immune interplay and the resistance of metastatic tumors to treatment. Our work outlines a set of questions and a tentative road map for further research in cancer ecology.

Quantifying Landscape-Flux via Single-Cell Transcriptomics Uncovers the Underlying Mechanism of Cell Cycle

Recent developments in single-cell sequencing technology enable the acquisition of entire transcriptome data. Understanding the underlying mechanism and identifying the driving force of transcriptional regulation governing cell function directly from these data remains challenging. This study reconstructs a continuous vector field of the cell cycle based on discrete single-cell RNA velocity to quantify the single-cell global nonequilibrium dynamic landscape-flux. It reveals that large fluctuations disrupt the global landscape and genetic perturbations alter landscape-flux, thus identifying key genes in maintaining cell cycle dynamics and predicting associated functional effects. Additionally, it quantifies the fundamental energy cost of the cell cycle initiation and unveils that sustaining the cell cycle requires curl flux and dissipation to maintain the oscillatory phase coherence. This study enables the inference of the cell cycle gene regulatory networks directly from the single-cell transcriptomic data, including the feedback mechanisms and interaction intensity. This provides a golden opportunity to experimentally verify the landscape-flux theory and also obtain its associated quantifications. It also offers a unique framework for combining the landscape-flux theory and single-cell high-through sequencing experiments for understanding the underlying mechanisms of the cell cycle and can be extended to other nonequilibrium biological processes, such as differentiation development and disease pathogenesis.

Deciphering cancer cell state plasticity with single-cell genomics and artificial intelligence

Cancer stem cell plasticity refers to the ability of tumour cells to dynamically switch between states-for example, from cancer stem cells to non-cancer stem cell states. Governed by regulatory processes, cells transition through a continuum, with this transition space often referred to as a cell state landscape. Plasticity in cancer cell states leads to divergent biological behaviours, with certain cell states, or state transitions, responsible for tumour progression and therapeutic response. The advent of single-cell assays means these features can now be measured for individual cancer cells and at scale. However, the high dimensionality of this data, complex relationships between genomic features, and a lack of precise knowledge of the genomic profiles defining cancer cell states have opened the door for artificial intelligence methods for depicting cancer cell state landscapes. The contribution of cell state plasticity to cancer phenotypes such as treatment resistance, metastasis, and dormancy has been masked by analysis of 'bulk' genomic data-constituted of the average signal from millions of cells. Single-cell technologies solve this problem by producing a high-dimensional cellular landscape of the tumour ecosystem, quantifying the genomic profiles of individual cells, and creating a more detailed model to investigate cancer plasticity (Genome Res 31:1719, 2021; Semin Cancer Biol 53: 48-58, 2018; Signal Transduct Target Ther 5:1-36, 2020). In conjunction, rapid development in artificial intelligence methods has led to numerous tools that can be employed to study cancer cell plasticity.

How important is EMT for cancer metastasis?

Epithelial-to-mesenchymal transition (EMT), a biological phenomenon of cellular plasticity initially reported in embryonic development, has been increasingly recognized for its importance in cancer progression and metastasis. Despite tremendous progress being made in the past 2 decades in our understanding of the molecular mechanism and functional importance of EMT in cancer, there are several mysteries around EMT that remain unresolved. In this Unsolved Mystery, we focus on the variety of EMT types in metastasis, cooperative and collective EMT behaviors, spatiotemporal characterization of EMT, and strategies of therapeutically targeting EMT. We also highlight new technical advances that will facilitate the efforts to elucidate the unsolved mysteries of EMT in metastasis.

KDM5B predicts temozolomide-resistant subclones in glioblastoma

Adaptive plasticity to the standard chemotherapeutic temozolomide (TMZ) leads to glioblastoma progression. Here, we examine early stages of this process in patient-derived cellular models, exposing the human lysine-specific demethylase 5B (KDM5B) as a prospective indicator for subclonal expansion. By integration of a reporter, we show its preferential activity in rare, stem-like ALDH1A1+ cells, immediately increasing expression upon TMZ exposure. Naive, genetically unmodified KDM5Bhigh cells phosphorylate AKT (pAKT) and act as slow-cycling persisters under TMZ. Knockdown of KDM5B reverses pAKT levels, simultaneously increasing PTEN expression and TMZ sensitivity. Pharmacological inhibition of PTEN rescues the effect. Interference with KDM5B subsequent to TMZ decreases cellular vitality, and clonal tracing with DNA barcoding demonstrates high individual levels of KDM5B to predict subclonal expansion already before TMZ exposure. Thus, KDM5Bhigh treatment-naive cells preferentially contribute to the dynamics of drug resistance under TMZ. These findings may serve as a cornerstone for future biomarker-assisted clinical trials.

Trellis tree-based analysis reveals stromal regulation of patient-derived organoid drug responses

Patient-derived organoids (PDOs) can model personalized therapy responses; however, current screening technologies cannot reveal drug response mechanisms or how tumor microenvironment cells alter therapeutic performance. To address this, we developed a highly multiplexed mass cytometry platform to measure post-translational modification (PTM) signaling, DNA damage, cell-cycle activity, and apoptosis in >2,500 colorectal cancer (CRC) PDOs and cancer-associated fibroblasts (CAFs) in response to clinical therapies at single-cell resolution. To compare patient- and microenvironment-specific drug responses in thousands of single-cell datasets, we developed "Trellis"-a highly scalable, tree-based treatment effect analysis method. Trellis single-cell screening revealed that on-target cell-cycle blockage and DNA-damage drug effects are common, even in chemorefractory PDOs. However, drug-induced apoptosis is rarer, patient-specific, and aligns with cancer cell PTM signaling. We find that CAFs can regulate PDO plasticity-shifting proliferative colonic stem cells (proCSCs) to slow-cycling revival colonic stem cells (revCSCs) to protect cancer cells from chemotherapy.

An oncogenic phenoscape of colonic stem cell polarization

Cancer cells are regulated by oncogenic mutations and microenvironmental signals, yet these processes are often studied separately. To functionally map how cell-intrinsic and cell-extrinsic cues co-regulate cell fate, we performed a systematic single-cell analysis of 1,107 colonic organoid cultures regulated by (1) colorectal cancer (CRC) oncogenic mutations, (2) microenvironmental fibroblasts and macrophages, (3) stromal ligands, and (4) signaling inhibitors. Multiplexed single-cell analysis revealed a stepwise epithelial differentiation phenoscape dictated by combinations of oncogenes and stromal ligands, spanning from fibroblast-induced Clusterin (CLU)+ revival colonic stem cells (revCSCs) to oncogene-driven LRIG1+ hyper-proliferative CSCs (proCSCs). The transition from revCSCs to proCSCs is regulated by decreasing WNT3A and TGF-β-driven YAP signaling and increasing KRASG12D or stromal EGF/Epiregulin-activated MAPK/PI3K flux. We find that APC loss and KRASG12D collaboratively limit access to revCSCs and disrupt stromal-epithelial communication-trapping epithelia in the proCSC fate. These results reveal that oncogenic mutations dominate homeostatic differentiation by obstructing cell-extrinsic regulation of cell-fate plasticity.

MAF amplification licenses ERα through epigenetic remodelling to drive breast cancer metastasis

MAF amplification increases the risk of breast cancer (BCa) metastasis through mechanisms that are still poorly understood yet have important clinical implications. Oestrogen-receptor-positive (ER+) BCa requires oestrogen for both growth and metastasis, albeit by ill-known mechanisms. Here we integrate proteomics, transcriptomics, epigenomics, chromatin accessibility and functional assays from human and syngeneic mouse BCa models to show that MAF directly interacts with oestrogen receptor alpha (ERα), thereby promoting a unique chromatin landscape that favours metastatic spread. We identify metastasis-promoting genes that are de novo licensed following oestrogen exposure in a MAF-dependent manner. The histone demethylase KDM1A is key to the epigenomic remodelling that facilitates the expression of the pro-metastatic MAF/oestrogen-driven gene expression program, and loss of KDM1A activity prevents this metastasis. We have thus determined that the molecular basis underlying MAF/oestrogen-mediated metastasis requires genetic, epigenetic and hormone signals from the systemic environment, which influence the ability of BCa cells to metastasize.

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.

Phenotypic maps for precision medicine: a promising systems biology tool for assessing therapy response and resistance at a personalized level

In this perspective we discuss how tumor heterogeneity and therapy resistance necessitate a focus on more personalized approaches, prompting a shift toward precision medicine. At the heart of the shift towards personalized medicine, omics-driven systems biology becomes a driving force as it leverages high-throughput technologies and novel bioinformatics tools. These enable the creation of systems-based maps, providing a comprehensive view of individual tumor's functional plasticity. We highlight the innovative PHENOSTAMP program, which leverages high-dimensional data to construct a visually intuitive and user-friendly map. This map was created to encapsulate complex transitional states in cancer cells, such as Epithelial-Mesenchymal Transition (EMT) and Mesenchymal-Epithelial Transition (MET), offering a visually intuitive way to understand disease progression and therapeutic responses at single-cell resolution in relation to EMT-related single-cell phenotypes. Most importantly, PHENOSTAMP functions as a reference map, which allows researchers and clinicians to assess one clinical specimen at a time in relation to their phenotypic heterogeneity, setting the foundation on constructing phenotypic maps for personalized medicine. This perspective argues that such dynamic predictive maps could also catalyze the development of personalized cancer treatment. They hold the potential to transform our understanding of cancer biology, providing a foundation for a future where therapy is tailored to each patient's unique molecular and cellular tumor profile. As our knowledge of cancer expands, these maps can be continually refined, ensuring they remain a valuable tool in precision oncology.

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

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

An introduction to representation learning for single-cell data analysis

Single-cell-resolved systems biology methods, including omics- and imaging-based measurement modalities, generate a wealth of high-dimensional data characterizing the heterogeneity of cell populations. Representation learning methods are routinely used to analyze these complex, high-dimensional data by projecting them into lower-dimensional embeddings. This facilitates the interpretation and interrogation of the structures, dynamics, and regulation of cell heterogeneity. Reflecting their central role in analyzing diverse single-cell data types, a myriad of representation learning methods exist, with new approaches continually emerging. Here, we contrast general features of representation learning methods spanning statistical, manifold learning, and neural network approaches. We consider key steps involved in representation learning with single-cell data, including data pre-processing, hyperparameter optimization, downstream analysis, and biological validation. Interdependencies and contingencies linking these steps are also highlighted. This overview is intended to guide researchers in the selection, application, and optimization of representation learning strategies for current and future single-cell research applications.

Multiscale geometric and topological analyses for characterizing and predicting immune responses from single cell data

Single cell genomics has revolutionized our ability to map immune heterogeneity and responses. With the influx of large-scale data sets from diverse modalities, the resolution achieved has supported the long-held notion that immune cells are naturally organized into hierarchical relationships, characterized at multiple levels. Such a multigranular structure corresponds to key geometric and topological features. Given that differences between an effective and ineffective immunological response may not be found at one level, there is vested interest in characterizing and predicting outcomes from such features. In this review, we highlight single cell methods and principles for learning geometric and topological properties of data at multiple scales, discussing their contributions to immunology. Ultimately, multiscale approaches go beyond classical clustering, revealing a more comprehensive picture of cellular heterogeneity.

Chromatin Organization and Transcriptional Programming of Breast Cancer Cell Identity

The advent of sequencing technologies for assessing chromosome conformations has provided a wealth of information on the organization of the 3-dimensional genome and its role in cancer progression. It is now known that changes in chromatin folding and accessibility can promote aberrant activation or repression of transcriptional programs that can drive tumorigenesis and progression in diverse cancers. This includes breast cancer, which comprises several distinct subtypes defined by their unique transcriptomes that dictate treatment response and patient outcomes. Of these, basal-like breast cancer is an aggressive subtype controlled by a pluripotency-enforcing transcriptome. Meanwhile, the more differentiated luminal subtype of breast cancer is driven by an estrogen receptor-dominated transcriptome that underlies its responsiveness to antihormone therapies and conveys improved patient outcomes. Despite the clear differences in molecular signatures, the genesis of each subtype from normal mammary epithelial cells remains unclear. Recent technical advances have revealed key distinctions in chromatin folding and organization between subtypes that could underlie their transcriptomic and, hence, phenotypic differences. These studies also suggest that proteins controlling particular chromatin states may be useful targets for treating aggressive disease. In this review, we explore the current state of understanding of chromatin architecture in breast cancer subtypes and its potential role in defining their phenotypic characteristics.

Inertial effect of cell state velocity on the quiescence-proliferation fate decision in breast cancer

Energy landscapes can provide intuitive depictions of population heterogeneity and dynamics. However, it is unclear whether individual cell behavior, hypothesized to be determined by initial position and noise, is faithfully recapitulated. Using the p21-/Cdk2-dependent quiescence-proliferation decision in breast cancer dormancy as a testbed, we examined single-cell dynamics on the landscape when perturbed by hypoxia, a dormancy-inducing stress. Combining trajectory-based energy landscape generation with single-cell time-lapse microscopy, we found that initial position on a p21/Cdk2 landscape did not fully explain the observed cell-fate heterogeneity under hypoxia. Instead, cells with higher cell state velocities prior to hypoxia, influenced by epigenetic parameters, tended to remain proliferative under hypoxia. Thus, the fate decision on this landscape is significantly influenced by "inertia", a velocity-dependent ability to resist directional changes despite reshaping of the underlying landscape, superseding positional effects. Such inertial effects may markedly influence cell-fate trajectories in tumors and other dynamically changing microenvironments.

Evolving cognition of the JAK-STAT signaling pathway: autoimmune disorders and cancer

The Janus kinase (JAK) signal transducer and activator of transcription (JAK-STAT) pathway is an evolutionarily conserved mechanism of transmembrane signal transduction that enables cells to communicate with the exterior environment. Various cytokines, interferons, growth factors, and other specific molecules activate JAK-STAT signaling to drive a series of physiological and pathological processes, including proliferation, metabolism, immune response, inflammation, and malignancy. Dysregulated JAK-STAT signaling and related genetic mutations are strongly associated with immune activation and cancer progression. Insights into the structures and functions of the JAK-STAT pathway have led to the development and approval of diverse drugs for the clinical treatment of diseases. Currently, drugs have been developed to mainly target the JAK-STAT pathway and are commonly divided into three subtypes: cytokine or receptor antibodies, JAK inhibitors, and STAT inhibitors. And novel agents also continue to be developed and tested in preclinical and clinical studies. The effectiveness and safety of each kind of drug also warrant further scientific trials before put into being clinical applications. Here, we review the current understanding of the fundamental composition and function of the JAK-STAT signaling pathway. We also discuss advancements in the understanding of JAK-STAT-related pathogenic mechanisms; targeted JAK-STAT therapies for various diseases, especially immune disorders, and cancers; newly developed JAK inhibitors; and current challenges and directions in the field.

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.

Epithelial-Mesenchymal Plasticity: Implications for Ferroptosis Vulnerability and Cancer Therapy

Cancers polarized to a mesenchymal or poorly differentiated state can often evade cell death induced by conventional therapies. The epithelial-mesenchymal transition is involved in lipid metabolism and increases polyunsaturated fatty acid levels in cancer cells, contributing to chemo- and radio-resistance. Altered metabolism in cancer enables invasion and metastasis but is prone to lipid peroxidation under oxidative stress. Cancers with mesenchymal rather than epithelial signatures are highly vulnerable to ferroptosis. Therapy-resistant persister cancer cells show a high mesenchymal cell state and dependence on the lipid peroxidase pathway, which can respond more sensitively to ferroptosis inducers. Cancer cells may survive under specific metabolic and oxidative stress conditions, and targeting this unique defense system can selectively kill only cancer cells. Therefore, this article summarizes the core regulatory mechanisms of ferroptosis in cancer, the relationship between ferroptosis and epithelial-mesenchymal plasticity, and the implications of epithelial-mesenchymal transition for ferroptosis-based cancer therapy.

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

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

CLINICAL HETEROGENEITY IN THE AGE OF BIG DATA, ADVANCED ANALYTICS, AND COMPLEXITY THEORY

Clinical heterogeneity remains a challenge in the practice of medicine and is an underlying motivation for much of biomedical research. Unfortunately, despite an abundance of technologies capable of producing millions of discrete data elements with information about a patient's health status or disease prognosis, our ability to translate those data into meaningful improvements in understanding of clinical heterogeneity is limited. To address this gap, we have applied newer approaches to manifold learning and developed additional and complementary techniques to interrogate and interpret complex, high dimensional omics data. The central premise is that there exist manifolds embedded in high dimensional data that represent fundamental biologic processes that may help address the challenges of clinical heterogeneity. Preliminary evidence from several real-world data sets suggests that these techniques can identify coherent and reproducible manifolds embedded in higher dimensional omics data. Work is currently ongoing to determine the clinical informativeness of these novel data structures.

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.