{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,3,4]],"date-time":"2026-03-04T07:04:25Z","timestamp":1772607865215,"version":"3.50.1"},"update-to":[{"DOI":"10.1371\/journal.pcbi.1009796","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2022,1,31]],"date-time":"2022-01-31T00:00:00Z","timestamp":1643587200000}}],"reference-count":65,"publisher":"Public Library of Science (PLoS)","issue":"1","license":[{"start":{"date-parts":[[2022,1,19]],"date-time":"2022-01-19T00:00:00Z","timestamp":1642550400000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["www.ploscompbiol.org"],"crossmark-restriction":false},"short-container-title":["PLoS Comput Biol"],"abstract":"<jats:p>The aim of this paper is two-fold. First, we propose a new computational method to investigate the particularities of evolution. Second, we apply this method to a model of gene regulatory networks (GRNs) and explore the evolution of mutational robustness and bistability. Living systems have developed their functions through evolutionary processes. To understand the particularities of this process theoretically, evolutionary simulation (ES) alone is insufficient because the outcomes of ES depend on evolutionary pathways. We need a reference system for comparison. An appropriate reference system for this purpose is an ensemble of the randomly sampled genotypes. However, generating high-fitness genotypes by simple random sampling is difficult because such genotypes are rare. In this study, we used the multicanonical Monte Carlo method developed in statistical physics to construct a reference ensemble of GRNs and compared it with the outcomes of ES. We obtained the following results. First, mutational robustness was significantly higher in ES than in the reference ensemble at the same fitness level. Second, the emergence of a new phenotype, bistability, was delayed in evolution. Third, the bistable group of GRNs contains many mutationally fragile GRNs compared with those in the non-bistable group. This suggests that the delayed emergence of bistability is a consequence of the mutation-selection mechanism.<\/jats:p>","DOI":"10.1371\/journal.pcbi.1009796","type":"journal-article","created":{"date-parts":[[2022,1,19]],"date-time":"2022-01-19T18:36:39Z","timestamp":1642617399000},"page":"e1009796","update-policy":"https:\/\/doi.org\/10.1371\/journal.pcbi.corrections_policy","source":"Crossref","is-referenced-by-count":14,"title":["Evolution enhances mutational robustness and suppresses the emergence of a new phenotype: A new computational approach for studying evolution"],"prefix":"10.1371","volume":"18","author":[{"ORCID":"https:\/\/orcid.org\/0000-0003-2375-5353","authenticated-orcid":true,"given":"Tadamune","family":"Kaneko","sequence":"first","affiliation":[]},{"ORCID":"https:\/\/orcid.org\/0000-0003-1582-8970","authenticated-orcid":true,"given":"Macoto","family":"Kikuchi","sequence":"additional","affiliation":[]}],"member":"340","published-online":{"date-parts":[[2022,1,19]]},"reference":[{"key":"pcbi.1009796.ref001","doi-asserted-by":"crossref","first-page":"828","DOI":"10.1038\/nrg1471","article-title":"Biological robustness","volume":"5","author":"H Kitano","year":"2004","journal-title":"Nat Rev Genet"},{"key":"pcbi.1009796.ref002","volume-title":"Robustness and Evolvability in Living Systems","author":"A Wagner","year":"2005"},{"key":"pcbi.1009796.ref003","doi-asserted-by":"crossref","first-page":"395","DOI":"10.1016\/j.tig.2009.07.005","article-title":"Robustness: mechanisms and consequences","volume":"25","author":"J Masel","year":"2009","journal-title":"Trends Genet"},{"key":"pcbi.1009796.ref004","doi-asserted-by":"crossref","first-page":"387","DOI":"10.1038\/nature00935","article-title":"Functional profiling of the Saccharomyces cerevisiae genome","volume":"418","author":"G Giaever","year":"2002","journal-title":"Nature"},{"key":"pcbi.1009796.ref005","doi-asserted-by":"crossref","DOI":"10.1038\/msb4100050","article-title":"Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection","volume":"2","author":"T Baba","year":"2006","journal-title":"Mol Syst Biol"},{"key":"pcbi.1009796.ref006","doi-asserted-by":"crossref","first-page":"617","DOI":"10.1038\/nbt.1628","article-title":"Analysis of a genome-wide set of gene deletions in the fission yeast Schizosaccharomyces pombe","volume":"28","author":"DU Kim","year":"2010","journal-title":"Nat Biotechnol"},{"key":"pcbi.1009796.ref007","doi-asserted-by":"crossref","first-page":"840","DOI":"10.1038\/nature06847","article-title":"Evolvability and hierarchy in rewired bacterial gene networks","volume":"452","author":"M Isaran","year":"2008","journal-title":"Nature"},{"key":"pcbi.1009796.ref008","doi-asserted-by":"crossref","first-page":"e15","DOI":"10.1371\/journal.pcbi.0030015","article-title":"Robustness can evolve gradually in complex regulatory gene networks with varying topology","volume":"3","author":"S Ciliberti","year":"2007","journal-title":"PLoS Comput Biol"},{"key":"pcbi.1009796.ref009","doi-asserted-by":"crossref","first-page":"13591","DOI":"10.1073\/pnas.0705396104","article-title":"Innovation and robustness in complex regulatory gene networks","volume":"104","author":"S Ciliberti","year":"2007","journal-title":"Proc Natl Acad Sci USA"},{"key":"pcbi.1009796.ref010","doi-asserted-by":"crossref","first-page":"279","DOI":"10.1098\/rspb.1994.0040","article-title":"From sequences to shapes and back: a case study in RNA secondary structures","volume":"255","author":"P Schuster","year":"1994","journal-title":"Proc R Soc Lond B"},{"key":"pcbi.1009796.ref011","doi-asserted-by":"crossref","first-page":"17263","DOI":"10.1073\/pnas.1109435108","article-title":"Motifs emerge from function in model gene regulatory networks","volume":"108","author":"Z Burda","year":"2011","journal-title":"Proc Natl Acad Sci USA"},{"key":"pcbi.1009796.ref012","doi-asserted-by":"crossref","first-page":"012727","DOI":"10.1103\/PhysRevE.87.012727","article-title":"Edge usage, motifs, and regulatory logic for cell cycling genetic networks","volume":"87","author":"M Zagorski","year":"2013","journal-title":"Phys Rev E"},{"key":"pcbi.1009796.ref013","doi-asserted-by":"crossref","first-page":"053037","DOI":"10.1088\/1367-2630\/15\/5\/053037","article-title":"Robustness leads close to the edge of chaos in coupled map networks: toward the understanding of biological networks","volume":"15","author":"N Saito","year":"2013","journal-title":"New J Phys"},{"key":"pcbi.1009796.ref014","doi-asserted-by":"crossref","first-page":"249","DOI":"10.1016\/0370-2693(91)91256-U","article-title":"Multicanonical algorithms for first order phase transitions","volume":"B267","author":"BA Berg","year":"1991","journal-title":"Phys Lett"},{"key":"pcbi.1009796.ref015","doi-asserted-by":"crossref","first-page":"9","DOI":"10.1103\/PhysRevLett.68.9","article-title":"Multicanonical ensemble: A new approach to simulate first-order phase transitions","volume":"68","author":"BA Berg","year":"1992","journal-title":"Phys Rev Lett"},{"key":"pcbi.1009796.ref016","doi-asserted-by":"crossref","first-page":"611","DOI":"10.1007\/s10463-014-0460-2","article-title":"Multicanonical MCMC for sampling rare events: an illustrative review","volume":"66","author":"Y Iba","year":"2014","journal-title":"Ann Inst Stat Math"},{"key":"pcbi.1009796.ref017","doi-asserted-by":"crossref","first-page":"031142","DOI":"10.1103\/PhysRevE.82.031142","article-title":"Multicanonical sampling of rare events in random matrices","volume":"82","author":"N Saito","year":"2010","journal-title":"Phys Rev E"},{"key":"pcbi.1009796.ref018","doi-asserted-by":"crossref","first-page":"223","DOI":"10.1016\/j.cpc.2010.06.039","article-title":"Probability of graphs with large spectral gap by multicanonical Monte Carlo","volume":"182","author":"N Saito","year":"2011","journal-title":"Comput Phys Commun"},{"key":"pcbi.1009796.ref019","doi-asserted-by":"crossref","first-page":"251","DOI":"10.1016\/j.cpc.2010.07.044","article-title":"Multicanonical sampling of rare trajectories in chaotic dynamical systems","volume":"182","author":"A Kitajima","year":"2011","journal-title":"Comput Phys Commun"},{"key":"pcbi.1009796.ref020","doi-asserted-by":"crossref","first-page":"e0 125062","DOI":"10.1371\/journal.pone.0125062","article-title":"Numerous but rare: An exploration of magic squares","volume":"10","author":"A Kitajima","year":"2015","journal-title":"PLoS ONE"},{"key":"pcbi.1009796.ref021","doi-asserted-by":"crossref","first-page":"e1007969","DOI":"10.1371\/journal.pcbi.1007969","article-title":"Emergence of cooperative bistability and robustness of gene regulatory networks","volume":"16","author":"S Nagata","year":"2020","journal-title":"PLoS Comput Biol"},{"key":"pcbi.1009796.ref022","doi-asserted-by":"crossref","first-page":"6840","DOI":"10.1073\/pnas.78.11.6840","article-title":"An amplified sensitivity arising from covalent modification in biological systems","volume":"78","author":"A Goldbeter","year":"1981","journal-title":"Proc Natl Acad Sci USA"},{"key":"pcbi.1009796.ref023","doi-asserted-by":"crossref","first-page":"496","DOI":"10.1016\/j.tibs.2014.08.003","article-title":"Ultrasensitivity part I: Michaelian responses and zero-order ultrasensitivity","volume":"39","author":"JE Ferrell","year":"2014","journal-title":"Trends Biochem Sci"},{"key":"pcbi.1009796.ref024","volume-title":"An Introduction to Systems Biology","author":"U Alon","year":"2020"},{"key":"pcbi.1009796.ref025","volume-title":"Genes & Signals","author":"M Ptashne","year":"2002"},{"key":"pcbi.1009796.ref026","volume-title":"A Genetic Switch","author":"M Ptashne","year":"2004"},{"key":"pcbi.1009796.ref027","doi-asserted-by":"crossref","first-page":"440","DOI":"10.1038\/msb.2010.96","article-title":"Lysogen stability is determined by the frequency of activity bursts from the fate-determining gene","volume":"6","author":"C Zong","year":"2010","journal-title":"Mol Syst Biol"},{"key":"pcbi.1009796.ref028","volume-title":"Molecular Biology of the Gene","author":"JD Watson","year":"2013","edition":"7"},{"key":"pcbi.1009796.ref029","doi-asserted-by":"crossref","first-page":"346","DOI":"10.1038\/ncb954","article-title":"Building a cell cycle oscillator: hysteresis and bistability in the activation of Cdc2","volume":"5","author":"JR Pomerening","year":"2003","journal-title":"Nat Cell Biol"},{"key":"pcbi.1009796.ref030","doi-asserted-by":"crossref","first-page":"975","DOI":"10.1073\/pnas.0235349100","article-title":"Hysteresis drives cell-cycle transitions in Xenopus laevis egg extracts","volume":"100","author":"W Sha","year":"2003","journal-title":"Proc Natl Acad Sci USA"},{"key":"pcbi.1009796.ref031","doi-asserted-by":"crossref","first-page":"263","DOI":"10.1016\/j.molcel.2011.01.012","article-title":"Ultrasensitivity in the regulation of Cdc25C by Cdk1","volume":"41","author":"NB Trunnell","year":"2011","journal-title":"Mol Cell"},{"key":"pcbi.1009796.ref032","doi-asserted-by":"crossref","first-page":"895","DOI":"10.1126\/science.280.5365.895","article-title":"The biochemical bases of an all-or-none cell fate switch in Xenopus oocytes","volume":"280","author":"JE Ferrell","year":"1998","journal-title":"Science"},{"key":"pcbi.1009796.ref033","doi-asserted-by":"crossref","first-page":"38002","DOI":"10.1209\/0295-5075\/124\/38002","article-title":"Cooperative reliable response from sloppy gene-expression dynamics","volume":"124","author":"M Inoue","year":"2018","journal-title":"Eutophys Lett"},{"key":"pcbi.1009796.ref034","doi-asserted-by":"crossref","first-page":"437","DOI":"10.1016\/0022-5193(69)90015-0","article-title":"Metabolic stability and epigenesis in randomly constructed genetic nets","volume":"22","author":"SA Kauffman","year":"1969","journal-title":"J Theor Biol"},{"key":"pcbi.1009796.ref035","doi-asserted-by":"crossref","first-page":"11","DOI":"10.1016\/S0022-5193(87)80029-2","article-title":"Towards a general theory of adaptive walks on rugged landscapes","volume":"128","author":"SA Kauffman","year":"1987","journal-title":"J Theor Biol"},{"key":"pcbi.1009796.ref036","doi-asserted-by":"crossref","DOI":"10.1093\/oso\/9780195079517.001.0001","volume-title":"The Origins of Order: Self-Organization and Selection in Evolution","author":"SA Kauffman","year":"1993"},{"key":"pcbi.1009796.ref037","doi-asserted-by":"crossref","first-page":"429","DOI":"10.1016\/S0022-5193(05)80391-1","article-title":"A connectionist model of development","volume":"152","author":"E Mjolsness","year":"1991","journal-title":"J Theor Biol"},{"key":"pcbi.1009796.ref038","doi-asserted-by":"crossref","first-page":"4387","DOI":"10.1073\/pnas.91.10.4387","article-title":"Evolution of gene networks by gene duplications: A mathematical model and its implications on genome organization","volume":"91","author":"A Wagner","year":"1994","journal-title":"Proc Natl Acad Sci USA"},{"key":"pcbi.1009796.ref039","doi-asserted-by":"crossref","first-page":"1008","DOI":"10.1111\/j.1558-5646.1996.tb02342.x","article-title":"Does evolutionary plasticity evolve","volume":"50","author":"A Wagner","year":"1996","journal-title":"Evolution"},{"key":"pcbi.1009796.ref040","doi-asserted-by":"crossref","first-page":"10528","DOI":"10.1073\/pnas.102303999","article-title":"Waddington\u2019s canalization revisited: developmental stability and evolution","volume":"99","author":"ML Siegal","year":"2002","journal-title":"Proc Natl Acad Sci USA"},{"key":"pcbi.1009796.ref041","doi-asserted-by":"crossref","first-page":"1106","DOI":"10.1111\/j.1420-9101.2004.00739.x","article-title":"Genetic assimilation can occur in the absence of selection for the assimilating phenotype, suggesting a role for the canalization heuristic","volume":"17","author":"J Masel","year":"2004","journal-title":"J Evol Biol"},{"key":"pcbi.1009796.ref042","doi-asserted-by":"crossref","first-page":"e434","DOI":"10.1371\/journal.pone.0000434","article-title":"Evolution of robustness to noise and mutation in gene expression dynamics","volume":"2","author":"K Kaneko","year":"2007","journal-title":"PLoS one"},{"key":"pcbi.1009796.ref043","doi-asserted-by":"crossref","first-page":"1284","DOI":"10.1111\/j.1420-9101.2011.02261.x","article-title":"Phenotypic robustness can increase phenotypic variability after nongenetic perturbations in gene regulatory circuits","volume":"24","author":"C Espinosa-Soto","year":"2011","journal-title":"J Evol Biol"},{"key":"pcbi.1009796.ref044","doi-asserted-by":"crossref","first-page":"e1003001","DOI":"10.1371\/journal.pcbi.1003001","article-title":"Cooperative adaptive responses in gene regulatory networks with many degrees of freedom","volume":"9","author":"M Inoue","year":"2013","journal-title":"PLoS Comput Biol"},{"key":"pcbi.1009796.ref045","doi-asserted-by":"crossref","first-page":"e3","DOI":"10.1371\/journal.pcbi.0040003","article-title":"A generic mechanism for adaptive growth rate regulation","volume":"4","author":"C Furusawa","year":"2008","journal-title":"PLoS Comput Biol"},{"key":"pcbi.1009796.ref046","doi-asserted-by":"crossref","first-page":"e1002865","DOI":"10.1371\/journal.pcbi.1002865","article-title":"The underlying molecular and network level mechanisms in the evolution of robustness in gene regulatory networks","volume":"9","author":"M Pujato","year":"2013","journal-title":"PLoS Comput Biol"},{"key":"pcbi.1009796.ref047","doi-asserted-by":"crossref","first-page":"141","DOI":"10.1007\/BF00339943","article-title":"\u201cNeural\u201d Computation of Decisionsin Optimization Problems","volume":"52","author":"JJ Hopfield","year":"1985","journal-title":"Biol Cybern"},{"key":"pcbi.1009796.ref048","doi-asserted-by":"crossref","first-page":"221","DOI":"10.1016\/S0955-0674(03)00017-6","article-title":"Sniffers, buzzers, toggles and blinkers: dynamics of regulatory and signaling pathways in the cell","volume":"15","author":"JJ Tyson","year":"2003","journal-title":"Curr Opin Cell Biol"},{"key":"pcbi.1009796.ref049","doi-asserted-by":"crossref","first-page":"64","DOI":"10.1038\/ng881","article-title":"Network motifs in the transcriptional regulation network of Escherichia coli","volume":"31","author":"SS Shen-Orr","year":"2002","journal-title":"Nat Genet"},{"key":"pcbi.1009796.ref050","doi-asserted-by":"crossref","first-page":"9716","DOI":"10.1073\/pnas.96.17.9716","article-title":"Neutral evolution of mutational robustness","volume":"17","author":"E van Nimwegen","year":"1999","journal-title":"Proc Natl Acad Sci USA"},{"key":"pcbi.1009796.ref051","doi-asserted-by":"crossref","first-page":"2321","DOI":"10.1111\/jeb.12959","article-title":"Selection for distinct gene expression properties favours the evolution of mutational robustness in gene regulatory networks","volume":"29","author":"E Espinosa-Soto","year":"2016","journal-title":"J Evol Biol"},{"key":"pcbi.1009796.ref052","doi-asserted-by":"crossref","first-page":"11980","DOI":"10.1073\/pnas.2133841100","article-title":"Struncure and function of the feed-forward loop network motif","volume":"100","author":"S Mangan","year":"2003","journal-title":"Proc Natl Acad Sci USA"},{"key":"pcbi.1009796.ref053","doi-asserted-by":"crossref","first-page":"227","DOI":"10.1063\/1.1349894","article-title":"Bistability in cell signaling: How to make continuous processes discontinuous, and reversible processes irreversible","volume":"11","author":"JE Ferrell","year":"2001","journal-title":"Chaos"},{"key":"pcbi.1009796.ref054","doi-asserted-by":"crossref","first-page":"140","DOI":"10.1016\/S0955-0674(02)00314-9","article-title":"Self-perpetuating states in signal transduction: positive feedback, doule-negative feedback and bistability","volume":"6","author":"JE Ferrell","year":"2002","journal-title":"Curr Open Chem Biol"},{"key":"pcbi.1009796.ref055","doi-asserted-by":"crossref","first-page":"1822","DOI":"10.1073\/pnas.0308265100","article-title":"Detection of multistability, bifurcations, and hysteresis in a large class of biological positive-feedback systems","volume":"101","author":"D Angeli","year":"2003","journal-title":"Proc Natl Acad Sci USA"},{"key":"pcbi.1009796.ref056","doi-asserted-by":"crossref","first-page":"046013","DOI":"10.1088\/1478-3975\/6\/4\/046013","article-title":"The combination of positive and negative feedback loops confers flexibility to biochemical switches","volume":"6","author":"B Pfeuty","year":"2009","journal-title":"Phys Biol"},{"key":"pcbi.1009796.ref057","doi-asserted-by":"crossref","first-page":"055003","DOI":"10.1088\/1478-3975\/9\/5\/055003","article-title":"Coupling between feedbak loops in autoregulatory networks affects bistability range, open-loop gain and switching times","volume":"9","author":"A Tiwari","year":"2012","journal-title":"Phys Biol"},{"key":"pcbi.1009796.ref058","doi-asserted-by":"crossref","first-page":"612","DOI":"10.1016\/j.tibs.2014.10.002","article-title":"Ultrasensitivity part III: cascades, bistable switches, and oscillators","volume":"39","author":"JE Ferrell","year":"2014","journal-title":"Trends Biochem Sci"},{"key":"pcbi.1009796.ref059","doi-asserted-by":"crossref","first-page":"695","DOI":"10.1016\/j.ydbio.2007.02.036","article-title":"Bifurcation dynamics in lineage-commitment in bipotent progenitor cells","volume":"305","author":"S Huang","year":"2007","journal-title":"Dev Biol"},{"key":"pcbi.1009796.ref060","doi-asserted-by":"crossref","first-page":"035007","DOI":"10.1088\/1478-3975\/aa6f90","article-title":"Operating principles of tristable circuits regulating cellular differentiation","volume":"14","author":"D Jia","year":"2017","journal-title":"Phys Biol"},{"key":"pcbi.1009796.ref061","doi-asserted-by":"crossref","first-page":"e1005456","DOI":"10.1371\/journal.pcbi.1005456","article-title":"Interrogating the topological robustness of gene regulatory circuits by randomization","volume":"13","author":"B Huang","year":"2017","journal-title":"PLoS Comput Biol"},{"key":"pcbi.1009796.ref062","doi-asserted-by":"crossref","first-page":"14086","DOI":"10.1073\/pnas.2334996100","article-title":"On the relation between fluctuation and response in biological systems","volume":"100","author":"K Sato","year":"2003","journal-title":"Proc Natl Acad Sci USA"},{"key":"pcbi.1009796.ref063","doi-asserted-by":"crossref","first-page":"211","DOI":"10.1103\/PhysRevLett.71.211","article-title":"New Monte Carlo Algorithm: Entropic Sampling","volume":"71","author":"J Lee","year":"1993","journal-title":"Phys Rev Lett"},{"key":"pcbi.1009796.ref064","doi-asserted-by":"crossref","first-page":"2050","DOI":"10.1103\/PhysRevLett.86.2050","article-title":"Efficient, Multiple-Range Random Walk Algorithm to Calculate the Density of State","volume":"86","author":"F Wang","year":"2001","journal-title":"Phys Rev Lett"},{"key":"pcbi.1009796.ref065","doi-asserted-by":"crossref","first-page":"056101","DOI":"10.1103\/PhysRevE.64.056101","article-title":"Determining the density of states for classical statistical models: A random walk algorithm to produce a flat histogram","volume":"64","author":"F Wang","year":"2001","journal-title":"Phys Rev E"}],"updated-by":[{"DOI":"10.1371\/journal.pcbi.1009796","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2022,1,31]],"date-time":"2022-01-31T00:00:00Z","timestamp":1643587200000}}],"container-title":["PLOS Computational Biology"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1009796","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2023,11,16]],"date-time":"2023-11-16T00:29:52Z","timestamp":1700094592000},"score":1,"resource":{"primary":{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1009796"}},"subtitle":[],"editor":[{"given":"Alexandre V.","family":"Morozov","sequence":"first","affiliation":[]}],"short-title":[],"issued":{"date-parts":[[2022,1,19]]},"references-count":65,"journal-issue":{"issue":"1","published-online":{"date-parts":[[2022,1,19]]}},"URL":"https:\/\/doi.org\/10.1371\/journal.pcbi.1009796","relation":{"new_version":[{"id-type":"doi","id":"10.1371\/journal.pcbi.1009796","asserted-by":"object"}]},"ISSN":["1553-7358"],"issn-type":[{"value":"1553-7358","type":"electronic"}],"subject":[],"published":{"date-parts":[[2022,1,19]]}}}