scholarly journals Predicting synchronized gene coexpression patterns from fibration symmetries in gene regulatory networks in bacteria

2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Ian Leifer ◽  
Mishael Sánchez-Pérez ◽  
Cecilia Ishida ◽  
Hernán A. Makse

Abstract Background Gene regulatory networks coordinate the expression of genes across physiological states and ensure a synchronized expression of genes in cellular subsystems, critical for the coherent functioning of cells. Here we address the question whether it is possible to predict gene synchronization from network structure alone. We have recently shown that synchronized gene expression can be predicted from symmetries in the gene regulatory networks described by the concept of symmetry fibrations. We showed that symmetry fibrations partition the genes into groups called fibers based on the symmetries of their ’input trees’, the set of paths in the network through which signals can reach a gene. In idealized dynamic gene expression models, all genes in a fiber are perfectly synchronized, while less idealized models—with gene input functions differencing between genes—predict symmetry breaking and desynchronization. Results To study the functional role of gene fibers and to test whether some of the fiber-induced coexpression remains in reality, we analyze gene fibrations for the gene regulatory networks of E. coli and B. subtilis and confront them with expression data. We find approximate gene coexpression patterns consistent with symmetry fibrations with idealized gene expression dynamics. This shows that network structure alone provides useful information about gene synchronization, and suggest that gene input functions within fibers may be further streamlined by evolutionary pressures to realize a coexpression of genes. Conclusions Thus, gene fibrations provide a sound conceptual tool to describe tunable coexpression induced by network topology and shaped by mechanistic details of gene expression.

2020 ◽  
Author(s):  
Ian Leifer ◽  
Mishael Sanchez ◽  
Cecilia Ishida ◽  
Hernan Makse

Abstract Background: Gene regulatory networks coordinate the expression of genes across physiological states and ensure a synchronized expression of genes in cellular subsystems, critical for the coherent functioning of cells. Here we address the questions whether it is possible to predict gene synchronization from network structure alone. We have recently shown that synchronized gene expression may be predicted from symmetries in the transcriptional regulatory networks (TRN) and described by the concept of symmetry fibrations. We showed that symmetry fibrations partition the genes into groups called fibers based on the symmetries of their 'input trees', the set of paths in the network through which signals can reach a gene. In idealized dynamic gene expression models, all genes in a fiber are perfectly synchronized, while less idealized models -- with gene input functions differencing between genes -- predict symmetry breaking and desynchronization. Results: To study the functional role of gene fibers and to test whether some of the fiber-induced coexpression remains in reality, we analyze gene fibrations for the transcription networks of E. coli and B. subtilis and confront them with expression data. We find approximate gene coexpression patterns consistent with symmetry fibrations with idealized gene expression dynamics. This shows that network structure alone provides useful information about gene synchronization, and suggest that gene input functions within fibers may be further streamlined by evolutionary pressures to realize a coexpression of genes. Conclusions: Thus, gene fibrations provides a sound conceptual tool to describe tunable coexpression induced by network topology and shaped by mechanistic details of gene expression.


2016 ◽  
Author(s):  
Dianbo Liu ◽  
Luca Albergante ◽  
Timothy J Newman

AbstractUsing a combination of mathematical modelling, statistical simulation and large-scale data analysis we study the properties of linear regulatory chains (LRCs) within gene regulatory networks (GRNs). Our modelling indicates that downstream genes embedded within LRCs are highly insulated from the variation in expression of upstream genes, and thus LRCs act as attenuators. This observation implies a progressively weaker functionality of LRCs as their length increases. When analysing the preponderance of LRCs in the GRNs of E. coli K12 and several other organisms, we find that very long LRCs are essentially absent. In both E. coli and M. tuberculosis we find that four-gene LRCs are intimately linked to identical feedback loops that are involved in potentially chaotic stress response, indicating that the dynamics of these potentially destabilising motifs are strongly restrained under homeostatic conditions. The same relationship is observed in a human cancer cell line (K562), and we postulate that four-gene LRCs act as “universal attenuators”. These findings suggest a role for long LRCs in dampening variation in gene expression, thereby protecting cell identity, and in controlling dramatic shifts in cell-wide gene expression through inhibiting chaos-generating motifs.In briefWe present a general principle that linear regulatory chains exponentially attenuate the range of expression in gene regulatory networks. The discovery of a universal interplay between linear regulatory chains and genetic feedback loops in microorganisms and a human cancer cell line is analysed and discussed.HighlightsWithin gene networks, linear regulatory chains act as exponentially strong attenuators of upstream variationBecause of their exponential behaviour, linear regulatory chains beyond a few genes provide no additional functionality and are rarely observed in gene networks across a range of different organismsNovel interactions between four-gene linear regulatory chains and feedback loops were discovered in E. coli, M. tuberculosis and human cancer cells, suggesting a universal mechanism of control.


2019 ◽  
Author(s):  
Zhang Zhang ◽  
Lifei Wang ◽  
Shuo Wang ◽  
Ruyi Tao ◽  
Jingshu Xiao ◽  
...  

SummaryReconstructing gene regulatory networks (GRNs) and inferring the gene dynamics are important to understand the behavior and the fate of the normal and abnormal cells. Gene regulatory networks could be reconstructed by experimental methods or from gene expression data. Recent advances in Single Cell RNA sequencing technology and the computational method to reconstruct trajectory have generated huge scRNA-seq data tagged with additional time labels. Here, we present a deep learning model “Neural Gene Network Constructor” (NGNC), for inferring gene regulatory network and reconstructing the gene dynamics simultaneously from time series gene expression data. NGNC is a model-free heterogenous model, which can reconstruct any network structure and non-linear dynamics. It consists of two parts: a network generator which incorporating gumbel softmax technique to generate candidate network structure, and a dynamics learner which adopting multiple feedforward neural networks to predict the dynamics. We compare our model with other well-known frameworks on the data set generated by GeneNetWeaver, and achieve the state of the arts results both on network reconstruction and dynamics learning.


2018 ◽  
Vol 47 (1) ◽  
pp. 447-467 ◽  
Author(s):  
David L. Shis ◽  
Matthew R. Bennett, ◽  
Oleg A. Igoshin

The ability of bacterial cells to adjust their gene expression program in response to environmental perturbation is often critical for their survival. Recent experimental advances allowing us to quantitatively record gene expression dynamics in single cells and in populations coupled with mathematical modeling enable mechanistic understanding on how these responses are shaped by the underlying regulatory networks. Here, we review how the combination of local and global factors affect dynamical responses of gene regulatory networks. Our goal is to discuss the general principles that allow extrapolation from a few model bacteria to less understood microbes. We emphasize that, in addition to well-studied effects of network architecture, network dynamics are shaped by global pleiotropic effects and cell physiology.


2021 ◽  
Author(s):  
Kenji Okubo ◽  
Kunihiko Kaneko

AbstractMendelian inheritance is a fundamental law of genetics. Considering two alleles in a diploid, a phenotype of a heterotype is dominated by a particular homotype according to the law of dominance. This picture is usually based on simple genotype-phenotype mapping in which one gene regulates one phenotype. However, in reality, some interactions between genes can result in deviation from Mendelian dominance.Here, by using the numerical evolution of diploid gene regulatory networks (GRNs), we discuss whether Mendelian dominance evolves beyond the classical case of one-to-one genotype-phenotype mapping. We examine whether complex genotype-phenotype mapping can achieve Mendelian dominance through the evolution of the GRN with interacting genes. Specifically, we extend the GRN model to a diploid case, in which two GRN matrices are added to give gene expression dynamics, and simulate evolution with meiosis and recombination. Our results reveal that Mendelian dominance evolves even under complex genotype-phenotype mapping. This dominance is achieved via a group of genotypes that differ from each other but have a common phenotype given by the expression of target genes. Calculating the degree of dominance shows that it increases through the evolution, correlating closely with the decrease in phenotypic fluctuations and the increase in robustness to initial noise. This evolution of Mendelian dominance is associated with phenotypic robustness against meiosis-induced genome mixing, whereas sexual recombination arising from the mixing of chromosomes from the parents further enhances dominance and robustness. Owing to this dominance, the robustness to genetic differences increases, while the optimal fitness is sustained up to a large difference between the two genomes. In summary, Mendelian dominance is achieved by groups of genotypes that are associated with the increase in phenotypic robustness to noise.Author summaryMendelian dominance is one of the most fundamental laws in genetics. When two conflicting characters occur in a single diploid, the dominant character is always chosen. Assuming that one gene makes one character, this law is simple to grasp. However, in reality, phenotypes are generated via interactions between several genes, which may alter Mendel’s dominance law. The evolution of robustness to noise and mutations has been investigated extensively using complex expression dynamics with gene regulatory networks. Here, we applied gene-expression dynamics with complex interactions to the case of a diploid and simulated the evolution of the gene regulatory network to generate the optimal phenotype given by a certain gene expression pattern. Interestingly, after evolution, Mendelian dominance is achieved via a group of genes. This group-based Mendelian dominance is shaped by phenotype insensitivity to genome mixing by meiosis and evolves concurrently with the robustness to noise. By focusing on the influence of phenotypic robustness, which has received considerable attention recently, our result provides a novel perspective as to why Mendel’s law of dominance is commonly observed.


2021 ◽  
Vol 7 (24) ◽  
pp. eabf8210
Author(s):  
Miki Tokuoka ◽  
Kazuki Maeda ◽  
Kenji Kobayashi ◽  
Atsushi Mochizuki ◽  
Yutaka Satou

In animal embryos, gene regulatory networks control the dynamics of gene expression in cells and coordinate such dynamics among cells. In ascidian embryos, gene expression dynamics have been dissected at the single-cell resolution. Here, we revealed mathematical functions that represent the regulatory logics of all regulatory genes expressed at the 32-cell stage when the germ layers are largely specified. These functions collectively explain the entire mechanism by which gene expression dynamics are controlled coordinately in early embryos. We found that regulatory functions for genes expressed in each of the specific lineages contain a common core regulatory mechanism. Last, we showed that the expression of the regulatory genes became reproducible by calculation and controllable by experimental manipulations. Thus, these regulatory functions represent an architectural design for the germ layer specification of this chordate and provide a platform for simulations and experiments to understand the operating principles of gene regulatory networks.


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Mika J. Välimäki ◽  
Robert S. Leigh ◽  
Sini M. Kinnunen ◽  
Alexander R. March ◽  
Ana Hernández de Sande ◽  
...  

AbstractBackgroundPharmacological modulation of cell fate decisions and developmental gene regulatory networks holds promise for the treatment of heart failure. Compounds that target tissue-specific transcription factors could overcome non-specific effects of small molecules and lead to the regeneration of heart muscle following myocardial infarction. Due to cellular heterogeneity in the heart, the activation of gene programs representing specific atrial and ventricular cardiomyocyte subtypes would be highly desirable. Chemical compounds that modulate atrial and ventricular cell fate could be used to improve subtype-specific differentiation of endogenous or exogenously delivered progenitor cells in order to promote cardiac regeneration.MethodsTranscription factor GATA4-targeted compounds that have previously shown in vivo efficacy in cardiac injury models were tested for stage-specific activation of atrial and ventricular reporter genes in differentiating pluripotent stem cells using a dual reporter assay. Chemically induced gene expression changes were characterized by qRT-PCR, global run-on sequencing (GRO-seq) and immunoblotting, and the network of cooperative proteins of GATA4 and NKX2-5 were further explored by the examination of the GATA4 and NKX2-5 interactome by BioID. Reporter gene assays were conducted to examine combinatorial effects of GATA-targeted compounds and bromodomain and extraterminal domain (BET) inhibition on chamber-specific gene expression.ResultsGATA4-targeted compounds 3i-1000 and 3i-1103 were identified as differential modulators of atrial and ventricular gene expression. More detailed structure-function analysis revealed a distinct subclass of GATA4/NKX2-5 inhibitory compounds with an acetyl lysine-like domain that contributed to ventricular cells (%Myl2-eGFP+). Additionally, BioID analysis indicated broad interaction between GATA4 and BET family of proteins, such as BRD4. This indicated the involvement of epigenetic modulators in the regulation of GATA-dependent transcription. In this line, reporter gene assays with combinatorial treatment of 3i-1000 and the BET bromodomain inhibitor (+)-JQ1 demonstrated the cooperative role of GATA4 and BRD4 in the modulation of chamber-specific cardiac gene expression.ConclusionsCollectively, these results indicate the potential for therapeutic alteration of cell fate decisions and pathological gene regulatory networks by GATA4-targeted compounds modulating chamber-specific transcriptional programs in multipotent cardiac progenitor cells and cardiomyocytes. The compound scaffolds described within this study could be used to develop regenerative strategies for myocardial regeneration.


2017 ◽  
Vol 15 (02) ◽  
pp. 1650045 ◽  
Author(s):  
Olga V. Petrovskaya ◽  
Evgeny D. Petrovskiy ◽  
Inna N. Lavrik ◽  
Vladimir A. Ivanisenko

Gene network modeling is one of the widely used approaches in systems biology. It allows for the study of complex genetic systems function, including so-called mosaic gene networks, which consist of functionally interacting subnetworks. We conducted a study of a mosaic gene networks modeling method based on integration of models of gene subnetworks by linear control functionals. An automatic modeling of 10,000 synthetic mosaic gene regulatory networks was carried out using computer experiments on gene knockdowns/knockouts. Structural analysis of graphs of generated mosaic gene regulatory networks has revealed that the most important factor for building accurate integrated mathematical models, among those analyzed in the study, is data on expression of genes corresponding to the vertices with high properties of centrality.


Sign in / Sign up

Export Citation Format

Share Document