Acknowledgments RP would like to thank Daniel Weissman, Amanda Casto and Robert Furrow for helpful discussions throughout the course of this work, Jeremy Van Cleve and Michael E. Palmer for technical support, and Isabel Gordo and Tiago Paixa˜o for helpful comments and insights. Special thanks to Manuel Irimia for all the guidance and help. References 1. Bornholdt S (2005) Systems biology. Less is more in modeling large genetic networks. Science 310: 449–451. 2. Bornholdt S (2008) Boolean network models of cellular regulation: prospects and limitations. J R Soc Interface 5: S85–S94. 3. Ku¨ rten KE (1988) Correspondence between neural threshold networks and Kauffman Boolean cellular automata. J Phys A: Math Gen 21: L615–L619. 4. Hopfield JJ (1982) Neural networks and physical systems with emergent collective computational abilities. Proc Natl Acad Sci U S A 79: 2554–2558. 5. Kauffman SA (1969) Metabolic stability and epigenesis in randomly constructed genetic nets. J Theor Biol 22: 437–467. 6. Aldana M, Coppersmith S, Kadanoff LP (2003) Boolean dynamics with random couplings. In: Kaplan E, Marsden JE, Sreenivasan KR, eds. Perspectives and Problems in Nonlinear Science. A Celebratory Volume in Honour of Lawerence Sirovid. Berlin: Springer. pp 23–89. 7. Drossel B (2008) Random Boolean Networks. arXiv:07063351v2 [cond-matstat- mech]. 8. Wagner A (1996) Does evolutionary plasticity evolve? Evolution 50: 1008–1023. 9. Siegal ML, Bergman A (2002) Waddington’s canalization revisited: develop- mental stability and evolution. Proc Natl Acad Sci U S A 99: 10528–10532. 10. Wagner A (1994) Evolution of gene networks by gene duplications: a mathematical model and its implications on genome organization. Proc Natl Acad Sci U S A 91: 4387–4391. 11. Ciliberti S, Martin OC, Wagner A (2007) Robustness can evolve gradually in complex regulatory gene networks with varying topology. PLoS Comput Biol 3(2): e15. doi:10.1371/journal.pcbi.0030015. 12. Huerta-Sanchez E, Durrett R (2007) Wagner’s canalization model. Theor Popul Biol 71: 121–130. 13. Draghi JA, Wagner GP (2009) The evolutionary dynamics of evolvability in a gene network model. J Evol Biol 22: 599–611. 14. Bornholdt S (2001) Modeling genetic networks and their evolution: a complex dynamical systems perspective. Biol Chem 382: 1289–1299. 15. Kaneko K (2007) Evolution of robustness to noise and mutation in gene expression dynamics. PLoS One 2: e434. 16. Azevedo RBR, Lohaus R, Srinivasan S, Dang KK, Burch CL (2006) Sexual reproduction selects for robustness and negative epistasis in artificial gene networks. Nature 440: 87–90. 17. Palmer ME, Feldman MW (2009) Dynamics of hybrid incompatibility in gene networks in a constant environment. Evolution 63: 418–431. 18. Li F, Long T, Lu Y, Ouyang Q, Tang C (2004) The yeast cell-cycle network is robustly designed. Proc Natl Acad Sci U S A 101: 4781–4786. 19. Masel J (2004) Genetic assimilation can occur in the absence of selection for the assimilating phenotype, suggesting a role for the canalization heuristic. J Evol Biol 17: 1106–1110. 20. McDonald D, Waterbury L, Knight R, Betterton MD (2008) Activating and inhibiting connections in biological network dynamics. Biol Direct 3: 49. 21. Ciliberti S, Martin OC, Wagner A (2007) Innovation and robustness in complex regulatory gene networks. Proc Natl Acad Sci U S A 104: 13591–13596. 22. Bergman A, Siegal ML (2003) Evolutionary capacitance as a general feature of complex gene networks. Nature 424: 604–607. 23. Drossel B, Mihaljev T, Greil F (2005) Number and length of attractors in a critical Kauffman model with connectivity one. Phys Rev Lett 94: 088701. 24. Rohlf T, Gulbahce N, Teuscher C (2007) Damage spreading and criticality in finite random dynamical networks. Phys Rev Lett 99: 248701. 25. Greil F, Drossel B (2005) Dynamics of critical Kauffman networks under asynchronous stochastic update. Phys Rev Lett 95: 3–6. 26. Greil F, Drossel B, Sattler J (2007) Critical Kauffman networks under deterministic asynchronous update. New J Phys 9: 373–373. 27. Klemm K, Bornholdt S (2005) Stable and unstable attractors in Boolean networks. Phys Rev E 72: 1–4. 28. Albert R, Othmer HG (2003) The topology of the regulatory interactions predicts the expression pattern of the segment polarity genes in Drosophila melanogaster. J Theor Biol 223: 1–18. Most Networks in Wagner’s Model Are Cycling Author Contributions Conceived and designed the experiments: RP EB MF. Performed the experiments: RP. Analyzed the data: RP. Wrote the paper: RP MF. 29. Espinosa-Soto C, Padilla-Longoria P, Alvarez-Buylla ER (2004) A gene regulatory network model for cell-fate determination during Arabidopsis thaliana flower development that is robust and recovers experimental gene expression profiles. Plant Cell 16: 2923–2939. 30. Borenstein E, Krakauer DC (2008) An end to endless forms: epistasis, phenotype distribution bias, and nonuniform evolution. PLoS Comput Biol 4: e1000202. 31. Sevim V, Rikvold PA (2008) Chaotic gene regulatory networks can be robust against mutations and noise. J Theor Biol 253: 323–332. 32. Leclerc RD (2008) Survival of the sparsest: robust gene networks are parsimonious. Mol Syst Biol 4: 213. 33. Martin OC, Wagner A (2009) Effects of recombination on complex regulatory circuits. Genetics 183: 673–684. 34. Siegal ML, Promislow DEL, Bergman A (2007) Functional and evolutionary inference in gene networks: does topology matter? Genetica 129: 83–103. 35. Rodriguez-Caso C, Corominas-Murtraa B, Sole´ RV (2009) On the basic computational structure of gene regulatory networks. Mol Biosyst 5: 1617–1629. 36. Rohlf T, Winkler CR (2009) Emergent network structure, evolvable robustness and non-linear effects of point mutations in an artificial genome model. arXiv:09083610v1 [q-bioMN]. 37. Burda Z, Zagorski M, Krzywicki A, Martin OC (2009) Sparse essential interactions in model networks of gene regulation. arXiv:09104077v1 [q- bioMN]. 38. Espinosa-Soto C, Wagner A (2010) Specialization can drive the evolution of modularity. PLoS Comput Biol 6: e1000719. 39. Kauffman SA, Peterson C, Samuelsson B, Troein C (2004) Genetic networks with canalyzing Boolean rules are always stable. Proc Natl Acad Sci U S A 101: 17102–17107. 40. Paul U, Kaufman V, Drossel B (2006) Properties of attractors of canalyzing random Boolean networks. Phys Rev E 73: 1–9. 41. Greil F, Drossel B (2007) Kauffman networks with threshold functions. Eur Phys J B 57: 109–113. 42. Ku¨ rten KE (1988) Critical phenomena in model neural networks. Phys Lett A 129: 157–160. 43. Rohlf T, Bornholdt S (2002) Criticality in random threshold networks: annealed approximation and beyond. Physica A 310: 245–259. 44. Anafi RC, Bates JHT (2010) Balancing robustness against the dangers of multiple attractors in a Hopfield-type model of biological attractors. PLoS One 5: e14413. 45. Gama-Castro S, Salgado H, Peralta-Gil M, Santos-Zavaleta A, Mun˜ iz Rascado L, et al. (2011) RegulonDB version 7.0: transcriptional regulation of Escherichia coli K-12 integrated within genetic sensory response units (Gensor Units). Nucleic Acids Res 39: D98–105. 46. Blattner FR, Plunkett G, 3rd, Bloch CA, Perna NT, Burland V, et al. (1997) The Complete Genome Sequence of Escherichia coli K-12. Science 277: 1453–1462. 47. Luscombe NM, Babu MM, Yu H, Snyder M, Teichmann SA, et al. (2004) Genomic analysis of regulatory network dynamics reveals large topological changes. Nature 431: 308–312. 48. Cherry JM, Ball C, Weng S, Juvik G, Schmidt R, et al. (1997) Genetic and physical maps of Saccharomyces cerevisiae. Nature 387: 67–73. 49. Thieffry D, Huerta AM, Pe E, Collado-Vides J (1998) From specific gene regulation to genomic networks: a global analysis of transcriptional regulation in Escherichia coli. BioEssays 20: 433–440. 50. Jeong H, Tombor B, Albert R, Oltvai ZN, Barabasi AL (2000) The large-scale organization of metabolic networks. Nature 407: 651–654. 51. Samuelsson B, Troein C (2003) Superpolynomial growth in the number of attractors in Kauffman networks. Phys Rev Lett 90: 098701. 52. Kauffman SA, Peterson C, Samuelsson B, Troein C (2003) Random Boolean network models and the yeast transcriptional network. Proc Natl Acad Sci U S A 100: 14796–14799. 53. Guelzim N, Bottani S, Bourgine P, Ke´pe`s F (2002) Topological and causal structure of the yeast transcriptional regulatory network. Nat Genet 31: 60–63. PLoS ONE | www.plosone.org 8 April 2012 | Volume 7 | Issue 4 | e34285