Abstract
Phenotypic plasticity of cancer cells can lead to complex cell state dynamics during tumor progression and acquired resistance. Highly plastic stem-like states may be inherently drug-resistant. Moreover, cell state dynamics in response to therapy allow a tumor to evade treatment. In both scenarios, quantifying plasticity is essential for identifying high-plasticity states or elucidating transition paths between states. Currently, methods to quantify plasticity tend to focus on 1) quantification of quasi-potential based on the underlying gene regulatory network dynamics of the system; or 2) inference of cell potency based on trajectory inference or lineage tracing in single-cell dynamics. Here, we explore both of these approaches and associated computational tools. We then discuss implications of each approach to plasticity metrics, and relevance to cancer treatment strategies.
1 Introduction
1.1 Overview
In the field of Network Physiology, cancer systems biology occupies an intriguing position. On the one hand, widespread research efforts advance data production from genes to patients and, in parallel, improving analytical methods for inferring molecular networks from these large datasets are providing insights that both leverage and go beyond reductionism-based knowledge. On the other hand, the built-in plasticity of heterogeneous cell states in tumors and the consequent lack of ground truths typically elucidated in physiological systems create profound uncertainty about structure and dynamics of cancer networks, whether inferred from top-down or bottom-up approaches.
In this review, we place studies on cancer cell plasticity and its underlying network dynamics in the context of broader studies on the regulation of cell plasticity in physiological self-organizing systems, as it occurs in, for example, brain or embryo development (Ivanov, 2021). In a nutshell, cancer has been understood as a disease in which regulation of the cell cycle is lost and cell proliferation has become a runaway process. This is an actionable perspective that has led to many advances in cancer treatment. However, in a larger sense, cancer is a disease of lost cell identity: tumors can be shrunk, slowed down, or almost eradicated, but in the vast majority of cases they relapse in a treatment-resistant or -tolerant state. Our current understanding of relapse is rooted on studies that unambiguously determined the heterogeneous nature of cancer cell states in a tumor. More recently, transitions among these states have been convincingly demonstrated. In this Introductory section, the evidence for tumor heterogeneity and cancer cell plasticity is first summarized, and the case is made for the key role of quantitative metrics for heterogeneity and plasticity.
In the rest of Introduction, theoretical frameworks for plasticity are recalled. The current noisy landscape of information formed by torrents of publications and a dataset tsunami can be overwhelming. We find it essential to grasp for theory as an anchor in reality, and a means for producing knowledge platforms that can be hardened and continuously improved upon (or falsified).
In later sections, current attempts at unveiling the mechanistic basis for plasticity are reviewed, with emphasis on the role of the dynamics of Gene Regulatory Networks (GRNs), and the dynamics of single-cell state transitions. This focus was motivated both by our direct experience in these areas, and by a broadening group of investigators that are collectively producing remarkable advances.
1.2 Cell heterogeneity and plasticity in cancer
Heterogeneity within tumors has been shown to be critical for acquired resistance to therapy in many cancer types (Altschuler and Wu, 2010; Calbo et al., 2011; Marusyk et al., 2012; Huang, 2013; Frick et al., 2015; Pisco and Huang, 2015; Jia et al., 2017; Lim et al., 2017; Su et al., 2019; Nath et al., 2021; Yabo et al., 2021). Several layers of heterogeneity exist and play a large role in cancer systems (Elowitz et al., 2002; Brock et al., 2009; Feinberg and Irizarry, 2010; Gupta et al., 2011; Pisco and Huang, 2015; Caiado et al., 2016; Kumar et al., 2019; Hayford et al., 2021). Genetic heterogeneity results from selection of mutants, each of which may respond differently to treatment. Non-genetic (or epigenetic) heterogeneity is dependent on epigenetic regulation of phenotype and can be attributed to several sources, including variability in chromatin accessibility, DNA methylation, and DNA-binding proteins that regulate transcription levels of genes. Finally, stochasticity arises from intrinsic sources, such as the probabilistic nature of biochemical reactions within a cell, or extrinsic sources, such as local fluctuations in chemical concentrations in the microenvironment (Swain et al., 2002). While transient, this variability can probabilistically drive transitions between phenotypes (Feinberg and Irizarry, 2010; Gupta et al., 2011; Liao et al., 2012; Hayford et al., 2021).
Together, these layers of heterogeneity—genetic, epigenetic, and stochastic—define the variability in phenotype. There is a critical need to quantify these levels of heterogeneity in cancer systems, as distinct phenotypes will presumably respond differentially to treatment, and changes in heterogeneity can lead to acquired resistance (Pisco and Huang, 2015; Brady et al., 2017; Jia et al., 2017; Mu et al., 2017; Zou et al., 2017; Jolly et al., 2018; Risom et al., 2018; Arozarena and Wellbrock, 2019; Su et al., 2019; Nath et al., 2021). Such dynamics of phenotype, or phenotypic plasticity, can lead to differential treatment response and/or resistance in several ways, including: 1) the existence of highly plastic, drug-resistant states; and/or, 2) cell state dynamics that evade treatment (Figure 1).
FIGURE 1
First, a particular phenotype may be intrinsically less susceptible to treatment, so transitions to this “drug-tolerant persister” phenotype can promote tumor survival (Sharma et al., 2010; Liau et al., 2017; Paudel et al., 2018; Risom et al., 2018; Jia et al., 2020; Cabanos and Hata, 2021; Oren et al., 2021). Often, such a phenotype has stem cell-like properties, suggesting the plastic potential of stem cell-like phenotypes is intrinsically tied to treatment response (Gupta et al., 2011; Chisholm et al., 2015; Liau et al., 2017; Wainwright and Scaffidi, 2017; Smith et al., 2018; Lytle et al., 2019; Neftel et al., 2019; Chan et al., 2021; Yabo et al., 2021). Lineage tracing analyses, as described in Section 3.3., that investigate the underlying mechanisms of the persister state and the transition paths towards it can point to strategies for reprogramming such states to sensitivity.
Second, cell state dynamics between various phenotypes can promote tumor survival by adaptation to treatment (Zhou and Li, 2016; Neftel et al., 2019; Wouters et al., 2020; Gay et al., 2021; Nath et al., 2021; Groves et al., 2022; Sutherland et al., 2022). In these cases, reprogramming cells towards a drug-sensitive state is infeasible, because the high degree of cell state transitions can allow for any cell state to become insensitive again. In this case, it would appear that plasticity itself should be the target.
In both scenarios, it is necessary to quantify the phenotypic plasticity of cancer cell states and the mechanisms underlying cell state dynamics, towards the goal of identifying therapeutic strategies that diminish the plastic capabilities of the tumor as a whole (Huang and Kauffman, 2013). Waddington’s landscape is a useful metaphor for understanding how cancer cells shift between phenotypes and can be quantified through the underlying gene regulatory network dynamics or via statistical mechanical modeling of cell state dynamics, such as those seen in single cell transcriptomics datasets.
1.3 Waddington’s landscape in cancer
In 1957, C.H. Waddington proposed the concept of an epigenetic landscape for understanding the regulation of phenotype in the context of biological differentiation (Waddington, 1957). In this analogy, cells roll downhill through canalized channels or “chreods” representing differentiation pathways. Cells at the top of the landscape are pluripotent stem cells, and as they travel down the landscape, they gradually become more committed to a particular cell fate. Thus, the epigenetic landscape could be thought of as a tool to uncover how epigenetic regulation in a cell (e.g., through chromatin accessibility or DNA-binding of transcripton factors) controls the cell’s phenotype. Waddington initially characterized this regulation as a complex system of interactions that he illustrated as strings pulling on and shaping the landscape from below (Waddington, 1957).
In normal development, cells are generally isogenic. In cancer, however, where the mutation rate is higher and multiple subclones may exist within a single tumor, genetic heterogeneity can be represented by a “fitness landscape” (Figure 2, bottom). In this landscape, mutants with higher fitness will be selected for via Darwinian evolution. For each location in the fitness landscape (each genome), an entire Waddington landscape of phenotypes exists (Figure 2, top). Similar to Waddington’s original conception, cells in the epigenetic landscape “fall downhill” towards the states with the lowest “potential.” These phenotypic transitions depend on the instability of each cell state, and a cell’s ability to transition can be defined by its plasticity.
FIGURE 2
1.4 System attractors, instability, and plasticity
The notion of plasticity goes hand in hand with the dynamical systems theoretical concept of instability. In dynamical systems, stability of a state requires more than stationarity; a stable state is one that is resilient to perturbations such that, after external influences such as changing microenvironmental conditions, the system returns to its original state (Huang, 2013). This idea is represented in the potential landscape, in which cells roll downhill toward local minima, as shown in Figure 3 (top). While there may be steady states throughout the landscape, such as the top of a flat hill or the bottom of a valley, a small push to a cell on top of a hill will cause it to roll down to a local minimum, far from its original starting state. On the contrary, a cell in a local minimum is resilient to small perturbations: it is in a stable “attractor” state of the landscape (Huang, 2009). The high-dimensional region around the attractor where a cell will roll towards the attractor is called the basin of attraction (Figure 3, bottom). Cell states with larger basins of attraction can withstand larger perturbations to their cell state, thereby demonstrating resilience of the system.
FIGURE 3
In dynamical systems theory, plasticity is a weaker kind of stability, in which a perturbed system neither returns to its original state nor escapes from it, but instead tracks the environmental change (Huang, 2013). However, in biology, plasticity and instability are often thought of as interchangeable: a more plastic cell state responds to an external perturbation by changing its state to a larger degree. In this view, cells with higher potential on the landscape are considered more plastic, as they are more likely to move through the landscape towards a lower-potential attractor. Quantifying both the quasi-potential of the landscape and potential trajectories through a landscape can elucidate the role of phenotype plasticity in cancer, such as plasticity in response to treatment. For example, by quantifying quasi-potential, one can identify highly plastic and/or stem-like cells that may be responsible for tumor propagation (Sharma et al., 2010; Grosse-Wilde et al., 2015; Chan et al., 2021; Gay et al., 2021). By identifying trajectories through the landscape of a tumor cell population, one can characterize the paths cells take epigenetically during tumor development, tumor metastasis, persistence and acquired resistance in response to treatment.
1.5 Quantifying plasticity as quasi-potential of Waddington’s landscape
While Waddington intended this picture purely as a metaphor, it has now been quantified in various ways, borrowing ideas from physics and dynamical systems theory to describe the underlying regulation of these processes (Wang et al., 2008; 2011; Zhou et al., 2012). The height of the landscape describes instability of each phenotype as a “quasi-potential,” a correlate of gravitational potential in a physical landscape (Figure 2, top). Quantification of this quasi-potential is informative for processes in which plasticity and instability plays a central role, including cancer systems (Huang et al., 2009; Huang and Kauffman, 2013; Hanahan, 2022). By modeling potential in an epigenetic landscape of phenotypically heterogeneous populations, one can better determine ways to control the permissivity of phenotype and prevent reprogramming of cell identity from a sensitive phenotype to an insensitive one, as often seen in acquired resistance.
Borrowing from physics, movement of cells in the landscape (i.e., changes in xi (x1,x2, … ,xN) over time, where x is the location of a cell in phenotype space) may be due to some “force” F(x), similar to the effect of gravity on movement through a physical landscape. A potential, U(x), can be defined such that the change in phenotype is equal to the gradient of this potential:
Cells will therefore “roll down” the gradient towards states with lower potential. It is worth noting that most high-dimensional, non-equilibrium biological systems are not simple gradient systems, and therefore the vector field is sometimes decomposed into two components: the gradient of some quasi-potential, and a remainder term (Wang et al., 2008). Still, the gradient term has been successfully used to understand pathways of transition through epigenetic landscapes, describing everything from differentiation to cell fate reprogramming (Wang et al., 2006; 2010; Zhou and Huang, 2010; Zhou et al., 2012; Wu and Wang, 2013a; 2013b; Li and Wang, 2014a; 2014b; Wang, 2015; Zhou and Li, 2016; Luo et al., 2017; Yan et al., 2019). Furthermore, the high dimensionality of complex biological systems can pose a problem for characterizing the structure of an interpretable, lower-dimensional epigenetic landscape. Recent work addressed this problem using a dimension reduction approach of the landscape (DRL), which projects high-dimensional landscapes into a lower-dimensional coordinate system based on variance in an associated probability density function (Kang and Li, 2021). This method was applied to cancer systems in the context of epithelial to mesenchymal transitions and metabolism (Kang and Li, 2021).
Several systems biology approaches have been developed to determine the driving force F(x) that shapes the epigenetic landscape and defines phenotypic plasticity (Huang, 2012; Zhou et al., 2012; Devaraj and Bose, 2020). Classical dynamical systems modeling of underlying gene regulatory networks is a bottom-up approach that can explain how phenotypic transitions are dependent on regulation of gene expression by transcription factors (TFs) (Bhattacharya et al., 2011; Wang et al., 2011; Joo et al., 2018). Alternatively, phenomenological top-down approaches based on analysis of large ‘omics’ datasets can approximate the potential landscape. For example, single-cell sequencing of transcriptomes samples the density of cells in the landscape and trajectory inference methods uncover transition paths between attractors, i.e., cell states (Saelens et al., 2019). These two orthogonal approaches are detailed in the following two sections.
2 Modeling plasticity in epigenetic landscapes via gene regulatory networks
2.1 Gene regulatory network structure and dynamics
To understand the driving force F(x) that defines the landscape quasi-potential, first it is important to understand gene regulatory networks (GRNs). A GRN is established by the fact that certain genes encode TFs which are capable of binding to DNA and regulate transcription of other genes into RNA. Because TFs can also control the transcription of other TFs (and sometimes themselves), a network of TFs and the genes they regulate can be constructed (Figure 4). The structure of the GRN for a particular cell is hardcore in the genome of a cell, as shown in Figure 5 (left), since each interaction in the network is a molecular interaction between a DNA-binding protein and the cis-regulatory loci (such as promoter and enhancer regions) for a particular gene (Huang, 2013). On the other hand, the dynamics of the network are described by the collective changes in gene expression over time. The dynamics of a GRN allow for various stable states dependent on the expression of genes in the network (Figure 5, right). Therefore, the state of the GRN, given by the expression levels of the genes within it, maps to a single location on the epigenetic landscape—the phenotypic state.
FIGURE 4
FIGURE 5
Quantifying the dynamics of TF binding can be calculated by adapting Hill kinetics to describe the rate at which a target gene is transcribed when regulated by TFs (Hill, 1913). The Hill equation is a sigmoidal function that describes activation (or repression) of a gene as dependent on the concentration of a regulator until it reaches saturation. This is a relatively realistic description of many gene control functions and can be derived directly from the binding of the TF to the promoter site. The dynamics, or the change over time, of each TF in the network can therefore be represented as a function of all “upstream” parent nodes in the network that influence its transcription. The system of such differential equations, where each TF in the network has a corresponding equation for its rate of change, defines the complete dynamics of the system. Based on this system of equations, GRN dynamics are equivalent to the driving force that pushes cells down the gradient of the potential in the landscape.
Several researchers have directly solved such systems of ODEs to quantify the plasticity of various systems (Wang et al., 2010; Wang et al., 2011; Zhou et al., 2012; Li and Wang, 2014b; Zhou and Li, 2016; Devaraj and Bose, 2020). Transition paths between stable states can then be calculated, such as by using a path-integral approach (Li et al., 2016; Lang et al., 2021). However, for high-dimensional systems, this system of equations often becomes intractable. Instead, several computation methods have been developed to approximate network interactions.
2.2 GRN simulations can model phenotypic transitions
In 1969, Stuart Kauffman introduced the idea of Boolean network models for biological systems (Kauffman S., 1969; Kauffman, S. A. 1969). Kauffman posited that, “while finely-graded intermediate levels of gene activity could occur,” genes tended to be very active or very inactive (Kauffman, 1971), consistent with switch-like Hill kinetics with a high Hill coefficient. Therefore, it is often useful to idealize the control of gene expression as a binary switch. Boolean logic determines the activity level of each gene given the binary states of its upstream regulating TFs by approximating the Hill equation, turning the smooth, monotonic function into a step function with activation (or repression) threshold of S (Kauffman, 1971; Glass and Kauffman, 1973; Thieffry and Thomas, 1998). Since Kauffman’s original idea, several studies have shown the utility of conceptualizing gene regulation as a set of binary genes coupled together through Boolean functions (Albert et al., 2008; Pomerance et al., 2009; Saadatpour and Albert, 2013; Masoudi-Nejad et al., 2015; Steinway et al., 2015; Zhou et al., 2016; Correia et al., 2018; Joo et al., 2018; Yachie-Kinoshita et al., 2018; Wooten et al., 2019).
While a Boolean approximation for transcriptional regulation is realistic for many biological systems, some genes are regulated by multiple TFs in a manner that does not use Boolean logic. For example, Kalir and Alon (2004) showed that gene regulation in an E. coli network of flagella biosynthesis follows a summation function (SUM), rather than Boolean logic gates (AND, OR, and NOT). Several studies have shown other functions, including complex functions with many inputs, are also possible (Yuh et al., 1998; Beer and Tavazoie, 2004; Istrail and Davidson, 2005). To model such complex systems, other types of networks must be used. One such approach is to adapt Boolean networks using probabilistic rules (Dorigo, 1994; Shmulevich et al., 2002; Li et al., 2007; Trairatphisan et al., 2013; Tercan et al., 2022). Probabilistic Boolean Networks (PBNs) find a set of functions for each node in the network, each with an associated probability of predicting the target node.
In order to understand systems with non-Boolean gene regulatory functions, other probabilistic methods of network inference, known as probabilistic graphical models (PGMs), may be used (Wang et al., 2005; Zou and Conzen, 2005; Li et al., 2011; Baba et al., 2014; Sanchez-Castillo et al., 2017). These models have multiple advantages over Boolean approaches. For example, they can infer non-linear relationships between TFs, so that the rule of interaction is not required a priori to have a particular form such as a Boolean function. One such PGM, known as a Bayesian network, considers a GRN to be a network (or graph) where each directed edge represents the probabilistic dependence among genes. PGMs are more phenomenological than Hill kinetics or Boolean network modeling, but they can mine information from transcriptomic data—for example, RNA-seq profiles for the TFs in the network—without assumptions, such as binariness, about the relationships between TFs (Li et al., 2007; Chai et al., 2014).
Several methods focus on predicting network structure alone (Margolin et al., 2004; Langfelder and Horvath, 2008; Huynh-Thu et al., 2010; Aibar et al., 2017; Chan et al., 2017; Moerman et al., 2018). Some of these approaches utilize similarity metrics on transcriptomic data, e.g., to identify co-expressed gene modules or find relationships between genes with high mutual information, such as WGCNA, GENIE3, ARACNE, or PIDC (Margolin et al., 2004; Langfelder and Horvath, 2008; Huynh-Thu et al., 2010; Chan et al., 2017). While these approaches have been successfully applied to several systems, including cancer, they can often find spurious relationships that do not correspond to physical interactions (cis-regulatory motifs, such as transcription factors binding to the promoter of a target gene). More recent methods can also incorporate this binding information to predict regulatory relationships. For example, SCENIC builds a network structure based on gene co-expression modules and transcription factor binding motif information from the RcisTarget database (Aibar et al., 2017). These tools identify network interactions that coordinate changes in cell identity, but do not predict single cell dynamics.
Many computational algorithms have also been developed to infer both network structure and single-cell dynamics based on Boolean, Bayesian, or other regulatory rules (Chan et al., 2017; Khan et al., 2017; Sanchez-Castillo et al., 2017; Chen and March 2018; Castro et al., 2019; Dunn et al., 2019; Wooten et al., 2019; Aalto et al., 2020; de Sande et al., 2020; Pratapa et al., 2020; Ramirez et al., 2020; Su et al., 2022; Hérault et al., 2023; Kamimoto et al., 2023). For example, BooleaBayes uses probabilistic Boolean rules to predict master regulators of heterogeneous phenotypes that, when perturbed, could destabilize particular phenotypes and therefore change the phenotypic composition of a tumor (Wooten et al., 2019; Olsen et al., 2021; Groves et al., 2022). SCODE models a GRN via ODEs using the gene expression matrix and associated pseudotime from a single-cell dataset and uses this GRN to reconstruct the expression dynamics (Matsumoto et al., 2017). CellOracle uses scRNA-seq and scATAC-seq to generate GRNs and simulate changes in gene expression following experimental perturbations (Kamimoto et al., 2023).
Regardless of the limitations or assumptions of network inference algorithms, these methods require biological data to fully characterize a system. Transcriptomics data are often used, sometimes in combination with other types of epigenomic or proteomic information (Margolin et al., 2004; Langfelder and Horvath, 2008; Liu et al., 2016; Duren et al., 2017; Ramirez et al., 2017; Wooten et al., 2019). Today, single-cell RNA-sequencing (scRNA-seq) is commonly used to obtain a more granular picture of transcriptional regulation and stable phenotypes in a system than bulk sequencing data can provide. Top-down, phenomenological approaches for modeling the epigenetic landscape can also utilize scRNA-seq data directly to find empirical patterns of expression. Because intratumoral heterogeneity and plasticity are relevant to acquire resistance in cancer, it is important to determine how cells change in phenotype in various contexts. These top-down approaches work towards the long-term goal of personalizing treatment by providing a framework for understanding plasticity in an individual patient’s tumor.
3 Modeling plasticity in epigenetic landscapes via single-cell dynamics
Single-cell sequencing methods have paved the way for data-driven approaches to quantifying plasticity. While classical dynamical systems modeling—i.e., modeling a GRN that determines a quasi-potential landscape—has the advantage of being predictive, it can be difficult or impossible to model the complete dynamics of a high-dimensional system. Alternatively, it is possible to use a data-driven, bottom-up approach by modeling single-cell dynamics as a Markovian process, which can identify transition paths heuristically from scRNA-seq data.
Borrowing once again from physics, a drift-diffusion equation can model the change in cell density for a given region of gene expression space (or, analogically, the phenotypic landscape):
where c is cell density of a given region of gene expression space, R describes the rate of accumulation and loss due to cell proliferation, death, and movement through the region, and v is the net average velocity (Weinreb et al., 2018). With additional assumptions, we can model the velocity as related to the deterministic average velocity field (due to the epigenetic landscape, for example,) and a stochastic component related to diffusion. This velocity field may be calculated heuristically from pseudo-temporal information using trajectory inference methods and can predict cell-state transitions in the epigenetic landscape (Qiu et al., 2022). Furthermore, drift-diffusion modeling of cell dynamics along a high-dimensional manifold in gene expression space can be used to infer dynamics through a Markov chain, with defined transition probabilities between cell states (Weinreb et al., 2018).
3.1 Trajectory inference and pseudotime as a measure of plasticity
Trajectory inference algorithms also aim to understand changes in cell density by ordering cells along a trajectory based on transcriptomic similarity, empirically determining transition paths in the system (Trapnell et al., 2014; Guo et al., 2016; Haghverdi et al., 2016; Welch et al., 2016; Qiu et al., 2017; Herring et al., 2018; Saelens et al., 2018; Setty et al., 2019; Wolf et al., 2019; Stassen et al., 2021). These trajectory inference algorithms tend to search for an underlying manifold of the data to delineate graph-based trajectories. By interrogating the structure of the single cell data in gene expression space, multifurcations, trees, and other graph structures can be identified. While these methods are unbiased and often unsupervised, they tend to require identification of a “root cell” to determine the directionality of transitions, as multiple trajectories could be explained by the same graph structure. Such a root cell, or source, can be thought of as having a high degree of plasticity, as defined by the quasi-potential of the underlying landscape. Therefore, these methods require a priori knowledge of the high-plasticity states of a system but are useful for identifying transition paths from these states.
Some methods utilize time-series data to determine directionality by optimal transport-based algorithms (Kimmel et al., 2019; Schiebinger et al., 2019; Marjanovic et al., 2020). Because scRNA-seq is a destructive method, the same single cell cannot be monitored and sequenced over time. Optimal transport-based methods overcome this experimental constraint by inferring “temporal couplings” across timepoints to determine the most likely phenotypic “descendants” of each cell at later timepoints. Ultimately, lineage tracing provides a benchmark for interrogating trajectories, as cell lineages across timepoints are identified via “barcodes,” thereby linking cell state in early timepoints to cell fate in later timepoints (Griffiths et al., 2018; Wagner and Klein, 2020; Wang et al., 2021).
3.2 RNA velocity-based measures of plasticity
In 2018, a novel approach to trajectory inference was developed based on RNA splicing dynamics (La Manno et al., 2018). By fitting an ordinary differential equation (ODE) model of RNA transcription, splicing, and degradation, La Manno et al. discovered that it was possible to infer short-term dynamics on a cell-by-cell basis (Figure 6). RNA velocity infers a steady-state ratio of unspliced to spliced counts of RNA on a gene-by-gene basis to fit the ODE model parameters, such as the degradation rate of the mRNA. As shown in Figure 6, an increase in RNA transcription from a particular gene results in a slow increase of unspliced RNA, followed by a delayed increase in spliced RNA. Therefore, by comparing the unspliced and spliced counts of a gene in each cell in this model, it is possible to determine the future state of each cell.
FIGURE 6
The timeframe for dynamic predictions is on the order of a few hours, similar to the average splicing rate. However, RNA velocity can be extrapolated to longer timeframes by considering the relationship between a cell’s velocity vector—i.e., the directionality and magnitude of its inferred change in gene expression—and the location of neighboring cells (Figure 6, right). These extrapolated dynamics can be used to make predictions about the future state of cells near the beginning of the trajectory. Because this method does not rely on multiple sampled timepoints or prior knowledge about the “root” cell of a trajectory, it is optimal for understanding the dynamics of systems for which a temporal series of samples is not possible, such as tumor dynamics from single biopsies. Together, these analysis methods can uncover an empirical epigenetic landscape by defining stable phenotypes and transition paths in scRNA-seq data sampled from various cancer systems, including human biopsies, to complement or replace quantification of GRN dynamics.
Since RNA velocity was first introduced, several methods have utilized the approach to quantify plasticity (Bergen et al., 2020; Gorin et al., 2020; Chen et al., 2022; Groves et al., 2022; Lange et al., 2022; Qiu et al., 2022). As described in the original paper, trajectory inference from RNA velocity translates the velocity vectors into a transition probability matrix, which can be used in a Markov Chain model of cell state dynamics (La Manno et al., 2018; Bergen et al., 2020). This approach assumes that movement of cells through a phenotypic landscape is a Markovian process in which cell fate depends only on the current state of the cell. With this assumption, we (Groves et al., 2022) and others (Weinreb et al., 2018; Bergen et al., 2020) proposed to quantify plasticity as a cell’s potential to move towards one or more attractors of the system, i.e., a cell’s ability to traverse a phenotypic landscape. Our metric, termed Cell Transport Potential (CTrP), quantifies the average distance traveled for each cell through a Markovian state transition graph, accounting for multiple possible cell fates (or absorbing states) (Groves et al., 2022). The gain with CTrP is an intuitive connection between cell landscape dynamics and transcriptomics: the higher the CTrP of a cell in a landscape, the larger the mobility in gene expression space expected for that cell. With CTrP, we identified highly plastic cell states across several small cell lung cancer (SCLC) human and mouse experimental models (Gay et al., 2021; Groves et al., 2022). For instance, in a circulating tumor cell-derived xenograft (CDX) model, we were able to determine that resistant tumors post-treatment originated from a small, high-CTrP cell state that arose after chemotherapy (Gay et al., 2021).
Other approaches have quantified plasticity as multipotency by defining possible cell fates for each cell state using RNA velocity-based transition probabilities. For example, CellRank builds on RNA velocity and trajectory inference models (such as pseudotime) to predict fate potentials given the stochastic nature of fate decisions (Lange et al., 2022). CellRank has been used to predict fate probabilities and reprogramming outcomes in several developmental systems (Hersbach et al., 2022; Lange et al., 2022; Van Bruggen et al., 2022; Bono et al., 2023; Matsushita et al., 2023), whereas applications to cancer systems have mainly focused on trajectory inference for the immune compartment rather than cancer cells (Xue et al., 2022; Friedrich et al., 2023; Jainarayanan et al., 2023). DeepVelo uses neural networks to learn transcriptomic dynamics, building a model that predicts trajectories, driver genes, and the effect of in silico perturbation on fate decisions (Chen et al., 2022); however, this approach has not yet been applied to cancer systems.
These approaches have been used to predict perturbations that can affect fate decisions. In cancer, these methods could identify treatment options for perturbing cells away from drug-resistant cell types (Wooten and Quaranta, 2017) or cancer attractors as a whole (Huang and Kauffman, 2013; Li et al., 2016).
3.3 Lineage tracing to understand cell state transitions
While trajectory inference of single cell sequencing data can provide high granularity for understanding phenotypic heterogeneity in cancer, such approaches to understand plasticity of cancer cells over time must account for the destructive nature of sequencing. Alternatively, lineage tracing methods have long been used to understand the temporal dynamics of cell state in cancer populations, particularly during tumor initiation and in response to treatment (Chakrabarti et al., 2018; Wagner et al., 2018; Weinreb et al., 2020; Wang et al., 2021; Singh and Saint-Antoine, 2023). Single cell time-lapse microscopy has been used to correlate cell state and fate, suggesting the existence of phenotype switching (Bhola and Simon, 2009; Spencer et al., 2009; Bertaux et al., 2014; Chakrabarti et al., 2018). For example, researchers used time-lapse microscopy to understand variability in the onset of apoptosis, finding that protein state gives rise to transient heritability between mother and daughter cells (Bhola and Simon, 2009; Spencer et al., 2009). This can be modeled mechanistically by considering stochastic fluctuations in protein levels (Bertaux et al., 2014). Together, these results connect cell state (assumed to be identical in sister cells) with cell fate (divergence in apoptotic response), and pave the way for understanding how subpopulations of a single tumor can have such different fates (drug sensitivity versus tolerance) in response to a single treatment.
More recently, several groups have used a modified Luria Delbrück fluctuation analysis to determine whether resistance to therapy is heritable or a result of transient reprogramming of phenotype (Shaffer et al., 2017; 2020; Russo et al., 2022; Singh and Saint-Antoine, 2023). For example, Shaffer et al. (2017) tested whether resistance to vemurafenib in BRAF-mutated melanoma was genetically heritable or transient. If the drug resistance was transient and due to epigenetic reprogramming, a Luria Delbrück fluctuation analysis would show a tighter distribution of resistant cell colony sizes, because all cells would be equally likely to form a resistant colony. The hypothesis of a transient pre-resistant state that could epigenetically reprogram to a stably resistant state (i.e., persister state) under drug was indeed supported by the results. This state was further characterized by a distinct transcriptional profile (including high expression of EGFR) and activation of transcription factors (JUN, AP-1, and TEAD). Shaffer et al. (2020) then expanded this work into a broadly applicable method, MemorySeq, that combines Luria-Delbrück fluctuation analysis and population-based RNA sequencing. Similarly, Russo et al. (2022) used fluctuation analysis to investigate drug-induced plasticity of colorectal cancer cells. In this cancer system, cell population dynamics were quantified with a mathematical model of transitions to a persister state, which was consistent with a drug-induced, rather than preexisting, persister state.
Together, these experiments and analyses have shown that drug resistance in cancer can arise from epigenetic reprogramming of transient, pre-resistant states, and that high degrees of transcriptional heterogeneity allow for rare cell populations to become stably resistant through plasticity. Importantly, these transitions to a stably resistant state are drug-induced rather than preexisting, solidifying the connection between treatment and plasticity of cancer cells.
4 Conclusion
The success of cancer therapies is often limited by mechanisms of cellular persistence and acquired resistance. Non-genetic plasticity has emerged as a major cause of treatment insensitivity or acquired resistance in several cancer types (Marjanovic et al., 2013; 2020; Pisco et al., 2013; Mu et al., 2017; Su et al., 2017; Zou et al., 2017; Qin et al., 2020; Quintanal-Villalonga et al., 2020; Chan et al., 2021; Hanahan, 2022). Targeting plasticity directly has been suggested as a possible treatment option for several cancers, including melanoma, breast cancer, and prostate cancer (Sáez-Ayala et al., 2013; Kemper et al., 2014; Ahmed and Haass, 2018; Risom et al., 2018; Arozarena and Wellbrock, 2019; Chapman et al., 2019; Boumahdi and de Sauvage, 2020; Yabo et al., 2021).
A few different methods for targeting plasticity can be envisioned. First, cell plasticity could be used advantageously to reprogram cells towards more drug-sensitive states (Yuan et al., 2019). For example, master TFs, identified through GRN analyses, could be controlled to direct phenotype switching to attractors that better respond to treatment, as shown in melanoma (Sáez-Ayala et al., 2013).
Second, preventing phenotype switching may be more desirable (Boumahdi and de Sauvage, 2020). Phenotypic plasticity is intrinsic to the epigenetic landscape: By shaping the landscape, GRN dynamics form transition paths and unused attractors, and cells may enter transition paths between stable attractors due to extrinsic perturbations or intrinsic stochasticity (Huang, 2013). The barrier to exit attractors may be lower in cancer than normal cells, with “de-canalized,” shallow valleys and attractor basins enabling cancer cells to stochastically sample the landscape and find new attractors that evade treatment (Jia et al., 2017). Targeting the mechanisms that allow for this stochastic search of drug-tolerant states in the landscape may lower plasticity and acquired resistance to therapy. For example, chromatin remodeling may be a key mechanism by which cells reprogram to other fates, as open chromatin has been shown to correlate with plasticity (Meshorer and Misteli, 2006; Giadrossi et al., 2007; Gaspar-Maia et al., 2011; Burdziak et al., 2023). In fact, a recent study on pancreatic ductal adenocarcinoma used this connection between plastic cells and accessible chromatin landscape to quantify plasticity as the entropy in prediction of transcriptomic fate based on chromatin accessibility (epigenomic state) (Burdziak et al., 2023). In cancers where plastic states with open chromatin landscapes exist, promoting repressive chromatin organization may be able to keep cells from transitioning between phenotypes during tumor progression or treatment evasion.
Modeling this plasticity through GRNs or single-cell dynamics can lead to new approaches to therapy. Development of strategies that target plasticity and systematically reprogram cell identity may ultimately enable to overcome persistence and acquired resistance in cancer. These goals should not be elusive, if they are rooted in our current understanding of mechanisms of GRN regulation. In a sense, it could be useful to start viewing cancer cells as driven by dysregulated GRNs, rather than by some mysterious “malignant” property (i.e., by a misguided absolute priority for self-preservation, as an invading virus or bacteria would do). Such perspective may engender a longer but perhaps more rewarding path to achieving a solution to this devastating disease.
Single-cell state transitions in response to perturbations are also broadly observed in other physiological systems. In fact, perfect adaptation is well-known in unicellular organisms, and adaptability is essential for tissue homeostasis. Thus, reproducible dynamics in physiological platforms, either spontaneous or under perturbation, can be used to place boundaries on cancer adaptive networks. Vice versa, cancer network studies can provide insights into the potential of physiological networks, e.g., in the context of evolution. We submit that bridging these areas of research will eventually broaden perspectives on network physiology.
Statements
Author contributions
SG contributed to investigation, writing, visualization, and funding acquisition. VQ contributed to writing (review and editing), supervision, project administration, and funding acquisition. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by funds from National Institutes of Health U54CA217450 (VQ and SG), 5T32CA009582-34 (SG) and from National Science Foundation DGE-1445197 (SG).
Conflict of interest
VQ is an Academic co-Founder and equity holder for Parthenon Therapeutics, Inc. and Duet BioSystems, Inc.
The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AaltoA.ViitasaariL.IlmonenP.MombaertsL.GonçalvesJ. (2020). Gene regulatory network inference from sparsely sampled noisy data. Nat. Commun.11, 3493. 10.1038/s41467-020-17217-1
2
AhmedF.HaassN. K. (2018). Microenvironment-driven dynamic heterogeneity and phenotypic plasticity as a mechanism of melanoma therapy resistance. Front. Oncol.8, 173. 10.3389/fonc.2018.00173
3
AibarS.González-BlasC. B.MoermanT.Huynh-ThuV. A.ImrichovaH.HulselmansG.et al (2017). Scenic: single-cell regulatory network inference and clustering. Nat. Methods14, 1083–1086. 10.1038/nmeth.4463
4
AlbertI.ThakarJ.LiS.ZhangR.AlbertR. (2008). Boolean network simulations for life scientists. Source Code Biol. Med.3, 16. 10.1186/1751-0473-3-16
5
AltschulerS. J.WuL. F. (2010). Cellular heterogeneity: do differences make a difference?10.1016/j.cell.2010.04.033
6
ArozarenaI.WellbrockC. (2019). Phenotype plasticity as enabler of melanoma progression and therapy resistance. Nat. Rev. Cancer19, 377–391. 10.1038/s41568-019-0154-4
7
BabaN.MohamadM. S.SallehA. H. M.HijaziM. H. A.ChaiL. E.ZainuddinM. M.et al (2014). “Continuous Dynamic Bayesian Network for gene regulatory network modelling,” in 2014 International Conference on Computational Science and Technology (ICCST). 10.1109/iccst.2014.7045200
8
BeerM. A.TavazoieS. (2004). Predicting gene expression from sequence. Cell117, 185–198. 10.1016/s0092-8674(04)00304-6
9
BergenV.LangeM.PeidliS.WolfF. A.TheisF. J. (2020). Generalizing RNA velocity to transient cell states through dynamical modeling. Nat. Biotechnol.38, 1408–1414. 10.1038/s41587-020-0591-3
10
BertauxF.StomaS.DrasdoD.BattG. (2014). Modeling dynamics of cell-to-cell variability in TRAIL-induced apoptosis explains fractional killing and predicts reversible resistance. PLoS Comput. Biol.10, e1003893. 10.1371/journal.pcbi.1003893
11
BhattacharyaS.ZhangQ.AndersenM. E. (2011). A deterministic map of Waddington’s epigenetic landscape for cell fate specification. Bmc Syst. Biol.5, 85. 10.1186/1752-0509-5-85
12
BholaP. D.SimonS. M. (2009). Determinism and divergence of apoptosis susceptibility in mammalian cells. J. Cell Sci.122, 4296–4302. 10.1242/jcs.055590
13
BonoC. D.LiuY.FerrenaA.ValentineA.ZhengD.MorrowB. E. (2023). Single-cell transcriptomics uncovers a non-autonomous Tbx1-dependent genetic program controlling cardiac neural crest cell development. Nat. Commun.14, 1551. 10.1038/s41467-023-37015-9
14
BoumahdiS.de SauvageF. J. (2020). The great escape: tumour cell plasticity in resistance to targeted therapy. Nat. Rev. Drug Discov.19, 39–56. 10.1038/s41573-019-0044-1
15
BradyS. W.McquerryJ. A.QiaoY.PiccoloS. R.ShresthaG.JenkinsD. F.et al (2017). Combating subclonal evolution of resistant cancer phenotypes. Nat. Commun.8, 1231. 10.1038/s41467-017-01174-3
16
BrockA.ChangH.HuangS. (2009). Non-genetic heterogeneity — a mutation-independent driving force for the somatic evolution of tumours. Nat. Rev. Genet.10, 336–342. 10.1038/nrg2556
17
BurdziakC.Alonso-CurbeloD.WalleT.ReyesJ.BarrigaF. M.HavivD.et al (2023). Epigenetic plasticity cooperates with cell-cell interactions to direct pancreatic tumorigenesis. Science380, eadd5327. 10.1126/science.add5327
18
CabanosH. F.HataA. N. (2021). Emerging insights into targeted therapy-tolerant persister cells in cancer. Cancers13, 2666. 10.3390/cancers13112666
19
CaiadoF.Silva‐SantosB.NorellH. (2016). Intra‐tumour heterogeneity – going beyond genetics. Febs J.283, 2245–2258. 10.1111/febs.13705
20
CalboJ.MontfortE. V.ProostN.DrunenE. V.BeverlooH. B.MeuwissenR.et al (2011). A functional role for tumor cell heterogeneity in a mouse model of small cell lung cancer. Cancer Cell19, 244–256. 10.1016/j.ccr.2010.12.021
21
CastroJ. C.ValdésI.Gonzalez-GarcíaL. N.DaniesG.CañasS.WinckF. V.et al (2019). Gene regulatory networks on transfer entropy (GRNTE): a novel approach to reconstruct gene regulatory interactions applied to a case study for the plant pathogen phytophthora infestans. Theor. Biol. Med. Model16, 7. 10.1186/s12976-019-0103-7
22
ChaiL. E.LohS. K.LowS. T.MohamadM. S.DerisS.ZakariaZ. (2014). A review on the computational approaches for gene regulatory network construction. Comput. Biol. Med.48, 55–65. 10.1016/j.compbiomed.2014.02.011
23
ChakrabartiS.PaekA. L.ReyesJ.LasickK. A.LahavG.MichorF. (2018). Hidden heterogeneity and circadian-controlled cell fate inferred from single cell lineages. Nat. Commun.9, 5372. 10.1038/s41467-018-07788-5
24
ChanT. E.StumpfM. P. H.BabtieA. C. (2017). Gene regulatory network inference from single-cell data using multivariate information measures. Cell Syst.5, 251–267. 10.1016/j.cels.2017.08.014
25
ChanJ. M.Quintanal-VillalongaÁ.GaoV. R.XieY.AllajV.ChaudharyO.et al (2021). Signatures of plasticity, metastasis, and immunosuppression in an atlas of human small cell lung cancer. Cancer Cell39, 1479–1496.e18. 10.1016/j.ccell.2021.09.008
26
ChapmanM. P.RisomT.AswaniA. J.LangerE. M.SearsR. C.TomlinC. J. (2019). Modeling differentiation-state transitions linked to therapeutic escape in triple-negative breast cancer. Plos Comput. Biol.15, e1006840. 10.1371/journal.pcbi.1006840
27
ChenS.MarJ. C. (2018). Evaluating methods of inferring gene regulatory networks highlights their lack of performance for single cell gene expression data. Bmc Bioinforma.19, 232. 10.1186/s12859-018-2217-z
28
ChenZ.KingW. C.HwangA.GersteinM.ZhangJ. (2022). DeepVelo: single-cell transcriptomic deep velocity field learning with neural ordinary differential equations. Sci. Adv.8, eabq3745. 10.1126/sciadv.abq3745
29
ChisholmR. H.LorenziT.LorzA.LarsenA. K.AlmeidaL. N. D.EscargueilA.et al (2015). Emergence of drug tolerance in cancer cell populations: an evolutionary outcome of selection, nongenetic instability, and stress-induced adaptation. Cancer Res.75, 930–939. 10.1158/0008-5472.can-14-2103
30
CorreiaR. B.GatesA. J.WangX.RochaL. M. (2018). CANA: a Python package for quantifying control and canalization in boolean networks. Front. Physiol.9, 1046. 10.3389/fphys.2018.01046
31
de SandeB. V.FlerinC.DavieK.WaegeneerM. D.HulselmansG.AibarS.et al (2020). A scalable SCENIC workflow for single-cell gene regulatory network analysis. Nat. Protoc.15, 2247–2276. 10.1038/s41596-020-0336-2
32
DevarajV.BoseB. (2020). The mathematics of phenotypic state transition: paths and potential. J. Indian I Sci.100, 451–464. 10.1007/s41745-020-00173-6
33
DorigoM. (1994). “Learning by probabilistic Boolean networks,” in Proceedings of 1994 IEEE International Conference on Neural Networks (ICNN'94), 887–891. 10.1109/icnn.1994.374297
34
DunnS.LiM. A.CarbogninE.SmithA.MartelloG. (2019). A common molecular logic determines embryonic stem cell self‐renewal and reprogramming. EMBO J.38, e100003. 10.15252/embj.2018100003
35
DurenZ.ChenX.JiangR.WangY.WongW. H. (2017). Modeling gene regulation from paired expression and chromatin accessibility data. Proc. Natl. Acad. Sci.114, E4914–E4923. 10.1073/pnas.1704553114
36
ElowitzM. B.LevineA. J.SiggiaE. D.SwainP. S. (2002). Stochastic gene expression in a single cell. Science297, 1183–1186. 10.1126/science.1070919
37
FeinbergA. P.IrizarryR. A. (2010). Evolution in health and medicine Sackler colloquium: stochastic epigenetic variation as a driving force of development, evolutionary adaptation, and disease. Proc. Natl. Acad. Sci.107, 1757–1764. 10.1073/pnas.0906183107
38
FrickP. L.PaudelB. B.TysonD. R.QuarantaV. (2015). Quantifying heterogeneity and dynamics of clonal fitness in response to perturbation. 10.1002/jcp.24888
39
FriedrichM. J.NeriP.KehlN.MichelJ.SteigerS.KilianM.et al (2023). The pre-existing T cell landscape determines the response to bispecific T cell engagers in multiple myeloma patients. Cancer Cell41, 711–725.e6. 10.1016/j.ccell.2023.02.008
40
Gaspar-MaiaA.AlajemA.MeshorerE.Ramalho-SantosM. (2011). Open chromatin in pluripotency and reprogramming. Nat. Rev. Mol. Cell Bio12, 36–47. 10.1038/nrm3036
41
GayC. M.StewartC. A.ParkE. M.DiaoL.GrovesS. M.HeekeS.et al (2021). Patterns of transcription factor programs and immune pathway activation define four major subtypes of SCLC with distinct therapeutic vulnerabilities. Cancer Cell39, 346–360.e7. 10.1016/j.ccell.2020.12.014
42
GiadrossiS.DvorkinaM.FisherA. G. (2007). Chromatin organization and differentiation in embryonic stem cell models. Curr. Opin. Genet. Dev.17, 132–138. 10.1016/j.gde.2007.02.012
43
GlassL.KauffmanS. A. (1973). The logical analysis of continuous, non-linear biochemical control networks. J. Theor. Biol.39, 103–129. 10.1016/0022-5193(73)90208-7
44
GorinG.SvenssonV.PachterL. (2020). Protein velocity and acceleration from single-cell multiomics experiments. Genome Biol.21, 39. 10.1186/s13059-020-1945-3
45
GriffithsJ. A.ScialdoneA.MarioniJ. C. (2018). Using single‐cell genomics to understand developmental processes and cell fate decisions. Mol. Syst. Biol.14, e8046. 10.15252/msb.20178046
46
Grosse-WildeA.d’HérouëlA. F.McIntoshE.ErtaylanG.SkupinA.KuestnerR. E.et al (2015). Stemness of the hybrid epithelial/mesenchymal state in breast cancer and its association with poor survival. Plos One10, e0126522. 10.1371/journal.pone.0126522
47
GrovesS. M.IldefonsoG. V.McAteeC. O.OzawaP. M. M.IrelandA. S.StaufferP. E.et al (2022). Archetype tasks link intratumoral heterogeneity to plasticity and cancer hallmarks in small cell lung cancer. Cell Syst.13, 690–710.e17. 10.1016/j.cels.2022.07.006
48
GuoM.BaoE. L.WagnerM.WhitsettJ. A.XuY. (2016). SLICE: determining cell differentiation and lineage based on single cell entropy. Nucleic Acids Res.45, e54. 10.1093/nar/gkw1278
49
GuptaP. B.FillmoreC. M.JiangG.ShapiraS. D.TaoK.KuperwasserC.et al (2011). Stochastic state transitions give rise to phenotypic equilibrium in populations of cancer cells. Cell146, 633–644. 10.1016/j.cell.2011.07.026
50
HaghverdiL.BüttnerM.WolfF. A.BuettnerF.TheisF. J. (2016). Diffusion pseudotime robustly reconstructs lineage branching. Nat. Methods13, 845–848. nmeth. 10.1038/nmeth.3971
51
HanahanD. (2022). Hallmarks of cancer: new dimensions. Cancer Discov.12, 31–46. 10.1158/2159-8290.cd-21-1059
52
HayfordC. E.TysonD. R.RobbinsC. J.FrickP. L.QuarantaV.HarrisL. A. (2021). An in vitro model of tumor heterogeneity resolves genetic, epigenetic, and stochastic sources of cell state variability. Plos Biol.19, e3000797. 10.1371/journal.pbio.3000797
53
HéraultL.PoplineauM.DuprezE.RemyÉ. (2023). A novel Boolean network inference strategy to model early hematopoiesis aging. Comput. Struct. Biotechnol. J.21, 21–33. 10.1016/j.csbj.2022.10.040
54
HerringC. A.BanerjeeA.McKinleyE. T.SimmonsA. J.PingJ.RolandJ. T.et al (2018). Unsupervised trajectory analysis of single-cell RNA-seq and imaging data reveals alternative tuft cell origins in the gut. Cell Syst.6, 37–51. 10.1016/j.cels.2017.10.012
55
HersbachB. A.FischerD. S.MasserdottiG.DeekshaMojžišováK.WaltzhöniT.et al (2022). Probing cell identity hierarchies by fate titration and collision during direct reprogramming. Mol. Syst. Biol.18, e11129. 10.15252/msb.202211129
56
HillA. V. (1913). The combinations of haemoglobin with oxygen and with carbon monoxide. I. Biochem. J.7, 471–480. 10.1042/bj0070471
57
HuangS.KauffmanS. (2013). How to escape the cancer attractor: rationale and limitations of multi-target drugs. Semin. Cancer Biol.23, 270–278. 10.1016/j.semcancer.2013.06.003
58
HuangS.ErnbergI.KauffmanS. (2009). Cancer attractors: a systems view of tumors from a gene network dynamics and developmental perspective. Semin. Cell Dev. Biol.20, 869–876. 10.1016/j.semcdb.2009.07.003
59
HuangS. (2009). Non-genetic heterogeneity of cells in development: more than just noise. Development136, 3853–3862. 10.1242/dev.035139
60
HuangS. (2012). The molecular and mathematical basis of Waddington’s epigenetic landscape: a framework for post‐darwinian biology?Bioessays34, 149–157. 10.1002/bies.201100031
61
HuangS. (2013). Genetic and non-genetic instability in tumor progression: link between the fitness landscape and the epigenetic landscape of cancer cells. Cancer Metast Rev.32, 423–448. 10.1007/s10555-013-9435-7
62
Huynh-ThuV. A.IrrthumA.WehenkelL.GeurtsP. (2010). Inferring regulatory networks from expression data using tree-based methods. Plos One5, 12776. 10.1371/journal.pone.0012776
63
IstrailS.DavidsonE. H. (2005). Logic functions of the genomic cis-regulatory code. P Natl. Acad. Sci. U. S. A.102, 4954–4959. 10.1073/pnas.0409624102
64
IvanovP. Ch (2021). The new field of network physiology: building the human physiolome. Front. Netw. Physiol.1, 711778. 10.3389/fnetp.2021.711778
65
JainarayananA.Mouroug-AnandN.Arbe-BarnesE. H.BushA. J.Bashford-RogersR.FramptonA.et al (2023). Pseudotime dynamics of T cells in pancreatic ductal adenocarcinoma inform distinct functional states within the regulatory and cytotoxic T cells. Iscience26, 106324. 10.1016/j.isci.2023.106324
66
JiaD.JollyM. K.KulkarniP.LevineH. (2017). Phenotypic plasticity and cell fate decisions in cancer: insights from dynamical systems theory. Cancers9, 70. 10.3390/cancers9070070
67
JiaD.PaudelB. B.HayfordC. E.HardemanK. N.LevineH.OnuchicJ. N.et al (2020). Drug-tolerant idling melanoma cells exhibit theory-predicted metabolic low-low phenotype. Front. Oncol.10, 1426. 10.3389/fonc.2020.01426
68
JollyM. K.KulkarniP.WeningerK.OrbanJ.LevineH. (2018). Phenotypic plasticity, bet-hedging, and androgen independence in prostate cancer: role of non-genetic heterogeneity. Front. Oncol.8, 50. 10.3389/fonc.2018.00050
69
JooJ. I.ZhouJ. X.HuangS.ChoK.-H. (2018). Determining relative dynamic stability of cell states using boolean network model. Sci. Rep-uk8, 12077. 10.1038/s41598-018-30544-0
70
KalirS.AlonU. (2004). Using a quantitative blueprint to reprogram the dynamics of the flagella gene network. Cell117, 713–720. 10.1016/j.cell.2004.05.010
71
KamimotoK.StringaB.HoffmannC. M.JindalK.Solnica-KrezelL.MorrisS. A. (2023). Dissecting cell identity via network inference and in silico gene perturbation. Nature614, 742–751. 10.1038/s41586-022-05688-9
72
KangX.LiC. (2021). A dimension reduction approach for energy landscape: identifying intermediate states in metabolism‐EMT network. Adv. Sci.8, 2003133. 10.1002/advs.202003133
73
KauffmanS. (1969). Homeostasis and differentiation in random genetic control networks. Nature224, 177–178. 10.1038/224177a0
74
KauffmanS. A. (1969). Metabolic stability and epigenesis in randomly constructed genetic nets. J. Theor. Biol.22, 437–467. 10.1016/0022-5193(69)90015-0
75
KauffmanS. (1971). Chapter 5 gene regulation networks: a theory for their global structure and behaviors. Curr. Top. Dev. Biol.6, 145–182. 10.1016/s0070-2153(08)60640-7
76
KemperK.de GoejeP. L.PeeperD. S.AmerongenR. van (2014). Phenotype switching: tumor cell plasticity as a resistance mechanism and target for therapy. Cancer Res.74, 5937–5941. 10.1158/0008-5472.can-14-1174
77
KhanF. M.MarquardtS.GuptaS. K.KnollS.SchmitzU.SpitschakA.et al (2017). Unraveling a tumor type-specific regulatory core underlying E2F1-mediated epithelial-mesenchymal transition to predict receptor protein signatures. Nat. Commun.8, 198. 10.1038/s41467-017-00268-2
78
KimmelJ. C.PenlandL.RubinsteinN. D.HendricksonD. G.KelleyD. R.RosenthalA. Z. (2019). Murine single-cell RNA-seq reveals cell-identity- and tissue-specific trajectories of aging. Genome Res.29, 2088–2103. 10.1101/gr.253880.119
79
KumarN.CramerG. M.DahajS. A. Z.SundaramB.CelliJ. P.KulkarniR. V. (2019). Stochastic modeling of phenotypic switching and chemoresistance in cancer cell populations. Sci. Rep-uk9, 10845. 10.1038/s41598-019-46926-x
80
La MannoG.SoldatovR.ZeiselA.BraunE.HochgernerH.PetukhovV.et al (2018). RNA velocity of single cells. Nature560, 494–498. 10.1038/s41586-018-0414-6
81
LangJ.NieQ.LiC. (2021). Landscape and kinetic path quantify critical transitions in epithelial-mesenchymal transition. Biophys. J.120, 4484–4500. 10.1016/j.bpj.2021.08.043
82
LangeM.BergenV.KleinM.SettyM.ReuterB.BakhtiM.et al (2022). CellRank for directed single-cell fate mapping. Nat. Methods19, 159–170. 10.1038/s41592-021-01346-6
83
LangfelderP.HorvathS. (2008). Wgcna: an R package for weighted correlation network analysis. Bmc Bioinforma.9, 559. 10.1186/1471-2105-9-559
84
LiC.WangJ. (2014a). Landscape and flux reveal a new global view and physical quantification of mammalian cell cycle. 10.1073/pnas.1408628111
85
LiC.WangJ. (2014b). Quantifying the underlying landscape and paths of cancer. J. Roy. Soc. Interface11, 20140774. 10.1098/rsif.2014.0774
86
LiP.ZhangC.PerkinsE. J.GongP.DengY.-EdwardJPerkinsE. J. P. (2007). Comparison of probabilistic Boolean network and dynamic Bayesian network approaches for inferring gene regulatory networks. Bmc Bioinforma.8, S13. 10.1186/1471-2105-8-s7-s13
87
LiH.WangN.GongP.PerkinsE. J.ZhangC. (2011). Learning the structure of gene regulatory networks from time series gene expression data. Bmc Genomics12, S13. 10.1186/1471-2164-12-s5-s13
88
LiC.HongT.NieQ. (2016). Quantifying the landscape and kinetic paths for epithelial–mesenchymal transition from a core circuit. Phys. Chem. Chem. Phys.18, 17949–17956. 10.1039/c6cp03174a
89
LiaoD.Estévez-SalmerónL.TlstyT. D. (2012). Generalized principles of stochasticity can be used to control dynamic heterogeneity. Phys. Biol.9, 065006. 10.1088/1478-3975/9/6/065006
90
LiauB. B.SieversC.DonohueL. K.GillespieS. M.FlavahanW. A.MillerT. E.et al (2017). Adaptive chromatin remodeling drives glioblastoma stem cell plasticity and drug tolerance. Cell Stem Cell20, 233–246. 10.1016/j.stem.2016.11.003
91
LimJ. S.IbasetaA.FischerM. M.CancillaB.O’YoungG.CristeaS.et al (2017). Intratumoural heterogeneity generated by Notch signalling promotes small-cell lung cancer. Nature545, 360–364. 10.1038/nature22323
92
LiuF.ZhangS.-W.GuoW.-F.WeiZ.-G.ChenL. (2016). Inference of gene regulatory network based on local bayesian networks. Plos Comput. Biol.12, e1005024. 10.1371/journal.pcbi.1005024
93
LuoX.XuL.HanB.WangJ. (2017). Funneled potential and flux landscapes dictate the stabilities of both the states and the flow: fission yeast cell cycle. Plos Comput. Biol.13, e1005710. 10.1371/journal.pcbi.1005710
94
LytleN. K.FergusonL. P.RajbhandariN.GilroyK.FoxR. G.DeshpandeA.et al (2019). A multiscale map of the stem cell state in pancreatic adenocarcinoma. Cell177, 572–586. 10.1016/j.cell.2019.03.010
95
MargolinA. A.NemenmanI.BassoK.WigginsC.StolovitzkyG.FaveraR. D.et al (2004). ARACNE: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context. Bmc Bioinforma.7, S7–S2105. 10.1186/1471-2105-7-s1-s7
96
MarjanovicN. D.WeinbergR. A.ChafferC. L. (2013). Cell plasticity and heterogeneity in cancer. Clin. Chem.59, 168–179. 10.1373/clinchem.2012.184655
97
MarjanovicN. D.HofreeM.ChanJ. E.CannerD.WuK.TrakalaM.et al (2020). Emergence of a high-plasticity cell state during lung cancer evolution. Cancer Cell38, 229–246. 10.1016/j.ccell.2020.06.012
98
MarusykA.AlmendroV.PolyakK. (2012). Intra-tumour heterogeneity: a looking glass for cancer?10.1038/nrc3261
99
Masoudi-NejadA.BidkhoriG.AshtianiS. H.NajafiA.BozorgmehrJ. H.WangE. (2015). Cancer systems biology and modeling: microscopic scale and multiscale approaches. Semin. Cancer Biol.30, 60–69. 10.1016/j.semcancer.2014.03.003
100
MatsumotoH.KiryuH.FurusawaC.KoM. S. H.KoS. B. H.GoudaN.et al (2017). Scode: an efficient regulatory network inference algorithm from single-cell RNA-seq during differentiation. Bioinformatics33, 2314–2321. 10.1093/bioinformatics/btx194
101
MatsushitaY.LiuJ.ChuA. K. Y.Tsutsumi-AraiC.NagataM.AraiY.et al (2023). Bone marrow endosteal stem cells dictate active osteogenesis and aggressive tumorigenesis. Nat. Commun.14, 2383. 10.1038/s41467-023-38034-2
102
MeshorerE.MisteliT. (2006). Chromatin in pluripotent embryonic stem cells and differentiation. Nat. Rev. Mol. Cell Bio7, 540–546. 10.1038/nrm1938
103
MoermanT.SantosS. A.González-BlasC. B.SimmJ.MoreauY.AertsJ.et al (2018). GRNBoost2 and arboreto: efficient and scalable inference of gene regulatory networks. Bioinformatics35, 2159–2161. 10.1093/bioinformatics/bty916
104
MuP.ZhangZ.BenelliM.KarthausW. R.HooverE.ChenC.-C.et al (2017). SOX2 promotes lineage plasticity and antiandrogen resistance in TP53- and RB1-deficient prostate cancer. Science355, 84–88. 10.1126/science.aah4307
105
NathA.CosgroveP. A.MirsafianH.ChristieE. L.PfliegerL.CopelandB.et al (2021). Evolution of core archetypal phenotypes in progressive high grade serous ovarian cancer. Nat. Commun.12, 3039. 10.1038/s41467-021-23171-3
106
NeftelC.LaffyJ.FilbinM. G.HaraT.ShoreM. E.RahmeG. J.et al (2019). An integrative model of cellular states, plasticity, and genetics for glioblastoma. Cell178, 835–849. 10.1016/j.cell.2019.06.024
107
OlsenR. R.IrelandA. S.KastnerD. W.GrovesS. M.SpainhowerK. B.PozoK.et al (2021). ASCL1 represses a SOX9+ neural crest stem-like state in small cell lung cancer. Gene Dev.35, 847–869. 10.1101/gad.348295.121
108
OrenY.TsabarM.CuocoM. S.Amir-ZilbersteinL.CabanosH. F.HütterJ.-C.et al (2021). Cycling cancer persister cells arise from lineages with distinct programs. Nature596, 576–582. 10.1038/s41586-021-03796-6
109
PaudelB. B.HarrisL. A.HardemanK. N.AbugableA. A.HayfordC. E.TysonD. R.et al (2018). A nonquiescent “idling” population state in drug-treated, BRAF-mutated melanoma. Biophys. J.114, 1499–1511. 10.1016/j.bpj.2018.01.016
110
PiscoA. O.HuangS. (2015). Non-genetic cancer cell plasticity and therapy-induced stemness in tumour relapse: ‘What does not kill me strengthens me. Brit J. Cancer112, 1725–1732. 10.1038/bjc.2015.146
111
PiscoA. O.BrockA.ZhouJ.MoorA.MojtahediM.JacksonD.et al (2013). Non-Darwinian dynamics in therapy-induced cancer drug resistance. Nat. Commun.4, 2467. 10.1038/ncomms3467
112
PomeranceA.OttE.GirvanM.LosertW. (2009). The effect of network topology on the stability of discrete state models of genetic control. P Natl. Acad. Sci. U. S. A.106, 8209–8214. 10.1073/pnas.0900142106
113
PratapaA.JalihalA. P.LawJ. N.BharadwajA.MuraliT. M. (2020). Benchmarking algorithms for gene regulatory network inference from single-cell transcriptomic data. Nat. Methods17, 147–154. 10.1038/s41592-019-0690-6
114
QinS.JiangJ.LuY.NiceE. C.HuangC.ZhangJ.et al (2020). Emerging role of tumor cell plasticity in modifying therapeutic response. Signal Transduct. Target Ther.5, 228. 10.1038/s41392-020-00313-5
115
QiuX.MaoQ.TangY.WangL.ChawlaR.PlinerH. A.et al (2017). Reversed graph embedding resolves complex single-cell trajectories. Nat. Methods14, 979–982. 10.1038/nmeth.4402
116
QiuX.ZhangY.Martin-RufinoJ. D.WengC.HosseinzadehS.YangD.et al (2022). Mapping transcriptomic vector fields of single cells. Cell185, 690–711.e45. 10.1016/j.cell.2021.12.045
117
Quintanal-VillalongaÁ.ChanJ. M.YuH. A.Pe’erD.SawyersC. L.SenT.et al (2020). Lineage plasticity in cancer: a shared pathway of therapeutic resistance. Nat. Rev. Clin. Oncol.17, 360–371. 10.1038/s41571-020-0340-z
118
RamirezR. N.El-AliN. C.MagerM. A.WymanD.ConesaA.MortazaviA. (2017). Dynamic gene regulatory networks of human myeloid differentiation. Cell Syst.4, 416–429. 10.1016/j.cels.2017.03.005
119
RamirezD.KoharV.LuM. (2020). Toward modeling context-specific EMT regulatory networks using temporal single cell RNA-seq data. Front. Mol. Biosci.7, 54. 10.3389/fmolb.2020.00054
120
RisomT.LangerE. M.ChapmanM. P.RantalaJ.FieldsA. J.BonifaceC.et al (2018). Differentiation-state plasticity is a targetable resistance mechanism in basal-like breast cancer. Nat. Commun.9, 3815. 10.1038/s41467-018-05729-w
121
RussoM.PompeiS.SogariA.CoriglianoM.CrisafulliG.PuliafitoA.et al (2022). A modified fluctuation-test framework characterizes the population dynamics and mutation rate of colorectal cancer persister cells. Nat. Genet.54, 976–984. 10.1038/s41588-022-01105-z
122
SaadatpourA.AlbertR. (2013). Boolean modeling of biological regulatory networks: a methodology tutorial. Methods62, 3–12. 10.1016/j.ymeth.2012.10.012
123
SaelensW.CannoodtR.SaeysY. (2018). A comprehensive evaluation of module detection methods for gene expression data. Nat. Commun.9, 1090. 10.1038/s41467-018-03424-4
124
SaelensW.CannoodtR.TodorovH.SaeysY. (2019). A comparison of single-cell trajectory inference methods. Nat. Biotechnol.37, 547–554. 10.1038/s41587-019-0071-9
125
Sáez-AyalaM.MontenegroM. F.Sánchez-del-CampoL.Fernández-PérezM. P.ChazarraS.FreterR.et al (2013). Directed phenotype switching as an effective antimelanoma strategy. Cancer Cell24, 105–119. 10.1016/j.ccr.2013.05.009
126
Sanchez-CastilloM.BlancoD.Tienda-LunaI. M.CarrionM. C.HuangY. (2017). A Bayesian framework for the inference of gene regulatory networks from time and pseudo-time series data. 10.1093/bioinformatics/btx605
127
SchiebingerG.ShuJ.TabakaM.ClearyB.SubramanianV.SolomonA.et al (2019). Optimal-transport analysis of single-cell gene expression identifies developmental trajectories in reprogramming. Cell176, 928–943. 10.1016/j.cell.2019.01.006
128
SettyM.KiseliovasV.LevineJ.GayosoA.MazutisL.Pe’erD. (2019). Characterization of cell fate probabilities in single-cell data with Palantir. Nat. Biotechnol.37, 451–460. 10.1038/s41587-019-0068-4
129
ShafferS. M.DunaginM. C.TorborgS. R.TorreE. A.EmertB.KreplerC.et al (2017). Rare cell variability and drug-induced reprogramming as a mode of cancer drug resistance. Nature546, 431–435. 10.1038/nature22794
130
ShafferS. M.EmertB. L.HuerosR. A. R.CoteC.HarmangeG.SchaffD. L.et al (2020). Memory sequencing reveals heritable single-cell gene expression programs associated with distinct cellular behaviors. Cell182, 947–959. 10.1016/j.cell.2020.07.003
131
SharmaS. V.LeeD. Y.LiB.QuinlanM. P.TakahashiF.MaheswaranS.et al (2010). A chromatin-mediated reversible drug-tolerant state in cancer cell subpopulations. Cell141, 69–80. 10.1016/j.cell.2010.02.027
132
ShmulevichI.DoughertyE. R.KimS.ZhangW. (2002). Probabilistic boolean networks: a rule-based uncertainty model for gene regulatory networks. Bioinformatics18, 261–274. 10.1093/bioinformatics/18.2.261
133
SinghA.Saint-AntoineM. (2023). Probing transient memory of cellular states using single-cell lineages. Front. Microbiol.13, 1050516. 10.3389/fmicb.2022.1050516
134
SmithB. A.BalanisN. G.NanjundiahA.SheuK. M.TsaiB. L.ZhangQ.et al (2018). A human adult stem cell signature marks aggressive variants across epithelial cancers. Cell Rep.24, 3353–3366. 10.1016/j.celrep.2018.08.062
135
SpencerS. L.GaudetS.AlbeckJ. G.BurkeJ. M.SorgerP. K. (2009). Non-genetic origins of cell-to-cell variability in TRAIL-induced apoptosis. Nature459, 428–432. 10.1038/nature08012
136
StassenS. V.YipG. G. K.WongK. K. Y.HoJ. W. K.TsiaK. K. (2021). Generalized and scalable trajectory inference in single-cell omics data with VIA. Nat. Commun.12, 5528. 10.1038/s41467-021-25773-3
137
SteinwayS. N.ZañudoJ. G. T.MichelP. J.FeithD. J.LoughranT. P.AlbertR. (2015). Combinatorial interventions inhibit TGFβ-driven epithelial-to-mesenchymal transition and support hybrid cellular phenotypes. Npj Syst. Biol. Appl.1, 15014. 10.1038/npjsba.2015.14
138
SuY.WeiW.RobertL.XueM.TsoiJ.Garcia-DiazA.et al (2017). Single-cell analysis resolves the cell state transition and signaling dynamics associated with melanoma drug-induced resistance. Proc. Natl. Acad. Sci.114, 13679–13684. 10.1073/pnas.1712064115
139
SuY.BintzM.YangY.RobertL.NgA. H. C.LiuV.et al (2019). Phenotypic heterogeneity and evolution of melanoma cells associated with targeted therapy resistance. Plos Comput. Biol.15, e1007034. 10.1371/journal.pcbi.1007034
140
SuK.KatebiA.KoharV.ClaussB.GordinD.QinZ. S.et al (2022). NetAct: a computational platform to construct core transcription factor regulatory networks using gene activity. Genome Biol.23, 270. 10.1186/s13059-022-02835-3
141
SutherlandK. D.IrelandA. S.OliverT. G. (2022). Killing SCLC: insights into how to target a shapeshifting tumor. Gene Dev.36, 241–258. 10.1101/gad.349359.122
142
SwainP. S.ElowitzM. B.SiggiaE. D. (2002). Intrinsic and extrinsic contributions to stochasticity in gene expression. Proc. Natl. Acad. Sci.99, 12795–12800. 10.1073/pnas.162041399
143
TercanB.AguilarB.HuangS.DoughertyE. R.ShmulevichI. (2022). Probabilistic boolean networks predict transcription factor targets to induce transdifferentiation. Iscience25, 104951. 10.1016/j.isci.2022.104951
144
ThieffryD.ThomasR. (1998). Qualitative analysis of gene networks. Pac. Symposium Biocomput., 77–88.
145
TrairatphisanP.MizeraA.PangJ.TantarA. A.SchneiderJ.SauterT. (2013). Recent development and biomedical applications of probabilistic Boolean networks. Cell Commun. Signal11, 46. 10.1186/1478-811x-11-46
146
TrapnellC.CacchiarelliD.GrimsbyJ.PokharelP.LiS.MorseM.et al (2014). The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells. Nat. Biotechnol.32, 381–386. 10.1038/nbt.2859
147
Van BruggenD.PohlF.LangsethC. M.KukanjaP.LeeH.AlbiachA. M.et al (2022). Developmental landscape of human forebrain at a single-cell level identifies early waves of oligodendrogenesis. Dev. Cell57, 1421–1436.e5. 10.1016/j.devcel.2022.04.016
148
WaddingtonC. H. (1957). The strategy of the genes. Available at: https://archive.org/details/in.ernet.dli.2015.547782/page/n1/mode/2up (Accessed February 16, 2022).
149
WagnerD. E.KleinA. M. (2020). Lineage tracing meets single-cell omics: opportunities and challenges. Nat. Rev. Genet.21, 410–427. 10.1038/s41576-020-0223-2
150
WagnerD. E.WeinrebC.CollinsZ. M.BriggsJ. A.MegasonS. G.KleinA. M. (2018). Single-cell mapping of gene expression landscapes and lineage in the zebrafish embryo. Science360, 981–987. 10.1126/science.aar4362
151
WainwrightE. N.ScaffidiP. (2017). Epigenetics and cancer stem cells: unleashing, hijacking, and restricting cellular plasticity. Trends Cancer3, 372–386. 10.1016/j.trecan.2017.04.004
152
WangJ.CheungL. W.-K.DelabieJ. (2005). New probabilistic graphical models for genetic regulatory networks studies. J. Biomed. Inf.38, 443–455. 10.1016/j.jbi.2005.04.003
153
WangJ.HuangB.XiaX.SunZ. (2006). Funneled landscape leads to robustness of cell networks: yeast cell cycle. Plos Comput. Biol.2, e147. 10.1371/journal.pcbi.0020147
154
WangJ.XuL.WangE. (2008). Potential landscape and flux framework of nonequilibrium networks: robustness, dissipation, and coherence of biochemical oscillations. Proc. Natl. Acad. Sci.105, 12271–12276. 10.1073/pnas.0800579105
155
WangJ.LiC.WangE. (2010). Potential and flux landscapes quantify the stability and robustness of budding yeast cell cycle network. Proc. Natl. Acad. Sci.107, 8195–8200. 10.1073/pnas.0910331107
156
WangJ.ZhangK.XuL.WangE. (2011). Quantifying the Waddington landscape and biological paths for development and differentiation. Proc. Natl. Acad. Sci.108, 8257–8262. 10.1073/pnas.1017017108
157
WangL.ZhangQ.QinQ.TrasanidisN.VinyardM.ChenH.et al (2021). Current progress and potential opportunities to infer single-cell developmental trajectory and cell fate. Curr. Opin. Syst. Biol.26, 1–11. 10.1016/j.coisb.2021.03.006
158
WangJ. (2015). Landscape and flux theory of non-equilibrium dynamical systems with application to biology. Adv. Phys.64, 1–137. 10.1080/00018732.2015.1037068
159
WeinrebC.WolockS.TusiB. K.SocolovskyM.KleinA. M. (2018). Fundamental limits on dynamic inference from single-cell snapshots. Proc. Natl. Acad. Sci.115, E2467-E2476. 10.1073/pnas.1714723115
160
WeinrebC.Rodriguez-FraticelliA.CamargoF. D.KleinA. M. (2020). Lineage tracing on transcriptional landscapes links state to fate during differentiation. Sci. N. Y. N. Y.367, eaaw3381. 10.1126/science.aaw3381
161
WelchJ. D.HarteminkA.erJ.PrinsJ. F. (2016). SLICER: inferring branched, nonlinear cellular trajectories from single cell RNA-seq data. Genome Biol.17, 106. 10.1186/s13059-016-0975-3
162
WolfF. A.HameyF. K.PlassM.SolanaJ.DahlinJ. S.GöttgensB.et al (2019). Paga: graph abstraction reconciles clustering with trajectory inference through a topology preserving map of single cells. Genome Biol.20, 59. 10.1186/s13059-019-1663-x
163
WootenD. J.GrovesS. M.TysonD. R.LiuQ.LimJ. S.AlbertR.et al (2019). Systems-level network modeling of Small Cell Lung Cancer subtypes identifies master regulators and destabilizers. Plos Comput. Biol.15, e1007343. 10.1371/journal.pcbi.1007343
164
WootenD. J.QuarantaV. (2017). Mathematical models of cell phenotype regulation and reprogramming: make cancer cells sensitive again!Biochim. Biophys. Acta1867. 10.1016/j.bbcan.2017.04.001
165
WoutersJ.Kalender-AtakZ.MinnoyeL.SpanierK. I.WaegeneerM. D.González-BlasC. B.et al (2020). Robust gene expression programs underlie recurrent cell states and phenotype switching in melanoma. Nat. Cell Biol.22, 986–998. 10.1038/s41556-020-0547-3
166
WuW.WangJ. (2013a). Landscape framework and global stability for stochastic reaction diffusion and general spatially extended systems with intrinsic fluctuations. J. Phys. Chem. B117, 12908–12934. 10.1021/jp402064y
167
WuW.WangJ. (2013b). Potential and flux field landscape theory. I. Global stability and dynamics of spatially dependent non-equilibrium systems. J. Chem. Phys.139, 121920. 10.1063/1.4816376
168
XueR.ZhangQ.CaoQ.KongR.XiangX.LiuH.et al (2022). Liver tumour immune microenvironment subtypes and neutrophil heterogeneity. Nature612, 141–147. 10.1038/s41586-022-05400-x
169
YaboY. A.NiclouS. P.GolebiewskaA. (2021). Cancer cell heterogeneity and plasticity: a paradigm shift in glioblastoma. Neuro-oncology24, 669–682. 10.1093/neuonc/noab269
170
Yachie‐KinoshitaA.OnishiK.OstblomJ.LangleyM. A.PosfaiE.RossantJ.et al (2018). Modeling signaling‐dependent pluripotency with Boolean logic to predict cell fate transitions. Mol. Syst. Biol.14, e7952. 10.15252/msb.20177952
171
YanH.LiB.WangJ. (2019). Non-equilibrium landscape and flux reveal how the central amygdala circuit gates passive and active defensive responses. J. Roy. Soc. Interface16, 20180756. 10.1098/rsif.2018.0756
172
YuanS.NorgardR. J.StangerB. Z. (2019). Cellular plasticity in cancer. Cancer Discov.9, 837–851. 10.1158/2159-8290.cd-19-0015
173
YuhC.-H.BolouriH.DavidsonE. H. (1998). Genomic cis-regulatory logic: experimental and computational analysis of a sea urchin gene. Science279, 1896–1902. 10.1126/science.279.5358.1896
174
ZhouJ. X.HuangS. (2010). Understanding gene circuits at cell-fate branch points for rational cell reprogramming. Trends Genet.27, 55–62. 10.1016/j.tig.2010.11.002
175
ZhouP.LiT. (2016). Construction of the landscape for multi-stable systems: potential landscape, quasi-potential, A-type integral and beyond. J. Chem. Phys.144, 094109. 10.1063/1.4943096
176
ZhouJ. X.AliyuM. D. S.AurellE.HuangS. (2012). Quasi-potential landscape in complex multi-stable systems. J. R. Soc. Interface R. Soc.9, 3539–3553. 10.1098/rsif.2012.0434
177
ZhouJ. X.SamalA.d’HérouëlA. F.PriceN. D.HuangS. (2016). Relative stability of network states in Boolean network models of gene regulation in development. Biosystems142, 15–24. 10.1016/j.biosystems.2016.03.002
178
ZouM.ConzenS. D. (2005). A new dynamic Bayesian network (DBN) approach for identifying gene regulatory networks from time course microarray data. Bioinformatics21, 71–79. 10.1093/bioinformatics/bth463
179
ZouM.ToivanenR.MitrofanovaA.FlochN.HayatiS.SunY.et al (2017). Transdifferentiation as a mechanism of treatment resistance in a mouse model of castration-resistant prostate cancer. Cancer Discov.7, 736–749. 10.1158/2159-8290.cd-16-1174
Summary
Keywords
plasticity, gene regulatory networks, dynamical systems, ScRNA-seq, cancer, network physiology
Citation
Groves SM and Quaranta V (2023) Quantifying cancer cell plasticity with gene regulatory networks and single-cell dynamics. Front. Netw. Physiol. 3:1225736. doi: 10.3389/fnetp.2023.1225736
Received
19 May 2023
Accepted
25 August 2023
Published
04 September 2023
Volume
3 - 2023
Edited by
Mohit Kumar Jolly, Indian Institute of Science (IISc), India
Reviewed by
Chunhe Li, Fudan University, China
Shaon Chakrabarti, National Centre for Biological Sciences, India
Mingyang Lu, Jackson Laboratory, United States
Xiaojun Tian, Arizona State University, United States
Stefan Semrau, New York Stem Cell Foundation, United States
Updates
Copyright
© 2023 Groves and Quaranta.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Vito Quaranta, vito.quaranta@vanderbilt.edu; Sarah M. Groves, smgroves@virginia.edu
† Present addresses: Sarah M. Groves, Department of Biomedical Engineering, University of Virginia, Charlottesville, VA, United States
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.