# Basic Properties and Qualitative Dynamics of a Vector-Borne Disease Model with Vector Stages and Vertical Transmission.

1. IntroductionVector-borne diseases have been the scourge of man and animals since the beginning of time [1]. Today, vector-borne diseases account for over 17% of all infectious diseases causing more than 1 million deaths annually, and their distribution is mainly determined by a complex dynamic of environmental and social factors [2]. In spite of all these inherent complexities, mathematical models have been used to translate assumptions concerning biological, environmental, and social aspects into mathematical structures, linking biological processes of transmission and dynamics of infection at population level. Such dynamic models have impacted both our understanding of epidemic spread and public health planning (for more details see [3-5] and references therein).

In this study our particular interest is in investigating qualitative properties of epidemic models of mosquito-borne diseases in which the vector is of genera Aedes. These mosquito species are known to transmit many vector-borne diseases of vast epidemiological importance including Dengue fever and Rift Valley fever (RVF), just to mention a few. These diseases are endemic in Sub-Saharan Africa with pronounced health and economic impacts on domestic animals and humans. In fact, losses due to RVF can reach millions of dollars during periods of disease outbreaks [6, 7]. An interesting phenomenon underlying many vector-borne diseases is their ability to persist year round, fluctuating seasonally but only falling to zero during some months. Hence, an important question is, how does the virus maintain itself in nature? For RVF it has been hypothesized that RVF virus (RVFV) is maintained through transovarial transmission in Aedes mosquito eggs [8]. Aedes eggs need be dry for several days before they can mature. After maturing, they hatch during the next flooding event large enough to cover them with water [9-11]. The eggs have high desiccation resistance and can survive dry conditions in a dormant form for months to years [12-15]. Thus, the emergence of adult Aedes mosquitoes from infected eggs can reintroduce RVF in livestock at the beginning of the rainy season, before other mosquitoes species amplify it further [16]. For instance, in eastern and southern Africa there is more and more evidence of disease activities between outbreaks [7,17-20] highlighting the role of vertical transmission for initial disease spread and endemicity.

In epidemiology disease spread and persistence are measured through quantities known as epidemic thresholds. Their derivation and characterization are one of the most important results of mathematical epidemic models. The basic reproductive number, [R.sup.0], is the most critical epidemic threshold given its applicability and suitability for deciding whether an outbreak will occur or fade out, making it essential for guiding disease control efforts. However, the derivation of this epidemic threshold in vector-borne disease models in particular suffers from a lack of uniqueness and it fails to give the correct average number of expected secondary infections produced by one infected individual [21]. This failure is more likely to compound itself when vertical transmission mode is included in the transmission model, since the resulting [R.sup.0] comes as the sum of the vertical and horizontal transmission components, if the next-generation method is used. Previous mathematical models have made a significant attempt in including vertical infection in modelling vector-borne diseases [22-28], but none of them discuss how the failure of [R.sup.0] compounds itself in the presence of vertical transmission. Therefore, the present work aims to discuss some relevant basic properties of vector-borne disease models when vertical infection is taken into account and their implications for disease control efforts. Further, our goal is to derive new epidemic thresholds useful for guiding control efforts in the settings of vector-borne disease models with vector stages that include vertical transmission mode.

We formulate two models, one simple but realistic and the other more complex with vector stages. The first is an extension of the one proposed by Ross [29] and popularized by Macdonald [30] and Anderson [31]. The model is used to discusssystem propertiessuchasthe asymmetric relationship between the host-to-vector and vector-to-host reproductive numbers. In addition, we highlight how to derive epidemic thresholds useful for guiding disease control efforts and discuss their relationships with vertical transmission efficiency. It is shown that the model has two model equilibria, namely, the disease-free and the endemic, and Lyapunov function theory is used to establish their global qualitative dynamics. The second is an extension of the basic model, where the dynamics of both aquatic and adult mosquitoes are modelled explicitly. In this model we let the populations of aquatic and adults vary with time but be limited by their respective carrying capacity. The inclusion of the explicit vector submodel allows for derivation of critical thresholds such as the reproductive number for both the vector population and the disease system. Then, these thresholds are used to determine global qualitative dynamics of both the disease-free and endemic equilibria.

The paper is set out as follows. In Section 2 we formulate and discuss the two model systems. In Section 3 we provide the epidemic threshold theorems regarding both the vector population and disease equilibria. We also discuss important model properties and how to derive and identify model epidemic thresholds useful for guiding disease control efforts. Furthermore, numerical simulations are carried out to investigate the influence of the key parameters on the spread of the disease (taking RVF as disease example), to support analytical analyses and conclusions and illustrate possible behavioural scenarios of the model with vector stages. Finally, in Section 4 we present a short discussion of the results and their biological implications.

2. Materials and Methods

For human and animal diseases, horizontal transmission typically occurs through direct or indirect physical contact with infectious hosts, or through disease vectors such as mosquitoes, ticks, or other biting insects. Among mosquito vectors vertical transmission is often through eggs. Of particular interest are female mosquitoes of genera Aedes which transmit the virus to their eggs. These eggs have some adaptive behaviour which allows them to stay dormant in nature for relatively long periods. Although vertical transmission also occurs among vectors involved in the transmission of Dengue disease, RVF is the disease for which the model is a good approximation. In particular, parameter values related to RVF are used to illustrate the dynamics of the disease numerically.

2.1. Host-Vector Basic Model with Vertical Transmission. Let [N.sub.h] = [S.sub.h] + [I.sub.h] + [R.sub.h] and [N.sub.v] = [S.sub.v] + [I.sub.v] denote the total host and vector populations sizes, respectively. We assume that individuals at each compartment mix homogeneously and each mosquito bites each individual host at a constant rate [alpha]/[N.sub.h], where [alpha] is the biting rate per unit time. Let [p.sub.hV] be a probability of successful infection transmission from an infected mosquito to a susceptible host and [p.sub.Vh] be a probability of successful infection transmission from an infected host to a susceptible mosquito per bite. Thus, the forces of infection are as follows, [mathematical expression not reproducible]. Hosts are recruited into the population at per capita rate [[mu].sub.h] which is proportional to the total population and leave each compartment through death. Noninfected mosquitoes join the susceptible compartment at rate [[mu].sub.v]([N.sub.v] - [q.sub.v] [I.sub.v]) while vertically infected mosquitoes join the infected class at rate [[mu].sub.v][q.sub.v][I.sub.v]. Assuming constant population sizes, that is, births equal to deaths, [R.sub.h] can be obtained when both [S.sub.h] and [I.sub.h] are known. In the same way, [S.sub.v] can be obtained when [I.sub.v] is known. Thus, the expressions for both [R.sub.h] and [S.sub.v] can be omitted and the system can be written in terms of proportions:

[mathematical expression not reproducible] (1)

where [q.sub.v] [member of] (0,1) is the probability of vertical infection, [[gamma].sub.h] the rate at which infected hosts recover from infection, and m denotes the ratio of female mosquitoes to hosts.

2.2. Model with Vector Stages and Vertical Transmission. Here we extend the basic model to include vector stages. Partial results of the resulting model without vertical infection have been obtained in [32]. Our aim is to extend their model and analysis by investigating global dynamics of all model equilibria and examine the extent to which vertical infection alters the dynamics of the system. The mosquito population is divided into aquatic (eggs, larvae, and pupae) and terrestrial (adults) subpopulations with [K.sub.1] and [K.sub.2] being their carrying capacity, respectively. The parameter [K.sub.1] represents the larval maximal capacity limited by the availability of breeding sites while [K.sub.2] is the maximal capacity of adult mosquitoes limited by factors conditioning their survival such as high altitudes and high temperatures. Further, the aquatic subpopulation is divided into epidemiological classes, susceptible ([L.sub.s]), and infected ([L.sub.i]) while adults are divided also into susceptible ([S.sub.v]) and infected ([I.sub.v]). The per capita oviposition rate is r(1-L/[K.sub.1]), where r is the intrinsic oviposition rate and L = [L.sub.s] + [L.sub.i]. Aquatic mosquitoes emerge as adults at a per capita rate [theta] where the proportion 1 - [q.sub.v] emerge noninfected while the remainder are infected. Disease transmission dynamics between vector and host populations remain the same as in the basic model. As a result the following nondimensional system of ordinary differential equations represents the model that governs the temporal evolution of the disease:

[mathematical expression not reproducible] (2)

Let [mathematical expression not reproducible] be the initial conditions of system (2). It is easy to check that the feasible region for (2) is the positive orthant of R6 and that the closed set

[mathematical expression not reproducible] (3)

is positively invariant for system (2).

2.2.1. Positivity of Solutions

Lemma 1 (see [33]). Let us denote [mathematical expression not reproducible] and consider the function f : [R.sub.+] x [R.sup.n] [right arrow] [R.sup.n] continuous with respect to t, u and Lipschitz with respect to u. If f(t, u) [greater than or equal to] 0 for (t, u) [member of] [R.sub.+] x [R.sup.n], with [u.sub.i] = 0, then, for every [u.sub.0] [member of] [R.sup.n.sub.+], there exists T > 0 such that the solution to

du/dt = f (t,u), u ([t.sub.0]) = [u.sub.0] (4)

exists and is unique and positive with value in [R.sup.n.sub.+] and defined on some interval [0, T). If T < [infinity], then

[mathematical expression not reproducible]. (5)

Theorem 2. The solution set {[L.sub.s](t), [L.sub.i](t), [S.sub.v](t), [I.sub.v](t), [S.sub.h](t), [I.sub.h](t)} of the model (2) exists and is unique and positive for t > 0.

Proof. Let X = ([L.sub.s](t), [L.sub.i](t), [S.sub.v](t), [I.sub.v](t), [S.sub.h](t), [I.sub.h](t)) and denote the function [mathematical expression not reproducible] such that

[mathematical expression not reproducible] (6)

Since the function f is continuous and Lipschitz continuous with respect to X, according to Picard's theorem, there exists [T.sub.0] >0 such that the solution to (2) exists locally at least on an interval of this form [0, [T.sub.0]]. Further, considering the initial condition [X.sub.1] = X([T.sub.0]) at [t.sub.0] = [T.sub.0] and using Picard's theorem, it follows that there exists [T.sub.0] [less than or equal to] [T.sub.1] [member of] [R.sub.+] such that the solution to (2) exists and is unique on [[T.sub.0], [T.sub.1]]. Since f is continuous and differentiable, the solution of (2) with a given initial condition is unique. Therefore, the solutions of (2) obtained on [0, [T.sub.0]] and on [[T.sub.0], [T.sub.1]] form the unique solution of (2) on [mathematical expression not reproducible] with the initial condition [X.sub.0] at [t.sub.0] = 0. Repeating this process again and again, we end up with the maximal forward interval of existence for the solutions of (2), say [mathematical expression not reproducible]. Furthermore, for X e [R.sup.6.sub.+],

[f.sub.i] ([L.sub.s] (t), [L.sub.i], (t), [S.sub.v] (t), [I.sub.v] (t), [S.sub.h] (t),0) [greater than or equal to] 0. (7)

Therefore, the solutions of (2) on [??] are positive, according to Lemma 1.

Finally, according to Theorem 2, the solutions to (2) are bounded on [0, [??]). In other words, they do not blow up on any finite interval of [R.sub.+]. It follows that, according to Lemma 1, the solution of (2) exists for all time. Hence for any initial condition in [R.sup.6.sub.+], system (2) possesses a unique and positive solution in [R.sup.6.sub.+].

3. Results

3.1. Analysis of the Basic Host-Vector Model

3.1.1. Model Equilibria and Stability Analysis. The basic host-vector model with vertical transmission exhibits two equilibria, namely, the disease-free [E.sup.0.sub.1] and the endemic [E.sup.*.sub.1], respectively. At the disease-free equilibrium,

[mathematical expression not reproducible], (8)

both vector and host populations persist but with no disease. The prevalence of the disease is denoted by

[mathematical expression not reproducible]

where [mathematical expression not reproducible], and [mathematical expression not reproducible] (10)

The Jacobian matrix of system (1) at [E.sup.0.sub.1] is given by

[mathematical expression not reproducible], (11)

such that the characteristic polynomial of matrix (11) is then given as

[mathematical expression not reproducible], (12)

with [a.sub.2] = V + c + [[mu].sub.h], [a.sub.1] = Vc(1 - [R.sup.e.sub.0]) + (V + c)[[mu].sub.h], and [a.sub.0] = [[mu].sub.h] Vc(1 - [R.sup.e.sub.0]). The coefficient [a.sub.2] > 0 and both [a.sub.1], [a.sub.0] are nonnegative if and only if [R.sup.e.sub.0] < 1. Hence, all Routh stability criteria are satisfied; that is, the three eigenvalues of matrix (11) are negative or have negative real parts. Furthermore, for [R.sup.e.sub.0] = 1, (9) becomes the disease-free equilibrium. Therefore, the following result holds.

Theorem 3. The disease-free equilibrium [E.sup.0.sub.1] = (1,0, 0) exists and it is globally asymptotically stable if [R.sup.e.sub.0] [less than or equal to] 1.

Alternatively, the global stability of the disease-free equilibrium [E.sup.0.sub.1] can be established using the following Lyapunov function:

[mathematical expression not reproducible] (13)

where [u.sup.1], [u.sub.2], [u.sub.3] are some positive constants. Calculating the derivative of V along the solutions of system (1), we obtain

[mathematical expression not reproducible] (14)

[mathematical expression not reproducible] (15)

[mathematical expression not reproducible] (16)

Choosing [u.sub.1] = [u.sub.2] = b/{[[gamma].sub.h] + [[mu].sub.h]) and [u.sub.3] = 1, (16) becomes

[mathematical expression not reproducible] (17)

Thus, V' is negative for [R.sup.e.sub.0] [less than or equal to] 1. Note also that V' = 0 if and only if [S.sub.h] = [S.sup.0.sub.h] and [I.sub.h] = [I.sub.v] = 0. Therefore, the largest invariant set for (1) is the singleton {[E.sup.0.sub.1]}. Hence, by LaSalle's invariance principle [34], [E.sup.0.sub.1] is globally asymptotically stable when [R.sup.e.sub.0] [less than or equal to] 1 and Theorem 3 is valid.

Remark 4. Clearly, the endemic equilibrium [E.sup.*.sub.1] exists and is unique for [R.sup.e.sub.0] > 1. This excludes the possibility of occurrence of backward bifurcation. This result is of great epidemiological significance in guiding efforts for disease control as it indicates that [R.sup.e.sub.0] = 1 is the critical epidemic threshold.

To establish the local stability of the endemic equilibrium we evaluate the Jacobian of the system at [E.sup.*.sub.1], which gives

[mathematical expression not reproducible]. (18)

The characteristic polynomial of matrix (18) is then given by

[P.sub.2] ([lambda]) = [[lambda].sup.3] + [b.sub.2] [[lambda].sup.2] + [b.sub.1] [lambda] + [b.sub.0], (19)

where

[mathematical expression not reproducible] (20)

The coefficient [b.sub.2] > 0 and both [b.sub.1], [b.sub.0] are nonnegative if and only if [R.sup.e.sub.0] > 1. Hence, all Routh criteria are satisfied; that is, the three eigenvalues of matrix (11) are negative or have negative real parts. Therefore, the following results holds.

Theorem 5. The endemic equilibrium [E.sup.*.sub.1] = ([S.sup.*.sub.h], [I.sup.*.sub.h], [I.sup.*.sub.v]) exists and it is locally asymptotically stable if [R.sup.e.sub.0] > 1.

A global stability result for the endemic equilibrium [E.sup.*.sub.1] of system (1) is given below.

Theorem 6. If [R.sup.e.sub.0] > 1, the endemic equilibrium [E.sup.*.sub.1] = ([S.sup.*.sub.h], [I.sub.*.sup.h], [I.sub.*.sup.v]) is globally asymptotically stable.

Proof. [mathematical expression not reproducible] be the components of the Lyapunov function

U([S.sub.h], [I.sub.h], [I.sub.v]) = [D.sub.1] + [D.sub.2] + [D.sub.3], (21)

where [v.sub.1], [v.sub.2], [v.sub.3] are some positive parameters to be chosen later. Differentiating U along the solutions of system (1), we obtain

[mathematical expression not reproducible] (22)

[mathematical expression not reproducible] (23)

for [v.sub.2] = [v.sub.1] and [[gamma].sub.h] + [[mu].sub.h] = a[S.sup.*.sub.h][I.sup.*.sub.v]/[I.sup.*.sub.h] at equilibrium.

[mathematical expression not reproducible] (24)

for [v.sub.3] = 1/b and [[mu].sub.v](1 - [q.sub.v]) = [I.sup.*.sub.h] (1 - [I.sup.*.sub.v])/[I.sup.*.sub.v] at equilibrium. Now all together,

[mathematical expression not reproducible] (25)

Using the inequality 1 - x + ln x [less than or equal to] 0for x > 0 with equality holding if and only if x =1 and the fact that the arithmetic mean is greater than or equal to the geometric mean, we obtain U!([S.sub.h], [I.sub.h], [I.sub.v]) [less than or equal to] 0 for all [S.sub.h], [I.sub.h], [I.sub.v] > 0. Furthermore, we obtain that U'([S.sub.h], [I.sub.h], [I.sub.v]) = 0 holds only when [S.sub.h] = [S.sub.*.sup.h], [I.sub.h] = [I.sup.*.sub.h], [I.sub.v] = [I.sup.*.sub.v] and that [E.sup.*.sub.1] is the only equilibrium state of these systems on this plane (line). Therefore, by LaSalle's invariance principle [34], the positive equilibrium [E.sup.*.sub.1] is globally asymptotically stable.

3.1.2. Epidemic Thresholds, Vertical Infection, and Basic Properties. One of the most important critical thresholds in epidemic models is the basic reproductive number, [R.sup.0], which is usually found using the next-generation method, as the dominant eigenvalue of the next-generation matrix [35, 36]. Following the method in [36] we write system (1) consisting only of infectious compartments as the difference between new infection and transfer rates and the resulting Jacobian matrices evaluated at the disease-free equilibrium [E.sub.0] = (1,0,0). For system (1) we have two infected classes, namely, [I.sub.v] and [I.sub.h]. It follows that the transmission and the transfer matrices F and F, respectively, are defined as

[mathematical expression not reproducible] (26)

Unlike in host-vector models without vertical transmission, the diagonal elements of the transmission matrix F are nonzero. This stems from the fact that in the presence of vertical transmission there is vector to vector transmission, which completely changes the nature of the basic reproductive number. Thus, the next-generation matrix, K, is then given by

[mathematical expression not reproducible] (27)

and the resulting dominant eigenvalue of the spectral radius [FV.sub.-1], which is the basic reproductive number, [R.sup.0], is given by

[mathematical expression not reproducible] (28)

When there is no vertical transmission, qV = 0, as in the case of malaria, [R.sup.0] = [R.sup.H.sub.0] is simply the geometric mean of the product of the number of new infections in hosts from one infected vector and the number of new infections in vectors from one infected host, in the limiting case that both populations are fully susceptible. The interpretation and epidemiological significance of [R.sup.H.sub.0] are well established. It is easy to see that the transmission of infection is increased with efficiency of vector biting a and probabilities of successful infection transmission, but it is hindered by high mosquito death rates and faster host recovery. The biting rate [alpha] appears as [[alpha].sup.2] because it enters twice in the transmission cycle [31]. An important parameter is the ratio of female mosquitoes to hosts m = [N.sub.v]/[N.sub.h], which is central for disease spread according to model settings. [R.sup.0] increases with the number (or density) of mosquitoes but decreases with the number (or density) of host population. This results from the asymmetry in the dependence of the vector's biting rate on the sizes of the host and vector populations such that when there are many more hosts compared to mosquitoes, sustained transmission maybe impossible. Therefore, in the absence of vertical transmission, for the infection to successfully spread and invade, the ratio of mosquitoes to hosts needs to be sufficiently large so that double bites are common [37]:

[mathematical expression not reproducible] (29)

where the critical ratio is given by [mathematical expression not reproducible]. Note that each mosquito could infect less than one host on average, and yet [R.sup.0] could still be more than unity. To elucidate this fact we write [R.sup.H.sub.0] as a product of each host type single-step reproductive number; that is,

[mathematical expression not reproducible]. (30)

where [R.sub.vh] represents the number of new infections in mosquitoes from a single infected host while [R.sub.hv], represents the number of new infections in hosts from a single infected mosquito. Clearly, [R.sup.H.sub.0] can be greater than unity even when one of these reproduction numbers is less than unity, and it can also be less than unity even if one of its components is greater than unity. In Figure 1(a) we depict contours plots corresponding to the overall [R.sup.0] in (31) along the plane ([R.sub.hv], [R.sub.vh]) in the presence of vertical transmission. Its effects in this asymmetric relationship between [R.sup.0] and its components is not very pronounced and the ratio of mosquitoes to hosts remains one of the leading factors when there is a large disparity between the sizes of the host and vector populations.

Clearly, the geometric mean is less than the average expected number of new infections per generation. This is the case where the next-generation method fails to produce the correct [R.sup.0] if transmission between hosts is intermediated by another host (for more discussion about the failure of this method see [21]). Instead, it gives the weighted average lying between the number of new infections each individual produces in the next infection event. Note that if the number of hosts is increased, this deficit is compounded. This can have serious implications for guiding disease control efforts as it fails to provide the actual severity of the infection. Therefore, in this study for application purpose we propose the use of another epidemic threshold. This is [R.sup.0] from (10) derived in Section 3.1.1, which is hereby referred to as the 'effective' reproductive number:

[R.sup.e.sub.0] = 1/1 - [q.sub.v] [p.sub.vh] [P.sub.hv][[alpha].sup.2]m/[[mu].sub.v] ([[gamma].sub.h] + [[mu].sub.h], (31)

as (1) it satisfies the property that the endemic equilibrium [E.sup.*.sub.1] only persists if [R.sup.e.sub.0] is greater than unity and (2) the endemic equilibrium exists without occurrence of backward bifurcation, meaning that [R.sup.e.sub.0] = 1 is the correct critical value. At this point, it is of particular interest to establish the critical ratio of mosquitoes to hosts in the settings of vertical transmission. From (31) we obtain

[mathematical expression not reproducible] (32,

that is, a new critical ratio [m.sup.*.sub.c] = (1 - [q.sub.v]) [m.sub.c] < [m.sub.c] for 0< [q.sub.v] <1. This result indicates that, in the presence of vertical transmission, the critical ratio of female mosquitoes to hosts reduces at a rate proportional to vertical transmission efficiency. This indicates an earlier occurrence of the point where an epidemic is just possible as above this level the equilibrium prevalence is expected to rapidly increase to its asymptotic value.

Others are epidemic thresholds that provide a direct measure of the control effort required for disease eradication. These measures are known as "type" reproductive numbers and are determined based on the next-generation matrix [38, 39] such that for our case we have host and vector type reproductive numbers [T.sup.h.sub.1] and [T.sup.v.sub.1], respectively. Note that in the absence of vertical transmission the host and vector type reproductive numbers are both equal to the square of the basic reproductive number, [R.sup.0] in (41), which turns out to be [R.sup.e.sub.0] in (31) for [q.sub.v] = 0. In the presence of vertical infection, [R.sup.2.sub.0] gives the expected number of secondary infections [R.sup.2.sub.0] after one average, complete (host-vector-host or vector-host-vector) transmission cycle but does not correspond to a specific population type [23]. For the case of host population, one infected host leads to some secondary host infections in the next host-vector-host transmission cycle. This results exactly from the horizontal transmission mechanism and it is given by [R.sup.H.sub.0]. However, further secondary host infections may also occur after any number of vector-vector transmission cycles as a result of transovarial transmission. From the expression of [R.sup.0] in (31) it can be seen that new infections resulting from vector-vector transmission cycle are given by 1/(1-[q.sub.v]) ;hence the host type reproductive number is then given by

[T.sup.h.sub.1] = 1/1 - [q.sub.v] [R.sup.H.sub.0] = [R.sup.e.sub.0]. (33)

From (33), [R.sub.H.sup.0] = [R.sub.e.sup.0] if 1/(1 - [q.sub.V])=1; that is, [q.sup.*.sub.v] =0 which is the boundary condition. An important question would be, what is the proportion of vertical infection that doubles [R.sup.H.sub.0]? For 1/(1 - [q.sub.v]) = 2 we obtain [q.sup.*.sub.v] = 0.5, which means that if about 50% of infected mosquitoes produce infected offspring, we expect [T.sup.h.sub.1] to double. Note that all secondary infections resulting from a single infected vector must occur in the next vector-host-vector transmission cycle; hence, the vector type reproductive number takes the form

[T.sup.V.sub.1] = [R.sup.H.sub.0] + [q.sub.v]. (34)

For details on the derivation of both types of reproductive numbers, see Supplementary Materials, Section A-1. A number of methods for deriving [R.sup.0] exist and in [21] it is demonstrated that the resulting [R.sup.0] is not unique but they may have a common threshold at [R.sup.0] = 1. Hence, their behaviour below and above unity may show significant differences. In Figure 1[C] we examine how vertical transmission efficiency contributes to the qualitative behaviour of [R.sup.2.sub.0], [T.sup.h.sub.1], and [T.sup.v.sub.1] above the unity threshold criterion. An important characteristic pattern of their relationship is that below a certain percentage of vertical transmission efficiency the three reproductive numbers are undistinguishable but beyond the relationship between vertical infection and host type reproductive number becomes exponential. This behaviour results from the fact that in each generation, the number of host infections is proportional to the number of infected mosquitoes, which decreases proportionally to the vertical infection efficiency. Hence, for low vertical infection efficiency, possibly the first generations make a significant contribution to subsequent transmission cycles (for more discussion see [23, 27]).

In Figure 1(b) we compare the basic reproductive number [R.sup.0] derived through the next-generation method with the one derived from the endemic equilibrium [R.sup.e.sub.0] with respect to vertical transmission efficiency. Clearly, the two give different relationships and it is easy to see that the effect of vertical infection is completely diluted on [R.sup.0]. Hence, this measure if used for disease control effort is more likely to provide misleading guidance. The relationship of the basic reproductive number [R.sup.e.sub.0] with mosquito to host ratio m and vertical infection [q.sub.v] shows some interesting patterns. For m > 1 the relationship between [R.sup.e.sub.0] and [q.sub.v] is exponential as described in Figure 1(c); however for m < 1 it is linear rendering [q.sub.v] having little effect if not negligible for [q.sub.v] [member of] (0,0.8) (see Figure 1(d) dashed lines). This suggests that when there are more hosts compared to mosquitoes, vertical infection efficiency becomes completely negligible as in the model settings, and if m < 1, sustained host-vector transmissions may be impossible. On the other hand, we observe that for [q.sub.v] [much greater than] 0.8 and even for m < 1, [R.sup.e.sub.0] becomes much larger highlighting the importance of the interplay between vertical infection and the ratio of mosquitoes to hosts.

Remark 7. In the settings of vector-borne diseases both host type reproductive number [T.sup.h.sub.1] or ([R.sup.e.sub.0]) and vector type reproductive number [T.sup.v.sub.1] are critical for providing useful informative indicator for guiding disease control efforts.

3.2. Effects of Ratio of Mosquitoes to Hosts and Vertical Transmission on Disease Prevalence. To fully capture effects of vertical transmission efficiency on vector-borne transmission models, we explore how this parameter influences the disease spread patterns.

Results are shown in Figure 2. The left figure explores how the first peak of an epidemic denoted by max([I.sub.h](t)) varies with vertical transmission efficiency, [q.sub.v]. The results show that max([I.sub.h](t)) increases linearly with [q.sub.v]. The righth-and figure plots the changes on time of occurrence of each epidemic peak at each value of vertical transmission efficiency. The results suggest that as vertical transmission efficiency increases the time position of each epidemic peak decreases linearly. This suggests that if vertical infection efficiency is increased, the first peak of an epidemic is likely to be attained earlier than the expected time. This suggests that while vertical transmission efficiency increases the epidemic size of an outbreak, it linearly decreases the time of the peak of an outbreak. The lower figures, (c) and (d), describes how the prevalence of infected hosts, [I.sub.h], varies with the ratio of female mosquitoes to hosts for different values of vertical transmission efficiency. Figure 2(c) indicates that the prevalence of the disease saturates at larger values of the ratio at earlier stage of the initial epidemic while Figure 2(d) indicates that the prevalence of the disease saturates at lower values of the ratio at later stage of the initial epidemic. This result stems from the fact that at earlier stage of the initial epidemic the size of the vector population is still very high compared to later stages of evolution of the epidemic.

3.3. Analysis of the Model with Vector Stages and Vertical Transmission

3.3.1. Model Equilibria. Model equilibria are obtained by setting the right-hand side of system (2) equal to zero. The model has three equilibria, two disease-free and one endemic equilibrium. Details about the derivation of the components of the endemic equilibrium are given in Section A-2 of Supplementary Materials. The first trivial equilibrium is both mosquito-free and disease-free:

[mathematical expression not reproducible], (35)

which corresponds only to the presence of hosts (this can be livestock for the case of RVF). However, this equilibrium is not biologically feasible or relevant in particular in tropical regions where vector-borne diseases such as Dengue and RVF are endemic. The second equilibrium corresponds to the coexistence of both vectors and hosts but without infection. This is the disease-free equilibrium:

[mathematical expression not reproducible], (36)

which can be used for computation of the basic reproduction number when using the next-generation method [36]. The nonzero components are

[mathematical expression not reproducible] (37)

Hence, [X.sup.1] is biologically significant if [Q.sub.0] > 1; that is, the mosquito population exists and can establish itself if on average each adult female mosquito produces at least one female mosquito at oviposition rate r during her life time period 1/[[mu].sub.v] after successfully surviving the aquatic stage at probability [theta]/([theta] + [micro]L].).

The endemic equilibrium which represents the prevalence of the disease in both host and vector populations is given by

[x.sup.*] = ([L.sup.*.sub.s], [L.sup.*.sub.i], [L.sup.*.sub.v], [L.sup.*.sub.h], [L.sup.*.sub.h]), (38)

[mathematical expression not reproducible] (39)

[mathematical expression not reproducible] (40)

where

[mathematical expression not reproducible] (41)

Similar to the discussion in Section 3.1.2 the above epidemic threshold in (41) is the "effective" or "control" reproductive number since it satisfies the property that the endemic equilibrium E* only persists if [R.sup.c.sub.0] is greater than unity and there is no backward bifurcation.

Remark 8. The relationship between the endemic equilibrium and [R.sub.0] which depends on the vector population threshold indicates that the existence of this equilibrium is governed by the norms of the vector population threshold. Therefore, this vector population threshold is an important parameter in vector control efforts as it provides a means for identifying key factors for reducing the vector population.

In the following section we derive [R.sup.0] using the next generation method.

3.3.2. The Basic Reproductive Number and Other Control Thresholds. The basic reproductive number, [R.sub.0], is a concept of great epidemiological significance. Following the method of [36], we write system (2) considering only disease compartments, [mathematical expression not reproducible], where

[mathematical expression not reproducible] (42)

Evaluating the Jacobian of above matrices at the disease-free equilibrium [X.sup.1] = (Z*,0, P*, 0,1, 0), the basic reproduction number is then the spectral radius of the next-generation matrix; that is, [R.sup.0] = [rho](F([X.sup.1])V [([X.sup.1]).sup.-1]),

[mathematical expression not reproducible] (43)

Following the discussion in Section 3.1.2 this epidemic threshold (43) can lead to an underestimate of the control effort required if used to guide disease control strategies. To overcome this inconsistency we use type reproductive numbers to estimate efforts required to provide informative indicators in control of vector-borne diseases. To avoid repetition we omit the steps required in the derivation process, such that the host type reproductive number is

[T.sup.h.sub.1] = [R.sup.H.sub.0]/1 - [q.sub.v] = [R.sup.c.sub.0] (44)

and the vector type reproductive number is then given by

[T.sup.v.sub.1] = [R.sup.H.sub.0] + [q.sub.v]. (45)

In the majority of mosquito-borne diseases, treatment or adequate vaccine is not likely to be feasible, either due to financial constraints or nonavailability, in particular in African communities with low resilience to economic challenges. In such cases the main preventive measures are individual protection against mosquito bites and the control of the proliferation of both the larval and adult mosquitoes. This can be achieved by lowering appropriate epidemic thresholds below unity, that is, forcing the reproductive number of the pathogen below 1. In the context of our model this implies reducing either of or both type reproductive numbers ([T.sup.h.sub.1] and/or [T.sup.v.sub.1]) below unity. Hence, we define the eradication effort to be the percentage reduction in the mosquito population size required to prevent endemic transmission of the pathogen [23], such that if using [T.sup.h.sub.1] the eradication effort is 100(1 - 1/[T.sup.h.sub.1]) = 100([R.sup.H.sub.0] + qv - 1)/[R.sup.H.sub.0] and for [T.sup.v.sub.1] the effort required to eradicate the disease is [mathematical expression not reproducible]. Therefore, vertical infection at [q.sub.v]% can be responsible for at most [q.sub.v]% of the required eradication effort, provided [R.sup.H.sub.0] > 1.

3.3.3. Stability of the Host-Only System. This fixed point [mathematical expression not reproducible] is an equilibrium point of the host-only system given by

d[S.sub.h]/dt = [[mu].sub.h] (1 - [S.sub.h)], (46)

which is local and globally asymptotically stable for all t > 0 as the derivative ([[mu].sub.h](1 - [S.sub.h]))' is negative.

3.3.4. Local Stability of the Disease-Free Equilibrium. The local stability of the disease-free equilibrium point [X.sup.1] = (l*,0,P*,0,1,0) is established by analysing the eigenvalues of the Jacobian matrix at [X.sup.1]:

[mathematical expression not reproducible] (47)

Direct computation shows that [[lambda].sub.1] = -[[mu].sub.h] is an eigenvalue of matrix (47) and the remaining are solutions of equation

[mathematical expression not reproducible] (48)

where

[mathematical expression not reproducible] (49)

By Descartes' rule of signs, since [c.sub.4] >0, [c.sub.3] > 0, [c.sub.2] > 0, and [c.sub.1] > 0, then all roots of (48) are negative or have a negative real part for [c.sub.0] > 0. The coefficient [c.sub.0] is nonnegative if and only if [R.sup.c.sub.0] < 1; hence the following result holds.

Theorem 9. The disease-free equilibrium point [X.sup.1] = (l*,0, P*, 0, 1, 0) is locally asymptotically stable whenever [R.sup.c.sub.0] <1.

3.3.5. Global Stability of the Disease-Free Equilibrium. To establish the global asymptotic stability (GAS) of [X.sup.1] we use results obtained by Kamgang and Sallet [40], which are the extension of some results in [36]. Around the disease-free equilibrium system (2) can be written as

[mathematical expression not reproducible] (50) where [x.sub.s] is the vector representing disease-free compartments ([L.sub.s], [S.sub.v], [S.sub.h]) and the vector [x.sub.I] represents the state of infected compartments ([L.sub.i],[I.sub.v], [I.sub.h]). This requires rewriting system equation around [X.sup.-1], as follows:

[mathematical expression not reproducible] (51)

Note that the latter arises from ([theta] + [[mu].sub.L])[L.sup.1.sub.s] = [delta]r(1 - [L.sup.1.sub.s]) [S.sup.1.sub.v]. Then, the following matrices are obtained,

[mathematical expression not reproducible] (52)

[mathematical expression not reproducible] (53)

From basic matrix properties and direct algebraic computation it follows that all eigenvalues of matrix [A.sub.1] are real and negative and both [A.sub.1] and [A.sub.2] are Metzler matrices. Hence, system [mathematical expression not reproducible] is GAS at the disease-free equilibrium [mathematical expression not reproducible]. To establish the global stability of overall system (50) at [X.sup.1] conditions of the following theorem in [40] must be satisfied.

Theorem 10. Let [PHI] [subset] U = [R.sup.3.sub.+] x [R.sup.3.sub.+]. System (50) is of class [C.sup.1], defined on U if

(1) U is positively invariant relative to (50).

(2) The system [mathematical expression not reproducible] is GAS at [x.sub.DPE,S].

(3) For any x [member of] [PHI], matrix [A.sub.2](x) is Metzler irreducible.

(4) There exists a matrix [[bar.A].sub.2], which is an upper bound of the set. M = {[A.sub.2](x) [member of] [M.sub.3],(R) | x [member of] [bar.[PHI]]}, with the property that if [[bar.A].sub.2] [member of] M, for any [bar.x] [member of] [[bar.sub.[PHI]], such that [A.sub.2] ([bar.x]) = [[bar.A].sub.2], then [bar.x] [member of] [R.sup.3] x {0}.

(5) The stability modulus of [mathematical expression not reproducible]), satisfies [alpha] ([bar.A].sub,.2) [less than or equal to] 0.

Then, [X.sup.1] is GAS in [bar.[PHI]].

Proof. Clearly, conditions (1)-(3) of the theorem have been satisfied. For all x [member of] [PHI], [A.sub.2] (x) is irreducible because [(I + |[A.sub.2](x)).sup.2] > 0. An upper bound of the set of matrices M, which is the matrix [[bar.A].sub.2], is given by matrix [A.sub.2]([bar.x]), where [mathematical expression not reproducible]. Similarly matrix [[bar.A].sub.2] is irreducible. Recall that from the Perron-Frobenius theorem for an irreducible matrix you get that one of the matrix eigenvalues is positive and greater than or equal to all others, that is, the dominant eigenvalue. Matrix [A.sub.2] is exactly the matrix used to compute the basic reproductive number, i.e., the dominant eigenvalue. For more details or proof in general settings see [40].

Now conditions (1)-(4) have been verified. To check the last condition, we make use of the following Lemma [40].

Lemma 11. Let H be a square Metzler matrix written in block form [mathematical expression not reproducible], with A and D square matrices. H is Metzler stable if and only if matrices A and D - [CA.sup.-1] B are Metzler stable.

Matrix [A.sub.2](x) in block matrices takes the following components:

[mathematical expression not reproducible] (54)

Clearly, A is a Metzler stable matrix and

[mathematical expression not reproducible] (55)

is also Metzler stable if [mu].sub.v] (1 - [q.sub.v])([[lambda].sub.h] + [[mu].sub.h]) - [[lambda].sub.hv] [[lambda].sub.h] [P.sup.*] [greater than or equal to] 0, that is, 1 - [R.sup.C.sub.0] [greater than or equal to] 0. Therefore, from Theorem 10 and Lemma 11 the following result holds.

Theorem 12. For [Q.sub.0] > 1 there exists a unique disease-free equilibrium [X.sup.1] = ([l.sup.*], 0, [P.sup.*],0,1,0), which is globally asymptotically stable whenever [R.sup.c.sub.0] [less than or equal to] 1.

Remark 13. The above result is of great epidemiological importance. It highlights two fundamental indicators: (1) Although vector control is central for control of vector-borne diseases, it does not mean eliminating all the vectors. Note that the effect of this include both infected and non-infected mosquitoes. (2) Disease eradication efforts are independent of the initial sizes of both vector and host populations; however, the ratio between the two populations is a key factor.

3.3.6. Local Stability of the Endemic Equilibrium. Results of the local stability of the disease-free equilibrium [X.sup.1] suggest that for [R.sup.c.sub.0] = 1 the Jacobian matrix (47) has zero eigenvalue while the remaining eigenvalues are negative or have a negative real part. Since the algebraic computation involved when establishing the stability of the endemic equilibrium through linearization is quite extensive, we employ the centre manifold theory [41]. This theory is used to examine existence of backward or forward bifurcation. The bifurcation occurs at [R.sup.c.sub.0] = 1, and choosing [q.sub.v] = [q.sup.*.sub.v] as abifurcation parameter, then

[mathematical expression not reproducible] (56)

The Jacobian matrix J([X.sup.1], [q.sup.*.sub.v]) is the same as matrix (47); hence, J([X.sup.1], [q.sup.*.sub.v]) has a simple zero eigenvalue [lambda] = 0 when [R.sup.c.sub.0] = 1 and others given as [lambda] = -[[mu].sub.h] and roots of the polynomial [[lambda].sup.4] + [c.sub.4] [[lambda].sup.3] + [c.sub.3] [[lambda].sup.2] + [c.sub.2][lambda] + [c.sub.1] = 0, where [c.sub.1], [c.sub.2], [c.sub.3], [c.sub.4] are as defined in Section 3.3.4. Let w = ([[omega].sub.1], [[omega].sub.2], [[omega].sub.3], [[omega].sub.4], [[omega].sub.5], [[omega].sub.6]) be the right eigenvector associated with zero eigenvalue of the Jacobian matrix J([X.sup.1], [q.sup.*.sub.v]). Its components are derived by solving J([X.sup.1], [q.sup.*.sub.v])x w = 0, which gives

[mathematical expression not reproducible] (57)

where

[mathematical expression not reproducible] (58)

Clearly, [F.sup.*] and [G.sup.*] are nonnegative if and only if [q.sup.*.sub.v] [greater than or equal to] 1/[Q.sub.0], [Q.sup.v.sub.0] > 1, and [Q.sub.0] > 1. Note that [Q.sup.v.sub.0] gives the average number of emerged infected female mosquitoes produced by one infected female mosquito. Hence, for [Q.sup.v.sub.0] >1 infected mosquitoes will persist in the population, leading to disease persistence. Let us denote v = ([v.sub.1], [v.sub.2], [v.sub.3], [v.sub.4], [v.sub.5], [v.sub.6]) as the left eigenvector associated with zero eigenvalue of the Jacobian matrix J ([X.sup.1], [q.sup.*.sub.v]). Its components are obtained when solving [J([X.sup.1], [q.sup.*.sub.v])].sup.T] x w = 0,(T represents transpose), which gives

[mathematical expression not reproducible] (59)

where [M.sup.*] = [theta][[mu].sub.v]([delta]r + [theta] + [[mu].sub.L])/[delta]r([theta] + [[mu].sub.L])([theta] + [delta][[mu].sub.v]).

Computation of Bifurcation Parameters a* and b*. Because [v.sup.1] = [v.sup.3] = [v.sup.5] = 0 the associated nonzero partial derivatives of [f.sub.1], [f.sub.3], [f.sub.5] are not considered and the remaining nonzero partial derivatives of system [mathematical expression not reproducible].

Thus, [mathematical expression not reproducible] after some algebraic simplification gives

[mathematical expression not reproducible] (60)

Given that the vital dynamics of the host are much more slower [mathematical expression not reproducible] and using the property that vw = 1, then a* < 0 for [w.sup.4] > 0. For b* the nonzero partial derivative of [mathematical expression not reproducible]. Using the property that [v.sub.w] = 1 we obtain that b* > 0. Since a* <0 for w4 > 0 and b* > 0, the model (2) exhibits a forward bifurcation at [R.sup.c.sub.0] = 1. Therefore, the following result holds.

Theorem 14. The endemic equilibrium point X* = ([L.sup.*.sub.s], [L.sup.*.sub.i], [S.sup.*.sub.v], [I.sup.*.sub.v], [S.sup.*.sub.h], [I.sup.*.sub.h]) exists and is locally asymptotically stable for [R.sup.c.sub.0] > 1.

Remark 15. The condition [Q.sup.v.sub.0] >1 is an indication that vertical transmission is vital for disease long-term persistence. A particular example is the case of Dengue fever and RVF. Upon failure of this condition, it is much more likely that results of Theorem 14 do not hold. Another important feature of [Q.sup.v.sub.0] is that vertical infection is proportional to [Q.sub.0], and hence vector control is still a viable control strategy even for vector-borne disease with transovarial transmission.

The above stability results and those in Theorems 9, 10, 12, and 14 can be summarized in Table 1 and in the bifurcation diagram in Figure 3.

3.4. Numerical Results of the Model with Vector Stages and Vertical Transmission. Numerical analysis using reasonable parameter values for RVF is carried out. In Figure 4 we plot time series of infected larval, adult mosquitoes and hosts, and their respective phase portraits. The latter are plotted for t [member of] [0 50] years for observing long-term dynamics. The parameter values used give [R.sup.c.sub.0] >1 and existence of a stable endemic equilibrium. These results not only confirm findings by analytical analysis (Theorems 9 and 12) but suggest that this equilibrium is also global asymptotic stable for the chosen parameter values. This also implies that the endemic equilibrium is stable not only close to the bifurcation point but also for Rc0 beyond unity (see inner figures, Figures 4(a)-4(c)).

4. Discussion and Conclusion

One of the main contributions of mathematical epidemic models is that they enhance our understanding of disease transmission and public health planning. This is possible because with models we are able to derive critical epidemic thresholds in terms of model parameters which carry important disease features and key players. One of these measures is the basic reproduction number [R.sup.0], which if greater than unity implies that the disease will persist in the population and fade out otherwise. However, a major challenge is to derive the correct [R.sup.0] for disease models in which infection transmission is intermediated by another host. This is the case of vector-borne diseases, where practical methods such as the next-generation matrix for deriving [R.sup.0] fails to give the correct value [21]. Rather, it gives the geometric mean of secondary infection per generation. Using the basic model we discussed alternative measures, known as host and vector type reproductive numbers ([T.sup.h.sub.1] and [T.sup.v.sub.1], respectively) [38, 39], to be used for guiding disease control efforts. Further, for both the basic and the complex model we showed that the basic reproductive number derived from the existence of endemic equilibrium is the correct epidemic measure and it is equivalent to the host type reproductive number [T.sup.h.sub.1]. In fact, this is an alternative method for finding [R.sup.0] when it is possible to derive the explicit components of the endemic equilibrium point [21, 42]. Then we established important relationships between the reproductive numbers [R.sup.2.sub.0], [T.sup.h.sub.1], [T.sup.v.sub.1] and vertical infection efficiency. It follows that there is a critical threshold of vertical infection from which the contribution of this transmission route becomes significant. Otherwise, if vertical infection efficiency is very low, its contribution to disease persistence may be insignificant. The work in [23, 27] discusses this feature in more detail and their findings suggest that these results stem from the fact that for low vertical infection efficiency the virus is rapidly lost unless there is a regular disease amplification in the host population. Recent entomological studies have confirmed this situation in particular for RVF where there has been more and more evidence of disease interepidemic activities [7,1720], underlined by sporadic outbreaks on livestock at low levels. During this period vertical infection plays a key role for providing protection against chance extinction, but only if it is efficient enough to span several mosquito generations [23], since at the end of unfavourable season the horizontally transmitted epidemic is declining.

Also, our results demonstrated that while the relationship between [R.sup.0] and vertical infection efficiency is linear, the relationship between the correct threshold [R.sup.e.sub.0] and vertical infection efficiency is exponential. This shows that the failure of [R.sup.0] as a critical epidemic threshold compounds itself when varying vertical transmission. Therefore, more caution is needed when deriving the basic reproductive number in vector-borne disease systems that include vertical transmission as one of the transmission modes. Another important property to highlight in our models is the female mosquitoes to host ratio which is more likely to be responsible for the asymmetric relationship between the host-to-vector and vector-to-host reproductive numbers as discussed in other studies [28, 43]. We found that in the presence of vertical infection the critical ratio of female mosquitoes to hosts decreases linearly when increasing vertical transmission efficiency. This result indicates an earlier occurrence of the point where an epidemic is just possible as above this level the equilibrium prevalence is expected to rapidly increase to its asymptotic value [37]. Further we analysed how vertical infection affects the prevalence of the disease at the first peak of an epidemic. The results showed that the prevalence at the first peak is positively correlated to vertical transmission efficiency while the time position of occurrence of each peak is negatively correlated to vertical transmission efficiency. These findings suggest that while vertical transmission efficiency increases the epidemic size of an outbreak, it reduces the duration of an outbreak. Additionally, the relationship between the prevalence of the disease and ratio of female mosquitoes to hosts was accessed for different values of vertical infection. Analysis showed that there is an optimal value of the ratio above which the prevalence of the disease decreases. However, such optimal value changes according to the stage of epidemic development, such that at early stage of the epidemic the optimal value is attained much later than that at later stage of the epidemic. Another interesting result is that after the optimal value of the ratio the effect of vertical transmission changes where higher values have less effect compared to lower values of vertical transmission efficiency.

Extending the basic model to include larval compartments with maximal capacities for both larval and adult populations we derived important critical thresholds for both vector population and disease system. We found that the model exhibits three boundary equilibria, namely, a mosquito-free system, a mosquito-host system but disease-free, and one where the disease persists in both populations. A detailed stability analysis is presented and numerical simulations are conducted using parameter values relevant to Rift Valley fever. The mathematical results are used to interpret the biological implications of the relationships between [R.sup.0], vertical infection, and the ratio of female mosquitoes to hosts for assessing viable control measures. Our analysis shows that the qualitative behaviours of the system are completely determined by three key quantities: [Q.sub.0], the vector population threshold; [Q.sup.v.sub.0], the average number of female mosquitoes produced by a single infected mosquito; and [R.sup.c.sub.0], the effective reproductive number (see the bifurcation diagram in Figure 3 and Table 1). These results provide important qualitative understanding of the interaction between vertical infection and ratio of female mosquitoes to hosts on the prevalence of the disease.

The mosquito population persists if each female adult mosquito produces at least one larval mosquito after surviving the larval stage. Then, if an infected individual is introduced in the community without any form of protection, the disease will invade the system leading to appearance of an endemic equilibrium if the epidemic criterion is satisfied ([R.sup.e.sub.0] > 1). The biological implication of this result is that the disease can only invade an area if already colonized by the mosquito population as expected. Hence, measures targeting critical features of the life cycle of the vector are viable disease control strategies. Furthermore, we found that the persistence of the disease is subjected to vertical transmission efficiency governed by the epidemic criterion [Q.sup.v.sub.0] > 1. This implies that for the disease to persist for each larval mosquito produced, 1/(1 - [q.sub.v]) should be infected. This explicit formulation of the threshold conditions in terms of parameters governing the infection transmission process and vector life history is of great epidemiological significance as it allows disease control efforts to be targeted at specific disease and vector stages. We have also shown that the disease-free equilibrium is locally asymptotically stable for [R.sup.c.sub.0] <1 and globally asymptotically stable if [R.sup.c.sub.0] [much less than] 1. The latter result is of outmost importance, because it shows that if at any time, through appropriate interventions (e.g., destruction of breeding sites, use of insecticides and repellents on the host, or even vaccination), we are able to lower [R.sup.c.sub.0] below 1, then the disease will disappear. When the vertical transmission route is possible, this requires also lowering the contribution from this transmission mode. That is, lowering [Q.sup.v.sub.0] which implies keeping qv very close to zero. Therefore, the design of control programs should take into account the implications of this mode of transmission in the case of vector-borne diseases that include vertical transmission in the vector population.

The goal of this study was to use mathematical models for discussing important properties of vector-borne disease systems and their implications in attempts for control strategies. Though we have used simple but realistic models, we believe that they remain valuable tools useful in the context of research and for providing qualitative understanding of complex processes underlying vector-borne disease transmission systems in particular in the context of transovarial transmission in the vector population. In conclusion, our results show that (1) including larval stages in the vector population when modelling vector-borne disease systems has benefits, as it allows determining conditions for colonization in terms of vector life cycle features; (2) the failure of [R.sup.0] derived from the next-generation method compounds itself in the presence of vertical transmission; (3) host type reproductive number gives the correct [R.sup.0] useful for guiding disease control strategies; (4) while vertical transmission efficiency is likely to increase the size of an epidemic, it decreases its duration; and (5) in the presence of vertical infection the critical ratio of female mosquitoes to hosts reduces linearly when increasing vertical transmission efficiency.

https://doi.org/10.1155/2018/2618985

Data Availability

The relevant data used to support the findings of this study are included within the article.

Conflicts of Interest

The author declares that he has no conflicts of interest.

Supplementary Materials

Supplementary material contains details on the analytical derivation of the basic reproductive number, type reproductive numbers, and model equilibria. Table 1: parameter values used in the simulations. (Supplementary Materials)

References

[1] D. J. Gubler, "Vector-borne diseases," Revue Scientifique et Technique de l'OIE, vol. 28, no. 2, pp. 583-588, 2009.

[2] World Health Organization, WHO--Vector-Borne Diseases, 2016.

[3] A. Huppert and G. Katriel, "Mathematical modelling and prediction in infectious disease epidemiology," Clinical Microbiology and Infection, vol. 19, no. 11, pp. 999-1005, 2013.

[4] F. Brauer, "Mathematical epidemiology: Past, present, and future," Infectious Disease Modelling, vol. 2, no. 2, pp. 113-127, 2017

[5] C. E. Walters, M. M. Mesle, and I. M. Hall, "Modelling the global spread of diseases: A review of current practice and capability," Epidemics, 2018.

[6] R. M. Murithi, P. Munyua, P M. Ithondeka et al., "Rift Valley fever in Kenya: History of epizootics and identification of vulnerable districts," Epidemiology and Infection, vol. 139, no. 3, pp. 372-380, 2011.

[7] C. Sindato, E. D. Karimuribo, D. U. Pfeiffer et al., "Spatial and temporal pattern of rift valley fever outbreaks in Tanzania; 1930 to 2007," PLoS ONE, vol. 9, no. 2, 2014.

[8] K. J. Linthicum, F. G. Davies, A. Kairo, and C. L. Bailey, "Rift Valley fever virus (family Bunyaviridae, genus Phlebovirus). Isolations from diptera collected during an inter-epizootic period in Kenya," Epidemiology & Infection, vol. 95, no. 1, pp. 197-209, 1985.

[9] W. Horsfall, H. Fowler, M. Lj, and L. Jr, Bionomics and Embryology of the Inland Flood Water Mosquito Aedes Vexans, University of Illinois Press, 1974.

[10] F. F. Luduena Almeida and D. E. Gorla, "The biology of Aedes (Ochlerotatus) albifasciatus Macquart, 1838 (Diptera:

Culicidae) in central Argentina," Memorias do Instituto Oswaldo Cruz, vol. 90, no. 4, pp. 463-468,1995.

[11] N. Becker, D. Petric, M. Zgomba et al., Mosquitoes and Their Control, Springer-Verlag Berlin Heidelberg 2003, 2010.

[12] K. J. Linthicum, T. M. Logan, C. L. Bailey, D. J. Dohm, and J. R. Moulton, "Transstadial and horizontal transmission of Rift Valley fever virus in Hyalomma truncatum," The American Journal of Tropical Medicine and Hygiene, vol. 41, no. 4, pp. 491-496, 1989.

[13] X. Guo, T. Zhao, Y. Dong, and B. Lu, "Survival and replication of dengue-2 virus in diapausing eggs of Aedes albopictus (diptera: Culicidae)," Journal of Medical Entomology, vol. 44, no. 3, pp. 492-497, 2007.

[14] B. Angel and V. Joshi, "Distribution and seasonality of vertically transmitted dengue viruses in Aedes mosquitoes in arid and semi-arid areas of Rajasthan, India," Journal of Vector Borne Diseases, vol. 45, no. 1, pp. 56-59, 2008.

[15] M. Pepin, M. Bouloy, B. H. Bird, A. Kemp, and J. Paweska, "Rift Valley fever virus (Bunyaviridae: Phlebovirus): An update on pathogenesis, molecular epidemiology, vectors, diagnostics and prevention," Veterinary Research, vol. 41, no. 6, 2010.

[16] Y. Ba, D. Diallo, C. M. F. Kebe, I. Dia, and M. Diallo, "Aspects of bioecology of two rift valley fever virus vectors in Senegal (West Africa): Aedes vexans and Culex poicilipes (Diptera: Culicidae)," Journal of Medical Entomology, vol. 42, no. 5, pp. 739-750, 2005.

[17] A. D. Labeaud, E. M. Muchiri, M. Ndzovu et al., "Interepidemic Rift Valley fever virus seropositivity, northeastern Kenya," Emerging Infectious Diseases, vol. 14, no. 8, pp. 1240-1246,2008.

[18] R. D. Sumaye, E. Geubbels, E. Mbeyela, and D. Berkvens, "Inter-epidemic transmission of rift valley fever in livestock in the Kilombero river valley, Tanzania: a cross-sectional survey," PLOS Neglected Tropical Diseases, vol. 7, no. 8, Article ID e2356, 2013.

[19] N. O. Owange, W. O. Ogara, H. Affognon et al., "Occurrence of rift valley fever in cattle in Ijara district, Kenya," Preventive Veterinary Medicine, vol. 117, no. 1, pp. 121-128, 2014.

[20] J. K. Lichoti, A. Kihara, A. A. Oriko et al., "Detection of rift valley fever virus interepidemic activity in some hotspot areas of kenya by sentinel animal surveillance, 2009-2012," Veterinary Medicine international, vol. 2014, Article ID 379010, 9 pages, 2014.

[21] J. Li, D. Blakeley, and R. J. Smith, "The failure of [R.sup.0]," Computational and Mathematical Methods in Medicine, vol. 2011, Article ID 527610, 17 pages, 2011.

[22] H. D. Gaff, D. M. Hartley, and N. P. Leahy, "An epidemiological model of rift valley fever," Electronic Journal of Differential Equations, vol. 115, pp. 1-12, 2007.

[23] B. Adams and M. Boots, "How important is vertical transmission in mosquitoes for the persistence of dengue? Insights from a mathematical model," Epidemics, vol. 2, no. 1, pp. 1-10, 2010.

[24] N. Chitnis, J. M. Hyman, and C. A. Manore, "Modelling vertical transmission in vector-borne diseases with applications to Rift Valley fever," Journal of Biological Dynamics, vol. 7, no. 1, pp. 1140, 2013.

[25] S. A. Pedro, S. Abelman, F. T. Ndjomatchoua, R. Sang, and H. E. Z. Tonnang, "Stability, bifurcation and chaos analysis of vectorborne disease model with application to rift valley fever," PLoS ONE, vol. 9, no. 10, Article ID e108172,2014.

[26] F. Chamchod, R. S. Cantrell, C. Cosner, A. N. Hassan, J. C. Beier, and S. Ruan, "A modeling approach to investigate epizootic outbreaks and enzootic maintenance of Rift Valley fever virus," Bulletin of Mathematical Biology, vol. 76, no. 8, pp. 2052-2072, 2014.

[27] S. A. Pedro, H. E. Tonnang, and S. Abelman, "Uncertainty and sensitivity analysis of a Rift Valley fever model," Applied Mathematics and Computation, vol. 279, pp. 170-186, 2016.

[28] S. A. Pedro, S. Abelman, and H. E. Z. Tonnang, "Predicting Rift Valley Fever Inter-epidemic Activities and Outbreak Patterns: Insights from a Stochastic Host-Vector Model," PLOS Neglected Tropical Diseases, vol. 10, no. 12, 2016.

[29] R. Ross, Murray London, 2nd edition, 1911.

[30] G. Macdonald, "The analysis of equilibrium in malaria," Tropical Diseases Bulletin, vol. 49, no. 9, pp. 813-829,1952.

[31] R. M. Anderson, Population Dynamics of Infectious Diseases: Theory and Applications, vol. 368,1982.

[32] N. A. Maidana and H. M. Yang, "Describing the geographic spread of dengue disease by traveling waves," Mathematical Biosciences, vol. 215, no. 1, pp. 64-77, 2008.

[33] H. R. Thieme, Mathematics in Population Biology, vol. 47, 2003.

[34] J. LaSalle, "The stability of dynamical systems," in Proceedings of the CBMS-NSF Regional Conference Series in Applied Mathematics, SIAM, Thailand, Philadelphia, USA, 1976.

[35] O. Diekmann and J. A. P. Heesterbeek, Mathematical Epidemiology of Infectious Diseases: Model Building, Analysis and Interpretation, John Wiley & Sons, 2000.

[36] P. van den Driessche and J. Watmough, "Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission," Mathematical Biosciences, vol. 180, pp. 29-48, 2002.

[37] M. J. Keeling and P. Rohani, Modeling Infectious Diseases in Humans and Animals, vol. 47, 2007

[38] M. G. Roberts and J. A. P. Heesterbeek, "A new method for estimating the effort required to control an infectious disease," Proceedings of the Royal Society B Biological Science,vol. 270, no. 1522, pp. 1359-1364, 2003.

[39] J. A. Heesterbeek and M. G. Roberts, "The type-reproduction number t in models for infectious disease control," Mathematical Biosciences, vol. 206, no. 1, pp. 3-10, 2007

[40] J. C. Kamgang and G. Sallet, "Global asymptotic stability for the disease free equilibrium for epidemiological models," Comptes Rendus Mathematique, vol. 341, no. 7, pp. 433-438, 2005.

[41] C. Castillo-Chavez and B. Song, "Dynamical models of tuberculosis and their applications," Mathematical Biosciences and Engineering, vol. 1, no. 2, pp. 361-404, 2004.

[42] J. M. Heffernan, R. J. Smith, and L. M. Wahl, "Perspectives on the basic reproductive ratio," Journal of the Royal Society Interface, vol. 2, no. 4, pp. 281-293, 2005.

[43] A. L. Lloyd, J. Zhang, and A. M. Root, "Stochasticity and heterogeneity in host-vector models," Journal of the Royal Society Interface, vol. 4, no. 16, pp. 851-863, 2007.

Sansao A. Pedro

Departamento de Matematica e Informatica, Universidade Eduardo Mondlane, 257, Maputo, Mozambique

Correspondence should be addressed to Sansao A. Pedro; sansaopedro@gmail.com

Received 7 July 2018; Revised 30 September 2018; Accepted 10 October 2018; Published 1 November 2018

Academic Editor: Zhidong Teng

Caption: Figure 1: Relationships among epidemic thresholds with respect to vertical transmission and ratio of mosquitoes to hosts. (a) The relationship between number of new infections in hosts and number of new infections in vectors. (b) Comparison between basic reproductive numbers [R.sup.0] and [R.sup.e.sub.0] derived from two different methods with respect to vertical infection. (c) Comparison of behaviours among the complete transmission cycle [R.sup.2.sub.0] and host type and vector type reproductive numbers with respect to vertical transmission efficiency away from the common epidemic criterion. (d) How the ratio of mosquitoes to hosts m and vertical infection [q.sub.V] contribute to patterns of [R.sup.e.sub.0]. Dashed lines indicated contours lines for [R.sup.e.sub.0] <1 while continuous lines are for [R.sup.e.sub.0] [greater than or equal to] 1.

Caption: Figure 2: The left figure (a) explores effects of vertical transmission on the prevalence of the disease on hosts, where ([I.sub.h]([q.sub.V])) is plotted at max([I.sub.h](t)) for every [q.sub.v]. The right figure (b) explores how the first peak of an epidemic (max([I.sub.h](t))) varies with vertical transmission efficiency ([q.sub.v]). The lower figures represent the dynamics of [I.sub.h](t) at time position (78) corresponding to day 27 (figure c) and at time position (150) corresponding to day 50 (figure d) for different values of vertical infection.

Caption: Figure 3: The bifurcation diagram in ([Q.sub.0], [R.sup.c.sub.0]) plane. The stability of equilibria in each region is given in Table 1.

Caption: Figure 4: Time series and phase diagrams plot. The inner figures describe the temporal evolution of [L.sub.i], [I.sub.v], [I.sub.h] after the transient period for [R.sup.c.sub.0] = 2.1002, while for the outer figures [R.sup.c.sub.0] = 1.0012. Parameter values used for [R.sup.c.sub.0] = 1.0012 are [[mu].sub.v] = 1/25, [[mu].sub.h] = 1/(8 * 360), [p.sub.hv] = 0.14, [p.sub.vh] = 0.35, [alpha] = 256/360, [[gamma].sub.h] = 1/4, m = 1.5, [q.sub.v] = 0.2, [theta] = 1/10, [[[mu].sub.L] = 1/18, r = 1/12, [delta] = 0.4 and parameter values used for [R.sup.c.sub.0] = 2.1002 are given in Table 1 of Section A-3 of Supplementary Materials.

Table 1: Existence and stability of model equilibria. [X.sup.0] [X.sup.1] if (0 < [g.sub.v] <1) I. [Q.sub.0] < 1 & GS DNE [R.sup.c.sub.0] < 1 II. [Q.sub.0] > 1 & DNE GS [R.sup.c.sub.0] < 1 III. [Q.sub.0] > 1 & DNE US [R.sup.c.sub.0] < 1 IV. [Q.sub.0] < 1 & US DNE [R.sup.c.sub.0] < 1 X*(if (0 < [g.sub.v] <1) I. DNE II. DNE III. LS IV. DNE DNE: does not exist, US: unstable, LS: locally stable, GS: globally stable.

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Sansao A. Pedro |

Publication: | Journal of Applied Mathematics |

Article Type: | Report |

Date: | Jan 1, 2018 |

Words: | 11458 |

Previous Article: | A Mathematical Model of Treatment and Vaccination Interventions of Pneumococcal Pneumonia Infection Dynamics. |

Next Article: | A Study of a Diseased Prey-Predator Model with Refuge in Prey and Harvesting from Predator. |

Topics: |