Engineering principles for rationally design therapeutic strategies against hepatocellular carcinoma
Alexis Hernández-Magaña1 
Antonio Bensussen2 
Juan Carlos Martínez-García2* 
Elena R. Álvarez-Buylla1,3*- 1Instituto de Ecología, Universidad Nacional Autónoma de México, Ciudad de México, Mexico
- 2Departamento de Control Automático, Cinvestav-IPN, Ciudad de México, Mexico
- 3Centro de Ciencias de la Complejidad (C3), Universidad Nacional Autónoma de México, Ciudad de México, Mexico
The search for new therapeutic strategies against cancer has favored the emergence of rationally designed treatments. These treatments have focused on attacking cell plasticity mechanisms to block the transformation of epithelial cells into cancerous cells. The aim of these approaches was to control particularly lethal cancers such as hepatocellular carcinoma. However, they have not been able to control the progression of cancer for unknown reasons. Facing this scenario, emerging areas such as systems biology propose using engineering principles to design and optimize cancer treatments. Beyond the possibilities that this approach might offer, it is necessary to know whether its implementation at a clinical level is viable or not. Therefore, in this paper, we will review the engineering principles that could be applied to rationally design strategies against hepatocellular carcinoma, and discuss whether the necessary elements exist to implement them. In particular, we will emphasize whether these engineering principles could be applied to fight hepatocellular carcinoma.
Introduction
Currently, cancer is one of the most important chronic-degenerative diseases, since it is responsible for 10 million deaths per year (Ferlay et al., 2021). However, not all cancers affect patients in the same way; it has been reported that some neoplasms such as hepatocellular carcinoma (HCC) are particularly lethal (Rumgay et al., 2022). Thus, it is a priority to identify the molecular mechanisms involved in the malignant transformation of somatic cells and determine the factors that increase the mortality rate of cancers like HCC. In this sense, experimental reports have shown that epithelial-to-mesenchymal transition (EMT) is crucial to regulate the development of HCC (Yuan et al., 2019). Broadly speaking, EMT is a reversible process by which epithelial cells lose cell-cell adhesion, change their morphology, and acquire properties of mesenchymal cells such as migratory ability (Figure 1A) (Cicchini et al., 2015; Wickramaratne et al., 2021).
Figure 1. Dynamical properties of biological systems. Panel (A) shows a schematic representation of cell plasticity and biological robustness. The blue arrows represent the probability of remaining in a specific phenotype, and gray arrows represent the probability of transitioning to other phenotypes. Panel (B) represents Waddington’s conceptual model of the epigenetic landscape applied to liver cells. This figure shows how terminal phenotypes such as hepatocytes and mesenchymal cells can be obtained from a totipotent cell.
Recently, a novel therapeutic approach that rationally targets EMT plasticity has been proposed (Yang et al., 2019; Huang et al., 2022). This approach uses drugs to either redirect or block the EMT in patients (Huang et al., 2022). Unfortunately, this strategy has not achieved notable improvements in cancer therapy (Shi et al., 2023). It is thought that the underperformance of this strategy is due to the lack of a comprehensive understanding of mechanisms that regulate tumor cell plasticity (Yang et al., 2019; Boumahdi and de Sauvage, 2020; Wickramaratne et al., 2021). In consequence, it is essential to explore better design techniques capable of dealing with complex problems such as the regulation of EMT plasticity. In this sense, a discipline that has historically specialized in solving complex problems is engineering.
In fact, abstract engineering principles have been successfully applied to create new metabolic pathways and artificial biological processes with broad spectrum applications (Choi et al., 2019). Furthermore, the methodology of engineering has been used to optimize already existing metabolic, genetic, and signaling pathways (He et al., 2016). Interestingly, recent research has focused on rationally designing cancer treatments (Wu et al., 2019; Shih and Chen, 2022). In the present work, we will examine the challenges that a rational design against HCC should be able to face and we will analyze the clinical viability of this approach in our current context.
Controlling biological systems according to engineering
The rational approach that engineering follows to solve problems begins with the creation of models that represent the system that is sought to be controlled (Ogata and Brewer, 1971; Mellodge, 2015). Then, control laws with inputs are proposed to regulate the system and guide it to a state or states that are of interest. Next, the parameters and materials necessary to implement the created control laws are analyzed. Subsequently, devices are built that implement the control laws with the proposed specifications, and finally, the performance of the solution is experimentally tested (Ogata and Brewer, 1971; Mellodge, 2015). This approach has passed integrally to metabolic engineering and synthetic biology, allowing the creation of highly complex de novo biological processes (Siuti et al., 2014; He et al., 2016). In these engineering applications, it is particularly important to know the dynamics of the flows of matter in the system, this in order to optimize already existing processes (Haellman et al., 2021).
To practically illustrate the principles mentioned above, Ajo-Franklin and colleagues noticed that some protein production systems in yeast gradually lost efficiency at the time the producing cells entered cell division (Ajo-Franklin et al., 2007). In consequence, this research team set out to design a protein production system with memory, that is, with the ability to maintain protein production regardless of cell division state. To do this, they conceptualized two genes, a controller gene and a process gene. The controller gene can only be induced with galactose as input, and once induced, it will express an artificial transcription factor fused to RFP, which activates the expression of a second gene. As a result of the above, the second gene expressed another artificial transcriptional factor fused to YFP, which is capable of inducing its own expression. The idea was that, when yeast transfected with this genetic system were induced with galactose, they would begin to produce YFP stably over time, regardless of whether they divided or not (Ajo-Franklin et al., 2007).
To ensure this behavior, the team represented the interactions of the system with a mathematical model and upon analyzing it, they found that the affinity constants of the transcriptional factors require certain specific values to achieve the desired behavior. After several tests in the laboratory, they created the transcriptional factors with the desired parametric values, and thus they were able to construct their production system with the properties of automatically producing a protein of interest over time, without the need to become induced (Ajo-Franklin et al., 2007). The above demonstrates that, through the use of mathematical models, it is possible to control a biological system at the molecular level so that it performs a desired behavior.
Therapeutic strategies to control hepatocellular carcinoma
In this direction, pioneering works have begun to implement engineering principles to create therapeutic strategies to eliminate HCC. From a control perspective, these treatments focus on destroying tumor cells and preventing metastasis (Wu et al., 2019). In particular, Han et al designed a series of genetic circuits with the ability to detect hallmark markers of HCC such as AFP activation, and interactions between YAP1 and 14-3-3s. The basic assumption was that these markers are notably increased in HCC cells compared to healthy cells, and by having high levels of these three input signals, their genetic circuits could express an antibody against VEGF in order to prevent tumor metastasis. This design was tested experimentally in mice, and it was found that there was a significant decrease in the size of HCC tumors in the diseased mice (Han et al., 2020).
Continuing with this line, Yang and Ding improved this idea and designed a genetic circuit capable of selectively eliminating HCC tumor cells (Yang and Ding, 2021). Furthermore, they added an extra layer of safety, coupling the expression of the anti-apoptotic Bcl-2 gene with the hBAX gene. The underlying motivation was to prevent premature activation of hBAX expression and accidently activate the apoptosis pathway, in case basal levels of YAP1 exist. The decisive test for this circuit was performed by mixing the IMR90 and PLC/PRF/5 cell lines, the first being from lung, while the second comes from liver carcinoma. In this experiment it was found that the circuit selectively killed PLC/PRF/5 cells, leaving lung cells undamaged (Yang and Ding, 2021). However, despite the promising results obtained in vitro, the circuit was not able to eliminate the entire population of cancer cells, which shows that it is necessary to improve the design before implementing in vivo applications (Yang and Ding, 2021). This same problem was found by Lin and colleagues, who designed a genetic circuit that used miR-196a and miR-126 as inputs to induce hBAX expression in mice with HCC (Lin et al., 2018). As proof-of-principle, these approaches demonstrate that it is possible to control HCC tumor cells. However, the way in which these control designs are conceptualized still needs to be improved, and to do so it could be necessary to delve deeper into the functioning of tumor cells.
Understanding the state space of a biological system
In the example of Ajo-Franklin and colleagues, it can be seen that they guided their genetic system toward a specific state, and the task was relatively simple because the system was small (Ajo-Franklin et al., 2007). However, controlling complex systems such as the regulatory networks that underlie the cell plasticity (Huang et al., 2022; Shi et al., 2023) is complicated due to the multistability of these biological systems, that is, the existence of several operating stable states (i.e., attractors) (Laurent and Kellershohn, 1999). In this regard, one way to represent the biological multistability was introduced by Conrad Hal Waddington, who proposed a conceptual model, known as the “epigenetic landscape”, that represents all the possible differentiation routes that a totipotent cell could take (Rajagopal and Stanger, 2016a) (Figure 1B).
In this model, a stem cell represented by a ball, placed at the top of the landscape would roll down and reach the bottom of valleys that are equivalent to stable cell phenotypes (e.g., cell types or cell fates), (Ladewig et al., 2013; Jia et al., 2017) (Figure 1B); these routes through the epigenetic landscape corresponds to cellular differentiation pathways. Although the epigenetic landscape metaphor was originally proposed to illustrate development, it is also useful to depict phenotypic transitions (i.e., cell plasticity) occurring in adult tissues during regeneration or cancer progression (Rajagopal and Stanger, 2016b; Pisco and Huang, 2015). It has been shown that differentiated cells can also move towards another differentiated state, a process called trans-differentiation, for instance, epithelial-mesenchymal transition (Jopling et al., 2011). Also, differentiated cells can transit towards less differentiated states (i.e., dedifferentiation) (Jopling et al., 2011), for example, hepatocytes can convert into liver progenitor cells in response to chronic injury (Li et al., 2020); even differentiated cells can transit to a pluripotent state (i.e., reprogramming) in response to overexpression of transcription factors such as OCT4, SOX2, KLF4, and MYC (Jopling et al., 2011). In this context, the heights of hills represent the difficulty of jumping from one attractor to another (Pisco and Huang, 2015). As discussed later, the enormous regenerative capacity of the liver depends significantly on trans-differentiation and dedifferentiation processes.
Notably, these phenotypic transitions are constrained by an underlying gene regulatory network (GRN). For example, various regulatory circuits controlling the EMT have been described experimentally (Saunders and McClay, 2014; Migault et al., 2022). In broad terms, epithelial markers (such as mir34a, mir200b, miR00c) and mesenchymal markers (such as SNAI1/2, TWIST1/2, ZEB1/2, FOXC2) form two distinct modules or “teams” (Hari et al., 2022). It has been observed that nodes within the same “team” activate each other, while inhibition is observed among members of different “teams”, enabling a switch behavior during EMT (Hari et al., 2022) and promoting the stability of phenotypes, either epithelial or mesenchymal, once reached.
Although GRNs underlying biological processes are extremely complex, it is possible to characterize their dynamical behavior and quantify their associated epigenetic landscape with the help of mathematical models. In this sense, mathematical abstractions, like differential equations, can be used to represent experimental information about the regulation of nodes (e.g., genes, proteins, miRNAs) and explore quantitively Waddington’s epigenetic landscape (Karlebach and Shamir, 2008; Álvarez-Buylla Roces et al., 2018). The specific methodology for deriving the epigenetic landscape from dynamical GRN models has been reviewed previously (Davila-velderrain et al., 2015; Álvarez-Buylla Roces et al., 2018). Thus, from a dynamical systems perspective, the epigenetic landscape graphically represents the topography of the state space (i.e., the set of possible states of the system) of a GRN, where stable states are at the bottom of valleys and unstable states at the top of hills (Huang, 2013). According to the notion of stability, a cell will tend to maintain its phenotype despite the presence of relatively slight environmental fluctuations or noise in gene expression (Huang, 2013).
Importantly, theoretical-experimental evidence has proven that Waddington´s epigenetic landscape can be quantified at the genome scale. Specifically, Shi and collaborators analyzed several independent sets of RNA-seq data; for integrating the transcriptomic profiles with a protein-protein interaction network created from interactions reported in databases and using specialized algorithms to quantify the levels of network entropy, it was observed the appearance of stable gene expression patterns that corresponded to the terminal phenotypes of somatic cells (Shi et al., 2018). In addition, they were able to trace all the differentiation routes that pluripotent cells follow to reach terminal phenotypes such as hepatocytes (Shi et al., 2018).
These results open the possibility of using the epigenetic landscape as a blueprint to control cell plasticity. It would be used to explore alterations in gene expression favoring that a specific cell fate can be attained, for example, senescence, apoptosis, or any other state that inhibits proliferation, which could be a powerful tool against diseases such as cancer, particularly in HCC.
Epigenetic landscape as a strategic blueprint
The next question to address is how to use information about the dynamics of biological systems as a design blueprint. The answer to this question may lie in recent advances in metabolic engineering. Specifically, Buldum and colleagues sought to improve the prediction of the yields of bioprocesses carried out by microorganisms (Buldum et al., 2020). To do this, they used as a case study the creation and optimization of a genetic circuit to induce the bacterial cellulose synthesis pathway in Escherichia coli. They used a GRN to integrate into the design all aspects of the biological multistability and molecular complexity of the pathway under natural conditions. With this approach, they were able to accurately predict and significantly improve the performance of transgenic bacteria (Buldum et al., 2020). Learning from the lessons given by metabolic engineering, it follows that the key to creating biological control systems focused against HCC is precisely to learn about the molecular mechanisms that control the appearance and progression of this disease. Therefore, all strategies that seek to improve the treatments against HCC must consider the diverse internal (e.g., GRNs, mutations) and external factors (e.g., immune system activity) that can influence cell plasticity and could hinder the design work. However, a comprehensive understanding of all these factors represents a major challenge.
Physiological context of the liver
Derived from the previous point, it is necessary to mention important aspects of liver physiology. In particular, this organ is composed of hepatocytes (epithelial cells that constitute more than 80% of the liver mass), cholangiocytes (biliary epithelial cells), endothelial cells, Kupffer cells, natural killer cells (NKs), and mesenchymal cells known as hepatic stellate cells (HSCs), whose chronic activation promotes liver fibrosis (Daoust and Cantero, 1959; Roskams, 2006; Higashi et al., 2017; Michalopoulos and Bhushan, 2021) (Figure 2A). Each cell type would correspond to a valley in the epigenetic landscape.
Figure 2. Homeostatic role of cell plasticity in the liver. Panel (A) shows a schematic representation of the liver under homeostatic conditions. Some of the main cell types are also shown in their anatomical location within the liver. Panel (B) shows the functioning and main features of the epigenetic landscape of normal liver cells.
Considering the phenotype of hepatocytes in more detail, most hepatocytes are in quiescence (G0) and show a particular gene expression pattern including transcription factors such as HNF4a and HNF1a (Berasain and Avila, 2015). However, the proliferation of hepatocytes can be triggered by factors like liver injury, diminished cell density, or cell-cell interaction loss (Loyer et al., 1996). It has been reported that the co-culture of hepatocytes with cholangiocytes inhibits the response to growth factors (Loyer et al., 1996). Interestingly, hepatocytes in culture suffer dedifferentiation; it was observed that a stiff extracellular matrix promotes EMT and apoptosis resistance through the activation of Akt and ERK1/2 signaling (Godoy et al., 2009). However, these changes were reversible when hepatocytes were cultured in soft gel collagen (Godoy et al., 2009). Another study has shown that the transcriptional profile of primary cultured hepatocytes is similar to those of liver diseases such as cirrhosis and HCC (Godoy et al., 2016). Particularly, HNF4a downregulation has been related to the loss of specific cell-type features in cultured hepatocytes and HCC (Godoy et al., 2016; Lau et al., 2018).
Remarkably, liver cells exhibit outstanding phenotypic plasticity in response to injury. From the epigenetic landscape approach, these phenotypic transitions can be represented as jumps between valleys. Examples of this include the differentiation of liver progenitor cells (LPCs) into hepatocytes or cholangiocytes, as well as the differentiation of hepatocytes and cholangiocytes into LPCs, and the trans-differentiation between hepatocytes and cholangiocytes (Ko et al., 2020; Li et al., 2020; So et al., 2020). It has also been observed that HSCs show mesenchymal stem cell properties and can also give rise to hepatocytes and cholangiocytes (Kordes et al., 2015; Kordes et al., 2021). In addition, hepatocytes can undergo partial EMT during liver regeneration (Oh et al., 2018). Notably, these phenotypic transitions enable the robustness of liver functions against various types of damage. For example, in conditions of chronic liver damage, in which the population of hepatocytes has been severely impaired, the conversion of cholangiocytes into hepatocytes can contribute significantly to liver regeneration (Deng et al., 2018; Li et al., 2020).
Thus, understanding HCC requires considering the intrinsic plasticity that liver cells possess. In this sense, epigenetic landscape modeling could be essential to map the stable phenotypes emerging from the GRN (Figure 2B).
Molecular bases of the EMT plasticity
Due to the relevance of EMT in the progression of HCC (Cicchini et al., 2015), we will focus our analysis on the regulation of this process. Molecularly, EMT plasticity is originated by several double negative feedback loops between epithelial and mesenchymal markers. For example, miR-200 and ZEB1, mir-34 as well as SNAIL, mir-200 and SNAIL (Garibaldi et al., 2012; Migault et al., 2022). Also, mutual negative regulation between SNAIL and HNF4A has also been identified in hepatocytes (Cicchini et al., 2015). Although these mutual inhibition interactions are related to the bistable behavior of EMT (Huang et al., 2017) they only partially control it. To understand exactly how EMT works, it is important to remember that this phenotypic transition is coupled to other cell processes like proliferation, senescence, stemness, and inflammation.
Multiple interactions have been reported among proliferation regulators such as β-catenin, YAP1, and c-Myc. Double negative feedback loops have been identified between Wnt- β-catenin signaling and HNF4A (Yang et al., 2013), as well as YAP1 and HNF4A (Cai et al., 2017). In addition, it has been found that the Wnt-β-catenin pathway positively regulates Snail (Jong et al., 2005) and mutual activation between YAP1 and SNAIL (Noce et al., 2019). It has also been reported that c-Myc positively regulates the SNAIL expression (Smith et al., 2009). HNF4A negatively regulates proliferation through repression of c-Myc and Cyclin D1 (Walesky and Apte, 2015). As mentioned above, the decreased expression of liver-specific genes such as HNF4A has been observed in cultured hepatocytes and HCC cells (Berasain and Avila, 2015). In this regard, over-expression of cell type-specific transcription factors decreased the efficiency of the dedifferentiation process induced by stem-related transcription factors in a neural-lineage cell line and hepatoblasts (Hikichi et al., 2013). In general, it is known that cell differentiation and the cell cycle are closely related through mutual inhibition between cell type-specific and cell cycle regulatory factors (Ruijtenberg and van den Heuvel, 2016; Zhao et al., 2020).
On the other hand, senescence has an influence on EMT and vice versa (Kishi et al., 2015). p53 activates the expression of miR200c and mir34a, which are associated with the conservation of the epithelial phenotype (Bommer et al., 2007; Chang et al., 2011). The loss of p53 decreases the expression of miR200c and mir34a and it induces EMT (Chang et al., 2011; Kim et al., 2011). At the post-transcriptional level, SNAIL and p53 negatively regulate each other (Lee et al., 2009; Lim et al., 2010; Ni et al., 2016). Moreover, it has been reported that SLUG suppresses the expression of p16 (Zhu et al., 2019).
The relationship between EMT and stemness has been also studied. Although OCT4 and NANOG activate SNAIL expression to induce EMT (Yin et al., 2015), the regulation is not simple since NANOG can also activate SLUG expression (Yao et al., 2016) and SNAIL can suppress NANOG expression promoting stem cell differentiation (Galvagni et al., 2015). In hepatocytes, HNF4 and SNAIL control the expression of stem markers through the regulation of miR-200 and miR-34a (Garibaldi et al., 2012). It is interesting that it has been reported that OCT/SOX2 activate miR200s family expression during fibroblast reprogramming (Wang et al., 2013) and that the miR-200 family negatively regulates SOX2 and KLF4 (Wellner et al., 2009; Pandey et al., 2015). It is known EMT plays a role in the emergence of cancer stem cells (Mani et al., 2008; Morel et al., 2008; May et al., 2011; Andriani et al., 2016; Jayachandran et al., 2016; Harner-Foreman et al., 2017; Jing et al., 2019) and in the reprogramming of somatic cells (Esteban et al., 2012; Shu and Pei, 2014).
Finally, experimental evidence shows that inflammation and EMT are also closely coordinated (López-Nouoa and Nieto, 2009). In particular, the IL-6/STAT3 pathway and HNF4A negatively regulate each other through miRNAs (Yan et al., 2018). Notably, a miRNA-mediated mutual inhibition circuit was also reported between NF-κB and HNF4A (Ning et al., 2014). Conversely, NF-κB activates the expression of mesenchymal markers such as SNAIL, SLUG, ZEB1, and TWIST (Min et al., 2008). Since HCC and other types of cancer arise and develop in a chronic inflammation context, the influence of inflammation on EMT and other processes is especially important (Taniguchi and Karin, 2018; Yan et al., 2018; Czauderna et al., 2019).
Collectively, these interactions highlight the tremendous complexity of the intracellular regulatory network underlying the EMT plasticity. In this context, mathematical-computational modeling has provided a framework to integrate the vast amount of available information and explain some aspects of cell plasticity, as we will see in the next sections.
Cancer stem cells are originated by plasticity
At this point it is clear that there is enough information to construct GRNs models. However, it is important to define and identify the types of tumor cells to be studied (Figure 3A). In this direction, tumors can be phenotypically heterogeneous, which means that tumor cells may exhibit variations in phenotypic traits such as gene expression, metabolism, motility, stemness, proliferation, and metastatic capacity (Marusyk and Polyak, 2010).
Figure 3. Emergence of CSCs due to cell plasticity alterations. Panel (A) shows a schematic representation of the anatomy of a tumor produced by HCC. This tumor shows the niche of the CSCs where they are protected from the external environment. Of importance is that the CSCs are heterogenous at phenotypic level, since they may present different expression levels of stem-cell markers (pink stars), hepatocyte markers (blue dots) or mesenchymal markers (orange squares). Panel (B) shows a representation of the effect of mutations and disturbances produced by the microenvironment on the topology of the epigenetic landscape. Collectively, those perturbations transform the topology of the epigenetic landscape to generate altered phenotypes such as CSCs (pink arrows).
In particular, tumor cells that present stem-like phenotypes, known as cancer stem cells (CSCs), are being intensively investigated. Notably, such cells present subpopulations that show different expression levels of stemness, epithelial and mesenchymal markers, and different proliferative capacities (Mohan et al., 2021). In this regard, phenotypic heterogeneity has important functional implications: quiescent CSCs or CSCs with a low cell division rate are related to therapy resistance and recurrence due to most chemotherapeutic agents target the proliferative cells, while CSCs with a hybrid epithelial and mesenchymal phenotype are highly prone to initiate metastasis (Mohan et al., 2021). Interestingly, the frequency of CSCs has been observed to increase according to the stage of tumor development (Yamashita and Wang, 2013).
Although several sources of intratumoral heterogeneity have been identified, including mutations, molecular epigenetic changes, and signals from the microenvironment (Marusyk et al., 2020), these factors are frequently studied separately, which may limit the understanding of intratumoral heterogeneity and its implications. In this sense, the epigenetic landscape approach could serve as an integrative framework.
Typically, somatic mutations are assumed to be the major source of intratumoral heterogeneity while network dynamic that controls phenotypic plasticity is ignored. In this regard, the integration of genetic and non-genetic mechanisms that control cell plasticity could give us a greater understanding of tumor cell plasticity and intratumoral heterogeneity (Pisco and Huang, 2015; Huang, 2020). It has been suggested that mutations may alter the topography of the epigenetic landscape, potentially changing the barriers between different valleys. This, in turn, could modify the probabilities of phenotypic transitions (Figure 3B) (Huang, 2013; Huang, 2020). Considering the transition into stem-like states, experimental observations might agree with this. For instance, it has also been observed that normal, immortalized, or transformed human breast epithelial cells can convert into CD44+ cells, which is a marker of stemness (Chaffer et al., 2011). However, transformed cells overexpressing RAS, telomerase, and a viral oncoprotein show an increased ability to become stem-like cells compared to non-transformed cells (Chaffer et al., 2011).
At least two mechanisms have been proposed for CSCs emergence in HCC: the transformation of liver progenitor cells and the differentiation of epithelial cells (Nio et al., 2017; Yagci et al., 2017; Wu et al., 2020). Indeed, several results indicate that the appearance of stem-cell markers is a result of EMT in HCC and other cancer types (Giannelli et al., 2016; Chen et al., 2017; Forte et al., 2017; Jing et al., 2019). For instance, Mani et al demonstrated that EMT can be induced either by TWIST, SNAIL ectopic expression or by TGF-β treatment, in immortalized human mammary cells, promoted the increase of the subpopulation CD44+. Remarkably, such cells were able to form mammospheres whereas CD44− cells were unable to do so (Mani et al., 2008).
On the other hand, environmental factors such as inflammation and hypoxia also promote the acquisition of stem cell traits (Figure 3B). IL-6 treatment induces EMT and the acquisition of stem cell properties in breast cancer cells (Xie et al., 2012). Similarly, IL-6 secreted by tumor-associated macrophages promoted the expansion of the CD44+ cell population in HCC (Wan et al., 2014). Also, hypoxia upregulates OCT4, NANOG, and BMI-1 in HCC and HIF-1α knockdown inhibited hypoxia-induced stem-like properties (Cui et al., 2017). It is of interest that both IL-6 and IL-1β upregulate HIF-1α expression (Wan et al., 2014; Zhang et al., 2018). Complementary observations obtained from modelling of EMT and stemness modules indicate that NF-κB might stabilize the epithelial-mesenchymal hybrid phenotype and increase the likelihood of gaining stem-like properties (Jolly et al., 2014).
On the whole, these data suggest that there are three main stable states: epithelial cells, mesenchymal cells, and stem-like cells, that must be considered to explore the phenotypic plasticity involved in CSC emergence. Using a quantitative epigenetic landscape model, one would identify how mutations and environmental signals specifically change phenotypic transition probabilities. Similarly, targets could be found that inhibit undesirable phenotypic transitions.
Cancer-associated senescence is a plastic state
In addition to the interactions and the characterization of the stable states of a network, it is necessary to know if there is evidence that proves that the system can leave or transit toward a state of interest. We will examine the case of senescence for its relevance in cancer progression.
Traditionally, senescence has been defined as a state of the permanent arrest of the cell cycle induced by several types of stress, such as severe DNA damage or oncogene activation (Herranz et al., 2018; Takasugi et al., 2022). Furthermore, senescent cells are unresponsive to growth factors, resistant to apoptosis, and secrete a mix of proinflammatory cytokines, chemokines, growth factors, and proteases called senescence-associated secretory phenotype (SASP) (Herranz et al., 2018; Paez Ribes et al., 2019). Senescence is induced in response to various stressors and prevents the proliferation of damaged cells, thus, it is considered a tumor suppressor mechanism. However, the accumulation of senescent cells in tumors may also favor a pro-inflammatory and pro-proliferative microenvironment via the SASP (Herranz et al., 2018). It suggests the manipulation of senescence must rely on context.
Although the senescent state is a very stable phenotype, the transition from senescent phenotype to proliferative phenotype has been observed experimentally (Lee and Schmitt, 2019; Saleh et al., 2019; Walcher et al., 2020). In particular, it has been reported that colon cancer cells with p53 and p16 inactivation were able to enter senescence (Carné Trécesson et al., 2011). Nevertheless, these senescent cells resumed proliferation after 4–5 weeks. Interestingly, such cells showed upregulation of vimentin and downregulation of E-cadherin, which is associated with EMT induction (Carné Trécesson et al., 2011).
In a similar way, senescent fibroblast attained escape from oncogene-induced senescence after a prolonged period in culture, because of telomerase reactivation (Patel et al., 2016). It is interesting that senescence escape of lymphoma cells has been associated with acquiring stem cell features (Milanovic et al., 2018). Remarkably, the reprogramming of senescent fibroblasts into induced pluripotent stem cells using OCT4, SOX2, KLF4, c-Myc, NANOG, and LIN28 has been experimentally demonstrated (Lapasset et al., 2011).
Together, these data indicate that even states like senescence can be reversible through manipulation (Figure 4). Importantly, the immune system is responsible for removing the senescent cells (Herranz et al., 2018), but it also plays a relevant role on EMT.
Figure 4. Observable changes produced by alterations of the epigenetic landscape. This figure shows a simplification of the epigenetic landscape in two dimensions. Under normal conditions, the valleys of senescent phenotypes are deeper than the valleys of proliferative phenotypes, this means that the proliferative cells tend to enter senescence more easily. However, mutations affect the topology of the landscape by increasing the valleys depth of proliferative phenotypes, which means that there will be more cells going into proliferation and leaving senescence.
Effect of perturbations produced by the immune system on EMT
Finally, it is necessary to put into perspective the effect of perturbations produced by external factors involved in the generation of cancer. In this direction, the immune system is essential to determine cancer progression. To be precise, recent evidence suggests that the main factor against tumor recovery is the time it takes for the immune system to control the neoplasia. For instance, it has been observed that the longer it takes NK cells to kill cancer cells, the more likely they are to survive (Ichise et al., 2022) (Figure 5A). Additionally, if the immune system takes longer to eliminate tumor cells, CCL2 levels increase in the tumor microenvironment (Qiu et al., 2018), which will favor macrophage infiltration. This is particularly important for HCC, since the immune cells with the greatest presence in the liver are precisely NK cells (Sajid et al., 2022).
Figure 5. First stages of HCC development. Panel (A) shows the early stages of HCC formation. NK cells are able to efficiently detect and eliminate transformed cells by apoptosis or by the release of Granzyme (B). However, if NK are slow to kill tumor cells, NK cells will be exhausted and allow the emergence of CSCs. Panel (B) shows the establishment of the tumor microenvironment. After the formation of CSCs, the internalization of macrophages is favored, which can still kill tumor cells. However, if they fail to kill the tumor in a short time, it is possible that the differentiation of M1 macrophages towards the TAM phenotype is induced, favoring the establishment of a pro-tumorigenic microenvironment.
In this scenario, the chances of newly infiltrating macrophages becoming tumor-associated macrophages (TAMs) are increased due to interactions with cancer cells and tumor microenvironment (Zhou et al., 2020) (Figure 5B). Concerning TAMs, they present two main variants, the M1-TAM phenotype and the M2-TAM phenotype (Arvanitakis et al., 2022). Like their non-tumor-associated counterparts, M1-TAMs promote tumor elimination while M2-TAMs produces an anti-inflammatory microenvironment, which promotes tumor progression and metastasis (Arvanitakis et al., 2022). TAMs make up more than 50% of the cells that make up the periphery of tumors (Azizi et al., 2018), and depending on the prevalence of the subtype of TAMs, the tumor may compromise the life of the patient. Specifically, a high prevalence of M2-TAMs is associated to bad prognosis for patients (Arvanitakis et al., 2022).
The importance of TAMs in the maintenance of tumors is mainly due to the fact that they produce an anti-inflammatory microenvironment rich in cytokines such as IL-10 and TGF-β, while these macrophages maintain low levels of proinflammatory cytokines such as TNF-α, IL-12, IL-6 and IL-1β (Feng et al., 2022). As a consequence, it has been reported that tumor microenvironment severely depletes the effector phenotypes of the immune system, causing both CD4+ T cells and CD8+ T cells to become anergic and express surface markers such as PD-1, facilitating their clearance (Ma et al., 2019; Di Blasi et al., 2020). On the other hand, TAMs also activate the growth pathways of tumor cells, maintaining an optimal microenvironment for their proliferation. Therefore, TAMs also induce EMT of epithelial cells, promoting their malignant transformation (Gifford et al., 2017; Ran et al., 2021). Furthermore, TAMs directly kill effector cells express the death ligands FasL and PD-L1 on their membranes, which induce apoptosis in effector cells (Feng et al., 2022) (Figure 6A).
Figure 6. Progression of HCC to metastasis. Panel (A) shows the conversion process of Th1 cells into Treg cells. Through the anti-inflammatory cytokines secreted by TAMs, Th1 cells lose their effectiveness in fighting the tumor, joining the tumor microenvironment as Treg cells. Panel (B) shows the last step before starting metastasis. In case that T lymphocytes cannot neutralize the tumor, the last cells to participate are the neutrophils, which, due to the action of the microenvironment formed by tumor cells, TAMs and Treg cells, end up joining the tumor microenvironment in the form of TANs. Hereafter, it is not known exactly what causes TANs to secrete their NETs, favoring the detachment of tumor cells to initiate metastasis.
Moreover, the sustained activity of TAMs converts effector immunological cells into regulatory phenotypes unable to inhibit the tumor. TAMs also promote the internalization of Myeloid inhibitory cells (MDSCs), which enhances the anti-inflammatory microenvironment and proliferation of the tumor through the activation of growth pathways (Pan et al., 2020). As a consequence of this prolonged stimulation, tumor cells increase the secretion of CXCL5 (Zhou et al., 2016), which is a localization signal for neutrophils. Once the neutrophils infiltrate in the affected tissue and interact with TAMs as well as the tumor microenvironment, they loss their effector phenotype and become tumor-associated neutrophils (TANs) (Wu and Zhang, 2020).
Analogous to TAMs, TANs remain at the periphery of the tumor, although the increase in its presence is negatively correlated with patient survival. The TANs, like the TAMs, also have N1-TAN and N2-TAN phenotypes, the latter being anti-inflammatory in nature and particularly harmful to the patient’s life (Arvanitakis et al., 2021). Recent research found a pernicious dynamic caused by TANs. First, it has been documented that N2-TANs positively regulate the appearance of regulatory phenotypes of CD4+ T lymphocytes by secreting CCL17 (Shaul et al., 2016). Likewise, it has been observed that TANs promote the appearance of TAMs and promote the growth of cancer cells, which generates an increase in CXCL5, proliferating the presence of TANs in the tumor with a positive feedback loop (Zhou et al., 2016). The increase in TANs produces the release of Neutrophil Extracellular Traps (NETs), which promote the detachment of tumor cells from the tissue (Masucci et al., 2020), and together with the increase in the expression of VEGF and other angiogenic signals by TAMs (Pan et al., 2020), gives rise to tumor metastasis (Figure 6B). In particular, it has been seen that the infiltration of TANs is very strong in the context of HCC, which associates these cells directly with the progression and worsening of the clinical condition of the patients (Zhou et al., 2016).
State of the art of cancer modeling
The data collected herein demonstrate that there is sufficient information to create models that adequately represent EMT. However, it remains to be determined whether there are tools capable of capturing all aspects needed to design therapeutic strategies against HCC. In this regard, a computational model provided by Chong et al proposed a GRN model in which they identified three attractors: homeostasis, cell cycle arrest, and senescence. It is of interest that the senescence phenotype was favored, as occurs experimentally (Chong et al., 2018). This model showed that p53 inactivation, in presence of DNA damage, modifies the epigenetic landscape structure and allows the emergence of a new attractor related to a cancer state (Chong et al., 2018). In the same line, Méndez-López and coworkers proposed a model of EMT in which they identified attractors that correspond to epithelial, senescent epithelial, and stem-like mesenchymal phenotypes. In this model, the loss of function of p16 abrogates the senescent epithelium attractor. Furthermore, the results suggested that senescent epithelial cells might exit the senescent state and acquire mesenchymal stem cell characteristics, as has been observed in spontaneous immortalization (Méndez-López et al., 2017).
Concerning to the immune system, computational models have been used to improve detection of tumors using IgG and IgM antibodies (Díaz-Zaragoza et al., 2015). Also, models of multi-scale agents have been created to simulate the interactions between CD4+ T cells and macrophages in the tumor microenvironment (Cess and Finley, 2020). In the same direction, models based on Boolean GRNs have been made to understand the fine details of macrophage differentiation and formation of TAMs in the context of the tumor microenvironment (Avila-Ponce de León et al., 2021). These models have shown that the phenotypic diversity of TAMs can be explained by a sensitive balance produced by transcriptional factors such as NF-κB, STAT3, STAT6, AP-1 and HIF-1α (Avila-Ponce de León et al., 2021). However, none of the aforementioned models has integrated by itself the complexity of cancer, in particular specialized aspects such as those seen in HCC.
Discussion
Currently, there is no single way to design therapeutic strategies against cancer, since all approaches, to the extent of their possibilities, may contribute against this disease. However, it is tempting to assume that it is possible to create rational therapeutic strategies. Especially if we consider the difference in mortality associated with different types of cancer, such as HCC (Rumgay et al., 2022). Furthermore, it would be extremely useful to be able to design personalized therapeutic strategies at low cost for each patient. In this direction, engineering shows that it is possible to create therapeutic strategies as long as two conditions are met: first, the behavior of the system to control can be visualized and second, the system must be controllable (Mellodge, 2015). Regarding the first condition, in the present work we have found ample experimental evidence that demonstrates that there is enough data to create models capable of visualizing the complexity of cancer, including specific types of cancer such as HCC. Concerning the second condition, the experimental evidence cited in this work proves that it is possible to guide cancer cells towards diverse phenotypes. In other words, it is theoretically possible to create therapeutic strategies using rational engineering principles to control cancer.
Consequently, a reasonable question is why this approach is not currently applied to the clinic. Although hypothetically this approach can be used, the truth is that there are several technical limitations that make its clinical implementation difficult. Firstly, even though there is sufficient information to model the complexity of cancer, the use of GRNs models requires a large number of biological components, which implies a high computational cost. Secondly, it is necessary to know all the possible side effects of the drugs used and genetic circuits in the design of therapeutic strategies. Regarding the first limitation, characterizing multistability requires excessive computational effort. To illustrate this point, suppose that a network of 30 nodes is being analyzed, this means that the number of states that the network can originate is
Nevertheless, there are advances in computer science that could solve the state space size problem. In particular, the methodology known as Process Hitting appears to be a viable solution to the problem, since it is capable of obtaining information from immensely large networks without being limited by the problem of the size of the state space (Folschette et al., 2015; Sheikh et al., 2019). Although this technique is still being experimented, it is not ruled out as a viable solution in the medium-term (Sheikh et al., 2019). Regarding the side effects of drugs, at the moment it is impossible to predict their behavior before clinically testing. However, databases of side effects such as SIDER (Kuhn et al., 2016) could help to select drugs to be considered in the design of therapeutic strategies.
In conclusion, from everything discussed up to this point, it is clear that the rational design of therapeutic strategies against liver cancer, as well as for other types of cancer, could be achieved with sustained efforts in the medium-term. As a fundamental condition for this, it is necessary to deal with the problem of the state space of the GRNs and in addition, the built models must necessarily consider the interactions of cancer cells with their microenvironment. This implies including models of the immune cells that constitute the periphery of the tumors, as well as the cancer cells in their different phenotypes inside the tumors.
Author contributions
AH: Conceptualization, Investigation, Methodology, Writing–original draft, Writing–review and editing. AB: Methodology, Visualization, Writing–original draft, Writing–review and editing. JM-G: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing–original draft, Writing–review and editing. EA-B: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. EA-B and JM-G acknowledge the support from UNAM-DGAPA PAPIIT IN211721 “Patrones genéricos y sistémicos de la diferenciación y la proliferación en los nichos de células troncales: Raíz de Arabidopsis thaliana como sistema de estudio teórico-experimental” and CONACYT-FORDECYT-PRONACES 194186/2020 “Biología matemática y computacional de sistemas médicos: modulación preventiva de la emergencia y progresión de enfermedades crónico-degenerativas,” respectively.
Acknowledgments
AH thanks to CONAHCYT for his doctoral fellowship. AB thanks to CONACYT and Cinvestav-IPN for his designation as “Investigador por México” and for the research support, respectively.
Conflict of interest
The authors declare 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
Ajo-Franklin, C. M., Drubin, D. A., Eskin, J. A., Gee, E. P. S., Landgraf, D., Phillips, I., et al. (2007). Rational design of memory in eukaryotic cells. Genes. Dev. 21 (18), 2271–2276. doi:10.1101/gad.1586107
Álvarez-Buylla Roces, M. E., Martínez-García, J. C., Dávila-Velderrain, J., Domínguez-Hüttinger, E., and Martínez-Sánchez, M. E. (2018) Modeling methods for medical systems biology, 1069, 258.
Andriani, F., Bertolini, G., Facchinetti, F., Baldoli, E., Moro, M., Casalini, P., et al. (2016). Conversion to stem-cell state in response to microenvironmental cues is regulated by balance between epithelial and mesenchymal features in lung cancer cells. Mol. Oncol. 10 (2), 253–271. doi:10.1016/j.molonc.2015.10.002
Arvanitakis, K., Koletsa, T., Mitroulis, I., and Germanidis, G. (2022). Tumor-associated macrophages in hepatocellular carcinoma pathogenesis, prognosis and therapy. Cancers 14, 226. doi:10.3390/cancers14010226
Arvanitakis, K., Mitroulis, I., and Germanidis, G. (2021). Tumor-associated neutrophils in hepatocellular carcinoma pathogenesis, prognosis, and therapy. Cancers 13, 2899–2917. doi:10.3390/cancers13122899
Avila-Ponce de León, U., Vázquez-Jiménez, A., Matadamas-Guzman, M., Pelayo, R., and Resendis-Antonio, O. (2021). Transcriptional and microenvironmental landscape of macrophage transition in cancer: a boolean analysis. Front. Immunol. 12 (642842), 1–17. doi:10.3389/fimmu.2021.642842
Azizi, E., Carr, A. J., Plitas, G., Cornish, A. E., Konopacki, C., Prabhakaran, S., et al. (2018). Single-cell map of diverse immune phenotypes in the breast tumor microenvironment. Cell. 174 (5), 1293–1308. doi:10.1016/j.cell.2018.05.060
Berasain, C., and Avila, M. A. (2015). Regulation of hepatocyte identity and quiescence. Cell. Mol. Life Sci. 72 (20), 3831–3851. doi:10.1007/s00018-015-1970-7
Bommer, G. T., Gerin, I., Feng, Y., Kaczorowski, A. J., Kuick, R., Love, R. E., et al. (2007). p53-mediated activation of miRNA34 candidate tumor-suppressor genes. Curr. Biol. 17, 1298–1307. doi:10.1016/j.cub.2007.06.068
Boumahdi, S., and de Sauvage, F. J. (2020). The great escape: tumour cell plasticity in resistance to targeted therapy. Nat. Rev. Drug Discov. 19 (1), 39–56. doi:10.1038/s41573-019-0044-1
Buldum, G., Tsipa, A., and Mantalaris, A. (2020). Linking engineered gene circuit kinetic modeling to cellulose biosynthesis prediction in Escherichia coli: toward bioprocessing of microbial cell factories. Ind. Eng. Chem. Res. 59 (10), 4659–4669. doi:10.1021/acs.iecr.9b05847
Cai, W. Y., Lin, L. Y., Hao, H., Zhang, S. M., Ma, F., Hong, X. X., et al. (2017). Yes-associated protein/TEA domain family member and hepatocyte nuclear factor 4-alpha (HNF4α) repress reciprocally to regulate hepatocarcinogenesis in rats and mice. Hepatology 65 (4), 1206–1221. doi:10.1002/hep.28911
Carné Trécesson, S. De, Guillemin, Y., Bélanger, A., Bernard, A. C., Preisser, L., Ravon, E., et al. (2011). Escape from p21-mediated oncogene-induced senescence leads to cell dedifferentiation and dependence on anti-apoptotic Bcl-xL and MCL1 proteins. J. Biol. Chem. 286 (15), 12825–12838. doi:10.1074/jbc.M110.186437
Cess, C. G., and Finley, S. D. (2020). Multi-scale modeling of macrophage—T cell interactions within the tumor microenvironment. PLoS Comput. Biol. 16 (12), e1008519–e1008535. doi:10.1371/journal.pcbi.1008519
Chaffer, C. L., Brueckmann, I., Scheel, C., Kaestli, A. J., Wiggins, P. A., Rodrigues, L. O., et al. (2011). Normal and neoplastic nonstem cells can spontaneously convert to a stem-like state. Proc. Natl. Acad. Sci. 108 (19), 7950–7955. doi:10.1073/pnas.1102454108
Chang, C. J., Chao, C. H., Xia, W., Yang, J. Y., Xiong, Y., Li, C. W., et al. (2011). P53 regulates epithelial-mesenchymal transition and stem cell properties through modulating miRNAs. Nat. Cell. Biol. 13 (3), 317–323. doi:10.1038/ncb2173
Chen, T., You, Y., Jiang, H., and Wang, Z. Z. (2017). Epithelial–mesenchymal transition (EMT): a biological process in the development, stem cell differentiation, and tumorigenesis. J. Cell. Physiol. 232 (12), 3261–3272. doi:10.1002/jcp.25797
Choi, K. R., Jang, W. D., Yang, D., Cho, J. S., Park, D., and Lee, S. Y. (2019). Systems metabolic engineering strategies: integrating systems and synthetic biology with metabolic engineering. Trends Biotechnol. 37 (8), 817–837. doi:10.1016/j.tibtech.2019.01.003
Chong, K. H., Zhang, X., and Zheng, J. (2018). Dynamical analysis of cellular ageing by modeling of gene regulatory network based attractor landscape. PLoS One 13 (6), 01978388–e197921. doi:10.1371/journal.pone.0197838
Cicchini, C., Amicone, L., Alonzi, T., Marchetti, A., Mancone, C., and Tripodi, M. (2015). Molecular mechanisms controlling the phenotype and the EMT/MET dynamics of hepatocyte. Liver Int. 35 (2), 302–310. doi:10.1111/liv.12577
Cui, C., Wong, C. C., Kai, A. K., Ho, D. W., Lau, E. Y., Tsui, Y., et al. (2017) SENP1 promotes hypoxia-induced cancer stemness by HIF-1 α deSUMOylation and SENP1/HIF-1 α positive feedback loop, 2149–2159.
Czauderna, C., Castven, D., Mahn, F. L., and Marquardt, J. U. (2019). Context-dependent role of NF-κB signaling in primary liver cancer—from tumor development to therapeutic implications. Cancers (Basel) 11 (8), 1053. doi:10.3390/cancers11081053
Daoust, R., and Cantero, A. (1959). The numerical proportions of cell types in rat liver during carcinogenesis by 4-dimethylaminoazobenzene (DAB). Cancer Res. 19 (7), 757–762.
Davila-velderrain, J., Martinez-garcia, J. C., and Alvarez-buylla, E. R. (2015). Modeling the epigenetic attractors landscape: toward a post-genomic mechanistic understanding of development. Front. Genet. 6 (April), 160–214. doi:10.3389/fgene.2015.00160
Deng, X., Zhang, X., Li, W., Feng, R.-X., Li, L., Yi, G.-R., et al. (2018). Chronic liver injury induces conversion of biliary epithelial cells into hepatocytes. Cell. Stem Cell. 23 (1), 114–122. doi:10.1016/j.stem.2018.05.022
Díaz-Zaragoza, M., Hernández-Ávila, R., Govezensky, T., Mendoza, L., Meneses-Ruíz, D. M., and Ostoa-Saloma, P. (2015). Comparison patterns of 4 T1 antigens recognized by humoral immune response mediated by IgG and IgM antibodies in female and male mice with breast cancer using 2D-immnunoblots. Immunobiology 220 (9), 1050–1058. doi:10.1016/j.imbio.2015.05.012
Di Blasi, D., Boldanova, T., Mori, L., Terracciano, L., Heim, M. H., and De Libero, G. (2020). Unique T-cell populations define immune-inflamed hepatocellular carcinoma. CMGH 9 (2), 195–218. doi:10.1016/j.jcmgh.2019.08.004
Esteban, M. A., Bao, X., Zhuang, Q., Zhou, T., Qin, B., and Pei, D. (2012). The mesenchymal-to-epithelial transition in somatic cell reprogramming. Curr. Opin. Genet. Dev. 22 (5), 423–428. doi:10.1016/j.gde.2012.09.004
Feng, Y., Ye, Z., Song, F., He, Y., and Liu, J. (2022). The role of TAMs in tumor microenvironment and new research progress. Stem Cells Int. 2022 (5775696), 5775696–5775711. doi:10.1155/2022/5775696
Ferlay, J., Colombet, M., Soerjomataram, I., Parkin, D. M., Piñeros, M., Znaor, A., et al. (2021). Cancer statistics for the year 2020: an overview. Int. J. Cancer. 149 (4), 1–12. doi:10.1002/ijc.33588
Folschette, M., Paulevé, L., Inoue, K., Magnin, M., and Roux, O. (2015). Identification of biological regulatory networks from Process Hitting models. Theor. Comput. Sci. 568 (C), 49–71. doi:10.1016/j.tcs.2014.12.002
Forte, E., Chimenti, I., Rosa, P., Angelini, F., Pagano, F., Calogero, A., et al. (2017). EMT/MET at the crossroad of stemness, regeneration and oncogenesis: the Ying-Yang equilibrium recapitulated in cell spheroids. Cancers (Basel) 9 (8), 98–15. doi:10.3390/cancers9080098
Galvagni, F., Lentucci, C., Neri, F., Dettori, D., De Clemente, C., Orlandini, M., et al. (2015). Snai1 promotes ESC exit from the pluripotency by direct repression of self-renewal genes. Stem Cells 33 (3), 742–750. doi:10.1002/stem.1898
Garibaldi, F., Cicchini, C., Conigliaro, A., Santangelo, L., Cozzolino, A. M., Grassi, G., et al. (2012). An epistatic mini-circuitry between the transcription factors Snail and HNF4α controls liver stem cell and hepatocyte features exhorting opposite regulation on stemness-inhibiting microRNAs. Cell. Death Differ. 19 (6), 937–946. doi:10.1038/cdd.2011.175
Giannelli, G., Koudelkova, P., Dituri, F., and Mikulits, W. (2016). Role of epithelial to mesenchymal transition in hepatocellular carcinoma. J. Hepatol. 65 (4), 798–808. doi:10.1016/j.jhep.2016.05.007
Gifford, G. B., Demitrack, E. S., Keeley, T. M., Tam, A., La Cunza, N., Dedhia, P. H., et al. (2017). Notch1 and Notch2 receptors regulate mouse and human gastric antral epithelial cell homoeostasis. Gut 66 (6), 1001–1011. doi:10.1136/gutjnl-2015-310811
Godoy, P., Hengstler, J. G., Ilkavets, I., Meyer, C., Bachmann, A., Müller, A., et al. (2009). Extracellular matrix modulates sensitivity of hepatocytes to fibroblastoid dedifferentiation and transforming growth factor beta-induced apoptosis. Hepatology 49 (6), 2031–2043. doi:10.1002/hep.22880
Godoy, P., Widera, A., Schmidt-Heck, W., Campos, G., Meyer, C., Cadenas, C., et al. (2016). Gene network activity in cultivated primary hepatocytes is highly similar to diseased mammalian liver tissue. Arch. Toxicol. 90 (10), 2513–2529. doi:10.1007/s00204-016-1761-4
Haellman, V., Saxena, P., Jiang, Y., and Fussenegger, M. (2021). Rational design and optimization of synthetic gene switches for controlling cell-fate decisions in pluripotent stem cells. Metab. Eng. 65, 99–110. doi:10.1016/j.ymben.2021.03.009
Han, X., Yang, J., Zeng, F., Weng, J., Zhang, Y., Peng, Q., et al. (2020). Programmable synthetic protein circuits for the identification and suppression of hepatocellular carcinoma. Mol. Ther. Oncolytics 17, 70–82. doi:10.1016/j.omto.2020.03.008
Hari, K., Ullanat, V., Balasubramanian, A., Gopalan, A., and Jolly, M. K. (2022). Landscape of Epithelial Mesenchymal Plasticity as an emergent property of coordinated teams in regulatory networks. Elife 11, 765355–e76631. doi:10.7554/eLife.76535
Harner-Foreman, N., Vadakekolathu, J., Laversin, S. A., Mathieu, M. G., Reeder, S., Pockley, A. G., et al. (2017). A novel spontaneous model of epithelial-mesenchymal transition (EMT) using a primary prostate cancer derived cell line demonstrating distinct stem-like characteristics. Sci. Rep. 7 (December 2016), 40633–40718. doi:10.1038/srep40633
He, F., Murabito, E., and Westerhoff, H. V. (2016). Synthetic biology and regulatory networks: where metabolic systems biology meets control engineering. J. R. Soc. Interface 13, 20151046–20151119. doi:10.1098/rsif.2015.1046
Herranz, N., Gil, J., Herranz, N., and Gil, J. (2018). Mechanisms and functions of cellular senescence Find the latest version: mechanisms and functions of cellular senescence. JCI J. Clin. Investig. 128 (4), 1238–1246. doi:10.1172/JCI95148
Higashi, T., Friedman, S. L., and Hoshida, Y. (2017). Hepatic stellate cells as key target in liver fibrosis. Adv. Drug Deliv. Rev. 121, 27–42. doi:10.1016/j.addr.2017.05.007
Hikichi, T., Matoba, R., Ikeda, T., Watanabe, A., Yamamoto, T., Yoshitake, S., et al. (2013). Transcription factors interfering with dedifferentiation induce cell type-specific transcriptional profiles. Proc. Natl. Acad. Sci. U. S. A. 110 (16), 6412–6417. doi:10.1073/pnas.1220200110
Huang, B., Lu, M., Jia, D., Ben-Jacob, E., Levine, H., and Onuchic, J. N. (2017). Interrogating the topological robustness of gene regulatory circuits by randomization. PLoS Comput. Biol. 13 (3), 10054566–e1005521. doi:10.1371/journal.pcbi.1005456
Huang, S. (2013). Genetic and non-genetic instability in tumor progression: link between the fitness landscape and the epigenetic landscape of cancer cells. Cancer Metastasis Rev. 32 (3–4), 423–448. doi:10.1007/s10555-013-9435-7
Huang, S. (2020). Reconciling non-genetic plasticity with somatic evolution in cancer. Trends Cancer 7 (4), 309–322. doi:10.1016/j.trecan.2020.12.007
Huang, Y., Hong, W., and Wei, X. (2022). The molecular mechanisms and therapeutic strategies of EMT in tumor progression and metastasis. J. Hematol. Oncol. 15, 129–227. doi:10.1186/s13045-022-01347-8
Ichise, H., Tsukamoto, S., Hirashima, T., Konishi, Y., Oki, C., Tsukiji, S., et al. Functional visualization of NK cell-mediated killing of metastatic single tumor cells. Elife. 2022;11:1–28. doi:10.7554/eLife.76269
Jayachandran, A., Dhungel, B., and Steel, J. C. (2016). Epithelial-to-mesenchymal plasticity of cancer stem cells: therapeutic targets in hepatocellular carcinoma. J. Hematol. Oncol. 9, 74–12. doi:10.1186/s13045-016-0307-9
Jia, D., Jolly, M. K., Kulkarni, P., and Levine, H. (2017). Phenotypic plasticity and cell fate decisions in cancer: insights from dynamical systems theory. Cancers (Basel) 9 (7), 70–19. doi:10.3390/cancers9070070
Jing, L., Sun, H., Wang, Y., Han, L., and Jiao, M. (2019) Epithelial – mesenchymal transition induced cancer - stem - cell - like characteristics in hepatocellular carcinoma, 1–11.
Jolly, M. K., Huang, B., Lu, M., Mani, S. A., Levine, H., Ben-jacob, E., et al. (2014) Towards elucidating the connection between epithelial – mesenchymal transitions and stemness.
Jong, I. Y., Li, X. Y., Ota, I., Fearon, E. R., and Weiss, S. J. (2005). Wnt-dependent regulation of the E-cadherin repressor snail. J. Biol. Chem. 280 (12), 11740–11748. doi:10.1074/jbc.M413878200
Jopling, C., Boue, S., and Belmonte, J. C. I. (2011). Dedifferentiation, transdifferentiation and reprogramming: three routes to regeneration. Nat. Rev. Mol. Cell. Biol. 12 (2), 79–89. doi:10.1038/nrm3043
Karlebach, G., and Shamir, R. (2008). Modelling and analysis of gene regulatory networks. Nat. Rev. Mol. Cell. Biol. 9 (10), 770–780. doi:10.1038/nrm2503
Kim, N. H., Kim, H. S., Li, X., Lee, I., Choi, H., Kang, S. E., et al. (2011). A p53/miRNA-34 axis regulates Snail1-dependent cancer cell epithelial–mesenchymal transition. J. Cell. Biol. 195 (3), 417–433. doi:10.1083/jcb.201103097
Kishi, S., Bayliss, P. E., and Hanai, J. I. (2015). A prospective epigenetic paradigm between cellular senescence and epithelial-mesenchymal transition in organismal development and aging. Transl. Res. 165 (1), 241–249. doi:10.1016/j.trsl.2014.05.007
Ko, S., Russell, J. O., Molina, L. M., and Monga, S. P. (2020) Liver progenitors and adult cell plasticity in hepatic injury and repair: knowns and unknowns.
Kordes, C., Bock, H. H., Reichert, D., May, P., and Häussinger, D. (2021). Hepatic stellate cells: current state and open questions. Biol. Chem. 402 (9), 1021–1032. doi:10.1515/hsz-2021-0180
Kordes, C., Sawitza, I., Schumacher, E., Diseases, I., and Häussinger, D. (2015). Beyond fibrosis: stellate cells as liver stem cells. Z. Gastroenterol. 53, 1425–1431. doi:10.1055/s-0035-1566895
Kuhn, M., Letunic, I., Jensen, L. J., and Bork, P. (2016). The SIDER database of drugs and side effects. Nucleic Acids Res. 44 (D1), D1075–D1079. doi:10.1093/nar/gkv1075
Ladewig, J., Koch, P., and Brüstle, O. (2013). Leveling Waddington: the emergence of direct programming and the loss of cell fate hierarchies. Nat. Publ. Gr 14, 225–236. doi:10.1038/nrm3543
Lapasset, L., Milhavet, O., Prieur, A., Besnard, E., Babled, A., Ät-Hamou, N., et al. (2011). Rejuvenating senescent and centenarian human cells by reprogramming through the pluripotent state. Genes. Dev. 25 (21), 2248–2253. doi:10.1101/gad.173922.111
Lau, H. H., Ng, N. H. J., Loo, L. S. W., Jasmen, J. B., and Teo, A. K. K. (2018). The molecular functions of hepatocyte nuclear factors – in and beyond the liver. J. Hepatol. 68 (5), 1033–1048. doi:10.1016/j.jhep.2017.11.026
Laurent, M., and Kellershohn, N. (1999). Multistability: a major means of differentiation and evolution in biological systems. Trends Biochem. Sci. 24 (11), 418–422. doi:10.1016/s0968-0004(99)01473-5
Lee, S., and Schmitt, C. A. (2019). The dynamic nature of senescence in cancer. Nat. Cell. Biol. 21 (January), 94–101. doi:10.1038/s41556-018-0249-2
Lee, S. H., Lee, S. J., Yeon, S. J., Xu, Y., Ho, S. K., Ha, N. C., et al. (2009). Blocking of p53-snail binding, promoted by oncogenic K-Ras, recovers p53 expression and function. Neoplasia 11 (1), 22–31. doi:10.1593/neo.81006
Li, W., Li, L., and Hui, L. (2020). Cell plasticity in liver regeneration. Trends Cell. Biol. 30 (4), 329–338. doi:10.1016/j.tcb.2020.01.007
Lim, S., Kim, H., and Jung, G. (2010). p53 inhibits tumor cell invasion via the degradation of snail protein in hepatocellular carcinoma. FEBS Lett. 584 (11), 2231–2236. doi:10.1016/j.febslet.2010.04.006
Lin, M. W., Tseng, Y. W., Shen, C. C., Hsu, M. N., Hwu, J. R., Chang, C. W., et al. (2018). Synthetic switch-based baculovirus for transgene expression control and selective killing of hepatocellular carcinoma cells. Nucleic Acids Res. 46 (15), e93–e15. doi:10.1093/nar/gky447
López-Nouoa, J. M., and Nieto, M. A. (2009). Inflammation and EMT: an alliance towards organ fibrosis and cancer progression. EMBO Mol. Med. 1 (6–7), 303–314. doi:10.1002/emmm.200900043
Loyer, P., Ilyin, G., Cariou, S., Glaise, D., Corlu, A., and Guguen-Guillouzo, C. (1996). Progression through G1 and S phases of adult rat hepatocytes. Prog. Cell. Cycle Res. 2, 37–47. doi:10.1007/978-1-4615-5873-6_4
Ma, J., Zheng, B., Goswami, S., Meng, L., Zhang, D., Cao, C., et al. (2019). PD1Hi CD8+ T cells correlate with exhausted signature and poor clinical outcome in hepatocellular carcinoma. J. Immunother. Cancer 7 (1), 331–415. doi:10.1186/s40425-019-0814-7
Mani, S. A., Guo, W., Liao, M. J., Eaton, E. N., Ayyanan, A., Zhou, A. Y., et al. (2008). The epithelial-mesenchymal transition generates cells with properties of stem cells. Cell. 133 (4), 704–715. doi:10.1016/j.cell.2008.03.027
Marusyk, A., Janiszewska, M., and Polyak, K. (2020). Intratumor heterogeneity: the rosetta stone of therapy resistance. Cancer Cell. 37 (4), 471–484. doi:10.1016/j.ccell.2020.03.007
Marusyk, A., and Polyak, K. (2010). Tumor heterogeneity: causes and consequences. Biochim. Biophys. Acta 1805, 105–117. doi:10.1016/j.bbcan.2009.11.002
Masucci, M. T., Minopoli, M., Del Vecchio, S., and Carriero, M. V. (2020). The emerging role of neutrophil extracellular Traps (NETs) in tumor progression and metastasis. Front. Immunol. 11, 1749–1757. doi:10.3389/fimmu.2020.01749
May, C. D., Sphyris, N., Evans, K. W., Werden, S. J., Guo, W., and Mani, S. A. (2011). Epithelial-mesenchymal transition and cancer stem cells: a dangerously dynamic duo in breast cancer progression. Breast Cancer Res. 13 (1), 202–210. doi:10.1186/bcr2789
Mellodge, P. (2015). A practical approach to dynamical systems for engineers. A Pract. Approach Dyn. Syst. Eng.
Méndez-López, L. F., Davila-Velderrain, J., Domínguez-Hüttinger, E., Enríquez-Olguín, C., Martínez-García, J. C., and Alvarez-Buylla, E. R. (2017). Gene regulatory network underlying the immortalization of epithelial cells. BMC Syst. Biol. 11 (1), 24. doi:10.1186/s12918-017-0393-5
Michalopoulos, G. K., and Bhushan, B. (2021). Liver regeneration: biological and pathological mechanisms and implications. Nat. Rev. Gastroenterol. Hepatol. 18 (1), 40–55. doi:10.1038/s41575-020-0342-4
Migault, M., Sapkota, S., and Bracken, C. P. (2022). Transcriptional and post-transcriptional control of epithelial-mesenchymal plasticity: why so many regulators? Cell. Mol. Life Sci. 79 (3), 182. doi:10.1007/s00018-022-04199-0
Milanovic, M., Fan, D. N. Y., Belenki, D., Däbritz, J. H. M., Zhao, Z., Yu, Y., et al. (2018). Senescence-associated reprogramming promotes cancer stemness. Nature 553 (7686), 96–100. doi:10.1038/nature25167
Min, C., Eddy, S. F., Sherr, D. H., and Sonenshein, G. E. (2008). NF- k B and epithelial to mesenchymal transition of. Cancer 744, 733–744. doi:10.1002/jcb.21695
Mohan, A., Raj Rajan, R., Mohan, G., Kollenchery Puthenveettil, P., and Maliekal, T. T. (2021). Markers and reporters to reveal the hierarchy in heterogeneous cancer stem cells. Front. Cell. Dev. Biol. 9 (June), 668851. doi:10.3389/fcell.2021.668851
Morel, A. P., Lièvre, M., Thomas, C., Hinkal, G., Ansieau, S., and Puisieux, A. (2008). Generation of breast cancer stem cells through epithelial-mesenchymal transition. PLoS One 3 (8), 28888–e2897. doi:10.1371/journal.pone.0002888
Ni, T., Li, X. Y., Lu, N., An, T., Liu, Z. P., Fu, R., et al. (2016). Snail1-dependent p53 repression regulates expansion and activity of tumour-initiating cells in breast cancer. Nat. Cell. Biol. 18 (11), 1221–1232. doi:10.1038/ncb3425
Ning, B., Ding, J., Liu, J., Yin, C., Xu, W., Cong, W., et al. (2014) Feedback circuit modulates liver cancer progression.
Nio, K., Yamashita, T., and Kaneko, S. (2017). The evolving concept of liver cancer stem cells. Mol. Cancer 16 (1), 4–12. doi:10.1186/s12943-016-0572-9
Noce, V., Battistelli, C., Cozzolino, A. M., Consalvi, V., Cicchini, C., Strippoli, R., et al. (2019). YAP integrates the regulatory Snail/HNF4α circuitry controlling epithelial/hepatocyte differentiation. Cell. Death Dis. 10 (10), 768. doi:10.1038/s41419-019-2000-8
Ogata, K., and Brewer, J. W. (1971). Modern control engineering. J. Dyn. Syst. Meas. Control 93, 63. doi:10.1115/1.3426465
Oh, S. H., Swiderska-Syn, M., Jewell, M. L., Premont, R. T., and Diehl, A. M. (2018). Liver regeneration requires Yap1-TGFβ-dependent epithelial-mesenchymal transition in hepatocytes. J. Hepatol. 69 (2), 359–367. doi:10.1016/j.jhep.2018.05.008
Paez Ribes, M., González-Gualda, E., Doherty, G. J., and Muñoz-Espín, D. (2019). Targeting senescent cells in translational medicine. EMBO Mol. Med. 11 (12), 102344–e10319. doi:10.15252/emmm.201810234
Pan, Y., Yu, Y., Wang, X., and Zhang, T. (2020). Tumor-associated macrophages in tumor immunity. Front. Immunol. 11, 1–9. doi:10.3389/fimmu.2020.583084
Pandey, A., Singh, P., Jauhari, A., Singh, T., Khan, F., Pant, A. B., et al. (2015). Critical role of the miR-200 family in regulating differentiation and proliferation of neurons. J. Neurochem. 133 (5), 640–652. doi:10.1111/jnc.13089
Patel, P. L., Suram, A., Mirani, N., Bischof, O., and Herbig, U. (2016). Derepression of hTERT gene expression promotes escape from oncogene-induced cellular senescence. Proc. Natl. Acad. Sci. U. S. A. 113 (34), E5024–E5033. doi:10.1073/pnas.1602379113
Pisco, A. O., and Huang, S. (2015). Non-genetic cancer cell plasticity and therapy-induced stemness in tumour relapse: what does not kill me strengthens me. Br. J. Cancer 112, 1725–1732. doi:10.1038/bjc.2015.146
Qiu, S. Q., Waaijer, S. J. H., Zwager, M. C., de Vries, E. G. E., van der Vegt, B., and Schröder, C. P. (2018). Tumor-associated macrophages in breast cancer: innocent bystander or important player? Cancer Treat. Rev. 70, 178–189. doi:10.1016/j.ctrv.2018.08.010
Rajagopal, J., and Stanger, B. Z. (2016a). Perspective plasticity in the adult: how should the Waddington diagram Be applied to regenerating tissues. Dev. Cell. 36 (2), 133–137. doi:10.1016/j.devcel.2015.12.021
Rajagopal, J., and Stanger, B. Z. (2016b). Perspective plasticity in the adult: how should the Waddington diagram Be applied to regenerating tissues ? Dev. Cell. 36, 133–137. doi:10.1016/j.devcel.2015.12.021
Ran, J., Li, Y., Liu, L., Zhu, Y., Ni, Y., Huang, H., et al. (2021). Apelin enhances biological functions in lung cancer A549 cells by downregulating exosomal miR-15a-5p. Carcinogenesis 42 (2), 243–253. doi:10.1093/carcin/bgaa089
Roskams, T. (2006). Different types of liver progenitor cells and their niches. J. Hepatol. 45 (1), 1–4. doi:10.1016/j.jhep.2006.05.002
Ruijtenberg, S., and van den Heuvel, S. (2016). Coordinating cell proliferation and differentiation: antagonism between cell cycle regulators and cell type-specific gene expression. Cell. Cycle 15 (2), 196–212. doi:10.1080/15384101.2015.1120925
Rumgay, H., Arnold, M., Ferlay, J., Lesi, O., Cabasag, C. J., Vignat, J., et al. (2022). Global burden of primary liver cancer in 2020 and predictions to 2040. J. Hepatol. 77 (6), 1598–1606. doi:10.1016/j.jhep.2022.08.021
Sajid, M., Liu, L., and Sun, C. (2022). The dynamic role of NK cells in liver cancers: role in HCC and HBV associated HCC and its therapeutic implications. Front. Immunol. 13, 887186–887213. doi:10.3389/fimmu.2022.887186
Saleh, T., Tyutyunyk-Massey, L., Murray, G. F., Alotaibi, M. R., Kawale, A. S., Elsayed, Z., et al. (2019). Tumor cell escape from therapy-induced senescence. Biochem. Pharmacol. 162 (December 2018), 202–212. doi:10.1016/j.bcp.2018.12.013
Saunders, L. R., and McClay, D. R. (2014). Sub-circuits of a gene regulatory network control a developmental epithelial-mesenchymal transition. Dev 141 (7), 1503–1513. doi:10.1242/dev.101436
Shaul, M. E., Levy, L., Sun, J., Mishalian, I., Singhal, S., Kapoor, V., et al. (2016). Tumor-associated neutrophils display a distinct N1 profile following TGFβ modulation: a transcriptomics analysis of pro-vs. antitumor TANs. Oncoimmunology 5 (11), 12322211–e1232314. doi:10.1080/2162402X.2016.1232221
Sheikh, I. A., Ahmad, J., Magnin, M., and Roux, O. (2019). Incorporating time delays in process hitting framework for dynamical modeling of large biological regulatory networks. Front. Physiol. 10 (FEB), 90–17. doi:10.3389/fphys.2019.00090
Shi, J., Teschendorff, A. E., Chen, W., Chen, L., and Li, T. (2018). Quantifying Waddington’s epigenetic landscape: a comparison of single-cell potency measures. Brief. Bioinform 21 (1), 248–261. doi:10.1093/bib/bby093
Shi, Z. D., Pang, K., Wu, Z. X., Dong, Y., Hao, L., Qin, J. X., et al. (2023). Tumor cell plasticity in targeted therapy-induced resistance: mechanisms and new strategies. Signal Transduct. Target Ther. 8 (1), 113–121. doi:10.1038/s41392-023-01383-x
Shih, R. M., and Chen, Y. Y. (2022). Engineering principles for synthetic biology circuits in cancer immunotherapy. Cancer Immunol. Res. 10 (1), 6–11. doi:10.1158/2326-6066.CIR-21-0769
Shu, X., and Pei, D. (2014). The function and regulation of mesenchymal-to-epithelial transition in somatic cell reprogramming. Curr. Opin. Genet. Dev. 28, 32–37. doi:10.1016/j.gde.2014.08.005
Siuti, P., Yazbek, J., and Lu, T. K. (2014). Engineering genetic circuits that compute and remember. Nat. Protoc. 9, 1292–1300. doi:10.1038/nprot.2014.089
Smith, A. P., Verrecchia, A., Fagà, G., Doni, M., Perna, D., Martinato, F., et al. (2009). A positive role for Myc in TGFbeta-induced Snail transcription and epithelial-to-mesenchymal transition. Oncogene 28 (3), 422–430. doi:10.1038/onc.2008.395
So, J., Kim, A., Lee, S. H., and Shin, D. (2020). Liver progenitor cell-driven liver regeneration. Exp. Mol. Med. 52 (8), 1230–1238. doi:10.1038/s12276-020-0483-0
Takasugi, M., Yoshida, Y., Hara, E., and Ohtani, N. (2022). The role of cellular senescence and SASP in tumour microenvironment. FEBS J. 290, 1348–1361. doi:10.1111/febs.16381
Taniguchi, K., and Karin, M. (2018). NF-B, inflammation, immunity and cancer: coming of age. Nat. Rev. Immunol. 18 (5), 309–324. doi:10.1038/nri.2017.142
Walcher, L., Kistenmacher, A. K., Suo, H., Kitte, R., Dluczek, S., Strauß, A., et al. (2020). Cancer stem cells—origins and biomarkers: perspectives for targeted personalized therapies. Front. Immunol. 11 (August), 1280–1333. doi:10.3389/fimmu.2020.01280
Walesky, C., and Apte, U. (2015). Role of hepatocyte nuclear factor 4α (HNF4α) in cell proliferation and cancer. Gene Expr. 16 (3), 101–108. doi:10.3727/105221615X14181438356292
Wan, S., Zhao, E., Kryczek, I., Vatan, L., Sadovskaya, A., Ludema, G., et al. (2014). Tumor-associated macrophages produce interleukin 6 and signal via STAT3 to promote expansion of human hepatocellular carcinoma stem cells. Gastroenterology 147 (6), 1393–1404. doi:10.1053/j.gastro.2014.08.039
Wang, G., Guo, X., Hong, W., Liu, Q., Wei, T., Lu, C., et al. (2013). Critical regulation of miR-200/ZEB2 pathway in Oct4/Sox2-induced mesenchymal-to-epithelial transition and induced pluripotent stem cell generation. Proc. Natl. Acad. Sci. U. S. A. 110 (8), 2858–2863. doi:10.1073/pnas.1212769110
Wellner, U., Schubert, J., Burk, U. C., Schmalhofer, O., Zhu, F., Sonntag, A., et al. (2009). The EMT-activator ZEB1 promotes tumorigenicity by repressing stemness-inhibiting microRNAs. Nat. Cell. Biol. 11 (12), 1487–1495. doi:10.1038/ncb1998
Wickramaratne, N., Li, R., Tian, T., Khoraki, J., Kang, H. S., Chmielewski, C., et al. (2021). Cholangiocyte Epithelial to Mesenchymal Transition (EMT) is a potential molecular mechanism driving ischemic cholangiopathy in liver transplantation. PLoS One 16 (7 July), 02469788–e247016. doi:10.1371/journal.pone.0246978
Wu, L., and Zhang, X. H. F. (2020). Tumor-associated neutrophils and macrophages—heterogenous but not chaotic. Front. Immunol. 11, 553967–554018. doi:10.3389/fimmu.2020.553967
Wu, M. R., Jusiak, B., and Lu, T. K. (2019). Engineering advanced cancer therapies with synthetic biology. Nat. Rev. Cancer 19 (4), 187–195. doi:10.1038/s41568-019-0121-0
Wu, Y., Zhang, J., Zhang, X., Zhou, H., Liu, G., and Li, Q. (2020). Cancer stem cells: a potential breakthrough in HCC-targeted therapy. Front. Pharmacol. 11 (March), 198–213. doi:10.3389/fphar.2020.00198
Xie, G., Yao, Q., Liu, Y., Du, S., Liu, A., Guo, Z., et al. (2012). IL-6-induced epithelial-mesenchymal transition promotes the generation of breast cancer stem-like cells analogous to mammosphere cultures. Int. J. Oncol. 40 (4), 1171–1179. doi:10.3892/ijo.2011.1275
Yagci, T., Cetin, M., and Ercin, P. B. (2017). Cancer stem cells in hepatocellular carcinoma. J. Gastrointest. Cancer 48 (3), 241–245. doi:10.1007/s12029-017-9960-7
Yamashita, T., and Wang, X. W. (2013). Cancer stem cells in the development of liver cancer Find the latest version: review series Cancer stem cells in the development of liver cancer. J. Clin. Investig. 123 (5), 1911–1918. doi:10.1172/JCI66024
Yan, L., Xu, F., and Dai, C. L. (2018). Relationship between epithelial-to-mesenchymal transition and the inflammatory microenvironment of hepatocellular carcinoma. J. Exp. Clin. Cancer Res. 37 (1), 203–213. doi:10.1186/s13046-018-0887-z
Yang, J., and Ding, S. (2021). Chimeric RNA-binding protein-based killing switch targeting hepatocellular carcinoma cells. Mol. Ther. Nucleic Acids 25, 683–695. doi:10.1016/j.omtn.2021.08.012
Yang, M., Li, S. N., Anjum, K. M., Gui, L. X., Zhu, S. S., Liu, J., et al. (2013). A double-negative feedback loop between Wnt-β-catenin signaling and HNF4α regulates epithelial-mesenchymal transition in hepatocellular carcinoma. J. Cell. Sci. 126 (24), 5692–5703. doi:10.1242/jcs.135053
Yang, X., Zhu, L., Wang, Y., Li, Y., and Wang, D. (2019). Current advance of therapeutic agents in clinical trials potentially targeting tumor plasticity. Front. Oncol. 9 (September), 887–910. doi:10.3389/fonc.2019.00887
Yao, C., Su, L., Shan, J., Zhu, C., Liu, L., Liu, C., et al. (2016). IGF/STAT3/NANOG/Slug signaling Axis simultaneously controls epithelial-mesenchymal transition and stemness maintenance in colorectal cancer. Stem Cells 34 (4), 820–831. doi:10.1002/stem.2320
Yin, X., Zhang, B., Zheng, S., Gao, D., Qiu, S., Wu, W., et al. Coexpression of gene Oct4 and Nanog initiates stem cell characteristics in hepatocellular carcinoma and promotes epithelial-mesenchymal transition through activation of Stat3/Snail signaling. 2015;1–13.
Yuan, S., Norgard, R. J., and Stanger, B. Z. (2019). Cellular plasticity in cancer. Cancer Discov. 9 (7), 837–851. doi:10.1158/2159-8290.CD-19-0015
Zhang, J., Zhang, Q., Lou, Y., Fu, Q., Chen, Q., Wei, T., et al. (2018). Hypoxia-inducible factor-1α/interleukin-1β signaling enhances hepatoma epithelial-mesenchymal transition through macrophages in a hypoxic-inflammatory microenvironment. Hepatology 67 (5), 1872–1889. doi:10.1002/hep.29681
Zhao, M. L., Rabiee, A., Kovary, K. M., Bahrami-Nejad, Z., Taylor, B., and Teruel, M. N. (2020). Molecular competition in G1 controls when cells simultaneously commit to terminally differentiate and exit the cell cycle. Cell. Rep. 31 (11), 107769. doi:10.1016/j.celrep.2020.107769
Zhou, J., Tang, Z., Gao, S., Li, C., Feng, Y., and Zhou, X. (2020). Tumor-associated macrophages: recent insights and therapies. Front. Oncol. 10, 188–213. doi:10.3389/fonc.2020.00188
Zhou, S. L., Zhou, Z. J., Hu, Z. Q., Huang, X. W., Wang, Z., Chen, E. B., et al. (2016). Tumor-associated neutrophils recruit macrophages and T-regulatory cells to promote progression of hepatocellular carcinoma and resistance to sorafenib. Gastroenterology 150 (7), 1646–1658. doi:10.1053/j.gastro.2016.02.040
Keywords: hepatocellular carcinoma, cancer, rational design, cell plasticity, medical systems biology
Citation: Hernández-Magaña A, Bensussen A, Martínez-García JC and Álvarez-Buylla ER (2024) Engineering principles for rationally design therapeutic strategies against hepatocellular carcinoma. Front. Mol. Biosci. 11:1404319. doi: 10.3389/fmolb.2024.1404319
Received: 20 March 2024; Accepted: 23 May 2024;
Published: 13 June 2024.
Edited by:
Valentina Audrito, Università del Piemonte Orientale, ItalyReviewed by:
Shilpi Minocha, Indian Institute of Technology Delhi, IndiaJinyu Zhang, East Tennessee State University, United States
Copyright © 2024 Hernández-Magaña, Bensussen, Martínez-García and Álvarez-Buylla. 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: Juan Carlos Martínez-García, juancarlos.martinez@cinvestav.mx; Elena R. Álvarez-Buylla, eabuylla@gmail.com