Abstract
This paper presents a modified random network model to illustrate how groups can form in the absence of evolutionary forces, assuming groups are collections of entities at any level of organization. This model is inspired by the Zero Force Evolutionary Law, which states that there is always a tendency for diversity and complexity to increase in any evolutionary system containing variation and heredity. That is, in the absence of evolutionary forces, the expectation is a continual increase in diversity and complexity at any level of biological hierarchy. I show that, when modeled, this expectation of increasing variation results not only in the formation of groups, but also in a higher probability of group formation than is found in a model that is purely random.
Similar content being viewed by others
Notes
This is often referred to in the graph theory and network literature as a clustering coefficient, and random networks are known for their extremely low clustering coefficients.
This example is also useful because it illustrates the point that there can be degrees of connectivity between entities. This is not apparent in the ZFEL model since it is designed such that entities are either connected or they are not. Weighted networks incorporate weighted edges, or different edge strengths, and I hope to address this issue in future research.
This denominator was also used when measuring 3-cycles and 4-cliques in the random model, except for changing the number of edges being chosen (see “Appendix 1” for details).
This denominator was also used for measuring 3-cycles and 4-cliques in the ZFEL model, except for changing the number of edges being chosen (see “Appendix 1” for details).
It is not necessary for each new node state to be unique; this method was just for ease of calculation and evaluation. A different model, where at time T 1 differentiated nodes change from state 0 to 1, and then at T 2 differentiated nodes change from either state 0 to 1 or 1 to 2, could also be used; however, the results in this paper would be the same under both methods because the only real constraint is that edges do not form between two nodes in state 0.
If this is still confusing, consider the following example: At time T 8 there are four nodes, A, B, C, and D in the states 0, 5, 3, and 7, respectively. There are also the following edges: {A, D} {A, C} and {C, B}. At time T 9 the four nodes, A, B, C, and D are in states 0, 5, 8, and 7, respectively. Because node C changed states, only the edge {A, D} is carried over to T 9 and every other possible pair of nodes has the same probability of forming a new edge.
A great amount of integration and interconnectedness among nodes is also known as high modularity in network theory. Although the definition of modularity varies in biology, the use of the term in network theory is consistent with an operational definition of modularity sometimes used in biology (see Dassow and Munro 1999, p. 312).
References
Albert, R., & Barabási, A.-L. (2002). Statistical mechanics of complex networks. Reviews of Modern Physics, 74(1), 47–97.
Alon, U. (2003). Biological networks: The tinkerer as an engineer. Science, 301, 1866–1867.
Bascompte, J. (2009). Disentangling the web of life. Science, 325, 416–419.
Beshers, S. N., & Fewell, J. H. (2001). Models of division of labor in social insects. Annual Review of Entomology, 46, 413–440.
Bonner, J. T. (1998). The origins of multicellularity. Integrative Biology, 1(1), 27–36.
Buss, L. W. (1987). The evolution of individuality. Princeton: Princeton University Press.
Croft, D. P., James, R., & Krause, J. (2007). Exploring animal social networks. Princeton, NJ: Princeton University Press.
Dassow, Gv., & Munro, E. M. (1999). Modularity in animal development and evolution: Elements of a conceptual framework for EvoDevo. Journal of Experimental Zoology, 285, 307–325.
Fewell, J. H. (2003). Social insect networks. Science, 301, 1867–1870.
Fewell, J. H., & Page, R. E., Jr. (1999). The emergence of division of labour in forced associations of normally solitary ant queens. Evolutionary Ecology Research, 1, 537–548.
Franks, D. W., James, R., Noble, J., & Ruxton, G. D. (2009). A foundation for developing a methodology for social network sampling. Behavioral Ecology and Sociobiology, 63, 1079–1088.
Godfrey-Smith, P. (2011). Darwinian populations and natural selection (224 p). USA: Oxford University Press.
Godfrey, S. S., Bull, C. M., James, R., & Murray, K. (2009). Network structure and parasite transmission in a group living lizard, the gidgee skink, Egernia stokesii. Behavioral Ecology and Sociobiology, 63, 1045–1056.
Grosberg, R. K., & Strathmann, R. R. (2007). The evolution of multicellularity: A minor major transition? Annual Review of Ecology, Evolution, and Systematics, 38, 621–654.
Henzi, S. P., Lusseau, D., Weingrill, T., van Schaik, C. P., & Barrett, L. (2009). Cyclicity in the structure of female baboon social networks. Behavioral Ecology and Sociobiology, 63, 1015–1021.
Jackson, M. O. (2008). Social and economic networks. Princeton, NJ: Princeton University Press.
Jan, N., Stauffer, D., & Moseley, L. (2000). A hypothesis for the evolution of sex. Theory in Biosciences, 119(2), 166–168.
Jeanson, R., Kukuk, P. F., & Fewell, J. H. (2005). Emergence of division of labour in halictine bees: contributions of social interactions and behavioural variance. Animal Behavior, 70, 1183–1193.
Kauffman, S. A. (1993). The origins of order: Self-organization and selection in evolution (Vol. xviii, 709 pp.). New York: Oxford University Press.
Kauffman, S. A. (1995). At home in the universe: The search for laws of self-organization and complexity (Vol. viii, 321 pp.). New York: Oxford University Press.
Krause, J., Lusseau, D., & James, R. (2009). Animal social networks: An introduction. Behavioral Ecology and Sociobiology, 63, 967–973.
Laughlin, S. B., & Sejnowski, J. (2003). Communication in neuronal networks. Science, 301, 1870–1874.
Lusseau, D. (2003). The emergent properties of a dolphin social network. Proceedings of the Royal Society of London B (supplement), 270, S186–S188.
Maynard Smith, J., & Szathmary, E. (1995). The major transitions in evolution. New York: Oxford University Press.
McAdams, H. H., & Shapiro, L. (2003). A bacterial cell-cycle regulatory network operating in time and space. Science, 301, 1874–1877.
McShea, D. W., & Brandon, R. N. (2010). Biology’s first law: The tendency for diversity and complexity to increase in evolutionary systems. Chicago: The University of Chicago Press.
Michod, R. E. (1995). Eros and evolution: A natural philosophy of sex (Vol xxi, 241 pp.). Reading, MA: Addison-Wesley.
Naug, D. (2009). Structure and resilience of the social network in an insect colony as a function of colony size. Behavioral Ecology and Sociobiology, 63, 1023–1028.
Newman, M. E. J. (2003). The structure and function of complex networks. Society for Industrial and Applied Mathematics, 45(2), 167–256.
Newman, S. A., & Muller, G. B. (2000). Epigenetic mechanisms of character origination. Journal of Experimental Zoology, 288, 304–317.
Nowak, M. A., Tarnita, C. E., & Wilson, E. O. (2010). The evolution of eusociality. Nature, 466, 1057–1062.
Okasha, S. (2007). Evolution and the levels of selection (Vol. xi, 263 pp.). Oxford/New York: Clarendon Press/Oxford University Press.
Örçal, B., Tüzel, E., Sevim, V., Naeem, J., & Erzan, A. (2000). Testing a hypothesis for the evolution of sex. International Journal of Modern Physics C, 11(5), 973–986.
Page, R. E. (1997). The evolution of insect societies. Endeavour, 21(3), 114–120.
Page, R. E, & Mitchell, S. D. (1991). Self organization and adaptation in insect societies. In PSA: Proceedings of the biennial meeting of the philosophy of science association 1990, Vol. 2: Symposia and invited papers (pp. 289–298).
Ramos-Fernández, G., Boyer, D., Aureli, F., & Vick, L. G. (2009). Association networks in spider monkeys (Ateles geoffroyi). Behavioral Ecology and Sociobiology, 63, 999–1013.
Sih, A., Hanser, S. F., & McHugh, K. A. (2009). Social network theory: New insights and issues for behavioral ecologists. Behavioral Ecology and Sociobiology, 63, 975–988.
Tüzel, E., Sevim, V., & Erzan, A. (2001a). Evolutionary route to diploidy and sex. Proceedings of the National Academy of Science USA, 98(24), 13774–13777.
Tüzel, E., Sevim, V., & Erzan, A. (2001b). Strategies for the evolution of sex. Physical Review E, 64(6), 061908/1–061908/13.
Acknowledgments
I am eternally grateful to Eric Bair for helping me program, and for running large numbers of simulations for me. I am also incredibly thankful to Tim Schwuchow for helping me create equations to represent my model as well as for giving me moral support. Thanks also goes to Robert Brandon, Daniel McShea, Alex Rosenberg, Carl Simpson and David McCandlish, and to the Philosophy of Biology reading group at Duke University for giving me comments on my ideas when they were in their first stages.
Author information
Authors and Affiliations
Corresponding author
Appendices
Appendix 1
In the following equations, n is the total number of nodes. In the ZFEL model, because two undifferentiated nodes cannot form an edge, x represents the number of differentiated nodes. Similar to cluster calculation (see “Methods”), to calculate the probability of a cycle of size 3 (3-cycle) forming in the random model, Eq. 3 divides the number of possible graphs producing a 3-cycle by the total number of possible graphs that can be formed with n nodes and three edges. To calculate the same probability in the ZFEL model, Eq. 4 does the same thing while taking into account the number of nodes that have differentiated, represented by x. To form a cycle in the ZFEL model either all three nodes are differentiated, represented by \( \left( {\begin{array}{*{20}c} x \\ 3 \\ \end{array} } \right) \), or only two of the three nodes must have differentiated. There are three ways to create a cycle with two differentiated nodes. For example, if we look at nodes A, B, and C again, with edges {A, B}, {B, C} and {A, C}, a cycle can have either A, B or C undifferentiated as long as the other two are differentiated. This is represented by \( (n - x) \cdot \left( {\begin{array}{*{20}c} x \\ 2 \\ \end{array} } \right) \), which assumes one undifferentiated node and chooses two differentiated nodes to make a cycle. The denominators of both cycle equations were explained in “Methods”.
Lastly, to calculate the probability of a clique of size 4 (4-clique) forming in the random model, Eq. 5 divides the number of possible graphs producing a 4-clique by the total number of possible graphs that can be formed with n nodes and six edges, and for the ZFEL model, Eq. 6 does the same thing while taking into account the number of nodes that have differentiated, represented by x. Because all nodes must be differentiated and able to form edges in a clique, for the ZFEL model, the four nodes must be chosen from the number of differentiated nodes, represented by \( \left( {\begin{array}{*{20}c} x \\ 4 \\ \end{array} } \right) \). The denominators of both clique equations were explained in “Methods”.
For all of the computations in this paper, the value of n was 100, and for the ZFEL model, the following values were input for x: 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, and 100. When x is 100 in the ZFEL model, it is equal to the random model, since in both cases all nodes are able to form edges.
Appendix 2
Cycles are more difficult to measure over time than clusters because of the quick shift from few or no cycles emerging (when edge numbers are low) to a somewhat complicated mess of many interconnected cycles (when edge numbers increase). Part of the problem is, unlike clusters, in graph theory, a node can be in more than one cycle at a time, meaning that when edge numbers increase only by a little, the number of cycles gets extremely difficult for the program to count. This was also a problem with cliques because a node can be in more than one clique at the same time as well. The ZFEL model had a much larger “window” in which cycles could be counted, however, most of those results could not be compared to the random model. Therefore, the most useful data comparing the two models was the longest cycle size over five different edge formation probabilities after ten generations (see Fig. 12). There are only two data points to display with the random model (the arrow is the second data point), however, the arrow represents the steep increase expected in the random model based on individual runs even though the 100-run simulation data could not be calculated. The trend in longest cycle length in the ZFEL model also increases rapidly, but at a much slower rate.
The minimum clique size measured was two nodes (an edge). In Fig. 13, which shows the average clique size in relation to the number of cliques, it is clear the difficulty of a 3-clique or larger forming in both models, even when the number of cliques is high. This means a 3-clique would likely form in an already very highly connected network. This can be seen in Fig. 14, where the largest clique size is compared to the number of cliques, and gives similar results to Fig. 13. Even at very large numbers of cliques, both models still barely produce cliques of size four or five. These numbers are the averages over 100 simulations. Unlike the other data, because clique size is so variable, it is quite possible that within a set of parameters, one simulation run could produce a clique of say, size 6, while the majority of runs only produce cliques of size 2. It is good to keep this in mind when interpreting the results. Figures 13 and 14 have logarithmic trend lines added.
Rights and permissions
About this article
Cite this article
Fleming, L. Network Theory and the Formation of Groups Without Evolutionary Forces. Evol Biol 39, 94–105 (2012). https://doi.org/10.1007/s11692-011-9140-9
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11692-011-9140-9