JAM Journal of Applied Mathematics 1687-0042 1110-757X Hindawi Publishing Corporation 10.1155/2015/275485 275485 Research Article A Time Scales Approach to Coinfection by Opportunistic Diseases Marvá Marcos 1 http://orcid.org/0000-0001-7215-5114 Venturino Ezio 2 Bravo de la Parra Rafael 1 Garira Winston 1 Departamento de Física y Matemáticas Universidad de Alcalá 28871 Alcalá de Henares Spain uah.es 2 Dipartimento di Matematica “Giuseppe Peano” Università di Torino Via Carlo Alberto 10 10123 Torino Italy unito.it 2015 21 4 2015 2015 23 10 2014 25 03 2015 21 4 2015 2015 Copyright © 2015 Marcos Marvá et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Traditional biomedical approaches treat diseases in isolation, but the importance of synergistic disease interactions is now recognized. As a first step we present and analyze a simple coinfection model for two diseases simultaneously affecting a population. The host population is affected by the primary disease, a long-term infection whose dynamics is described by a SIS model with demography, which facilitates individuals acquiring a second disease, secondary (or opportunistic) disease. The secondary disease is instead a short-term infection affecting only the primary infected individuals. Its dynamics is also represented by a SIS model with no demography. To distinguish between short- and long-term infection the complete model is written as a two-time-scale system. The primary disease acts at the slow time scale while the secondary disease does at the fast one, allowing a dimension reduction of the system and making its analysis tractable. We show that an opportunistic disease outbreak might change drastically the outcome of the primary epidemic process, although it does among the outcomes allowed by the primary disease. We have found situations in which either acting on the opportunistic disease transmission or recovery rates or controlling the susceptible and infected population size allows eradicating/promoting disease endemicity.

1. Introduction

Coinfection is the simultaneous infection of a host by multiple pathogen species. The global incidence of coinfection among humans is huge  and supposed to be more common than single infection. The interactions between pathogen species within their host can have either positive or negative effects on each other. The net effect of coinfection on human health is thought to be negative .

The case of positive parasite interactions falls into the concept of syndemic: aggregation of two or more diseases in a population in which there is some level of positive biological interaction that exacerbates the negative health effects of any or all of the diseases . From the point of view of prevention and treatment of disease it is the opposite case that is important, sometimes called counter syndemic: disease interactions that yield a lower whole effect than the sum effects of the individual diseases involved. An example of counter syndemic is that of human immunodeficiency virus (HIV) transiently suppressed during acute measles infections. A broadly extended syndemic involves tuberculosis (TB) and HIV . The World Health Organization  reports that people living with HIV are around 30 times more likely to develop TB than persons without HIV and also that TB is the most common occurring illness among people living with HIV. Other syndemics involving infectious diseases have been described in the literature: HIV and malaria syndemic ; the helminthic infections, malaria, and HIV/AIDS syndemic ; the pertussis, influenza, and tuberculosis syndemic ; and the HIV and sexually transmitted disease (STD) syndemic .

In this work we deal with a particular, but very common, type of coinfection. We consider the interactions of two diseases, the first one of the type called primary disease and the second one of the opportunistic disease type. Only relatively few pathogen species cause disease in otherwise healthy individuals . Those few are called primary pathogens. The diseases that they cause, primary diseases, are the result of their only activity within a healthy host. An opportunistic disease, on the other hand, is characterized  as a serious, usually progressive infection by a microorganism that has limited (or no) pathogenic capacity under ordinary circumstances, but which has been able to cause serious disease as a result of the predisposing effect of another disease or of its treatment.

The importance of opportunistic diseases for public health [2, 12, 13] is underrepresented in the mathematical modelling literature. A reason for that is that models of coinfection usually result in large dimensional systems which are difficult to be studied analytically. The main aim of this work is to settle a model describing the interaction between both, the primary and the secondary diseases. The model that we present in this work tries to capture the basic features of a coinfection model using for it the least possible number of variables. The dynamics of the primary disease is represented by means of a S I S model. All individuals affected by the primary infection are assumed to be susceptible of being infected by the opportunistic disease. As the dynamics of the opportunistic disease is also described in terms of a S I S model, we only distinguish three types of individuals in the population: individuals with no infection, susceptible, individuals infected by the primary disease but not by the opportunistic disease, primary infected, and individuals infected by both diseases, coinfected.

Specifically, we want to know whether or not the coinfection by a secondary disease produces epidemiological scenarios not allowed by the primary disease submodel. In the latter case, it is of interest to assess if coinfection has any influence on the actual outcome of the model, even if it is only among those allowed by the primary disease submodel. On the other hand, and in any case, we look for identifying mechanisms to modulate the epidemiological outcome.

A primary disease enabling secondary infections has typically a long illness period. It must produce a persistent alteration of the immune response which weakens the body’s ability to clear secondary diseases. On the other hand, a compromised immune system presents an opportunity that a secondary pathogen must rapidly take advantage of. As a simplified approximation of the general case we suppose that the primary disease is a long-term infection that evolves slowly compared to the opportunistic disease which has a rapid evolution and, thus, can be considered a short-term infection. This difference in the acting speed of both infections is reflected in our model in two different issues. Firstly, we assume that demography has an impact in the primary disease, due to its slow evolution, whereas it is negligible for the opportunistic disease which evolves in a short period of time. Secondly, the system of differential equations, in terms of which we express our model, possesses two time scales: the slow one encompassing the demography and the primary disease evolution and the fast one associated with the opportunistic disease evolution.

The inclusion of two time scales in the system has the advantage of allowing its reduction. The asymptotic behavior of the solutions of the initial three-dimensional system can be studied through a planar system. The reduction of the system is undertaken with the help of aggregation methods . The general aim of these methods is studying the relationships between a large class of complex systems, in which many variables are involved, and their corresponding reduced or aggregated systems, governed by a few global variables. The idea behind the reduction of the system in our model is considering the evolution of the secondary infection as instantaneous in relation to that of the primary one. Obviously this is but an approximation which, on the other hand, can be precisely treated with the help of the aggregation method. The steady state rapidly, almost instantaneously, reached by the opportunistic disease serves to merging in one single variable those variables corresponding to primary infected and coinfected individuals. The result is a S I S type model where the effect of the opportunistic disease is reflected in its parameters.

The model is presented in Section 2. In this section the reduction of the system is also included. Section 3 is devoted to the analysis of the reduced system. This analysis allows a discussion of the permanence of the population as well as of the influence of the final size of the opportunistic disease on the outcome of the primary epidemic. This discussion is the content of Section 4.

2. The Model

We build up in this section a model of coinfection that describes the interaction between two diseases, one of primary type whereas the second one is of opportunistic type. Only the individuals infected by the primary disease are susceptible of being infected by the opportunistic disease. Moreover, the interaction of both diseases occurs at different time scales, the evolution of the opportunistic disease being much faster than that of the primary one. The model is written in terms of a slow-fast ordinary differential equations model. After building the slow-fast model, the separation of time scales allows us to apply approximate aggregation techniques [14, 15] to get a smaller dimensional system. For the convenience of readers nonfamiliar with it, the reduction procedure is sketched in Section 2.4. The section finishes describing which kind of information about the slow-fast system can be retrieved from the reduced system.

2.1. The Primary Disease Submodel

The primary disease dynamics is described by a S I S model with demographic effects. In a S I S model individuals are divided into susceptible ( S ) and infected ( I ). The latter return to the susceptible class on recovery because the disease confers no immunity against reinfection . It is appropriate for most diseases transmitted by bacterial or helminth agents and most sexually transmitted diseases. Concerning transmission, there are two extreme traditional forms:  the density-dependent transmission (DDT) and the frequency-dependent transmission (FDT). In DDT the rate of contact between susceptible and infected individuals increases with host density while in FDT this rate of contact is independent of host density. The fact that the primary disease acts together with demography at the same time scale leads us to assume it does with density-dependent transmission. On the other hand, in the case of the opportunistic disease which turns out to evolve at a faster time scale we consider that it does with frequency-dependent transmission [19, 20].

We denote by γ the recovery rate and by β the constant transmission rate. The parameter μ describes the additional disease-induced mortality.

We consider demographic effects with only horizontal transmission of the disease. In many mathematical models, from a demographic point of view, the differences between susceptible and infected individuals are reduced to an additional disease-related death rate or disease-induced reduction in fecundity . However, there are experimental evidences of the influence of disease on host competitive abilities  which have already been introduced in ecoepidemiological models . We adopt this last approach. The intrinsic per capita fertility rate of uninfected individuals is given by r . The reduction on intrinsic per capita fertility rate of infected individuals is represented by the parameter a ( 0,1 ) . The natural death rate is denoted by m . The effects of intraspecific competition reducing population growth are introduced in the model by means of parameters c S S , c S I , c I S , and c I I . To be precise, the parameters c S S and c I I represent intraclass competition between susceptible and infected individuals, respectively, whereas the parameters c S I and c I S introduce the interclass impact of infected on susceptible individuals and of susceptible on infected individuals, respectively.

The primary disease submodel is given by the equations (1) d S d t = r S + a r I - m S - c S S S + c S I I S - β S I + γ I , d I d t = - m I density independent growth - c I S S + c I I I I competition + β S I - γ I - μ I transmission , recovery , and disease mortality . As mentioned in the introduction, to our knowledge, the primary disease submodel (1) has not been previously analyzed. However, we postpone its analysis to Section 3, once we have described the full model and the aforementioned reduction process.

2.2. The Opportunistic Disease Submodel

The opportunistic disease spreads only through the individuals infected by the primary disease. We consider that the opportunistic disease dynamics is also described by a S I S model. Individuals infected by the primary disease are further classified into those not infected by the opportunistic disease ( U ), primary infected, and those infected by both diseases ( V ), coinfected.

The fast evolution of the opportunistic disease, compared to primary disease and demography, suggests not including demographic effects and choosing the frequency-dependent transmission form. Let λ and δ be, respectively, the constant transmission and recovery rates.

The opportunistic disease submodel is represented by the equations (2) d U d τ = - λ U V U + V + δ V , d V d τ = λ U V U + V - δ V . We use τ to denote the time variable for the fast time scale. It is related to variable time t in system (1) as t = ε τ , where ε is a small positive constant representing the ratio between time scales.

2.3. The Full Two-Time-Scale Model

Finally, we construct the model encompassing both diseases. It has the form of a system with three state variables: susceptible S , primary infected U , and coinfected V individuals. It is a system with two time scales that is expressed in terms of the fast time variable τ . The terms associated with the slow time scale, demography and primary disease dynamics, appear multiplied by ε in (3). The fast part of system (3), the opportunistic disease dynamics, coincides with system (2).

In the slow part of system (3) we have to define different rates for primary infected and coinfected individuals. We denote by β U and β V the constant primary disease transmission rates due to primary infected and coinfected individuals, respectively. We assume that there is no direct connection between the susceptible and coinfected stages. A susceptible individual must first acquire the primary disease and later be infected by the opportunistic one. On the other hand, a coinfected individual must first recover from the opportunistic disease and then, being just primary infected, can also recover from the primary one. The primary disease recovery rate is still denoted by γ . Parameters μ U and μ V describe the additional primary disease-induced mortality in primary infected and coinfected individuals, respectively.

Concerning the part of demography, we keep the same intrinsic per capita fertility rate of uninfected individuals r and the individuals natural death rate m as in system (1). We include different coefficients of reduction on intrinsic per capita fertility rate for primary infected and coinfected individuals: a U and a V . We assume them to verify 0 < a V < a U < 1 supposing that coinfected individuals participate in reproduction though at a smaller rate. To distinguish the effects of intraspecific competition among the three stages we need to introduce nine parameters c S S , c S U , c S V , c U S , c U U , c U V , c V S , c V U , and c V V . They represent the competition, either intraclass or interclass, between the two stages in each of the nine different interaction pairs.

The complete two-time-scale system reads as follows: (3) d S d τ = ε r S + a U r U + a V r V - m S - c S S S + c S U U + c S V V S - β U S U - β V S V + γ U , d U d τ = - λ U V U + V + δ V + ε - m U - c U S S + c U U U + c U V V U + β U S U + β V S V - γ U - μ U U , d V d τ = λ U V U + V - δ V + ε - m V - c V S S + c V U U + c V V V V - μ V V .

2.4. Reduction of the Model

In this section we take advantage of the two time scales to reduce the dimension of the complete system (3). In the next section, as a consequence of this reduction, we perform the analysis of the model by means of a planar system. The reduction follows the technique called approximate aggregation method [14, 15]. The first step is writing the system in the so-called slow-fast form. This is easily done in system (3) using the change of variables ( S , U , V ) ( S , I , V ) , where I = U + V represents all the infected individuals, both primary infected and coinfected. Consider (4) d S d τ = ε r S + a U r I - V + a V r V - m S - c S S S + c S U I - V + c S V V S - β U S I - V - β V S V + γ I - V , d I d τ = ε - m I - c U S S + c U U I - V + c U V V I - V + β U S I - V + β V S V - γ I - V - μ U I - V - c V S S + c V U I - V + c V V V V - μ V V , d V d τ = λ I - V V I - δ V + ε - m V - c V S S + c V U I - V + c V V V V - μ V V . The key point of the new form of system (3) is making it visible that variables S and I are slow (the right-hand side terms of their equations have ε as a factor) in the sense that they almost do not change at the fast time scale. The fast dynamics is concentrated in the first terms without ε in the equation for V . The approximation that the aggregation method proposes consists in separating both dynamics. Firstly, the nonslow variables are calculated in terms of the slow ones by assuming that they are the equilibria (called fast equilibria) determining the long-term behaviour of the fast dynamics. Secondly, these obtained values of the nonslow variables are substituted into the equations of the slow ones yielding a reduced system for the latter. In this reduced or aggregated system the fast dynamics is summarized in its parameters. In the particular case of system (4), the only nonslow variable is V and the fast dynamics reduces to the equation (5) d V d τ = λ I - V V I - δ V . Assuming the slow variable I to be constant, the analysis of (5) for positive values of V ( 0 ) gives (6) lim τ V τ = ν I = 0 if δ λ 1 - δ λ I if δ < λ , which corresponds to the results of a classical S I S model without demography and frequency-dependent transmission . If the recovery rate is larger than the transmission rate, the disease disappears since the number of coinfected individuals tends rapidly to zero. On the other hand, if the recovery rate is smaller than the transmission rate, the disease becomes endemic with a stable fraction ν = 1 - δ / λ of the infected population I remaining coinfected.

The fast equilibria V = ν I found in (5) are the values to be substituted into the equations for the slow variables S and I to obtain the following reduced system: (7) d S d t = r S + a - r I - m S - c S S S + c - S I I S - β - S I + γ - I , d I d t = - m I - c - I S S + c - I I I I + β - S I - γ - I - μ - I , which has the same form as the primary disease submodel (1). In its parameters the effect of fast dynamics, the opportunistic disease, is implicit through ν (8) a - = 1 - ν a U + ν a V , c - S I = 1 - ν c S U + ν c S V , c - I S = 1 - ν c U S + ν c V S , c - I I = 1 - ν 2 c U U + 1 - ν ν c U V + ν 1 - ν c V U + ν 2 c V V , β - = 1 - ν β U + ν β V , γ - = 1 - ν γ , μ - = 1 - ν μ U + ν μ V . The reduced system (7) is useful to analyze the asymptotic behaviour of the solutions of the complete system (3) [14, 15]. In particular, the existence of a hyperbolic asymptotically stable equilibria ( S , I ) of system (7) ensures the existence, for ε small enough, of an equilibria of system (3) with the same characteristics and a form very close to ( S , ( 1 - ν ) I , ν I ) . In the next section we carry out the analysis of the stability of equilibria of system (7) obtaining thus the corresponding results for the complete model (3).

Note that the reduced system (7) and the primary disease submodel (1) are the same, the only difference being the values of the respective coefficients. Indeed, when δ > λ the opportunistic disease cannot invade the population and, in this case, the coefficients (8) of systems (7) and (1) are exactly the same.

3. Analysis of the Reduced System

We proceed in this section to analyze the reduced system (7).

We first note that E 0 = ( 0,0 ) is an equilibrium point, the positive S semiaxis, { ( S , 0 ) : S > 0 } , is invariant, and on the positive I semiaxis, { ( 0 , I ) : I > 0 } , the vector field associated with system (7) points to the interior of the positive quadrant. We then have that the closed positive quadrant R + 2 = { ( S , I ) : S 0 , I 0 } is positively invariant.

In the next result we prove that, as expected, if the susceptible fertility rate r is not strictly larger than the natural death rate m the population gets extinct.

Proposition 1.

If r m then any solution ( S ( t ) , I ( t ) ) of system (7) with nonnegative initial conditions ( S ( 0 ) , I ( 0 ) ) tends to E 0 .

Proof.

Let us call W = S + I . Summing up both equations in system (7) and now choosing c = min { c S S , c - S I , c - I S , c - I I } we have d W / d t - c W 2 that, by integration, yields (9) 0 W t W ( 0 ) 1 + W 0 c t t 0 since R + 2 is positively invariant.

Henceforth, we assume that r > m . This assumption prevents the population from extinction. The linearization of system (7) at the equilibrium E 0 has the matrix (10) r - m a - r + γ - 0 - m + γ - + μ - with one positive and one negative eigenvalues. The unstable manifold of E 0 , associated with r - m , is included in the S axis, while the stable manifold, associated with - ( m + γ - + μ - ) , is tangent at 0 to the eigenvector a - r + γ - , - ( r + γ - + μ - ) and lies completely outside the interior of the positive quadrant.

Assuming r > m the only nonnegative solution tending to 0 is E 0 itself. We prove next that all nonnegative solutions are forward bounded.

Proposition 2.

Let r > m . If ( S ( t ) , I ( t ) ) is any solution of system (7) with nonnegative initial conditions ( S ( 0 ) , I ( 0 ) ) then it is bounded on [ 0 , ) .

Proof.

Calling W = S + I and letting c = min { c S S , c - S I , c - I S , c - I I } we have (11) d W d t + r - m W 2 r - m W - c W 2 . Function g ( W ) = 2 ( r - m ) W - c W 2 attains its maximum on [ 0 , ) at W = ( r - m ) / c , so that (12) d W d t + r - m W r - m 2 c . Multiplying both sides of (12) by e ( r - m ) t and rearranging terms yield (13) d d t e r - m t W e ( r - m ) t r - m 2 c , which implies, integrating (13) on [ 0 , t ] , (14) e r - m t W t - W 0 r - m 2 c r - m e r - m t - 1 . Rearranging terms in expression (14) leads to (15) W ( t ) W ( 0 ) e - r - m t + r - m c 1 - e - r - m t and, finally, we have that W ( t ) max { W ( 0 ) , ( r - m ) / c } for every t [ 0 , ) .

In addition to the trivial equilibrium E 0 , system (7) possesses a disease-free equilibrium E 1 = ( S 1 , 0 ) , where (16) S 1 = r - m c S S that represents the stable size of the population in case of no infection. The growth of the population in the absence of infections is logistic and S 1 represents its carrying capacity.

Proposition 3.

Let r > m . The equilibrium point E 1 = ( S 1 , 0 ) of system (7) verifies

if   S 1 ( β - - c - I S ) > m + γ - + μ - then E 1 is a saddle point, the stable manifold of which coincides with the positive S semiaxis;

if   S 1 ( β - - c - I S ) < m + γ - + μ - then E 1 is locally asymptotically stable;

if   β - - c - I S 0 then the basin of attraction of E 1 includes R + 2 - { 0 } .

Proof.

To prove the two first items it suffices to calculate the matrix of the linearization of system (7) at E 1 : (17) - r - m a - r + γ - - S 1 c - S I + β - 0 S 1 β - - c - I S - m + γ - + μ - . One of the eigenvalues, - ( r - m ) , is negative while stability depends on the other one, S 1 ( β - - c - I S ) - ( m + γ - + μ - ) , being positive or negative: E 1 is a saddle or (locally) asymptotically stable, respectively.

To prove the last assertion we first note that there exist no interior equilibria because the right-hand side of the I equation is always negative for positive S and I . Now the Poincaré-Bendixson theorem implies that there is no closed orbit in the interior of the positive quadrant and therefore that all positive solutions that are forward bounded must tend to the unique nonnegative equilibrium point, E 1 .

Up to now we have obtained the condition of nonextinction of the population, r > m , and a sufficient condition, β - - c - I S 0 , for a disease-free scenario in the long term. This last condition says that if the competition coefficient c - I S , representing the impact of susceptible on infected individuals, is larger than the transmission rate β - then the infection disappears independently of the initial conditions. More significant cases exist when the simple competitive pressure of susceptible on infected individuals is not enough to compensate transmission.

From now on we are also assuming that β - - c - I S > 0 . In this case the conditions of local stability of the equilibrium E 1 can be expressed in terms of the parameter (18) A - = m + γ - + μ - β - - c - I S . Thus, Proposition 3 can be restated as follows: for r > m and β - > c - I S , if S 1 > A - or S 1 < A - , then the equilibrium E 1 is a saddle point or locally asymptotically stable, respectively. The parameter A - represents a threshold population size allowing or not the increase of the infection when it is rare. If the susceptible population is close to its carrying capacity, S 1 , a few infected individuals are able to spread the disease if the number of susceptible individuals is large enough, S 1 > A - . On the other hand, the infection disappears if the susceptible population is under the threshold A - .

In the next results we search for conditions ensuring the endemicity of the infection. To express them in a simpler form we define another parameter (19) B - = a - r + γ - c - S I + β - , which can be interpreted through the terms depending on I in the first equation of system (7). This equation can be written in the following form: (20) d S d t = r - m S - c S S S 2 + a - r + γ - - c - S I + β - S I , where we note that depending on whether a - r + γ - - ( c - S I + β - ) S ( t ) is positive or negative, the existence of infected individuals makes the susceptible growth rate increase or decrease, respectively. The size of the susceptible population determines if the infection has a positive or a negative effect on its growing. If S ( t ) < B - , the more the infected individuals the larger the susceptible population growth rate, while S ( t ) > B - yields a larger decrease of the susceptible population growth rate whenever there is a larger infected population.

Using parameters S 1 , A - , and B - , system (7) can be expressed as follows: (21) d S d t = c S S S S 1 - S + c - S I + β - B - - S I , d I d t = β - - c - I S S - A - I - c - I I I 2 .

The equation of the S -nullcline of system (21), for S 1 B - , is (22) I = Φ S = c S S c - S I + β - · S ( S 1 - S ) S - B - . We are interested in the part included in the positive quadrant. This is for S 1 < B - an increasing branch going from the point E 1 = ( S 1 , 0 ) to the asymptote S = B - (see Figure 1, left panel) and for S 1 > B - a decreasing branch going from the asymptote S = B - to the point E 1 (see Figure 1, right panel). In the case S 1 = B - the S -nullcline in the positive quadrant reduces to the line S = B - (see Figure 1, center panel).

Possible profiles of the S -nullcline I = Φ ( S ) and the I -nullcline I = Ψ ( S ) and the S -nullcline I = Φ ( S ) (in black, their intersection with the positive cone).

The I -nullcline of system (21) in the open positive quadrant is the line (23) I = Ψ S = β - - c - I S c - I I S - A - . It is immediate to prove that if A - max { S 1 , B - } there are no interior equilibria of the system (see panels in Figure 1) and thus, applying again the Poincaré-Bendixson theorem, we get that all positive solutions tend to E 1 . We gather these results in the next proposition.

Proposition 4.

Let r > m and β - - c - I S > 0 . If A - max { S 1 , B - } , then system (7) possesses a unique nonnegative equilibrium point E 1 = ( S 1 , 0 ) that is asymptotically stable and attracts every positive solution.

Condition A - max { S 1 , B - } tell us, on the one hand, that the infection cannot invade due to A - S 1 and, on the other hand, that infected individuals cannot help in attaining the invasion threshold because A - B - . The consequence is that infection disappears.

There are two situations for the infection to become endemic. The first one is allowing invasion, that is, S 1 > A - , that is treated in Proposition 5. The second one does not allow infection invasion for a low number of infected individuals, S 1 < A - , but larger numbers of infected individuals might help the susceptible population growing, A - B - , so as to maintain this latter over the invasion threshold. In Proposition 6 are detailed sufficient conditions to meet this second situation.

Proposition 5.

Let r > m and β - - c - I S > 0 . If S 1 > A - then system (7) possesses a unique interior equilibrium point E + = ( S + , I + ) that is locally asymptotically stable. If, in addition, S 1 < B - , then E + attracts every positive solution.

Proof.

The assumptions on parameters yield the existence of unique interior equilibrium (see Figure 2). The asymptotic stability that follows can be proved by linearization.

Note that the condition S 1 > A - might not ensure that all positive solutions tend to the interior equilibrium. Due to the Poincaré-Bendixson theorem, it might happen that some of these solutions tend to a limit cycle included in the open positive quadrant surrounding the equilibrium.

Condition S 1 < B - excludes the existence of any limit cycle because the region A = ( S , I ) ; Φ ( S ) S Ψ ( S ) , 0 I is an invariant (“trapping”) subregion A R + 2 such that E + A , the boundary of A , and A R + 2 (see Figure 2). That is, any orbit surrounding E + must enter A but cannot leave from there.

The invariant region A = ( S , I ) ; Φ ( S ) S Ψ ( S ) , 0 I mentioned in the proof of Proposition 5.

In any case, what condition S 1 > A - ensures is the endemicity of the infection. In the Proposition 6 we state conditions leading the population to the disease-free state or towards conditional endemicity related with a bistable scenario. By conditional endemicity we mean that the outcome of the model can be either disease-free (see Figure 3(a)) or with an endemic disease depending on the initial amount of susceptible and infected individuals; see Figure 3(c).

Related to Proposition 6, possible configurations of the nullclines when S 1 < A - < B - . (a) S 1 is global attractor. (b) S 1 is an attractor. (c) The bistable case: S 1 and E + 2 are locally asymptotically stable while E + 1 is a saddle.

Indeed, we introduce R that appears later in the corresponding proof and drives the epidemic outcome. This quantity depends on the parameters of the model and allows discriminating whether conditional disease endemicity is allowed or cannot occur.

Proposition 6.

Let r > m and β - - c - I S > 0 . If S 1 < A - < B - then system (7) possesses the asymptotically stable equilibrium E 1 = ( S 1 , 0 ) . Furthermore, let us define (24) R = c - I I r - m + a - r + γ - β - - c - I S + c - S I + β - μ - + γ - 2 c S S c - I I + c - S I + β - β - - c - I S a - r + γ - μ - + γ - . We have:

if R < 1 then there is no interior equilibrium point and the basin of attraction of E 1 includes R + 2 0 ;

if R = 1 then there is only one interior equilibrium point E + = ( S + , I + ) that is unstable; the equilibrium E 1 = ( S 1 , 0 ) attracts every solution with initial values in the interior of   R + 2 E + ;

if R > 1 then there are two interior equilibrium points E + 1 = ( S + 1 , I + 1 ) and E + 2 = ( S + 2 , I + 2 ) , with S + 1 < S + 2 and I + 1 < I + 2 . E + 1 is a saddle and E + 2 is locally asymptotically stable.

Proof.

Note that the asymptotic stability of E 1 follows directly from Proposition 3 since S 1 < A - .

Next, we focus on showing the relation between R and the existence of equilibrium points of system (7). Equating the nullclines I = Φ ( S ) and I = Ψ ( S ) of system (7) yields (25) c S S c - S I + β - · S S 1 - S S - B - = β - - c - I S c - I I S - A - . Keeping in mind the definition of A - and B - , the previous expression is equivalent to (26) 0 = c S S c - I I + c - S I + β - β - - c - I S S 2 - c - I I r - m + a - r + γ - β - - c - I S + c - S I + β - μ - + γ - S + a - r + γ - μ - + γ - . Now, direct calculations lead to the fact that R being smaller than, equal to, or larger than 1, is equivalent to the discriminant of (26) being negative, zero, or positive. This yields the number of equilibrium points.

Note that when R < 1 there is no interior equilibrium point and a direct application of the Poincaré-Bendixson theorem yields statement 1.

Concerning statement ( 2 ) , direct calculations show that when R = 1 the equilibrium E is not hyperbolic so that we cannot use the linearization criterion. Note that the region B = ( S , I ) ; Ψ ( S ) S Φ ( S ) , 0 I is an invariant region such that E + B and B R + 2 . Indeed, any solution with initial values in the interior of B converges to E 0 since Ψ ( S ) < S < Φ ( S ) , which implies that E + is unstable. The nonexistence of periodic orbits can be argued as done in the proof of Proposition 5.

We now assume R > 1 and analyze the stability of the equilibria E = ( S , I ) by means of the well known trace-determinant criterion. Let us consider the Jacobian matrix of the flow of system (7) at the equilibrium point which, taking into account the fact that Φ ( S ) = I = Ψ ( S ) , simplifies to (27) J = - c S S S - a - r + γ - I S a - r + γ - - C - S I + β - S β - - c - I S I - c - I I I . This immediately yields tr F < 0 . Furthermore, a direct calculation leads to (28) det J I = c S S c - I I + C - S I + β - β - - c - I S S 2 + a - r + γ - c - I I I - β - - c - I S S · S - 1 . Using again the fact that I = Ψ ( S ) , if and only if c - I I I = ( β - - c - I S ) S - ( μ - + γ - ) , we have (29) det J I = k c - I I + C - S I + β - β - - c - I S S 2 - a - r + γ - μ - + γ - S , so that det J > 0 is equivalent to (30) S > a - r + γ - μ - + γ - k c - I I + C - S I + β - β - - c - I S , which entails local stability. On the contrary, E is unstable if (31) S < a - r + γ - μ - + γ - k c - I I + C - S I + β - β - - c - I S . The S component of the equilibrium points E + 1 and E + 2 can be explicitly calculated from (26). Direct calculations show that S + 1 fulfills condition (31) while condition (30) holds for S + 2 .

4. Discussion

We have set up a model aimed at ascertaining the impact of an opportunistic disease outbreak in a population already affected by a primary disease by assuming that both diseases evolve within different time scales. For the discussion of results, let us remember the two main aims stated in the Introduction. On the one hand, we wanted to know whether the coinfection by a secondary disease produces epidemiological scenarios not allowed by the primary disease submodel or not. In the latter case, it is of interest to determine if coinfection has any influence on the actual outcome of the model, even if just among those allowed by the primary disease submodel.

The answer to the first question is negative, as we have pointed out at the end of Section 2. Thus, the catalog of possible qualitative epidemic behaviors remains unchanged by the influence of a secondary disease under the assumptions considered here. We can restate this fact by saying that there is neither functional nor dynamical emergence .

Nevertheless, the effect of the opportunistic disease must be taken into account. In Section 3 we have found that A - , B - , and S 1 , as well as R , are key parameters to describe the outcome of the model. All of them, but S 1 , depend on ν , the fraction of coinfected individuals which, in turn, depends on λ / δ , the ratio of the parameters describing the opportunistic disease dynamics. It means that the opportunistic disease can be decisive in the long-term behavior of the slow-fast model. Therefore, the interest lies on how A - , B - , and R vary with the quotient δ / λ . Unfortunately, such a dependence is, in general, not simple (just see (5) and (8)) and we resort to numerical tools to illustrate the effect of varying δ / λ in the outcome of the model. Figure 4 displays the different outcomes of the aggregated model as function of δ and λ : in yellow, values of δ and λ leading to the disease-free scenario; in orange, values leading to an endemic primary infection outcome; in gray, conditional coinfection, meaning values leading to either disease-free or endemic coinfection scenario, depending on the initial amount of susceptible and infected individuals; in red, values leading to disease endemic coinfection outcome.

Possible epidemic outcomes. Yellow: disease-free. Orange: endemic primary infection. Gray: disease-free or endemic coinfection depending on initial values. Red: endemic coinfection. (a) Parameter values: a U = 0 . 9 , a V = 0 . 7 , r = 26 , m = 12 , μ U = 0 . 3 , μ V = 0 . 5 , c S S = 3 . 8 , c S U = 0 . 5 , c S V = 0 . 5 , c U S = 2 . 6 , c U U = 0 . 1 , c U V = 1 , c V S = 0 . 5 , c V U = 4 , c V V = 4 , β U = 4 , β V = 8 , γ = 0 . 2 . The parameter values in (b) are the same as before but m = 17 , c S S = 2 . 8 , β V = 4 .

In Figure 4(b), the epidemiological outcome changes from disease-free to endemic coinfection as the ratio δ / λ increases and crosses the threshold δ / λ = 1 . Instead, in Figure 4(a), note that the disease-free region (in the parameters space) overlaps the region λ > δ (above the dotted line) where the opportunistic disease would be able to invade if there were primary infected individuals in the population. As the ratio δ / λ increases, the epidemiological outcome changes from disease-free to conditional coinfection and a further increase leads to endemic coinfection. The only difference between the parameter values used in each figure is on c S S and β V . And this fact leads us to another interesting finding: there is a delicate interplay between competition coefficients and infection parameters, captured by the definition of A - , B - , S 1 , and R , which must be taken into account. Although we could not derive general results of such an interdependence, we have shown that it must be taken into account.

Summing up, both the irruption of an opportunistic disease and the competitive pressure of individuals being in different epidemiological state may affect the evolution of the primary disease outbreak. The effect can be determined by means of the parameters A - , B - , and R on δ / λ . And it leads us to the second objective of this work.

Related to our second objective, our results point out two different kinds of mechanisms to modulate the outcome of the model, each of them feasible within certain ranges of the parameter values.

On the one hand, having control on parameters δ and λ may allow certain leeway to reverse/promote epidemic outbreaks or infection/coinfection scenarios. Indeed, once the actual parameter values of the model are known one can compute (the equivalent of) Figure 4 and get enhanced comprehension on the epidemiological context as well as ascertaining the effect on the epidemic outbreak of changes on δ or λ . In this sense, it is interesting to note that any action or measure taken to modify the secondary infection recovery rate δ or transmission rate λ such that δ / λ remains constant is completely ineffective. In addition, controlling individuals competitive pressure may be relevant for the epidemiological outcome.

On the other hand, the results in Propositions 3 and 6 suggest that acting on the susceptible/infected individuals population size in order to keep the population above/below certain threshold allows having control on epidemic outbreaks. In particular, according to Proposition 3, A - is a susceptible population size threshold allowing or not the increase of the infection when it is rare. Therefore, introducing/culling (removing) susceptible individuals to bring the population above/below this threshold may certainly modify the outcome. Besides, under the hypotheses of Proposition 6, we show that when R > 1 whether the disease establishes itself or not depends on the initial amount of susceptible and infected individuals. From a mathematical point of view, this scenario is characterized by the fact that the disease-free equilibrium E 1 and an endemic disease (interior) equilibrium E + 2 coexist and are locally asymptotically stable. There is also a saddle node interior equilibrium E + 1 . The stable manifold of E + 1 separates the basins of attraction of the disease-free and the endemic disease steady states. This stable manifold cannot be calculated straightforward but can be computed using, for instance, the results in [24, 25].

A final comment has to do with the selection of the transmission form of the opportunistic disease. Preliminary calculations show that considering DDT instead of FDT leads to equivalent results. This means that even if the nullclines are different, the possible outcomes (say the dynamical scenarios) of the corresponding aggregated model are the same.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

M. Marvá and R. Bravo de la Parra are partially supported by Ministerio de Ciencia e Innovación (Spain), Projects MTM2011-24321 and MTM2011-25238. E. Venturino is partially supported by the Project “Metodi numerici in teoria delle popolazioni” of the Dipartimento di Matematica “Giuseppe Peano.”

Cox F. E. G. Concomitant infections, parasites and immune responses Parasitology 2001 122 supplement S23 S38 10.1017/s003118200001698x 2-s2.0-0034966996 Griffiths E. C. Pedersen A. B. P. Fenton A. Petchey O. L. The nature and consequences of coinfection in humans Journal of Infection 2011 63 3 200 206 10.1016/j.jinf.2011.06.005 2-s2.0-80052266824 Merrill S. Introducing Syndemics: A Critical Systems Approach to Public and Community Health 2009 Wiley Kwan C. Ernst J. D. HIV and tuberculosis: a deadly human syndemic Clinical Microbiology Reviews 2011 24 2 351 376 10.1128/cmr.00042-10 2-s2.0-79955059280 World Health Organization WHO Report 2011: Global Tuberculosis Control 2011 Abu-Raddad L. J. Patnaik P. Kublin J. G. Dual infection with HIV and malaria fuels the spread of both diseases in Sub-Saharan Africa Science 2006 314 5805 1603 1606 10.1126/science.1132338 2-s2.0-33845386531 Ivan E. Crowther N. J. Mutimura E. Osuwat L. O. Janssen S. Grobusch M. P. Helminthic infections rates and malaria in HIV-infected pregnant women on anti-retroviral therapy in Rwanda PLoS Neglected Tropical Diseases 2013 7 8 e2380 10.1371/journal.pntd.0002380 2-s2.0-84883325650 Herring D. A. Sattenspiel L. Social contexts, syndemics, and infectious disease in North Aboriginal populations American Journal of Human Biology 2007 19 2 190 202 10.1002/ajhb.20618 2-s2.0-33947410848 Eaton D. K. Lowry R. Brener N. D. Kann L. Romero L. Wechsler H. Trends in human immunodeficiency virus- and sexually transmitted disease-related risk behaviors among U.S. high school students, 1991–2009 American Journal of Preventive Medicine 2011 40 4 427 433 10.1016/j.amepre.2010.12.010 2-s2.0-79952561416 Baron S. Medical Microbiology 1996 4th Galveston, Tex, USA The University of Texas Medical Branch at Galveston Symmers W. S. Opportunistic infections. The concept of ‘opportunistic infections’ Proceedings of the Royal Society of Medicine 1965 58 341 346 Koster F. T. Curlin G. C. Aziz K. N. A. Haque A. Synergistic impact of measles and diarrhoea on nutrition and mortality in Bangladesh Bulletin of the World Health Organization 1981 59 6 901 908 2-s2.0-0019843231 Zlamy M. Kofler S. Orth D. Würzner R. Heinz-Erian P. Streng A. Prelog M. The impact of Rotavirus mass vaccination on hospitalization rates, nosocomial Rotavirus gastroenteritis and secondary blood stream infections BMC Infectious Diseases 2013 13 1, article 112 10.1186/1471-2334-13-112 2-s2.0-84874397445 Auger P. Bravo de la Parra R. Poggiale J.-C. Sánchez E. Nguyen-Huu T. Magal P. Ruan S. Aggregation of variables and applications to population dynamics Structured Population Models in Biology and Epidemiology 2008 1936 Berlin, Germany Springer 209 263 Lecture Notes in Mathematics 10.1007/978-3-540-78273-5_5 MR2433578 Auger P. Poggiale J. C. Sánchez E. A review on spatial aggregation methods involving several time scales Ecological Complexity 2012 10 12 25 10.1016/j.ecocom.2011.09.001 2-s2.0-84860456660 Marvá M. Bravo de la Parra R. Poggiale J. C. Approximate aggregation of a two time scales periodic multi-strain SIS epidemic model: A patchy environment with fast migrations Ecological Complexity 2012 10 1 34 41 2-s2.0-84860452464 10.1016/j.ecocom.2011.09.002 Brauer F. Castillo-Chávez C. Mathematical Models in Population Biology and Epidemiology 2001 40 New York, NY, USA Springer Texts in Applied Mathematics 10.1007/978-1-4757-3516-1 MR1822695 Kang Y. Castillo-Chavez C. Dynamics of SI models with both horizontal and vertical transmissions as well as Allee effects Mathematical Biosciences 2014 248 97 116 10.1016/j.mbs.2013.12.006 MR3162646 2-s2.0-84893475422 McCallum H. Barlow N. Hone J. How should pathogen transmission be modelled? Trends in Ecology & Evolution 2001 16 6 295 300 10.1016/s0169-5347(01)02144-9 2-s2.0-0035358852 Begon M. Bennett M. Bowers R. G. French N. P. Hazel S. M. Turner J. A clarification of transmission terms in host-microparasite models: numbers, densities and areas Epidemiology & Infection 2002 129 1 147 153 10.1017/s0950268802007148 2-s2.0-0036677204 Brauer F. Brauer F. van den Driessche P. Wu J. Compartmental models in epidemiology Mathematical Epidemiology 2008 1945 Berlin, Germany Springer 19 79 Lecture Notes in Mathematics 10.1007/978-3-540-78911-6_2 MR2428372 2-s2.0-42249097285 Bedhomme S. Agnew P. Vital Y. Sidobre C. Michalakis Y. Prevalence-dependent costs of parasite virulence PLoS Biology 2005 3 8, article e262 10.1371/journal.pbio.0030262 2-s2.0-23944451633 Sieber M. Malchow H. Hilker F. M. Disease-induced modification of prey competition in eco-epidemiological models Ecological Complexity 2014 18 74 82 10.1016/j.ecocom.2013.06.002 2-s2.0-84900424129 Cavoretto R. Chaudhuri S. de Rossi A. Menduni E. Moretti F. Rodi M. C. Venturino E. Simos T. Psihoyios G. Tsitouras C. Anastassi Z. Approximation of dynamical system's separatrix curves 1389 International Conference on Numerical Analysis and Applied Mathematics (ICNAAM '11) September 2011 Halkidiki, Greece 1220 1223 AIP Conference Proceedings 10.1063/1.3637836 Cavoretto R. de Rossi A. Perracchione E. Venturino E. Reliable approximation of separatrix manifolds in competition models with safety niches International Journal of Computer Mathematics 2014 10.1080/00207160.2013.867955 2-s2.0-84896574146