Open Access is an initiative that aims to make scientific research freely available to all. To date our community has made over 100 million downloads. It’s based on principles of collaboration, unobstructed discovery, and, most importantly, scientific progression. As PhD students, we found it difficult to access the research we needed, so we decided to create a new Open Access publisher that levels the playing field for scientists across the world. How? By making research easy to access, and puts the academic needs of the researchers before the business interests of publishers.
We are a community of more than 103,000 authors and editors from 3,291 institutions spanning 160 countries, including Nobel Prize winners and some of the world’s most-cited researchers. Publishing on IntechOpen allows authors to earn citations and find new collaborators, meaning more people see your work not only from your own field of study, but from other related fields too.
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.
chapter and author info
Meksianis Z. Ndii*
Department of Mathematics, Faculty of Sciences and Engineering, University of Nusa Cendana, Indonesia
Eti D. Wiraningsih
Department of Mathematics, Faculty of Mathematics and Natural Sciences, State University of Jakarta, Indonesia
Department of Mathematics, Faculty of Mathematics and Natural Sciences, Padjadjaran University, Indonesia
Asep K. Supriatna
Department of Mathematics, Faculty of Mathematics and Natural Sciences, Padjadjaran University, Indonesia
*Address all correspondence to: email@example.com
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 . It is estimated that around 390 million cases happen each year . 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 . It is estimated that around 500,000 individuals get severe dengue and require hospitalisation. Of these, about 2.5% die . Without a proper treatment, the fatality rate can reach 20% . Dengue is also a substantial public health and economic burden .
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 . One of the strategies against dengue is by the use of Wolbachiabacterium. There are two Wolbachiastrains used in the experiments: WMelPopand WMel. WMelPopstrain can reduce the mosquito lifespan of more than 50% and almost 20% reduction in fecundity . WMelstrain reduces the lifespan of around 10% and only small reduction in fecundity . Wolbachiacan reduce the level of virus in the salivary glands. Wolbachiagives reproductive advantage for Wolbachia-carrying female mosquitoes known as cytoplasmic incompatibility (CI). The Wolbachia-carrying female mosquitoes can reproduce when mating with both non-Wolbachiaand Wolbachia-carrying male mosquitoes. Non-Wolbachiafemale mosquitoes can reproduce when mating with non-Wolbachiamales . Field experiments showed that Wolbachia-carrying mosquitoes have established and dominated the population . When Wolbachia-carrying mosquitoes persist in the field, the Wolbachiaintervention 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 , 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 Wolbachiaintervention, give examples of the mathematical models, and show several numerical simulations to illustrate the model’s solutions.
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 . The model is based on a standard SIR model where the human population is divided into susceptible (S), infected (I), and recovered (R) . 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:
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 Wolbachiaand dengue
In this section, we review existing mathematical model of Wolbachia-carrying mosquito population dynamics and dengue with Wolbachiaintervention.
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.  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.  formulated a mathematical model to assess the best strategies for releasing Wolbachia-carrying mosquitoes. They found that initial quantities of non-Wolbachiaand Wolbachia-carrying mosquitoes and augmentation methods (timing, quantity, and order of frequency) determine the success of the Wolbachiaintervention. 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-Wolbachiaand 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  used reaction diffusion approach and incorporated slow and fast dispersal mode to assess the dynamics of the Wolbachiaspread. They found that temperature affects the wavespeed of the Wolbachia-carrying Aedes aegypti, that is, Wolbachiainvasion for Aedes aegyptiincreases when the temperature decreases within the optimal temperature rate for mosquito survival. Hancock et al.  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 Wolbachiainto population. In their other paper , 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.  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 Wolbachiaeffectiveness in reducing dengue transmission. They found that Wolbachiacan reduce primary and secondary dengue infections with higher reduction in secondary infections. Hughes and Britton  found that Wolbachiacan reduce dengue transmission in areas where the basic reproduction number is not too high. This implies that Wolbachiacan reduce dengue transmission in areas with low to moderate transmission settings, which is similar to that found Ndii et al.  and Ferguson et al. . Supriatna et al.  showed that Wolbachiacan reduce the value of basic reproduction number. In their other paper, they showed that the predatory and Wolbachiacan reduce primary and secondary infections . Furthermore, Supriatna et al.  investigated the use of vaccine and Wolbachiaon dengue transmission dynamics  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 Wolbachiaintervention 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 (and ), male (and ) and female (and ) mosquitoes. Note that the aquatic compartment consists of eggs, larvae and pupae, which are grouped into one compartment. Furthermore, the subscripts and denote non-Wolbachiaand Wolbachia-carrying mosquito population.
The effect of CI is captured by the following expression. The non-Wolbachiafemale mosquitoes reproduce when mating with non-Wolbachiamales, which is governed by the following equations:
and the Wolbachia-carrying females reproduce when mating with non-Wolbachiaand Wolbachia-carrying males, which is governed by the following equations:
Note that the population growth is limited by carrying capacity . The maternal transmission is not perfect . This means that not all Wolbachia-carrying aquatic mature to be Wolbachia-carrying adult. There is a proportion of that mature to be non-Wolbachiaadults that is . Note that the ratio of male and female mosquitoes is denoted by (). The model is governed by the following systems of differential equations:
where 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 , , , and .
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. .
Figure 1 shows the numerical solutions of the model using the parameter values given in Table 1, but the Wolbachiaadult mosquito death rate is . This reflects the WMelPop Wolbachiastrain which reduces the mosquito lifespan by a half. Figure 1 shows that the non-Wolbachiamosquitoes 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 Wolbachiamosquito death rate is which reflects the WMel Wolbachiastrain. This strain reduces the mosquito lifespan by around 10%. It shows that the Wolbachia-carrying mosquitoes dominate the population. This means that WMelstrain can be used in the Wolbachiaintervention. Figure 3 shows the simulation results using WMelparameter values with initial conditions of , , , and . It shows that the non-Wolbachiamosquitoes 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. . The model consists of human, non-Wolbachiaand Wolbachia-carrying mosquito population. The human population is divided into susceptible (), exposed to serotype (), infected to serotype (), temporary immunity to the serotype (), recovered class (), susceptible, exposed and infected to strain (, respectively). The superscript means individuals that were previously infected by serotype and currently infected by serotype . The mosquito population is divided into susceptible (and ), exposed to serotype (and ) and infected to serotype (and ). The subscript and is for non-Wolbachiaand Wolbachia-carrying mosquitoes.
The model is governed by the following system of differential equations:
Model for non-Wolbachiamosquito population
Model for Wolbachia-carrying mosquito population
where the force of infections are
where is the antibody-dependent enhancement factor for serotype . Note that the susceptible human becomes exposed to dengue after being bitten by non-Wolbachiaand 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 Wolbachiacan 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 Wolbachiacan 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.
The use of Wolbachiabacterium 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 Wolbachiabacterium. 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 Wolbachiain 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.  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 Wolbachiastrains. Furthermore, the initial conditions also affect the persistence of Wolbachia-carrying mosquito populations. It is shown that Wolbachiacan potentially reduce the primary and secondary infections with higher reduction in secondary infections. Results suggest that using Wolbachiacan potentially reduce the transmission of dengue and hence minimise the public health and economic burden.
The results showed that the Wolbachiacan persist in the population. When mosquitoes are infected with the WMelstrain of Wolbachia. For dengue mathematical models with Wolbachia, it shows that the Wolbachiacan potentially reduce the primary and secondary infections. This means that using Wolbachiacan be an alternative strategy against dengue.
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).
Meksianis Z. Ndii, Eti D. Wiraningsih, Nursanti Anggriani and Asep K. Supriatna (November 5th 2018). Mathematical Model as a Tool for the Control of Vector-Borne Diseases: Wolbachia Example, Dengue Fever - a Resilient Threat in the Face of Innovation, Jorge Abelardo Falcón-Lezama, Miguel Betancourt-Cravioto and Roberto Tapia-Conyer, IntechOpen, DOI: 10.5772/intechopen.79754. Available from:
Over 21,000 IntechOpen readers like this topic
Help us write another book on this subject and reach those readers
We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.