Dongya Jia

Implications of the Hybrid Epithelial/Mesenchymal Phenotype in Metastasis

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.

HIF-1-Mediated Suppression of Acyl-CoA Dehydrogenases and Fatty Acid Oxidation Is Critical for Cancer Progression

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.

Tumor Budding: The Name is EMT. Partial EMT.

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.

Stability of the hybrid epithelial/mesenchymal phenotype

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.

Coupling the modules of EMT and stemness: A tunable ‘stemness window’ model

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.

Interrogating the topological robustness of gene regulatory circuits by randomization.

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.

Phosphorylation-induced conformational dynamics in an intrinsically disordered protein and potential role in phenotypic heterogeneity

Significance The onset of androgen resistance is a major impediment in treating prostate cancer (PCa). However, the underlying molecular mechanisms are not fully understood. Here, we integrate multiple biophysical approaches that report conformational preferences of the intrinsically disordered protein (IDP) Prostate-Associated Gene 4 (PAGE4) with human cancer biology and nonlinear dynamics. Based on our biophysical, biochemical, and cellular data, a mathematical model is presented that suggests a mechanism by which phosphorylation-induced conformational changes in PAGE4 can switch from an androgen-dependent phenotype to an androgen-independent phenotype in PCa cells. The study underscores how IDPs that engage in multiple “promiscuous” interactions due to their conformational dynamics when overexpressed or aberrantly expressed can stochastically orchestrate phenotypic heterogeneity in PCa. Intrinsically disordered proteins (IDPs) that lack a unique 3D structure and comprise a large fraction of the human proteome play important roles in numerous cellular functions. Prostate-Associated Gene 4 (PAGE4) is an IDP that acts as a potentiator of the Activator Protein-1 (AP-1) transcription factor. Homeodomain-Interacting Protein Kinase 1 (HIPK1) phosphorylates PAGE4 at S9 and T51, but only T51 is critical for its activity. Here, we identify a second kinase, CDC-Like Kinase 2 (CLK2), which acts on PAGE4 and hyperphosphorylates it at multiple S/T residues, including S9 and T51. We demonstrate that HIPK1 is expressed in both androgen-dependent and androgen-independent prostate cancer (PCa) cells, whereas CLK2 and PAGE4 are expressed only in androgen-dependent cells. Cell-based studies indicate that PAGE4 interaction with the two kinases leads to opposing functions. HIPK1-phosphorylated PAGE4 (HIPK1-PAGE4) potentiates c-Jun, whereas CLK2-phosphorylated PAGE4 (CLK2-PAGE4) attenuates c-Jun activity. Consistent with the cellular data, biophysical measurements (small-angle X-ray scattering, single-molecule fluorescence resonance energy transfer, and NMR) indicate that HIPK1-PAGE4 exhibits a relatively compact conformational ensemble that binds AP-1, whereas CLK2-PAGE4 is more expanded and resembles a random coil with diminished affinity for AP-1. Taken together, the results suggest that the phosphorylation-induced conformational dynamics of PAGE4 may play a role in modulating changes between PCa cell phenotypes. A mathematical model based on our experimental data demonstrates how differential phosphorylation of PAGE4 can lead to transitions between androgen-dependent and androgen-independent phenotypes by altering the AP-1/androgen receptor regulatory circuit in PCa cells.

Modeling the Genetic Regulation of Cancer Metabolism: Interplay between Glycolysis and Oxidative Phosphorylation

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 The Cancer Genome Atlas patient transcriptomics data of multiple cancer types and single-cell RNA-seq data of lung adenocarcinoma, we confirmed an anticorrelation 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. Cancer Res; 77(7); 1564-74. ©2017 AACR.

The GRHL2/ZEB Feedback Loop—A Key Axis in the Regulation of EMT in Breast Cancer

More than 90% of cancer‐related deaths are caused by metastasis. Epithelial‐to‐Mesenchymal Transition (EMT) causes tumor cell dissemination while the reverse process, Mesenchymal‐to‐Epithelial Transition (MET) allows cancer cells to grow and establish a potentially deadly metastatic lesion. Recent evidence indicates that in addition to E and M, cells can adopt a stable hybrid Epithelial/Mesenchymal (E/M) state where they can move collectively leading to clusters of Circulating Tumor Cells—the “bad actors” of metastasis. EMT is postulated to occur in all four major histological breast cancer subtypes. Here, we identify a set of genes strongly correlated with CDH1 in 877 cancer cell lines, and differentially expressed genes in cell lines overexpressing ZEB1, SNAIL, and TWIST. GRHL2 and ESRP1 appear in both these sets and also correlate with CDH1 at the protein level in 40 breast cancer specimens. Next, we find that GRHL2 and CD24 expression coincide with an epithelial character in human mammary epithelial cells. Further, we show that high GRHL2 expression is highly correlated with worse relapse‐free survival in all four subtypes of breast cancer. Finally, we integrate CD24, GRHL2, and ESRP1 into a mathematical model of EMT regulation to validate the role of these players in EMT. Our data analysis and modeling results highlight the relationships among multiple crucial EMT/MET drivers including ZEB1, GRHL2, CD24, and ESRP1, particularly in basal‐like breast cancers, which are most similar to triple‐negative breast cancer (TNBC) and are considered the most dangerous subtype. J. Cell. Biochem. 118: 2559–2570, 2017.

Phenotypic Plasticity and Cell Fate Decisions in Cancer: Insights from Dynamical Systems Theory

Waddington’s epigenetic landscape, a famous metaphor in developmental biology, depicts how a stem cell progresses from an undifferentiated phenotype to a differentiated one. The concept of “landscape” in the context of dynamical systems theory represents a high-dimensional space, in which each cell phenotype is considered as an “attractor” that is determined by interactions between multiple molecular players, and is buffered against environmental fluctuations. In addition, biological noise is thought to play an important role during these cell-fate decisions and in fact controls transitions between different phenotypes. Here, we discuss the phenotypic transitions in cancer from a dynamical systems perspective and invoke the concept of “cancer attractors”—hidden stable states of the underlying regulatory network that are not occupied by normal cells. Phenotypic transitions in cancer occur at varying levels depending on the context. Using epithelial-to-mesenchymal transition (EMT), cancer stem-like properties, metabolic reprogramming and the emergence of therapy resistance as examples, we illustrate how phenotypic plasticity in cancer cells enables them to acquire hybrid phenotypes (such as hybrid epithelial/mesenchymal and hybrid metabolic phenotypes) that tend to be more aggressive and notoriously resilient to therapies such as chemotherapy and androgen-deprivation therapy. Furthermore, we highlight multiple factors that may give rise to phenotypic plasticity in cancer cells, such as (a) multi-stability or oscillatory behaviors governed by underlying regulatory networks involved in cell-fate decisions in cancer cells, and (b) network rewiring due to conformational dynamics of intrinsically disordered proteins (IDPs) that are highly enriched in cancer cells. We conclude by discussing why a therapeutic approach that promotes “recanalization”, i.e., the exit from “cancer attractors” and re-entry into “normal attractors”, is more likely to succeed rather than a conventional approach that targets individual molecules/pathways.