Transitions between epithelial and mesenchymal phenotypes – the epithelial to ­mesenchymal transition (EMT) and its reverse the mesenchymal to epithelial transition (MET) – are hallmarks of cancer metastasis. While transitioning between the epithelial and mesenchymal phenotypes, cells can also attain a hybrid epithelial/mesenchymal (E/M) (i.e., partial or intermediate EMT) phenotype. Cells in this phenotype have mixed epithelial (e.g., adhesion) and mesenchymal (e.g., migration) properties, thereby allowing them to move collectively as clusters. If these clusters reach the bloodstream intact, they can give rise to clusters of circulating tumor cells (CTCs), as have often been seen experimentally. Here, we review the operating principles of the core regulatory network for EMT/MET that acts as a “three-way” switch giving rise to three distinct phenotypes – E, M and hybrid E/M – and present a theoretical framework that can elucidate the role of many other players in regulating epithelial plasticity. Furthermore, we highlight recent studies on partial EMT and its association with drug resistance and tumor-initiating potential; and discuss how cell–cell communication between cells in a partial EMT phenotype can enable the formation of clusters of CTCs. These clusters can be more apoptosis-resistant and have more tumor-initiating potential than singly moving CTCs with a wholly mesenchymal (complete EMT) phenotype. Also, more such clusters can be formed under inflammatory conditions that are often generated by various therapies. Finally, we discuss the multiple advantages that the partial EMT or hybrid E/M phenotype have as compared to a complete EMT phenotype and argue that these collectively migrating cells are the primary “bad actors” of metastasis.

Epithelial-to-Mesenchymal Transition (EMT) and its reverse – Mesenchymal to Epithelial Transition (MET) – are hallmarks of cellular plasticity during embryonic development and cancer metastasis. During EMT, epithelial cells lose cell-cell adhesion and gain migratory and invasive traits either partially or completely, leading to a hybrid epithelial/mesenchymal (hybrid E/M) or a mesenchymal phenotype respectively. Mesenchymal cells move individually, but hybrid E/M cells migrate collectively as observed during gastrulation, wound healing, and the formation of tumor clusters detected as Circulating Tumor Cells (CTCs). Typically, the hybrid E/M phenotype has largely been tacitly assumed to be transient and ‘metastable’. Here, we identify certain ‘phenotypic stability factors’ (PSFs) such as GRHL2 that couple to the core EMT decision-making circuit (miR-200/ZEB) and stabilize hybrid E/M phenotype. Further, we show that H1975 lung cancer cells can display a stable hybrid E/M phenotype and migrate collectively, a behavior that is impaired by knockdown of GRHL2 and another previously identified PSF - OVOL. In addition, our computational model predicts that GRHL2 can also associate hybrid E/M phenotype with high tumor-initiating potential, a prediction strengthened by the observation that the higher levels of these PSFs may be predictive of poor patient outcome. Finally, based on these specific examples, we deduce certain network motifs that can stabilize the hybrid E/M phenotype. Our results suggest that partial EMT, i.e. a hybrid E/M phenotype, need not be ‘metastable’, and strengthen the emerging notion that partial EMT, but not necessarily a complete EMT, is associated with aggressive tumor progression.

Tumor budding is a histological phenomenon encountered in various cancers, whereby individual malignant cells and/or small clusters of malignant cells are seen in the tumor stroma. Postulated to be mirror epithelial-mesenchymal transition, tumor budding has been associated with poor cancer outcomes. However, the vast heterogeneity in its exact definition, methodology of assessment, and patient stratification need to be resolved before it can be routinely used as a standardized prognostic feature. Here, we discuss the heterogeneity in defining and assessing tumor budding, its clinical significance across multiple cancer types, and its prospective implementation in clinical practice. Next, we review the emerging evidence about partial, rather than complete, epithelial-mesenchymal phenotype at the tumor bud level, and its connection with tumor proliferation, quiescence, and stemness. Finally, based on recent literature, indicating a co-expression of epithelial and mesenchymal markers in many tumor buds, we posit tumor budding to be a manifestation of this hybrid epithelial/mesenchymal phenotype displaying collective cell migration.

One of the most important roles of cells is performing their cellular tasks properly for survival. Cells usually achieve robust functionality, for example, cell-fate decision-making and signal transduction, through multiple layers of regulation involving many genes. Despite the combinatorial complexity of gene regulation, its quantitative behavior has been typically studied on the basis of experimentally verified core gene regulatory circuitry, composed of a small set of important elements. It is still unclear how such a core circuit operates in the presence of many other regulatory molecules and in a crowded and noisy cellular environment. Here we report a new computational method, named random circuit perturbation (RACIPE), for interrogating the robust dynamical behavior of a gene regulatory circuit even without accurate measurements of circuit kinetic parameters. RACIPE generates an ensemble of random kinetic models corresponding to a fixed circuit topology, and utilizes statistical tools to identify generic properties of the circuit. By applying RACIPE to simple toggle-switch-like motifs, we observed that the stable states of all models converge to experimentally observed gene state clusters even when the parameters are strongly perturbed. RACIPE was further applied to a proposed 22-gene network of the Epithelial-to-Mesenchymal Transition (EMT), from which we identified four experimentally observed gene states, including the states that are associated with two different types of hybrid Epithelial/Mesenchymal phenotypes. Our results suggest that dynamics of a gene circuit is mainly determined by its topology, not by detailed circuit parameters. Our work provides a theoretical foundation for circuit-based systems biology modeling. We anticipate RACIPE to be a powerful tool to predict and decode circuit design principles in an unbiased manner, and to quantitatively evaluate the robustness and heterogeneity of gene expression.

Hypoxia-inducible factor 1 (HIF-1) mediates a metabolic switch that blocks the conversion of pyruvate to acetyl-CoA in cancer cells. Here, we report that HIF-1α also inhibits fatty acid β-oxidation (FAO), another major source of acetyl-CoA. We identified a PGC-1β-mediated pathway by which HIF-1 inhibits the medium- and long-chain acyl-CoA dehydrogenases (MCAD and LCAD), resulting in decreased reactive oxygen species levels and enhanced proliferation. Surprisingly, we further uncovered that blocking LCAD, but not MCAD, blunts PTEN expression and dramatically affects tumor growth in vivo. Analysis of 158 liver cancer samples showed that decreased LCAD expression predicts patient mortality. Altogether, we have identified a previously unappreciated mechanism by which HIF-1 suppresses FAO to facilitate cancer progression.

Metabolic plasticity enables cancer cells to switch their metabolism phenotypes between glycolysis and oxidative phosphorylation (OXPHOS) during tumorigenesis and metastasis. However, it is still largely unknown how cancer cells orchestrate gene regulation to balance their glycolysis and OXPHOS activities. Previously, by modeling the gene regulation of cancer metabolism we have reported that cancer cells can acquire a stable hybrid metabolic state in which both glycolysis and OXPHOS can be used. Here, to comprehensively characterize cancer metabolic activity, we establish a theoretical framework by coupling gene regulation with metabolic pathways. Our modeling results demonstrate a direct association between the activities of AMPK and HIF-1, master regulators of OXPHOS and glycolysis, respectively, with the activities of three major metabolic pathways: glucose oxidation, glycolysis, and fatty acid oxidation. Our model further characterizes the hybrid metabolic state and a metabolically inactive state where cells have low activity of both glycolysis and OXPHOS. We verify the model prediction using metabolomics and transcriptomics data from paired tumor and adjacent benign tissue samples from a cohort of breast cancer patients and RNA-sequencing data from The Cancer Genome Atlas. We further validate the model prediction by in vitro studies of aggressive triple-negative breast cancer (TNBC) cells. The experimental results confirm that TNBC cells can maintain a hybrid metabolic phenotype and targeting both glycolysis and OXPHOS is necessary to eliminate their metabolic plasticity. In summary, our work serves as a platform to symmetrically study how tuning gene activity modulates metabolic pathway activity, and vice versa.

Abnormal metabolism is a hallmark of cancer, yet its regulation remains poorly understood. Cancer cells were considered to utilize primarily glycolysis for ATP production, referred to as the Warburg effect. However, recent evidence suggests that oxidative phosphorylation (OXPHOS) plays a crucial role during cancer progression. Here we utilized a systems biology approach to decipher the regulatory principle of glycolysis and OXPHOS. Integrating information from literature, we constructed a regulatory network of genes and metabolites from which we extracted a core circuit containing HIF-1, AMPK, and ROS. Our circuit analysis showed that while normal cells have an oxidative state and a glycolytic state, cancer cells can access a hybrid state with both metabolic modes coexisting. This was due to higher ROS production and/or oncogene activation, such as RAS, MYC, and c-SRC. Guided by the model, we developed two signatures consisting of AMPK and HIF-1 downstream genes, respectively, to quantify the activity of glycolysis and OXPHOS. By applying the AMPK and HIF-1 signatures to TCGA patient transcriptomics data of multiple cancer types and single-cell RNA-seq data of lung adenocarcinoma, we confirmed an anti-correlation between AMPK and HIF-1 activities and the association of metabolic states with oncogenes. We propose that the hybrid phenotype contributes to metabolic plasticity, allowing cancer cells to adapt to various microenvironments. Using model simulations, our theoretical framework of metabolism can serve as a platform to decode cancer metabolic plasticity and design cancer therapies targeting metabolism.

Metastasis of carcinoma involves migration of tumor cells to distant organs and initiate secondary tumors. Migration requires a complete or partial Epithelial-to-Mesenchymal Transition (EMT), and tumor-initiation requires cells possessing stemness. Epithelial cells (E) undergoing a complete EMT to become mesenchymal (M) have been suggested to be more likely to possess stemness. However, recent studies suggest that stemness can also be associated with cells undergoing a partial EMT (hybrid E/M phenotype). Therefore, the correlation between EMT and stemness remains elusive. Here, using a theoretical framework that couples the core EMT and stemness modules (miR-200/ZEB and LIN28/let-7), we demonstrate that the positioning of ‘stemness window’ on the ‘EMT axis’ need not be universal; rather it can be fine-tuned. Particularly, we present OVOL as an example of a modulating factor that, due to its coupling with miR-200/ZEB/LIN28/let-7 circuit, fine-tunes the EMT-stemness interplay. Coupling OVOL can inhibit the stemness likelihood of M and elevate that of the hybrid E/M (partial EMT) phenotype, thereby pulling the ‘stemness window’ away from the M end of ‘EMT axis’. Our results unify various apparently contradictory experimental findings regarding the interconnection between EMT and stemness, corroborate the emerging notion that partial EMT associates with stemness, and offer new testable predictions.