View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by CERN Document Server Neutral Mutations and Punctuated Equilibrium in Evolving Genetic Networks Stefan Bornholdta and Kim Sneppenb a Institut f¨ur Theoretische Physik, Universit¨at Kiel, Leibnizstr. 15, D-24098 Kiel, Germany [email protected] b NORDITA, Blegdamsvej 17, DK-2100 Copenhagen, Denmark [email protected] (August 23, 1997) Abstract We consider evolution of Boolean networks, and demonstrate how require- ment of continuity in evolution may leads to punctuated equilibrium. We discuss evolution of genetic flexibility and how this may reconcile the expo- nential distribution of species lifetimes seen by Van Valen with the power law distribution of genera lifetimes obtained when one averages over the fossil record. PACS number(s): 87.10.+e, 02.70.Lq, 05.40.+j Typeset using REVTEX 1 Evolution of life is presumably a random process with selection [1]. It has been discussed whether this process can be viewed as some hill climbing process [2], or whether evolution mostly happens as a random walk where most changes have no influence on the phenotype, and thus may be considered as neutral [3]. The case of evolution as adaptation in some externally imposed fitness landscapes have been originally proposed by Sewall Wright [2], and later formed the basis for models of punctuated equilibrium by Newman [4] and Lande [5]. The case for neutral evolution has been presented by Kimura [3], and is experimentally supported on microlevel by observation of the many functionally identically variants there are of most of the important macromolecules of life. The observation of punctuated equilibrium in the fossil record, recently discussed by Gould and Eldredge [6], may be taken as an indication that evolution of a species consists of exaptations of jumping from one hill top to another nearby in some fitness landscape. Naturally such jumps will be rare, separated by large time intervals where species are located at a fitness peak, and the resulting evolutionary pattern will show punctuations as indeed seen in the fossil record. This picture of single species evolution in a given fixed landscape has been modeled explicitly by Newman and Lande [4,5]. However, also neutral evolution may show punctuations: as, for example, might be visualized by finding the exit in a labyrinth or from finding a golf hole by random walk in a flat landscape. The picture here is that genetic changes always take place, but that the phenotypic changes only rarely occur. This has recently been demonstrated by modeling the evolution of RNA secondary structure by P. Schuster and collaborators [7]. For these molecules most base pair exchanges do not induce any changes in their secondary structure, and thus they may be considered as neutral. Occasionally, however, one mutation may suddenly lead to a complete readjustment of the structure, and presumably therefore to a major change in its functionality. In any case, as demonstrated by Bak and Sneppen [8], punctuated equilibrium on the organism level might be connected to the episodic punctuations observed on the ecosys- tem level. The crucial element of such an extrapolation is that the environment of each 2 species depends on species which are ecological neighbors, thereby allowing punctuations to propagate across the ecosystem. In the present paper we introduce a model of single species evolution without inter-species competition. The evolutionary steps in the model consists of random mutations combined with selection of mutants preserving the phenotype with respect to a given environment. Thus, the only requirement is continuity in phenotype. Other changes in genotype are allowed, creating a path of neutral mutations. We will discuss how this requirement of continuity in evolution may constrain and guide the evolution of an individual species in the face of a constantly changing environment. Our fundamental constituents are the genes of the organism, and the evolution we con- sider is on the genetic network level. Although genetic networks consist of biochemical switches [9], it has been proposed that the on-off nature of these switches can be well approximated by Boolean functions [10–12]. Thus we here consider networks of random Boolean functions. The functionality we test for is attractors of these networks [11]. We model continuity by testing for reaching a given attractor on subsequent steps, but allowing changes that modify attractors that are not tested from the actual initial condition. In sub- sequent steps, the initial condition (modeling the environment) assumes new random values such that previous neutral mutations may then surface in the phenotype. The model may be viewed as coarse grained in time, demanding continuity in fulfilling present demands on a given environment, but opening for silent changes in how the organism may respond to changed environments. We note that the continuity requirement also means that we only evolve species for which at least one neighbour genotype has quite similar dynamics. Thus our continuity requirement favours evolution of robust networks where changing of a few rules does not leads to overall changes of expression. In the simplest form we start the simulation by selecting a genetic network with N genes and assign each a Boolean variable σi = −1 or 1. Further for each of the N genes we define an updating matrix in the form of a lookup table which determines its output for each of the 2N input states from the N genes in the system. Finally we define which gene is actively 3 connected to which, by a matrix wi,j that defines the input to state i from gene number j as wi,jσj. The entry value of the connectivity matrix wi,j may take values −1, 0 and 1. A value 0 corresponds to no coupling from gene i to gene j, and is defined to enter the updating matrix by the entry value −1 always. Typically only a fraction of the connectivity matrix entries are in use, and the average number of active inputs per gene is called the connectivity K. It varies between 0 and N, meaning that K includes self couplings. Thus K =0means that all is fixed to the output state specified by input corresponding to (−1, ..., −1) to all states. Initially we start with a low but finite connectivity. In the following, an initial average connectivity of K = 1 coupling per site will be used. Boolean networks can exhibit a rich dynamical behavior, including fixed points, periodic attractors and long transients to reach these. Further the number of attractors, their length, and the length of the transients to reach these are known to depend strongly on the coordi- nation number [13]. In this paper we do not address any questions connected with the time scale of these attractors. Instead we consider a longer evolutionary time scale where the geometry of the network may change. This reminds of evolutionary time scales of biological species as compared with the much shorter lifetimes of organisms. The evolutionary time step of the system is: 1. Select a random initial state of the system {σi}. Let your mother system evolve from this state until a final attractor is determined. 2. Create a daughter network by a) adding, b) removing, or c) adding and removing a weight in the coupling matrix at random, each option occurring with probability =1/3. Let you daughter system evolve from the same initial state as that selected for the mother and test whether it reaches the same attractor as the mother system did. In case it does then replace mother with daughter network and go to step 3. In case another attractor is reached, keep mother network and go to step 3. 3. Then finally one random bit of the total N ∗ 2N lookup table entries is flipped to 4 another value. This allows for a convenient self averaging of the system, and in fact represents a very slow change. Iterating these steps makes an evolutionary algorithm that represents the sole require- ment of continuity in evolution and how this may proceed under an environment that fluctu- ates. No selective pressure is applied. However, by not testing the entire basin structure of the network we allow for silent mutations in the organism that may surface in a later changed environment. In our model this new environment is represented by the new input vector selected in the next time step. In practice this new environment may well represent either physical changes in the environment or maybe changes in neighbor species in an ecosystem, that we do not attempt to model here. In figure 1a we show how the connectivity of this system evolves with time in a network of size N = 16. One observes that the typical connectivity of the network is confined to lower values than for random networks. This is further quantified in Figure 2 where the distributions of average connectivities are displayed in the statistically stationary state. All data are taken beginning after an equilibration time of ten percent of the length of the total run. Notice that there are two distributions, one counting the frequency of connectivities for all new “species” and one counting the time averaged distribution of connectivities. These two distributions diverge strongly for high connectivities because the few species with high connectivity have very long lifetimes, i.e., it is very difficult to find mutations which do not change the activity pattern of the networks completely. In our case, the activity pattern consists of the transient and the final periodic attractor following the given initial state. The time scale of these patterns becomes large for networks with high connectivity, making it more difficult to keep the exact dynamic pattern under the mutation of a weight.