Eukaryotic cells are dynamically ordered or critical but not chaotic.

Two important theoretical approaches have been developed to generically characterize the relationship between the structure and function of large genetic networks: the continuous approach, based on reaction-kinetics differential equations, and the Boolean approach, based on difference equations and discrete logical rules. These two approaches do not always coincide in their predictions for the same system. Nonetheless, both of them predict that the highly nonlinear dynamics exhibited by genetic regulatory systems can be characterized into two broad regimes, to wit, an ordered regime where the system is robust against perturbations, and a chaotic regime where the system is extremely sensitive to perturbations. It has been a plausible and long-standing hypothesis that genomic regulatory networks of real cells operate in the ordered regime or at the border between order and chaos. This hypothesis is indirectly supported by the robustness and stability observed in the phenotypic traits of living organisms under genetic perturbations. However, there has been no systematic study to determine whether the gene-expression patterns of real cells are compatible with the dynamically ordered regimes predicted by theoretical models. Using the Boolean approach, here we show what we believe to be the first direct evidence that the underlying genetic network of HeLa cells appears to operate either in the ordered regime or at the border between order and chaos but does not appear to be chaotic.

[1]  A. L. Ruhkin Testing Randomness: A Suite of Statistical Procedures , 2001 .

[2]  S. Kauffman Metabolic stability and epigenesis in randomly constructed genetic nets. , 1969, Journal of theoretical biology.

[3]  Gabriel S. Eichler,et al.  Cell fates as high-dimensional attractor states of a complex gene regulatory network. , 2005, Physical review letters.

[4]  D. Murray,et al.  A genomewide oscillation in transcription gates DNA replication and cell cycle. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[5]  P. Cluzel,et al.  A natural class of robust networks , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[6]  C. Espinosa-Soto,et al.  A Gene Regulatory Network Model for Cell-Fate Determination during Arabidopsis thaliana Flower Development That Is Robust and Recovers Experimental Gene Expression Profilesw⃞ , 2004, The Plant Cell Online.

[7]  Jerrold E. Marsden,et al.  Perspectives and Problems in Nonlinear Science , 2003 .

[8]  Denis Thieffry,et al.  Genetic control of flower morphogenesis in Arabidopsis thaliana: a logical analysis , 1999, Bioinform..

[9]  Ilya Shmulevich,et al.  Binary analysis and optimization-based normalization of gene expression data , 2002, Bioinform..

[10]  H. Othmer,et al.  The topology of the regulatory interactions predicts the expression pattern of the segment polarity genes in Drosophila melanogaster. , 2003, Journal of theoretical biology.

[11]  Jaakko Astola,et al.  Microarray quality control , 2004 .

[12]  Trevor Hastie,et al.  The Elements of Statistical Learning , 2001 .

[13]  L. Glass,et al.  Ordered and disordered dynamics in random networks , 1998 .

[14]  M. Aldana Boolean dynamics of networks with scale-free topology , 2003 .

[15]  S. Kauffman,et al.  Activities and sensitivities in boolean network models. , 2004, Physical review letters.

[16]  Carsten Peterson,et al.  Random Boolean network models and the yeast transcriptional network , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[17]  Abraham Lempel,et al.  On the Complexity of Finite Sequences , 1976, IEEE Trans. Inf. Theory.

[18]  E. Winzeler,et al.  Genomics, gene expression and DNA arrays , 2000, Nature.

[19]  Michael A. Savageau,et al.  Effects of alternative connectivity on behavior of randomly constructed Boolean networks , 2002 .

[20]  E. Álvarez-Buylla,et al.  Dynamics of the genetic regulatory network for Arabidopsis thaliana flower morphogenesis. , 1998, Journal of theoretical biology.

[21]  Stuart A. Kauffman,et al.  The origins of order , 1993 .

[22]  Jason H. Moore,et al.  A global view of epistasis , 2005, Nature Genetics.

[23]  Vladimir D. Gusev,et al.  On the complexity measures of genetic sequences , 1999, Bioinform..

[24]  P. Lio’,et al.  Periodic gene expression program of the fission yeast cell cycle , 2004, Nature Genetics.

[25]  C. Ball,et al.  Identification of genes periodically expressed in the human cell cycle and their expression in tumors. , 2002, Molecular biology of the cell.

[26]  Nicola J. Rinaldi,et al.  Transcriptional Regulatory Networks in Saccharomyces cerevisiae , 2002, Science.

[27]  A. Wagner Distributed robustness versus redundancy as causes of mutational robustness. , 2005, BioEssays : news and reviews in molecular, cellular and developmental biology.

[28]  A. Zeng,et al.  An extended transcriptional regulatory network of Escherichia coli and analysis of its hierarchical structure and network motifs. , 2004, Nucleic acids research.

[29]  S. Huang,et al.  Shape-dependent control of cell growth, differentiation, and apoptosis: switching between attractors in cell regulatory networks. , 2000, Experimental cell research.

[30]  J. Collado-Vides,et al.  Identifying global regulators in transcriptional regulatory networks in bacteria. , 2003, Current opinion in microbiology.

[31]  M Villani,et al.  Genetic network models and statistical properties of gene expression data in knock-out experiments. , 2004, Journal of theoretical biology.

[32]  Q. Ouyang,et al.  The yeast cell-cycle network is robustly designed. , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[33]  L. Glass,et al.  Common Chaos in Arbitrarily Complex Feedback Networks , 1997 .

[34]  Jaakko Astola,et al.  The role of certain Post classes in Boolean network models of genetic networks , 2003, Proceedings of the National Academy of Sciences of the United States of America.