1
|
Eidi Z, Khorasani N, Sadeghi M. Correspondence between multiple signaling and developmental cellular patterns: a computational perspective. Front Cell Dev Biol 2024; 12:1310265. [PMID: 39139453 PMCID: PMC11319269 DOI: 10.3389/fcell.2024.1310265] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Accepted: 07/02/2024] [Indexed: 08/15/2024] Open
Abstract
The spatial arrangement of variant phenotypes during stem cell division plays a crucial role in the self-organization of cell tissues. The patterns observed in these cellular assemblies, where multiple phenotypes vie for space and resources, are largely influenced by a mixture of different diffusible chemical signals. This complex process is carried out within a chronological framework of interplaying intracellular and intercellular events. This includes receiving external stimulants, whether secreted by other individuals or provided by the environment, interpreting these environmental signals, and incorporating the information to designate cell fate. Here, given two distinct signaling patterns generated by Turing systems, we investigated the spatial distribution of differentiating cells that use these signals as external cues for modifying the production rates. By proposing a computational map, we show that there is a correspondence between the multiple signaling and developmental cellular patterns. In other words, the model provides an appropriate prediction for the final structure of the differentiated cells in a multi-signal, multi-cell environment. Conversely, when a final snapshot of cellular patterns is given, our algorithm can partially identify the signaling patterns that influenced the formation of the cellular structure, provided that the governing dynamic of the signaling patterns is already known.
Collapse
Affiliation(s)
- Zahra Eidi
- School of Biological Sciences, Institute for Research in Fundamental Sciences (IPM), Tehran, Iran
| | - Najme Khorasani
- School of Biological Sciences, Institute for Research in Fundamental Sciences (IPM), Tehran, Iran
| | - Mehdi Sadeghi
- National Institute of Genetic Engineering and Biotechnology (NIGEB), Tehran, Iran
| |
Collapse
|
2
|
Nou XA, Voigt CA. Sentinel cells programmed to respond to environmental DNA including human sequences. Nat Chem Biol 2024; 20:211-220. [PMID: 37770697 DOI: 10.1038/s41589-023-01431-1] [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: 08/08/2022] [Accepted: 08/31/2023] [Indexed: 09/30/2023]
Abstract
Monitoring environmental DNA can track the presence of organisms, from viruses to animals, but requires continuous sampling of transient sequences from a complex milieu. Here we designed living sentinels using Bacillus subtilis to report the uptake of a DNA sequence after matching it to a preencoded target. Overexpression of ComK increased DNA uptake 3,000-fold, allowing for femtomolar detection in samples dominated by background DNA. This capability was demonstrated using human sequences containing single-nucleotide polymorphisms (SNPs) associated with facial features. Sequences were recorded with high efficiency and were protected from nucleases for weeks. The SNP could be determined by sequencing or in vivo using CRISPR interference to turn on reporter expression in response to a specific base. Multiple SNPs were recorded by one cell or through a consortium in which each member recorded a different sequence. Sentinel cells could surveil for specific sequences over long periods of time for applications spanning forensics, ecology and epidemiology.
Collapse
Affiliation(s)
- Xuefei Angelina Nou
- Synthetic Biology Center, Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Christopher A Voigt
- Synthetic Biology Center, Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA, USA.
| |
Collapse
|
3
|
Le Priol C, Azencott CA, Gidrol X. Detection of genes with differential expression dispersion unravels the role of autophagy in cancer progression. PLoS Comput Biol 2023; 19:e1010342. [PMID: 36893104 PMCID: PMC9997931 DOI: 10.1371/journal.pcbi.1010342] [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: 07/01/2022] [Accepted: 02/09/2023] [Indexed: 03/10/2023] Open
Abstract
The majority of gene expression studies focus on the search for genes whose mean expression is different between two or more populations of samples in the so-called "differential expression analysis" approach. However, a difference in variance in gene expression may also be biologically and physiologically relevant. In the classical statistical model used to analyze RNA-sequencing (RNA-seq) data, the dispersion, which defines the variance, is only considered as a parameter to be estimated prior to identifying a difference in mean expression between conditions of interest. Here, we propose to evaluate four recently published methods, which detect differences in both the mean and dispersion in RNA-seq data. We thoroughly investigated the performance of these methods on simulated datasets and characterized parameter settings to reliably detect genes with a differential expression dispersion. We applied these methods to The Cancer Genome Atlas datasets. Interestingly, among the genes with an increased expression dispersion in tumors and without a change in mean expression, we identified some key cellular functions, most of which were related to catabolism and were overrepresented in most of the analyzed cancers. In particular, our results highlight autophagy, whose role in cancerogenesis is context-dependent, illustrating the potential of the differential dispersion approach to gain new insights into biological processes and to discover new biomarkers.
Collapse
Affiliation(s)
- Christophe Le Priol
- Univ. Grenoble Alpes, INSERM, CEA-IRIG, Biomics, Grenoble, France
- * E-mail: (CLP); (XG)
| | - Chloé-Agathe Azencott
- Center for Computational Biology, Mines ParisTech, PSL Research University, Paris, France
- Institut Curie, Paris, France
- INSERM U900, Paris, France
| | - Xavier Gidrol
- Univ. Grenoble Alpes, INSERM, CEA-IRIG, Biomics, Grenoble, France
- * E-mail: (CLP); (XG)
| |
Collapse
|
4
|
Khorasani N, Sadeghi M. A computational model of stem cells' decision-making mechanism to maintain tissue homeostasis and organization in the presence of stochasticity. Sci Rep 2022; 12:9167. [PMID: 35654903 PMCID: PMC9163052 DOI: 10.1038/s41598-022-12717-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2021] [Accepted: 05/10/2022] [Indexed: 11/09/2022] Open
Abstract
The maintenance of multi-cellular developed tissue depends on the proper cell production rate to replace the cells destroyed by the programmed process of cell death. The stem cell is the main source of producing cells in a developed normal tissue. It makes the stem cell the lead role in the scene of a fully formed developed tissue to fulfill its proper functionality. By focusing on the impact of stochasticity, here, we propose a computational model to reveal the internal mechanism of a stem cell, which generates the right proportion of different types of specialized cells, distribute them into their right position, and in the presence of intercellular reactions, maintain the organized structure in a homeostatic state. The result demonstrates that the spatial pattern could be harassed by the population geometries. Besides, it clearly shows that our model with progenitor cells able to recover the stem cell presence could retrieve the initial pattern appropriately in the case of injury. One of the fascinating outcomes of this project is demonstrating the contradictory roles of stochasticity. It breaks the proper boundaries of the initial spatial pattern in the population. While, on the flip side of the coin, it is the exact factor that provides the demanded non-genetic diversity in the tissue. The remarkable characteristic of the introduced model as the stem cells' internal mechanism is that it could control the overall behavior of the population without need for any external factors.
Collapse
Affiliation(s)
- Najme Khorasani
- School of Biological Sciences, Institute for Research in Fundamental Sciences (IPM), Tehran, Iran.
| | - Mehdi Sadeghi
- National Institute of Genetic Engineering and Biotechnology (NIGEB), Tehran, Iran
| |
Collapse
|
5
|
Perkins ML, Gandara L, Crocker J. A synthetic synthesis to explore animal evolution and development. Philos Trans R Soc Lond B Biol Sci 2022; 377:20200517. [PMID: 35634925 PMCID: PMC9149795 DOI: 10.1098/rstb.2020.0517] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Identifying the general principles by which genotypes are converted into phenotypes remains a challenge in the post-genomic era. We still lack a predictive understanding of how genes shape interactions among cells and tissues in response to signalling and environmental cues, and hence how regulatory networks generate the phenotypic variation required for adaptive evolution. Here, we discuss how techniques borrowed from synthetic biology may facilitate a systematic exploration of evolvability across biological scales. Synthetic approaches permit controlled manipulation of both endogenous and fully engineered systems, providing a flexible platform for investigating causal mechanisms in vivo. Combining synthetic approaches with multi-level phenotyping (phenomics) will supply a detailed, quantitative characterization of how internal and external stimuli shape the morphology and behaviour of living organisms. We advocate integrating high-throughput experimental data with mathematical and computational techniques from a variety of disciplines in order to pursue a comprehensive theory of evolution. This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.
Collapse
Affiliation(s)
- Mindy Liu Perkins
- Developmental Biology Unit, European Molecular Biology Laboratory, 69117 Heidelberg, Germany
| | - Lautaro Gandara
- Developmental Biology Unit, European Molecular Biology Laboratory, 69117 Heidelberg, Germany
| | - Justin Crocker
- Developmental Biology Unit, European Molecular Biology Laboratory, 69117 Heidelberg, Germany
| |
Collapse
|
6
|
Ziv N, Brenes LR, Johnson A. Multiple molecular events underlie stochastic switching between 2 heritable cell states in fungi. PLoS Biol 2022; 20:e3001657. [PMID: 35594297 PMCID: PMC9162332 DOI: 10.1371/journal.pbio.3001657] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Revised: 06/02/2022] [Accepted: 05/04/2022] [Indexed: 02/07/2023] Open
Abstract
Eukaryotic transcriptional networks are often large and contain several levels of feedback regulation. Many of these networks have the ability to generate and maintain several distinct transcriptional states across multiple cell divisions and to switch between them. In certain instances, switching between cell states is stochastic, occurring in a small subset of cells of an isogenic population in a seemingly homogenous environment. Given the scarcity and unpredictability of switching in these cases, investigating the determining molecular events is challenging. White-opaque switching in the fungal species Candida albicans is an example of stably inherited cell states that are determined by a complex transcriptional network and can serve as an experimentally accessible model system to study characteristics important for stochastic cell fate switching in eukaryotes. In standard lab media, genetically identical cells maintain their cellular identity (either "white" or "opaque") through thousands of cell divisions, and switching between the states is rare and stochastic. By isolating populations of white or opaque cells, previous studies have elucidated the many differences between the 2 stable cell states and identified a set of transcriptional regulators needed for cell type switching and maintenance of the 2 cell types. Yet, little is known about the molecular events that determine the rare, stochastic switching events that occur in single cells. We use microfluidics combined with fluorescent reporters to directly observe rare switching events between the white and opaque states. We investigate the stochastic nature of switching by beginning with white cells and monitoring the activation of Wor1, a master regulator and marker for the opaque state, in single cells and throughout cell pedigrees. Our results indicate that switching requires 2 stochastic steps; first an event occurs that predisposes a lineage of cells to switch. In the second step, some, but not all, of those predisposed cells rapidly express high levels of Wor1 and commit to the opaque state. To further understand the rapid rise in Wor1, we used a synthetic inducible system in Saccharomyces cerevisiae into which a controllable C. albicans Wor1 and a reporter for its transcriptional control region have been introduced. We document that Wor1 positive autoregulation is highly cooperative (Hill coefficient > 3), leading to rapid activation and producing an "all or none" rather than a graded response. Taken together, our results suggest that reaching a threshold level of a master regulator is sufficient to drive cell type switching in single cells and that an earlier molecular event increases the probability of reaching that threshold in certain small lineages of cells. Quantitative molecular analysis of the white-opaque circuit can serve as a model for the general understanding of complex circuits.
Collapse
Affiliation(s)
- Naomi Ziv
- Department of Microbiology and Immunology, University of California, San Francisco, San Francisco, California, United States of America
- * E-mail: (NZ); (AJ)
| | - Lucas R. Brenes
- Department of Microbiology and Immunology, University of California, San Francisco, San Francisco, California, United States of America
| | - Alexander Johnson
- Department of Microbiology and Immunology, University of California, San Francisco, San Francisco, California, United States of America
- * E-mail: (NZ); (AJ)
| |
Collapse
|
7
|
Zhao N, Liu H, Yan F. Oscillation dynamic mechanism driven by time delays in the competent gene regulatory circuit of B. subtilis. INT J BIOMATH 2021. [DOI: 10.1142/s1793524522500176] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Bacillus subtilis with competent states absorbs DNA and may improve the growth of bacteria by integrating new genetic material. Therefore, it is important to clarify how the genes interact in the circuit so that cells enter into a competent state or return to a vegetative state. The gene regulatory circuit consists of two positive feedback loops and one negative feedback loop. In this paper, a mathematical model is developed by considering transcription time delays to further study dynamic behavior of the B. subtilis competent gene regulatory network. Combined with theoretical calculation and numerical simulation, it is verified that the time delay in indirect transcription inhibition indeed has the effect of inducing the periodic oscillation of the B. subtilis competent system. In addition, some important chemical reaction rates can also regulate system dynamic behavior. However, under the control of time delay, the effects of the important chemical reaction rates have changed significantly. In particular, the time delay can advance critical value of the important chemical reaction rates where vibration occurs and can also weaken or even eliminate the effect of the important chemical reaction rates. These results will help us to analyze the competent state of B. subtilis.
Collapse
Affiliation(s)
- Na Zhao
- Department of Mathematics, Yunnan Normal University, Kunming 650500, P. R. China
- Key Laboratory of Complex System Modeling and Application, for Universities in Yunnan, Kunming 650500, P. R. China
| | - Haihong Liu
- Department of Mathematics, Yunnan Normal University, Kunming 650500, P. R. China
- Key Laboratory of Complex System Modeling and Application, for Universities in Yunnan, Kunming 650500, P. R. China
| | - Fang Yan
- Department of Mathematics, Yunnan Normal University, Kunming 650500, P. R. China
- Key Laboratory of Complex System Modeling and Application, for Universities in Yunnan, Kunming 650500, P. R. China
| |
Collapse
|
8
|
Xiong K, Gerstein M, Masel J. Differences in evolutionary accessibility determine which equally effective regulatory motif evolves to generate pulses. Genetics 2021; 219:6358726. [PMID: 34740240 DOI: 10.1093/genetics/iyab140] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2021] [Accepted: 08/17/2021] [Indexed: 01/02/2023] Open
Abstract
Transcriptional regulatory networks (TRNs) are enriched for certain "motifs." Motif usage is commonly interpreted in adaptationist terms, i.e., that the optimal motif evolves. But certain motifs can also evolve more easily than others. Here, we computationally evolved TRNs to produce a pulse of an effector protein. Two well-known motifs, type 1 incoherent feed-forward loops (I1FFLs) and negative feedback loops (NFBLs), evolved as the primary solutions. The relative rates at which these two motifs evolve depend on selection conditions, but under all conditions, either motif achieves similar performance. I1FFLs generally evolve more often than NFBLs. Selection for a tall pulse favors NFBLs, while selection for a fast response favors I1FFLs. I1FFLs are more evolutionarily accessible early on, before the effector protein evolves high expression; when NFBLs subsequently evolve, they tend to do so from a conjugated I1FFL-NFBL genotype. In the empirical S. cerevisiae TRN, output genes of NFBLs had higher expression levels than those of I1FFLs. These results suggest that evolutionary accessibility, and not relative functionality, shapes which motifs evolve in TRNs, and does so as a function of the expression levels of particular genes.
Collapse
Affiliation(s)
- Kun Xiong
- Department of Molecular and Cellular Biology, University of Arizona, Tucson, AZ 85721, USA.,Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT 06520, USA
| | - Mark Gerstein
- Department of Molecular Biophysics and Biochemistry, Yale University, New Haven, CT 06520, USA.,Program in Computational Biology and Bioinformatics, Yale University, New Haven, CT 06520, USA.,Department of Computer Science, Yale University, New Haven, CT 06520, USA.,Department of Statistics and Data Science, Yale University, New Haven, CT 06520, USA
| | - Joanna Masel
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson,AZ 85721, USA
| |
Collapse
|
9
|
Song YM, Hong H, Kim JK. Universally valid reduction of multiscale stochastic biochemical systems using simple non-elementary propensities. PLoS Comput Biol 2021; 17:e1008952. [PMID: 34662330 PMCID: PMC8562860 DOI: 10.1371/journal.pcbi.1008952] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2021] [Revised: 11/02/2021] [Accepted: 10/04/2021] [Indexed: 11/19/2022] Open
Abstract
Biochemical systems consist of numerous elementary reactions governed by the law of mass action. However, experimentally characterizing all the elementary reactions is nearly impossible. Thus, over a century, their deterministic models that typically contain rapid reversible bindings have been simplified with non-elementary reaction functions (e.g., Michaelis-Menten and Morrison equations). Although the non-elementary reaction functions are derived by applying the quasi-steady-state approximation (QSSA) to deterministic systems, they have also been widely used to derive propensities for stochastic simulations due to computational efficiency and simplicity. However, the validity condition for this heuristic approach has not been identified even for the reversible binding between molecules, such as protein-DNA, enzyme-substrate, and receptor-ligand, which is the basis for living cells. Here, we find that the non-elementary propensities based on the deterministic total QSSA can accurately capture the stochastic dynamics of the reversible binding in general. However, serious errors occur when reactant molecules with similar levels tightly bind, unlike deterministic systems. In that case, the non-elementary propensities distort the stochastic dynamics of a bistable switch in the cell cycle and an oscillator in the circadian clock. Accordingly, we derive alternative non-elementary propensities with the stochastic low-state QSSA, developed in this study. This provides a universally valid framework for simplifying multiscale stochastic biochemical systems with rapid reversible bindings, critical for efficient stochastic simulations of cell signaling and gene regulation. To facilitate the framework, we provide a user-friendly open-source computational package, ASSISTER, that automatically performs the present framework.
Collapse
Affiliation(s)
- Yun Min Song
- Department of Mathematical Sciences, Korea Advanced Institute of Science and Technology, Daejeon, Republic of Korea
- Biomedical Mathematics Group, Institute for Basic Science, Daejeon, Republic of Korea
| | - Hyukpyo Hong
- Department of Mathematical Sciences, Korea Advanced Institute of Science and Technology, Daejeon, Republic of Korea
- Biomedical Mathematics Group, Institute for Basic Science, Daejeon, Republic of Korea
| | - Jae Kyoung Kim
- Department of Mathematical Sciences, Korea Advanced Institute of Science and Technology, Daejeon, Republic of Korea
- Biomedical Mathematics Group, Institute for Basic Science, Daejeon, Republic of Korea
| |
Collapse
|
10
|
Wei L, Li S, Zhang P, Hu T, Zhang MQ, Xie Z, Wang X. Characterizing microRNA-mediated modulation of gene expression noise and its effect on synthetic gene circuits. Cell Rep 2021; 36:109573. [PMID: 34433047 DOI: 10.1016/j.celrep.2021.109573] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2020] [Revised: 07/13/2021] [Accepted: 07/28/2021] [Indexed: 10/20/2022] Open
Abstract
MicroRNAs (miRNAs) have been shown to modulate gene expression noise, but less is known about how miRNAs with different properties may regulate noise differently. Here, we investigate the role of competing RNAs and the composition of miRNA response elements (MREs) in modulating noise. We find that weak competing RNAs could introduce lower noise than strong competing RNAs. In comparison with a single MRE, both repetitive and composite MREs can reduce the noise at low expression, but repetitive MREs can elevate the noise remarkably at high expression. We further observed the behavior of a synthetic cell-type classifier with miRNAs as inputs and find that miRNAs and MREs that could introduce higher noise tend to enhance cell state transition. These results provide a systematic and quantitative understanding of the function of miRNAs in controlling gene expression noise and the utilization of miRNAs to modulate the behavior of synthetic gene circuits.
Collapse
Affiliation(s)
- Lei Wei
- Ministry of Education Key Laboratory of Bioinformatics; Center for Synthetic and Systems Biology; Bioinformatics Division, Beijing National Research Center for Information Science and Technology; Department of Automation, Tsinghua University, Beijing 100084, China
| | - Shuailin Li
- School of Life Sciences, Tsinghua University, Beijing 100084, China
| | - Pengcheng Zhang
- Ministry of Education Key Laboratory of Bioinformatics; Center for Synthetic and Systems Biology; Bioinformatics Division, Beijing National Research Center for Information Science and Technology; Department of Automation, Tsinghua University, Beijing 100084, China
| | - Tao Hu
- Ministry of Education Key Laboratory of Bioinformatics; Center for Synthetic and Systems Biology; Bioinformatics Division, Beijing National Research Center for Information Science and Technology; Department of Automation, Tsinghua University, Beijing 100084, China
| | - Michael Q Zhang
- Ministry of Education Key Laboratory of Bioinformatics; Center for Synthetic and Systems Biology; Bioinformatics Division, Beijing National Research Center for Information Science and Technology; Department of Automation, Tsinghua University, Beijing 100084, China; Department of Basic Medical Sciences, School of Medicine, Tsinghua University, Beijing 100084, China; Department of Molecular and Cell Biology, Center for Systems Biology, The University of Texas, Richardson, TX 75080-3021, USA
| | - Zhen Xie
- Ministry of Education Key Laboratory of Bioinformatics; Center for Synthetic and Systems Biology; Bioinformatics Division, Beijing National Research Center for Information Science and Technology; Department of Automation, Tsinghua University, Beijing 100084, China
| | - Xiaowo Wang
- Ministry of Education Key Laboratory of Bioinformatics; Center for Synthetic and Systems Biology; Bioinformatics Division, Beijing National Research Center for Information Science and Technology; Department of Automation, Tsinghua University, Beijing 100084, China.
| |
Collapse
|
11
|
Alves R, Salvadó B, Milo R, Vilaprinyo E, Sorribas A. Maximization of information transmission influences selection of native phosphorelay architectures. PeerJ 2021; 9:e11558. [PMID: 34178454 PMCID: PMC8199921 DOI: 10.7717/peerj.11558] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2020] [Accepted: 05/12/2021] [Indexed: 01/28/2023] Open
Abstract
Phosphorelays are signal transduction circuits that sense environmental changes and adjust cellular metabolism. Five different circuit architectures account for 99% of all phosphorelay operons annotated in over 9,000 fully sequenced genomes. Here we asked what biological design principles, if any, could explain selection among those architectures in nature. We began by studying kinetically well characterized phosphorelays (Spo0 of Bacillus subtilis and Sln1 of Saccharomyces cerevisiae). We find that natural circuit architecture maximizes information transmission in both cases. We use mathematical models to compare information transmission among the architectures for a realistic range of concentration and parameter values. Mapping experimentally determined phosphorelay protein concentrations onto that range reveals that the native architecture maximizes information transmission in sixteen out of seventeen analyzed phosphorelays. These results suggest that maximization of information transmission is important in the selection of native phosphorelay architectures, parameter values and protein concentrations.
Collapse
Affiliation(s)
- Rui Alves
- Ciències Mèdiques Bàsiques, Universitat de Lleida, Lleida, Spain
| | - Baldiri Salvadó
- Ciències Mèdiques Bàsiques, Universitat de Lleida, Lleida, Spain
| | - Ron Milo
- Plant and Environmental Science, Weizmann Institute of Science, Rehovot, Israel
| | - Ester Vilaprinyo
- Ciències Mèdiques Bàsiques, Universitat de Lleida, Lleida, Spain
| | - Albert Sorribas
- Ciències Mèdiques Bàsiques, Universitat de Lleida, Lleida, Spain
| |
Collapse
|
12
|
Das S, Barik D. Scaling of intrinsic noise in an autocratic reaction network. Phys Rev E 2021; 103:042403. [PMID: 34006004 DOI: 10.1103/physreve.103.042403] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2020] [Accepted: 03/16/2021] [Indexed: 11/07/2022]
Abstract
Biochemical reactions in living cells often produce stochastic trajectories due to the fluctuations of the finite number of the macromolecular species present inside the cell. A significant number of computational and theoretical studies have previously investigated stochasticity in small regulatory networks to understand its origin and regulation. At the systems level regulatory networks have been determined to be hierarchical resembling social networks. In order to determine the stochasticity in networks with hierarchical architecture, here we computationally investigated intrinsic noise in an autocratic reaction network in which only the upstream regulators regulate the downstream regulators. We studied the effects of the qualitative and quantitative nature of regulatory interactions on the stochasticity in the network. We established an unconventional scaling of noise with average abundance in which the noise passes through a minimum indicating that the network can be noisy both in the low and high abundance regimes. We determined that the bursty kinetics of the trajectories are responsible for such scaling. The scaling of noise remains intact for a mixed network that includes democratic subnetworks within the autocratic network.
Collapse
Affiliation(s)
- Soutrick Das
- School of Chemistry, University of Hyderabad, Gachibowli, 500046, Hyderabad, India
| | - Debashis Barik
- School of Chemistry, University of Hyderabad, Gachibowli, 500046, Hyderabad, India
| |
Collapse
|
13
|
Deshmukh S, Saini S. Phenotypic Heterogeneity in Tumor Progression, and Its Possible Role in the Onset of Cancer. Front Genet 2020; 11:604528. [PMID: 33329751 PMCID: PMC7734151 DOI: 10.3389/fgene.2020.604528] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2020] [Accepted: 11/10/2020] [Indexed: 12/20/2022] Open
Abstract
Heterogeneity among isogenic cells/individuals has been known for at least 150 years. Even Mendel, working on pea plants, realized that not all tall plants were identical. However, Mendel was more interested in the discontinuous variation between genetically distinct individuals. The concept of environment dictating distinct phenotypes among isogenic individuals has since been shown to impact the evolution of populations in numerous examples at different scales of life. In this review, we discuss how phenotypic heterogeneity and its evolutionary implications exist at all levels of life, from viruses to mammals. In particular, we discuss how a particular disease condition (cancer) is impacted by heterogeneity among isogenic cells, and propose a potential role that phenotypic heterogeneity might play toward the onset of the disease.
Collapse
Affiliation(s)
- Saniya Deshmukh
- Department of Chemical Engineering, Indian Institute of Technology Bombay, Mumbai, India
| | - Supreet Saini
- Department of Chemical Engineering, Indian Institute of Technology Bombay, Mumbai, India
| |
Collapse
|
14
|
Flux variability analysis reveals a tragedy of commons in cancer cells. SN APPLIED SCIENCES 2020. [DOI: 10.1007/s42452-020-03762-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022] Open
|
15
|
Guinn MT, Wan Y, Levovitz S, Yang D, Rosner MR, Balázsi G. Observation and Control of Gene Expression Noise: Barrier Crossing Analogies Between Drug Resistance and Metastasis. Front Genet 2020; 11:586726. [PMID: 33193723 PMCID: PMC7662081 DOI: 10.3389/fgene.2020.586726] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2020] [Accepted: 09/29/2020] [Indexed: 12/12/2022] Open
Affiliation(s)
- Michael Tyler Guinn
- Biomedical Engineering Department, Stony Brook University, Stony Brook, NY, United States.,Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, United States.,Stony Brook Medical Scientist Training Program, Stony Brook, NY, United States
| | - Yiming Wan
- Biomedical Engineering Department, Stony Brook University, Stony Brook, NY, United States.,Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, United States
| | - Sarah Levovitz
- Biomedical Engineering Department, Stony Brook University, Stony Brook, NY, United States.,Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, United States
| | - Dongbo Yang
- Ben May Department for Cancer Research, The University of Chicago, Chicago, IL, United States
| | - Marsha R Rosner
- Ben May Department for Cancer Research, The University of Chicago, Chicago, IL, United States
| | - Gábor Balázsi
- Biomedical Engineering Department, Stony Brook University, Stony Brook, NY, United States.,Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, United States
| |
Collapse
|
16
|
Synthetic gene-regulatory networks in the opportunistic human pathogen Streptococcus pneumoniae. Proc Natl Acad Sci U S A 2020; 117:27608-27619. [PMID: 33087560 PMCID: PMC7959565 DOI: 10.1073/pnas.1920015117] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Streptococcus pneumoniae is a major human pathogen responsible for enormous global morbidity and mortality. Despite this, the pneumococcus makes up part of the commensal nasopharyngeal flora. How the pneumococcus switches from this commensal to pathogenic state and causes disease is unclear and very likely involves variability in expression of its virulence factors. Here, we used synthetic biology approaches to generate complex gene-regulatory networks such as logic gates and toggle switches. We show that these networks are functional in vivo to control capsule production in an influenza-superinfection model. This opens the field of systematically testing the role of phenotypic variation in pneumococcal virulence. The approaches used here may serve as an example for synthetic biology projects in unrelated organisms. Streptococcus pneumoniae can cause disease in various human tissues and organs, including the ear, the brain, the blood, and the lung, and thus in highly diverse and dynamic environments. It is challenging to study how pneumococci control virulence factor expression, because cues of natural environments and the presence of an immune system are difficult to simulate in vitro. Here, we apply synthetic biology methods to reverse-engineer gene expression control in S. pneumoniae. A selection platform is described that allows for straightforward identification of transcriptional regulatory elements out of combinatorial libraries. We present TetR- and LacI-regulated promoters that show expression ranges of four orders of magnitude. Based on these promoters, regulatory networks of higher complexity are assembled, such as logic AND gates and IMPLY gates. We demonstrate single-copy genome-integrated toggle switches that give rise to bimodal population distributions. The tools described here can be used to mimic complex expression patterns, such as the ones found for pneumococcal virulence factors. Indeed, we were able to rewire gene expression of the capsule operon, the main pneumococcal virulence factor, to be externally inducible (YES gate) or to act as an IMPLY gate (only expressed in absence of inducer). Importantly, we demonstrate that these synthetic gene-regulatory networks are functional in an influenza A virus superinfection murine model of pneumonia, paving the way for in vivo investigations of the importance of gene expression control on the pathogenicity of S. pneumoniae.
Collapse
|
17
|
Khorasani N, Sadeghi M, Nowzari-Dalini A. A computational model of stem cell molecular mechanism to maintain tissue homeostasis. PLoS One 2020; 15:e0236519. [PMID: 32730297 PMCID: PMC7392222 DOI: 10.1371/journal.pone.0236519] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2020] [Accepted: 07/07/2020] [Indexed: 11/24/2022] Open
Abstract
Stem cells, with their capacity to self-renew and to differentiate to more specialized cell types, play a key role to maintain homeostasis in adult tissues. To investigate how, in the dynamic stochastic environment of a tissue, non-genetic diversity and the precise balance between proliferation and differentiation are achieved, it is necessary to understand the molecular mechanisms of the stem cells in decision making process. By focusing on the impact of stochasticity, we proposed a computational model describing the regulatory circuitry as a tri-stable dynamical system to reveal the mechanism which orchestrate this balance. Our model explains how the distribution of noise in genes, linked to the cell regulatory networks, affects cell decision-making to maintain homeostatic state. The noise effect on tissue homeostasis is achieved by regulating the probability of differentiation and self-renewal through symmetric and/or asymmetric cell divisions. Our model reveals, when mutations due to the replication of DNA in stem cell division, are inevitable, how mutations contribute to either aging gradually or the development of cancer in a short period of time. Furthermore, our model sheds some light on the impact of more complex regulatory networks on the system robustness against perturbations.
Collapse
Affiliation(s)
- Najme Khorasani
- School of Mathematics, Statistics and Computer Science, College of Science, University of Tehran, Tehran, Iran
| | - Mehdi Sadeghi
- National Institute of Genetic Engineering and Biotechnology (NIGEB), Tehran, Iran.,School of Biological Sciences, Institute for Research in Fundamental Sciences (IPM), Tehran, Iran
| | - Abbas Nowzari-Dalini
- School of Mathematics, Statistics and Computer Science, College of Science, University of Tehran, Tehran, Iran
| |
Collapse
|
18
|
Das S, Barik D. Qualitative and quantitative nature of mutual interactions dictate chemical noise in a democratic reaction network. Phys Rev E 2020; 101:042407. [PMID: 32422814 DOI: 10.1103/physreve.101.042407] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2019] [Accepted: 03/16/2020] [Indexed: 06/11/2023]
Abstract
The functions of a living cell rely on a complex network of biochemical reactions that allow it to respond against various internal and external cues. The outcomes of these chemical reactions are often stochastic due to intrinsic and extrinsic noise leading to population heterogeneity. The majority of calculations of stochasticity in reaction networks have focused on small regulatory networks addressing the role of timescales, feedback regulations, and network topology in propagation of noise. Here we computationally investigated chemical noise in a network with democratic architecture where each node is regulated by all other nodes in the network. We studied the effects of the qualitative and quantitative nature of mutual interactions on the propagation of both intrinsic and extrinsic noise in the network. We show that an increased number of inhibitory signals lead to ultrasensitive switching of average and that leads to sharp transition of intrinsic noise. The intrinsic noise exhibits a biphasic power-law scaling with the average, and the scaling coefficients strongly correlate with the strength of inhibitory signal. The noise strength critically depends on the strength of the interactions, where negative interactions attenuate both intrinsic and extrinsic noise.
Collapse
Affiliation(s)
- Soutrick Das
- School of Chemistry, University of Hyderabad, Gachibowli, 500046 Hyderabad, India
| | - Debashis Barik
- School of Chemistry, University of Hyderabad, Gachibowli, 500046 Hyderabad, India
| |
Collapse
|
19
|
Dynamic Analysis of the Time-Delayed Genetic Regulatory Network Between Two Auto-Regulated and Mutually Inhibitory Genes. Bull Math Biol 2020; 82:46. [DOI: 10.1007/s11538-020-00722-1] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2019] [Accepted: 03/16/2020] [Indexed: 01/14/2023]
|
20
|
Rodrigo G. Insights about collective decision-making at the genetic level. Biophys Rev 2019; 12:19-24. [PMID: 31845181 DOI: 10.1007/s12551-019-00608-0] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2019] [Accepted: 12/05/2019] [Indexed: 01/08/2023] Open
Abstract
By living in a collective, individuals can share and aggregate information to base their decisions on the many rather than on the one, thereby increasing accuracy. But a collective can also be defined at the molecular level. In the following, we reason that genes, by working collectively, share fundamental features with social organisms, which ends, without invoking cognition, in wiser responses. For that, we compile into a single picture the terms redundancy, stochastic resonance, intrinsic and extrinsic noise, and cross-regulation.
Collapse
Affiliation(s)
- Guillermo Rodrigo
- Institute for Integrative Systems Biology (I2SysBio), CSIC - U. Valencia, 46980, Paterna, Spain.
| |
Collapse
|
21
|
Guillemin A, Duchesne R, Crauste F, Gonin-Giraud S, Gandrillon O. Drugs modulating stochastic gene expression affect the erythroid differentiation process. PLoS One 2019; 14:e0225166. [PMID: 31751364 PMCID: PMC6872177 DOI: 10.1371/journal.pone.0225166] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2019] [Accepted: 10/30/2019] [Indexed: 12/30/2022] Open
Abstract
To better understand the mechanisms behind cells decision-making to differentiate, we assessed the influence of stochastic gene expression (SGE) modulation on the erythroid differentiation process. It has been suggested that stochastic gene expression has a role in cell fate decision-making which is revealed by single-cell analyses but studies dedicated to demonstrate the consistency of this link are still lacking. Recent observations showed that SGE significantly increased during differentiation and a few showed that an increase of the level of SGE is accompanied by an increase in the differentiation process. However, a consistent relation in both increasing and decreasing directions has never been shown in the same cellular system. Such demonstration would require to be able to experimentally manipulate simultaneously the level of SGE and cell differentiation in order to observe if cell behavior matches with the current theory. We identified three drugs that modulate SGE in primary erythroid progenitor cells. Both Artemisinin and Indomethacin decreased SGE and reduced the amount of differentiated cells. On the contrary, a third component called MB-3 simultaneously increased the level of SGE and the amount of differentiated cells. We then used a dynamical modelling approach which confirmed that differentiation rates were indeed affected by the drug treatment. Using single-cell analysis and modeling tools, we provide experimental evidence that, in a physiologically relevant cellular system, SGE is linked to differentiation.
Collapse
Affiliation(s)
- Anissa Guillemin
- Laboratoire de biologie et modélisation de la cellule. LBMC - Ecole Normale Supérieure - Lyon, Université Claude Bernard Lyon 1, Centre National de la Recherche Scientifique: UMR5239, Institut National de la Santé et de la Recherche Médicale: U1210 - Ecole Normale Supérieure de Lyon 46 allée d’Italie 69007 Lyon, France
| | - Ronan Duchesne
- Laboratoire de biologie et modélisation de la cellule. LBMC - Ecole Normale Supérieure - Lyon, Université Claude Bernard Lyon 1, Centre National de la Recherche Scientifique: UMR5239, Institut National de la Santé et de la Recherche Médicale: U1210 - Ecole Normale Supérieure de Lyon 46 allée d’Italie 69007 Lyon, France
- Inria Dracula, Villeurbanne, France
| | - Fabien Crauste
- Inria Dracula, Villeurbanne, France
- Univ. Bordeaux, CNRS, Bordeaux INP, IMB, UMR 5251, F-33400, Talence, France
| | - Sandrine Gonin-Giraud
- Laboratoire de biologie et modélisation de la cellule. LBMC - Ecole Normale Supérieure - Lyon, Université Claude Bernard Lyon 1, Centre National de la Recherche Scientifique: UMR5239, Institut National de la Santé et de la Recherche Médicale: U1210 - Ecole Normale Supérieure de Lyon 46 allée d’Italie 69007 Lyon, France
| | - Olivier Gandrillon
- Laboratoire de biologie et modélisation de la cellule. LBMC - Ecole Normale Supérieure - Lyon, Université Claude Bernard Lyon 1, Centre National de la Recherche Scientifique: UMR5239, Institut National de la Santé et de la Recherche Médicale: U1210 - Ecole Normale Supérieure de Lyon 46 allée d’Italie 69007 Lyon, France
- Inria Dracula, Villeurbanne, France
| |
Collapse
|
22
|
VÁZQUEZ-JIMÉNEZ AARÓN, SANTILLÁN MOISÉS, RODRÍGUEZ-GONZÁLEZ JESÚS. CHARACTERIZATION OF INTRINSIC AND EXTRINSIC NOISE EFFECTS IN POSITIVELY REGULATED GENES. J BIOL SYST 2019. [DOI: 10.1142/s0218339019500165] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Gene regulation is fundamental for cell survival. This regulation must be both robust to noise and sensitive enough to external stimuli to elicit the proper responses. In this work, we study, through stochastic numerical simulations, how a gene regulatory network with a positive feedback loop responds to environmental changes in the presence of intrinsic and extrinsic noises. Noise effects were characterized by measuring the statistical differences between two protein time series resulting from identical systems subject to the same source of extrinsic noise. A robust analysis was implemented by modifying the kinetic system parameters. We found that the common source of time-varying extrinsic fluctuations leads to a correlation in the systems it affects. The correlation and the extrinsic and intrinsic noise components are modulated by the update period and noise intensity parameters. Our results suggest that noise perception is controlled through the parameters associated with the response time: degradation rates and promoter dissociation constant.
Collapse
Affiliation(s)
- AARÓN VÁZQUEZ-JIMÉNEZ
- Centro de Investigación y Estudios Avanzados del IPN, Unidad Monterrey, Vía del Conocimiento 201, Parque PIIT, CP 66600 Apodaca NL, México
| | - MOISÉS SANTILLÁN
- Centro de Investigación y Estudios Avanzados del IPN, Unidad Monterrey, Vía del Conocimiento 201, Parque PIIT, CP 66600 Apodaca NL, México
| | - JESÚS RODRÍGUEZ-GONZÁLEZ
- Centro de Investigación y Estudios Avanzados del IPN, Unidad Monterrey, Vía del Conocimiento 201, Parque PIIT, CP 66600 Apodaca NL, México
| |
Collapse
|
23
|
Yan J, Hilfinger A, Vinnicombe G, Paulsson J. Kinetic Uncertainty Relations for the Control of Stochastic Reaction Networks. PHYSICAL REVIEW LETTERS 2019; 123:108101. [PMID: 31573304 DOI: 10.1103/physrevlett.123.108101] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/23/2019] [Revised: 06/27/2019] [Indexed: 06/10/2023]
Abstract
Nonequilibrium stochastic reaction networks are commonly found in both biological and nonbiological systems, but have remained hard to analyze because small differences in rate functions or topology can change the dynamics drastically. Here, we conjecture exact quantitative inequalities that relate the extent of fluctuations in connected components, for various network topologies. Specifically, we find that regardless of how two components affect each other's production rates, it is impossible to suppress fluctuations below the uncontrolled equivalents for both components: one must increase its fluctuations for the other to be suppressed. For systems in which components control each other in ringlike structures, it appears that fluctuations can only be suppressed in one component if all other components instead increase fluctuations, compared to the case without control. Even the general N-component system-with arbitrary connections and parameters-must have at least one component with increased fluctuations to reduce fluctuations in others. In connected reaction networks it thus appears impossible to reduce the statistical uncertainty in all components, regardless of the control mechanisms or energy dissipation.
Collapse
Affiliation(s)
- Jiawei Yan
- Department of Systems Biology, Harvard University, 200 Longwood Avenue, Boston, Massachusetts 02115, USA
- Department of Molecular and Cellular Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Andreas Hilfinger
- Department of Chemical & Physical Sciences, University of Toronto, Mississauga, Ontario L5L 1C6, Canada
- Department of Mathematics, University of Toronto, 40 St. George Street, Toronto, Ontario M5S 2E4, Canada
- Department of Cell & Systems Biology, University of Toronto, 25 Harbord Street, Toronto, Ontario M5S 3G5, Canada
| | - Glenn Vinnicombe
- Department of Engineering, University of Cambridge, Cambridge CB2 1PZ, United Kingdom
| | - Johan Paulsson
- Department of Systems Biology, Harvard University, 200 Longwood Avenue, Boston, Massachusetts 02115, USA
| |
Collapse
|
24
|
Farquhar KS, Charlebois DA, Szenk M, Cohen J, Nevozhay D, Balázsi G. Role of network-mediated stochasticity in mammalian drug resistance. Nat Commun 2019; 10:2766. [PMID: 31235692 PMCID: PMC6591227 DOI: 10.1038/s41467-019-10330-w] [Citation(s) in RCA: 63] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2018] [Accepted: 05/03/2019] [Indexed: 11/11/2022] Open
Abstract
A major challenge in biology is that genetically identical cells in the same environment can display gene expression stochasticity (noise), which contributes to bet-hedging, drug tolerance, and cell-fate switching. The magnitude and timescales of stochastic fluctuations can depend on the gene regulatory network. Currently, it is unclear how gene expression noise of specific networks impacts the evolution of drug resistance in mammalian cells. Answering this question requires adjusting network noise independently from mean expression. Here, we develop positive and negative feedback-based synthetic gene circuits to decouple noise from the mean for Puromycin resistance gene expression in Chinese Hamster Ovary cells. In low Puromycin concentrations, the high-noise, positive-feedback network delays long-term adaptation, whereas it facilitates adaptation under high Puromycin concentration. Accordingly, the low-noise, negative-feedback circuit can maintain resistance by acquiring mutations while the positive-feedback circuit remains mutation-free and regains drug sensitivity. These findings may have profound implications for chemotherapeutic inefficiency and cancer relapse. The role of gene expression noise in the evolution of drug resistance in mammalian cells is unclear. Here, by uncoupling noise from mean expression of a drug resistance gene in CHO cells the authors show that noisy expression aids adaptation to high drug levels, but delays it at low drug levels.
Collapse
Affiliation(s)
- Kevin S Farquhar
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, 11794, USA.,Genetics and Epigenetics Graduate Program, The University of Texas MD Anderson Cancer Center, UT Health Graduate School of Biomedical Sciences, Houston, TX, 77030, USA
| | - Daniel A Charlebois
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, 11794, USA.,Department of Physics, University of Alberta, Edmonton, AB, 4-181 CCIS, T6G-2E1, Canada
| | - Mariola Szenk
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, 11794, USA.,Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, 11794, USA
| | - Joseph Cohen
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, 11794, USA
| | - Dmitry Nevozhay
- School of Biomedicine, Far Eastern Federal University, 8 Sukhanova Street, Vladivostok, 690950, Russia.,Department of Systems Biology, The University of Texas MD Anderson Cancer Center, Houston, TX, 77030, USA
| | - Gábor Balázsi
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, 11794, USA. .,Genetics and Epigenetics Graduate Program, The University of Texas MD Anderson Cancer Center, UT Health Graduate School of Biomedical Sciences, Houston, TX, 77030, USA. .,Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, 11794, USA. .,Department of Systems Biology, The University of Texas MD Anderson Cancer Center, Houston, TX, 77030, USA.
| |
Collapse
|
25
|
Firman T, Amgalan A, Ghosh K. Maximum Caliber Can Build and Infer Models of Oscillation in a Three-Gene Feedback Network. J Phys Chem B 2019; 123:343-355. [PMID: 30507199 DOI: 10.1021/acs.jpcb.8b07465] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Single-cell protein expression time trajectories provide rich temporal data quantifying cellular variability and its role in dictating fitness. However, theoretical models to analyze and fully extract information from these measurements remain limited for three reasons: (i) gene expression profiles are noisy, rendering models of averages inapplicable, (ii) experiments typically measure only a few protein species while leaving other molecular actors-necessary to build traditional bottom-up models-unnoticed, and (iii) measured data are in fluorescence, not particle number. We recently addressed these challenges in an alternate top-down approach using the principle of Maximum Caliber (MaxCal) to model genetic switches with one and two protein species. In the present work we address scalability and broader applicability of MaxCal by extending to a three-gene (A, B, C) feedback network that exhibits oscillation, commonly known as the repressilator. We test MaxCal's inferential power by using synthetic data of noisy protein number time traces-serving as a proxy for experimental data-generated from a known underlying model. We notice that the minimal MaxCal model-accounting for production, degradation, and only one type of symmetric coupling between all three species-reasonably infers several underlying features of the circuit such as the effective production rate, degradation rate, frequency of oscillation, and protein number distribution. Next, we build models of higher complexity including different levels of coupling between A, B, and C and rigorously assess their relative performance. While the minimal model (with four parameters) performs remarkably well, we note that the most complex model (with six parameters) allowing all possible forms of crosstalk between A, B, and C slightly improves prediction of rates, but avoids ad hoc assumption of all the other models. It is also the model of choice based on Bayesian information criteria. We further analyzed time trajectories in arbitrary fluorescence (using synthetic trajectories) to mimic realistic data. We conclude that even with a three-protein system including both fluorescence noise and intrinsic gene expression fluctuations, MaxCal can faithfully infer underlying details of the network, opening future directions to model other network motifs with many species.
Collapse
|
26
|
Firman T, Amgalan A, Ghosh K. Maximum Caliber Can Build and Infer Models of Oscillation in a Three-Gene Feedback Network. J Phys Chem A 2018. [DOI: 10.1021/acs.jpca.8b07465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
27
|
Arbel-Goren R, Di Patti F, Fanelli D, Stavans J. Noise⁻Seeded Developmental Pattern Formation in Filamentous Cyanobacteria. Life (Basel) 2018; 8:life8040058. [PMID: 30423937 PMCID: PMC6316479 DOI: 10.3390/life8040058] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2018] [Revised: 10/24/2018] [Accepted: 11/05/2018] [Indexed: 11/16/2022] Open
Abstract
Under nitrogen-poor conditions, multicellular cyanobacteria such as Anabaena sp. PCC 7120 undergo a process of differentiation, forming nearly regular, developmental patterns of individual nitrogen-fixing cells, called heterocysts, interspersed between intervals of vegetative cells that carry out photosynthesis. Developmental pattern formation is mediated by morphogen species that can act as activators and inhibitors, some of which can diffuse along filaments. We survey the limitations of the classical, deterministic Turing mechanism that has been often invoked to explain pattern formation in these systems, and then, focusing on a simpler system governed by birth-death processes, we illustrate pedagogically a recently proposed paradigm that provides a much more robust description of pattern formation: stochastic Turing patterns. We emphasize the essential role that cell-to-cell differences in molecular numbers—caused by inevitable fluctuations in gene expression—play, the so called demographic noise, in seeding the formation of stochastic Turing patterns over a much larger region of parameter space, compared to their deterministic counterparts.
Collapse
Affiliation(s)
- Rinat Arbel-Goren
- Department of Physics of Complex Systems, Weizmann Institute of Science, Rehovot 7610001, Israel.
| | - Francesca Di Patti
- Consorzio Interuniversitario Nazionale per la Scienza e Tecnologia dei Materiali, Dip. di Chimica, Università degli Studi di Firenze, Via della Lastruccia 3-13, 50019 Sesto Fiorentino, Firenze, Italy.
- Istituto Nazionale di Fisica Nucleare, Sezione di Firenze, via G. Sansone 1, 50019 Sesto Fiorentino, Firenze, Italy.
- Centro Interdipartimentale per lo Studio delle Dinamiche Complesse, via G. Sansone 1, 50019 Sesto Fiorentino, Firenze, Italy.
| | - Duccio Fanelli
- Istituto Nazionale di Fisica Nucleare, Sezione di Firenze, via G. Sansone 1, 50019 Sesto Fiorentino, Firenze, Italy.
- Centro Interdipartimentale per lo Studio delle Dinamiche Complesse, via G. Sansone 1, 50019 Sesto Fiorentino, Firenze, Italy.
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, 50019 Sesto Fiorentino, Firenze, Italy.
| | - Joel Stavans
- Department of Physics of Complex Systems, Weizmann Institute of Science, Rehovot 7610001, Israel.
| |
Collapse
|
28
|
Role of noise and parametric variation in the dynamics of gene regulatory circuits. NPJ Syst Biol Appl 2018; 4:40. [PMID: 30416751 PMCID: PMC6218471 DOI: 10.1038/s41540-018-0076-x] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2018] [Revised: 10/14/2018] [Accepted: 10/16/2018] [Indexed: 12/21/2022] Open
Abstract
Stochasticity in gene expression impacts the dynamics and functions of gene regulatory circuits. Intrinsic noises, including those that are caused by low copy number of molecules and transcriptional bursting, are usually studied by stochastic simulations. However, the role of extrinsic factors, such as cell-to-cell variability and heterogeneity in the microenvironment, is still elusive. To evaluate the effects of both the intrinsic and extrinsic noises, we develop a method, named sRACIPE, by integrating stochastic analysis with random circuit perturbation (RACIPE) method. RACIPE uniquely generates and analyzes an ensemble of models with random kinetic parameters. Previously, we have shown that the gene expression from random models form robust and functionally related clusters. In sRACIPE we further develop two stochastic simulation schemes, aiming to reduce the computational cost without sacrificing the convergence of statistics. One scheme uses constant noise to capture the basins of attraction, and the other one uses simulated annealing to detect the stability of states. By testing the methods on several synthetic gene regulatory circuits and an epithelial-mesenchymal transition network in squamous cell carcinoma, we demonstrate that sRACIPE can interpret the experimental observations from single-cell gene expression data. We observe that parametric variation (the spread of parameters around a median value) increases the spread of the gene expression clusters, whereas high noise merges the states. Our approach quantifies the robustness of a gene circuit in the presence of noise and sheds light on a new mechanism of noise-induced hybrid states. We have implemented sRACIPE as an R package.
Collapse
|
29
|
Hao L, Yang Z, Turcotte M. Time-scale separation and stochasticity conspire to impact phenotypic dynamics in the canonical and inverted Bacillus subtilis core genetic regulation circuits. QUANTITATIVE BIOLOGY 2018. [DOI: 10.1007/s40484-018-0151-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
30
|
Álvarez-Buylla Roces ME, Martínez-García JC, Dávila-Velderrain J, Domínguez-Hüttinger E, Martínez-Sánchez ME. Medical Systems Biology. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2018; 1069:1-33. [PMID: 30076565 DOI: 10.1007/978-3-319-89354-9_1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
The aim of this volume is to encourage the use of systems-level methodologies to contribute to the improvement of human-health . We intend to motivate biomedical researchers to complement their current theoretical and empirical practice with up-to-date systems biology conceptual approaches. Our perspective is based on the deep understanding of the key biomolecular regulatory mechanisms that underlie health, as well as the emergence and progression of human-disease . We strongly believe that the contemporary systems biology perspective opens the door to the effective development of novel methodologies to the improvement of prevention . This requires a deeper and integrative understanding of the involved underlying systems-level mechanisms. In order to explain our proposal in a simple way, in this chapter we privilege the conceptual exposition of our chosen framework over formal considerations. The formal exposition of our proposal will be expanded and discussed later in the next chapters.
Collapse
|
31
|
Wu F, Chen X, Zheng Y, Duan J, Kurths J, Li X. Lévy noise induced transition and enhanced stability in a gene regulatory network. CHAOS (WOODBURY, N.Y.) 2018; 28:075510. [PMID: 30070513 DOI: 10.1063/1.5025235] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/08/2018] [Accepted: 05/09/2018] [Indexed: 06/08/2023]
Abstract
We investigate a quantitative bistable two-dimensional model (MeKS network) of gene expression dynamics describing the competence development in the Bacillus subtilis under the influence of Lévy as well as Brownian motions. To analyze the transitions between the vegetative and the competence regions therein, two dimensionless deterministic quantities, the mean first exit time (MFET) and the first escape probability, are determined from a microscopic perspective, as well as their averaged versions from a macroscopic perspective. The relative contribution factor λ, the ratio of non-Gaussian and Gaussian noise strengths, is adopted to identify an optimum choice in these transitions. Additionally, we use a recent geometric concept, the stochastic basin of attraction (SBA), to exhibit a pictorial comprehension about the influence of the Lévy motion on the basin stability of the competence state. Our main results indicate that (i) the transitions between the vegetative and the competence regions can be induced by the noise intensities, the relative contribution factor λ and the Lévy motion index α; (ii) a higher noise intensity and a larger α with smaller jump magnitude make the MFET shorter, and the MFET as a function of λ exhibits one maximum value, which is a signature of the noise-enhanced stability phenomenon for the vegetative state; (iii) a larger α makes the transition from the vegetative to the adjacent competence region to occur at the highest probability. The Lévy motion index α0≈0.5 (a larger jump magnitude with a lower frequency) is an ideal choice to implement the transition to the non-adjacent competence region; (iv) there is an expansion in SBA when α decreases.
Collapse
Affiliation(s)
- Fengyan Wu
- Center for Mathematical Sciences & School of Mathematics and Statistics, Huazhong University of Science and Technology, Wuhan 430074, China
| | - Xiaoli Chen
- Center for Mathematical Sciences & School of Mathematics and Statistics, Huazhong University of Science and Technology, Wuhan 430074, China
| | - Yayun Zheng
- Center for Mathematical Sciences & School of Mathematics and Statistics, Huazhong University of Science and Technology, Wuhan 430074, China
| | - Jinqiao Duan
- Wuhan National Laboratory for Optoelectronics, Wuhan 430074, China
| | - Jürgen Kurths
- Department of Applied Mathematics, Illinois Institute of Technology, Chicago, Illinois 60616, USA
| | - Xiaofan Li
- Wuhan National Laboratory for Optoelectronics, Wuhan 430074, China
| |
Collapse
|
32
|
Rosenthal AZ, Qi Y, Hormoz S, Park J, Li SHJ, Elowitz MB. Metabolic interactions between dynamic bacterial subpopulations. eLife 2018; 7:33099. [PMID: 29809139 PMCID: PMC6025961 DOI: 10.7554/elife.33099] [Citation(s) in RCA: 51] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2017] [Accepted: 05/21/2018] [Indexed: 01/08/2023] Open
Abstract
Individual microbial species are known to occupy distinct metabolic niches within multi-species communities. However, it has remained largely unclear whether metabolic specialization can similarly occur within a clonal bacterial population. More specifically, it is not clear what functions such specialization could provide and how specialization could be coordinated dynamically. Here, we show that exponentially growing Bacillus subtilis cultures divide into distinct interacting metabolic subpopulations, including one population that produces acetate, and another population that differentially expresses metabolic genes for the production of acetoin, a pH-neutral storage molecule. These subpopulations exhibit distinct growth rates and dynamic interconversion between states. Furthermore, acetate concentration influences the relative sizes of the different subpopulations. These results show that clonal populations can use metabolic specialization to control the environment through a process of dynamic, environmentally-sensitive state-switching. The chemical reactions that occur within a living organism are collectively referred to as its metabolism. Many metabolic reactions produce byproducts that will poison the cells if they are not dealt with: fermenting bacteria, for example, release harmful organic acids and alcohols. How the bacteria respond to these toxins has been most studied at the level of entire microbial populations, meaning the activities of individual cells are effectively “averaged” together. Yet, even two bacteria with the same genes and living in the same environment can behave in different ways. This raises the question: do bacterial populations specialize into distinct subpopulations that play distinct roles when dealing with metabolic products, or do all cells in the community act in unison? Rosenthal et al. set out to answer this question for a community of Bacillus subtilis, a bacterium that is commonly studied in the laboratory and used for the industrial production of enzymes. The analysis focused on genes involved in fundamental metabolic processes, known as the TCA cycle, which the bacteria use to generate energy and build biomass. The experiments revealed that, even when all the cells are genetically identical, different Bacillus subtilis cells do indeed specialize into metabolic subpopulations with distinct growth rates. Time-lapse movies of bacteria that made fluorescent markers of different colors whenever certain metabolic genes became active showed cells switching different colors on and off, indicating that they switch between metabolic subpopulations. Further biochemical studies and measures of gene activity revealed that the different subpopulations produce and release distinct metabolic products, including toxic byproducts. Notably, the release of these metabolites by one subpopulation appeared to activate other subpopulations within the community. This example of cells specializing into unique interacting metabolic subpopulations provides insight into several fundamental issues in microbiology and beyond. It is relevant to evolutionary biologists, since the fact that fractions of the population can switch in and out of a metabolic state, instead of evolving into several inflexible specialists, may provide an evolutionary advantage in fluctuating natural environments by reducing the risk of extinction. It also has implications for industrial fermentation processes and metabolic engineering, and may help biotechnologists design more efficient ways to harness bacterial metabolism to produce useful products.
Collapse
Affiliation(s)
- Adam Z Rosenthal
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, United States.,Department of Applied Physics, California Institute of Technology, Pasadena, United States
| | - Yutao Qi
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, United States.,Department of Applied Physics, California Institute of Technology, Pasadena, United States
| | - Sahand Hormoz
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, United States.,Department of Applied Physics, California Institute of Technology, Pasadena, United States
| | - Jin Park
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, United States.,Department of Applied Physics, California Institute of Technology, Pasadena, United States
| | - Sophia Hsin-Jung Li
- Department of Molecular Biology, Princeton University, Princeton, United States
| | - Michael B Elowitz
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, United States.,Department of Applied Physics, California Institute of Technology, Pasadena, United States.,Howard Hughes Medical Institute, Pasadena, United States
| |
Collapse
|
33
|
Abstract
Engineering synthetic gene regulatory circuits proceeds through iterative cycles of design, building, and testing. Initial circuit designs must rely on often-incomplete models of regulation established by fields of reductive inquiry—biochemistry and molecular and systems biology. As differences in designed and experimentally observed circuit behavior are inevitably encountered, investigated, and resolved, each turn of the engineering cycle can force a resynthesis in understanding of natural network function. Here, we outline research that uses the process of gene circuit engineering to advance biological discovery. Synthetic gene circuit engineering research has not only refined our understanding of cellular regulation but furnished biologists with a toolkit that can be directed at natural systems to exact precision manipulation of network structure. As we discuss, using circuit engineering to predictively reorganize, rewire, and reconstruct cellular regulation serves as the ultimate means of testing and understanding how cellular phenotype emerges from systems-level network function.
Collapse
Affiliation(s)
- Caleb J. Bashor
- Institute for Medical Engineering and Science, Department of Biological Engineering, and Synthetic Biology Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA;,
| | - James J. Collins
- Institute for Medical Engineering and Science, Department of Biological Engineering, and Synthetic Biology Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA;,
- Harvard–MIT Program in Health Sciences and Technology, Cambridge, Massachusetts 02139, USA
- Broad Institute of MIT and Harvard, Cambridge, Massachusetts 02142, USA
- Wyss Institute for Biologically Inspired Engineering, Harvard University, Boston, Massachusetts 02115, USA
| |
Collapse
|
34
|
Wierenga H, Ten Wolde PR, Becker NB. Quantifying fluctuations in reversible enzymatic cycles and clocks. Phys Rev E 2018; 97:042404. [PMID: 29758603 DOI: 10.1103/physreve.97.042404] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2018] [Indexed: 11/07/2022]
Abstract
Biochemical reactions are fundamentally noisy at a molecular scale. This limits the precision of reaction networks, but it also allows fluctuation measurements that may reveal the structure and dynamics of the underlying biochemical network. Here, we study nonequilibrium reaction cycles, such as the mechanochemical cycle of molecular motors, the phosphorylation cycle of circadian clock proteins, or the transition state cycle of enzymes. Fluctuations in such cycles may be measured using either of two classical definitions of the randomness parameter, which we show to be equivalent in general microscopically reversible cycles. We define a stochastic period for reversible cycles and present analytical solutions for its moments. Furthermore, we associate the two forms of the randomness parameter with the thermodynamic uncertainty relation, which sets limits on the timing precision of the cycle in terms of thermodynamic quantities. Our results should prove useful also for the study of temporal fluctuations in more general networks.
Collapse
Affiliation(s)
| | | | - Nils B Becker
- DKFZ, Bioquant, Im Neuenheimer Feld 267, 69120 Heidelberg, Germany
| |
Collapse
|
35
|
Hortsch SK, Kremling A. Characterization of noise in multistable genetic circuits reveals ways to modulate heterogeneity. PLoS One 2018; 13:e0194779. [PMID: 29579101 PMCID: PMC5868838 DOI: 10.1371/journal.pone.0194779] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2017] [Accepted: 03/11/2018] [Indexed: 02/04/2023] Open
Abstract
Random fluctuations in the amount of cellular components like mRNA and protein molecules are inevitable due to the stochastic and discrete nature of biochemical reactions. If large enough, this so-called “cellular noise” can lead to random transitions between the expression states of a multistable genetic circuit. That way, heterogeneity within isogenic populations is created. Our aim is to understand which dynamical features of a simple autoregulatory system determine its intrinsic noise level, and how they can be modified in order to regulate state-transitions. To that end, novel mathematical methods for the state-dependent characterization and prediction of noise in multistable systems are developed. First, we introduce the hybrid LNA, a modified version of the Linear Noise Approximation. It yields good predictions on variances of mRNA and protein fluctuations, even for reaction systems comprising low-copy-number components (e.g. mRNA) and highly nonlinear reaction rates. Furthermore, the temporal structure of fluctuations and the skewness of the protein distribution are characterized via state-dependent protein burst sizes and burst frequencies. Based on this mathematical framework, we develop graphical methods which support the intuitive design of regulatory circuits with a desired noise pattern. The methods are then used to predict how overall noise in the system can be adapted, and how state-specific noise modifications are possible that allow, e.g., the generation of unidirectional transitions. Our considerations are validated by stochastic simulations. This way, a design of genetic circuits is possible that takes population heterogeneity into account and is valuable in applications of synthetic biology and biotechnology. Moreover, natural phenomena like the bimodal development of genetic competence can be studied.
Collapse
Affiliation(s)
- Sayuri Katharina Hortsch
- Systems Biotechnology, Faculty of Mechanical Engineering, Technical University of Munich, Garching, Germany
- * E-mail:
| | - Andreas Kremling
- Systems Biotechnology, Faculty of Mechanical Engineering, Technical University of Munich, Garching, Germany
| |
Collapse
|
36
|
Vigouroux A, Oldewurtel E, Cui L, Bikard D, van Teeffelen S. Tuning dCas9's ability to block transcription enables robust, noiseless knockdown of bacterial genes. Mol Syst Biol 2018; 14:e7899. [PMID: 29519933 PMCID: PMC5842579 DOI: 10.15252/msb.20177899] [Citation(s) in RCA: 74] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2017] [Revised: 02/08/2018] [Accepted: 02/14/2018] [Indexed: 02/06/2023] Open
Abstract
Over the past few years, tools that make use of the Cas9 nuclease have led to many breakthroughs, including in the control of gene expression. The catalytically dead variant of Cas9 known as dCas9 can be guided by small RNAs to block transcription of target genes, in a strategy also known as CRISPRi. Here, we reveal that the level of complementarity between the guide RNA and the target controls the rate at which RNA polymerase "kicks out" dCas9 from the target and completes transcription. We use this mechanism to precisely and robustly reduce gene expression by defined relative amounts. Alternatively, tuning repression by changing dCas9 concentration is noisy and promoter-strength dependent. We demonstrate broad applicability of this method to the study of genetic regulation and cellular physiology. First, we characterize feedback strength of a model auto-repressor. Second, we study the impact of amount variations of cell-wall synthesizing enzymes on cell morphology. Finally, we multiplex the system to obtain any combination of fractional repression of two genes.
Collapse
Affiliation(s)
- Antoine Vigouroux
- Synthetic Biology Laboratory, Institut Pasteur, Paris, France
- Microbial Morphogenesis and Growth Laboratory, Institut Pasteur, Paris, France
| | - Enno Oldewurtel
- Microbial Morphogenesis and Growth Laboratory, Institut Pasteur, Paris, France
| | - Lun Cui
- Synthetic Biology Laboratory, Institut Pasteur, Paris, France
| | - David Bikard
- Synthetic Biology Laboratory, Institut Pasteur, Paris, France
| | - Sven van Teeffelen
- Microbial Morphogenesis and Growth Laboratory, Institut Pasteur, Paris, France
| |
Collapse
|
37
|
Guillemin A, Richard A, Gonin-Giraud S, Gandrillon O. Automated cell cycle and cell size measurements for single-cell gene expression studies. BMC Res Notes 2018; 11:92. [PMID: 29391045 PMCID: PMC5796519 DOI: 10.1186/s13104-018-3195-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2017] [Accepted: 01/23/2018] [Indexed: 11/10/2022] Open
Abstract
OBJECTIVES Recent rise of single-cell studies revealed the importance of understanding the role of cell-to-cell variability, especially at the transcriptomic level. One of the numerous sources of cell-to-cell variation in gene expression is the heterogeneity in cell proliferation state. In order to identify how cell cycle and cell size influences gene expression variability at the single-cell level, we provide an universal and automatic toxic-free label method, compatible with single-cell high-throughput RT-qPCR. The method consists of isolating cells after a double-stained, analyzing their morphological parameters and performing a transcriptomic analysis on the same identified cells. RESULTS This led to an unbiased gene expression analysis and could be also used for improving single-cell tracking and imaging when combined with cell isolation. As an application for this technique, we showed that cell-to-cell variability in chicken erythroid progenitors was negligibly influenced by cell size nor cell cycle.
Collapse
Affiliation(s)
- Anissa Guillemin
- Laboratoire de biologie et modélisation de la cellule. LBMC-Ecole Normale Supérieure-Lyon, Université Claude Bernard Lyon 1, Institut National de la Santé et de la Recherche Médicale: U1210-Ecole Normale Supérieure de Lyon, Centre National de la Recherche Scientifique: UMR5239, 46 Allée d’Italie, 69007 Lyon, France
| | - Angélique Richard
- Laboratoire de biologie et modélisation de la cellule. LBMC-Ecole Normale Supérieure-Lyon, Université Claude Bernard Lyon 1, Institut National de la Santé et de la Recherche Médicale: U1210-Ecole Normale Supérieure de Lyon, Centre National de la Recherche Scientifique: UMR5239, 46 Allée d’Italie, 69007 Lyon, France
| | - Sandrine Gonin-Giraud
- Laboratoire de biologie et modélisation de la cellule. LBMC-Ecole Normale Supérieure-Lyon, Université Claude Bernard Lyon 1, Institut National de la Santé et de la Recherche Médicale: U1210-Ecole Normale Supérieure de Lyon, Centre National de la Recherche Scientifique: UMR5239, 46 Allée d’Italie, 69007 Lyon, France
| | - Olivier Gandrillon
- Laboratoire de biologie et modélisation de la cellule. LBMC-Ecole Normale Supérieure-Lyon, Université Claude Bernard Lyon 1, Institut National de la Santé et de la Recherche Médicale: U1210-Ecole Normale Supérieure de Lyon, Centre National de la Recherche Scientifique: UMR5239, 46 Allée d’Italie, 69007 Lyon, France
- Inria Dracula, 69100 Villeurbanne, France
| |
Collapse
|
38
|
Crocker J, Ilsley GR. Using synthetic biology to study gene regulatory evolution. Curr Opin Genet Dev 2017; 47:91-101. [DOI: 10.1016/j.gde.2017.09.001] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2017] [Revised: 09/06/2017] [Accepted: 09/11/2017] [Indexed: 12/21/2022]
|
39
|
Dey A, Barik D. Parallel arrangements of positive feedback loops limit cell-to-cell variability in differentiation. PLoS One 2017; 12:e0188623. [PMID: 29186174 PMCID: PMC5706692 DOI: 10.1371/journal.pone.0188623] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2017] [Accepted: 11/10/2017] [Indexed: 02/04/2023] Open
Abstract
Cellular differentiations are often regulated by bistable switches resulting from specific arrangements of multiple positive feedback loops (PFL) fused to one another. Although bistability generates digital responses at the cellular level, stochasticity in chemical reactions causes population heterogeneity in terms of its differentiated states. We hypothesized that the specific arrangements of PFLs may have evolved to minimize the cellular heterogeneity in differentiation. In order to test this we investigated variability in cellular differentiation controlled either by parallel or serial arrangements of multiple PFLs having similar average properties under extrinsic and intrinsic noises. We find that motifs with PFLs fused in parallel to one another around a central regulator are less susceptible to noise as compared to the motifs with PFLs arranged serially. Our calculations suggest that the increased resistance to noise in parallel motifs originate from the less sensitivity of bifurcation points to the extrinsic noise. Whereas estimation of mean residence times indicate that stable branches of bifurcations are robust to intrinsic noise in parallel motifs as compared to serial motifs. Model conclusions are consistent both in AND- and OR-gate input signal configurations and also with two different modeling strategies. Our investigations provide some insight into recent findings that differentiation of preadipocyte to mature adipocyte is controlled by network of parallel PFLs.
Collapse
Affiliation(s)
- Anupam Dey
- School of Chemistry, University of Hyderabad, Central University P.O., Hyderabad, Telangana, India
| | - Debashis Barik
- School of Chemistry, University of Hyderabad, Central University P.O., Hyderabad, Telangana, India
- * E-mail:
| |
Collapse
|
40
|
Mann M, Mehta A, Zhao JL, Lee K, Marinov GK, Garcia-Flores Y, Lu LF, Rudensky AY, Baltimore D. An NF-κB-microRNA regulatory network tunes macrophage inflammatory responses. Nat Commun 2017; 8:851. [PMID: 29021573 PMCID: PMC5636846 DOI: 10.1038/s41467-017-00972-z] [Citation(s) in RCA: 174] [Impact Index Per Article: 24.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2017] [Accepted: 08/10/2017] [Indexed: 01/02/2023] Open
Abstract
The innate inflammatory response must be tightly regulated to ensure effective immune protection. NF-κB is a key mediator of the inflammatory response, and its dysregulation has been associated with immune-related malignancies. Here, we describe a miRNA-based regulatory network that enables precise NF-κB activity in mouse macrophages. Elevated miR-155 expression potentiates NF-κB activity in miR-146a-deficient mice, leading to both an overactive acute inflammatory response and chronic inflammation. Enforced miR-155 expression overrides miR-146a-mediated repression of NF-κB activation, thus emphasizing the dominant function of miR-155 in promoting inflammation. Moreover, miR-155-deficient macrophages exhibit a suboptimal inflammatory response when exposed to low levels of inflammatory stimuli. Importantly, we demonstrate a temporal asymmetry between miR-155 and miR-146a expression during macrophage activation, which creates a combined positive and negative feedback network controlling NF-κB activity. This miRNA-based regulatory network enables a robust yet time-limited inflammatory response essential for functional immunity. MicroRNAs (miR) are important regulators of gene transcription, with miR-155 and miR-146a both implicated in macrophage activation. Here the authors show that NF-κB signalling, miR-155 and miR-146a form a complex network of cross-regulations to control gene transcription in macrophages for modulating inflammatory responses.
Collapse
Affiliation(s)
- Mati Mann
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, 91125, USA.
| | - Arnav Mehta
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, 91125, USA.,David Geffen School of Medicine, University of California, Los Angeles, CA, 90095, USA
| | - Jimmy L Zhao
- Department of Medicine, New York Presbyterian Hospital, Weill Cornell Medical College, 525 E 68th Street, New York, NY, 10065, USA.,Division of Hematology Oncology, Memorial Sloan Kettering Cancer Center, 1275 York Avenue, New York, NY, 10065, USA
| | - Kevin Lee
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, 91125, USA
| | - Georgi K Marinov
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, 91125, USA
| | - Yvette Garcia-Flores
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, 91125, USA
| | - Li-Fan Lu
- Division of Biological Sciences, University of California, La Jolla, San Diego, CA, 92093, USA.,Moores Cancer Center, University of California, La Jolla, San Diego, CA, 92093, USA.,Center for Microbiome Innovation, University of California, La Jolla, San Diego, CA, 92093, USA
| | - Alexander Y Rudensky
- Howard Hughes Medical Institute and Immunology Program, Ludwig Center at Memorial Sloan-Kettering Cancer Center, Memorial Sloan-Kettering Cancer Center, New York, NY, 10065, USA
| | - David Baltimore
- Division of Biology and Biological Engineering, California Institute of Technology, Pasadena, CA, 91125, USA.
| |
Collapse
|
41
|
Rodrigo G, Bajić D, Elola I, Poyatos JF. Deconstructing a multiple antibiotic resistance regulation through the quantification of its input function. NPJ Syst Biol Appl 2017; 3:30. [PMID: 29018569 PMCID: PMC5630622 DOI: 10.1038/s41540-017-0031-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2017] [Revised: 09/01/2017] [Accepted: 09/12/2017] [Indexed: 12/21/2022] Open
Abstract
Many essential bacterial responses present complex transcriptional regulation of gene expression. To what extent can the study of these responses substantiate the logic of their regulation? Here, we show how the input function of the genes constituting the response, i.e., the information of how their transcription rates change as function of the signals acting on the regulators, can serve as a quantitative tool to deconstruct the corresponding regulatory logic. To demonstrate this approach, we consider the multiple antibiotic resistance (mar) response in Escherichia coli. By characterizing the input function of its representative genes in wild-type and mutant bacteria, we recognize a dual autoregulation motif as main determinant of the response, which is further adjusted by the interplay with other regulators. We show that basic attributes, like its reaction to a wide range of stress or its moderate expression change, are associated with a strong negative autoregulation, while others, like the buffering of metabolic signals or the lack of memory to previous stress, are related to a weak positive autoregulation. With a mathematical model of the input functions, we identify some constraints fixing the molecular attributes of the regulators, and also notice the relevance of the bicystronic architecture harboring the dual autoregulation that is unique in E. coli. The input function emerges then as a tool to disentangle the rationale behind most of the attributes defining the mar phenotype. Overall, the present study supports the value of characterizing input functions to deconstruct the complexity of regulatory architectures in prokaryotic and eukaryotic systems. Many cellular responses result from the integration of numerous regulatory signals. To deconstruct the regulation of one of these responses, which enables resistance to multiple antibiotics in Escherichia coli, a team led by Juan F. Poyatos at the National Center for Biotechnology in Madrid studied the response input function by combining theoretical models and experiments. This function quantifies the rate of transcription of the genes constituting the response with respect to the signals acting on its cognate regulators. By examining how the shape of the function changes in different situations, e.g., when a given regulator is mutated, the team identified the implications for the specificity and dynamics of the response of a dual autoregulation at the core of the control architecture. The use of input functions as quantitative tools allows us to reverse engineer the complex regulations that dictate essential physiological functions in both prokaryotic and eukaryotic cells.
Collapse
Affiliation(s)
- Guillermo Rodrigo
- Instituto de Biología Molecular y Celular de Plantas, CSIC-UPV, 46022 Valencia, Spain
| | - Djordje Bajić
- Logic of Genomic Systems Laboratory, CNB-CSIC, 28049 Madrid, Spain.,Present Address: Department of Ecology and Evolutionary Biology, Yale University, New Haven, USA
| | - Ignacio Elola
- Logic of Genomic Systems Laboratory, CNB-CSIC, 28049 Madrid, Spain
| | - Juan F Poyatos
- Logic of Genomic Systems Laboratory, CNB-CSIC, 28049 Madrid, Spain
| |
Collapse
|
42
|
Brophy JAN, LaRue T, Dinneny JR. Understanding and engineering plant form. Semin Cell Dev Biol 2017; 79:68-77. [PMID: 28864344 DOI: 10.1016/j.semcdb.2017.08.051] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2017] [Revised: 08/25/2017] [Accepted: 08/28/2017] [Indexed: 11/18/2022]
Abstract
A plant's form is an important determinant of its fitness and economic value. Here, we review strategies for producing plants with altered forms. Historically, the process of changing a plant's form has been slow in agriculture, requiring iterative rounds of growth and selection. We discuss modern techniques for identifying genes involved in the development of plant form and tools that will be needed to effectively design and engineer plants with altered forms. Synthetic genetic circuits are highlighted for their potential to generate novel plant forms. We emphasize understanding development as a prerequisite to engineering and discuss the potential role of computer models in translating knowledge about single genes or pathways into a more comprehensive understanding of development.
Collapse
Affiliation(s)
- Jennifer A N Brophy
- Carnegie Institution for Science, Department of Plant Biology, Stanford, CA 94305, USA
| | - Therese LaRue
- Stanford University, Department of Biology, Stanford, CA 94305, USA
| | - José R Dinneny
- Carnegie Institution for Science, Department of Plant Biology, Stanford, CA 94305, USA.
| |
Collapse
|
43
|
A kinetic model of multiple phenotypic states for breast cancer cells. Sci Rep 2017; 7:9890. [PMID: 28852133 PMCID: PMC5574983 DOI: 10.1038/s41598-017-10321-1] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2017] [Accepted: 08/07/2017] [Indexed: 12/31/2022] Open
Abstract
Quantitative modeling of microscopic genes regulatory mechanisms in an individual cell is a crucial step towards understanding various macroscopic physiological phenomena of cell populations. Based on the regulatory mechanisms of genes zeb1 and cdh1 in the growth and development of breast cancer cells, we propose a kinetic model at the level of single cell. By constructing the effective landscape of underlying stationary probability for the genes expressions, it is found that (i) each breast cancer cell has three phenotypic states (i.e., the stem-like, basal, and luminal states) which correspond to three attractions of the probability landscape. (ii) The interconversions between phenotypic states can be induced by the noise intensity and the property of phenotypic switching is quantified by the mean first-passage time. (iii) Under certain conditions, the probabilities of each cancer cell appearing in the three states are consistent with the macroscopic phenotypic equilibrium proportions in the breast cancer SUM159 cell line. (iv) Our kinetic model involving the TGF-β signal can also qualitatively explain several macroscopic physiological phenomena of breast cancer cells, such as the "TGF-β paradox" in tumor therapy, the five clinical subtypes of breast cancer cells, and the effects of transient TGF-β on breast cancer metastasis.
Collapse
|
44
|
Kim JK, Sontag ED. Reduction of multiscale stochastic biochemical reaction networks using exact moment derivation. PLoS Comput Biol 2017; 13:e1005571. [PMID: 28582397 PMCID: PMC5481150 DOI: 10.1371/journal.pcbi.1005571] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2017] [Revised: 06/22/2017] [Accepted: 05/12/2017] [Indexed: 12/14/2022] Open
Abstract
Biochemical reaction networks (BRNs) in a cell frequently consist of reactions with disparate timescales. The stochastic simulations of such multiscale BRNs are prohibitively slow due to high computational cost for the simulations of fast reactions. One way to resolve this problem uses the fact that fast species regulated by fast reactions quickly equilibrate to their stationary distribution while slow species are unlikely to be changed. Thus, on a slow timescale, fast species can be replaced by their quasi-steady state (QSS): their stationary conditional expectation values for given slow species. As the QSS are determined solely by the state of slow species, such replacement leads to a reduced model, where fast species are eliminated. However, it is challenging to derive the QSS in the presence of nonlinear reactions. While various approximation schemes for the QSS have been developed, they often lead to considerable errors. Here, we propose two classes of multiscale BRNs which can be reduced by deriving an exact QSS rather than approximations. Specifically, if fast species constitute either a feedforward network or a complex balanced network, the reduced model based on the exact QSS can be derived. Such BRNs are frequently observed in a cell as the feedforward network is one of fundamental motifs of gene or protein regulatory networks. Furthermore, complex balanced networks also include various types of fast reversible bindings such as bindings between transcriptional factors and gene regulatory sites. The reduced models based on exact QSS, which can be calculated by the computational packages provided in this work, accurately approximate the slow scale dynamics of the original full model with much lower computational cost.
Collapse
Affiliation(s)
- Jae Kyoung Kim
- Department of Mathematical Sciences, Korea Advanced Institute of Science and Technology, Daejeon, Korea
- * E-mail: (JKK); , (EDS)
| | - Eduardo D. Sontag
- Department of Mathematics and Center for Quantitative Biology, Rutgers University, New Brunswick, New Jersey, United States of America
- * E-mail: (JKK); , (EDS)
| |
Collapse
|
45
|
Haerter JO, Díaz-Guilera A, Serrano MÁ. Noise-induced polarization switching in complex networks. Phys Rev E 2017; 95:042305. [PMID: 28505785 DOI: 10.1103/physreve.95.042305] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2016] [Indexed: 11/07/2022]
Abstract
The combination of bistability and noise is ubiquitous in complex systems, from biology to social interactions, and has important implications for their functioning and resilience. Here we use a simple three-state dynamical process, in which nodes go from one pole to another through an intermediate state, to show that noise can induce polarization switching in bistable systems if dynamical correlations are significant. In large, fully connected networks, where dynamical correlations can be neglected, increasing noise yields a collapse of bistability to an unpolarized configuration where the three possible states of the nodes are equally likely. In contrast, increased noise induces abrupt and irreversible polarization switching in sparsely connected networks. In multiplexes, where each layer can have a different polarization tendency, one layer is dominant and progressively imposes its polarization state on the other, offsetting or promoting the ability of noise to switch its polarization. Overall, we show that the interplay of noise and dynamical correlations can yield discontinuous transitions between extremes, which cannot be explained by a simple mean-field description.
Collapse
Affiliation(s)
- Jan O Haerter
- Departament de Física de la Matèria Condensada, Universitat de Barcelona, Martí i Franquès 1, 08028 Barcelona, Spain.,Niels Bohr Institutet, Københavns Universitet, Blegdamsvej 17, 2100 København Ø, Denmark
| | - Albert Díaz-Guilera
- Departament de Física de la Matèria Condensada, Universitat de Barcelona, Martí i Franquès 1, 08028 Barcelona, Spain.,Universitat de Barcelona Institute of Complex Systems (UBICS), Universitat de Barcelona, Barcelona, Spain
| | - M Ángeles Serrano
- Departament de Física de la Matèria Condensada, Universitat de Barcelona, Martí i Franquès 1, 08028 Barcelona, Spain.,Universitat de Barcelona Institute of Complex Systems (UBICS), Universitat de Barcelona, Barcelona, Spain.,ICREA, Pg. Lluís Companys 23, 08010 Barcelona, Spain
| |
Collapse
|
46
|
Forment M, Rodrigo G. Molecular noise can minimize the collective sensitivity of a clonal heterogeneous cell population. J Theor Biol 2016; 416:38-44. [PMID: 28043818 DOI: 10.1016/j.jtbi.2016.12.023] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2016] [Revised: 12/19/2016] [Accepted: 12/28/2016] [Indexed: 11/17/2022]
Abstract
It is now widely accepted that molecular noise, rather than be always detrimental, introduces in many circumstances the required boost to reach fundamental cellular activities or strategies otherwise unattainable. In threshold-like genetic systems, molecular noise serves to generate heterogeneous responses in a clonal population, also in a tissue, due to cell-to-cell variability. Here, we derived a mathematical framework from which we could study in detail this effect. We focused on a minimal decision-making gene circuit implemented as a transcriptional positive-feedback loop. We evidenced that when the individual responses of each cell within the population are averaged, a sort of collective behavior, the resulting dose-response curve is linearized. In other words, the population is less sensitive than the individuals, which otherwise enhances the information transfer from signal to response. We found that the distance to the ideal linear response of the cell population is minimized for a particular noise level, and also characterized the interplay between intrinsic and extrinsic noise. Overall, our results highlight how cells could, by acting as a collective, entangle their genetic systems with their environments by adjusting the intracellular noise levels.
Collapse
Affiliation(s)
- Marzo Forment
- Instituto de Biología Molecular y Celular de Plantas, CSIC - UPV, 46022 Valencia, Spain
| | - Guillermo Rodrigo
- Instituto de Biología Molecular y Celular de Plantas, CSIC - UPV, 46022 Valencia, Spain.
| |
Collapse
|
47
|
Richard A, Boullu L, Herbach U, Bonnafoux A, Morin V, Vallin E, Guillemin A, Papili Gao N, Gunawan R, Cosette J, Arnaud O, Kupiec JJ, Espinasse T, Gonin-Giraud S, Gandrillon O. Single-Cell-Based Analysis Highlights a Surge in Cell-to-Cell Molecular Variability Preceding Irreversible Commitment in a Differentiation Process. PLoS Biol 2016; 14:e1002585. [PMID: 28027290 PMCID: PMC5191835 DOI: 10.1371/journal.pbio.1002585] [Citation(s) in RCA: 128] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2016] [Accepted: 09/22/2016] [Indexed: 12/31/2022] Open
Abstract
In some recent studies, a view emerged that stochastic dynamics governing the switching of cells from one differentiation state to another could be characterized by a peak in gene expression variability at the point of fate commitment. We have tested this hypothesis at the single-cell level by analyzing primary chicken erythroid progenitors through their differentiation process and measuring the expression of selected genes at six sequential time-points after induction of differentiation. In contrast to population-based expression data, single-cell gene expression data revealed a high cell-to-cell variability, which was masked by averaging. We were able to show that the correlation network was a very dynamical entity and that a subgroup of genes tend to follow the predictions from the dynamical network biomarker (DNB) theory. In addition, we also identified a small group of functionally related genes encoding proteins involved in sterol synthesis that could act as the initial drivers of the differentiation. In order to assess quantitatively the cell-to-cell variability in gene expression and its evolution in time, we used Shannon entropy as a measure of the heterogeneity. Entropy values showed a significant increase in the first 8 h of the differentiation process, reaching a peak between 8 and 24 h, before decreasing to significantly lower values. Moreover, we observed that the previous point of maximum entropy precedes two paramount key points: an irreversible commitment to differentiation between 24 and 48 h followed by a significant increase in cell size variability at 48 h. In conclusion, when analyzed at the single cell level, the differentiation process looks very different from its classical population average view. New observables (like entropy) can be computed, the behavior of which is fully compatible with the idea that differentiation is not a “simple” program that all cells execute identically but results from the dynamical behavior of the underlying molecular network. A single-cell transcriptomics analysis offers a new dynamical view of the differentiation process, involving an increase in between-cell variability prior to commitment. The differentiation process has classically been seen as a stereotyped program leading from one progenitor toward a functional cell. This vision was based upon cell population-based analyses averaged over millions of cells. However, new methods have recently emerged that allow interrogation of the molecular content at the single-cell level, challenging this view with a new model suggesting that cell-to-cell gene expression stochasticity could play a key role in differentiation. We took advantage of a physiologically relevant avian cellular model to analyze the expression level of 92 genes in individual cells collected at several time-points during differentiation. We first observed that the process analyzed at the single-cell level is very different and much less well ordered than the population-based average view. Furthermore, we showed that cell-to-cell variability in gene expression peaks transiently before strongly decreasing. This rise in variability precedes two key events: an irreversible commitment to differentiation, followed by a significant increase in cell size variability. Altogether, our results support the idea that differentiation is not a “simple” series of well-ordered molecular events executed identically by all cells in a population but likely results from dynamical behavior of the underlying molecular network.
Collapse
Affiliation(s)
- Angélique Richard
- Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, Laboratory of Biology and Modelling of the Cell, 46 allée d’Italie Site Jacques Monod, F-69007, Lyon, France
| | - Loïs Boullu
- Inria Team Dracula, Inria Center Grenoble Rhône-Alpes, France
- Université de Lyon, Université Lyon 1, CNRS UMR 5208, Institut Camille Jordan 43 blvd du 11 novembre 1918, F-69622 Villeurbanne-Cedex, France
- Département de Mathématiques et de statistiques de l’Université de Montréal, Pavillon André-Aisenstadt, 2920, chemin de la Tour, Montréal (Québec) H3T 1J4 Canada
| | - Ulysse Herbach
- Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, Laboratory of Biology and Modelling of the Cell, 46 allée d’Italie Site Jacques Monod, F-69007, Lyon, France
- Inria Team Dracula, Inria Center Grenoble Rhône-Alpes, France
- Université de Lyon, Université Lyon 1, CNRS UMR 5208, Institut Camille Jordan 43 blvd du 11 novembre 1918, F-69622 Villeurbanne-Cedex, France
| | - Arnaud Bonnafoux
- Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, Laboratory of Biology and Modelling of the Cell, 46 allée d’Italie Site Jacques Monod, F-69007, Lyon, France
- Inria Team Dracula, Inria Center Grenoble Rhône-Alpes, France
- The CoSMo company. 5 passage du Vercors – 69007 LYON – France
| | - Valérie Morin
- Univ Lyon, Univ Claude Bernard, CNRS UMR 5310 - INSERM U1217, Institut NeuroMyoGène, F-69622 Villeurbanne-Cedex, France
| | - Elodie Vallin
- Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, Laboratory of Biology and Modelling of the Cell, 46 allée d’Italie Site Jacques Monod, F-69007, Lyon, France
| | - Anissa Guillemin
- Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, Laboratory of Biology and Modelling of the Cell, 46 allée d’Italie Site Jacques Monod, F-69007, Lyon, France
| | - Nan Papili Gao
- Institute for Chemical and Bioengineering, ETH Zurich, Zurich, Switzerland
- Swiss Institute of Bioinformatics, Quartier Sorge - Batiment Genopode, 1015 Lausanne Switzerland
| | - Rudiyanto Gunawan
- Institute for Chemical and Bioengineering, ETH Zurich, Zurich, Switzerland
- Swiss Institute of Bioinformatics, Quartier Sorge - Batiment Genopode, 1015 Lausanne Switzerland
| | - Jérémie Cosette
- Genethon – Institut National de la Santé et de la Recherche Médicale – INSERM, Université d’Evry-Val-d’Essone – 1 rue de l’internationale 91000 Evry, France
| | - Ophélie Arnaud
- RIKEN - Center for Life Science Technologies (Division of Genomic Technologies)—CLST (DGT), 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| | | | - Thibault Espinasse
- Université de Lyon, Université Lyon 1, CNRS UMR 5208, Institut Camille Jordan 43 blvd du 11 novembre 1918, F-69622 Villeurbanne-Cedex, France
| | - Sandrine Gonin-Giraud
- Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, Laboratory of Biology and Modelling of the Cell, 46 allée d’Italie Site Jacques Monod, F-69007, Lyon, France
| | - Olivier Gandrillon
- Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, Laboratory of Biology and Modelling of the Cell, 46 allée d’Italie Site Jacques Monod, F-69007, Lyon, France
- Inria Team Dracula, Inria Center Grenoble Rhône-Alpes, France
- * E-mail:
| |
Collapse
|
48
|
Veliz-Cuba A, Gupta C, Bennett MR, Josić K, Ott W. Effects of cell cycle noise on excitable gene circuits. Phys Biol 2016; 13:066007. [PMID: 27902489 DOI: 10.1088/1478-3975/13/6/066007] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
We assess the impact of cell cycle noise on gene circuit dynamics. For bistable genetic switches and excitable circuits, we find that transitions between metastable states most likely occur just after cell division and that this concentration effect intensifies in the presence of transcriptional delay. We explain this concentration effect with a three-states stochastic model. For genetic oscillators, we quantify the temporal correlations between daughter cells induced by cell division. Temporal correlations must be captured properly in order to accurately quantify noise sources within gene networks.
Collapse
|
49
|
Bracken CP, Scott HS, Goodall GJ. A network-biology perspective of microRNA function and dysfunction in cancer. Nat Rev Genet 2016; 17:719-732. [DOI: 10.1038/nrg.2016.134] [Citation(s) in RCA: 468] [Impact Index Per Article: 58.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
50
|
Rodrigo G, Bajic D, Elola I, Poyatos JF. Antagonistic autoregulation speeds up a homogeneous response in Escherichia coli. Sci Rep 2016; 6:36196. [PMID: 27796341 PMCID: PMC5086920 DOI: 10.1038/srep36196] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2016] [Accepted: 10/10/2016] [Indexed: 01/15/2023] Open
Abstract
By integrating positive and negative feedback loops, biological systems establish intricate gene expression patterns linked to multistability, pulsing, and oscillations. This depends on the specific characteristics of each interlinked feedback, and thus one would expect additional expression programs to be found. Here, we investigate one such program associated with an antagonistic positive and negative transcriptional autoregulatory motif derived from the multiple antibiotic resistance (mar) system of Escherichia coli. We studied the dynamics of the system by combining a predictive mathematical model with high-resolution experimental measures of the response both at the population and single-cell level. We show that in this motif the weak positive autoregulation does not slow down but rather enhances response speedup in combination with a strong negative feedback loop. This balance of feedback strengths anticipates a homogeneous population phenotype, which we corroborate experimentally. Theoretical analysis also emphasized the specific molecular properties that determine the dynamics of the mar phenotype. More broadly, response acceleration could provide a rationale for the presence of weak positive feedbacks in other biological scenarios exhibiting these interlinked regulatory architectures.
Collapse
Affiliation(s)
- Guillermo Rodrigo
- Instituto de Biología Molecular y Cellular de Plantas, CSIC–UPV, 46022 Valencia, Spain
| | - Djordje Bajic
- Logic of Genomic Systems Laboratory, CNB–CSIC, 28049 Madrid, Spain
| | - Ignacio Elola
- Logic of Genomic Systems Laboratory, CNB–CSIC, 28049 Madrid, Spain
| | - Juan F. Poyatos
- Logic of Genomic Systems Laboratory, CNB–CSIC, 28049 Madrid, Spain
| |
Collapse
|