1
|
Daniels BC, Wang Y, Page RE, Amdam GV. Identifying a developmental transition in honey bees using gene expression data. PLoS Comput Biol 2023; 19:e1010704. [PMID: 37733808 PMCID: PMC10547183 DOI: 10.1371/journal.pcbi.1010704] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2022] [Revised: 10/03/2023] [Accepted: 09/05/2023] [Indexed: 09/23/2023] Open
Abstract
In many organisms, interactions among genes lead to multiple functional states, and changes to interactions can lead to transitions into new states. These transitions can be related to bifurcations (or critical points) in dynamical systems theory. Characterizing these collective transitions is a major challenge for systems biology. Here, we develop a statistical method for identifying bistability near a continuous transition directly from high-dimensional gene expression data. We apply the method to data from honey bees, where a known developmental transition occurs between bees performing tasks in the nest and leaving the nest to forage. Our method, which makes use of the expected shape of the distribution of gene expression levels near a transition, successfully identifies the emergence of bistability and links it to genes that are known to be involved in the behavioral transition. This proof of concept demonstrates that going beyond correlative analysis to infer the shape of gene expression distributions might be used more generally to identify collective transitions from gene expression data.
Collapse
Affiliation(s)
- Bryan C. Daniels
- School of Complex Adaptive Systems, Arizona State University, Tempe, Arizona, United States of America
| | - Ying Wang
- Banner Health Corporation, Phoenix, Arizona, United States of America
| | - Robert E. Page
- School of Life Sciences, Arizona State University, Tempe, Arizona, United States of America
- Department of Entomology and Nematology, University of California Davis, Davis, California, United States of America
| | - Gro V. Amdam
- School of Life Sciences, Arizona State University, Tempe, Arizona, United States of America
- Department of Ecology and Natural Resource Management, Norwegian University of Life Sciences, Aas, Norway
| |
Collapse
|
2
|
Li G, LeFebre R, Starman A, Chappell P, Mugler A, Sun B. Temporal signals drive the emergence of multicellular information networks. Proc Natl Acad Sci U S A 2022; 119:e2202204119. [PMID: 36067282 PMCID: PMC9477235 DOI: 10.1073/pnas.2202204119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2022] [Accepted: 07/18/2022] [Indexed: 11/18/2022] Open
Abstract
Coordinated responses to environmental stimuli are critical for multicellular organisms. To overcome the obstacles of cell-to-cell heterogeneity and noisy signaling dynamics within individual cells, cells must effectively exchange information with peers. However, the dynamics and mechanisms of collective information transfer driven by external signals are poorly understood. Here we investigate the calcium dynamics of neuronal cells that form confluent monolayers and respond to cyclic ATP stimuli in microfluidic devices. Using Granger inference to reconstruct the underlying causal relations between the cells, we find that the cells self-organize into spatially decentralized and temporally stationary networks to support information transfer via gap junction channels. The connectivity of the causal networks depends on the temporal profile of the external stimuli, where short periods, or long periods with small duty fractions, lead to reduced connectivity and fractured network topology. We build a theoretical model based on communicating excitable units that reproduces our observations. The model further predicts that connectivity of the causal network is maximal at an optimal communication strength, which is confirmed by the experiments. Together, our results show that information transfer between neuronal cells is externally regulated by the temporal profile of the stimuli and internally regulated by cell-cell communication.
Collapse
Affiliation(s)
- Guanyu Li
- Department of Physics, Oregon State University, Corvallis, OR 97331
| | - Ryan LeFebre
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA 15260
| | - Alia Starman
- Department of Biomedical Sciences, Carlson College of Veterinary Medicine, Oregon State University, Corvallis, OR 97331
| | - Patrick Chappell
- Department of Biomedical Sciences, Carlson College of Veterinary Medicine, Oregon State University, Corvallis, OR 97331
| | - Andrew Mugler
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA 15260
| | - Bo Sun
- Department of Physics, Oregon State University, Corvallis, OR 97331
| |
Collapse
|
3
|
Emergent phenomena in living systems: A statistical mechanical perspective. J Biosci 2022. [DOI: 10.1007/s12038-021-00247-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
4
|
Bose I. Tipping the Balance: A Criticality Perspective. ENTROPY 2022; 24:e24030405. [PMID: 35327916 PMCID: PMC8947304 DOI: 10.3390/e24030405] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/17/2022] [Revised: 03/10/2022] [Accepted: 03/12/2022] [Indexed: 01/02/2023]
Abstract
Cell populations are often characterised by phenotypic heterogeneity in the form of two distinct subpopulations. We consider a model of tumour cells consisting of two subpopulations: non-cancer promoting (NCP) and cancer-promoting (CP). Under steady state conditions, the model has similarities with a well-known model of population genetics which exhibits a purely noise-induced transition from unimodality to bimodality at a critical value of the noise intensity σ2. The noise is associated with the parameter λ representing the system-environment coupling. In the case of the tumour model, λ has a natural interpretation in terms of the tissue microenvironment which has considerable influence on the phenotypic composition of the tumour. Oncogenic transformations give rise to considerable fluctuations in the parameter. We compute the λ−σ2 phase diagram in a stochastic setting, drawing analogies between bifurcations and phase transitions. In the region of bimodality, a transition from a state of balance to a state of dominance, in terms of the competing subpopulations, occurs at λ = 0. Away from this point, the NCP (CP) subpopulation becomes dominant as λ changes towards positive (negative) values. The variance of the steady state probability density function as well as two entropic measures provide characteristic signatures at the transition point.
Collapse
Affiliation(s)
- Indrani Bose
- Department of Physics, Bose Institute, 93/1, A. P. C. Road, Kolkata 700009, India
| |
Collapse
|
5
|
Vennettilli M, Erez A, Mugler A. Multicellular sensing at a feedback-induced critical point. Phys Rev E 2020; 102:052411. [PMID: 33327087 DOI: 10.1103/physreve.102.052411] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2020] [Accepted: 10/06/2020] [Indexed: 11/07/2022]
Abstract
Feedback in sensory biochemical networks can give rise to bifurcations in cells' behavioral response. These bifurcations share many properties with thermodynamic critical points. Evidence suggests that biological systems may operate near these critical points, but the functional benefit of doing so remains poorly understood. Here we investigate a simple biochemical model with nonlinear feedback and multicellular communication to determine if criticality provides a functional benefit in terms of the ability to gain information about a stochastic chemical signal. We find that when signal fluctuations are slow, the mutual information between the signal and the intracellular readout is maximized at criticality, because the benefit of high signal susceptibility outweighs the detriment of high readout noise. When cells communicate, criticality gives rise to long-range correlations in readout molecule number among cells. Consequently, we find that communication increases the mutual information between a given cell's readout and the spatial average of the signal across the population. Finally, we find that both with and without communication, the sensory benefits of criticality compete with critical slowing down, such that the information rate, as opposed to the information itself, is minimized at the critical point. Our results reveal the costs and benefits of feedback-induced criticality for multicellular sensing.
Collapse
Affiliation(s)
- Michael Vennettilli
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, USA
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Amir Erez
- Department of Molecular Biology, Princeton University, Princeton, New Jersey 08544, USA
- The Racah Institute of Physics, Hebrew University, Jerusalem 91904, Israel
| | - Andrew Mugler
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, USA
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| |
Collapse
|
6
|
Erez A, Byrd TA, Vennettilli M, Mugler A. Cell-to-Cell Information at a Feedback-Induced Bifurcation Point. PHYSICAL REVIEW LETTERS 2020; 125:048103. [PMID: 32794792 DOI: 10.1103/physrevlett.125.048103] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/31/2019] [Accepted: 07/06/2020] [Indexed: 06/11/2023]
Abstract
A ubiquitous way that cells share information is by exchanging molecules. Yet, the fundamental ways that this information exchange is influenced by intracellular dynamics remain unclear. Here we use information theory to investigate a simple model of two interacting cells with internal feedback. We show that cell-to-cell molecule exchange induces a collective two-cell critical point and that the mutual information between the cells peaks at this critical point. Information can remain large far from the critical point on a manifold of cellular states but scales logarithmically with the correlation time of the system, resulting in an information-correlation time trade-off. This trade-off is strictly imposed, suggesting the correlation time as a proxy for the mutual information.
Collapse
Affiliation(s)
- Amir Erez
- Department of Molecular Biology, Princeton University, Princeton, New Jersey 08544, USA
| | - Tommy A Byrd
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, USA
| | - Michael Vennettilli
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, USA
| | - Andrew Mugler
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, USA
| |
Collapse
|
7
|
Byrd TA, Erez A, Vogel RM, Peterson C, Vennettilli M, Altan-Bonnet G, Mugler A. Critical slowing down in biochemical networks with feedback. Phys Rev E 2019; 100:022415. [PMID: 31574667 PMCID: PMC8499154 DOI: 10.1103/physreve.100.022415] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2019] [Indexed: 06/10/2023]
Abstract
Near a bifurcation point, the response time of a system is expected to diverge due to the phenomenon of critical slowing down. We investigate critical slowing down in well-mixed stochastic models of biochemical feedback by exploiting a mapping to the mean-field Ising universality class. We analyze the responses to a sudden quench and to continuous driving in the model parameters. In the latter case, we demonstrate that our class of models exhibits the Kibble-Zurek collapse, which predicts the scaling of hysteresis in cellular responses to gradual perturbations. We discuss the implications of our results in terms of the tradeoff between a precise and a fast response. Finally, we use our mapping to quantify critical slowing down in T cells, where the addition of a drug is equivalent to a sudden quench in parameter space.
Collapse
Affiliation(s)
- Tommy A. Byrd
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, USA
| | - Amir Erez
- Department of Molecular Biology, Princeton University, Princeton, New Jersey 08544, USA
| | - Robert M. Vogel
- IBM T. J. Watson Research Center, Yorktown Heights, New York 10598, USA
| | - Curtis Peterson
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, USA
- Department of Physics and School of Mathematical and Statistical Sciences, Arizona State University, Tempe, Arizona 85287, USA
| | - Michael Vennettilli
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, USA
| | - Grégoire Altan-Bonnet
- Immunodynamics Group, Cancer and Inflammation Program, National Cancer Institute, National Institutes of Health, Bethesda, Maryland 20814, USA
| | - Andrew Mugler
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, USA
| |
Collapse
|