Mutation rules and the evolution of sparseness and modularity in biological systems

Revisiting the evolution of bow-tie architecture in signaling networks

Bow-tie architecture is a layered network structure that has a narrow middle layer with multiple inputs and outputs. Such structures are widely seen in the molecular networks in cells, suggesting that a universal evolutionary mechanism underlies the emergence of bow-tie architecture. The previous theoretical studies have implemented evolutionary simulations of the feedforward network to satisfy a given input-output goal and proposed that the bow-tie architecture emerges when the ideal input-output relation is given as a rank-deficient matrix with mutations in network link intensities in a multiplicative manner. Here, we report that the bow-tie network inevitably appears when the link intensities representing molecular interactions are small at the initial condition of the evolutionary simulation, regardless of the rank of the goal matrix. Our dynamical system analysis clarifies the mechanisms underlying the emergence of the bow-tie structure. Further, we demonstrate that the increase in the input-output matrix reduces the width of the middle layer, resulting in the emergence of bow-tie architecture, even when evolution starts from large link intensities. Our data suggest that bow-tie architecture emerges as a side effect of evolution rather than as a result of evolutionary adaptation.

Modularity in Biological Networks

Network modeling, from the ecological to the molecular scale has become an essential tool for studying the structure, dynamics and complex behavior of living systems. Graph representations of the relationships between biological components open up a wide variety of methods for discovering the mechanistic and functional properties of biological systems. Many biological networks are organized into a modular structure, so methods to discover such modules are essential if we are to understand the biological system as a whole. However, most of the methods used in biology to this end, have a limited applicability, as they are very specific to the system they were developed for. Conversely, from the statistical physics and network science perspective, graph modularity has been theoretically studied and several methods of a very general nature have been developed. It is our perspective that in particular for the modularity detection problem, biology and theoretical physics/network science are less connected than they should. The central goal of this review is to provide the necessary background and present the most applicable and pertinent methods for community detection in a way that motivates their further usage in biological research.

The relation between crosstalk and gene regulation form revisited

Genes differ in the frequency at which they are expressed and in the form of regulation used to control their activity. In particular, positive or negative regulation can lead to activation of a gene in response to an external signal. Previous works proposed that the form of regulation of a gene correlates with its frequency of usage: positive regulation when the gene is frequently expressed and negative regulation when infrequently expressed. Such network design means that, in the absence of their regulators, the genes are found in their least required activity state, hence regulatory intervention is often necessary. Due to the multitude of genes and regulators, spurious binding and unbinding events, called “crosstalk”, could occur. To determine how the form of regulation affects the global crosstalk in the network, we used a mathematical model that includes multiple regulators and multiple target genes. We found that crosstalk depends non-monotonically on the availability of regulators. Our analysis showed that excess use of regulation entailed by the formerly suggested network design caused high crosstalk levels in a large part of the parameter space. We therefore considered the opposite ‘idle’ design, where the default unregulated state of genes is their frequently required activity state. We found, that ‘idle’ design minimized the use of regulation and thus minimized crosstalk. In addition, we estimated global crosstalk of S. cerevisiae using transcription factors binding data. We demonstrated that even partial network data could suffice to estimate its global crosstalk, suggesting its applicability to additional organisms. We found that S. cerevisiae estimated crosstalk is lower than that of a random network, suggesting that natural selection reduces crosstalk. In summary, our study highlights a new type of protein production cost which is typically overlooked: that of regulatory interference caused by the presence of excess regulators in the cell. It demonstrates the importance of whole-network descriptions, which could show effects missed by single-gene models.

Genes differ in the frequency at which they are expressed and in the form of regulation used to control their activity. The basic level of regulation is mediated by different types of DNA-binding proteins, where each type regulates particular gene(s). We distinguish between two basic forms of regulation: positive—if a gene is activated by the binding of its regulatory protein, and negative—if it is active unless bound by its regulatory protein. Due to the multitude of genes and regulators, spurious binding and unbinding events, called “crosstalk”, could occur. How does the form of regulation, positive or negative, affect the extent of regulatory crosstalk? To address this question, we used a mathematical model integrating many genes and many regulators. As intuition suggests, we found that in most of the parameter space, crosstalk increased with the availability of regulators. We propose, that crosstalk is usually reduced when networks are designed such that minimal regulation is needed, which we call the ‘idle’ design. In other words: a frequently needed gene will use negative regulation and conversely, a scarcely needed gene will employ positive regulation. In both cases, the requirement for the regulators is minimized. In addition, we demonstrate how crosstalk can be calculated from available datasets and discuss the technical challenges in such calculation, specifically data incompleteness.

Genotype Components as Predictors of Phenotype in Model Gene Regulatory Networks

Models of gene regulatory networks (GRN) have proven useful for understanding many aspects of the highly complex behavior of biological control networks. Randomly generated non-Boolean networks were used in experimental simulations to generate data on dynamic phenotypes as a function of several genotypic parameters. We found that predictive relationships between some phenotypes and quantitative genotypic parameters such as number of network genes, interaction density, and initial condition could be derived depending on the strength of the topological (positional) genotype on specific phenotypes. We quantitated the strength of the topological genotype effect (TGE) on a number of phenotypes in multi-gene networks. For phenotypes with a low influence of topological genotype, derived and empirical relationships using quantitative genotype parameters were accurate in phenotypic outcomes. We found a number of dynamic network properties, including oscillation behaviors, that were largely dependent on genotype topology, and for which no such general quantitative relationships were determinable. It remains to be determined if these results are applicable to biological gene regulatory networks.

Quantitative Operating Principles of Yeast Metabolism during Adaptation to Heat Stress

Microorganisms evolved adaptive responses to survive stressful challenges in ever-changing environments. Understanding the relationships between the physiological/metabolic adjustments allowing cellular stress adaptation and gene expression changes being used by organisms to achieve such adjustments may significantly impact our ability to understand and/or guide evolution. Here, we studied those relationships during adaptation to various stress challenges in Saccharomyces cerevisiae, focusing on heat stress responses. We combined dozens of independent experiments measuring whole-genome gene expression changes during stress responses with a simplified kinetic model of central metabolism. We identified alternative quantitative ranges for a set of physiological variables in the model (production of ATP, trehalose, NADH, etc.) that are specific for adaptation to either heat stress or desiccation/rehydration. Our approach is scalable to other adaptive responses and could assist in developing biotechnological applications to manipulate cells for medical, biotechnological, or synthetic biology purposes.

Evolutionary trade-offs and the structure of polymorphisms

Populations of organisms show genetic differences called polymorphisms. Understanding the effects of polymorphisms is important for biology and medicine. Here, we ask which polymorphisms occur at high frequency when organisms evolve under trade-offs between multiple tasks. Multiple tasks present a problem, because it is not possible to be optimal at all tasks simultaneously and hence compromises are necessary. Recent work indicates that trade-offs lead to a simple geometry of phenotypes in the space of traits: phenotypes fall on the Pareto front, which is shaped as a polytope: a line, triangle, tetrahedron etc. The vertices of these polytopes are the optimal phenotypes for a single task. Up to now, work on this Pareto approach has not considered its genetic underpinnings. Here, we address this by asking how the polymorphism structure of a population is affected by evolution under trade-offs. We simulate a multi-task selection scenario, in which the population evolves to the Pareto front: the line segment between two archetypes or the triangle between three archetypes. We find that polymorphisms that become prevalent in the population have pleiotropic phenotypic effects that align with the Pareto front. Similarly, epistatic effects between prevalent polymorphisms are parallel to the front. Alignment with the front occurs also for asexual mating. Alignment is reduced when drift or linkage is strong, and is replaced by a more complex structure in which many perpendicular allele effects cancel out. Aligned polymorphism structure allows mating to produce offspring that stand a good chance of being optimal multi-taskers in at least one of the locales available to the species.This article is part of the theme issue 'Self-organization in cell biology'.

The genetics of phenotypic plasticity. XVI. Interactions among traits and the flow of information

Although the environment varies, adaptive trait plasticity is uncommon, which can be due to either costs or limitations. Currently there is little evidence for costs of plasticity; limitations are a more promising explanation, including information reliability. A possible cause for a decrease in information reliability is the channeling of environmental information through one trait that then affects the phenotype of a second trait, the information path. Using an individual-based simulation model, I explored the ways in which configurations of trait interactions and patterns of environmental variation in space and time affect the evolution of phenotypic plasticity. I found that genotypes and phenotypes evolved to shorten and simplify the information path from the environment to fitness. A shortened path was characterized by a decrease in the amount of plasticity for traits that had a less direct connection between the environment of development and fitness. A simplified path was characterized by a decrease in the amount of plasticity for traits that had multiple paths between the environment and their phenotype. These results suggest that an eighth proposition be added to the theory of the evolution of phenotypic plasticity: Trait plasticity will evolve to minimize the information path between the environment and fitness.

On the role of sparseness in the evolution of modularity in gene regulatory networks

Modularity is a widespread property in biological systems. It implies that interactions occur mainly within groups of system elements. A modular arrangement facilitates adjustment of one module without perturbing the rest of the system. Therefore, modularity of developmental mechanisms is a major factor for evolvability, the potential to produce beneficial variation from random genetic change. Understanding how modularity evolves in gene regulatory networks, that create the distinct gene activity patterns that characterize different parts of an organism, is key to developmental and evolutionary biology. One hypothesis for the evolution of modules suggests that interactions between some sets of genes become maladaptive when selection favours additional gene activity patterns. The removal of such interactions by selection would result in the formation of modules. A second hypothesis suggests that modularity evolves in response to sparseness, the scarcity of interactions within a system. Here I simulate the evolution of gene regulatory networks and analyse diverse experimentally sustained networks to study the relationship between sparseness and modularity. My results suggest that sparseness alone is neither sufficient nor necessary to explain modularity in gene regulatory networks. However, sparseness amplifies the effects of forms of selection that, like selection for additional gene activity patterns, already produce an increase in modularity. That evolution of new gene activity patterns is frequent across evolution also supports that it is a major factor in the evolution of modularity. That sparseness is widespread across gene regulatory networks indicates that it may have facilitated the evolution of modules in a wide variety of cases.

Modular systems have performance and design advantages over non-modular systems. Thus, modularity is very important for the development of a wide range of new technological or clinical applications. Moreover, modularity is paramount to evolutionary biology since it allows adjusting one organismal function without disturbing other previously evolved functions. But how does modularity itself evolve? Here I analyse the structure of regulatory networks and follow simulations of network evolution to study two hypotheses for the origin of modules in gene regulatory networks. The first hypothesis considers that sparseness, a low number of interactions among the network genes, could be responsible for the evolution of modular networks. The second, that modules evolve when selection favours the production of additional gene activity patterns. I found that sparseness alone is neither sufficient nor necessary to explain modularity in gene regulatory networks. However, it enhances the effects of selection for multiple gene activity patterns. While selection for multiple patterns may be decisive in the evolution of modularity, that sparseness is widespread across gene regulatory networks suggests that its contributions should not be neglected.

Optimal pathways for the assembly of the Apaf-1·cytochrome c complex into apoptosome

The formation of a heptameric apoptosome is a crucial event in the intrinsic cell death pathway. Considerable progress has been made towards unraveling the constituents and the structure of the apoptosome as well as the mechanism of apoptosome-mediated caspase-9 activation. However, a significant gap remains in the understanding of this process, i.e., how seven Apaf-1·cytochrome c complexes stepwisely assemble into an apoptosome. Here, we construct a biophysical model that incorporates current biochemical knowledge about the formation of apoptosome. We propose 11 elementary routes and enumerate all 2047 possible assembly pathways from the Apaf-1·cytochrome c complex to the heptameric apoptosome. By combining mathematical analysis and numerical simulation, we find that two elementary routes are the most favorable biochemical reaction routes and there are 52 optimal assembly pathways which are economical and relatively fast. Our study yields the first comprehensive analysis of apoptosome assembly and provides insights into complex assembly pathways.

Endocytosis as a stabilizing mechanism for tissue homeostasis

Cells in tissues communicate by secreted growth factors (GF) and other signals. An important function of cell circuits is tissue homeostasis: maintaining proper balance between the amounts of different cell types. Homeostasis requires negative feedback on the GFs, to avoid a runaway situation in which cells stimulate each other and grow without control. Feedback can be obtained in at least two ways: endocytosis in which a cell removes its cognate GF by internalization and cross-inhibition in which a GF down-regulates the production of another GF. Here we ask whether there are design principles for cell circuits to achieve tissue homeostasis. We develop an analytically solvable framework for circuits with multiple cell types and find that feedback by endocytosis is far more robust to parameter variation and has faster responses than cross-inhibition. Endocytosis, which is found ubiquitously across tissues, can even provide homeostasis to three and four communicating cell types. These design principles form a conceptual basis for how tissues maintain a healthy balance of cell types and how balance may be disrupted in diseases such as degeneration and fibrosis.

Genetic code asymmetry supports diversity through experimentation with posttranslational modifications

Protein N-glycosylation has been identified in all three domains of life presumably conserved for its early role in glycoprotein folding. However, the N-glycans added to proteins in the secretory pathway of multicellular organisms are remodeling in the Golgi, increasing structural diversity exponentially and adding new layers of functionality in immunity, metabolism and other systems. The branching and elongation of N-glycan chains found on cell surface receptors generates a gradation of affinities for carbohydrate-binding proteins, the galectin, selectin and siglec families. These interactions adapt cellular responsiveness to environmental conditions, but their complexity presents a daunting challenge to drug design. To gain further insight, I review how N-glycans biosynthesis and biophysical properties provide a selective advantage in the form of tunable and ultrasensitive stimulus-response relationships. In addition, the N-glycosylation motif favors step-wise mutational experimentation with sites. Glycoproteins display accelerated evolution during vertebrate radiation, and the encoding asymmetry of NXS/T(X≠P) has left behind phylogenetic evidence suggesting that the genetic code may have been selected to optimize diversity in part through emerging posttranslational modifications.

How evolution learns to generalise: Using the principles of learning theory to understand the evolution of developmental organisation

One of the most intriguing questions in evolution is how organisms exhibit suitable phenotypic variation to rapidly adapt in novel selective environments. Such variability is crucial for evolvability, but poorly understood. In particular, how can natural selection favour developmental organisations that facilitate adaptive evolution in previously unseen environments? Such a capacity suggests foresight that is incompatible with the short-sighted concept of natural selection. A potential resolution is provided by the idea that evolution may discover and exploit information not only about the particular phenotypes selected in the past, but their underlying structural regularities: new phenotypes, with the same underlying regularities, but novel particulars, may then be useful in new environments. If true, we still need to understand the conditions in which natural selection will discover such deep regularities rather than exploiting ‘quick fixes’ (i.e., fixes that provide adaptive phenotypes in the short term, but limit future evolvability). Here we argue that the ability of evolution to discover such regularities is formally analogous to learning principles, familiar in humans and machines, that enable generalisation from past experience. Conversely, natural selection that fails to enhance evolvability is directly analogous to the learning problem of over-fitting and the subsequent failure to generalise. We support the conclusion that evolving systems and learning systems are different instantiations of the same algorithmic principles by showing that existing results from the learning domain can be transferred to the evolution domain. Specifically, we show that conditions that alleviate over-fitting in learning systems successfully predict which biological conditions (e.g., environmental variation, regularity, noise or a pressure for developmental simplicity) enhance evolvability. This equivalence provides access to a well-developed theoretical framework from learning theory that enables a characterisation of the general conditions for the evolution of evolvability.

A striking feature of evolving organisms is their ability to acquire novel characteristics that help them adapt in new environments. The origin and the conditions of such ability remain elusive and is a long-standing question in evolutionary biology. Recent theory suggests that organisms can evolve designs that help them generate novel features that are more likely to be beneficial. Specifically, this is possible when the environments that organisms are exposed to share common regularities. However, the organisms develop robust designs that tend to produce what had been selected in the past and might be inflexible for future environments. The resolution comes from a recent theory introduced by Watson and Szathmáry that suggests a deep analogy between learning and evolution. Accordingly, here we utilise learning theory to explain the conditions that lead to more evolvable designs. We successfully demonstrate this by equating evolvability to the way humans and machines generalise to previously-unseen situations. Specifically, we show that the same conditions that enhance generalisation in learning systems have biological analogues and help us understand why environmental noise and the reproductive and maintenance costs of gene-regulatory connections can lead to more evolvable designs.

Evolution of Modular Networks Under Selection for Non-Linearly Denoising

Many biological and artifact networks often represent modular structures in which the network can be decomposed into several subnetworks. Here, we propose a simple model for the modular network evolution based on the nonlinear denoising in node activities. This model suggests that modular networks can evolve under certain conditions — if the stipulated goals for the networks or the input and target output pairs involve modular features, or if the signal transfer in a node is carried out in a nonlinear manner with respect to the saturation at the upper and lower bounds. Our model highlights the positive role played by noise in modular network evolution.

Habitat variability does not generally promote metabolic network modularity in flies and mammals

The evolution of species habitat range is an important topic over a wide range of research fields. In higher organisms, habitat range evolution is generally associated with genetic events such as gene duplication. However, the specific factors that determine habitat variability remain unclear at higher levels of biological organization (e.g., biochemical networks). One widely accepted hypothesis developed from both theoretical and empirical analyses is that habitat variability promotes network modularity; however, this relationship has not yet been directly tested in higher organisms. Therefore, I investigated the relationship between habitat variability and metabolic network modularity using compound and enzymatic networks in flies and mammals. Contrary to expectation, there was no clear positive correlation between habitat variability and network modularity. As an exception, the network modularity increased with habitat variability in the enzymatic networks of flies. However, the observed association was likely an artifact, and the frequency of gene duplication appears to be the main factor contributing to network modularity. These findings raise the question of whether or not there is a general mechanism for habitat range expansion at a higher level (i.e., above the gene scale). This study suggests that the currently widely accepted hypothesis for habitat variability should be reconsidered.

Evolution of sparsity and modularity in a model of protein allostery

The sequence of a protein is not only constrained by its physical and biochemical properties under current selection, but also by features of its past evolutionary history. Understanding the extent and the form that these evolutionary constraints may take is important to interpret the information in protein sequences. To study this problem, we introduce a simple but physical model of protein evolution where selection targets allostery, the functional coupling of distal sites on protein surfaces. This model shows how the geometrical organization of couplings between amino acids within a protein structure can depend crucially on its evolutionary history. In particular, two scenarios are found to generate a spatial concentration of functional constraints: high mutation rates and fluctuating selective pressures. This second scenario offers a plausible explanation for the high tolerance of natural proteins to mutations and for the spatial organization of their least tolerant amino acids, as revealed by sequence analysis and mutagenesis experiments. It also implies a faculty to adapt to new selective pressures that is consistent with observations. The model illustrates how several independent functional modules may emerge within the same protein structure, depending on the nature of past environmental fluctuations. Our model thus relates the evolutionary history of proteins to the geometry of their functional constraints, with implications for decoding and engineering protein sequences.

Evolution of bow-tie architectures in biology

Bow-tie or hourglass structure is a common architectural feature found in many biological systems. A bow-tie in a multi-layered structure occurs when intermediate layers have much fewer components than the input and output layers. Examples include metabolism where a handful of building blocks mediate between multiple input nutrients and multiple output biomass components, and signaling networks where information from numerous receptor types passes through a small set of signaling pathways to regulate multiple output genes. Little is known, however, about how bow-tie architectures evolve. Here, we address the evolution of bow-tie architectures using simulations of multi-layered systems evolving to fulfill a given input-output goal. We find that bow-ties spontaneously evolve when the information in the evolutionary goal can be compressed. Mathematically speaking, bow-ties evolve when the rank of the input-output matrix describing the evolutionary goal is deficient. The maximal compression possible (the rank of the goal) determines the size of the narrowest part of the network—that is the bow-tie. A further requirement is that a process is active to reduce the number of links in the network, such as product-rule mutations, otherwise a non-bow-tie solution is found in the evolutionary simulations. This offers a mechanism to understand a common architectural principle of biological systems, and a way to quantitate the effective rank of the goals under which they evolved.

Many biological systems show bow-tie (also called hourglass) architecture. A bow-tie means that a large number of inputs are converted to a small number of intermediates, which then fan out to generate a large number of outputs. For example, cells use a wide variety of nutrients; process them into 12 metabolic precursors, which are then used to make all of the cells biomass. Similar principles exist in biological signaling and in the information processing in the visual system. Despite the ubiquity of bow-tie structures in biology, there is no explanation of how they evolved. Here, we find that bow-ties spontaneously evolve when the information in the evolutionary goal they evolved to satisfy can be compressed. Mathematically, this means that the matrix representing the goal has deficient rank. The maximal compression possible determines the width of the bow-tie—the narrowest part in the network (equal to the rank of the goal matrix). This offers a mechanism to understand a common architectural principle of biological systems, and a way to quantitate the rank of the goals under which they evolved.

Climatic seasonality may affect ecological network structure: food webs and mutualistic networks

Ecological networks exhibit non-random structural patterns, such as modularity and nestedness, which determine ecosystem stability with species diversity and connectance. Such structure-stability relationships are well known. However, another important perspective is less well understood: the relationship between the environment and structure. Inspired by theoretical studies that suggest that network structure can change due to environmental variability, we collected data on a number of empirical food webs and mutualistic networks and evaluated the effect of climatic seasonality on ecological network structure. As expected, we found that climatic seasonality affects ecological network structure. In particular, an increase in modularity due to climatic seasonality was observed in food webs; however, it is debatable whether this occurs in mutualistic networks. Interestingly, the type of climatic seasonality that affects network structure differs with ecosystem type. Rainfall and temperature seasonality influence freshwater food webs and mutualistic networks, respectively; food webs are smaller, and more modular, with increasing rainfall seasonality. Mutualistic networks exhibit a higher diversity (particularly of animals) with increasing temperature seasonality. These results confirm the theoretical prediction that the stability increases with greater perturbation. Although these results are still debatable because of several limitations in the data analysis, they may enhance our understanding of environment-structure relationships.