Journal of Applied Mathematics

Volume 2015 (2015), Article ID 275485, 10 pages

http://dx.doi.org/10.1155/2015/275485

## A Time Scales Approach to Coinfection by Opportunistic Diseases

^{1}Departamento de Física y Matemáticas, Universidad de Alcalá, 28871 Alcalá de Henares, Spain^{2}Dipartimento di Matematica “Giuseppe Peano”, Università di Torino, Via Carlo Alberto 10, 10123 Torino, Italy

Received 23 October 2014; Accepted 25 March 2015

Academic Editor: Winston Garira

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.

#### Abstract

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 [1] 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 [2].

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 [3]. 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 [4]. The World Health Organization [5] 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 [6]; the helminthic infections, malaria, and HIV/AIDS syndemic [7]; the pertussis, influenza, and tuberculosis syndemic [8]; and the HIV and sexually transmitted disease (STD) syndemic [9].

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 [10]. 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 [11] 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 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 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 [14–16]. 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 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 model with demographic effects. In a model individuals are divided into susceptible () and infected (). The latter return to the susceptible class on recovery because the disease confers no immunity against reinfection [17]. 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: [18] 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 [21]. However, there are experimental evidences of the influence of disease on host competitive abilities [22] which have already been introduced in ecoepidemiological models [23]. We adopt this last approach. The intrinsic per capita fertility rate of uninfected individuals is given by . The reduction on intrinsic per capita fertility rate of infected individuals is represented by the parameter . The natural death rate is denoted by . The effects of intraspecific competition reducing population growth are introduced in the model by means of parameters , , , and . To be precise, the parameters and represent intraclass competition between susceptible and infected individuals, respectively, whereas the parameters and introduce the interclass impact of infected on susceptible individuals and of susceptible on infected individuals, respectively.

The primary disease submodel is given by the equationsAs 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 model. Individuals infected by the primary disease are further classified into those not infected by the opportunistic disease (), primary infected, and those infected by both diseases (), 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 equationsWe use to denote the time variable for the fast time scale. It is related to variable time in system (1) as , 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 , primary infected , and coinfected 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 and 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 and 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 and the individuals natural death rate as in system (1). We include different coefficients of reduction on intrinsic per capita fertility rate for primary infected and coinfected individuals: and . We assume them to verify 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 , , , , , , , , and . 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:

##### 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 , where represents all the infected individuals, both primary infected and coinfected. ConsiderThe key point of the new form of system (3) is making it visible that variables and 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 . 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 and the fast dynamics reduces to the equationAssuming the slow variable to be constant, the analysis of (5) for positive values of giveswhich corresponds to the results of a classical model without demography and frequency-dependent transmission [17]. 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 of the infected population remaining coinfected.

The fast equilibria found in (5) are the values to be substituted into the equations for the slow variables and to obtain the following reduced system: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 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 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 . 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 is an equilibrium point, the positive semiaxis, , is invariant, and on the positive semiaxis, , the vector field associated with system (7) points to the interior of the positive quadrant. We then have that the closed positive quadrant is positively invariant.

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

Proposition 1. *If then any solution of system (7) with nonnegative initial conditions tends to .*

*Proof. *Let us call . Summing up both equations in system (7) and now choosing we have that, by integration, yieldssince is positively invariant.

Henceforth, we assume that . This assumption prevents the population from extinction. The linearization of system (7) at the equilibrium has the matrix with one positive and one negative eigenvalues. The unstable manifold of , associated with , is included in the axis, while the stable manifold, associated with , is tangent at to the eigenvector and lies completely outside the interior of the positive quadrant.

Assuming the only nonnegative solution tending to is itself. We prove next that all nonnegative solutions are forward bounded.

Proposition 2. *Let . If is any solution of system (7) with nonnegative initial conditions then it is bounded on .*

*Proof. *Calling and letting we have Function attains its maximum on at , so that Multiplying both sides of (12) by and rearranging terms yield which implies, integrating (13) on ,Rearranging terms in expression (14) leads toand, finally, we have that for every .

In addition to the trivial equilibrium , system (7) possesses a disease-free equilibrium , wherethat 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 represents its carrying capacity.

Proposition 3. *Let . The equilibrium point of system (7) verifies*(1)*if then is a saddle point, the stable manifold of which coincides with the positive semiaxis;*(2)*if then is locally asymptotically stable;*(3)*if then the basin of attraction of includes .*

*Proof. *To prove the two first items it suffices to calculate the matrix of the linearization of system (7) at :One of the eigenvalues, , is negative while stability depends on the other one, , being positive or negative: 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 equation is always negative for positive and . 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, .

*Up to now we have obtained the condition of nonextinction of the population, , and a sufficient condition, , for a disease-free scenario in the long term. This last condition says that if the competition coefficient , 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 . In this case the conditions of local stability of the equilibrium can be expressed in terms of the parameterThus, Proposition 3 can be restated as follows: for and , if or , then the equilibrium is a saddle point or locally asymptotically stable, respectively. The parameter 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, , a few infected individuals are able to spread the disease if the number of susceptible individuals is large enough, . On the other hand, the infection disappears if the susceptible population is under the threshold .*

*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 which can be interpreted through the terms depending on in the first equation of system (7). This equation can be written in the following form: where we note that depending on whether 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 , the more the infected individuals the larger the susceptible population growth rate, while yields a larger decrease of the susceptible population growth rate whenever there is a larger infected population.*

*Using parameters , , and , system (7) can be expressed as follows:*

*The equation of the -nullcline of system (21), for , is We are interested in the part included in the positive quadrant. This is for an increasing branch going from the point to the asymptote (see Figure 1, left panel) and for a decreasing branch going from the asymptote to the point (see Figure 1, right panel). In the case the -nullcline in the positive quadrant reduces to the line (see Figure 1, center panel).*