1
|
Bergman DR, Karikomi MK, Yu M, Nie Q, MacLean AL. Modeling the effects of EMT-immune dynamics on carcinoma disease progression. Commun Biol 2021; 4:983. [PMID: 34408236 PMCID: PMC8373868 DOI: 10.1038/s42003-021-02499-y] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2020] [Accepted: 07/27/2021] [Indexed: 02/07/2023] Open
Abstract
During progression from carcinoma in situ to an invasive tumor, the immune system is engaged in complex sets of interactions with various tumor cells. Tumor cell plasticity alters disease trajectories via epithelial-to-mesenchymal transition (EMT). Several of the same pathways that regulate EMT are involved in tumor-immune interactions, yet little is known about the mechanisms and consequences of crosstalk between these regulatory processes. Here we introduce a multiscale evolutionary model to describe tumor-immune-EMT interactions and their impact on epithelial cancer progression from in situ to invasive disease. Through simulation of patient cohorts in silico, the model predicts that a controllable region maximizes invasion-free survival. This controllable region depends on properties of the mesenchymal tumor cell phenotype: its growth rate and its immune-evasiveness. In light of the model predictions, we analyze EMT-inflammation-associated data from The Cancer Genome Atlas, and find that association with EMT worsens invasion-free survival probabilities. This result supports the predictions of the model, and leads to the identification of genes that influence outcomes in bladder and uterine cancer, including FGF pathway members. These results suggest new means to delay disease progression, and demonstrate the importance of studying cancer-immune interactions in light of EMT.
Collapse
Affiliation(s)
- Daniel R. Bergman
- grid.266093.80000 0001 0668 7243Department of Mathematics, University of California, Irvine, CA USA
| | - Matthew K. Karikomi
- grid.266093.80000 0001 0668 7243Department of Mathematics, University of California, Irvine, CA USA
| | - Min Yu
- grid.42505.360000 0001 2156 6853USC Norris Comprehensive Cancer Center, Keck School of Medicine of the University of Southern California, Los Angeles, CA, USA ,grid.42505.360000 0001 2156 6853Department of Stem Cell Biology and Regenerative Medicine, Keck School of Medicine, University of Southern California, Los Angeles, CA USA
| | - Qing Nie
- grid.266093.80000 0001 0668 7243Department of Mathematics, University of California, Irvine, CA USA ,grid.266093.80000 0001 0668 7243Department of Cell and Developmental Biology, University of California, Irvine, CA USA
| | - Adam L. MacLean
- grid.266093.80000 0001 0668 7243Department of Mathematics, University of California, Irvine, CA USA ,grid.42505.360000 0001 2156 6853USC Norris Comprehensive Cancer Center, Keck School of Medicine of the University of Southern California, Los Angeles, CA, USA ,grid.42505.360000 0001 2156 6853Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, CA, USA
| |
Collapse
|
2
|
Qiu Y, Fung L, Schilling TF, Nie Q. Multiple morphogens and rapid elongation promote segmental patterning during development. PLoS Comput Biol 2021; 17:e1009077. [PMID: 34161317 PMCID: PMC8259987 DOI: 10.1371/journal.pcbi.1009077] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2020] [Revised: 07/06/2021] [Accepted: 05/13/2021] [Indexed: 12/21/2022] Open
Abstract
The vertebrate hindbrain is segmented into rhombomeres (r) initially defined by distinct domains of gene expression. Previous studies have shown that noise-induced gene regulation and cell sorting are critical for the sharpening of rhombomere boundaries, which start out rough in the forming neural plate (NP) and sharpen over time. However, the mechanisms controlling simultaneous formation of multiple rhombomeres and accuracy in their sizes are unclear. We have developed a stochastic multiscale cell-based model that explicitly incorporates dynamic morphogenetic changes (i.e. convergent-extension of the NP), multiple morphogens, and gene regulatory networks to investigate the formation of rhombomeres and their corresponding boundaries in the zebrafish hindbrain. During pattern initiation, the short-range signal, fibroblast growth factor (FGF), works together with the longer-range morphogen, retinoic acid (RA), to specify all of these boundaries and maintain accurately sized segments with sharp boundaries. At later stages of patterning, we show a nonlinear change in the shape of rhombomeres with rapid left-right narrowing of the NP followed by slower dynamics. Rapid initial convergence improves boundary sharpness and segment size by regulating cell sorting and cell fate both independently and coordinately. Overall, multiple morphogens and tissue dynamics synergize to regulate the sizes and boundaries of multiple segments during development. In segmental pattern formation, chemical gradients control gene expression in a concentration-dependent manner to specify distinct gene expression domains. Despite the stochasticity inherent to such biological processes, precise and accurate borders form between segmental gene expression domains. Previous work has revealed synergy between gene regulation and cell sorting in sharpening borders that are initially rough. However, it is still poorly understood how size and boundary sharpness of multiple segments are regulated in a tissue that changes dramatically in its morphology as the embryo develops. Here we develop a stochastic multiscale cell-base model to investigate these questions. Two novel strategies synergize to promote accurate segment formation, a combination of long- and short-range morphogens plus rapid tissue convergence, with one responsible for pattern initiation and the other enabling pattern refinement.
Collapse
Affiliation(s)
- Yuchi Qiu
- Department of Mathematics, University of California, Irvine, California, United States of America
| | - Lianna Fung
- Department of Developmental and Cell Biology, University of California, Irvine, California, United States of America
- The NSF-Simons Center for Multiscale Cell Fate Research, University of California, Irvine, California, United States of America
| | - Thomas F. Schilling
- Department of Developmental and Cell Biology, University of California, Irvine, California, United States of America
- The NSF-Simons Center for Multiscale Cell Fate Research, University of California, Irvine, California, United States of America
- * E-mail: (TFS); (QN)
| | - Qing Nie
- Department of Mathematics, University of California, Irvine, California, United States of America
- Department of Developmental and Cell Biology, University of California, Irvine, California, United States of America
- The NSF-Simons Center for Multiscale Cell Fate Research, University of California, Irvine, California, United States of America
- * E-mail: (TFS); (QN)
| |
Collapse
|
3
|
Sha Y, Wang S, Zhou P, Nie Q. Inference and multiscale model of epithelial-to-mesenchymal transition via single-cell transcriptomic data. Nucleic Acids Res 2020; 48:9505-9520. [PMID: 32870263 PMCID: PMC7515733 DOI: 10.1093/nar/gkaa725] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2020] [Revised: 07/19/2020] [Accepted: 08/20/2020] [Indexed: 12/17/2022] Open
Abstract
Rapid growth of single-cell transcriptomic data provides unprecedented opportunities for close scrutinizing of dynamical cellular processes. Through investigating epithelial-to-mesenchymal transition (EMT), we develop an integrative tool that combines unsupervised learning of single-cell transcriptomic data and multiscale mathematical modeling to analyze transitions during cell fate decision. Our approach allows identification of individual cells making transition between all cell states, and inference of genes that drive transitions. Multiscale extractions of single-cell scale outputs naturally reveal intermediate cell states (ICS) and ICS-regulated transition trajectories, producing emergent population-scale models to be explored for design principles. Testing on the newly designed single-cell gene regulatory network model and applying to twelve published single-cell EMT datasets in cancer and embryogenesis, we uncover the roles of ICS on adaptation, noise attenuation, and transition efficiency in EMT, and reveal their trade-off relations. Overall, our unsupervised learning method is applicable to general single-cell transcriptomic datasets, and our integrative approach at single-cell resolution may be adopted for other cell fate transition systems beyond EMT.
Collapse
Affiliation(s)
- Yutong Sha
- Department of Mathematics, University of California, Irvine, Irvine, CA 92697, USA.,The NSF-Simons Center for Multiscale Cell Fate Research, University of California, Irvine, Irvine, CA 92697, USA
| | - Shuxiong Wang
- Department of Mathematics, University of California, Irvine, Irvine, CA 92697, USA
| | - Peijie Zhou
- Department of Mathematics, University of California, Irvine, Irvine, CA 92697, USA
| | - Qing Nie
- Department of Mathematics, University of California, Irvine, Irvine, CA 92697, USA.,The NSF-Simons Center for Multiscale Cell Fate Research, University of California, Irvine, Irvine, CA 92697, USA.,Department of Developmental and Cell Biology, University of California, Irvine, Irvine, CA 92697, USA
| |
Collapse
|
4
|
Qiu Y, Chen W, Nie Q. STOCHASTIC DYNAMICS OF CELL LINEAGE IN TISSUE HOMEOSTASIS. DISCRETE AND CONTINUOUS DYNAMICAL SYSTEMS-SERIES B 2019; 24:3971-3994. [PMID: 32269502 PMCID: PMC7141575 DOI: 10.3934/dcdsb.2018339] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
During epithelium tissue maintenance, lineages of cells differentiate and proliferate in a coordinated way to provide the desirable size and spatial organization of different types of cells. While mathematical models through deterministic description have been used to dissect role of feedback regulations on tissue layer size and stratification, how the stochastic effects influence tissue maintenance remains largely unknown. Here we present a stochastic continuum model for cell lineages to investigate how both layer thickness and layer stratification are affected by noise. We find that the cell-intrinsic noise often causes reduction and oscillation of layer size whereas the cell-extrinsic noise increases the thickness, and sometimes, leads to uncontrollable growth of the tissue layer. The layer stratification usually deteriorates as the noise level increases in the cell lineage systems. Interestingly, the morphogen noise, which mixes both cell-intrinsic noise and cell-extrinsic noise, can lead to larger size of layer with little impact on the layer stratification. By investigating different combinations of the three types of noise, we find the layer thickness variability is reduced when cell-extrinsic noise level is high or morphogen noise level is low. Interestingly, there exists a tradeoff between low thickness variability and strong layer stratification due to competition among the three types of noise, suggesting robust layer homeostasis requires balanced levels of different types of noise in the cell lineage systems.
Collapse
Affiliation(s)
- Yuchi Qiu
- Department of Mathematics, University of California, Irvine Irvine, CA 92697, USA
| | - Weitao Chen
- Department of Mathematics, University of California, Riverside Riverside, CA 92507, USA
| | | |
Collapse
|
5
|
Wang S, Moise AR. Recent insights on the role and regulation of retinoic acid signaling during epicardial development. Genesis 2019; 57:e23303. [PMID: 31066193 PMCID: PMC6682438 DOI: 10.1002/dvg.23303] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2019] [Revised: 04/23/2019] [Accepted: 04/24/2019] [Indexed: 12/18/2022]
Abstract
The vitamin A metabolite, retinoic acid, carries out essential and conserved roles in vertebrate heart development. Retinoic acid signals via retinoic acid receptors (RAR)/retinoid X receptors (RXRs) heterodimers to induce the expression of genes that control cell fate specification, proliferation, and differentiation. Alterations in retinoic acid levels are often associated with congenital heart defects. Therefore, embryonic levels of retinoic acid need to be carefully regulated through the activity of enzymes, binding proteins and transporters involved in vitamin A metabolism. Here, we review evidence of the complex mechanisms that control the fetal uptake and synthesis of retinoic acid from vitamin A precursors. Next, we highlight recent evidence of the role of retinoic acid in orchestrating myocardial compact zone growth and coronary vascular development.
Collapse
Affiliation(s)
- Suya Wang
- Department of Cardiology, Boston Children's Hospital, 300 Longwood Avenue, Boston, MA, 02115, USA
| | - Alexander R. Moise
- Medical Sciences Division, Northern Ontario School of Medicine, Sudbury, ON P3E 2C6, Canada
- Departments of Chemistry and Biochemistry, and Biology and Biomolecular Sciences Program, Laurentian University, Sudbury, ON, P3E 2C6 Canada
- Department of Pharmacology and Toxicology, School of Pharmacy, University of Kansas, Lawrence, KS, 66045, USA
| |
Collapse
|
6
|
Zhornikova P, Golyandina N, Spirov A. Noise model estimation with application to gene expression. J Bioinform Comput Biol 2019; 17:1950009. [PMID: 31057070 DOI: 10.1142/s0219720019500094] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Algorithms for the estimation of noise level and the detection of noise model are proposed. They are applied to gene expression data for Drosophila embryos. The 2D data on gene expression and the extracted 1D profiles are considered. Since the 1D data contain processing errors, an algorithm for separation of these processing errors is constructed to estimate the biological noise level. An approach to discrimination between the additive and multiplicative models is suggested for the 1D and 2D cases. Singular spectrum analysis and its 2D extension are exploited for the pattern extraction. The algorithms are tested on artificial data similar to the real data. Comparison of the results, which are obtained by the 1D and 2D methods, is performed for Krüppel and giant genes.
Collapse
Affiliation(s)
- Polina Zhornikova
- * Faculty of Mathematics and Mechanics, St. Petersburg State University, Universitetskaya Nab. 7/9, 199034 St. Petersburg, Russia
| | - Nina Golyandina
- * Faculty of Mathematics and Mechanics, St. Petersburg State University, Universitetskaya Nab. 7/9, 199034 St. Petersburg, Russia
| | - Alexander Spirov
- † The Sechenov Institute of Evolutionary Physiology and Biochemistry Russian Academy of Sciences, Torez Pr. 44, 194223 St. Petersburg, Russia
| |
Collapse
|
7
|
Viner H, Nitsan I, Sapir L, Drori S, Tzlil S. Mechanical Communication Acts as a Noise Filter. iScience 2019; 14:58-68. [PMID: 30927696 PMCID: PMC6441679 DOI: 10.1016/j.isci.2019.02.030] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2018] [Revised: 12/29/2018] [Accepted: 02/26/2019] [Indexed: 11/15/2022] Open
Abstract
Cells can communicate mechanically by responding to mechanical deformations generated by their neighbors. Here, we describe a new role for mechanical communication by demonstrating that mechanical coupling between cells acts as a signaling cue that reduces intrinsic noise in the interacting cells. We measure mechanical interaction between beating cardiac cells cultured on a patterned flexible substrate and find that beat-to-beat variability decays exponentially with coupling strength. To demonstrate that such noise reduction is indeed a direct consequence of mechanical coupling, we reproduce the exponential decay in an assay where a beating cell interacts mechanically with an artificial stochastic ‘mechanical cell’. The mechanical cell consists of a probe that mimics the deformations generated by a stochastically beating neighboring cardiac cell. We show that noise reduction through mechanical coupling persists long after stimulation stops and identify microtubule integrity, NOX2, and CaMKII as mediators of noise reduction. Mechanical communication reduces intrinsic noise in interacting cells Cardiac cell beating noise decays exponentially with the strength of mechanical coupling Identical exponential decay length is obtained using a stochastic mechanical cell NOX2, ROS, and CaMKII are involved in mechanical communication-induced noise reduction
Collapse
Affiliation(s)
- Hen Viner
- Faculty of Mechanical Engineering, Technion - Israel Institute of Technology, Haifa 32000, Israel
| | - Ido Nitsan
- Faculty of Mechanical Engineering, Technion - Israel Institute of Technology, Haifa 32000, Israel
| | - Liel Sapir
- Faculty of Mechanical Engineering, Technion - Israel Institute of Technology, Haifa 32000, Israel
| | - Stavit Drori
- Faculty of Mechanical Engineering, Technion - Israel Institute of Technology, Haifa 32000, Israel
| | - Shelly Tzlil
- Faculty of Mechanical Engineering, Technion - Israel Institute of Technology, Haifa 32000, Israel.
| |
Collapse
|
8
|
Abstract
The transition of epithelial cells into a mesenchymal state (epithelial-to-mesenchymal transition or EMT) is a highly dynamic process implicated in various biological processes. During EMT, cells do not necessarily exist in 'pure' epithelial or mesenchymal states. There are cells with mixed (or hybrid) features of the two, which are termed as the intermediate cell states (ICSs). While the exact functions of ICS remain elusive, together with EMT it appears to play important roles in embryogenesis, tissue development, and pathological processes such as cancer metastasis. Recent single cell experiments and advanced mathematical modeling have improved our capability in identifying ICS and provided a better understanding of ICS in development and disease. Here, we review the recent findings related to the ICS in/or EMT and highlight the challenges in the identification and functional characterization of ICS.
Collapse
Affiliation(s)
- Yutong Sha
- Department of Mathematics, University of California, Irvine, CA 92697, United States of America
- Co-first authors
| | - Daniel Haensel
- Department of Biological Chemistry, School of Medicine, University of California, Irvine, CA 92697, United States of America
- Co-first authors
| | - Guadalupe Gutierrez
- Department of Biological Chemistry, School of Medicine, University of California, Irvine, CA 92697, United States of America
| | - Huijing Du
- Department of Mathematics, University of Nebraska-Lincoln, Lincoln, NE 68588, United States of America
| | - Xing Dai
- Department of Biological Chemistry, School of Medicine, University of California, Irvine, CA 92697, United States of America
| | - Qing Nie
- Department of Mathematics, University of California, Irvine, CA 92697, United States of America
- Department of Development and Cell Biology, University of California, Irvine, CA 92697, United States of America
| |
Collapse
|