Plan
Comptes Rendus

Threshold conditions for infection persistence in complex host-vectors interactions
[Conditions seuil de persistence d’infection dans un système complexe d’interactions hôte–vecteur.]
Comptes Rendus. Biologies, Volume 325 (2002) no. 11, pp. 1073-1084.

Résumés

As classically defined by Macdonald in the early 1950s, for the case of diseases with one vector and one host, the Basic Reproduction Number, R0, is defined as the number of secondary infections caused by a single infective of the same type (vector or host) during its infectiousness period in an entirely susceptible population. In the case of a disease which has one vector and one host, it is easy to show that R0 coincides with the threshold for the establishment of an endemic state: if R0>1<1, the disease can invade (cannot invade) the host population. In this paper we examine various epidemic situations in which there are more than one vector and/or host. We show that in those more complex systems it is not possible to deduce a single R0 but rather a threshold for infection persistence which is a composite of several quantities closely related to the classical expression of R0. Another definition of R0 given by Diekmann, Heesterbeek and Metz, and denoted in this paper RNGO0 is discussed and applied as an alternative to calculate the thresholds for infection establishment.

Dans le cas d’une infection à un seul vecteur et à un seul hôte, le taux de reproduction de base R0 a été défini par Macdonald dans les années 1950 par le nombre d’infections secondaires causées par un seul individu infecté (hôte ou vecteur) pendant toute la durée de sa période infectueuse, au contact d’une population entièrement saine. Dans le cas d’une infection à un seul hôte et un seul vecteur, il est facile de montrer que R0 correspond à un seuil d’installation d’un état endémique: Si R0>1<1, l’épidémie envahit (disparaît de) la population hôte. Dans cet article, nous étudions diverses situations d’épidémies à plusieurs hôtes et/ou vecteurs. Nous montrons que, dans ces systèmes plus complexes, il n’est plus possible de définir un seul R0, mais plutôt un seuil de persistence de l’infection, qui est un agrégat de plusieurs quantités fortement connectées à l’expression classique de R0. Nous discutons une autre définition du R0, initialement proposée par Diekmann, Heesterbeck et Metz, notée RNGO0, qui est proposée comme une alternative pour calculer le seuil d’installation d’une infection.

Métadonnées
Reçu le :
Accepté le :
Publié le :
DOI : 10.1016/S1631-0691(02)01534-2
Keywords: basic reproductive number, vector-borne infections, yellow fever
Mot clés : taux de reproduction de base, infection transmise par un vecteur, fièvre jaune

Luiz Fernandes Lopez 1 ; Francisco Antonio Bezerra Coutinho 1 ; Marcelo Nascimento Burattini 1 ; Eduardo Massad 1

1 School of Medicine, The University of São Paulo and LIM01/HCFMUSP, Av. Dr. Arnaldo 455, São Paulo 01246-903, SP, Brazil
@article{CRBIOL_2002__325_11_1073_0,
     author = {Luiz Fernandes Lopez and Francisco Antonio Bezerra Coutinho and Marcelo Nascimento Burattini and Eduardo Massad},
     title = {Threshold conditions for infection persistence in complex host-vectors interactions},
     journal = {Comptes Rendus. Biologies},
     pages = {1073--1084},
     publisher = {Elsevier},
     volume = {325},
     number = {11},
     year = {2002},
     doi = {10.1016/S1631-0691(02)01534-2},
     language = {en},
}
TY  - JOUR
AU  - Luiz Fernandes Lopez
AU  - Francisco Antonio Bezerra Coutinho
AU  - Marcelo Nascimento Burattini
AU  - Eduardo Massad
TI  - Threshold conditions for infection persistence in complex host-vectors interactions
JO  - Comptes Rendus. Biologies
PY  - 2002
SP  - 1073
EP  - 1084
VL  - 325
IS  - 11
PB  - Elsevier
DO  - 10.1016/S1631-0691(02)01534-2
LA  - en
ID  - CRBIOL_2002__325_11_1073_0
ER  - 
%0 Journal Article
%A Luiz Fernandes Lopez
%A Francisco Antonio Bezerra Coutinho
%A Marcelo Nascimento Burattini
%A Eduardo Massad
%T Threshold conditions for infection persistence in complex host-vectors interactions
%J Comptes Rendus. Biologies
%D 2002
%P 1073-1084
%V 325
%N 11
%I Elsevier
%R 10.1016/S1631-0691(02)01534-2
%G en
%F CRBIOL_2002__325_11_1073_0
Luiz Fernandes Lopez; Francisco Antonio Bezerra Coutinho; Marcelo Nascimento Burattini; Eduardo Massad. Threshold conditions for infection persistence in complex host-vectors interactions. Comptes Rendus. Biologies, Volume 325 (2002) no. 11, pp. 1073-1084. doi : 10.1016/S1631-0691(02)01534-2. https://comptes-rendus.academie-sciences.fr/biologies/articles/10.1016/S1631-0691(02)01534-2/

Version originale du texte intégral

1 Introduction

Often unnoticed by practicing physicians in the temperate zone, arthropod-borne diseases account for a huge proportion of the spectrum of human maladies worldwide, and the problem appears to be growing 〚1〛. Despite the enormous effort of the medical and scientific community, controlling disease agents transmitted by arthropod vectors has proven to be difficult. The list of emerging and re-emerging infections is enormous but it is worth citing just a few: dengue, malaria, yellow fever, various mosquito-borne encephalitis, leishmaniasis, and Lyme disease.

Most of the techniques used for the control and eradication of vector-borne diseases were developed in the early 20th century. Rules for source reduction, insecticides, biological control, vaccination, chemotherapy and personal protection were all laid down nearly a century ago 〚2〛. Many of these techniques are still effective, others succeeded initially but failed later for a variety of reasons. Investigators must now incorporate new approaches that will allow them to move to the next level of control to alleviate the effects of vector-borne diseases on human and animal health 〚2〛.

The central parameter related to the intensity of transmission of infection is the so called basic reproduction number (R0), defined by Macdonald 〚3〛 as the number of secondary infections produced by a single infective in an entirely susceptible population (see next section). Originally applied in the context of malaria, R0 is a function of the vector population density as related to the host population, m, the average daily biting rate of the vector, a, the host susceptibility, b, the vector mortality rate, μ, the parasite extrinsic incubation period in days, n, and the parasitemia recovery rate, r, according to the (now) historical equation:

R0=ma2b exp -μnrμ
(actually, Macdonald denoted R0 as z0 in his original paper). From the definition of the basic reproduction number it can be demonstrated that if R0 is not greater than one, that is, when an index case (the first infective individual) is not able to generate at least one new infection, the disease dies out. Hence, in the original Macdonald analysis, R0 coincides with the threshold for the infection persistence. For an interesting historical account of R0 see 〚4〛.

This paper is organized as follows. In section 2 we revisit the Macdonald analysis of malaria transmission and deduce the general expression of R0 for any one vector – one host epidemic system. We show, using a dynamical system approach, that R0 coincides with the threshold for the infection persistence.

In section 3 we examine the definition of the basic reproduction number given by Dieckmann, Heesterbeek and Metz 〚5〛 as compared with the classical Mcdonald analysis.

In section 4 we examine more complex vector-hosts systems, exemplifying with a 2 vectors–2 hosts system, and show that there are four basic reproduction numbers of the Macdonald type, one for each kind of index case. Obviously, those R0's do not coincide with the threshold for epidemic persistence. However, we show that the threshold can be expressed in terms of those basic reproduction numbers. There is, however, only one basic reproduction number of the type defined by Dieckmann, Heesterbeek and Metz 〚5〛. As we shall see the same threshold can be deduced from it.

In section 5 we discuss an even more complex epidemic system, namely the yellow fever case, which comprises three different kinds of hosts and two vectors. In this case there are five R0's, one for each kind of index case. Also in this particular case, those R0's do not coincide with the threshold for epidemic persistence. Again, the threshold is given by a combination of them.

Finally, in the discussion section we summarize our findings.

2 The classical Macdonald analysis

In his 1952 seminal paper, Macdonald 〚3〛 addressed the problem of a system involving one vector (Anopheles mosquitoes) and one host (man). As mentioned above, his definition of R0 is the number of secondary infections in the first generation, that is, produced by a single infectee along his entire infectiousness period. We shall deduce an explicit expression for R0 from an intuitive perspective to show that it coincides with the threshold for the establishment of the disease. We do this because, as shown in the next section for more complex systems, this approach does not work in such a simple way.

Let us begin by assuming that the index case is a human host. The question to be answered is how many human secondary infections this index case produces in his/her entire infectiousness period.

Let Nm be the number of female mosquitoes. Let a be the average daily biting rate female anophelines inflict in the human population. The number of bites in the human population per unit of time is, therefore, Nma. Let Nh be the number of humans and r be the rate of recovery from parasitemia in the human cases. Therefore, the index case produces NmaNhRchm infected mosquitoes, where ch→m is the probability that a mosquito gets the infection after biting an infective human. These NmaNhrchm infected mosquitoes, in turn, produce aNmaNhrchm1μbmhe-μδ new human cases in the first generation, where 1μ is the average life expectancy of mosquitoes, bm→h is the probability that a human gets the infection after being bitten by an infective mosquito and e-μδ is the fraction of the infected mosquito population that survives through the extrinsic incubation period δ of the parasite. Note that, once infective a mosquito is assumed to remain so for life. Therefore, the expression for R0 is 〚3〛:

R 0 = a N m a N h r c h m 1 μ b m h e - μ δ

Similarly, if we begin with an infective mosquito as an index case, and compute the number of infected mosquitoes this index case produces in the first generation we get the same expression.

Let us now see how this deduction can be performed by a dynamical system approach.

Let Yh be the number of infected humans, and Yv the number of infected vectors. We can write

dYhdt=YvaNhbvhSh-rYhdYvdt=Svt-δaNhchme-μδYht-δ-μYv
where Sh and Sv are the number of susceptible humans and vectors, respectively.

To deduce the threshold for the disease to establish in the human population we analyse the stability of the trivial solution Sh=Nh, Sv=Nv, Yv=Yh=0, that is, the solution representing the absence of the infection. Linearising the system (3) around the trivial solution, we get

dyhdt=yvabvh-ryhdyvdt=NvaNhchme-μδyht-δ-μyv
where yv and yh are small deviations from zero. From the system (4), assuming solutions of the type yh=Aeλ𝑡 and yv=Beλ𝑡, we get the following characteristic equation for λ:

-λ+rabvhNvaNhchme-μδe-λδ-λ+μ
or

λ 2 + μ + r λ + μ r - N v a N h c h m e - μ δ e - λ δ a b v h = 0

It follows (see Appendix) that the roots of equation (5) or (6) have negative real parts if

μ r - N v a N h a b v h c h m e - μ δ > 0

The above result is the same as that obtained by the intuitive McDonald’s approach.

This still holds true for slightly more complex systems, such as those with one vector and two hosts populations or two vectors with one host populations. In these cases, the expression for R0 is partitioned in a sum with the individual terms of each component of the transmission chain 〚6〛.

3 The next generation operator

In a classical paper, Dieckman et al. 〚5〛 propose a new definition of the basic reproduction number for infections and we now study how it compares with the classical Macdonald definition described above.

Those authors define R0 as being the greatest eigenvalue of an operator which they call ‘the next generation operator’ (NGO). The case of vector-transmitted infections was analysed in a recent book by Dieckman and Heeterbeek 〚7〛.

In this section we give the next generation operator for the case of one-vector/one-host, exemplified by malaria. In this case, the next generation operator reduces to a two-by-two matrix

NGO = A v v A v h A h v A h h

The elements have the following interpretation. The element Av→h, for instance, means the number of infected humans generated by a single infected vector during its infectious period. Therefore, we have:

A v v = 0 A v h = a 1 μ b m h e - μ δ A h h = 0 A h v = N m a N h r c h m

Essentially, the NGO is the mean infectious output over all possible progressions of the infection within the host individual and therefore contain any process that influences output of infectious material to others. In particular, the NGO was derived without appealing to the dynamical system (3).

In this case the greatest eigenvalue of the NGO matrix, that is, RNGO0, is

R0 NGO =aNmaNhrchm1μbmhe-μδ
which is the square root of the Macdonald R0. It follows from the general theory of the Next Generation Operator 〚5〛 that if R0 NGO <1R0 NGO >1 the disease cannot (can) invade the host population.

4 Difficulties with more complex systems

The simple one vector/one host system characteristic of a great number of vector-borne diseases such as malaria, however, is not unique in the biology of the diseases transmitted by arthropods 〚2〛.

In some cases, more than one vector, or more than one host, and in some cases, more than one parasite may be involved. Examples of the complexity of the vector-borne infections abound (see, for instance any zoonoses textbooks, for instance Palmer, Soulsby and Simpson 〚8〛).

In some situations, several hosts may act as reservoirs of infection, being, therefore, intermediate steps in the chain of transmission. In the dengue/yellow fever complex, we have two (Haemagogus and Sabethes mosquitoes) or three (Haemagogus, Sabethes and Aedes mosquitoes) vectors for yellow fever, one common vector for dengue and yellow fever (Aedes mosquitoes), and at least two hosts (one human and one primate reservoir). As we demonstrate in the following sections, in such complex situations, a single threshold in the Macdonald’s R0 style simply does not exist.

4.1 The intuitive approach

Let us consider now a more complex system, namely one with two vectors and two host populations. This hypothetical system comprises an interaction between a human (Nh) and an animal (Nm) population, mediated by two species of vectors, referred to as (NA) and (NH). We shall try to deduce an explicit expression for R0 from the intuitive perspective used by Macdonald 〚3〛.

Let us begin by assuming that there is an index case in the human population. The question to be answered is how many human and animal secondary infections this index case produces in his/her entire infectiousness period.

Let aA be the average daily biting rate female mosquitoes of type A inflict in both host populations. The number of bites in the human population per unit of time is, therefore, NAaANhNh+Nm. Let rh be the rate of recovery from parasitemia in the human cases. Therefore, the index case produces NAaANhNh+Nm1Nh1rhchA infected mosquitoes, where 1Nh is the fraction of the bites inflicted in the index case, and ch→A is the probability that a mosquito gets the infection after biting an infective human. Those NAaANhNh+Nm1Nh1rhchA infected mosquitoes, in turn, produce

NAaANhNh+Nm1Nh1rhchAaA1μAbAhe-μAδA
new human cases in the first generation, where 1μA is the average life expectancy of mosquitoes, bA→h is the probability that a human gets the infection after being bitten by an infective mosquito and e-μAδA is the fraction of the infected mosquito population that survives through the extrinsic incubation period δA of the parasite. Note that, once infective, a mosquito is assumed to remain so for life.

By the same token, the same index case also produces

NHaHNhNh+Nm1Nh1rhchHaA1μHbHhe-μHδH
new human cases in the first generation, through the vector population H.

On the other hand, the same human index case also produces

NAaANhNh+Nm1Nh1rhchAaA1μAbAme-μAδA
secondary animal cases due to vector population of type A, and

NHaHNhNh+Nm1Nh1rhchHaA1μHbHme-μHδH
secondary animal cases due to vector population of type H.

Moreover, an animal index case also produces a similar amount of secondary animal and human infections due to each vector population.

Therefore, the human index case generates:

R hAh +R hHh
where

R hAh = N A a A N h N h + N m 1 N h 1 r h c h A a A 1 μ A b A h e - μ A δ A R hHh = N H a H N h N h + N m 1 N h 1 r h c h H a A 1 μ H b H h e - μ H δ H

human cases through the mosquitoes vectors A and H, respectively.

Moreover, the human index case also generates:

R hAm +R hHm
where

R hAm =NAaANhNh+Nm1Nh1rhchAaA1μAbAme-μAδAR hHm =NHaHNhNh+Nm1Nh1rhchHaA1μHbHme-μHδH
animal cases through the mosquitoes vectors A and H, respectively.

In the case when the index case is an animal, we find that it generates:

R mAh +R mHh
where

R mAh =NAaANmNh+Nm1Nm1rmcmAaA1μAbAhe-μAδAR mHh =NHaHNmNh+Nm1Nm1rmcmHaA1μHbHhe-μHδH
human cases through the mosquitoes vectors A and H, respectively.

It also generates:

R mAm +R mHm
where

R mAm =NAaANmNh+Nm1Nm1rmcmAaA1μAbAme-μAδAR mHm =NHaHNmNh+Nm1Nm1rmcmHaA1μHbHme-μHδH
animal cases through the mosquitoes vectors A and H, respectively.

Since, strictly speaking, the expression for the basic reproduction number should be calculated starting from an index case and counting the first generation number of cases of the same kind of the index case we have two basic reproduction numbers, depending on the kind of index case.

If the index case is a human being we tentatively set

R 0 h = R hAh + R hHh

On the other hand, if the index case is an animal, we tentatively set

R 0 m = R mAm + R mHm

However, we still have to investigate the case where the index case is a mosquito.

Assume that the index case is a mosquito of type A. It produces

aANhNh+Nm1μAbAh
infected humans and

aANmNh+Nm1μAbAm
animal cases. Those, in turn, produce

R AhA +R AmA
where

R AhA =NAaANhNh+Nm1Nm+Nh1rhchAaA1μAbAhe-μAδAR AmA =NAaANhNh+Nm1Nm+Nh1rhcmAaA1μAbAme-μAδA
infected mosquitoes of type A.

They also produce

R AhH +R AmH
where

R AhH =NHaHNhNh+Nm1Nm+Nh1rhchHaA1μAbAhe-μ𝐴δ𝐴R AmH =NHaHNhNh+Nm1Nm+Nh1rhcmHaA1μAbAme-μ𝐴δ𝐴
infected mosquitoes of type H.

So, if we stick to the convention that R0 is the number of secondary cases of the same type of the index case, we have two more expressions:

R0A=R AhA +R AmA
and

R 0 H = R HhH + R HmH

We, therefore, have a problem. We have four tentative R0 and we do not know if they are somehow linked with the threshold condition. Due to this fact we turn now to the dynamical system approach for the problem. First, however, we summarize and simplify our notation by defining the following transmission coefficients

ϐ Ah = a A N h N h + N m b A h ϐ Hh = a H N h N h + N m b H h ϐ Am = a A N m N h + N m b A m ϐ Hm = a H N m N h + N m b H m ϐ hA = a A N A N h + N m c h A e - μ A δ A ϐ mA = a A N A N h + N m c m A e - μ A δ A ϐ hH = a H N A N h + N m c h H e - μ H δ H ϐ mH = a H N A N h + N m c m H e - μ H δ H

Also, for future use, let us summarize the following parameters as:

R hAh = R AhA = ϐ hA ϐ Ah r h μ A R hHh = R HhH = ϐ hH ϐ Hh r h μ H R mAm = R AmA = ϐ mA ϐ Am r m μ A R mHm = R HmH = ϐ mH ϐ Hm r m μ H R hHm = ϐ hH ϐ Hm r h μ H R mAh = ϐ mA ϐ Ah r m μ A R mHh = ϐ mH ϐ Hh r m μ H R hAm = ϐ hA ϐ Am r h μ A R AhH = ϐ Ah ϐ hH r m μ H R HmA = ϐ Hm ϐ mA r m μ H R AmH = ϐ Am ϐ mH r m μ A R HhA = ϐ Hh ϐ hA r h μ H

Note also that we have

R hHm R mAh = R AhH R HmA R mHh R hAm = R AmH R HhA

Those parameters have a clear biological significance. For instance, RhAh is the basic reproduction number of a human index case who generates further human infections through the infective mosquito A. Note that in RhAh the mosquito population H and the population of animals m exist but are supposed not to contribute to RhAh. However, the animal population appear in the parameter RmHm through the factors NHNH+Nm and NANh+Nm.

4.2 The dynamical system approach

The dynamical system associated with such a situation is:

{ d Y h d t = Y A a A S h N h + N m b A h + Y H a H S h N h N m b H h - r h Y h d Y m d t = Y A a A S m N h + N m b A m + Y H a H S m N h N m b H m - r m Y m d Y A d t = N A a A Y h t - δ A N h + N m c h A e - μ A δ A + N A a A Y m t - δ A N h + N m c m A e - μ A δ A - μ A Y A dY H d t = N H a H Y h t - δ H N h + N m c h H e - μ H δ H + N H a H Y m t - δ H N h + N m c m H e - μ H δ H - μ H Y H

Linearising the system (13) around the solution without disease, namely Sh=Nh, Sm=Nm, SA=NA, SH=NH, Yh=Ym=YA=YH=0, we get the following system:

{dyhdt=yAϐ Ah +yHϐ Hh -rhyhdymdt=yAϐ Am +yHϐ Hm -rmymdyAdt=yht-δAϐ hA +ymt-δAϐ mA -μAyAdyHdt=yht-δHϐ hH +ymt-δHϐ mH -μHyH
where yii=h,m,A,H are small deviations from zero. From the system (14) we get the following characteristic equation:

- λ + r h 0 ϐ Ah ϐ Hh 0 - λ + r m ϐ Am ϐ Hm ϐ hA e - λ δ A ϐ mA e - λ δ A - λ + μ A 0 ϐ hH e - λ δ H ϐ mH e - λ δ H 0 - λ + μ H = 0

The trivial solution will then be stable if all the roots of the characteristic equation (15) have negative real parts. However, it is almost hopeless to study such a complex characteristic equation as (15), and a sensible procedure should be to neglect all the terms e-λδii=A,H, in which case equation (15) would became a fourth order algebraic equation which can be analysed using the Routh–Hurwitz criteria 〚9〛. This can be justified very laboriously using perturbation theory as is done for the simple case of one-host/one-parasite in the appendix. However, it is easier to apply the Next Generation Operator approach to calculate the threshold conditions for the establishment of the epidemic. This is the subject of the next section.

4.3 The next generation operator approach

For the reasons explained in the Appendix, we can neglect the time delays from the very beginning. In this case, the linearised system (14) becomes:

{dyhdt=yAϐ Ah +yHϐ Hh -rhyhdymdt=yAϐ Am +yHϐ Hm -rmymdyAdt=yhϐ hA +ymϐ mA -μAyAdyHdt=yhϐ hH +ymϐ mH -μHyH
which can be written in the form

d𝐲¯dt=𝐓-𝐃𝐲¯
where T is a positive matrix and D is a diagonal positive matrix (see Dieckman and Heesterbeck 〚7〛, p. 105). The next generation operator in this case is given by a four by four matrix:

NGO = TD - 1 = 0 0 ϐ hA μ A ϐ hH μ H 0 0 ϐ mA μ A ϐ mH μ H ϐ Ah r h ϐ Am r m 0 0 ϐ Hh r h ϐ Hh r m 0 0

Although in this case we derived the NGO from the linearised form of the system (14), given by equation (16), this is not necessary. The NGO can be derived independently and, in particular, it does not depend on the assumption made in going from (14) to (16) of neglecting the time delays.

According to the Dieckmann et al. 〚5〛 (see also 〚10〛) theorem, the R0 given by the Next Generation Operator, RNGO0, is the largest eigenvalue of matrix (18). We then have

R0 NGO =12B+B2+4A
where

B=R AhA +R HhH +R AmA +R HmH =R hAh +R hHh +R mAm +R mHm
and

A = R hHm R mAh + R mHh R hAm - R mHm R hAh - R mAm R hHh

If RNGO0<1, then the solution without disease of system (13) is stable. Therefore, the threshold condition is given, then, by

R0 NGO =12B+B2+4A1
which implies

A+B1
or

T = 1 - R hAh - R hHh - R mAm - R mHm - R hHm R mAh - R mHh R hAm + R mHm R hAh + R mAm R hHh 0

The same conclusion can be obtained, more laboriously, by applying the Routh–Hurwitz criteria to equation (15), and taking the exponential terms e-λδi (i = A, H) as 1.

Some particular cases can be analysed. Suppose that all RiJi are zero except one. Then, T becomes negative if this non-zero RiJi is greater than 1, and the disease can establish itself in the host population. Another more interesting case occurs when, for instance, RhAh=RhHh=0. One might think, at first sight, that in this case a human index case cannot trigger a human epidemic. However, it can be seen that the threshold, which determines whether T is positive or negative is given by

T=1-R mAm -R mHm -R hHm R mAh -R mHh R hAm
which has a simple interpretation. In this case the fourth term in the brackets, RhHmRmAh, describes the basic reproduction number for a human index who transmit the infection to the animal population through the mosquito type H. The infected animal then passes the disease, through mosquito type A back to the human hosts. Therefore, even when RhAh=RhHh=0, a human index case can trigger a human epidemic. For this, it suffices that RhHmRmAh>1. The other terms in the expression (20) refer to transmission among the animal population by all the possible ways and all can be less than 1 without preventing the disease establishing itself.

By the same token, if we make RmAmRmHm=0 (which apparently implies that an animal index case cannot trigger an animal epidemic), then

T=1-R hAh -R hHh -R mAh R hHm -R hAm R mHh
which can be negative when RmAhRhHm>1, even when the remaining components RiJi are all less than 1.

Another attempt to interpret our condition is as follows. Suppose that when the index case is h, and RmAm=RmHm=0. Then, we have that T reduces to (21). Suppose now that when the index case is m, and RhAh=RhHh=0. In this case T reduces to

T = 1 - R mAm - R mHm - R mAh R hHm - R mHh R hAm

To investigate the cases where the index cases are mosquitoes, we rewrite T using the relations given by equations (11) and (12). We get:

T = 1 - R AhA - R HhH - R AmA - R HmH - R AhH R HmA - R AmH R HhA + R HmH R AhA + R AmA R HhH

If we now assume that the index cases is a mosquito of the type A, we put RHhH=RHmH=0 to obtain

T = 1 - R AhA - R AmA - R AhH R HmA - R AmH R HhA

Similarly, when the index case is a mosquito of type H we put RAhA=RAmA=0 to obtain

T = 1 - R HhH - R HmH - R HmA R AhH - R HhA R AmH

5 The yellow fever model

Urban yellow fever is transmitted from person to person by peridomestic Aedes aegypti mosquitoes 〚11〛. By contrast, jungle yellow fever is a zoonosis, transmitted from monkeys to humans by mosquitoes that breed in tree-holes of the genuses Haemagogus and Sabethes in the rain forest ecosystem of South America. The jungle form is only partly controlled by vaccination of rural residents and provides a source of infection to population centres infested with Ae aegypti. In the early 20th century, when it was discovered that the yellow fever virus was transmitted in its urban cycle by Aedes aegypti, measures of control were introduced leading to its disappearance. Progressive neglect of the disease, however, led to a new outbreak in Africa in 1927 〚12〛 during which the etiological agent was isolated; some years later a vaccine was discovered and yellow fever disappeared again. Unfortunately, reinfestation with the Aedes vector, which began in the 1970s, is now virtually complete, and vector control is substantially more difficult than before. The threat of urban yellow fever is greatest in towns near forests, but improved transport links increase the likelihood of spread by viremic people to non-endemic areas 〚13〛.

Classified as one of the viral hemorrhagic fevers, yellow fever is unique in its severity, in particular because of its hepatic impairment. Yellow fever is currently endemic and epidemic in tropical areas of the Americas and Africa 〚14–16〛.

From the point of view of its dynamics, yellow fever differs from other vector borne infection by involving three host populations and two vector populations and, therefore, as we shall see, the threshold condition for non-existence of the infection is quite complicated. We recently calculated the risk of urban yellow fever in a dengue infested area 〚17〛.

5.1 The dynamical system approach

Let Yct be the number of infected human individuals living in cities, Yft the number of infected ‘fishermen’ (those individuals who either live part of time in forest areas or who eventually frequent those areas for leisure or other purposes), Ymt the number of infected non-human primates, YAt the number of infected mosquitoes of the Aedes genus, and YHt the number of infected mosquitoes of the Haemagogus or Sabetes genuses. Let also Ni i=c,f,m,A,H be the total number of individuals in each population considered.

The linearised equations for Yi around the solution without the disease are:

{dYcdt=YAaANcNc+NfbAc-rcYcdYfdt=YAaANfNc+NfbAf+YHaHNfNc+NfbHf-rfYfdYmdt=YHaHNmNm+NfbHm-rmYmdYAdt=NAaAYct-δANc+NfgcAe-μAδA+NAaAYft-δANm+NfgfAe-μAδA-μAYAdYHdt=NHaHYft-δHNm+NfgfHe-μHδH+NHaHYmt-δHNm+NfgmHe-μHδH-μHYH
where aii=A,H is the average number of bites per unit time that the vectors inflict on the hosts; bij (i = A, H and j = cf) is the fraction of the potentially infective bites aiYi per unit time, which are actually infective; gji (i = A, H and j = cf) is the fraction of the bites inflicted on infective hosts that are really infective to the vectors. Finally, rj (j = c, f, m) is the rate with which the hosts leave the infective compartment (either by death or by recovery from the infection), δi (i = A, H) is the incubation period of the virus in the mosquito population (that is, the extrinsic incubation period, the time elapsed between infection and infectiousness in the vectors) and e-μiδi is the natural mortality rate of the vectors (therefore μi (i = A, H) is the fraction of mosquitoes which survives throughout the extrinsic incubation period).

A few words about the system (26). Take for instance the second equation: in it, the term YAaA represents the total number of potentially infective bites Aedes mosquitoes inflict in the fraction NfNc+Nf of the ‘fishermen’ hosts, and bA→f is the fraction of those bites which are actually infective for the host. Therefore, the complete term YAaANfNc+NfbAf is the number of ‘fishermen’ who get the infection per unit of time from the Aedes population. Note that we have linearised the original equations. In its complete form, the fraction NfNc+Nf would be SfNc+Nf, where Sf is the number of susceptible ‘fishermen’. By the same token, the second term YHaHNfNc+NfbHf is the number of ‘fishermen’ who get the infection per unit of time from the Haemagogus population. And finally, the last term rfYf represents the number of ‘fishermen’ removed from the infectious state by death or recovery.

Similarly, the equations for the mosquitoes, for instance, the fourth equation could be explained in words as follows. The first term NAaA represents the total number of bites Aedes mosquitoes inflict in the fraction Yct-δANc+Nf of infected citizens t-δA days ago, gc→A is the fraction of those bites that are actually infective for the mosquitoes, and e-μAδA is the fraction of mosquitoes that survives throughout the extrinsic incubation period, therefore becoming infective to the hosts. Therefore, the complete term NAaAYct-δANc+NfgcAe-μAδA represents the number of Aedes per unit time which gets the infection from the ‘citizens’ hosts. Note that the original non-linearised system should be SAt-δAaAYct-δANc+NfgcAe-μAδA, where SAt-δA is the number of susceptible Aedes mosquitoes (t-δA) days ago. The second term of the fourth equation NAaAYft-δANm+NfgfAe-μAδA, by the same reasoning, represents the number of Aedes per unit time which gets the infection from the ‘fishermen’ hosts. And finally, the last term μAYA represent Aedes removed from the infective condition by mortality.

Calculating the threshold condition for yellow fever

In order to simplify the notation let us define the following transmission coefficients

ϐ Ac = a A N c N c + N f b A c ϐ Af = a A N f N c + N f b A f ϐ Hf = a H N f N c + N f b H f ϐ Hm = a H N m N m + N f b H m ϐ cA = a A N A N c + N f g c A e - μ A δ A ϐ fA = a A N A N m + N f g f A e - μ A δ A ϐ fH = a H N H N m + N f g f H e - μ H δ H ϐ mH = a H N H N m + N f g m H e - μ H δ H

To calculate the threshold condition we examine the stability of the trivial solution of system (26). We get the following characteristic equation, written as a determinant:

- r c + λ 0 0 ϐ Ac 0 0 - r f + λ 0 ϐ Af ϐ Hf 0 0 - r m + λ 0 ϐ Hm ϐ cA e - λ δ A ϐ fA e - λ δ A 0 - μ A + λ 0 0 ϐ fH e - λ δ H ϐ mH e - λ δ H 0 - μ H + λ = 0

Again, the trivial solution will then be stable if all the roots of the characteristic equation (27) have negative real parts. As mentioned above, it is almost hopeless to study such a complex characteristic equation as (27), and a sensible procedure should be to neglect all the terms e-λδi (i = A, H), in which case equation (27) would became a fifth order algebraic equation which can be analysed using the Routh–Hurwitz criteria 〚9〛. This can be justified by using perturbation theory as done in the appendix. However, once again, it is easier to apply the Next Generation Operator approach to calculate the threshold conditions for the establishment of the epidemic.

5.2 The next generation operator approach

The next generation operator approach for the above model of yellow fever is given by

NGO = 0 0 0 ϐ Ac μ A 0 0 0 0 ϐ Af μ A ϐ Hf μ H 0 0 0 0 ϐ Hm μ H ϐ cA e - μ A δ A r c ϐ fA e - μ A δ A r f 0 0 0 0 ϐ fH e - μ H δ H r f ϐ m e - μ H δ H r m 0 0

According to the Dieckmann et al. theorem, the R0 given by the Next Generation Operator, RNGO0, is the largest eigenvalue of matrix (28). We then have

R0 NGO =12B+B2+4A
where

B=R fHf +R fAf +R mHm +R cAc
and

A = - R fAf R mHm - R cAc R mHm - R cAc R fHf

The threshold condition is given, then, by

R0 NGO =12B+B2+4A1
which implies

A+B1
or,

T = 1 - R cAc + R fAf + R fHf + R mHm - R cAc R fHf - R cAc R mHm - R mHm R fAf 0

Some particular cases can be of interest. Suppose that transmission to fishermen is zero, that is, RfAf and RfHf are equal to zero. Then:

T = 1 - R cAc + R mHm - R cAc R mHm 0

Therefore, the threshold occurs when either RmHm or RcAc is greater than one. However, if both RmHm and RcAc are greater than one, then T may become negative. Therefore we have to apply the Routh–Hurwitz 〚9〛 condition to each particular case. However, there will obviously be an epidemic because the population of humans living in cities becomes decoupled from the population of non-human primates.

Suppose now that we have transmission only to fishermen. In this case, the threshold becomes

T = 1 - R fAf + R fHf 0

It is obvious that the threshold occurs when T becomes negative and is given by

R fAf <1R fHf <1
but

R fAf +R fHf >1
which is intuitively very reasonable.

6 Summary and discussion

The basic reproduction number, as defined by Macdonald 〚3〛 is the average number of secondary cases produced by an index case during its infectiousness period. Note that the Macdonald definition implies secondary cases of the same kind as the index case. Furthermore, the basic reproduction number in a disease which involves only one host and one vector coincides with the threshold that breaks the stability of the trivial solution, as shown in this paper. Another definition of R0 is the largest eigenvalue of the Next Generation Operator. The dominant eigenvalue of the NGO gives the average multiplication factor with which successive generations grow. It is therefore precisely the average number of cases in the next generation of infecteds produced by a typical case in the present generation. Thus, a threshold is immediately apparent: the disease cannot invade the population if the dominant eigenvalue of the NGO is less than 1.

Suppose now an indirectly transmitted disease which involves more than one type of host and/or vector. Then, as we have shown in this paper, several basic reproduction numbers can be defined. Notwithstanding, there is only one threshold for the epidemic persistence that can be written as a function of the many R0's and some other quantities, which are not R0, but are similar. In fact, they are the number of secondary cases of host(s)/vector(s) of a different kind of the index case.

In section 2, we revisited the classical Macdonald model, calculated the basic reproduction number and showed, by using a dynamical system approach, that the threshold and the R0 are the same quantity.

In section 3, we discussed in detail the case of an hypothetical infection system with two hosts (denoted h and m) and two vectors (denoted A and H). We showed that this system has four R0, namely

R hAh = R AhA R hHh = R HhH R mAm = R AmA R mHm = R HmH

We also calculated the threshold for the epidemic persistence, which is given by a quantity T (see equation (19)), namely:

T = 1 - R hAh - R hHh - R mAm - R mHm - R hHm R mAh - R mHh R hAm + R mHm R hAh + R mAm R hHh

If T is greater than zero there is no epidemic persistence. Note that T involves terms that are R0 as defined above (like RhAh or RhHh), but it involves as well terms that are not R0 but refers to the number of human or mosquitoes secondary infections produced by a single infective mosquitoe or human, respectively (like RhHm or RmAh).

Acknowledgements

The authors would like to thank Professor Hans Heesterbeek for his helpful comments on the manuscript. Eventual remaining errors are our own responsibility.

Appendix

We write equation (6) as:

λ2+a1λ+a21-a22e-λδ=0
where

a 1 = μ + ra 21 = μ ra 22 = N v a N h c h m e - μ δ a b v h

Assume a root

z = x + i y

We can separate equation (34) in its real and imaginary parts. We get:

x2-y2+a1x+a21-a22e𝑥δ cos yδ=0
and

2 x y + a 1 y - a 22 e 𝑥 δ sin y δ = 0

Let us solve (35) and (36) perturbatively around the solution obtained with δ=0. We write

x=x0+x1δ+x2δ2+...
and

y = y 0 + y 1 δ + y 2 δ 2 + . . .

Replacing (37) and (38) in (35) and (36), we get, keeping only terms of zero order in δ:

x02+a1x0+a21-a22=02x0y0+a1y0=0
whose solutions are

y0=0x0=-a1±a12-4a21-a222
and

x 0 = - a 1 2 y 0 = ± 4 a 21 - a 22 - a 1 2 2

However, the second solution (41) represents a root with negative real part, and therefore is uninteresting. We can see that in the zero order in δ a root crosses the imaginary axis through the real axis when a21-a22→0.

The terms in first order in δ give

2x0x1-2y0y1+a1x1-a22x0=0
and

2 x 0 y 1 + 2 x 1 y 0 + a 1 y 1 - a 22 y 0 = 0

Replacing y0=0 in equations (42) and (43), we get:

2x0x1+a1x1-a22x0=02x0y1+a1y1=0
whose solution is

y 1 = 0 x 1 = a 22 x 0 2 x 0 + a 1

Then we see that when a21-a220, x10, together with x0.

The terms in second order in δ, after substituting y0=0 and y1=0, give

x12+2x0x2+a1x2-a22x1-12a22x02=02x0y2+a1y2=0
whose solution in

y 2 = 0 x 2 = a 22 x 1 - x 1 2 + 1 2 a 22 x 0 2 2 x 0 + a 1

Again we see that when a21-a220, x20 together with x0 and x1.

So, up to second order in δ, a root crosses the imaginary axis through the real axis when a21-a220. It is not difficult to see that all terms in the expansion vanish when a21-a220, showing that a21-a22=0 is the threshold for the infection to persist.


Bibliographie

[〚1〛] J. Goddard Infectious Diseases and Arthropods, Humana Press Inc, Totowa, New Jersey, 2000

[〚2〛] B.J. Beaty; W.C. Marquardt The Biology of Diseases Vectors, University Press of Colorado, Niwot, Colorado, USA, 1996

[〚3〛] G. Macdonald The analysis of equilibrium in malaria, Trop. Dis. Bull., Volume 49 (1952), pp. 813-828

[〚4〛] K. Dietz The estimation of the basic reproduction number for infectious diseases, Stat. Meth. M, Volume 2 (1993), pp. 23-41

[〚5〛] O. Diekmann; J.P. Heesterbeek; J.A.J. Metz On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations, J. Math. Biol. (1990), pp. 365-382

[〚6〛] M.N. Burattini; F.A.B. Coutinho; L.F. Lopez; E. Massad Modelling the dynamics of leishmaniasis considering human, animal host and vector populations, J. Biol. Syst., Volume 6 (1998) no. 4, pp. 337-356

[〚7〛] O. Diekmann; J.P. Heesterbeek Mathematical epidemiology of infectious diseases, John Wiley & Son Ltd, Chichester, 2000

[〚8〛] S.R. Palmer; L. Soulsby; D.I.H. Simpson Zoonoses. Biology, Clinical Practices and Public Health Control, Oxford University Press, Oxford, 1998

[〚9〛] J.D. Murray Mathematical Biology, Springer Verlag, Berlin, 1993

[〚10〛] A. Nold Heterogeneity in disease-transmission modeling, Math. Biosci., Volume 52 (1980), pp. 227-240

[〚11〛] T.P. Monath Yellow Fever (S.A. Plotkin; W.A. Orestein, eds.), Vaccines, W.B. Sounders Company, Philadelphia, USA, 1999, pp. 815-880

[〚12〛] J.P. Digoutte Present status of an arbovirus infection: yellow fever, its natural history of hemorrhagic fever, Rift Valley fever, Bull. Soc. Pathol. Exot., Volume 92 (1999) no. 5, pp. 343-348

[〚13〛] T.P. Monath Facing up to re-emergence of urban yellow fever, The Lancet, Volume 353 (1999) no. 9164, p. 1541

[〚14〛] T.P. Monath Yellow Fever (K.S. Warren; A.A.F. Mahmoud, eds.), Tropical and Geographical Medicine, McGraw-Hill Information Services Company, New York, 1990, pp. 661-674

[〚15〛] E.J. Sanders; A.A. Marfin; P.M. Tukei; G. Kuria; G. Ademba; N.N. Agata; J.O. Ouma; C.B. Cropp; N. Karabatsos; P. Reiter; P.S. Moore; D.J. Gubler First recorded outbreak of yellow fever in Kenya, 1992–1993. I. Epidemiologic investigations, Am. J. Trop. Med. Hyg., Volume 59 (1998) no. 4, pp. 644-649

[〚16〛] P. Reiter; R. Cordellier; J.O. Ouma; C.B. Cropp; H.M. Savage; E.J. Sanders; A.A. Marfin; P.M. Tukei; N.N. Agata; L.G. Gitau; B.A. Rapuoda; D.J. Gubler First recorded outbreak of yellow fever in Kenya, 1992–1993. II. Entomologic investigations, Am. J. Trop. Med. Hyg., Volume 59 (1998) no. 4, pp. 650-656

[〚17〛] E. Massad; F.A.B. Coutinho; M.N. Burattini; L.F. Lopez The risk of yellow fever in a dengue-infested area, Trans. R. Soc. Trop. Med. Hyg., Volume 95 (2001) no. 3, pp. 370-374


Commentaires - Politique