Identical populations that diverge, since they are not allowed to interbreed (Figure 1). For the details of the simulation we refer the reader to the Methods Bay 41-4109 web section and Additional file 1.Population fitness increases stochasticallyResults Our model simulates evolution and hybrid formation for diploid organisms, using gene regulatory networks and their response to environmental inputs. On the individual level, we simulated the genetic and regulatory characteristics of each diploid individual as a Boolean network. Each phenotype is defined as the regulatory steady state that is reached by a network in response to inputs from outside the network (environmental factors). We calculate fitness as the similarity of the network phenotype (dynamically stable state) to the ideal phenotype under certain environmental conditions.At the beginning of the simulation the population fitness is close to zero, because the PubMed ID:https://www.ncbi.nlm.nih.gov/pubmed/28607003 network dynamics of na e individuals are entirely random. Whilst we generate networks to display certain quasi-biological characteristics in their topology – such as a scale-free-like degree distribution [43] and clustering [44] – the state of each node is just as likely to be the correct response to a given set of environmental factors as it is to be the opposite. Since fitness is averaged over nodes responding correctly and nodes responding incorrectly, they tend to cancel out in a network that has not undergone adaptation, leading to the observed fitness measure around zero. However, after mutations through evolutionary time, individuals are produced with a response closer to the optimal one. These individuals are selected to produce offspring, which increases the average fitness of the population over time. As can be seen in Figure 2, the population fitness increases through evolutionary time, with a characteristic pattern of punctuated equilibria. This feature of the graph represents the fact that improvements in fitness are brought about by rare but important mutations that produce sudden increases, alternated by periods of more or less constant fitness where only neutral or deleterious mutations occur. The latter are quickly removed from the population by selection, causing them to have little effect on the fitness of the overall population. In the very beginning, all individuals in the population are identical, but during the adaptation period, the diversity increases. This causes the population fitness to become more variable over time in some runs of the simulation.Hybrids display higher fitness than intra-population crossesAfter the networks have reached a fitness plateau with only infrequent increases in fitness (typically after 1000 generations), we duplicate the population and continueEmmrich et al. BMC Evolutionary Biology (2015) 15:Page 4 ofFigure 1 Scheme of the evolutionary algorithm and the hybrid formation. The initial network (A) is duplicated to form the initial population (B), which evolves over hundreds of generations by mutation, reassortment and selection, to form an adapted population (C). This population is duplicated, generating two initially identical populations (D). Both populations carry on evolving, but do not interbreed, causing them to diverge (indicated by the orange and blue colors) (E). At regular intervals, individuals are hybridized between the populations and the fitness of the hybrids is measured (F). The hybrid individuals are then discarded and are not used for further breeding. Over time, th.