ABSTRACT
Maintaining sustainable ecosystems are important for all the inhabitants of earth. Also, an important component of sustainable ecosystems is the maintenance of healthy coexistence of consumers and their resources which can include diseases in the species involved. We formulate a model, where the resources are plants, to explore how consumer-resource coexistence could of itself limit the spread of infectious diseases. The important mathematical features of the model are discussed using the basic reproduction number and the consumption number. The results show an association between species coexistence and a decrease in ecosystem resource disease. The possible importance of these results are discussed.
1. Introduction
Ecosystem sustainability is of obvious benefit to overall regional ecology and geography as well as benefiting humans, economically and socially. There are also less obvious benefits. This study aims to illustrate that one of these possible hidden benefits could be how ecosystem sustainability, through species coexistence, could affect the spread of diseases. Throughout the manuscript coexistence is used to describe a natural system where consumers and their resources both persist (survive). Due to the fact that consumers depend on their resources the persistence of consumers indicates this coexistence. Thus, when consumers are absent this coexistence has broken down.
In human and some animal populations control of infectious diseases can be done by applying appropriate control measures such as vaccination, treatment, quarantine, etc. For instance, the spread of rabies to wildlife can be reduced by oral vaccination of raccoons [Citation41]. In natural ecosystems it is not easy to apply these types of control and other methods of reducing the spread of infectious diseases are needed. For example, it has been shown that maintaining biodiversity can reduce disease transmission [Citation29]. That infectious diseases can reduce ecosystem sustainability is reasonably obvious and understood [Citation15]. On the other hand, the possibility that ecosystem sustainability can reduce disease is less intuitive. This study explores how maintaining ecosystem sustainability could potential reduce the spread of infectious diseases in an ecosystem. We consider one aspect contributing to sustainable ecosystems, the ability of species to coexist, in this case consumers and their resources.
The field of mathematical epidemiology, in which disease progression is modelled, originated from work by Daniel Bernoulli in the 1700s [Citation6] and was developed further by a number of people from ideas of Ronald Ross [Citation42]. These developments included the influential papers of McKendrick and Kermack [Citation6]. Robert M. May and Roy M. Anderson made important contributions including the concept of the basic reproduction number defined as the average number of infections caused by susceptible individuals [Citation1]. Epidemiological consumer-resource models with resources (prey or plants) and consumers (predators or herbivores) have since been given considerable attention [Citation2,Citation31,Citation40,Citation43,Citation47]. Less of these studies consider infection in the resources, as used here, where these are plants [Citation36].
The relationship between a consumer and infected resource can be complex with existing evidence of either increases or decreases in resource infections in response to consumers [Citation18,Citation25,Citation27,Citation39]. For the first scenario, the idea is that consumption helps to eliminate infectious individuals thus preventing the spread of disease [Citation25,Citation39] and in the second scenario predators facilitate the spread of disease [Citation18,Citation27].
In the presence of a consumer and diseased resource, the overall system can show different complex dynamical behaviours, such as bistability, quasi-periodicity and chaos [Citation10,Citation32]. Here we use the classic Rosenzweig–MacArthur model approach for consumer-resource dynamics which incorporates reasonable biology in the form of logistic resource growth and a saturating functional response of resource consumption, and can produce steady-state or cyclical dynamics [Citation28].
Epidemiologically, the basic reproduction number is a threshold quantity that describes a condition under which an epidemic can be eradicated from a system [Citation8,Citation44,Citation46]. On the other hand, the consumption number denoted by can be understood as the parameter combination ensuring the resource consumption required for survival or coexistence of the consumer with its resource [Citation13,Citation17].
For dynamical systems, a bifurcation is understood as an abrupt topological change in the dynamics associated with a smooth change in parameter values (the bifurcation parameters). Epidemic models can have a variety of bifurcation behaviours at [Citation6]; see [Citation5] for examples. In particular, for a forward bifurcation, this change does not change the number of stable states. In contrast, a backward bifurcation occurs when there is a change from one stable state to two stable states (bistability). For a forward bifurcation at , the equilibrium of the infected resource remains zero so long as and then increases continuously as increases [Citation5,Citation44]. For a forward bifurcation reducing this parameter below results in a return to an infected resource of zero. In contrast, for a backward bifurcation it can be necessary to bring well below 1 to eradicate the disease [Citation3,Citation19,Citation22,Citation37]. Therefore, in understanding infection dynamics, it is helpful to consider if a backward bifurcation is possible. Another type of bifurcation is a change from a steady-state solution to cyclical dynamics, called a Hopf bifurcation, which can also be important for disease progression [Citation4,Citation11].
Using our model and these threshold concepts, we show the types of bifurcations and illustrate the possibility of how strengthening the coexistence of species could lead to eradication of certain diseases from the system where not expected, or speed up the expected eradication.
2. Models and methods
A simple mathematical model methodology is used to investigate the dynamics of an ecosystem with disease in the resources. The model is formulated using the following assumptions. The key players are resources, consumers and pathogens. All variables are measured as densities of biomass. Resources are made up of a susceptible (uninfected) group S and an infected group I. These resources are consumed by herbivores, the consumers (Y ). P is the density of pathogen in the environment [Citation12,Citation14,Citation35,Citation44]. The pathogen attack the plants with the capacity of causing damage which prevents consumers feeding on them. The nature of interactions between S and Y as well as S and P are non-linear [Citation34]. Thus, a non-linear function, a hyperbolic functional response, is used to model interactions between S and Y. This function is chosen as it has been shown to be an appropriate model for the interaction of consumers and resources, where the resources are plants [Citation45]. The force of infection is the rate at which susceptible resources become infected from the pathogen and is given by , where δ is the pathogen concentration at which the force of infection is at half of its upper limit [Citation12,Citation35]. The growth of susceptible (uninfected) resources are modelled as logistic [Citation36]. The pathogen density P increases based on the amount of infected resource at a rate ν and dies at a rate ξ. Based on these assumptions and formulations, we obtain the model (1) (1) where S and I denote densities of susceptible and infected resources at any particular moment in time. Similarly, Y and P denote densities of consumers and pathogen in the environment. For all the numerical calculations used in this study these densities are in . The meanings of parameters are given in Table .
Table 1. Variables/parameters for the model.
Analyses of this model reveal that it has many equilibrium points. However, the most interesting of these is when all variables are positive. In particular, we are interested in situations where the disease can be eradicated with the survival of consumers and resources.
The basic reproduction number is calculated because this is where any potential bifurcation will occur with eradication of the disease for . However, the possibility of a backward bifurcation is also investigated because then decreasing below 1 is not sufficient to ensure disease eradication. The basic reproduction number for model (Equation1(1) (1) ) is computed using the next generation matrix method [Citation46] and is (2) (2) The consumption number gives a condition under which the equilibrium points of the system are stable [Citation13,Citation17]. It has been shown that the threshold requirement for consumer survival is that for a broad range of consumer resource models [Citation17]. Thus, this quantity provides insight into the strength of coexistence of resources and consumers because for values further above 1 the likelihood of being forced below 1 by external drivers is reduced. Another way of explaining this is that for a change in parameters such that will require greater effort. is calculated using the approach in [Citation13] and is (3) (3) The consumer-resource-disease dynamical system is investigated using the concepts of the basic reproduction number and the consumption number using mathematical analyses and numerical results.
3. Results
3.1. Model analyses
The system (Equation1(1) (1) ) has many equilibrium points: (4) (4) (5) (5) (6) (6) (7) (7) (8) (8) where and satisfy the equation: (9) (9) and . Algebraically, equation (Equation9(9) (9) ) becomes (10) (10) Note: exists for , exists for and exists for and .
For the equilibrium points and the consumer is absent (i.e. a lack of coexistence) while for and the consumer is present (i.e coexistence). From these analyses a necessary condition for coexistence is that . Also, is simply a trivial equilibrium point. The remaining equilibrium points coincide (i.e. ) if . This suggests that and are bifurcation quantities. We present a bifurcation analyses of the model (Equation1(1) (1) ) to explore this possibility further. First we consider stability about the equilibrium points.
3.1.1. Stability analyses
The short and long term dynamics of the model can be described by stability conditions about its equilibrium points [Citation33] and these are summarized in the following theorem. To discuss the impact of consumers on the dynamics of the system, the stability results are presented for the following cases: consumers absent (a lack of coexistence) and consumers present (coexistence). For and stability for the general cases are difficult to prove. Where possible we do investigate this stability analytically here and support these results by further numerical simulations.
Theorem 3.1
no coexistence
The trivial equilibrium point is unstable.
The equilibrium point is stable if and .
For some the equilibrium point is stable provided .
Proof.
The Jacobian of model (Equation1(1) (1) ) evaluated at the trivial equilibrium point has the following eigenvalues: (11) (11) Thus, is unstable since .
The Jacobian of model (Equation1(1) (1) ) evaluated at the equilibrium point has the following eigenvalues: (12) (12) (13) (13) Clearly, and . Thus, the equilibrium point is stable if and .
To simplify our analyses, we make the following assumptions: . Based on these assumptions, and the characteristic polynomial of the Jacobian of model (Equation1(1) (1) ) evaluated at the equilibrium point becomes (14) (14) where . A necessary and sufficient condition for to be stable is the Routh Hurwitz stability criterion: for i=0,1,2,3,4 and and . By algebraic manipulations we establish that these conditions are satisfied as follows: and Thus, for some , the equilibrium point is stable provided .
Ecologically, for a lack of coexistence, Theorem 3.1 implies the following: (i) consumers and pathogens cannot all go extinct (ii) consumers and infected resources can be eradicated when and (iii) for and , infected resources can persist in the system (i.e. an endemic steady state).
Theorem 3.2
coexistence
For :
the equilibrium point is stable if .
the equilibrium point is stable for some
(i.e. ).
Proof.
The Jacobian of model (Equation1(1) (1) ) evaluated at the equilibrium point has the following eigenvalues: (15) (15) (16) (16) where . Clearly, and .
. Given that , . Observe . By elementary algebraic calculations we have that Thus, is stable if and .
From (a), is stable if and . Since , then for some (i.e. ). Taking into account that , we also conclude that for some (i.e. ), is stable if .
Note that for and , the equilibrium point might be unstable. Using the analytical expression of the Jacobian evaluated at the equilibrium point results in overly complicated eigenvalues. Thus, to show where might be unstable we consider special cases. Various sets of parameter values such that and , resulted in the eigenvalues having a positive real part. For instance, considering parameter values such that and which is greater than corresponds to eigenvalues: . This suggest that for and , is likely to be unstable.
Ecologically, for coexistence of the consumer and its resources, Theorem 3.2 implies the following: (i) for , (a) infected resources are eradicated when and (b) infected resources can also be eradicated for some (for to exist must be positive and thus (see above) and then for (which is ecologically essential) will be larger than unity and this then implies the coexistence of resources and the consumer) (ii) for and , is unstable. We explore this possibility further, numerically, below.
3.1.2. Bifurcation analyses
Bifurcation analyses are presented here to improve in an understanding of the dynamics of model (Equation1(1) (1) ). In particular, bifurcation near the equilibrium points is explored (local bifurcations).
Theorem 3.3
In the absence of consumers, model (Equation1(1) (1) ) has a forward bifurcation about .
Model (Equation1(1) (1) ) does not have a backward bifurcation about for .
A Hopf Bifurcation occurs at , about the equilibrium point .
Proof.
For the case with no consumers, we consider the disease-free equilibrium () and endemic equilibrium () ( is trivial). Forward bifurcation is where the disease-free equilibrium loses its stability and a stable endemic equilibrium appears as increases through 1 [Citation7,Citation19,Citation22]. is a stable disease-free equilibrium for model (Equation1(1) (1) ) when (see Theorem 3.1 for the proof). When , become a stable endemic equilibrium for model (Equation1(1) (1) ) (see Theorem 3.1 for the proof). Note that when , the two equilibrium points coincide (i.e. ). Based on these arguments, model (Equation1(1) (1) ) has a forward bifurcation at . Thus, in the absence of consumers, we expect disease eradication for model (Equation1(1) (1) ) when .
Backward bifurcation is a type of bifurcation where a stable endemic equilibrium co-exists with a stable disease-free equilibrium when [Citation7,Citation19,Citation22]. Considering (Equation10(10) (10) ) the following points hold for equation (Equation9(9) (9) ), there is
a unique positive solution for ,
no positive solution for ,
no positive solution for .
Thus, for equilibrium point there is no endemic equilibrium when and so it does not have a backward bifurcation about .
A Hopf bifurcation occurs when a pair of eigenvalues crosses the imaginary axis. From the stability analyses in Theorem 3.1, when the Jacobian evaluated at has negative real eigenvalues with zero imaginary parts. Also, when the Jacobian evaluated at has a pair of complex conjugate (pure imaginary) eigenvalues. Next, when , the Jacobian evaluated at has a pair of complex conjugate eigenvalues with positive real part, provided . Therefore, increasing through , a pair of complex conjugate eigenvalues emerges and crosses the imaginary axis from left to right. Thus, a Hopf bifurcation occurs at about the equilibrium point .
3.2. Numerical examples
Of interested is disease eradication and it is already known that diseases are generally eradicated when the basic reproduction number is less than one and persist when the basic reproduction number is greater than one. Thus, we perform numerical simulation for these two categories using the parameter values given in Table : and .
Table 2. Parameters values used for model simulation with their source (σ and onwards were estimated in order to yield reasonable model behaviour).
Firstly, for , using the parameter values (Table ) results in and the dynamics given in Figure . For Figure (a–d), σ is varied to obtain a range of such that: (a) for () it takes about 200 years for the infected resources to be eradicated, (b) for () it takes about 150 years for infected resources to be eradicated, (c) for () it takes about 100 years for infected resources to be eradicated, and (d) for ( ) it take about 50 years for infected resources to be eradicated. So, for fixed, increasing results in faster eradication of the disease from the ecosystem.
Secondly, for , again using the parameter values in Table , changing β results in and the dynamics given in Figure . Again, Figure (a–d) was obtained by varying σ to get a range of while keeping fixed. The results from Figure are: (a) for () the infected resources increase to a maximum biomass of about 30 g/m2, (b) for () the infected resources increase to a maximum biomass of about 10 g/m2, (c) for () the infected resources decrease to a minimum biomass of about 1.0 g/m2, (d) for () the infected resources decrease to zero biomass (total eradication of infected population). Each of these dynamics occur over approximately 200 years except for the last one that is over a period of 100 years (Figure (d)).
3.3. Generalized model
To show how the results could be extended to more complex models we consider an ecosystem which consists of multiple patches with resources and consumers together with pathogens causing infectious diseases. We make the simplifying assumption that the disease is transmitted through direct contact with the pathogen and does not spread from one patch to another. Thus, while movement between and within patches is important, this is ignored. Also, the pathogen is assumed to affect each patch equally. Based on these assumptions we obtain the generalized model: (17) (17) The meaning of variables and parameters are the same as given in Table except for the subscript which represent the patch for the resources or consumer.
3.3.1. Generalized model: Preliminary analysis
The basic reproduction number for the generalized model (Equation17(17) (17) ) is also computed using the next generation matrix method [Citation46] and is (18) (18) where .
The consumption number for the generalized model (Equation17(17) (17) ) is also calculated using the approach in [Citation13,Citation46] and is (19) (19) where .
Thus, because the basic reproduction and the consumption numbers for this more general model can also be calculated it should also be possible to use them to investigate the system dynamics as illustrated using the simplified version of the model. This possibility needs further analyses and will be considered in future work.
4. Discussion
In summary, a consumer-resource model is presented here where resources are exposed to disease. The basic reproduction number and the consumption number are calculated and shown to govern the dynamics in agreement with fundamental results in the literature [Citation8,Citation13,Citation17,Citation44,Citation46].
As explained is the introduction understanding the bifurcation dynamics at can be important [Citation6]. For our system we show that there is no backward bifurcation and so going through in either direction does not affect the qualitative behaviour of the dynamics. We show that a Hopf bifurcation occurs for a prescribed combination of parameter values and so cyclic dynamics are possible as found in other epidemiological models [Citation4,Citation10,Citation26]. These cycles can be important as they can create further cyclical dynamics in the disease itself [Citation9,Citation28]. Also, if species densities are cycling then they can be vulnerable to environmental stochasticity [Citation20].
It is well known in epidemiology that when the basic reproduction number is less than unity there is a greater chance of eradicating a disease. Using our model we show that forward bifurcation is possible and the disease is eradicated when as expected [Citation8,Citation44,Citation46]. However, we also show here that increasing the consumption number results in faster eradication of the disease from the system. Thus, increasing the strength of coexistence between the consumer and the resources results in faster disease elimination. As far as we know this effect on the rate of disease elimination has not been demonstrated previously and could be important in situations where stochastic events prevent elimination if the process is slow [Citation21].
When , the model can show disease persistence, as expected. However, here we show that consumers can change this result. In fact, mathematical analyses and numerical examples show how with consumers present even with the disease can be eradicated for some large enough. As explained in the introduction, larger values of will improve the likelihood of coexistence of resources and consumers. These results suggest an association between strengthening the coexistence of species and the eradication of certain diseases from a system. This association highlights the possible importance of future research on the extent to which interventions that support coexistence might reduce the spread of infectious diseases (including the prevalence and consequences of such diseases). Moreover, if future research indicted a causal relation between species coexistence and disease prevalence then a possible consequence could be that consumer extinction (through external factors) would increase diseases, further intensifying pressure on plant diversity as well.
One possible mechanism for this type of process is that consumption eliminates infectious individuals thus preventing the spread of disease [Citation25,Citation39]. Also, in a predator prey situation this process can be exacerbated by the possible weakening of prey by disease making them more vulnerable to predation [Citation25].
Infectious diseases in natural resources are common and inevitable. This study provides epidemiological reproduction and consumption numbers for understanding how such a disease might operate in a consumer resource system. These are biologically conceivable quantities and in principle could be measured [Citation26]. Thus, field studies together with the models could be used to predict disease progression. For example, for a particular ecosystem, field estimates for variables such as carrying capacity, conversion and consumption rates of resource to consumer biomass, and the other variables that constitute could be used together with measures of the parameters in to determine the possible disease dynamics. In particular, eradication or not of the disease can then be assessed. The resilience of the ecosystem can then be explored under scenarios of pressure from factors such climatic variability, human settlement or poaching. Also, if it is ascertained that consumers (predators or herbivores) are important for disease dynamics this knowledge could be important when managing culling or hunting programmes.
The original model described here is simple but gives the fundamental characteristics of possible consumer-resource systems. However, it still displays complex behaviour and its simplicity allows a more complete analysis [Citation20]. Many possible extensions can be used to effectively describe actual situations. We use one example of a more complex model, with multiple patches and consumers, resources and disease, and the same results are likely to hold. Also, this model is made with a particular set of assumptions that could be altered to better suit specific ”real world” systems, and these alternative models could be worth pursuing in future research efforts.
It is understood to a certain extent that in ecosystems the occurrence of infectious diseases can have a negative effect in reducing ecosystem sustainability by affecting the coexistence of species. Alternatively, it has been shown that ecosystem biodiversity can affect disease transmission rates [Citation30]. Thus, it appears that ecosystem sustainability could affect levels of infections diseases. Here we show how coexistence of consumers and their resources could contribute to reduced levels of disease. This possibility could be important for more complicated systems and the overall persistence of ecosystems. A further advantage is that if disease is reduced in natural systems then this could also reduce infections in domesticated species in which there is contact, such as food crops, with possible economic and public health cost benefits [Citation24]. Lastly, if incidences of disease increase as consumers are driven extinct then eventual declines in resource diversity are also a possible consequence.
Disclosure statement
No potential conflict of interest was reported by the authors.
Additional information
Funding
References
- R.M. Anderson and R.M. May, The population dynamics of microparasites and their invertebrate hosts, Phil. Trans. R. Soc. Lond. B 291(1054) (1981), pp. 451–524.
- P. Auger, R. Mchich, T. Chowdhury, G. Sallet, M. Tchuente, and J. Chattopadhyay, Effects of a disease affecting a predator on the dynamics of a predator–prey system, J. Theor. Biol. 258(3) (2009), pp. 344–351.
- M. Barfield, M. Martcheva, N. Tuncer, and R.D. Holt, Backward bifurcation and oscillations in a nested immuno-eco-epidemiological model, J. Biol. Dyn. 12(1) (2018), pp. 51–88.
- A.M. Bate and F.M. Hilker, Complex dynamics in an eco-epidemiological model, Bull Math. Biol.75(11) (2013), pp. 2059–2078.
- F. Brauer, Backward bifurcations in simple vaccination models, J. Math. Anal. Appl. 298(2) (2004), pp. 418–431.
- F. Brauer, Mathematical epidemiology: Past, present, and future, Infect. Dis. Model. 2(2) (2017), pp. 113–127.
- C. Castillo-Chavez, K. Cooke, W. Huang, and S.A. Levin, Results on the dynamics for models for the sexual transmission of the human immunodeficiency virus, Appl. Math. Lett. 2(4) (1989), pp. 327–331.
- C. Castillo-Chavez, Z. Feng, and W. Huang, On the computation of R0 and its role on global stability in Mathematical Approaches for Emerging and Re-emerging Infectious Diseases: An Introduction, IMA 125, Springer-Verlag; 2002
- R.D. Cavanagh, X. Lambin, T. Ergon, M. Bennett, I.M. Graham, D. van Soolingen, and M. Begon, Disease dynamics in cyclic populations of field voles (Microtus agrestis): Cowpox virus and vole tuberculosis (Mycobacterium microti), Proc. R. Soc. B Biol. Sci. 271(1541) (2004), pp. 859–867.
- S. Chakraborty, B.W. Kooi, B. Biswas, and J. Chattopadhyay, Revealing the role of predator interference in a predatorprey system with disease in prey population, Ecol. Complex. 21 (2015), pp. 100–111.
- J. Chattopadhyay and O. Arino, A predator-prey model with disease in the prey, Nonlinear Anal. Theory Methods Appl. 36 (6), (1999), pp. 747–766.
- C.T. Codeo, Endemic and epidemic dynamics of cholera: The role of the aquatic reservoir, BMC Infectious Diseases 1 (1), (2001), pp. 209.
- O.C. Collins and K.J. Duffy, Consumption threshold used to investigate stability and ecological dominance in consumer-resource dynamics, Ecol. Model. 319 (2016), pp. 155–162.
- O.C. Collins and K.J. Duffy, Optimal control of maize foliar diseases using the plants population dynamics, Acta Agric. Scand B Soil Plant Sci. 66 (2016), pp. 20–26.
- P. Daszak, A.A. Cunningham, and A.D. Hyatt, Emerging infectious diseases of wildlife–threats to biodiversity and human health, Science 287 (2000), pp. 443–449.
- K.J. Duffy, Simulations to investigate animal movement effects on population, Nat. Resour. Model.24 (2001), pp. 48–60.
- K.J. Duffy and O.C Collins, Identifying stability conditions and Hopf bifurcations in a consumer resource model using a consumption threshold, Ecol. Complex. 28 (2016), pp. 212–217.
- M.A. Duffy, J.M. Housley, R.M. Penczykowski, C.E. Caceres, and S.E. Hall, Unhealthy herds: Indirect effects of predators enhance two drivers of disease spread, Funct. Ecol. 25(5) (2011), pp. 945–953.
- J. Dushoff, W. Huang, and C. Castillo-Chavez, Backwards bifurcations and catastrophe in simple models of fatal diseases, J. Math. Biol. 36(3) (1998), pp. 227–248.
- J.M. Fryxell, C. Packer, K. McCann, E.J. Solberg, and B-E. Sther, Resource management cycles and the sustainability of harvested wildlife populations, Science 328 (5980), (2010), pp. 903–906.
- S. Greenhalgh, A.P. Galvani, and J. Medlock, Disease elimination and re-emergence in differential-equation models, J. Theor. Biol. 387 (2015), pp. 174–180.
- A.B. Gumel, Causes of backward bifurcations in some epidemiological models, J. Math. Anal. Appl.395(1) (2012), pp. 355–365.
- M. Haque and E. Venturino, The role of transmissible diseases in the Holling-Tanner predator-prey model, Theor. Popul. Biol. 70(3) (2006), pp. 273–288.
- M.J. Hatcher, J.T. Dick, and A.M. Dunn, How parasites affect interactions between competitors and predators, Ecol. Lett. 9 (11), (2006), pp. 1253–1271.
- H.W. Hethcote, W. Wang, L. Han, and Z. Ma, A predatorprey model with infected prey, Theor. Popul. Biol. 66(3) (2004), pp. 259–268.
- F.M. Hilker and K. Schmitz, Disease-induced stabilization of predator-prey oscillations, J. Theor. Biol.255(3) (2008), pp. 299–306.
- R.D. Holt and M. Roy, Predation can increase the prevalence of infectious diseases, Am. Nat. 169(5) (2007), pp. 690–699.
- P.J. Hurtado, S.R. Hall, and S.P. Ellner, Infectious disease in consumer populations: Dynamic consequences of resource-mediated transmission and infectiousness, Theor. Ecol. 7 (2014), pp. 163–179.
- P.T. Johnson, D.L. Preston, J.T. Hoverman, and K.L. Richgels, Biodiversity decreases disease through predictable changes in host community competence, Nature 494 (2013), pp. 230–233.
- F. Keesing, L.K. Belden, P. Daszak, A. Dobson, C.D. Harvell, R.D. Holt, P. Hudson, A. Jolles, K.E.Jones, C.E. Mitchell, and S.S. Myers, Impacts of biodiversity on the emergence and transmission of infectious diseases, Nature 468 (2010), pp. 647–652.
- B.W. Kooi, G.A.K. van Voorn, and K. pada Das, Stabilisation and complex dynamics in a predator-prey model with predator suffering from an infectious disease, Ecol. Complex. 8 (2011), pp. 113–122.
- B.W. Kooi and E. Venturino, Ecoepidemic predator–prey model with feeding satiation, prey herd behavior and abandoned infected prey, Math. Biosci. 274 (2016), pp. 58–72.
- S. Liao and J. Wang, Stability analysis and application of a mathematical cholera model, Math. Biosci. Eng. 8 (2011), pp. 733–752.
- H. McCallum, N. Barlow, and J. Hone, How should pathogen transmission be modelled?, Trends Ecol. Evol. 16(6) (2001), pp. 295–300.
- Z. Mukandavire, S. Liao, J. Wang, H. Gaff, D.L. Smith, and J.G. Morris, Estimating the reproductive numbers for the 2008–2009 cholera outbreaks in Zimbabwe, Proc. R. Acad. Sci. 108 (21), (2011), pp. 8767–8772.
- T. Nakazawa, T. Yamanaka, and S. Urano, Model analysis for plant disease dynamics co-mediated by herbivory and herbivore-borne phytopathogens, Biol. Lett. 8 (2012), pp. 685–688.
- N. Chitnis, J.M. Cushing, J.M. Hyman, Bifurcation analysis of a mathematical model for malaria transmission. SIAM J. Appl. Math. 67(1) (2006), pp. 24–45.
- N. Owen-Smith, Functional heterogeneity in resources within landscapes and herbivore population dynamics, Landscape Ecol. 19 (2004), pp. 761–771.
- C. Packer, R.D. Holt, P.J. Hudson, K.D. Lafferty, and A.P. Dobson, Keeping the herds healthy and alert: Implications of predator control for infectious disease, Ecol. Lett. 6(9) (2003), pp. 797–802.
- K. Pada Das, K. Kundu, and J.A. Chattopadhyay, A predator-prey mathematical model with both the populations affected by diseases, Ecol. Complex. 8 (2011), pp. 68–80.
- A.H. Robbins, M.D. Borden, B.S. Windmiller, M. Niezgoda, L.C. Marcus, S.M. O'Brien, S.M. Kreindel, M.W. McGuill, A. DeMariaJr, C.E. Rupprecht, and S. Rowell, Prevention of the spread of rabies to wildlife by oral vaccination of raccoons in Massachusetts, J. Am. Vet. Med. Assoc. 213(10) (1998), pp. 1407–1412.
- R. Ross, The Prevention of Malaria, 2nd ed. John Murray, London, UK, 1911.
- J.J. Tewa, V.Y. Djeumen, and S. Bowong, Predator-prey model with Holling response function of type II and SIS infectious disease, Appl. Math. Model. 37 (2013), pp. 4825–4841.
- J.H. Tien and D.J.D. Earn, Multiple transmission pathways and disease dynamics in a waterborne pathogen model, Bull. Math. Biol. 72 (2010), pp. 1506–1533.
- P. Turchin, Complex Population Dynamics: A Theoretical/Empirical Synthesis, Princeton University Press, Princeton, NJ, 2003.
- P. van den Driessche and J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci. 180 (2002), pp. 29–48.
- E. Venturino, Ecoepidemiology: a more comprehensive view of population interactions, Math. Model Nat. Phenom. 11(1) (2016), pp. 49–90.