1
|
Jia C, Grima R. Holimap: an accurate and efficient method for solving stochastic gene network dynamics. Nat Commun 2024; 15:6557. [PMID: 39095346 PMCID: PMC11297302 DOI: 10.1038/s41467-024-50716-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2024] [Accepted: 07/13/2024] [Indexed: 08/04/2024] Open
Abstract
Gene-gene interactions are crucial to the control of sub-cellular processes but our understanding of their stochastic dynamics is hindered by the lack of simulation methods that can accurately and efficiently predict how the distributions of gene product numbers vary across parameter space. To overcome these difficulties, here we present Holimap (high-order linear-mapping approximation), an approach that approximates the protein or mRNA number distributions of a complex gene regulatory network by the distributions of a much simpler reaction system. We demonstrate Holimap's computational advantages over conventional methods by applying it to predict the stochastic time-dependent dynamics of various gene networks, including transcriptional networks ranging from simple autoregulatory loops to complex randomly connected networks, post-transcriptional networks, and post-translational networks. Holimap is ideally suited to study how the intricate network of gene-gene interactions results in precise coordination and control of gene expression.
Collapse
Affiliation(s)
- Chen Jia
- Applied and Computational Mathematics Division, Beijing Computational Science Research Center, Beijing, China
| | - Ramon Grima
- School of Biological Sciences, University of Edinburgh, Edinburgh, UK.
| |
Collapse
|
2
|
Nemsick S, Hansen AS. Molecular models of bidirectional promoter regulation. Curr Opin Struct Biol 2024; 87:102865. [PMID: 38905929 DOI: 10.1016/j.sbi.2024.102865] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2023] [Revised: 03/30/2024] [Accepted: 05/27/2024] [Indexed: 06/23/2024]
Abstract
Approximately 11% of human genes are transcribed by a bidirectional promoter (BDP), defined as two genes with <1 kb between their transcription start sites. Despite their evolutionary conservation and enrichment for housekeeping genes and oncogenes, the regulatory role of BDPs remains unclear. BDPs have been suggested to facilitate gene coregulation and/or decrease expression noise. This review discusses these potential regulatory functions through the context of six prospective underlying mechanistic models: a single nucleosome free region, shared transcription factor/regulator binding, cooperative negative supercoiling, bimodal histone marks, joint activation by enhancer(s), and RNA-mediated recruitment of regulators. These molecular mechanisms may act independently and/or cooperatively to facilitate the coregulation and/or decreased expression noise predicted of BDPs.
Collapse
Affiliation(s)
- Sarah Nemsick
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA; The Gene Regulation Observatory, Broad Institute of MIT and Harvard, Cambridge, MA 02139, USA; Koch Institute for Integrative Cancer Research, Cambridge, MA 02139, USA
| | - Anders S Hansen
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA; The Gene Regulation Observatory, Broad Institute of MIT and Harvard, Cambridge, MA 02139, USA; Koch Institute for Integrative Cancer Research, Cambridge, MA 02139, USA.
| |
Collapse
|
3
|
Lin WH, Opoc FG, Liao CW, Roy K, Steinmetz L, Leu JY. Histone deacetylase Hos2 regulates protein expression noise by potentially modulating the protein translation machinery. Nucleic Acids Res 2024; 52:7556-7571. [PMID: 38783136 PMCID: PMC11260488 DOI: 10.1093/nar/gkae432] [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: 02/02/2024] [Revised: 05/05/2024] [Accepted: 05/09/2024] [Indexed: 05/25/2024] Open
Abstract
Non-genetic variations derived from expression noise at transcript or protein levels can result in cell-to-cell heterogeneity within an isogenic population. Although cells have developed strategies to reduce noise in some cellular functions, this heterogeneity can also facilitate varying levels of regulation and provide evolutionary benefits in specific environments. Despite several general characteristics of cellular noise having been revealed, the detailed molecular pathways underlying noise regulation remain elusive. Here, we established a dual-fluorescent reporter system in Saccharomyces cerevisiae and performed experimental evolution to search for mutations that increase expression noise. By analyzing evolved cells using bulk segregant analysis coupled with whole-genome sequencing, we identified the histone deacetylase Hos2 as a negative noise regulator. A hos2 mutant down-regulated multiple ribosomal protein genes and exhibited partially compromised protein translation, indicating that Hos2 may regulate protein expression noise by modulating the translation machinery. Treating cells with translation inhibitors or introducing mutations into several Hos2-regulated ribosomal protein genes-RPS9A, RPS28B and RPL42A-enhanced protein expression noise. Our study provides an effective strategy for identifying noise regulators and also sheds light on how cells regulate non-genetic variation through protein translation.
Collapse
Affiliation(s)
- Wei-Han Lin
- Doctoral Program in Microbial Genomics, National Chung Hsing University and Academia Sinica, Taiwan
- Institute of Molecular Biology, Academia Sinica, Taipei 115, Taiwan
| | - Florica J G Opoc
- Institute of Molecular Biology, Academia Sinica, Taipei 115, Taiwan
| | - Chia-Wei Liao
- Institute of Molecular Biology, Academia Sinica, Taipei 115, Taiwan
| | - Kevin R Roy
- Stanford Genome Technology Center, Stanford University, Palo Alto, CA 94304, USA
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA
| | - Lars M Steinmetz
- Stanford Genome Technology Center, Stanford University, Palo Alto, CA 94304, USA
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA
- European Molecular Biology Laboratory (EMBL), Genome Biology Unit, Heidelberg 69117, Germany
| | - Jun-Yi Leu
- Doctoral Program in Microbial Genomics, National Chung Hsing University and Academia Sinica, Taiwan
- Institute of Molecular Biology, Academia Sinica, Taipei 115, Taiwan
| |
Collapse
|
4
|
Pal S, Dhar R. Living in a noisy world-origins of gene expression noise and its impact on cellular decision-making. FEBS Lett 2024; 598:1673-1691. [PMID: 38724715 DOI: 10.1002/1873-3468.14898] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Revised: 03/23/2024] [Accepted: 03/27/2024] [Indexed: 07/23/2024]
Abstract
The expression level of a gene can vary between genetically identical cells under the same environmental condition-a phenomenon referred to as gene expression noise. Several studies have now elucidated a central role of transcription factors in the generation of expression noise. Transcription factors, as the key components of gene regulatory networks, drive many important cellular decisions in response to cellular and environmental signals. Therefore, a very relevant question is how expression noise impacts gene regulation and influences cellular decision-making. In this Review, we summarize the current understanding of the molecular origins of expression noise, highlighting the role of transcription factors in this process, and discuss the ways in which noise can influence cellular decision-making. As advances in single-cell technologies open new avenues for studying expression noise as well as gene regulatory circuits, a better understanding of the influence of noise on cellular decisions will have important implications for many biological processes.
Collapse
Affiliation(s)
- Sampriti Pal
- Department of Bioscience and Biotechnology, IIT Kharagpur, India
| | - Riddhiman Dhar
- Department of Bioscience and Biotechnology, IIT Kharagpur, India
| |
Collapse
|
5
|
Wang L, Zhang K, Xu L, Wang J. Understanding underlying physical mechanism reveals early warning indicators and key elements for adaptive infections disease networks. PNAS NEXUS 2024; 3:pgae237. [PMID: 39035039 PMCID: PMC11259140 DOI: 10.1093/pnasnexus/pgae237] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/01/2024] [Accepted: 06/03/2024] [Indexed: 07/23/2024]
Abstract
The study of infectious diseases holds significant scientific and societal importance, yet current research on the mechanisms of disease emergence and prediction methods still face challenging issues. This research uses the landscape and flux theoretical framework to reveal the non-equilibrium dynamics of adaptive infectious diseases and uncover its underlying physical mechanism. This allows the quantification of dynamics, characterizing the system with two basins of attraction determined by gradient and rotational flux forces. Quantification of entropy production rates provides insights into the system deviating from equilibrium and associated dissipative costs. The study identifies early warning indicators for the critical transition, emphasizing the advantage of observing time irreversibility from time series over theoretical entropy production and flux. The presence of rotational flux leads to an irreversible pathway between disease states. Through global sensitivity analysis, we identified the key factors influencing infectious diseases. In summary, this research offers valuable insights into infectious disease dynamics and presents a practical approach for predicting the onset of critical transition, addressing existing research gaps.
Collapse
Affiliation(s)
- Linqi Wang
- Center of Theoretical Physics, College of Physics, Jilin University, Changchun, Jilin, 130012, China
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin, 130022, China
| | - Kun Zhang
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin, 130022, China
| | - Li Xu
- State Key Laboratory of Electroanalytical Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin, 130022, China
| | - Jin Wang
- Department of Chemistry, Physics and Astronomy, State University of New York at Stony Brook, Stony Brook, NY 11794, USA
| |
Collapse
|
6
|
Zhang Z, Zabaikina I, Nieto C, Vahdat Z, Bokes P, Singh A. Stochastic Gene Expression in Proliferating Cells: Differing Noise Intensity in Single-Cell and Population Perspectives. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.06.28.601263. [PMID: 38979195 PMCID: PMC11230457 DOI: 10.1101/2024.06.28.601263] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/10/2024]
Abstract
Random fluctuations (noise) in gene expression can be studied from two complementary perspectives: following expression in a single cell over time or comparing expression between cells in a proliferating population at a given time. Here, we systematically investigated scenarios where both perspectives lead to different levels of noise in a given gene product. We first consider a stable protein, whose concentration is diluted by cellular growth, and the protein inhibits growth at high concentrations, establishing a positive feedback loop. For a stochastic model with molecular bursting of gene products, we analytically predict and contrast the steady-state distributions of protein concentration in both frameworks. Although positive feedback amplifies the noise in expression, this amplification is much higher in the population framework compared to following a single cell over time. We also study other processes that lead to different noise levels even in the absence of such dilution-based feedback. When considering randomness in the partitioning of molecules between daughters during mitosis, we find that in the single-cell perspective, the noise in protein concentration is independent of noise in the cell cycle duration. In contrast, partitioning noise is amplified in the population perspective by increasing randomness in cell-cycle time. Overall, our results show that the commonly used single-cell framework that does not account for proliferating cells can, in some cases, underestimate the noise in gene product levels. These results have important implications for studying the inter-cellular variation of different stress-related expression programs across cell types that are known to inhibit cellular growth.
Collapse
Affiliation(s)
- Zhanhao Zhang
- Department of Electrical and Computer Engineering, University of Delaware. Newark, DE 19716, USA
| | - Iryna Zabaikina
- Department of Applied Mathematics and Statistics, Comenius University, Bratislava 84248, Slovakia
| | - César Nieto
- Department of Electrical and Computer Engineering, University of Delaware. Newark, DE 19716, USA
| | - Zahra Vahdat
- Department of Electrical and Computer Engineering, University of Delaware. Newark, DE 19716, USA
| | - Pavol Bokes
- Department of Applied Mathematics and Statistics, Comenius University, Bratislava 84248, Slovakia
| | - Abhyudai Singh
- Department of Electrical and Computer Engineering, University of Delaware. Newark, DE 19716, USA
| |
Collapse
|
7
|
Hong CKY, Ramu A, Zhao S, Cohen BA. Effect of genomic and cellular environments on gene expression noise. Genome Biol 2024; 25:137. [PMID: 38790076 PMCID: PMC11127367 DOI: 10.1186/s13059-024-03277-9] [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: 12/07/2022] [Accepted: 05/13/2024] [Indexed: 05/26/2024] Open
Abstract
BACKGROUND Individual cells from isogenic populations often display large cell-to-cell differences in gene expression. This "noise" in expression derives from several sources, including the genomic and cellular environment in which a gene resides. Large-scale maps of genomic environments have revealed the effects of epigenetic modifications and transcription factor occupancy on mean expression levels, but leveraging such maps to explain expression noise will require new methods to assay how expression noise changes at locations across the genome. RESULTS To address this gap, we present Single-cell Analysis of Reporter Gene Expression Noise and Transcriptome (SARGENT), a method that simultaneously measures the noisiness of reporter genes integrated throughout the genome and the global mRNA profiles of individual reporter-gene-containing cells. Using SARGENT, we perform the first comprehensive genome-wide survey of how genomic locations impact gene expression noise. We find that the mean and noise of expression correlate with different histone modifications. We quantify the intrinsic and extrinsic components of reporter gene noise and, using the associated mRNA profiles, assign the extrinsic component to differences between the CD24+ "stem-like" substate and the more "differentiated" substate. SARGENT also reveals the effects of transgene integrations on endogenous gene expression, which will help guide the search for "safe-harbor" loci. CONCLUSIONS Taken together, we show that SARGENT is a powerful tool to measure both the mean and noise of gene expression at locations across the genome and that the data generatd by SARGENT reveals important insights into the regulation of gene expression noise genome-wide.
Collapse
Affiliation(s)
- Clarice K Y Hong
- The Edison Family Center for Genome Sciences and Systems Biology, School of Medicine, Washington University in St. Louis, Saint Louis, MO, 63110, USA
- Department of Genetics, School of Medicine, Washington University in St. Louis, Saint Louis, MO, 63110, USA
| | - Avinash Ramu
- The Edison Family Center for Genome Sciences and Systems Biology, School of Medicine, Washington University in St. Louis, Saint Louis, MO, 63110, USA
- Department of Genetics, School of Medicine, Washington University in St. Louis, Saint Louis, MO, 63110, USA
| | - Siqi Zhao
- The Edison Family Center for Genome Sciences and Systems Biology, School of Medicine, Washington University in St. Louis, Saint Louis, MO, 63110, USA
- Department of Genetics, School of Medicine, Washington University in St. Louis, Saint Louis, MO, 63110, USA
| | - Barak A Cohen
- The Edison Family Center for Genome Sciences and Systems Biology, School of Medicine, Washington University in St. Louis, Saint Louis, MO, 63110, USA.
- Department of Genetics, School of Medicine, Washington University in St. Louis, Saint Louis, MO, 63110, USA.
| |
Collapse
|
8
|
Stevanovic M, Teuber Carvalho JP, Bittihn P, Schultz D. Dynamical model of antibiotic responses linking expression of resistance genes to metabolism explains emergence of heterogeneity during drug exposures. Phys Biol 2024; 21:036002. [PMID: 38412523 PMCID: PMC10988634 DOI: 10.1088/1478-3975/ad2d64] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2023] [Revised: 01/25/2024] [Accepted: 02/27/2024] [Indexed: 02/29/2024]
Abstract
Antibiotic responses in bacteria are highly dynamic and heterogeneous, with sudden exposure of bacterial colonies to high drug doses resulting in the coexistence of recovered and arrested cells. The dynamics of the response is determined by regulatory circuits controlling the expression of resistance genes, which are in turn modulated by the drug's action on cell growth and metabolism. Despite advances in understanding gene regulation at the molecular level, we still lack a framework to describe how feedback mechanisms resulting from the interdependence between expression of resistance and cell metabolism can amplify naturally occurring noise and create heterogeneity at the population level. To understand how this interplay affects cell survival upon exposure, we constructed a mathematical model of the dynamics of antibiotic responses that links metabolism and regulation of gene expression, based on the tetracycline resistancetetoperon inE. coli. We use this model to interpret measurements of growth and expression of resistance in microfluidic experiments, both in single cells and in biofilms. We also implemented a stochastic model of the drug response, to show that exposure to high drug levels results in large variations of recovery times and heterogeneity at the population level. We show that stochasticity is important to determine how nutrient quality affects cell survival during exposure to high drug concentrations. A quantitative description of how microbes respond to antibiotics in dynamical environments is crucial to understand population-level behaviors such as biofilms and pathogenesis.
Collapse
Affiliation(s)
- Mirjana Stevanovic
- Department of Microbiology and Immunology, Geisel School of Medicine at Dartmouth, Hanover, NH, United States of America
| | - João Pedro Teuber Carvalho
- Department of Microbiology and Immunology, Geisel School of Medicine at Dartmouth, Hanover, NH, United States of America
| | - Philip Bittihn
- Max Planck Institute for Dynamics and Self-Organization, Göttingen, Germany
- Institute for the Dynamics of Complex Systems, University of Göttingen, Göttingen, Germany
| | - Daniel Schultz
- Department of Microbiology and Immunology, Geisel School of Medicine at Dartmouth, Hanover, NH, United States of America
| |
Collapse
|
9
|
Giudice L, Mohamed A, Malm T. StellarPath: Hierarchical-vertical multi-omics classifier synergizes stable markers and interpretable similarity networks for patient profiling. PLoS Comput Biol 2024; 20:e1012022. [PMID: 38607982 PMCID: PMC11042724 DOI: 10.1371/journal.pcbi.1012022] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2023] [Revised: 04/24/2024] [Accepted: 03/25/2024] [Indexed: 04/14/2024] Open
Abstract
The Patient Similarity Network paradigm implies modeling the similarity between patients based on specific data. The similarity can summarize patients' relationships from high-dimensional data, such as biological omics. The end PSN can undergo un/supervised learning tasks while being strongly interpretable, tailored for precision medicine, and ready to be analyzed with graph-theory methods. However, these benefits are not guaranteed and depend on the granularity of the summarized data, the clarity of the similarity measure, the complexity of the network's topology, and the implemented methods for analysis. To date, no patient classifier fully leverages the paradigm's inherent benefits. PSNs remain complex, unexploited, and meaningless. We present StellarPath, a hierarchical-vertical patient classifier that leverages pathway analysis and patient similarity concepts to find meaningful features for both classes and individuals. StellarPath processes omics data, hierarchically integrates them into pathways, and uses a novel similarity to measure how patients' pathway activity is alike. It selects biologically relevant molecules, pathways, and networks, considering molecule stability and topology. A graph convolutional neural network then predicts unknown patients based on known cases. StellarPath excels in classification performances and computational resources across sixteen datasets. It demonstrates proficiency in inferring the class of new patients described in external independent studies, following its initial training and testing phases on a local dataset. It advances the PSN paradigm and provides new markers, insights, and tools for in-depth patient profiling.
Collapse
Affiliation(s)
- Luca Giudice
- A.I. Virtanen Institute for Molecular Sciences, University of Eastern Finland, Kuopio, Finland
| | - Ahmed Mohamed
- A.I. Virtanen Institute for Molecular Sciences, University of Eastern Finland, Kuopio, Finland
| | - Tarja Malm
- A.I. Virtanen Institute for Molecular Sciences, University of Eastern Finland, Kuopio, Finland
| |
Collapse
|
10
|
Murugesan A, Alshagrawi RA, Thiyagarajan R, Kandhavelu M. A dual fluorescence protein expression system detects cell cycle dependent protein noise. Int J Biol Macromol 2024; 263:130262. [PMID: 38378117 DOI: 10.1016/j.ijbiomac.2024.130262] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2023] [Revised: 02/14/2024] [Accepted: 02/15/2024] [Indexed: 02/22/2024]
Abstract
Inherently identical cells exhibit significant phenotypic variation. It can be essential for many biological processes and is known to arise from stochastic, 'noisy', gene expression that is determined by intrinsic and extrinsic components. It is now obvious that the noise varies as a function of inducer concentration. However, its fluctuation over the cell cycle is limited. Applying dual colour fluorescence protein reporter system, Cyan Fluorescent Protein (CFP) and Yellow fluorescent protein (YFP) tagged multi-copy plasmids, we determine variation of the noise components over the phases in lac promoter induced by Isopropyl β-D-1-thiogalactopyranoside (IPTG) and in presence of additional Magnesium, Mg2+ ion. We, also, estimate the how such system deviates from observations of single-copy plasmid. Found 25 % difference between multi-copy system and single-copy system clarifies that observed noise is considerable and estimates population behaviour during the cell cycle. We show that total variation in cells induced with IPTG is determined by higher extrinsic than intrinsic noise. It increases from Lag to Exponential phase and decreases from Retardation to Stationary phase. By observing slow and fast dividing cells, we show that 5 mM Mg2+ increases population homogeneity compared to 2.5 mM Mg2+ in the environment. The experimental data obtained using dual colour fluorescence protein reporter system demonstrates that protein expression noise, depending on intra cellular ionic concentration, is tightly controlled by phase of the cell.
Collapse
Affiliation(s)
- Akshaya Murugesan
- Molecular Signaling Group, Faculty of Medicine and Health Technology, Tampere University and BioMediTech, P.O. Box 553, 33101 Tampere, Finland.; Department of Biotechnology, Lady Doak College, Madurai Kamaraj University, Thallakulam, Madurai 625002, India
| | - Reshod A Alshagrawi
- Department of Food Science and Nutrition, College of Food Science and Agriculture, King Saud University, Riyadh, Saudi Arabia
| | - Ramesh Thiyagarajan
- Department of Basic Medical Sciences, College of Medicine, Prince Sattam Bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
| | - Meenakshisundaram Kandhavelu
- Molecular Signaling Group, Faculty of Medicine and Health Technology, Tampere University and BioMediTech, P.O. Box 553, 33101 Tampere, Finland..
| |
Collapse
|
11
|
Ilker E, Hinczewski M. Bioenergetic costs and the evolution of noise regulation by microRNAs. Proc Natl Acad Sci U S A 2024; 121:e2308796121. [PMID: 38386708 PMCID: PMC10907262 DOI: 10.1073/pnas.2308796121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2023] [Accepted: 01/14/2024] [Indexed: 02/24/2024] Open
Abstract
Noise control, together with other regulatory functions facilitated by microRNAs (miRNAs), is believed to have played important roles in the evolution of multicellular eukaryotic organisms. miRNAs can dampen protein fluctuations via enhanced degradation of messenger RNA (mRNA), but this requires compensation by increased mRNA transcription to maintain the same expression levels. The overall mechanism is metabolically expensive, leading to questions about how it might have evolved in the first place. We develop a stochastic model of miRNA noise regulation, coupled with a detailed analysis of the associated metabolic costs. Additionally, we calculate binding free energies for a range of miRNA seeds, the short sequences which govern target recognition. We argue that natural selection may have fine-tuned the Michaelis-Menten constant [Formula: see text] describing miRNA-mRNA affinity and show supporting evidence from analysis of experimental data. [Formula: see text] is constrained by seed length, and optimal noise control (minimum protein variance at a given energy cost) is achievable for seeds of 6 to 7 nucleotides in length, the most commonly observed types. Moreover, at optimality, the degree of noise reduction approaches the theoretical bound set by the Wiener-Kolmogorov linear filter. The results illustrate how selective pressure toward energy efficiency has potentially shaped a crucial regulatory pathway in eukaryotes.
Collapse
Affiliation(s)
- Efe Ilker
- Max Planck Institute for the Physics of Complex Systems, Dresden01187, Germany
| | - Michael Hinczewski
- Department of Physics, Case Western Reserve University, Cleveland, OH44106
| |
Collapse
|
12
|
Fromm B, Sorger T. Rapid adaptation of cellular metabolic rate to the MicroRNA complements of mammals and its relevance to the evolution of endothermy. iScience 2024; 27:108740. [PMID: 38327773 PMCID: PMC10847693 DOI: 10.1016/j.isci.2023.108740] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2022] [Revised: 09/13/2023] [Accepted: 12/12/2023] [Indexed: 02/09/2024] Open
Abstract
The metabolic efficiency of mammalian cells depends on the attenuation of intrinsic translation noise by microRNAs. We devised a metric of cellular metabolic rate (cMR), rMR/Mexp optimally fit to the number of microRNA families (mirFam), that is robust to variation in mass and sensitive to body temperature (Tb), consistent with the heat dissipation limit theory of Speakman and Król (2010). Using mirFam as predictor, an Ornstein-Uhlenbeck process of stabilizing selection, with an adaptive shift at the divergence of Boreoeutheria, accounted for 95% of the variation in cMR across mammals. Branchwise rates of evolution of cMR, mirFam and Tb concurrently increased 6- to 7-fold at the divergence of Boreoeutheria, independent of mass. Cellular MR variation across placental mammals was also predicted by the sum of model conserved microRNA-target interactions, revealing an unexpected degree of integration of the microRNA-target apparatus into the energy economy of the mammalian cell.
Collapse
Affiliation(s)
- Bastian Fromm
- The Arctic University Museum of Norway, UiT- The Arctic University of Norway, Tromsø, Norway
| | - Thomas Sorger
- Department of Biology, Roger Williams University, Bristol, RI 02809, USA
| |
Collapse
|
13
|
Liao Y. Emerging tools for uncovering genetic and transcriptomic heterogeneities in bacteria. Biophys Rev 2024; 16:109-124. [PMID: 38495445 PMCID: PMC10937887 DOI: 10.1007/s12551-023-01178-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2023] [Accepted: 12/11/2023] [Indexed: 03/19/2024] Open
Abstract
Bacterial communities display an astonishing degree of heterogeneities among their constituent cells across both the genomic and transcriptomic levels, giving rise to diverse social interactions and stress-adaptation strategies indispensable for proliferating in the natural environment (Ackermann in Nat Rev Microbiol 13:497-508, 2015). Our knowledge about bacterial heterogeneities and their physiological ramifications critically depends on our ability to unambiguously resolve the genetic and phenotypic states of the individual cells that make up the population. In this short review, I highlight several recently developed methods for studying bacterial heterogeneities, primarily focusing on single-cell techniques based on advanced sequencing and microscopy technologies. I will discuss the working principle of each technique as well as the types of problems each technique is best positioned to address. With significant improvements in resolution and throughput, these emerging tools together offer unprecedented and complementary views of various types of heterogeneities found within bacterial populations, paving the way for mechanistic dissections and systematic interventions in laboratory and clinical settings.
Collapse
Affiliation(s)
- Yi Liao
- Division of Life Science, Hong Kong University of Science and Technology, Hong Kong SAR, China
| |
Collapse
|
14
|
Khorasani N, Sadeghi M. A computational model of stem cells' internal mechanism to recapitulate spatial patterning and maintain the self-organized pattern in the homeostasis state. Sci Rep 2024; 14:1528. [PMID: 38233402 PMCID: PMC10794714 DOI: 10.1038/s41598-024-51386-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2023] [Accepted: 01/04/2024] [Indexed: 01/19/2024] Open
Abstract
The complex functioning of multi-cellular tissue development relies on proper cell production rates to replace dead or differentiated specialized cells. Stem cells are critical for tissue development and maintenance, as they produce specialized cells to meet the tissues' demands. In this study, we propose a computational model to investigate the stem cell's mechanism, which generates the appropriate proportion of specialized cells, and distributes them to their correct position to form and maintain the organized structure in the population through intercellular reactions. Our computational model focuses on early development, where the populations overall behavior is determined by stem cells and signaling molecules. The model does not include complicated factors such as movement of specialized cells or outside signaling sources. The results indicate that in our model, the stem cells can organize the population into a desired spatial pattern, which demonstrates their ability to self-organize as long as the corresponding leading signal is present. We also investigate the impact of stochasticity, which provides desired non-genetic diversity; however, it can also break the proper boundaries of the desired spatial pattern. We further examine the role of the death rate in maintaining the system's steady state. Overall, our study sheds light on the strategies employed by stem cells to organize specialized cells and maintain proper functionality. Our findings provide insight into the complex mechanisms involved in tissue development and maintenance, which could lead to new approaches in regenerative medicine and tissue engineering.
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
|
15
|
Nakamura YT, Himeoka Y, Saito N, Furusawa C. Evolution of hierarchy and irreversibility in theoretical cell differentiation model. PNAS NEXUS 2024; 3:pgad454. [PMID: 38205032 PMCID: PMC10776358 DOI: 10.1093/pnasnexus/pgad454] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/03/2023] [Accepted: 12/18/2023] [Indexed: 01/12/2024]
Abstract
The process of cell differentiation in multicellular organisms is characterized by hierarchy and irreversibility in many cases. However, the conditions and selection pressures that give rise to these characteristics remain poorly understood. By using a mathematical model, here we show that the network of differentiation potency (differentiation diagram) becomes necessarily hierarchical and irreversible by increasing the number of terminally differentiated states under certain conditions. The mechanisms generating these characteristics are clarified using geometry in the cell state space. The results demonstrate that the hierarchical organization and irreversibility can manifest independently of direct selection pressures associated with these characteristics, instead they appear to evolve as byproducts of selective forces favoring a diversity of differentiated cell types. The study also provides a new perspective on the structure of gene regulatory networks that produce hierarchical and irreversible differentiation diagrams. These results indicate some constraints on cell differentiation, which are expected to provide a starting point for theoretical discussion of the implicit limits and directions of evolution in multicellular organisms.
Collapse
Affiliation(s)
- Yoshiyuki T Nakamura
- Department of Physics, The University of Tokyo, Bunkyo-ku 113-0033, Japan
- Universal Biology Institute, The University of Tokyo, Bunkyo-ku 113-0033, Japan
- Center for Biosystems Dynamics Research, RIKEN, Suita 565-0874, Japan
| | - Yusuke Himeoka
- Universal Biology Institute, The University of Tokyo, Bunkyo-ku 113-0033, Japan
| | - Nen Saito
- Graduate School of Integrated Sciences for Life, Hiroshima University, Higashihiroshima 739-8526, Japan
- Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences, Okazaki 444-8787, Japan
| | - Chikara Furusawa
- Department of Physics, The University of Tokyo, Bunkyo-ku 113-0033, Japan
- Universal Biology Institute, The University of Tokyo, Bunkyo-ku 113-0033, Japan
- Center for Biosystems Dynamics Research, RIKEN, Suita 565-0874, Japan
| |
Collapse
|
16
|
Mellini M, Letizia M, Caruso L, Guiducci A, Meneghini C, Heeb S, Williams P, Cámara M, Visca P, Imperi F, Leoni L, Rampioni G. RsaL-driven negative regulation promotes heterogeneity in Pseudomonas aeruginosa quorum sensing. mBio 2023; 14:e0203923. [PMID: 37843294 PMCID: PMC10746200 DOI: 10.1128/mbio.02039-23] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2023] [Accepted: 09/05/2023] [Indexed: 10/17/2023] Open
Abstract
IMPORTANCE Single-cell analyses can reveal that despite experiencing identical physico-chemical conditions, individual bacterial cells within a monoclonal population may exhibit variations in gene expression. Such phenotypic heterogeneity has been described for several aspects of bacterial physiology, including QS activation. This study demonstrates that the transition of non-quorate cells to the quorate state is a graded process that does not occur at a specific cell density and that subpopulations of non-quorate cells also persist at high cell density. Here, we provide a mechanistic explanation for this phenomenon, showing that a negative feedback regulatory loop integrated into the las system has a pivotal role in promoting cell-to-cell variation in the QS activation state and in limiting the transition of non-quorate cells to the quorate state in P. aeruginosa.
Collapse
Affiliation(s)
- Marta Mellini
- Department of Science, University Roma Tre, Rome, Italy
| | | | | | | | | | - Stephan Heeb
- National Biofilms Innovation Centre, Biodiscovery Institute and School of Life Sciences, University of Nottingham, Nottingham, United Kingdom
| | - Paul Williams
- National Biofilms Innovation Centre, Biodiscovery Institute and School of Life Sciences, University of Nottingham, Nottingham, United Kingdom
| | - Miguel Cámara
- National Biofilms Innovation Centre, Biodiscovery Institute and School of Life Sciences, University of Nottingham, Nottingham, United Kingdom
| | - Paolo Visca
- Department of Science, University Roma Tre, Rome, Italy
- NBFC, National Biodiversity Future Center, Palermo, Italy
- IRCCS Fondazione Santa Lucia, Rome, Italy
| | - Francesco Imperi
- Department of Science, University Roma Tre, Rome, Italy
- NBFC, National Biodiversity Future Center, Palermo, Italy
- IRCCS Fondazione Santa Lucia, Rome, Italy
| | - Livia Leoni
- Department of Science, University Roma Tre, Rome, Italy
| | - Giordano Rampioni
- Department of Science, University Roma Tre, Rome, Italy
- IRCCS Fondazione Santa Lucia, Rome, Italy
| |
Collapse
|
17
|
Feng H, Li F, Wang T, Xing XH, Zeng AP, Zhang C. Deep-learning-assisted Sort-Seq enables high-throughput profiling of gene expression characteristics with high precision. SCIENCE ADVANCES 2023; 9:eadg5296. [PMID: 37939173 PMCID: PMC10631719 DOI: 10.1126/sciadv.adg5296] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/04/2023] [Accepted: 10/06/2023] [Indexed: 11/10/2023]
Abstract
Owing to the nondeterministic and nonlinear nature of gene expression, the steady-state intracellular protein abundance of a clonal population forms a distribution. The characteristics of this distribution, including expression strength and noise, are closely related to cellular behavior. However, quantitative description of these characteristics has so far relied on arrayed methods, which are time-consuming and labor-intensive. To address this issue, we propose a deep-learning-assisted Sort-Seq approach (dSort-Seq) in this work, enabling high-throughput profiling of expression properties with high precision. We demonstrated the validity of dSort-Seq for large-scale assaying of the dose-response relationships of biosensors. In addition, we comprehensively investigated the contribution of transcription and translation to noise production in Escherichia coli, from which we found that the expression noise is strongly coupled with the mean expression level. We also found that the transcriptional interference caused by overlapping RpoD-binding sites contributes to noise production, which suggested the existence of a simple and feasible noise control strategy in E. coli.
Collapse
Affiliation(s)
- Huibao Feng
- MOE Key Laboratory for Industrial Biocatalysis, Institute of Biochemical Engineering, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Fan Li
- MOE Key Laboratory for Industrial Biocatalysis, Institute of Biochemical Engineering, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
| | - Tianmin Wang
- Tsinghua-Peking Center for Life Sciences, School of Medicine, Tsinghua University, Beijing 100084, China
- School of Life Science and Technology, ShanghaiTech University, Shanghai 201210, China
| | - Xin-hui Xing
- MOE Key Laboratory for Industrial Biocatalysis, Institute of Biochemical Engineering, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
- Shenzhen International Graduate School, Tsinghua University, Shenzhen 518055, China
| | - An-ping Zeng
- Institute of Bioprocess and Biosystems Engineering, Hamburg University of Technology, Hamburg 21073, Germany
- Center of Synthetic Biology and Integrated Bioengineering, School of Engineering, Westlake University, Hangzhou 310024, China
| | - Chong Zhang
- MOE Key Laboratory for Industrial Biocatalysis, Institute of Biochemical Engineering, Department of Chemical Engineering, Tsinghua University, Beijing 100084, China
- Center for Synthetic and Systems Biology, Tsinghua University, Beijing 100084, China
| |
Collapse
|
18
|
Keseroglu K, Zinani OQH, Keskin S, Seawall H, Alpay EE, Özbudak EM. Stochastic gene expression and environmental stressors trigger variable somite segmentation phenotypes. Nat Commun 2023; 14:6497. [PMID: 37838784 PMCID: PMC10576776 DOI: 10.1038/s41467-023-42220-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2023] [Accepted: 10/04/2023] [Indexed: 10/16/2023] Open
Abstract
Mutations of several genes cause incomplete penetrance and variable expressivity of phenotypes, which are usually attributed to modifier genes or gene-environment interactions. Here, we show stochastic gene expression underlies the variability of somite segmentation defects in embryos mutant for segmentation clock genes her1 or her7. Phenotypic strength is further augmented by low temperature and hypoxia. By performing live imaging of the segmentation clock reporters, we further show that groups of cells with higher oscillation amplitudes successfully form somites while those with lower amplitudes fail to do so. In unfavorable environments, the number of cycles with high amplitude oscillations and the number of successful segmentations proportionally decrease. These results suggest that individual oscillation cycles stochastically fail to pass a threshold amplitude, resulting in segmentation defects in mutants. Our quantitative methodology is adaptable to investigate variable phenotypes of mutant genes in different tissues.
Collapse
Affiliation(s)
- Kemal Keseroglu
- Division of Developmental Biology, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, 45229, USA
| | - Oriana Q H Zinani
- Division of Developmental Biology, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, 45229, USA
- Molecular and Developmental Biology Graduate Program, University of Cincinnati, College of Medicine, Cincinnati, OH, 45229, USA
| | - Sevdenur Keskin
- Allergy and Immunology, University of Arkansas for Medical Science and Arkansas Children's Hospital, Little Rock, AR, 72202, USA
| | - Hannah Seawall
- Division of Developmental Biology, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, 45229, USA
| | - Eslim E Alpay
- Division of Developmental Biology, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, 45229, USA
| | - Ertuğrul M Özbudak
- Division of Developmental Biology, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, 45229, USA.
- Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati, OH, 45229, USA.
| |
Collapse
|
19
|
Lewin LE, Daniels KG, Hurst LD. Genes for highly abundant proteins in Escherichia coli avoid 5' codons that promote ribosomal initiation. PLoS Comput Biol 2023; 19:e1011581. [PMID: 37878567 PMCID: PMC10599525 DOI: 10.1371/journal.pcbi.1011581] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2023] [Accepted: 10/09/2023] [Indexed: 10/27/2023] Open
Abstract
In many species highly expressed genes (HEGs) over-employ the synonymous codons that match the more abundant iso-acceptor tRNAs. Bacterial transgene codon randomization experiments report, however, that enrichment with such "translationally optimal" codons has little to no effect on the resultant protein level. By contrast, consistent with the view that ribosomal initiation is rate limiting, synonymous codon usage following the 5' ATG greatly influences protein levels, at least in part by modifying RNA stability. For the design of bacterial transgenes, for simple codon based in silico inference of protein levels and for understanding selection on synonymous mutations, it would be valuable to computationally determine initiation optimality (IO) scores for codons for any given species. One attractive approach is to characterize the 5' codon enrichment of HEGs compared with the most lowly expressed genes, just as translational optimality scores of codons have been similarly defined employing the full gene body. Here we determine the viability of this approach employing a unique opportunity: for Escherichia coli there is both the most extensive protein abundance data for native genes and a unique large-scale transgene codon randomization experiment enabling objective definition of the 5' codons that cause, rather than just correlate with, high protein abundance (that we equate with initiation optimality, broadly defined). Surprisingly, the 5' ends of native genes that specify highly abundant proteins avoid such initiation optimal codons. We find that this is probably owing to conflicting selection pressures particular to native HEGs, including selection favouring low initiation rates, this potentially enabling high efficiency of ribosomal usage and low noise. While the classical HEG enrichment approach does not work, rendering simple prediction of native protein abundance from 5' codon content futile, we report evidence that initiation optimality scores derived from the transgene experiment may hold relevance for in silico transgene design for a broad spectrum of bacteria.
Collapse
Affiliation(s)
- Loveday E. Lewin
- The Milner Centre for Evolution, Department of Life Sciences, University of Bath, Bath, United Kingdom
| | - Kate G. Daniels
- The Milner Centre for Evolution, Department of Life Sciences, University of Bath, Bath, United Kingdom
| | - Laurence D. Hurst
- The Milner Centre for Evolution, Department of Life Sciences, University of Bath, Bath, United Kingdom
| |
Collapse
|
20
|
Stevanovic M, Carvalho JPT, Bittihn P, Schultz D. Dynamical model of antibiotic responses linking expression of resistance to metabolism explains emergence of heterogeneity during drug exposures. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.09.22.558994. [PMID: 37790326 PMCID: PMC10542528 DOI: 10.1101/2023.09.22.558994] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
Antibiotic responses in bacteria are highly dynamic and heterogeneous, with sudden exposure of bacterial colonies to high drug doses resulting in the coexistence of recovered and arrested cells. The dynamics of the response is determined by regulatory circuits controlling the expression of resistance genes, which are in turn modulated by the drug's action on cell growth and metabolism. Despite advances in understanding gene regulation at the molecular level, we still lack a framework to describe how feedback mechanisms resulting from the interdependence between expression of resistance and cell metabolism can amplify naturally occurring noise and create heterogeneity at the population level. To understand how this interplay affects cell survival upon exposure, we constructed a mathematical model of the dynamics of antibiotic responses that links metabolism and regulation of gene expression, based on the tetracycline resistance tet operon in E. coli. We use this model to interpret measurements of growth and expression of resistance in microfluidic experiments, both in single cells and in biofilms. We also implemented a stochastic model of the drug response, to show that exposure to high drug levels results in large variations of recovery times and heterogeneity at the population level. We show that stochasticity is important to determine how nutrient quality affects cell survival during exposure to high drug concentrations. A quantitative description of how microbes respond to antibiotics in dynamical environments is crucial to understand population-level behaviors such as biofilms and pathogenesis.
Collapse
Affiliation(s)
- Mirjana Stevanovic
- Department of Microbiology and Immunology, Geisel School of Medicine at Dartmouth, Hanover, NH, United States
| | - João Pedro Teuber Carvalho
- Department of Microbiology and Immunology, Geisel School of Medicine at Dartmouth, Hanover, NH, United States
| | - Philip Bittihn
- Max Planck Institute for Dynamics and Self-Organization, Göttingen, Germany
- Institute for the Dynamics of Complex Systems, University of Göttingen, Göttingen, Germany
| | - Daniel Schultz
- Department of Microbiology and Immunology, Geisel School of Medicine at Dartmouth, Hanover, NH, United States
| |
Collapse
|
21
|
Zhang J, Han X, Ma L, Xu S, Lin Y. Deciphering a global source of non-genetic heterogeneity in cancer cells. Nucleic Acids Res 2023; 51:9019-9038. [PMID: 37587722 PMCID: PMC10516630 DOI: 10.1093/nar/gkad666] [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: 09/13/2022] [Revised: 07/09/2023] [Accepted: 08/09/2023] [Indexed: 08/18/2023] Open
Abstract
Cell-to-cell variability within a clonal population, also known as non-genetic heterogeneity, has created significant challenges for intervening with diseases such as cancer. While non-genetic heterogeneity can arise from the variability in the expression of specific genes, it remains largely unclear whether and how clonal cells could be heterogeneous in the expression of the entire transcriptome. Here, we showed that gene transcriptional activity is globally modulated in individual cancer cells, leading to non-genetic heterogeneity in the global transcription rate. Such heterogeneity contributes to cell-to-cell variability in transcriptome size and displays both dynamic and static characteristics, with the global transcription rate temporally modulated in a cell-cycle-coupled manner and the time-averaged rate being distinct between cells and heritable across generations. Additional evidence indicated the role of ATP metabolism in this heterogeneity, and suggested its implication in intrinsic cancer drug tolerance. Collectively, our work shed light on the mode, mechanism, and implication of a global but often hidden source of non-genetic heterogeneity.
Collapse
Affiliation(s)
- Jianhan Zhang
- Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- The MOE Key Laboratory of Cell Proliferation and Differentiation, School of Life Sciences, Peking University, Beijing 100871, China
| | - Xu Han
- Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- The MOE Key Laboratory of Cell Proliferation and Differentiation, School of Life Sciences, Peking University, Beijing 100871, China
| | - Liang Ma
- Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- The MOE Key Laboratory of Cell Proliferation and Differentiation, School of Life Sciences, Peking University, Beijing 100871, China
| | - Shuhui Xu
- Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- The MOE Key Laboratory of Cell Proliferation and Differentiation, School of Life Sciences, Peking University, Beijing 100871, China
| | - Yihan Lin
- Center for Quantitative Biology and Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- The MOE Key Laboratory of Cell Proliferation and Differentiation, School of Life Sciences, Peking University, Beijing 100871, China
| |
Collapse
|
22
|
Damour A, Slaninova V, Radulescu O, Bertrand E, Basyuk E. Transcriptional Stochasticity as a Key Aspect of HIV-1 Latency. Viruses 2023; 15:1969. [PMID: 37766375 PMCID: PMC10535884 DOI: 10.3390/v15091969] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2023] [Revised: 09/16/2023] [Accepted: 09/18/2023] [Indexed: 09/29/2023] Open
Abstract
This review summarizes current advances in the role of transcriptional stochasticity in HIV-1 latency, which were possible in a large part due to the development of single-cell approaches. HIV-1 transcription proceeds in bursts of RNA production, which stem from the stochastic switching of the viral promoter between ON and OFF states. This switching is caused by random binding dynamics of transcription factors and nucleosomes to the viral promoter and occurs at several time scales from minutes to hours. Transcriptional bursts are mainly controlled by the core transcription factors TBP, SP1 and NF-κb, the chromatin status of the viral promoter and RNA polymerase II pausing. In particular, spontaneous variability in the promoter chromatin creates heterogeneity in the response to activators such as TNF-α, which is then amplified by the Tat feedback loop to generate high and low viral transcriptional states. This phenomenon is likely at the basis of the partial and stochastic response of latent T cells from HIV-1 patients to latency-reversing agents, which is a barrier for the development of shock-and-kill strategies of viral eradication. A detailed understanding of the transcriptional stochasticity of HIV-1 and the possibility to precisely model this phenomenon will be important assets to develop more effective therapeutic strategies.
Collapse
Affiliation(s)
- Alexia Damour
- MFP UMR 5234 CNRS, Université de Bordeaux, 33076 Bordeaux, France;
| | - Vera Slaninova
- IGH UMR 9002 CNRS, Université de Montpellier, 34094 Montpellier, France;
| | - Ovidiu Radulescu
- LPHI, UMR 5294 CNRS, University of Montpellier, 34095 Montpellier, France;
| | - Edouard Bertrand
- IGH UMR 9002 CNRS, Université de Montpellier, 34094 Montpellier, France;
| | - Eugenia Basyuk
- MFP UMR 5234 CNRS, Université de Bordeaux, 33076 Bordeaux, France;
| |
Collapse
|
23
|
Kell B, Ripsman R, Hilfinger A. Noise properties of adaptation-conferring biochemical control modules. Proc Natl Acad Sci U S A 2023; 120:e2302016120. [PMID: 37695915 PMCID: PMC10515136 DOI: 10.1073/pnas.2302016120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2023] [Accepted: 06/12/2023] [Indexed: 09/13/2023] Open
Abstract
A key goal of synthetic biology is to develop functional biochemical modules with network-independent properties. Antithetic integral feedback (AIF) is a recently developed control module in which two control species perfectly annihilate each other's biological activity. The AIF module confers robust perfect adaptation to the steady-state average level of a controlled intracellular component when subjected to sustained perturbations. Recent work has suggested that such robustness comes at the unavoidable price of increased stochastic fluctuations around average levels. We present theoretical results that support and quantify this trade-off for the commonly analyzed AIF variant in the idealized limit with perfect annihilation. However, we also show that this trade-off is a singular limit of the control module: Even minute deviations from perfect adaptation allow systems to achieve effective noise suppression as long as cells can pay the corresponding energetic cost. We further show that a variant of the AIF control module can achieve significant noise suppression even in the idealized limit with perfect adaptation. This atypical configuration may thus be preferable in synthetic biology applications.
Collapse
Affiliation(s)
- Brayden Kell
- Department of Physics, University of Toronto, Toronto, ONM5S 1A7, Canada
- Department of Chemical and Physical Sciences, University of Toronto, Mississauga, ONL5L 1C6, Canada
- Department of Molecular Biosciences, Northwestern University, Evanston, IL60208
- National Science Foundation-Simons Center for Quantitative Biology, Northwestern University, Evanston, IL60208
| | - Ryan Ripsman
- Department of Physics, University of Toronto, Toronto, ONM5S 1A7, Canada
- Department of Laboratory Medicine and Pathobiology, University of Toronto, Toronto, ONM5S 1A8, Canada
| | - Andreas Hilfinger
- Department of Physics, University of Toronto, Toronto, ONM5S 1A7, Canada
- Department of Chemical and Physical Sciences, University of Toronto, Mississauga, ONL5L 1C6, Canada
- Department of Mathematics, University of Toronto, Toronto, ONM5S 2E4, Canada
- Department of Cell and Systems Biology, University of Toronto, Toronto, ONM5S 3G5, Canada
| |
Collapse
|
24
|
Douaihy M, Topno R, Lagha M, Bertrand E, Radulescu O. BurstDECONV: a signal deconvolution method to uncover mechanisms of transcriptional bursting in live cells. Nucleic Acids Res 2023; 51:e88. [PMID: 37522372 PMCID: PMC10484743 DOI: 10.1093/nar/gkad629] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2023] [Revised: 06/23/2023] [Accepted: 07/17/2023] [Indexed: 08/01/2023] Open
Abstract
Monitoring transcription in living cells gives access to the dynamics of this complex fundamental process. It reveals that transcription is discontinuous, whereby active periods (bursts) are separated by one or several types of inactive periods of distinct lifetimes. However, decoding temporal fluctuations arising from live imaging and inferring the distinct transcriptional steps eliciting them is a challenge. We present BurstDECONV, a novel statistical inference method that deconvolves signal traces into individual transcription initiation events. We use the distribution of waiting times between successive polymerase initiation events to identify mechanistic features of transcription such as the number of rate-limiting steps and their kinetics. Comparison of our method to alternative methods emphasizes its advantages in terms of precision and flexibility. Unique features such as the direct determination of the number of promoter states and the simultaneous analysis of several potential transcription models make BurstDECONV an ideal analytic framework for live cell transcription imaging experiments. Using simulated realistic data, we found that our method is robust with regards to noise or suboptimal experimental designs. To show its generality, we applied it to different biological contexts such as Drosophila embryos or human cells.
Collapse
Affiliation(s)
- Maria Douaihy
- LPHI, University of Montpellier and CNRS, Place Eugène Bataillon, Montpellier 34095, France
- IGMM, University of Montpellier and CNRS, 1919 Rte de Mende, Montpellier 34090, France
| | - Rachel Topno
- LPHI, University of Montpellier and CNRS, Place Eugène Bataillon, Montpellier 34095, France
- IGH, University of Montpellier and CNRS, 141 Rue de la Cardonille, Montpellier 34094, France
| | - Mounia Lagha
- IGMM, University of Montpellier and CNRS, 1919 Rte de Mende, Montpellier 34090, France
| | - Edouard Bertrand
- IGH, University of Montpellier and CNRS, 141 Rue de la Cardonille, Montpellier 34094, France
| | - Ovidiu Radulescu
- LPHI, University of Montpellier and CNRS, Place Eugène Bataillon, Montpellier 34095, France
| |
Collapse
|
25
|
Haerinck J, Goossens S, Berx G. The epithelial-mesenchymal plasticity landscape: principles of design and mechanisms of regulation. Nat Rev Genet 2023; 24:590-609. [PMID: 37169858 DOI: 10.1038/s41576-023-00601-0] [Citation(s) in RCA: 18] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/30/2023] [Indexed: 05/13/2023]
Abstract
Epithelial-mesenchymal plasticity (EMP) enables cells to interconvert between several states across the epithelial-mesenchymal landscape, thereby acquiring hybrid epithelial/mesenchymal phenotypic features. This plasticity is crucial for embryonic development and wound healing, but also underlies the acquisition of several malignant traits during cancer progression. Recent research using systems biology and single-cell profiling methods has provided novel insights into the main forces that shape EMP, which include the microenvironment, lineage specification and cell identity, and the genome. Additionally, key roles have emerged for hysteresis (cell memory) and cellular noise, which can drive stochastic transitions between cell states. Here, we review these forces and the distinct but interwoven layers of regulatory control that stabilize EMP states or facilitate epithelial-mesenchymal transitions (EMTs) and discuss the therapeutic potential of manipulating the EMP landscape.
Collapse
Affiliation(s)
- Jef Haerinck
- Molecular and Cellular Oncology Laboratory, Department of Biomedical Molecular Biology, Ghent University, Ghent, Belgium
- Cancer Research Institute Ghent (CRIG), Ghent, Belgium
| | - Steven Goossens
- Cancer Research Institute Ghent (CRIG), Ghent, Belgium
- Unit for Translational Research in Oncology, Department of Diagnostic Sciences, Ghent University, Ghent, Belgium
| | - Geert Berx
- Molecular and Cellular Oncology Laboratory, Department of Biomedical Molecular Biology, Ghent University, Ghent, Belgium.
- Cancer Research Institute Ghent (CRIG), Ghent, Belgium.
| |
Collapse
|
26
|
Das S, Singh A, Shah P. Evaluating single-cell variability in proteasomal decay. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.08.22.554358. [PMID: 37662347 PMCID: PMC10473619 DOI: 10.1101/2023.08.22.554358] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/05/2023]
Abstract
Gene expression is a stochastic process that leads to variability in mRNA and protein abundances even within an isogenic population of cells grown in the same environment. This variation, often called gene-expression noise, has typically been attributed to transcriptional and translational processes while ignoring the contributions of protein decay variability across cells. Here we estimate the single-cell protein decay rates of two degron GFPs in Saccharomyces cerevisiae using time-lapse microscopy. We find substantial cell-to-cell variability in the decay rates of the degron GFPs. We evaluate cellular features that explain the variability in the proteasomal decay and find that the amount of 20s catalytic beta subunit of the proteasome marginally explains the observed variability in the degron GFP half-lives. We propose alternate hypotheses that might explain the observed variability in the decay of the two degron GFPs. Overall, our study highlights the importance of studying the kinetics of the decay process at single-cell resolution and that decay rates vary at the single-cell level, and that the decay process is stochastic. A complex model of decay dynamics must be included when modeling stochastic gene expression to estimate gene expression noise.
Collapse
Affiliation(s)
| | - Abhyudai Singh
- Department of Electrical and Computer Engineering, Biomedical Engineering, University of Delaware
| | | |
Collapse
|
27
|
Gholami S, Ilie S. Quantifying Parameter Interdependence in Stochastic Discrete Models of Biochemical Systems. ENTROPY (BASEL, SWITZERLAND) 2023; 25:1168. [PMID: 37628198 PMCID: PMC10452982 DOI: 10.3390/e25081168] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/03/2023] [Revised: 07/31/2023] [Accepted: 08/03/2023] [Indexed: 08/27/2023]
Abstract
Stochastic modeling of biochemical processes at the cellular level has been the subject of intense research in recent years. The Chemical Master Equation is a broadly utilized stochastic discrete model of such processes. Numerous important biochemical systems consist of many species subject to many reactions. As a result, their mathematical models depend on many parameters. In applications, some of the model parameters may be unknown, so their values need to be estimated from the experimental data. However, the problem of parameter value inference can be quite challenging, especially in the stochastic setting. To estimate accurately the values of a subset of parameters, the system should be sensitive with respect to variations in each of these parameters and they should not be correlated. In this paper, we propose a technique for detecting collinearity among models' parameters and we apply this method for selecting subsets of parameters that can be estimated from the available data. The analysis relies on finite-difference sensitivity estimations and the singular value decomposition of the sensitivity matrix. We illustrated the advantages of the proposed method by successfully testing it on several models of biochemical systems of practical interest.
Collapse
Affiliation(s)
- Samaneh Gholami
- Department of Mathematics, Toronto Metropolitan University, Toronto, ON M5B 2K3, Canada;
| | | |
Collapse
|
28
|
Torres A, Cockerell S, Phillips M, Balázsi G, Ghosh K. MaxCal can infer models from coupled stochastic trajectories of gene expression and cell division. Biophys J 2023; 122:2623-2635. [PMID: 37218129 PMCID: PMC10397576 DOI: 10.1016/j.bpj.2023.05.017] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2023] [Revised: 05/03/2023] [Accepted: 05/18/2023] [Indexed: 05/24/2023] Open
Abstract
Gene expression is inherently noisy due to small numbers of proteins and nucleic acids inside a cell. Likewise, cell division is stochastic, particularly when tracking at the level of a single cell. The two can be coupled when gene expression affects the rate of cell division. Single-cell time-lapse experiments can measure both fluctuations by simultaneously recording protein levels inside a cell and its stochastic division. These information-rich noisy trajectory data sets can be harnessed to learn about the underlying molecular and cellular details that are often not known a priori. A critical question is: How can we infer a model given data where fluctuations at two levels-gene expression and cell division-are intricately convoluted? We show the principle of maximum caliber (MaxCal)-integrated within a Bayesian framework-can be used to infer several cellular and molecular details (division rates, protein production, and degradation rates) from these coupled stochastic trajectories (CSTs). We demonstrate this proof of concept using synthetic data generated from a known model. An additional challenge in data analysis is that trajectories are often not in protein numbers, but in noisy fluorescence that depends on protein number in a probabilistic manner. We again show that MaxCal can infer important molecular and cellular rates even when data are in fluorescence, another example of CST with three confounding factors-gene expression noise, cell division noise, and fluorescence distortion-all coupled. Our approach will provide guidance to build models in synthetic biology experiments as well as general biological systems where examples of CSTs are abundant.
Collapse
Affiliation(s)
- Andrew Torres
- Department of Physics and Astronomy, University of Denver, Denver, Colorado
| | - Spencer Cockerell
- Department of Physics and Astronomy, University of Denver, Denver, Colorado
| | - Michael Phillips
- Department of Physics and Astronomy, University of Denver, Denver, Colorado
| | - Gábor Balázsi
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York
| | - Kingshuk Ghosh
- Molecular and Cellular Biophysics, University of Denver, Denver, Colorado; Department of Physics and Astronomy, University of Denver, Denver, Colorado.
| |
Collapse
|
29
|
Ramu A, Cohen BA. Transcription factor fluctuations underlie cell-to-cell variability in a signaling pathway response. Genetics 2023; 224:iyad094. [PMID: 37226217 PMCID: PMC10691749 DOI: 10.1093/genetics/iyad094] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2023] [Revised: 04/25/2023] [Accepted: 05/05/2023] [Indexed: 05/26/2023] Open
Abstract
Stochastic differences among clonal cells can initiate cell fate decisions in development or cause cell-to-cell differences in the responses to drugs or extracellular ligands. One hypothesis is that some of this phenotypic variability is caused by stochastic fluctuations in the activities of transcription factors (TFs). We tested this hypothesis in NIH3T3-CG cells using the response to Hedgehog signaling as a model cellular response. Here, we present evidence for the existence of distinct fast- and slow-responding substates in NIH3T3-CG cells. These two substates have distinct expression profiles, and fluctuations in the Prrx1 TF underlie some of the differences in expression and responsiveness between fast and slow cells. Our results show that fluctuations in TFs can contribute to cell-to-cell differences in Hedgehog signaling.
Collapse
Affiliation(s)
- Avinash Ramu
- The Edison Family Center for Genome Sciences and Systems Biology, School of Medicine, Washington University in St. Louis, Saint Louis, MO 63110, USA
- Department of Genetics, School of Medicine, Washington University in St. Louis, Saint Louis, MO 63110, USA
| | - Barak A Cohen
- The Edison Family Center for Genome Sciences and Systems Biology, School of Medicine, Washington University in St. Louis, Saint Louis, MO 63110, USA
- Department of Genetics, School of Medicine, Washington University in St. Louis, Saint Louis, MO 63110, USA
| |
Collapse
|
30
|
Rosales-Alvarez RE, Rettkowski J, Herman JS, Dumbović G, Cabezas-Wallscheid N, Grün D. VarID2 quantifies gene expression noise dynamics and unveils functional heterogeneity of ageing hematopoietic stem cells. Genome Biol 2023; 24:148. [PMID: 37353813 PMCID: PMC10290360 DOI: 10.1186/s13059-023-02974-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2022] [Accepted: 05/18/2023] [Indexed: 06/25/2023] Open
Abstract
Variability of gene expression due to stochasticity of transcription or variation of extrinsic signals, termed biological noise, is a potential driving force of cellular differentiation. Utilizing single-cell RNA-sequencing, we develop VarID2 for the quantification of biological noise at single-cell resolution. VarID2 reveals enhanced nuclear versus cytoplasmic noise, and distinct regulatory modes stratified by correlation between noise, expression, and chromatin accessibility. Noise levels are minimal in murine hematopoietic stem cells (HSCs) and increase during differentiation and ageing. Differential noise identifies myeloid-biased Dlk1+ long-term HSCs in aged mice with enhanced quiescence and self-renewal capacity. VarID2 reveals noise dynamics invisible to conventional single-cell transcriptome analysis.
Collapse
Affiliation(s)
- Reyna Edith Rosales-Alvarez
- Würzburg Institute of Systems Immunology, Max Planck Research Group at the Julius-Maximilians-Universität Würzburg, Würzburg, Germany
- International Max Planck Research School for Immunobiology, Epigenetics, and Metabolism (IMPRS-IEM), Freiburg, Germany
- Faculty of Biology, University of Freiburg, Freiburg, Germany
| | - Jasmin Rettkowski
- Max Planck Institute of Immunobiology and Epigenetics, Freiburg, Germany
- Spemann Graduate School of Biology and Medicine (SGBM), Freiburg, Germany
| | - Josip Stefan Herman
- Würzburg Institute of Systems Immunology, Max Planck Research Group at the Julius-Maximilians-Universität Würzburg, Würzburg, Germany
| | - Gabrijela Dumbović
- Würzburg Institute of Systems Immunology, Max Planck Research Group at the Julius-Maximilians-Universität Würzburg, Würzburg, Germany
| | - Nina Cabezas-Wallscheid
- Max Planck Institute of Immunobiology and Epigenetics, Freiburg, Germany
- CIBSS-Centre for Integrative Biological Signaling Studies, University of Freiburg, Freiburg, Germany
| | - Dominic Grün
- Würzburg Institute of Systems Immunology, Max Planck Research Group at the Julius-Maximilians-Universität Würzburg, Würzburg, Germany.
- Helmholtz Institute for RNA-Based Infection Research (HIRI), Helmholtz-Center for Infection Research (HZI), Würzburg, Germany.
| |
Collapse
|
31
|
Tobler S, Köhler K, Sinha T, Hafen E, Kapur M. Understanding Randomness on a Molecular Level: A Diagnostic Tool. CBE LIFE SCIENCES EDUCATION 2023; 22:ar17. [PMID: 36862800 PMCID: PMC10228260 DOI: 10.1187/cbe.22-05-0097] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/23/2022] [Revised: 01/24/2023] [Accepted: 02/03/2023] [Indexed: 06/02/2023]
Abstract
Undergraduate biology students' molecular-level understanding of stochastic (also referred to as random or noisy) processes found in biological systems is often limited to those examples discussed in class. Therefore, students frequently display little ability to accurately transfer their knowledge to other contexts. Furthermore, elaborate tools to assess students' understanding of these stochastic processes are missing, despite the fundamental nature of this concept and the increasing evidence demonstrating its importance in biology. Thus, we developed the Molecular Randomness Concept Inventory (MRCI), an instrument composed of nine multiple-choice questions based on students' most prevalent misconceptions, to quantify students' understanding of stochastic processes in biological systems. The MRCI was administered to 67 first-year natural science students in Switzerland. The psychometric properties of the inventory were analyzed using classical test theory and Rasch modeling. Moreover, think-aloud interviews were conducted to ensure response validity. Results indicate that the MRCI yields valid and reliable estimations of students' conceptual understanding of molecular randomness in the higher educational setting studied. Ultimately, the performance analysis sheds light on the extent and the limitations of students' understanding of the concept of stochasticity on a molecular level.
Collapse
Affiliation(s)
- Samuel Tobler
- Professorship for Learning Sciences and Higher Education and
| | - Katja Köhler
- Department of Biology, ETH Zurich, 8092 Zurich, Switzerland
| | - Tanmay Sinha
- Professorship for Learning Sciences and Higher Education and
| | - Ernst Hafen
- Department of Biology, ETH Zurich, 8092 Zurich, Switzerland
| | - Manu Kapur
- Professorship for Learning Sciences and Higher Education and
| |
Collapse
|
32
|
Aubé S, Nielly-Thibault L, Landry CR. Evolutionary trade-off and mutational bias could favor transcriptional over translational divergence within paralog pairs. PLoS Genet 2023; 19:e1010756. [PMID: 37235586 DOI: 10.1371/journal.pgen.1010756] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2023] [Accepted: 04/21/2023] [Indexed: 05/28/2023] Open
Abstract
How changes in the different steps of protein synthesis-transcription, translation and degradation-contribute to differences of protein abundance among genes is not fully understood. There is however accumulating evidence that transcriptional divergence might have a prominent role. Here, we show that yeast paralogous genes are more divergent in transcription than in translation. We explore two causal mechanisms for this predominance of transcriptional divergence: an evolutionary trade-off between the precision and economy of gene expression and a larger mutational target size for transcription. Performing simulations within a minimal model of post-duplication evolution, we find that both mechanisms are consistent with the observed divergence patterns. We also investigate how additional properties of the effects of mutations on gene expression, such as their asymmetry and correlation across levels of regulation, can shape the evolution of paralogs. Our results highlight the importance of fully characterizing the distributions of mutational effects on transcription and translation. They also show how general trade-offs in cellular processes and mutation bias can have far-reaching evolutionary impacts.
Collapse
Affiliation(s)
- Simon Aubé
- Département de biochimie, de microbiologie et de bio-informatique, Faculté des sciences et de génie, Université Laval, Québec, Québec, Canada
- Institut de Biologie Intégrative et des Systèmes, Université Laval, Québec, Québec, Canada
- PROTEO, Le regroupement québécois de recherche sur la fonction, l'ingénierie et les applications des protéines, Université Laval, Québec, Québec, Canada
- Centre de Recherche en Données Massives, Université Laval, Québec, Québec, Canada
| | - Lou Nielly-Thibault
- Institut de Biologie Intégrative et des Systèmes, Université Laval, Québec, Québec, Canada
- PROTEO, Le regroupement québécois de recherche sur la fonction, l'ingénierie et les applications des protéines, Université Laval, Québec, Québec, Canada
- Centre de Recherche en Données Massives, Université Laval, Québec, Québec, Canada
- Département de biologie, Faculté des sciences et de génie, Université Laval, Québec, Québec, Canada
| | - Christian R Landry
- Département de biochimie, de microbiologie et de bio-informatique, Faculté des sciences et de génie, Université Laval, Québec, Québec, Canada
- Institut de Biologie Intégrative et des Systèmes, Université Laval, Québec, Québec, Canada
- PROTEO, Le regroupement québécois de recherche sur la fonction, l'ingénierie et les applications des protéines, Université Laval, Québec, Québec, Canada
- Centre de Recherche en Données Massives, Université Laval, Québec, Québec, Canada
- Département de biologie, Faculté des sciences et de génie, Université Laval, Québec, Québec, Canada
| |
Collapse
|
33
|
Gopalan V, Hannenhalli S. Towards a Synthesis of the Non-Genetic and Genetic Views of Cancer in Understanding Pancreatic Ductal Adenocarcinoma Initiation and Prevention. Cancers (Basel) 2023; 15:cancers15072159. [PMID: 37046820 PMCID: PMC10093726 DOI: 10.3390/cancers15072159] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Revised: 03/22/2023] [Accepted: 03/27/2023] [Indexed: 04/14/2023] Open
Abstract
While much of the research in oncogenesis and cancer therapy has focused on mutations in key cancer driver genes, more recent work suggests a complementary non-genetic paradigm. This paradigm focuses on how transcriptional and phenotypic heterogeneity, even in clonally derived cells, can create sub-populations associated with oncogenesis, metastasis, and therapy resistance. We discuss this complementary paradigm in the context of pancreatic ductal adenocarcinoma. A better understanding of cellular transcriptional heterogeneity and its association with oncogenesis can lead to more effective therapies that prevent tumor initiation and slow progression.
Collapse
Affiliation(s)
- Vishaka Gopalan
- Cancer Data Science Laboratory, National Cancer Institute, National Institutes of Health, Bethesda, MD 20814, USA
| | - Sridhar Hannenhalli
- Cancer Data Science Laboratory, National Cancer Institute, National Institutes of Health, Bethesda, MD 20814, USA
| |
Collapse
|
34
|
Abstract
Organ development and homeostasis involve dynamic interactions between individual cells that collectively regulate tissue architecture and function. To ensure the highest tissue fidelity, equally fit cell populations are continuously renewed by stochastic replacement events, while cells perceived as less fit are actively removed by their fitter counterparts. This renewal is mediated by surveillance mechanisms that are collectively known as cell competition. Recent studies have revealed that cell competition has roles in most, if not all, developing and adult tissues. They have also established that cell competition functions both as a tumour-suppressive mechanism and as a tumour-promoting mechanism, thereby critically influencing cancer initiation and development. This Review discusses the latest insights into the mechanisms of cell competition and its different roles during embryonic development, homeostasis and cancer.
Collapse
|
35
|
Polizzi S, Marzi T, Matteuzzi T, Castellani G, Bazzani A. Random Walk Approximation for Stochastic Processes on Graphs. ENTROPY (BASEL, SWITZERLAND) 2023; 25:394. [PMID: 36981283 PMCID: PMC10047168 DOI: 10.3390/e25030394] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/12/2023] [Revised: 02/13/2023] [Accepted: 02/18/2023] [Indexed: 06/18/2023]
Abstract
We introduce the Random Walk Approximation (RWA), a new method to approximate the stationary solution of master equations describing stochastic processes taking place on graphs. Our approximation can be used for all processes governed by non-linear master equations without long-range interactions and with a conserved number of entities, which are typical in biological systems, such as gene regulatory or chemical reaction networks, where no exact solution exists. For linear systems, the RWA becomes the exact result obtained from the maximum entropy principle. The RWA allows having a simple analytical, even though approximated, form of the solution, which is global and easier to deal with than the standard System Size Expansion (SSE). Here, we give some theoretically sufficient conditions for the validity of the RWA and estimate the order of error calculated by the approximation with respect to the number of particles. We compare RWA with SSE for two examples, a toy model and the more realistic dual phosphorylation cycle, governed by the same underlying process. Both approximations are compared with the exact integration of the master equation, showing for the RWA good performances of the same order or better than the SSE, even in regions where sufficient conditions are not met.
Collapse
Affiliation(s)
- Stefano Polizzi
- Department of Physics and Astronomy A. Righi, University of Bologna, 40127 Bologna, Italy
| | - Tommaso Marzi
- Department of Physics and Astronomy A. Righi, University of Bologna, 40127 Bologna, Italy
| | - Tommaso Matteuzzi
- Department of Physics and Astronomy, University of Florence, 50019 Sesto Fiorentino, Italy
| | - Gastone Castellani
- Department of Experimental, Diagnostic and Specialty Medicine, University of Bologna, 40138 Bologna, Italy
| | - Armando Bazzani
- Department of Physics and Astronomy A. Righi, University of Bologna, 40127 Bologna, Italy
| |
Collapse
|
36
|
Dal Co A, Ackermann M, van Vliet S. Spatial self-organization of metabolism in microbial systems: A matter of enzymes and chemicals. Cell Syst 2023; 14:98-108. [PMID: 36796335 DOI: 10.1016/j.cels.2022.12.009] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2022] [Revised: 08/14/2022] [Accepted: 12/21/2022] [Indexed: 02/17/2023]
Abstract
Most bacteria live in dense, spatially structured communities such as biofilms. The high density allows cells to alter the local microenvironment, whereas the limited mobility can cause species to become spatially organized. Together, these factors can spatially organize metabolic processes within microbial communities so that cells in different locations perform different metabolic reactions. The overall metabolic activity of a community depends both on how metabolic reactions are arranged in space and on how they are coupled, i.e., how cells in different regions exchange metabolites. Here, we review mechanisms that lead to the spatial organization of metabolic processes in microbial systems. We discuss factors that determine the length scales over which metabolic activities are arranged in space and highlight how the spatial organization of metabolic processes affects the ecology and evolution of microbial communities. Finally, we define key open questions that we believe should be the main focus of future research.
Collapse
Affiliation(s)
- Alma Dal Co
- Department of Computational Biology, University of Lausanne, 1015 Lausanne, Switzerland
| | - Martin Ackermann
- Department of Environmental Systems Science, ETH Zurich, 8092 Zurich, Switzerland; Department of Environmental Microbiology, Eawag, 8600 Duebendorf, Switzerland.
| | | |
Collapse
|
37
|
Samanta T, Kar S. Efficient Quantification of Extrinsic Fluctuations via Stochastic Simulations. Methods Mol Biol 2023; 2634:153-165. [PMID: 37074578 DOI: 10.1007/978-1-0716-3008-2_7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/20/2023]
Abstract
At the molecular level, all the biological processes are exposed to fluctuations emanating from various sources in and around the cellular system. Often these fluctuations dictate the outcome of a cell-fate decision-making event. Thus, having an accurate estimate of these fluctuations for any biological network is extremely important. There are well-established theoretical and numerical methods to quantify the intrinsic fluctuation present within a biological network arising due to the low copy numbers of cellular components. Unfortunately, the extrinsic fluctuations arising due to cell division events, epigenetic regulation, etc. have received very little attention. However, recent studies demonstrate that these extrinsic fluctuations significantly affect the transcriptional heterogeneity of certain important genes. Herein, we propose a new stochastic simulation algorithm to efficiently estimate these extrinsic fluctuations for experimentally constructed bidirectional transcriptional reporter systems along with the intrinsic variability. We use the Nanog transcriptional regulatory network and its variants to illustrate our numerical method. Our method reconciled experimental observations related to Nanog transcription, made exciting predictions, and can be applied to quantify intrinsic and extrinsic fluctuations for any similar transcriptional regulatory network.
Collapse
Affiliation(s)
| | - Sandip Kar
- Department of Chemistry, IIT Bombay, Mumbai, India.
| |
Collapse
|
38
|
Boe RH, Ayyappan V, Schuh L, Raj A. Allelic correlation is a marker of trade-offs between barriers to transmission of expression variability and signal responsiveness in genetic networks. Cell Syst 2022; 13:1016-1032.e6. [PMID: 36450286 PMCID: PMC9811561 DOI: 10.1016/j.cels.2022.10.008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2021] [Revised: 06/28/2022] [Accepted: 10/28/2022] [Indexed: 12/03/2022]
Abstract
Genetic networks should respond to signals but prevent the transmission of spontaneous fluctuations. Limited data from mammalian cells suggest that noise transmission is uncommon, but systematic claims about noise transmission have been limited by the inability to directly measure it. Here, we build a mathematical framework modeling allelic correlation and noise transmission, showing that allelic correlation and noise transmission correspond across model parameters and network architectures. Limiting noise transmission comes with the trade-off of being unresponsive to signals, and within responsive regimes, there is a further trade-off between response time and basal noise transmission. Analysis of allele-specific single-cell RNA-sequencing data revealed that genes encoding upstream factors in signaling pathways and cell-type-specific factors have higher allelic correlation than downstream factors, suggesting they are more subject to regulation. Overall, our findings suggest that some noise transmission must result from signal responsiveness, but it can be minimized by trading off for a slower response. A record of this paper's transparent peer review process is included in the supplemental information.
Collapse
Affiliation(s)
- Ryan H Boe
- Genetics and Epigenetics, Cell and Molecular Biology Graduate Group, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Vinay Ayyappan
- Department of Bioengineering, School of Engineering and Applied Sciences, University of Pennsylvania, Philadelphia, PA, USA
| | - Lea Schuh
- Institute of AI for Health, Helmholtz Zentrum München, German Research Center for Environmental Health, 85764 Neuherberg, Germany; Department of Mathematics, Technical University of Munich, Garching 85748, Germany
| | - Arjun Raj
- Department of Bioengineering, School of Engineering and Applied Sciences, University of Pennsylvania, Philadelphia, PA, USA; Department of Genetics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA.
| |
Collapse
|
39
|
Sharmeen N, Law C, Wu C. Polarization and cell-fate decision facilitated by the adaptor Ste50p in Saccharomyces cerevisiae. PLoS One 2022; 17:e0278614. [PMID: 36538537 PMCID: PMC9767377 DOI: 10.1371/journal.pone.0278614] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2022] [Accepted: 11/18/2022] [Indexed: 12/24/2022] Open
Abstract
In response to pheromone, many proteins localize on the plasma membrane of yeast cell to reform it into a polarized shmoo structure. The adaptor protein Ste50p, known as a pheromone signal enhancer critical for shmoo polarization, has never been explored systematically for its localization and function in the polarization process. Time-lapse single-cell imaging and quantitation shown here characterizes Ste50p involvement in the establishment of cell polarity. We found that Ste50p patches on the cell cortex mark the point of shmoo initiation, these patches could move, and remain associated with the growing shmoo tip in a pheromone concentration time-dependent manner until shmoo maturation. A Ste50p mutant impaired in patch localization suffers a delay in polarization. By quantitative analysis we show that polarization correlates with the rising levels of Ste50p, enabling rapid cell responses to pheromone that correspond to a critical level of Ste50p at the initial G1 phase. We exploited the quantitative differences in the pattern of Ste50p expression to correlate with the cell-cell phenotypic heterogeneity, showing Ste50p involvement in the cellular differentiation choice. Taken together, these findings present Ste50p to be part of the early shmoo development phase, suggesting that Ste50p may be involved with the polarisome in the initiation of polarization, and plays a role in regulating the polarized growth of shmoo during pheromone response.
Collapse
Affiliation(s)
- Nusrat Sharmeen
- Division of Experimental Medicine, Department of Medicine, McGill University, Montreal, Quebec, Canada
- * E-mail:
| | - Chris Law
- Centre for Microscopy and Cellular Imaging, Department of Biology, Concordia University, Montreal, Quebec, Canada
| | - Cunle Wu
- Division of Experimental Medicine, Department of Medicine, McGill University, Montreal, Quebec, Canada
- Human Health Therapeutics Research Centre, National Research Council Canada, Montreal, Quebec, Canada
| |
Collapse
|
40
|
Abstract
The ability of bacteria to respond to changes in their environment is critical to their survival, allowing them to withstand stress, form complex communities, and induce virulence responses during host infection. A remarkable feature of many of these bacterial responses is that they are often variable across individual cells, despite occurring in an isogenic population exposed to a homogeneous environmental change, a phenomenon known as phenotypic heterogeneity. Phenotypic heterogeneity can enable bet-hedging or division of labor strategies that allow bacteria to survive fluctuating conditions. Investigating the significance of phenotypic heterogeneity in environmental transitions requires dynamic, single-cell data. Technical advances in quantitative single-cell measurements, imaging, and microfluidics have led to a surge of publications on this topic. Here, we review recent discoveries on single-cell bacterial responses to environmental transitions of various origins and complexities, from simple diauxic shifts to community behaviors in biofilm formation to virulence regulation during infection. We describe how these studies firmly establish that this form of heterogeneity is prevalent and a conserved mechanism by which bacteria cope with fluctuating conditions. We end with an outline of current challenges and future directions for the field. While it remains challenging to predict how an individual bacterium will respond to a given environmental input, we anticipate that capturing the dynamics of the process will begin to resolve this and facilitate rational perturbation of environmental responses for therapeutic and bioengineering purposes.
Collapse
|
41
|
Kumar MS, Slud EV, Hehnly C, Zhang L, Broach J, Irizarry RA, Schiff SJ, Paulson JN. Differential richness inference for 16S rRNA marker gene surveys. Genome Biol 2022; 23:166. [PMID: 35915508 PMCID: PMC9344657 DOI: 10.1186/s13059-022-02722-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Accepted: 06/28/2022] [Indexed: 12/24/2022] Open
Abstract
Background Individual and environmental health outcomes are frequently linked to changes in the diversity of associated microbial communities. Thus, deriving health indicators based on microbiome diversity measures is essential. While microbiome data generated using high-throughput 16S rRNA marker gene surveys are appealing for this purpose, 16S surveys also generate a plethora of spurious microbial taxa. Results When this artificial inflation in the observed number of taxa is ignored, we find that changes in the abundance of detected taxa confound current methods for inferring differences in richness. Experimental evidence, theory-guided exploratory data analyses, and existing literature support the conclusion that most sub-genus discoveries are spurious artifacts of clustering 16S sequencing reads. We proceed to model a 16S survey’s systematic patterns of sub-genus taxa generation as a function of genus abundance to derive a robust control for false taxa accumulation. These controls unlock classical regression approaches for highly flexible differential richness inference at various levels of the surveyed microbial assemblage: from sample groups to specific taxa collections. The proposed methodology for differential richness inference is available through an R package, Prokounter. Conclusions False species discoveries bias richness estimation and confound differential richness inference. In the case of 16S microbiome surveys, supporting evidence indicate that most sub-genus taxa are spurious. Based on this finding, a flexible method is proposed and is shown to overcome the confounding problem noted with current approaches for differential richness inference. Package availability: https://github.com/mskb01/prokounter Supplementary Information The online version contains supplementary material available at 10.1186/s13059-022-02722-x.
Collapse
|
42
|
Koraichi MB, Touzel MP, Mazzolini A, Mora T, Walczak AM. NoisET: Noise Learning and Expansion Detection of T-Cell Receptors. J Phys Chem A 2022; 126:7407-7414. [PMID: 36178325 DOI: 10.1021/acs.jpca.2c05002] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
High-throughput sequencing of T- and B-cell receptors makes it possible to track immune repertoires across time, in different tissues, in acute and chronic diseases and in healthy individuals. However, quantitative comparison between repertoires is confounded by variability in the read count of each receptor clonotype due to sampling, library preparation, and expression noise. We review methods for accounting for both biological and experimental noise and present an easy-to-use python package NoisET that implements and generalizes a previously developed Bayesian method. It can be used to learn experimental noise models for repertoire sequencing from replicates, and to detect responding clones following a stimulus. We test the package on different repertoire sequencing technologies and data sets. We review how such approaches have been used to identify responding clonotypes in vaccination and disease data. Availability: NoisET is freely available to use with source code at github.com/statbiophys/NoisET.
Collapse
Affiliation(s)
- Meriem Bensouda Koraichi
- Laboratoire de physique de l' École normale supérieure, CNRS, PSL University, Sorbonne Université, and Université de Paris, Paris75005, France
| | | | - Andrea Mazzolini
- Laboratoire de physique de l' École normale supérieure, CNRS, PSL University, Sorbonne Université, and Université de Paris, Paris75005, France
| | - Thierry Mora
- Laboratoire de physique de l' École normale supérieure, CNRS, PSL University, Sorbonne Université, and Université de Paris, Paris75005, France
| | - Aleksandra M Walczak
- Laboratoire de physique de l' École normale supérieure, CNRS, PSL University, Sorbonne Université, and Université de Paris, Paris75005, France
| |
Collapse
|
43
|
Martin CS, Jubelin G, Darsonval M, Leroy S, Leneveu-Jenvrin C, Hmidene G, Omhover L, Stahl V, Guillier L, Briandet R, Desvaux M, Dubois-Brissonnet F. Genetic, physiological, and cellular heterogeneities of bacterial pathogens in food matrices: Consequences for food safety. Compr Rev Food Sci Food Saf 2022; 21:4294-4326. [PMID: 36018457 DOI: 10.1111/1541-4337.13020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2022] [Revised: 07/12/2022] [Accepted: 07/14/2022] [Indexed: 01/28/2023]
Abstract
In complex food systems, bacteria live in heterogeneous microstructures, and the population displays phenotypic heterogeneities at the single-cell level. This review provides an overview of spatiotemporal drivers of phenotypic heterogeneity of bacterial pathogens in food matrices at three levels. The first level is the genotypic heterogeneity due to the possibility for various strains of a given species to contaminate food, each of them having specific genetic features. Then, physiological heterogeneities are induced within the same strain, due to specific microenvironments and heterogeneous adaptative responses to the food microstructure. The third level of phenotypic heterogeneity is related to cellular heterogeneity of the same strain in a specific microenvironment. Finally, we consider how these phenotypic heterogeneities at the single-cell level could be implemented in mathematical models to predict bacterial behavior and help ensure microbiological food safety.
Collapse
Affiliation(s)
- Cédric Saint Martin
- MICALIS Institute, Université Paris-Saclay, INRAE, AgroParisTech, Jouy-en-Josas, France.,Université Clermont Auvergne, INRAE, UMR454 MEDIS, Clermont-Ferrand, France
| | - Grégory Jubelin
- Université Clermont Auvergne, INRAE, UMR454 MEDIS, Clermont-Ferrand, France
| | - Maud Darsonval
- MICALIS Institute, Université Paris-Saclay, INRAE, AgroParisTech, Jouy-en-Josas, France
| | - Sabine Leroy
- Université Clermont Auvergne, INRAE, UMR454 MEDIS, Clermont-Ferrand, France
| | - Charlène Leneveu-Jenvrin
- MICALIS Institute, Université Paris-Saclay, INRAE, AgroParisTech, Jouy-en-Josas, France.,Association pour le Développement de l'Industrie de la Viande (ADIV), Clermont-Ferrand, France
| | - Ghaya Hmidene
- Risk Assessment Department, ANSES, Maisons-Alfort, France
| | - Lysiane Omhover
- Aerial, Technical Institute of Agro-Industry, Illkirch, France
| | - Valérie Stahl
- Aerial, Technical Institute of Agro-Industry, Illkirch, France
| | | | - Romain Briandet
- MICALIS Institute, Université Paris-Saclay, INRAE, AgroParisTech, Jouy-en-Josas, France
| | - Mickaël Desvaux
- Université Clermont Auvergne, INRAE, UMR454 MEDIS, Clermont-Ferrand, France
| | | |
Collapse
|
44
|
Gupta A, Martin-Rufino JD, Jones TR, Subramanian V, Qiu X, Grody EI, Bloemendal A, Weng C, Niu SY, Min KH, Mehta A, Zhang K, Siraj L, Al' Khafaji A, Sankaran VG, Raychaudhuri S, Cleary B, Grossman S, Lander ES. Inferring gene regulation from stochastic transcriptional variation across single cells at steady state. Proc Natl Acad Sci U S A 2022; 119:e2207392119. [PMID: 35969771 PMCID: PMC9407670 DOI: 10.1073/pnas.2207392119] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2022] [Accepted: 07/20/2022] [Indexed: 12/24/2022] Open
Abstract
Regulatory relationships between transcription factors (TFs) and their target genes lie at the heart of cellular identity and function; however, uncovering these relationships is often labor-intensive and requires perturbations. Here, we propose a principled framework to systematically infer gene regulation for all TFs simultaneously in cells at steady state by leveraging the intrinsic variation in the transcriptional abundance across single cells. Through modeling and simulations, we characterize how transcriptional bursts of a TF gene are propagated to its target genes, including the expected ranges of time delay and magnitude of maximum covariation. We distinguish these temporal trends from the time-invariant covariation arising from cell states, and we delineate the experimental and technical requirements for leveraging these small but meaningful cofluctuations in the presence of measurement noise. While current technology does not yet allow adequate power for definitively detecting regulatory relationships for all TFs simultaneously in cells at steady state, we investigate a small-scale dataset to inform future experimental design. This study supports the potential value of mapping regulatory connections through stochastic variation, and it motivates further technological development to achieve its full potential.
Collapse
Affiliation(s)
- Anika Gupta
- Broad Institute of MIT and Harvard, Cambridge, MA 02142
- Department of Biomedical Informatics, Harvard Medical School, Boston, MA 02115
| | - Jorge D. Martin-Rufino
- Broad Institute of MIT and Harvard, Cambridge, MA 02142
- Division of Hematology/Oncology, Boston Children’s Hospital, Boston, MA 02115
- Dana-Farber Cancer Institute, Boston, MA 02215
| | | | | | - Xiaojie Qiu
- Whitehead Institute for Biomedical Research, Cambridge, MA 02142
- HHMI, Massachusetts Institute of Technology, Cambridge, MA 02139
| | | | | | - Chen Weng
- Broad Institute of MIT and Harvard, Cambridge, MA 02142
- Division of Hematology/Oncology, Boston Children’s Hospital, Boston, MA 02115
- Dana-Farber Cancer Institute, Boston, MA 02215
- Whitehead Institute for Biomedical Research, Cambridge, MA 02142
| | | | - Kyung Hoi Min
- Whitehead Institute for Biomedical Research, Cambridge, MA 02142
- Department of Electrical Engineering and Computer Science, Massachusetts Institute of Technology, Cambridge, MA 02139
| | - Arnav Mehta
- Broad Institute of MIT and Harvard, Cambridge, MA 02142
- Dana-Farber Cancer Institute, Boston, MA 02215
- Department of Medicine, Massachusetts General Hospital, Boston, MA 02114
| | - Kaite Zhang
- Broad Institute of MIT and Harvard, Cambridge, MA 02142
| | - Layla Siraj
- Broad Institute of MIT and Harvard, Cambridge, MA 02142
| | | | - Vijay G. Sankaran
- Broad Institute of MIT and Harvard, Cambridge, MA 02142
- Division of Hematology/Oncology, Boston Children’s Hospital, Boston, MA 02115
- Dana-Farber Cancer Institute, Boston, MA 02215
| | - Soumya Raychaudhuri
- Broad Institute of MIT and Harvard, Cambridge, MA 02142
- Department of Biomedical Informatics, Harvard Medical School, Boston, MA 02115
- Center for Data Sciences, Brigham and Women’s Hospital, Boston, MA 02115
| | - Brian Cleary
- Broad Institute of MIT and Harvard, Cambridge, MA 02142
| | | | - Eric S. Lander
- Broad Institute of MIT and Harvard, Cambridge, MA 02142
- Department of Biology, Massachusetts Institute of Technology, Cambridge, MA 02142
- Department of Systems Biology, Harvard Medical School, Boston, MA 02115
| |
Collapse
|
45
|
Loell K, Wu Y, Staller MV, Cohen B. Activation domains can decouple the mean and noise of gene expression. Cell Rep 2022; 40:111118. [PMID: 35858548 PMCID: PMC9912357 DOI: 10.1016/j.celrep.2022.111118] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2021] [Revised: 01/18/2022] [Accepted: 06/28/2022] [Indexed: 11/03/2022] Open
Abstract
Regulatory mechanisms set a gene's average level of expression, but a gene's expression constantly fluctuates around that average. These stochastic fluctuations, or expression noise, play a role in cell-fate transitions, bet hedging in microbes, and the development of chemotherapeutic resistance in cancer. An outstanding question is what regulatory mechanisms contribute to noise. Here, we demonstrate that, for a fixed mean level of expression, strong activation domains (ADs) at low abundance produce high expression noise, while weak ADs at high abundance generate lower expression noise. We conclude that differences in noise can be explained by the interplay between a TF's nuclear concentration and the strength of its AD's effect on mean expression, without invoking differences between classes of ADs. These results raise the possibility of engineering gene expression noise independently of mean levels in synthetic biology contexts and provide a potential mechanism for natural selection to tune the noisiness of gene expression.
Collapse
Affiliation(s)
- Kaiser Loell
- Department of Genetics, Washington University School of Medicine in St. Louis, St. Louis, MO 63108, USA,The Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine in St. Louis, St. Louis, MO 63108, USA
| | - Yawei Wu
- Department of Genetics, Washington University School of Medicine in St. Louis, St. Louis, MO 63108, USA,The Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine in St. Louis, St. Louis, MO 63108, USA
| | - Max V. Staller
- Center for Computational Biology, University of California, Berkeley, Berkeley, CA 94720, USA
| | - Barak Cohen
- Department of Genetics, Washington University School of Medicine in St. Louis, St. Louis, MO 63108, USA; The Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine in St. Louis, St. Louis, MO 63108, USA.
| |
Collapse
|
46
|
Grimaudo R, Lazzari P, Solidoro C, Valenti D. Effects of solar irradiance noise on a complex marine trophic web. Sci Rep 2022; 12:12163. [PMID: 35842433 PMCID: PMC9288520 DOI: 10.1038/s41598-022-16236-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2022] [Accepted: 07/06/2022] [Indexed: 11/09/2022] Open
Abstract
The analysis of experimental data of the solar irradiance, collected on the marine surface, clearly highlights the intrinsic stochasticity of such an environmental parameter. Given this result, effects of randomly fluctuating irradiance on the population dynamics of a marine ecosystem are studied on the basis of the stochastic 0-dimensional biogeochemical flux model. The noisy fluctuations of the irradiance are formally described as a multiplicative Ornstein-Uhlenbeck process, that is a self-correlated Gaussian noise. Nonmonotonic behaviours of the variance of the marine populations' biomass are found with respect to the intensity and the autocorrelation time of the noise source, manifesting a noise-induced transition of the ecosystem to an out-of-equilibrium steady state. Moreover, evidence of noise-induced effects on the organic carbon cycling processes underlying the food web dynamics are highlighted. The reported results clearly show the profound impact the stochastic environmental variables can have on both the populations and the biogeochemistry at the basis of a marine trophic network.
Collapse
Affiliation(s)
- Roberto Grimaudo
- Dipartimento di Fisica e Chimica Emilio Segré, Universitá degli Studi di Palermo, Viale delle Scienze, Ed. 18, I-90128, Palermo, Italy
| | - Paolo Lazzari
- National Institute of Oceanography and Applied Geophysics - OGS, via Beirut 2, I-34014, Trieste, Italy.
| | - Cosimo Solidoro
- National Institute of Oceanography and Applied Geophysics - OGS, via Beirut 2, I-34014, Trieste, Italy
| | - Davide Valenti
- Dipartimento di Fisica e Chimica Emilio Segré, Universitá degli Studi di Palermo, Viale delle Scienze, Ed. 18, I-90128, Palermo, Italy
| |
Collapse
|
47
|
Das AK. Stochastic gene transcription with non-competitive transcription regulatory architecture. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2022; 45:61. [PMID: 35831727 DOI: 10.1140/epje/s10189-022-00213-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/08/2022] [Accepted: 06/22/2022] [Indexed: 06/15/2023]
Abstract
The transcription factors, such as activators and repressors, can interact with the promoter of gene either in a competitive or non-competitive way. In this paper, we construct a stochastic model with non-competitive transcriptional regulatory architecture and develop an analytical theory that re-establishes the experimental results with an improved data fitting. The analytical expressions in the theory allow us to study the nature of the system corresponding to any of its parameters and hence, enable us to find out the factors that govern the regulation of gene expression for that architecture. We notice that, along with transcriptional reinitiation and repressors, there are other parameters that can control the noisiness of this network. We also observe that, the Fano factor (at mRNA level) varies from sub-Poissonian regime to super-Poissonian regime. In addition to the aforementioned properties, we observe some anomalous characteristics of the Fano factor (at mRNA level) and that of the variance of protein at lower activator concentrations in the presence of repressor molecules. This model is useful to understand the architecture of interactions which may buffer the stochasticity inherent to gene transcription.
Collapse
|
48
|
Zreika S, Fourneaux C, Vallin E, Modolo L, Seraphin R, Moussy A, Ventre E, Bouvier M, Ozier-Lafontaine A, Bonnaffoux A, Picard F, Gandrillon O, Gonin-Giraud S. Evidence for close molecular proximity between reverting and undifferentiated cells. BMC Biol 2022; 20:155. [PMID: 35794592 PMCID: PMC9258043 DOI: 10.1186/s12915-022-01363-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2022] [Accepted: 06/27/2022] [Indexed: 11/28/2022] Open
Abstract
Background According to Waddington’s epigenetic landscape concept, the differentiation process can be illustrated by a cell akin to a ball rolling down from the top of a hill (proliferation state) and crossing furrows before stopping in basins or “attractor states” to reach its stable differentiated state. However, it is now clear that some committed cells can retain a certain degree of plasticity and reacquire phenotypical characteristics of a more pluripotent cell state. In line with this dynamic model, we have previously shown that differentiating cells (chicken erythrocytic progenitors (T2EC)) retain for 24 h the ability to self-renew when transferred back in self-renewal conditions. Despite those intriguing and promising results, the underlying molecular state of those “reverting” cells remains unexplored. The aim of the present study was therefore to molecularly characterize the T2EC reversion process by combining advanced statistical tools to make the most of single-cell transcriptomic data. For this purpose, T2EC, initially maintained in a self-renewal medium (0H), were induced to differentiate for 24H (24H differentiating cells); then, a part of these cells was transferred back to the self-renewal medium (48H reverting cells) and the other part was maintained in the differentiation medium for another 24H (48H differentiating cells). For each time point, cell transcriptomes were generated using scRT-qPCR and scRNAseq. Results Our results showed a strong overlap between 0H and 48H reverting cells when applying dimensional reduction. Moreover, the statistical comparison of cell distributions and differential expression analysis indicated no significant differences between these two cell groups. Interestingly, gene pattern distributions highlighted that, while 48H reverting cells have gene expression pattern more similar to 0H cells, they are not completely identical, which suggest that for some genes a longer delay may be required for the cells to fully recover. Finally, sparse PLS (sparse partial least square) analysis showed that only the expression of 3 genes discriminates 48H reverting and 0H cells. Conclusions Altogether, we show that reverting cells return to an earlier molecular state almost identical to undifferentiated cells and demonstrate a previously undocumented physiological and molecular plasticity during the differentiation process, which most likely results from the dynamic behavior of the underlying molecular network. Supplementary Information The online version contains supplementary material available at 10.1186/s12915-022-01363-7.
Collapse
|
49
|
Zinani OQ, Keseroğlu K, Dey S, Ay A, Singh A, Özbudak EM. Gene copy number and negative feedback differentially regulate transcriptional variability of segmentation clock genes. iScience 2022; 25:104579. [PMID: 35789861 PMCID: PMC9250017 DOI: 10.1016/j.isci.2022.104579] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2021] [Revised: 03/11/2022] [Accepted: 06/07/2022] [Indexed: 10/26/2022] Open
|
50
|
Kolbe N, Hexemer L, Bammert LM, Loewer A, Lukáčová-Medvid’ová M, Legewie S. Data-based stochastic modeling reveals sources of activity bursts in single-cell TGF-β signaling. PLoS Comput Biol 2022; 18:e1010266. [PMID: 35759468 PMCID: PMC9269928 DOI: 10.1371/journal.pcbi.1010266] [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: 07/30/2021] [Revised: 07/08/2022] [Accepted: 05/30/2022] [Indexed: 11/30/2022] Open
Abstract
Cells sense their surrounding by employing intracellular signaling pathways that transmit hormonal signals from the cell membrane to the nucleus. TGF-β/SMAD signaling encodes various cell fates, controls tissue homeostasis and is deregulated in diseases such as cancer. The pathway shows strong heterogeneity at the single-cell level, but quantitative insights into mechanisms underlying fluctuations at various time scales are still missing, partly due to inefficiency in the calibration of stochastic models that mechanistically describe signaling processes. In this work we analyze single-cell TGF-β/SMAD signaling and show that it exhibits temporal stochastic bursts which are dose-dependent and whose number and magnitude correlate with cell migration. We propose a stochastic modeling approach to mechanistically describe these pathway fluctuations with high computational efficiency. Employing high-order numerical integration and fitting to burst statistics we enable efficient quantitative parameter estimation and discriminate models that assume noise in different reactions at the receptor level. This modeling approach suggests that stochasticity in the internalization of TGF-β receptors into endosomes plays a key role in the observed temporal bursting. Further, the model predicts the single-cell dynamics of TGF-β/SMAD signaling in untested conditions, e.g., successfully reflects memory effects of signaling noise and cellular sensitivity towards repeated stimulation. Taken together, our computational framework based on burst analysis, noise modeling and path computation scheme is a suitable tool for the data-based modeling of complex signaling pathways, capable of identifying the source of temporal noise. Fluctuations in molecular networks give rise to heterogeneity in cellular behavior and therefore promote the diversification of tissues. For a better understanding of cellular decision making, it is important to identify sources of molecular fluctuations and to quantitatively describe them by predictive mathematical models. In this work, we focused on temporal fluctuations of the TGF-β signaling pathway that is important for controlling cell division and migration. We characterized a single-cell dataset comprising hundreds of cells using time series analysis and a large-scale stochastic model. By fitting several model variants to the data, we identified the stochastic internalization of cell surface receptors into endosomes as a main source of temporal fluctuations (’bursts’) in the signaling pathway. The corresponding model accurately predicted novel experimental data, and provided insights into the long-term memory of signaling fluctuations. In summary, we propose a modeling approach to quantitatively describe heterogeneous behavior in large-scale single-cell datasets and to identify the underlying biological mechanisms.
Collapse
Affiliation(s)
- Niklas Kolbe
- Institute of Geometry and Practical Mathematics, RWTH Aachen University, Aachen, Germany
- * E-mail: (NK); (SL)
| | - Lorenz Hexemer
- Institute of Molecular Biology (IMB), Mainz, Germany
- Department of Systems Biology, Institute for Biomedical Genetics (IBMG), University of Stuttgart, Stuttgart, Germany
| | | | - Alexander Loewer
- Systems Biology of the Stress Response, Department of Biology, Technical University of Darmstadt, Darmstadt, Germany
| | | | - Stefan Legewie
- Department of Systems Biology, Institute for Biomedical Genetics (IBMG), University of Stuttgart, Stuttgart, Germany
- Stuttgart Research Center for Systems Biology (SRCSB), University of Stuttgart, Stuttgart, Germany
- * E-mail: (NK); (SL)
| |
Collapse
|