Robustness Of Gene Regulatory Networks Research Articles

Mechanisms of robustness in gene regulatory networks involved in neural development.

The functions of living organisms are affected by different kinds of perturbation, both internal and external, which in many cases have functional effects and phenotypic impact. The effects of these perturbations become particularly relevant for multicellular organisms with complex body patterns and cell type heterogeneity, where transcriptional programs controlled by gene regulatory networks determine, for example, the cell fate during embryonic development. Therefore, an essential aspect of development in these organisms is the ability to maintain the functionality of their genetic developmental programs even in the presence of genetic variation, changing environmental conditions and biochemical noise, a property commonly termed robustness. We discuss the implication of different molecular mechanisms of robustness involved in neurodevelopment, which is characterized by the interplay of many developmental programs at a molecular, cellular and systemic level. We specifically focus on processes affecting the function of gene regulatory networks, encompassing transcriptional regulatory elements and post-transcriptional processes such as miRNA-based regulation, but also higher order regulatory organization, such as gene network topology. We also present cases where impairment of robustness mechanisms can be associated with neurodevelopmental disorders, as well as reasons why understanding these mechanisms should represent an important part of the study of gene regulatory networks driving neural development.

Emergence of cooperative bistability and robustness of gene regulatory networks.

Gene regulatory networks (GRNs) are complex systems in which many genes regulate mutually to adapt the cell state to environmental conditions. In addition to function, the GRNs possess several kinds of robustness. This robustness means that systems do not lose their functionality when exposed to disturbances such as mutations or noise, and is widely observed at many levels in living systems. Both function and robustness have been acquired through evolution. In this respect, GRNs utilized in living systems are rare among all possible GRNs. In this study, we explored the fitness landscape of GRNs and investigated how robustness emerged in highly-fit GRNs. We considered a toy model of GRNs with one input gene and one output gene. The difference in the expression level of the output gene between two input states, "on" and "off", was considered as fitness. Thus, the determination of the fitness of a GRN was based on how sensitively it responded to the input. We employed the multicanonical Monte Carlo method, which can sample GRNs randomly in a wide range of fitness levels, and classified the GRNs according to their fitness. As a result, the following properties were found: (1) Highly-fit GRNs exhibited bistability for intermediate input between "on" and "off". This means that such GRNs responded to two input states by using different fixed points of dynamics. This bistability emerges necessarily as fitness increases. (2) These highly-fit GRNs were robust against noise because of their bistability. In other words, noise robustness is a byproduct of high fitness. (3) GRNs that were robust against mutations were not extremely rare among the highly-fit GRNs. This implies that mutational robustness is readily acquired through the evolutionary process. These properties are universal irrespective of the evolutionary pathway, because the results do not rely on evolutionary simulation.

Changes of signal transductivity and robustness of gene regulatory network in the carcinogenesis of leukemic subtypes via microarray sample data.

Mutation accumulation and epigenetic alterations in genes are important for carcinogenesis. Because leukemogenesis-related signal pathways have been investigated and microarray sample data have been produced in acute myeloid leukemia (AML), myelodysplastic syndromes (MDS) and normal cells, systems analysis in coupling pathways becomes possible.Based on system modeling and identification, we could construct the coupling pathways and their associated gene regulatory networks using microarray sample data. By applying system theory to the estimated system model in coupling pathways, we can then obtain transductivity sensitivity, basal sensitivity and error sensitivity of each protein to identify the potential impact of genetic mutations, epigenetic alterations and the coupling of other pathways from the perspective of energy, respectively. By comparing the results in AML, MDS and normal cells, we investigated the potential critical genetic mutations and epigenetic alterations that activate or repress specific cellular functions to promote MDS or AML leukemogenesis. We suggested that epigenetic modification of β-catenin and signal integration of CSLs, AP-2α, STATs, c-Jun and β-catenin could contribute to cell proliferation at AML and MDS. Epigenetic regulation of ERK and genetic mutation of p53 could lead to the repressed apoptosis, cell cycle arrest and DNA repair in leukemic cells. Genetic mutation of JAK, epigenetic regulation of ERK, and signal integration of C/EBPα could result in the promotion of MDS cell differentiation. According to the results, we proposed three drugs, decitabine, genistein, and monorden for preventing AML leukemogenesis, while three drugs, decitabine, thalidomide, and geldanamycin, for preventing MDS leukemogenesis.

Selection for distinct gene expression properties favours the evolution of mutational robustness in gene regulatory networks.

Mutational robustness is a genotype's tendency to keep a phenotypic trait with little and few changes in the face of mutations. Mutational robustness is both ubiquitous and evolutionarily important as it affects in different ways the probability that new phenotypic variation arises. Understanding the origins of robustness is specially relevant for systems of development that are phylogenetically widespread and that construct phenotypic traits with a strong impact on fitness. Gene regulatory networks are examples of this class of systems. They comprise sets of genes that, through cross-regulation, build the gene activity patterns that define cellular responses, different tissues or distinct cell types. Several empirical observations, such as a greater robustness of wild-type phenotypes, suggest that stabilizing selection underlies the evolution of mutational robustness. However, the role of selection in the evolution of robustness is still under debate. Computer simulations of the dynamics and evolution of gene regulatory networks have shown that selection for any gene activity pattern that is steady and self-sustaining is sufficient to promote the evolution of mutational robustness. Here, I generalize this scenario using a computational model to show that selection for different aspects of a gene activity phenotype increases mutational robustness. Mutational robustness evolves even when selection favours properties that conflict with the stationarity of a gene activity pattern. The results that I present support an important role for stabilizing selection in the evolution of robustness in gene regulatory networks.

Model checking the evolution of gene regulatory networks

The behaviour of gene regulatory networks (GRNs) is typically analysed using simulation-based statistical testing-like methods. In this paper, we demonstrate that we can replace this approach by a formal verification-like method that gives higher assurance and scalability. We focus on Wagner’s weighted GRN model with varying weights, which is used in evolutionary biology. In the model, weight parameters represent the gene interaction strength that may change due to genetic mutations. For a property of interest, we synthesise the constraints over the parameter space that represent the set of GRNs satisfying the property. We experimentally show that our parameter synthesis procedure computes the mutational robustness of GRNs—an important problem of interest in evolutionary biology—more efficiently than the classical simulation method. We specify the property in linear temporal logic. We employ symbolic bounded model checking and SMT solving to compute the space of GRNs that satisfy the property, which amounts to synthesizing a set of linear constraints on the weights.

Evolving robust gene regulatory networks.

Design and implementation of robust network modules is essential for construction of complex biological systems through hierarchical assembly of ‘parts’ and ‘devices’. The robustness of gene regulatory networks (GRNs) is ascribed chiefly to the underlying topology. The automatic designing capability of GRN topology that can exhibit robust behavior can dramatically change the current practice in synthetic biology. A recent study shows that Darwinian evolution can gradually develop higher topological robustness. Subsequently, this work presents an evolutionary algorithm that simulates natural evolution in silico, for identifying network topologies that are robust to perturbations. We present a Monte Carlo based method for quantifying topological robustness and designed a fitness approximation approach for efficient calculation of topological robustness which is computationally very intensive. The proposed framework was verified using two classic GRN behaviors: oscillation and bistability, although the framework is generalized for evolving other types of responses. The algorithm identified robust GRN architectures which were verified using different analysis and comparison. Analysis of the results also shed light on the relationship among robustness, cooperativity and complexity. This study also shows that nature has already evolved very robust architectures for its crucial systems; hence simulation of this natural process can be very valuable for designing robust biological systems.

What are memories made of? How Polycomb and Trithorax proteins mediate epigenetic memory

In any biological system with memory, the state of the system depends on its history. Epigenetic memory maintains gene expression states through cell generations without a change in DNA sequence and in the absence of initiating signals. It is immensely powerful in biological systems - it adds long-term stability to gene expression states and increases the robustness of gene regulatory networks. The Polycomb group (PcG) and Trithorax group (TrxG) proteins can confer long-term, mitotically heritable memory by sustaining silent and active gene expression states, respectively. Several recent studies have advanced our understanding of the molecular mechanisms of this epigenetic memory during DNA replication and mitosis.

The underlying molecular and network level mechanisms in the evolution of robustness in gene regulatory networks.

Gene regulatory networks show robustness to perturbations. Previous works identified robustness as an emergent property of gene network evolution but the underlying molecular mechanisms are poorly understood. We used a multi-tier modeling approach that integrates molecular sequence and structure information with network architecture and population dynamics. Structural models of transcription factor-DNA complexes are used to estimate relative binding specificities. In this model, mutations in the DNA cause changes on two levels: (a) at the sequence level in individual binding sites (modulating binding specificity), and (b) at the network level (creating and destroying binding sites). We used this model to dissect the underlying mechanisms responsible for the evolution of robustness in gene regulatory networks. Results suggest that in sparse architectures (represented by short promoters), a mixture of local-sequence and network-architecture level changes are exploited. At the local-sequence level, robustness evolves by decreasing the probabilities of both the destruction of existent and generation of new binding sites. Meanwhile, in highly interconnected architectures (represented by long promoters), robustness evolves almost entirely via network level changes, deleting and creating binding sites that modify the network architecture.

Computational Systems Biology

Computational Systems Biology

Robustness of the hypoxic response: Another job for miRNAs?

Living organisms are constantly exposed to environmental and genetic perturbations. Biological robustness enables these organisms to maintain their functional stability in the presence of external or internal changes. It has been proposed that microRNAs (miRNAs), small non-coding regulatory RNAs, contribute to robustness of gene regulatory networks. The hypoxic response is a major and well-characterized example of a cellular and systemic response to environmental stress that needs to be robust. miRNAs regulate the response to hypoxia, both at the level of the main transcription factor that mediates this response, the hypoxia-inducible factor (HIF), and at the level of one of the most important systemic outcomes of the response: angiogenesis. In this review, we will take the hypoxic response as a paradigm of miRNAs participating in circuits that provide robustness to biological responses.

Mutational Robustness of Gene Regulatory Networks

Mutational robustness of gene regulatory networks refers to their ability to generate constant biological output upon mutations that change network structure. Such networks contain regulatory interactions (transcription factor – target gene interactions) but often also protein-protein interactions between transcription factors. Using computational modeling, we study factors that influence robustness and we infer several network properties governing it. These include the type of mutation, i.e. whether a regulatory interaction or a protein-protein interaction is mutated, and in the case of mutation of a regulatory interaction, the sign of the interaction (activating vs. repressive). In addition, we analyze the effect of combinations of mutations and we compare networks containing monomeric with those containing dimeric transcription factors. Our results are consistent with available data on biological networks, for example based on evolutionary conservation of network features. As a novel and remarkable property, we predict that networks are more robust against mutations in monomer than in dimer transcription factors, a prediction for which analysis of conservation of DNA binding residues in monomeric vs. dimeric transcription factors provides indirect evidence.

In SilicoEvolution of Gene Cooption in Pattern-Forming Gene Networks

Gene recruitment or cooption occurs when a gene, which may be part of an existing gene regulatory network (GRN), comes under the control of a new regulatory system. Such re-arrangement of pre-existing networks is likely more common for increasing genomic complexity than the creation of new genes. Using evolutionary computations (EC), we investigate how cooption affects the evolvability, outgrowth and robustness of GRNs. We use a data-driven model of insect segmentation, for the fruit fly Drosophila, and evaluate fitness by robustness to maternal variability—a major constraint in biological development. We compare two mechanisms of gene cooption: a simpler one with gene Introduction and Withdrawal operators; and one in which GRN elements can be altered by transposon infection. Starting from a minimal 2-gene network, insufficient for fitting the Drosophila gene expression patterns, we find a general trend of coopting available genes into the GRN, in order to better fit the data. With the transposon mechanism, we find co-evolutionary oscillations between genes and their transposons. These oscillations may offer a new technique in EC for overcoming premature convergence. Finally, we comment on how a differential equations (in contrast to Boolean) approach is necessary for addressing realistic continuous variation in biochemical parameters.

On the Interplay between Entropy and Robustness of Gene Regulatory Networks

The interplay between entropy and robustness of gene network is a core mechanism of systems biology. The entropy is a measure of randomness or disorder of a physical system due to random parameter fluctuation and environmental noises in gene regulatory networks. The robustness of a gene regulatory network, which can be measured as the ability to tolerate the random parameter fluctuation and to attenuate the effect of environmental noise, will be discussed from the robust H∞ stabilization and filtering perspective. In this review, we will also discuss their balancing roles in evolution and potential applications in systems and synthetic biology.

Modeling stochasticity and robustness in gene regulatory networks

Motivation: Understanding gene regulation in biological processes and modeling the robustness of underlying regulatory networks is an important problem that is currently being addressed by computational systems biologists. Lately, there has been a renewed interest in Boolean modeling techniques for gene regulatory networks (GRNs). However, due to their deterministic nature, it is often difficult to identify whether these modeling approaches are robust to the addition of stochastic noise that is widespread in gene regulatory processes. Stochasticity in Boolean models of GRNs has been addressed relatively sparingly in the past, mainly by flipping the expression of genes between different expression levels with a predefined probability. This stochasticity in nodes (SIN) model leads to over representation of noise in GRNs and hence non-correspondence with biological observations.Results: In this article, we introduce the stochasticity in functions (SIF) model for simulating stochasticity in Boolean models of GRNs. By providing biological motivation behind the use of the SIF model and applying it to the T-helper and T-cell activation networks, we show that the SIF model provides more biologically robust results than the existing SIN model of stochasticity in GRNs.Availability: Algorithms are made available under our Boolean modeling toolbox, GenYsis. The software binaries can be downloaded from http://si2.epfl.ch/∼garg/genysis.html.Contact: abhishek.garg@epfl.ch

Don’t Try This at Home

The most common source of genetic variation that can lead to evolution is certainly single-nucleotide mutations. But less frequent events like gene duplication are likely also important, as they allow more dramatic changes that result in rewiring of genetic regulatory networks. Isalan et al . explored in bacteria the repercussions of duplication events that would allow splicing of a regulatory promoter sequence of one gene with the open reading frame of another gene. They engineered almost 600 such events in genes encoding transcription factors and monitored effects on survival, growth, and evolvability of the bacterial strains. The remarkable robustness of gene regulatory networks was confirmed by the fact that over 95% of these changes did not prevent viability of the bacteria. Only 16% of the strains showed any change in growth, and some grew even better than wild-type bacteria. The authors also subjected the strains to selective pressures in the lab over many generations and found that even their relatively small (in evolutionary terms) sampling of 600 new networks yielded a dozen or so with apparent evolutionary advantages in fitness. Although some of the rewiring events caused new positive or negative transcriptional feedback loops, these did not cause changes in the amounts of green fluorescent protein–tagged proteins produced from the experimental genes. This suggests that the regulatory networks themselves can compensate for such a change and emphasize that posttranscriptional mechanisms should not be underestimated in comparison with transcriptional control mechanisms. Bennett and Hasty add valuable commentary, noting that you would not want to try such rewiring of processing units in your home computer (analogies between biological control mechanism and systems design in modern technology go only so far) and that advocates of "intelligent design" who claim that the complexity of biological systems would cause them to fail in the presence of randomly generated changes will have to take into account these new data to the contrary. M. Isalan, C. Lemerle, K. Michalodimitrakis, C. Horn, P. Beltrao, E. Raineri, M. Garriga-Canut, L. Serrano, Evolvability and hierarchy in rewired bacterial gene networks. Nature 452 , 840-845 (2008). [PubMed] M. R. Bennett, J. Hasty, Genome rewired. Nature 452 , 824-825 (2008). [PubMed]

Threshold-dominated regulation hides genetic variation in gene expression networks

BackgroundIn dynamical models with feedback and sigmoidal response functions, some or all variables have thresholds around which they regulate themselves or other variables. A mathematical analysis has shown that when the dose-response functions approach binary or on/off responses, any variable with an equilibrium value close to one of its thresholds is very robust to parameter perturbations of a homeostatic state. We denote this threshold robustness. To check the empirical relevance of this phenomenon with response function steepnesses ranging from a near on/off response down to Michaelis-Menten conditions, we have performed a simulation study to investigate the degree of threshold robustness in models for a three-gene system with one downstream gene, using several logical input gates, but excluding models with positive feedback to avoid multistationarity. Varying parameter values representing functional genetic variation, we have analysed the coefficient of variation (CV) of the gene product concentrations in the stable state for the regulating genes in absolute terms and compared to the CV for the unregulating downstream gene. The sigmoidal or binary dose-response functions in these models can be considered as phenomenological models of the aggregated effects on protein or mRNA expression rates of all cellular reactions involved in gene expression.ResultsFor all the models, threshold robustness increases with increasing response steepness. The CVs of the regulating genes are significantly smaller than for the unregulating gene, in particular for steep responses. The effect becomes less prominent as steepnesses approach Michaelis-Menten conditions. If the parameter perturbation shifts the equilibrium value too far away from threshold, the gene product is no longer an effective regulator and robustness is lost. Threshold robustness arises when a variable is an active regulator around its threshold, and this function is maintained by the feedback loop that the regulator necessarily takes part in and also is regulated by. In the present study all feedback loops are negative, and our results suggest that threshold robustness is maintained by negative feedback which necessarily exists in the homeostatic state.ConclusionThreshold robustness of a variable can be seen as its ability to maintain an active regulation around its threshold in a homeostatic state despite external perturbations. The feedback loop that the system necessarily possesses in this state, ensures that the robust variable is itself regulated and kept close to its threshold. Our results suggest that threshold regulation is a generic phenomenon in feedback-regulated networks with sigmoidal response functions, at least when there is no positive feedback. Threshold robustness in gene regulatory networks illustrates that hidden genetic variation can be explained by systemic properties of the genotype-phenotype map.