Received 26 February 2016; accepted 10 May 2016; published 13 May 2016
In recent years, the concern for mutualism is growing, since almost all species have mutualistic relationship with other species  -  . Microbial species influence on the abundances and ecological functions of related species   . Many bacterial species coexist in a syntrophic association; that is, one species lives off the products of another species. Here, we pay attention to obligate mutualism between a pair of species; one species cannot survive without the other. The systems of obligate mutualism usually have some mechanism to avoid a sudden increase (or decrease) of one species  . We discuss the reason why such a mechanism is necessary.
The relationship of obligate mutualism has the similar behaviour as the male-female relationship  -  : female (or male) cannot get offspring without the partner. The sex ratio of many animals is nearly one-to-one. Fisher (1930) first explained the reason why the 1:1 sex ratio is optimal  . Later his argument was recognized as in the framework of evolutionarily stable strategy (ESS)   . However, Fisher’s theory has some problems   . For example, his theory may conflict with real data. In most animals, the sex ratio () at birth is slightly biased ()  . Consider a large population (resident) which is slightly biased (say). According to Fisher’s theory, the resident is easily beaten by all mutants which satisfy; in particular, the mutant of asexual reproduction () is optimal  .
It is known that ESS is the most powerful tool to obtain the optimal strategy which beats the other strategies  . To win such competitions may be one of most important driving forces in evolution, but it is not a necessary condition for evolution. For instance, the most sustainable strategy is advantageous for evolution  -  . Yoshimura et al. have explained the evolutionary origins of periodical cicadas in North America not by ESS but by sustainability   . If the competition is dominant, real cicadas (13-and 17-year cycles) might be beaten by the other phenotypes with shorter cycles. The sustainability also explains the fact that the life span of species is not so long  . Similarly, Tainaka et al. (2006) have explained the evolution of 1:1 sex ratio on the basis of sustainability  . In the present paper, we extend this idea to obligate mutualism.
Lattice models have been applied to ecological problems, where simulations have been carried out by two methods: local and global interactions   . In the former, interactions occur between neighboring sites, whereas in the latter they occur between any pair of lattice sites. The simulation of global interaction is often called “lattice gas model”    ; individuals contact like the collision of gas molecules. In most cases, the dynamics of global interaction can be represented by mean-field theory. However, for the local interaction, there are usually no theories; simulations are thus necessary to obtain various results.
The most famous model of population dynamics in ecology is Lotka-Volterra equations (LVEs)   . However, they have a flaw, when we apply them to mutualism   . For example, the population sizes of both species may increase infinitely (“divergence problem”)  . Moreover, LVEs never predict Allee effect for obligate mutualism  . Recently, Iwata et al.  introduced a simple population model for mutualism, applying the lattice gas model. They assume that all interactions occur between any pair of cells. In their model, the divergence has never occurred, and Allee effect has been derived in a simple form. The present work is the local-interaction version of Iwata’s model. Especially, we explore equilibrium densities to know the sustainability of populations.
In the next section, we explain our model and method of both local and global interactions. In Section 3, the prediction by mean-field theory is reported. We briefly describe the results obtained by Iwata et al.  , and we add new results, such as equilibrium densities. Section 4 is devoted to report simulation results. In the case of local interaction, the dynamics is similar to those for global interaction, such as Allee effect. However, there is a large difference between local and global interactions. In the former case, the sustainable range in parameter space is very narrow, compared to global interaction. Namely, the proportion (ratio) between both species is strongly restricted for their survival. In the final section, we present an idea of “ratio selection” which is an extension of the sex-ratio selection  -  .
2. Model and Methods
Consider a system consisting of two mutualistic species X and Y (Figure 1). Each lattice site is labeled by X, Y or O, where O means the empty site. The empty cell is introduced to prohibit the divergence of population sizes  . The reactions are defined by
Figure 1. Lattice model for obligate mutualism. The cell X (Y) indicates the occupation site (individual) of species X (Y), and O is empty. The empty cell is introduced to avoid the divergence of population sizes  .
The reactions (1a) and (1c) respectively denote the birth and death processes of species X, and () denotes the reproduction (mortality) rate of species X. Similarly, the reactions (1b) and (1d) mean the birth and death processes of species Y, respectively.
Simulations of lattice model are usually carried out by either local or global interaction   . First, we explain the simulation procedure of local interaction. Reaction processes are performed in the following steps:
1) Initially, we distribute X, Y on a square lattice. Each cell is one of three sites: X, Y, and O (see Figure 1).
2) To update, we choose a target site randomly.
a) Reaction (1c) or (1d): If the target cell is X (or Y), then it becomes O by the rate (or).
b) Reaction (1a) or (1b): If the target site is O, we choose another cell from the neighboring 8 sites (Moore neighborhood) around the target cell. When the second site is X (or Y), then the target cell become X (or Y) by the rate (or).
3) We repeat step 2) by times, where is the lattice size and is the total number of sites (). This step is a unit time called the Monte Carlo step.
4) We further continue the updates, until the system reaches stationary state.
The reproduction rates and for local interaction are assumed that
where we call () the birth rate of species X (Y), and () the neighboring density of species X (Y). For instance, if there are five X cells in the Moore neighborhood, we put.
Next, we explain the method of global interaction (lattice gas model). Almost all procedures are the same as local interaction. However, in the global simulations, the birth process [reaction (1a) or (1b)] occurs between any pair of lattice sites. In mean-field limit, the reproduction rates are replaced as follows:
where () is the overall density of species X (Y).
3. Prediction of Global Interaction
If the reaction (1a) occurs between any pair of lattice sites, and if the lattice size is sufficiently large, then the mean-field theory holds:
where the factor in the right hand sides denotes the density of empty cell. The first and second terms in Equation (4) come from birth and death processes, respectively. Inserting Equations (3) into (4), we have
Note that this equation is almost the same as in the male-female system   . Iwata et al.  proved that the Equation (5) led to one of two phases: extinction or Allee-effect phase (Figure 2). In the former case, both species always go extinct. The latter is called the survival/extinction (S/E) phase, because this phase has two stable equilibriums which means survival and extinction. The condition for the Allee-effect phase has been represented by
In summary, two conditions are necessary for survival; one is the relation (6), and the other is that initial densities of two species are higher than the separatrix which is schematically shown in Figure 2(b).
We explicitly obtain equilibrium density in stable state, when
Figure 2. Schematic illustration of population dynamics. (a) extinction phase, and (b) Allee effect phase  .
From Equations (5) and (7), we get
From Equation (8), we can find the following relation: for. This means that the equilibrium density satisfies
Inserting Equations (7) and (10) into (5), we have
Equation (11) denotes a typical Allee-effect equation  . According to the sign of D, the dynamics is classified into two phases: extinction phase for and Allee-effect phase for. So long as initial densities of both male and female are high enough, the steady-state densities in survival phase is given by
It is therefore found in the mean-field limit that 1) the maximum value of total density () is achieved at, 2) a positive stable equilibrium (Allee-effect phase) exists, when or
This inequality is the same as equation (6) under the condition of.
4. Simulation Results
We report simulation results of individual-based (lattice) model. In the case of global interaction, the mean-field theory [Equation (5)] predicts that the dynamics has two phases. One is an extinction phase: both species always go extinct. The other is Allee-effect phase. A typical example of the latter phase is shown in Figure 3 (), where (1) represents simulation result and (2) is the prediction by mean-field theory. All orbits asymptotically reach one of two stable equilibria (filled circles). Both equilibria locate on the dotted line represented by Equation (10). The densities of both species at survival equilibrium (right filled circle) are expressed by Equation (13a). Unless both densities are considerably high, both species goes extinct. The red curve denotes the separatrix which determines whether the final equilibrium is survival or extinction state. The unstable equilibrium (open circle) locates on the separatrix. We confirm the results of global simulation are well predicted by mean-field theory [Equation (5)].
In the case of local interaction, the dynamics are similar to the predictions of lattice gas theory. For instance, the Allee effect can be observed: unless both densities of species are considerably high, the population goes extinct. However, in the case of local interaction, simulation results exhibit a distinct difference from those of global interaction. In Figure 4, typical spatial distributions in stationary state are displayed (), where (1) and (2) are the patterns of global simulation, and (3) and (4) are those of local simulation. The value of is optimal at; the total population size takes the maximum value at [see Figure 4(a) and Figure 4(c)]. In contrast, both Figure 4(b) and Figure 4(d) display at. It is found from Figure 4(d) that the total population size largely decreases, when the value of slightly deviates from the optimality.
Figure 3. Typical cases of population dynamics in Allee-effect phase (). (a) simulation, and (b) theory (Equation (5)). We use (). Each orbit asymptotically reaches one of two stable equilibriums (filled circles). All equilibriums locate on the dotted line defined by Equation (10). Red curve denotes the separatrix which is obtained by simulation.
(a) (b)(c) (d)
Figure 4. Typical spatial patterns in stationary state (). (a) and (b) Global interaction. (c) and (d) Local interaction. The value of is 0.5 for (a) and (c), and 0.53 for (b) and (d). Blue and red cells denote species X and Y, respectively.
Figure 5. Steady-state densities at stable equilibrium. (a) Global interaction, but (b) local interaction. The densities at stable equilibrium are plotted against, where we use,. Each plot is obtained by the average over for global interaction, but for local interaction. The initial condition is set to be the random distribution without empty cell in order to avoid extinction. The total density of both species has a peak at (dotted line).
ratio of densities. However, it is necessary for the sustainability in local interaction that the value of should be very close to 0.5 [see Figure 5(b)]. In other word, both species must have 1:1 ratio for the sustainability.
When, the sustainable range becomes biased. Figure 6 is the same figure as Figure 5, but the mortality rate of X is much higher than that of Y (). In global simulation, almost all ratios between species X and Y are sustainable. However, in local simulation, both species can coexist near where the total density takes the maximum value. The ratio () of both birth rates is largely X-biased; namely, the birth rate () of species X must be higher than for the survival of both species. The survival condition is strongly restricted in local simulation.
We have developed a lattice population of obligate mutualism. Local interaction strongly effects on the populations dynamics. The survival condition is extremely limited for local simulation. When both species have the same mortality rate (), it is necessary for the sustainability in local interaction that both species must have 1:1 ratio (see Figure 5). When the mortality rate of species X is larger than that of Y (), the birth rate of X must be larger than that of Y (see Figure 6). In the case of local interaction, a severe condition is required. Individuals of both species must live together in close proximity. On the basis of these results, we present a theory of “ratio selection” for obligate mutualism. When two species interact locally (long-range interaction is prohibited), both species require a specific (proper) ratio. The proper ratio may be determined by various factors, such as mutual body sizes, mutual shapes and obligate requirements of both species.
The ratio selection is a generalization of sex-ratio selection. In fact, our model is also applicable to the male- female system. Equation (5) is almost the same as obtained in male-female system   , where X (Y) corresponds to male (female) and means the sex ratio at birth. Fisher first explained the evolution of 1:1 sex ratio  . However, as described before, his theory (ESS) cannot explain real data of male-biased sex ratio    . On the other hand, the present work well explains the male-biased sex ratio (). As shown in Figure 6(b), the biased sex ratio is explained by the fact that the mortality rate of male is slightly higher than that of female (). Since males rapidly die than females, the inequality is necessary to keep mutual encounter (close proximity). According to ESS, the strategy of is optimal. However, this strategy cannot survive, because its steady-state density is zero [see Figure 6(b)].
Many empirical data suggest the ratio selection; two-species systems of obligate mutualism usually have some mechanism to avoid an abrupt increase of one species  . 1) Coral and algae: the ratio between cell number in stomach (“gastrodermis”) of coral and zooxanthellae number is just 1:1  . When the density of algae is too high, the excess algae are excluded from coral. 2) Chlorella and Paramecium: the number of Chlorella (Chlorella) inside each body of Paramecium (Paramecium bursaria) is about 400 individuals  . 3) Yucca and
Figure 6. Same as Figure 5, but for; we use and. The total density of both species has a peak, when the value of is larger than as indicated by dotted line.
yucca moth: too many yucca moths become harmful for yucca plant  . 4) Chloroplast and mitochondria: the number of these organelles in a host cell is not so arbitrary  . The ratio selection comes from the requirement that all individuals of one species need a high local density of the other species. We consider such requirement evolutionarily leads to endosymbionts from free-living ancestors  . This is because endosymbiosis is more stable than free living to keep their proper ratio. Similarly, the ratio selection may be a driving force for the evolution from endosymbionts to organelles   .
The authors sincerely thank to professors Jin Yoshimura for valuable comments. They also thank to Mr. Keiji Amemiya for the help of simulations.
 Bashary, R. and Bronstein, J.L. (2004) Game Structures in Mutualisms: What Can the Evidence Tell us about the Kind of Models We Need? Advances in the Study of Behavior, 34, 59-101.
 Pellmyr, O. and Huth, C.J. (2002) Evolutionary Stability of Mutualism between Yuccas and Yucca Moths. Nature, 372, 257-260.
 Yokoi, H., Uehara, T., Kawai, T., Tateoka, Y. and Tainaka, K. (2014) Lattice and Lattice Gas Models for Commensalism: Two Shellfishes in Intertidal Zone. Open Journal of Ecology, 4, 671-677.
 Goto, R., Okamoto, T., Kiers, E.T., Kawakita, A. and Kato, M. (2010) Selective Flower Abortion Maintains Moth Cooperation in a Newly Discovered Pollination Mutualism. Ecology Letters, 13, 321-329.
 Berec, L. Boukal, D.S. and Berec, M. (2001) Linking the Allee Effect, Sexual Reproduction and Temperature-Dependent Sex Determination via Spatial Dynamics. American Naturalist, 157, 217-230.
 Tainaka, K., Hayashi, T. and Yoshimura, J. (2006) Sustainable Sex Ratio in Lattice Populations. Europhysics Letters, 74, 554-559.
 Iwata, S., Kobayashi, K., Higa, S., Yoshimura, J. and Tainaka, K. (2011) A Simple Population Theory for Mutualism by the Use of Lattice Gas Model. Ecological Modelling, 222, 2042-2048.
 Fisher, R.A. (1930) The Genetical Theory of Natural Selection. Clarendon Press, Oxford.
 Smith, J.M. (1968) Mathematical Ideas in Biology. Cambridge University Press, Cambridge.
 Shaw, R.F. and Mohler, J.D. (1953) The Selective Advantage of the Sex Ratio. American Naturalist, 87, 337-342.
 Cowgill, U.M. and Hutchinson, G.E. (1963) Differential Mortality among the Sexes in Childhood and Its Possible Significance in Human Evolution. Proceedings of the National Academy of Sciences of the United States of America, 49, 425-429.
 Yoshimura, J. (1997) The Evolutionary Origins of Periodical Cicadas during Ice Ages. The American Naturalist, 149, 112-124.
 Tanaka, Y., Yoshimura, J., Simon, C., Cooley, J.R. and Tainaka, K. (2009) The Allee Effect in the Selection for Prime-Numbered Cycles in Periodical Cicadas. Proceedings of the National Academy of Sciences of the United States of America, 106, 8975-8979.
 Tainaka, K. (2003) Perturbation Expansion and Optimized Death Rate in a Lattice Ecosystem. Ecological Modelling, 163, 73-85.
 Tainaka, K. (1988) Lattice Model for the Lotka-Volterra System. Journal of the Physical Society of Japan, 57, 2588-2590.
 Tainaka, K. and Araki, N. (1999) Press Perturbation in Lattice Ecosystems: Parity Law and Optimum Strategy. Journal of Theoretical Biology, 197, 1-13.
 Frisch, U., Hasslacher, B. and Pomeau. Y. (1986) Lattice-Gas Automata for the Navier-Stokes Equation. Physical Review Letters, 56, 1505-1508.
 Hagiwara, T., Ushimaru, T., Tainaka, K., Kurachi, H. and Yoshimura, J. (2011) Apoptosis at Inflection Point in Liquid Culture of Budding Yeasts. PLoS ONE, 6, e19224.
 Hofbauer, J. and Sigmund, K. (1998) Evolutionary Games and Population Dynamics. Cambridge University Press, Cambridge.
 Hardy, I.C.W. (2002) Sex Ratios: Concepts and Research Methods. Cambridge University Press, Cambridge.
 Hamilton, W.D. (1967) Extraordinary Sex Ratios. Science, 156, 477-488.
 Rosenberg, E., Koren, O., Reshef, L., Efrony, R. and Zilber-Rosenberg, I. (2007) The Role of Microorganisms in Coral Health, Disease and Evolution. Nature Reviews Microbiology, 5, 355-362.
 Miwa, I. (2009) Regulation of Circadian Rhythms of Paramecium Bursaria by Symbiotic Chlorella Species. In: Fujishima, M., Ed., Endosymbionts in Paramecium, Springer-Verlag, Berlin, 83-110.
 Fiegna, F., Yu, Y.-T.N., Kadam, S.V. and Velicer, G.J. (2006). Evolution of an Obligate Social Cheater to a Superior Cooperator. Nature, 441, 310-314.
 Milo, R. and Philips, R. (2014) Cell Biology by the Number. http://book.bionumbers.org/
 Hosokawa, T., Ishii, Y., Nikoh, N., Fujie, M., Satoh, N. and Fukatsu, T. (2016) Obligate Bacterial Mutualists Evolving from Environmental Bacteria in Natural Insect Populations. Nature Microbiology, 1, Article Number: 15011.
 Sagan, L. (1967) On the Origin of Mitosing Cells. Journal of Theoretical Biology, 14, 255-274.