Animal and Plant Species in Antagonistic and Mutualistic Networks , a Multispecies Mean Field Model

We construct a dynamic system model of a web of A N animal species interacting with P N plant species using a system of coupled differential equations. The model has a parameter which represents the effect of animal on plant species that can be positive (mutualistic interaction) or negative (antagonistic interaction). We work a Multispecies Mean Field Model MMFM in which the full set of A N and P N variables are reduced to a couple of an average animal and plant species. We study the conditions for existence of the MMFM and relate the result to the difference between animal and plant species in the network. We compare our results with empirical data from pollinators (mutualistic) and herbivorous insect (antagonistic) networks. We conclude by combining analytical and empirical work that antagonistic networks present fewer animal species in relation to plant species than mutualistic ones.


Introduction
Since the seminal paper [1], the stability complexity debate has evolved to a hot topic in theoretical ecology [2] [3].From a mathematical perspective, the stability complexity problem can be posed in at least two different contexts: firstly, the traditional approach of population dynamics, a set of coupled differential equations representing the density of a set of many interacting species [1] [4] [5] [6].In this framework, the main idea is to work the stability of the fixed points of the equations describing the time evolution of the populations [7].On the other side, the stability complexity problem can be formulated using a more contemporary framework, namely the network theory [8] [9] [10].In this approach, the focus is on the number, distribution, or assortativity of interactions, in one word: the pattern of connections in the network [11] [12].
In this work, we explore both mathematical approaches.Firstly, the coupled dynamic system is a flexible schema that permits a modeling of complex systems using phenomenological parameters such as the intrinsic growth rate of a population, its carrying capacity, or the interaction strength among species.All cited parameters can, in principle, be empirically estimated [13].In addition, most empirical data available to compare with dynamic system models are records of food webs or pollinator networks which indicate the presence of interaction among a set of species.Confident time series records of populations interacting in a community are not available, but network inventories of interaction among species are common in the literature.We cite as examples web interaction inventories [14], pollinators flower networks [15] host parasite networks [16], and herbivory networks [17].
In this manuscript, we work with a dynamic multi-species model representing two different situations: mutualistic and antagonistic interactions.Pollination webs are formed by species which both profit from the interaction and then should be favored by adaptive forces.On the other hand, Herbivory (trophic) webs are formed by bipartite webs in which only herbivores profit.This is differently from pollination where in this case the plants are injured in the interaction and adaptive forces involved in the interaction are opposite to what found in pollination.We consider a set of A N animal species and a set of P N plant species that interact with each other, but do not interact among themselves.The objective of this paper is to construct a Multispecies Mean Field Model MMFM reducing the original equations to two equations that represent the average animal and plant species.In this approach, the set of coupled dynamic equations may have an intricate topology with many fixed points, limit cycles, strange attractors, or even chaos because of the intrinsic non-linearity [18].The reduced MMFM, as we shall see, has just one single fixed point of interest.
Despite the strong simplification in this methodology, we extract useful hints about network properties, espeacially regarding the network pattern related to the stability of interacting systems, a key point in the complexity stability debate.
We are aware that the MMFM that we present here lacks a rigorous mathematical foundation.However, we persist in this theoretic field because of the promising results we have found until now.
In this paper, we consider a MMFM of a dynamic system of coupled Lotka Volterra like equations.The MMFM is a simplified two-dimension version of a multiple dimension system.The Lotka Volterra equations we deal with in our manuscript correspond to linear functional responses in the interaction term [19] [20], while we discuss more involved non-linear responses in the text.The analysis of the conditions for existence of the MMFM reveals some properties of the complex many dimensional dynamic network underlying the complex dynamics.
The main result of the cited paper of May [1] is that ecological networks are sparse.This conclusion is somewhat unexpected, we cite MacArthur that believed that large number of interaction will stabilize networks and not the opposite [21].Our work also uses a coupled of differential equations to model ecological webs.In contrast, we use different technique and our main findings concern differences between antagonistic and mutualistic networks.Our manuscript shows a theoretical model that explains the difference between the number of plant and animals in herbivorous and pollinator networks.The rest of the paper is organized as follows: in Section 2, we present the full multiple dimension dynamic system used to describe the model and the associated MMFM.In Section 3 we separately compare the existing conditions for mutualistic and antagonistic networks with empirical data.Finally, we present our main findings in Section 4 and discuss the number of animal and plants in the community.

Methodology
This section is divided in two parts.Initially, we expose the biological data used to test the methodology, namely the set of binary and quantitative interaction networks.The second part of the methodology is devoted to describe the mean field model that is in the center of our analysis.

The Biological Data Set
In this work, we employed a set of quantitative webs extracted from the literature.We use a total of total 44 n = quantitative webs divided into two categories: pollinators ( ) .Pollinator webs are a typical example of a mutualistic networks formed by flowering plants and pollinator animal species, [22] [23].Herbivory webs are a well studied example of antagonistic interaction, these networks are formed by herbivorous insect species and their respective host plants [24].
The empirical objects used in this work are ecological interaction webs which essentially are matrices of size A N versus P N corresponding to two groups of interacting animal and plant species, respectively.In the case of an interaction between species i of the first set and species j of the second the matrix element is , 1 i j a = , otherwise , 0 i j a = , for the binary case.In the more general case of quantitative matrices , i j a can assume any positive value.The connec- tivity of an animal species is defined as = ∑ and for a plant species

The Mean Field Approximation
We employ a model that is a close extension of the multispecies model developed in the reference [26].We extend the previous work to a model that encompasses both mutualistic and antagonistic interactions.The theory is based on a two-species model proposed in the reference [2] for the mutualistic case but it is even older in the antagonist situation [27].The evolution of the animal j A and plant i P abundances in time is achieved by:

A q A T A P A t A q A T A P A t P r P S P P A t P r P S P P A t
where i P and j A represent the abundances of plant i and animal j ; i r and j q are the growth rates of plants i and j while i S and j T are the respective intraspecific competition coefficients.In addition, ij α and ij β are the per-capita effects of animal j on plant i , and vice-versa.The positive signal of ij α models a mutualistic dynamic, the effect of animals on plants is positive, so it implies in an increase in plant species abundance.On the other hand, a negative signal models an antagonistic interaction in which animal species restrict plant species growth.
We adopt a mean field approach that simplifies the model allowing further analytical insight.The reference [26] used a similar reduction technique to approximate modeling.Indeed, such reductions are common in theoretical ecology [27].We assume that all plant species are equivalent, which means i r r = ; and the same average equivalence is applied to animals: Moreover, we assume in the mean field approach that , ; a similar estimation is performed on the equations of plant species.In addition, we assume that the approximations A A ≈ and are valid in the case differences inside the groups j A and i P are not significant.This con- dition is assumed in the mean field approximation.Following these estimations we obtain: Finally, we use tot We are aware that the mean field statement employed in the analytical treatment that reduces the dimension of the system from There are four solutions for this pair of coupled system: the trivial null solution ( ) 0, 0 , a solution with absence of plants ( ) 0, q T , another with absence of animals ( ) ,0 r s and a non-trivial solution with A * and P * positives given by: We note that any solution is biologically feasible only for A * and P * posi- tives.Bellow, we develop this point in detail for mutualistic and antagonistic networks.We chose to drop the ± symbol in the last equation, instead we adopt 0 α > for mutualistic interactions and 0 α < for herbivory.

Discussion of a Nonlinear Model
We may also consider more sophisticated mutualistic models that take into account a non-linear saturation process as in the reference [19] [20].We remark that the nonlinear model, as any ordinary differential equation model, is a mean field approach; to expand our methodology to a non mean field regime we would have to explicitly include the space by using differential partial equations or cellular automata models [12].The linear mutualistic model presents an explosive behavior that can be solved using saturation in the interaction response between animal-plant.In reference [20] the mutualistic interaction for one animal and one plant species A and P is modeled as: However, for low handling time h , the non-linear equations reduce to the standard linear model.The main point in model choice is perhaps: how far are the linear model predictions compared to the non-linear one.The result of reference [20] is interesting, as it shows that the critical points of the non-linear model are in a broader interval which include the linear one, and the two models coincide in the limit of small non-linearity.In addition, the analytical treatment of the non-linear model is mathematically cumbersome and less intuitive than the linear case.Moreover, the linear model allows us to explicitly make analytical inferences and to test them.Despite its simplification, the linear model produces sound theoretical results that can be compared with empirical data.

Condition for Existence of the Mean Field Model
In our work, we study network properties and analyze what the MMFM can tell about network patterns that produce viable communities.In this way, we do not focus on species parameter q , r , S and T but on the interaction network properties: α and β at one side and network parameters: We distinguish two situations in the study: the mutualistic interactions and the antagonistic interactions.The first is characterized by a positive effect of the plant animal encounter on plant population growth 0 α > whereas the second shows the opposite tendency that plant-animal interactions decrease plant growth 0 α < .In the following, we distinguish these two situations to study MMFM conditions for the existence of viable solutions.

Mutualism
The condition for non-trivial existence of the system formed by Equations ( 7) and ( 8) assumes that both A * and P * should be positive.Since the numera- tor of the right side of equation is always positive, in order to fulfill both conditions it is necessary that the denominator will be positive.Then, we have: [22].Moreover, the condition of small product αβ was explored in reference [28], where it is interpreted as the effect of asymmetric dependence in mutualistic networks.The limitations of the linear approach in this context was discussed in reference [20].

Antagonism
In the case of antagonistic interactions, we perform the transformation α α → − .In this new situation, the denominator of Equations ( 7) and (8) will always be positive, and as a consequence, Equation ( 7) is trivially positive.The fixed point related to Equation ( 8) still has an additional condition that should be fulfilled: This new condition imposed to antagonistic systems implies that A CN should be small, which means that the overall number of links in an antagonistic web should not be large.The simplest way to guarantee this condition is that the number of animals in the web, A N , should be small.We call attention that this condition is not present in mutualistic networks, therefore antagonistic webs should have a comparatively small number of animal species.We compare the number of animal and plant species for mutualistic and antagonistic cases in the next subsection.

Data Comparison
The condition for existence (11) which is exclusive of antagonistic networks claims that A N should be small.In this way, we expect that the number of animals in antagonistic networks should be more restricted than the number of plants.
The relative difference between animals and plants is presented in Figure 1.In this picture we show a bar-plot comparing the ratio

Discussion
In this paper, we work a multidimensional dynamic model that encompasses coupled equations to two equations that describe the dynamic evolution of a paradigmatic animal and plant of the community.At one hand, the MMFM implicates in a loss of information of the individual species and their interactions.On the other hand, the MMFM allows for computing the existence conditions of the system, and as a consequence, to make inference about these conditions of the full multispecies interaction network.
The essence of Mean Field approach in theoretical science consists in erasing individual differences to capture overall features of the system.In a mean field approximation, the interaction details among individuals in a network are less important than the average properties of the system dynamics.The challenge of the Mean Field Theory consists in not taking into account superfluous heterogeneities of the problem, but instead focusing on overall aspects of the phenomenon.The objective of the Mean Field theory is neither to give answers to individual differences, nor to look at small differences in the system.In this perspective, the MMFM still gives answers that can be empirically verified, despite obvious limitations and strong reductionism.The dependence asymmetry and the sparsity of the networks are positive results of the Mean Field approach and have already been discussed in the literature [26].The difference between animal and plant species has not previously explored, and our work is the first to point out this theoretical result of the MMFM.In summary, the Mean Field Theory should be judged a posteriori.If the theoretical model gives a positive answer to empirical data, only then can it be taken seriously.
The idea that antagonistic interactions are more modular like, whereas mutualistic ones are more nested was extensively explored in the literature [15] [30].The antagonistic relations tend to follow an arms race co-evolutionary strategy that produces a formation of cliques in the network, while the mutualistic interactions point to a generalist plus sub-generalist co-evolutionary strategy [31].Our manuscript works out another, not so explored, difference between antagonistic and mutualistic networks: the number of animal and plants in their composition.In fact, an alternative and appealing title for our paper could be: "Why pollinator networks have more animals than plants in their structures than herbivorous networks?"This paper essays an explanation for this question, in fact, this is the objective of our work.
It is well known that biodiversity of animal species is larger than plant species [32] [33] [34].Following this simple rule, we should generally expect that the number of animal species should surpass the number of plant species in interacting networks.This expectation could be used as a null model for plant and animal species occupancy in an interacting network.Our theoretical results suggest that it is not valid for antagonistic interaction networks, which is in agreement with empirical results.In this way, as stated in Figure 1, pollinator networks are closer to a null model of animal and plant occupancy than antagonistic networks which present a strong misbalance in their composition.
Finally, we believe that this work is important for the discussion of the complexity stability debate.We are conscious of the strong simplification of our mathematical framework.Initially we propose a model of coupled differential equations that is already a mean field approach because it does not take into account the spatial structure of the landscape.Over this model we proceed to perform a further simplification, reducing the number of variables from A P N N + to 2. In this way, our theoretic approach proceeds over a double mean field approximation: first, we exclude the spatiality, and secondly we reduce the number of degrees of freedom.A comparison of our results with empirical data suggest that despite the strong simplification of our theoretical approach, the mean field strategy is indeed a valuable tool.
The empirical data used in this work is the same of the reference [25], where we show the full characteristic of the interacting networks.The herbivorous and pollinator webs used in this work are the totality of available webs in the literature at the time of the computation analysis.That means, we have not performed any exclusion criterion in the statistical analysis.The indices of the network we use are: the number of plant P N and animal species A N , the total number of connections (without weight) among plant and animal species L , and the connectance means, we average the individual elements of the adjacency matrix by the connectance of the matrix.In addition, we use the notation tot 1 equations' average for the set of animal and plant species in the following way.We add all d d

N
, C , α or β should be small.The condition of small product A P N N × and connectance C is related to sparsity and it is a well known property of ecological networks[15]

Figure 1 .
Figure 1.The figure shows the ratio A P N N , the number of animal species divided by plant species in the network.The figure shows a larger relative number of animals in pollinator networks as compared to herbivorous networks.We depict a constant dashed line 1 A P N N = to help the visualization.