Open access peer-reviewed chapter

# Mathematical Model as a Tool for the Control of Vector-Borne Diseases: Wolbachia Example

Written By

Meksianis Z. Ndii, Eti D. Wiraningsih, Nursanti Anggriani and Asep K. Supriatna

Submitted: March 26th, 2018 Reviewed: June 26th, 2018 Published: November 5th, 2018

DOI: 10.5772/intechopen.79754

From the Edited Volume

## Dengue Fever

Edited by Jorge Abelardo Falcón-Lezama, Miguel Betancourt-Cravioto and Roberto Tapia-Conyer

Chapter metrics overview

View Full Metrics

## Abstract

Dengue is a vector-borne disease that risks two-thirds of the world’s population particularly in tropical and subtropical regions. Strategies have been implemented, but they are only effective in the short term. A new innovative and promising strategy against dengue is by the use of Wolbachia bacterium. This requires that Wolbachia-carrying mosquitoes should persist in the population. To assess the persistence of Wolbachia-carrying mosquitoes and its effects on dengue, a number of mathematical models have been formulated and analysed. In this chapter, we review the existing mathematical models of Wolbachia-carrying mosquito population dynamics and dengue with Wolbachia intervention and provide examples of the mathematical models. Simulations of the models are presented to illustrate the model’s solutions.

### Keywords

• Wolbachia
• dengue
• mathematical model

## 1. Introduction

Dengue is a vector-borne disease caused by four distinct serotypes (DEN1–DEN4), and is endemic in most countries particularly in tropical and subtropical areas [1]. It is estimated that around 390 million cases happen each year [2]. Individuals obtain lifelong immunity to the serotype that they are infected with, but have a higher chance to get the most severe form of dengue in the subsequent infection [1]. It is estimated that around 500,000 individuals get severe dengue and require hospitalisation. Of these, about 2.5% die [3]. Without a proper treatment, the fatality rate can reach 20% [3]. Dengue is also a substantial public health and economic burden [4].

A number of strategies have been implemented, but they are generally effective in the short term. Although some progresses have been made for dengue antiviral treatment, dengue control strategies still depend on vector control [5]. One of the strategies against dengue is by the use of Wolbachia bacterium. There are two Wolbachia strains used in the experiments: WMelPop and WMel. WMelPop strain can reduce the mosquito lifespan of more than 50% and almost 20% reduction in fecundity [6]. WMel strain reduces the lifespan of around 10% and only small reduction in fecundity [6]. Wolbachia can reduce the level of virus in the salivary glands. Wolbachia gives reproductive advantage for Wolbachia-carrying female mosquitoes known as cytoplasmic incompatibility (CI). The Wolbachia-carrying female mosquitoes can reproduce when mating with both non-Wolbachia and Wolbachia-carrying male mosquitoes. Non-Wolbachia female mosquitoes can reproduce when mating with non-Wolbachia males [7]. Field experiments showed that Wolbachia-carrying mosquitoes have established and dominated the population [8]. When Wolbachia-carrying mosquitoes persist in the field, the Wolbachia intervention can be implemented. The question that arises is that to what extend this intervention can reduce dengue transmission? To answer the above question, a number of mathematical models have been formulated and analysed. Mathematical model is a useful tool to understand complex phenomena. This can be used to understand population dynamics [9], disease transmission dynamics [10, 11], and others [12, 13]. A number of mathematical models have been developed to examine the persistence and spread of Wolbachia-carrying mosquitoes and its effects on dengue transmission dynamics. In this chapter, we review the existing mathematical models of Wolbachia-carrying mosquito population dynamics and dengue with Wolbachia intervention, give examples of the mathematical models, and show several numerical simulations to illustrate the model’s solutions.

## 2. Mathematical modelling

This section presents background on mathematical modelling of infectious diseases. Mathematical modelling is a useful tool to understand complex phenomena including disease transmission dynamics and their control strategies. There are several types of modelling that are generally used: deterministic, stochastic, statistical, agent-based modelling, and the others. A deterministic model is mostly used because it is easily solved and can include many parameters or variables. The model is in the form of system of differential equations.

Many mathematical models have been developed to investigate disease transmission dynamics including vector-borne diseases [14]. The model is based on a standard SIR model where the human population is divided into susceptible (S), infected (I), and recovered (R) [15]. The susceptible individuals become exposed after being contacted with infected individuals at a rate β . They then recover at a rate γ . The model is written in the following system of differential equations:

dS dt = β SI N , dI dt = β SI N γI , dR dt = γI . E1

The model can then be extended to include other compartments and parameters depending on the characteristics of diseases. For example, if the disease has long incubation period, we can add exposed compartment. If the disease is transmitted via vector, we can add another system of equations describing vector dynamics. When one aims to investigate the effects of vaccination, vaccinated compartment can be included. The important principles in modelling are to know characteristics of studied phenomena and the purpose of the research. The principles have been applied when we formulate mathematical models for Wolbachia-carrying mosquito population dynamics and dengue with Wolbachia.

## 3. Overview of mathematical models of Wolbachia and dengue

In this section, we review existing mathematical model of Wolbachia-carrying mosquito population dynamics and dengue with Wolbachia intervention.

Many (spatial and non-spatial) mathematical models have been formulated to analyse the persistence and spread or dispersal of Wolbachia-carrying mosquitoes in the populations [9, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29]. The general aim is to understand the underlying factors required for the persistence and spread of Wolbachia-carrying mosquitoes.

A number of nonspatial mathematical model for Wolbachia-carrying mosquito population dynamics have been developed. Ndii et al. [19] developed a mathematical model for Wolbachia-carrying mosquito population dynamics and assessed the persistence of Wolbachia-carrying mosquito populations. They found that Wolbachia-carrying mosquitoes persist in the population given that the death rate is not too high. Zhang et al. [30] formulated a mathematical model to assess the best strategies for releasing Wolbachia-carrying mosquitoes. They found that initial quantities of non-Wolbachia and Wolbachia-carrying mosquitoes and augmentation methods (timing, quantity, and order of frequency) determine the success of the Wolbachia intervention. They also formulated birth-pulse model with different density dependent death rate functions. They found that for condition with a strong density dependent death rate, the initial ratio of non-Wolbachia and Wolbachia-carrying mosquitoes should exceed a critical threshold for Wolbachia-carrying mosquitoes to dominate the population.

The spatial mathematical models have been developed to assess the Wolbachia-carrying mosquitoes’ dispersal. Chan and Kim [9] used reaction diffusion approach and incorporated slow and fast dispersal mode to assess the dynamics of the Wolbachia spread. They found that temperature affects the wavespeed of the Wolbachia-carrying Aedes aegypti, that is, Wolbachia invasion for Aedes aegypti increases when the temperature decreases within the optimal temperature rate for mosquito survival. Hancock et al. [17] developed a metapopulation model to assess the spatial dynamics of Wolbachia. They found that spatial variation in the density-dependent competition experienced by juvenile host insects can influence the spread of Wolbachia into population. In their other paper [16], they found a new expression for the threshold which takes into account the main aspects of insects’ life history. They showed that constant or pulsed immigrations affect the spread of Wolbachia-carrying mosquitoes.

Mathematical models for Wolbachia-carrying mosquitoes’ populations consider several important aspects. They are cytoplasmic incompatibility (CI), the maternal transmission, Wolbachia-carrying mosquito death rate, release strategies of Wolbachia-carrying mosquitoes [9, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 31, 32]. These are expressed in the parameters, variables, or simulations.

A number of mathematical models have been developed to understand dengue transmission mathematical models [33, 34]. However, little mathematical models have been developed to investigate the efficacy of Wolbachia-intervention [35, 36, 37, 38, 39, 40] in reducing dengue transmission. Hancock et al. [39] developed a mathematical model and investigated the strategies for releasing Wolbachia-carrying mosquitoes and its effects on dengue transmission dynamics. They found that male-biased releases can substantially reduce the dengue transmission. Furthermore, male-biased release can be an effective strategy that results in the persistence of Wolbachia-carrying mosquitoes. Ndii et al. [36, 41, 42] formulated single and two serotype dengue mathematical models to investigate the Wolbachia effectiveness in reducing dengue transmission. They found that Wolbachia can reduce primary and secondary dengue infections with higher reduction in secondary infections. Hughes and Britton [35] found that Wolbachia can reduce dengue transmission in areas where the basic reproduction number is not too high. This implies that Wolbachia can reduce dengue transmission in areas with low to moderate transmission settings, which is similar to that found Ndii et al. [36] and Ferguson et al. [37]. Supriatna et al. [40] showed that Wolbachia can reduce the value of basic reproduction number. In their other paper, they showed that the predatory and Wolbachia can reduce primary and secondary infections [43]. Furthermore, Supriatna et al. [44] investigated the use of vaccine and Wolbachia on dengue transmission dynamics [44] and showed that the optimal dengue control is determined by the epidemiological parameters and economic factors. Furthermore, they found that introducing too many Wolbachia-carrying mosquitoes would be counter-productive.

## 4. Examples and numerical simulations of mathematical models

In this section, we present examples of mathematical models of Wolbachia-carrying mosquito population dynamics and dengue with Wolbachia intervention and their numerical simulations.

### 4.1. Mathematical model of Wolbachia-carrying mosquito population dynamics and numerical simulations

#### 4.1.1. Mathematical model of Wolbachia-carrying mosquito population dynamics

Here, we present an example of the mathematical model of the Wolbachia-carrying mosquito population dynamics. We present the model by Ndii et al. [19, 45] and show several numerical simulations. The mosquito population is divided into aquatic ( A N and A W ), male ( M N and M W ) and female ( F N and F W ) mosquitoes. Note that the aquatic compartment consists of eggs, larvae and pupae, which are grouped into one compartment. Furthermore, the subscripts M and W denote non-Wolbachia and Wolbachia-carrying mosquito population.

The effect of CI is captured by the following expression. The non-Wolbachia female mosquitoes reproduce when mating with non-Wolbachia males, which is governed by the following equations:

ρ N M N F N M N + F N + M W + F W E2
and the Wolbachia-carrying females reproduce when mating with non-Wolbachia and Wolbachia-carrying males, which is governed by the following equations:
ρ W F W M N + M W M N + F N + M W + F W E3

Note that the population growth is limited by carrying capacity K . The maternal transmission is not perfect [46]. This means that not all Wolbachia-carrying aquatic mature to be Wolbachia-carrying adult. There is a proportion of 1 α that mature to be non-Wolbachia adults that is ϵ NW 1 α . Note that the ratio of male and female mosquitoes is denoted by ϵ ( ϵ N , ϵ W , ϵ NW ). The model is governed by the following systems of differential equations:

dA N dt = ρ N M N F N P 1 A N + A W K μ NA A N γ N A N , dM N dt = ϵ N γ N A N μ N M N + ϵ NW 1 α W γ W A W , dF N dt = 1 ϵ N γ N A N μ N F N + 1 ϵ NW 1 α W γ W A W , dA W dt = ρ W F W M W + M N P 1 A N + A W K μ WA A W γ W A W , dM W dt = ϵ W α W γ W A W μ W M W , dF W dt = 1 ϵ W α W γ W A W μ W F W . E4

where P = M N + F N + M W + F W is the total population.

#### 4.1.2. Numerical simulations

In this section, numerical simulations are conducted to illustrate the solutions of the model. The parameter values used are given in Table 1. The initial conditions are A N 0 = 0 , F N 0 = M N 0 = 7253 , A W 0 = 0 , and M W 0 = F W 0 = 14200 .

Symbol Description Value Unit Source
ρ N Non-Wolbachia reproductive rate 1.25 day 1 [19]
μ NA Non-Wolbachia aquatic death rate 1/7.78 day 1 [47]
γ N Non-Wolbachia maturation rate 1/6.67 day 1 [48]
ϵ N The proportion of non-Wolbachia adult male offspring 0.5 Proportion [49]
μ N Non-Wolbachia adult death rate 1/14 day 1 [47]
μ WA Wolbachia aquatic death rate 1/7.78 day 1 [47]
μ W Wolbachia adult death rate 1/7 day 1 [46]
ρ W Wolbachia reproductive rate 0.95 ρ N day 1 [19]
γ W Wolbachia maturation rate 1/6.67 day 1 [46]
ϵ W The proportion of Wolbachia-infected male adults 0.5 N/A Assumed
ϵ NW The rate of uninfected males hatched from a Wolbachia-infected mother 0.5 N/A Assumed
α W The proportion of Wolbachia-infected offspring from a Wolbachia-infected mother 0.9 N/A [7, 46, 50, 51]
K Carrying capacity 300,000 [48]

### Table 1.

Parameters, description, values and sources for the model of Wolbachia-carrying mosquitoes.

Symbol Description Value Unit Source
α Maternal transmission 0.9 N/A [19, 46, 52]
B Human birth rate 1 / 70 × 365 day 1 [3]
b N Biting rate of non-W mosquitoes 0.63 day 1 [53]
b W Biting rate of W mosquitoes 0.95 b N day 1 [54]
γ H Progression rate from exposed to infectious 1/5.5 day 1 [1]
γ N Progression from exposed to infectious class of non-W mosquitoes 1/10 day 1 [55]
γ W Progression from exposed to infectious class of W mosquitoes 1/10 day 1 [55]
K Carrying capacity 3 × N H N/A [55]
λ Force of infection Eqs. (21)(23)
μ N Adult mosquito death rate (non-W) 1/14 day 1 [47]
μ H Human death rate 1/ 70 × 365 day 1 [3]
μ NA Death rate of aquatic non-W mosquitoes 1/14 day 1 [47]
μ W Adult aquatic death rate 1.1 μ N day 1 [46, 51]
μ WA Death rate of W mosquitoes 1/14 day 1 [47]
ϕ ADE 1.1 N/A [56]
ρ N Reproductive rate of non-W mosquitoes 1.25 day 1 [19]
ρ W Reproductive rate of W-mosquitoes 0.95 ρ N day 1 [46]
σ Recovery rate 1/5 day 1 [1]
T N Transmission probability from non-W mosquitoes to human 0.5 N/A [36]
T HW Transmission probability from W mosquitoes to human 0.5 T N N/A [36, 57]
θ Progression rate from temporary immunity class to susceptible class 1 / 0.5 × 365 day 1 [58]
τ N Maturation rate of non-W mosquitoes 1/10 day 1 [47]
τ W Maturation rate of W mosquitoes 1/10 day 1 [47]

### Table 2.

Parameter descriptions, values, and sources. Note that W and N are used to indicate Wolbachia-carrying and non-Wolbachia mosquitoes in the parameter descriptions, respectively. N H = 10 5 .

Figure 1 shows the numerical solutions of the model using the parameter values given in Table 1, but the Wolbachia adult mosquito death rate is 2 × μ N . This reflects the WMelPop Wolbachia strain which reduces the mosquito lifespan by a half. Figure 1 shows that the non-Wolbachia mosquitoes dominate the population. This means that this strain cannot be used as a strategy to reduce dengue transmission. Figure 2 shows the numerical solutions of the model using the parameter values given in Table 1. The Wolbachia mosquito death rate is 1.1 × μ N which reflects the WMel Wolbachia strain. This strain reduces the mosquito lifespan by around 10%. It shows that the Wolbachia-carrying mosquitoes dominate the population. This means that WMel strain can be used in the Wolbachia intervention. Figure 3 shows the simulation results using WMel parameter values with initial conditions of A N 0 = 0 , F N 0 = M N 0 = 7253 , A W 0 = 0 , and M W 0 = F W 0 = 145 . It shows that the non-Wolbachia mosquitoes dominate the populations. It implies that the initial conditions also affects the persistence of Wolbachia-carrying mosquitoes.

### 4.2. Dengue mathematical model and numerical simulations

#### 4.2.1. Dengue mathematical model in the presence of Wolbachia

In this section, we give example of two-serotype dengue mathematical model. We present the model by Ndii et al. [42]. The model consists of human, non-Wolbachia and Wolbachia-carrying mosquito population. The human population is divided into susceptible ( S H ), exposed to serotype i ( E H i ), infected to serotype i ( I H i ), temporary immunity to the serotype i ( X H i ), recovered class ( R H ), susceptible, exposed and infected to j strain ( S H ji , E H ji , I H ji , respectively). The superscript ji means individuals that were previously infected by serotype i and currently infected by serotype j . The mosquito population is divided into susceptible ( S N and S W ), exposed to serotype i ( E N i and E W i ) and infected to serotype i ( I N i and I W i ). The subscript N and W is for non-Wolbachia and Wolbachia-carrying mosquitoes.

The model is governed by the following system of differential equations:

dS H dt = BN H i = 1 2 λ H i S H μ H S H , E5
dE H i dt = λ H i S H γ H E H i μ H E H i , E6
dI H i dt = γ H E H i σ I H i μ H I H i , E7
dX H i dt = σ I H i θ i X H i μ H X H i , E8
dS H ji dt = λ H j S H ji μ H S H ji , E9
dE H ji dt = λ H j S H ji γ H E H ji μ H E H ji , E10
dI H ji dt = γ H E H ji σ I H ji μ H + d I H ji , E11
dR H dt = j = 1 , j i 2 σ I H ji μ H R H , E12

Model for non-Wolbachia mosquito population

dA N dt = ρ N F N 2 2 F N + F W 1 A N + A W K τ N A N μ NA A N , E13
dS N dt = τ N A N 2 + 1 α τ W A W 2 i = 1 2 λ N i S N μ N S N , E14
dE N i dt = λ N i S N γ N E N i μ N E N i , E15
dI N i dt = γ N E N i μ N I N i . E16

Model for Wolbachia-carrying mosquito population

dA W dt = ρ W F W 2 1 A N + A W K τ W A W μ WA A W , E17
dS W dt = ατ W A W 2 i = 1 2 λ W i S W μ W S W , E18
dE W i dt = λ W i S W γ W E W i μ W E W i , E19
dI W dt = γ W E W i μ W I W i , E20

where the force of infections are

λ H i = b N T i I N i N H + b W T HW i I W i N H , E21
λ N i = b N T i I H i N H + ϕ i b N T i I H ij N H , E22
λ W i = b W T i I H i N H + ϕ i b W T i I H ij N H , E23

where ϕ i is the antibody-dependent enhancement factor for serotype i . Note that the susceptible human becomes exposed to dengue after being bitten by non-Wolbachia and Wolbachia-infected mosquitoes, which then becomes infected and have temporary immunity. After a certain period in temporary immunity class, they become susceptible to the other dengue serotype. They will have secondary infection after being bitten by infected mosquitoes carrying different dengue serotype to that they are previously infected.

#### 4.2.2. Numerical simulations

This section presents numerical simulations of the model. Figures 4 and 5 show the numerical simulations of primary and secondary infections in the absence and presence of Wolbachia, respectively.

Figures 4 and 5 show that Wolbachia can reduce dengue transmission. The number of infections in the presence of Wolbachia-carrying mosquitoes (see Figure 5) is smaller than that in the absence of Wolbachia-carrying mosquitoes (see Figure 4). This means that the Wolbachia can potentially be used to break the cycle of dengue transmission. Note that the parameter values are largely uncertain. Therefore, large data set is needed to validate the model against data.

## 5. Discussion and conclusions

The use of Wolbachia bacterium has been proposed as a new innovative strategy against dengue. A lot of research have been conducted to look at the persistence of Wolbachia-carrying mosquitoes and the potential reduction in the number of dengue cases by the use of Wolbachia bacterium. One of the approaches is by the use of mathematical model. It can be seen that mathematical model can provide insights into the persistence and the effectiveness of the Wolbachia in reducing dengue transmission dynamics.

One of the important steps in modelling is model’s validation. The model can be validated against the real data. Although several parameters can be obtained from literature, it is important to estimate the influential parameters such as transmission rate against the real data. Ferguson et al. [37] validated their model against the real data. Furthermore, most parameters are strongly uncertain, which indicate that sensitivity analysis is strongly required. This aims to find the most important parameters which can guide us in collecting appropriate data to be estimated.

Models presented in this work do not take into account the environmental factors such as temperature and rainfall. These may affect the dynamics of mosquito population and hence dengue transmission dynamics. Furthermore, in our work, the ratio of male and female mosquitoes is equal, which possibly affects the mosquito’s population dynamics. It is important to consider sex-biased ratio to determine its effects on the persistence of Wolbachia-carrying mosquitoes and dengue reduction.

In this paper, we review existing mathematical models of Wolbachia-carrying mosquitoes’ population dynamics and dengue with Wolbachia. Examples of the mathematical models are given. It shows that Wolbachia-carrying mosquitoes can persist in the population depending on the Wolbachia strains. Furthermore, the initial conditions also affect the persistence of Wolbachia-carrying mosquito populations. It is shown that Wolbachia can potentially reduce the primary and secondary infections with higher reduction in secondary infections. Results suggest that using Wolbachia can potentially reduce the transmission of dengue and hence minimise the public health and economic burden.

The results showed that the Wolbachia can persist in the population. When mosquitoes are infected with the WMel strain of Wolbachia. For dengue mathematical models with Wolbachia, it shows that the Wolbachia can potentially reduce the primary and secondary infections. This means that using Wolbachia can be an alternative strategy against dengue.

## Acknowledgments

MZN acknowledges funding from Ministry of Research, Technology and Higher Education of Indonesia through Penelitian Pascadoctor scheme 2018 (Grant No: 70/UN15.19/LT/2018). EDW received funding from Ministry of Research, Technology and Higher Education of Indonesia (Grant No: 3/E/KPT/2018). NA received Academic Leadership Grant from Padjadjaran University (Grant No: 2297 /UN6.D/KS/2018) and Penelitian Dasar Unggulan Perguruan Tinggi (Grant No: 1126/UN6.D/LT/2018).

## References

1. 1. Gubler DJ. Dengue and dengue hemorrhagic fever. Clinical Microbiology Reviews. 1998;11:480-496
2. 2. Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, et al. The global distribution and burden of dengue. Nature. 2013;496:1476-4687
3. 3. World Health Organisation. Dengue and Severe Dengue, Fact Sheet No. 117 [Online]. February 2018. Available from: http://www.who.int/news-room/fact-sheets/detail/dengue-and-severe-dengue [Accessed: June 1, 2018]
4. 4. Shepard DS, Undurraga EA, Halasa YA, Stanaway JD. The global economic burden of dengue: A systematic analysis. The Lancet Infectious Diseases. 2016;16(8):935-941
5. 5. Dorigatti I, McCormack C, Nedjati-Gilani G, Ferguson NM. Using Wolbachia for dengue control: Insights from modelling. Trends in Parasitology. 2018;34(2):102-113
6. 6. McMeniman CJ, Lane RV, Cass BN, Fong AWC, Sidhu M, Wang Y-F, O’Neill SL. Stable introduction of a life-shortening Wolbachia infection into the mosquito Aedes aegypti. Science. 2009;323(5910):141-144
7. 7. Sinkins SP. Wolbachia and cytoplasmic incompatibility in mosquitoes. Insect Biochemistry and Molecular Biology. 2004;34(7):723-729
8. 8. Hoffmann AA, Iturbe-Ormaetxe I, Callahan AG, Phillips BL, Billington K, Axford JK, Montgomery B, Turley AP, O’Neill SL. Stability of the WMel Wolbachia infection following invasion into Aedes aegypti populations. PLoS Neglected Tropical Diseases. 2014;8:e3115
9. 9. Chan MHT, Kim PS. Modelling a Wolbachia invasion using a slow–fast dispersal reaction–diffusion approach. Bulletin of Mathematical Biology. 2013;75(9):1501-1523
10. 10. Ndii MZ, Supriatna AK. Stochastic mathematical models in epidemiology. Information. 2017;20:6185-6196
11. 11. Ndii MZ, Amarti Z, Wiraningsih ED, Supriatna AK. Rabies epidemic model with uncertainty in parameters: Crisp and fuzzy approaches. IOP Conference Series: Materials Science and Engineering. 2018;332(1):012031
12. 12. Ndii MZ, Carnia E, Supriatna AK. Mathematical models for the spread of rumors: A review. In: Gaol FL, Hutagalung F, Chew FP, editors. Issues and Trends in Interdisciplinary Behavior and Social Science. USA: CRC Press; 2018. pp. 266-290
13. 13. Tambaru D, Djahi BS, Ndii MZ. The effects of hard water consumption on kidney function: Insights from mathematical modelling. AIP Conference Proceedings. 2018;1937(1):020020
14. 14. Reiner RC, Perkins T, Barker CM, et al. A systematic review of mathematical models of mosquito-borne pathogen transmission: 1970–2010. Journal of the Royal Society Interface. 2013;10(81):20120921
15. 15. Keeling MJ, Rohani P. Modeling Infectious Diseases in Humans and Animals. USA: Princeton University Press; 2008
16. 16. Hancock PA, Sinkins SP, Godfray HCJ. Population dynamic models of the spread of Wolbachia. The American Naturalist. 2011;177(3):323-333
17. 17. Hancock PA, Godfray HCJ. Modelling the spread of Wolbachia in spatially heterogeneous environments. Journal of Royal Society Interface. 2012;9(76):3045-3054
18. 18. Hancock PA, White VL, Callahan AG, Godfray CHJ, Hoffmann AA, Ritchie SA. Density dependent population dynamics in Aedes aegypti slow the spread of wMel Wolbachia. Journal of Applied Ecology. 2016;53(3):785-793
19. 19. Ndii MZ, Hickson RI, Mercer GN. Modelling the introduction of Wolbachia into Aedes aegypti to reduce dengue transmission. The ANZIAM Journal. 2012;53:213-227
20. 20. Li Y, Liu X. An impulsive model for Wolbachia infection control of mosquito-borne diseases with general birth and death rate functions. Nonlinear Analysis: Real World Applications. 2017;37:412-432
21. 21. Turelli M, Barton NH. Deploying dengue-suppressing Wolbachia: Robust models predict slow but effective spatial spread in Aedes aegypti. Theoretical Population Biology. 2017;115:45-60
22. 22. Xue L, Manore CA, Thongsripong P, Hyman JM. Two-sex mosquito model for the persistence of Wolbachia. Journal of Biological Dynamics. 2017;11(Suppl. 1):216-237
23. 23. Zhang X, Tang S, Liu Q, Cheke RA, Zhu H. Models to assess the effects of non-identical sex ratio augmentations of Wolbachia-carrying mosquitoes on the control of dengue disease. Mathematical Biosciences. 2018;299:58-72
24. 24. Guevara-Souza M, Vallejo EE. A computer simulation model of Wolbachia invasion for disease vector population modification. BMC Bioinformatics. 2015;16(1):317
25. 25. Strau JF, Telschow A. Modeling the indirect effect of Wolbachia on the infection dynamics of horizontally transmitted viruses. Frontiers in Microbiology. 2015;6:378
26. 26. Bliman PA, Aronna MS, Coelho FC, Moacyr A, da Silva HB. Ensuring successful introduction of Wolbachia in natural populations of Aedes aegypti by means of feedback control. Journal of Mathematical Biology. 2018;76(5):1269-1300
27. 27. Li Y, Liu X. A sex-structured model with birth pulse and release strategy for the spread of Wolbachia in mosquito population. Journal of Theoretical Biology. 2018;448:53-65
28. 28. Huang M, Luo J, Hu L, Zheng B, Yu J. Assessing the efficiency of Wolbachia driven aedes mosquito suppression by delay differential equations. Journal of Theoretical Biology. 2018;440:1-11
29. 29. Zheng B, Guo W, Hu L, Huang M, Yu J. Complex Wolbachia infection dynamics in mosquitoes with imperfect maternal transmission. Mathematical Biosciences and Engineering. 2018;15:523-541
30. 30. Zhang X, Tang S, Cheke RA. Models to assess how best to replace dengue virus vectors with Wolbachia-infected mosquito populations. Mathematical Biosciences. 2015;269:164-177
31. 31. Keeling MJ, Jiggins FM, Read JM. The invasion and coexistence of competing Wolbachia strains. Heredity. 2003;91:382-388
32. 32. Rasgon JL. Using Predictive Models to Optimize Wolbachia-Based Strategies for Vector-Borne Disease Control. New York, NY: Springer New York; 2008. pp. 114-125
33. 33. Supriatna AK, Soewono E, van Gils SA. A two-age-classes dengue transmission model. Mathematical Biosciences. 2008;216(1):114-121
34. 34. Anggriani N, Supriatna AK, Soewono E. A critical protection level derived from dengue infection mathematical model considering asymptomatic and symptomatic classes. Journal of Physics: Conference Series. 2013;423(1):012056
35. 35. Hughes H, Britton NF. Modelling the use of Wolbachia to control dengue fever transmission. Bulletin of Mathematical Biology. 2013;75(5):796-818
36. 36. Ndii MZ, Hickson RI, Allingham D, Mercer GN. Modelling the transmission dynamics of dengue in the presence of Wolbachia. Mathematical Biosciences. 2015;262:157-166
37. 37. Ferguson NM, Kien DTH, Clapham H, Aguas R, Trung VT, Chau TNB, Popovici J, Ryan PA, O’Neill SL, McGraw EA, Long VT, Dui LT, Nguyen HL, Van Vinh Chau N, Wills B, Simmons CP. Modeling the impact on virus transmission of Wolbachia-mediated blocking of dengue virus infection of Aedes aegypti. Science Translational Medicine. 2015;7(279):279ra37
38. 38. Qu Z, Xue L, Hyman JM. Modeling the transmission of Wolbachia in mosquitoes for controlling mosquito-borne diseases. SIAM Journal on Applied Mathematics. 2018;78(2):826-852
39. 39. Hancock PA, Sinkins SP, Godfray HCJ. Strategies for introducing Wolbachia to reduce transmission of mosquito-borne diseases. PLoS Neglected Tropical Diseases. 2011;5(4):1-10
40. 40. Supriatna AK, Anggriani N. System dynamics model of Wolbachia infection in dengue transmission. Procedia Engineering. 2012;50:12-18
41. 41. Ndii MZ, Allingham D, Hickson RI, Glass K. The effect of Wolbachia on dengue outbreaks when dengue is repeatedly introduced. Theoretical Population Biology. 2016;111:9-15
42. 42. Ndii MZ, Allingham D, Hickson RI, Glass K. The effect of Wolbachia on dengue dynamics in the presence of two serotypes of dengue: Symmetric and asymmetric epidemiological characteristics. Epidemiology and Infection. 2016;144(13):2874-2882
43. 43. Anggriani N, Supriatna AK. The effect of Wolbachia introduction and predatory in dengue disease transmission. Advanced Science, Engineering and Medicine. 2015;7(10):864-868
44. 44. Supriatna AK, Anggriani N, Nurulputri L, Wulantini R, Aldila D. The optimal release strategy of Wolbachia infected mosquitoes to control dengue disease. Advanced Science, Engineering and Medicine. 2014;6(8):831-837
45. 45. Ndii MZ. Analysis of a Wolbachia intervention to reduce dengue spread [Master’s thesis]. The Australian National University; 2011
46. 46. Walker T, Johnson PH, Moreira LA, Iturbe-Ormaetxe I, Frentiu FD, et al. The WMel Wolbachia strain blocks dengue and invades caged Aedes aegypti populations. Nature. 2011;476:450-453
47. 47. Yang HM, Macoris MLG, Galvani KC, Andrighetti MTM, Wanderley DMV. Assessing the effects of temperature on the population of Aedes aegypti, the vector of dengue. Epidemiology and Infection. 2009;137(8):1188-1202
48. 48. Maidana NA, Yang HM. Describing the geographic spread of dengue disease by traveling waves. Mathematical Biosciences. 2008;215(1):64-77
49. 49. Ruang-areerate T, Kittayapong P. Wolbachia transinfection in aedes aegypti: A potential gene driver of dengue vectors. Proceedings of the National Academy of Sciences. 2006;103(33):12534-12539
50. 50. Xi Z, Khoo CCH, Dobson SL. Wolbachia establishment and invasion in an Aedes aegypti laboratory population. Science. 2005;310:326-328
51. 51. Yeap HL, Mee P, Walker T, Weeks AR, O’Neill SL, et al. Dynamics of the “Popcorn” Wolbachia infection in outbred Aedes aegypti informs prospects for mosquito vector control. Genetics. 2011;187:583-595
52. 52. Hoffmann AA, Turelli M, Harshman LG. Factors affecting the distribution of cytoplasmic incompatibility in Drosophila simulans. Genetics. 1990;126:933-948
53. 53. Scott TW, Amerasinghe PH, Morrison AC, Lorenz LH, Clark GG, Strickman D, Kittayapong P, Edman JD. Longitudinal studies of Aedes aegypti (Diptera: Culicidae) in Thailand and Puerto Rico: Blood feeding frequency. Journal of Medical Entomology. 2000;37:89-101
54. 54. Turley AP, Moreira LA, O’Neill SL, McGraw EA. Wolbachia infection reduces blood–feeding success in the dengue fever mosquito, Aedes aegypti. PLoS Neglected Tropical Diseases. 2009;3:e516
55. 55. Chowel G, Diaz-Duenas P, Miller JC, Velazco AA, Hyman JM, Fenimore PW, Castillo-Chaves C. Estimation of the reproduction number of dengue fever from spatial epidemic data. Mathematical Biosciences. 2007;208:571-589
56. 56. Hu K, Thoens C, Bianco S, Edlund S, Davis M, Douglas J, Kaufman JH. The effect of antibody-dependent enhancement, cross immunity, and vector population on the dynamics of dengue fever. Journal of Theoretical Biology. 2013;319:62-74
57. 57. Bian G, Xu Y, Lu P, Xie Y, Xi Z. The endosymbiotic bacterium Wolbachia induces resistance to dengue virus in Aedes aegypti. PLoS Pathogen. 2010;6:e1000833
58. 58. Wearing HJ, Rohani P. Ecological and immunological determinants of dengue epidemics. Proceedings of the National Academy of Sciences. 2006;103:11802-11807

Written By

Meksianis Z. Ndii, Eti D. Wiraningsih, Nursanti Anggriani and Asep K. Supriatna

Submitted: March 26th, 2018 Reviewed: June 26th, 2018 Published: November 5th, 2018