A model of dengue fever
 M Derouich^{1},
 A Boutayeb^{1, 2}Email author and
 EH Twizell^{2}
DOI: 10.1186/1475925X24
© Derouich et al; licensee BioMed Central Ltd. 2003
Received: 29 November 2002
Accepted: 19 February 2003
Published: 19 February 2003
Abstract
Background
Dengue is a disease which is now endemic in more than 100 countries of Africa, America, Asia and the Western Pacific. It is transmitted to the man by mosquitoes (Aedes) and exists in two forms: Dengue Fever and Dengue Haemorrhagic Fever. The disease can be contracted by one of the four different viruses. Moreover, immunity is acquired only to the serotype contracted and a contact with a second serotype becomes more dangerous.
Methods
The present paper deals with a succession of two epidemics caused by two different viruses. The dynamics of the disease is studied by a compartmental model involving ordinary differential equations for the human and the mosquito populations.
Results
Stability of the equilibrium points is given and a simulation is carried out with different values of the parameters. The epidemic dynamics is discussed and illustration is given by figures for different values of the parameters.
Conclusion
The proposed model allows for better understanding of the disease dynamics. Environment and vaccination strategies are discussed especially in the case of the succession of two epidemics with two different viruses.
Background
With medical research achievements in terms of vaccination, antibiotics and improvement of life conditions from the second half of the 20^{th} century, it was expected that infectious diseases were going to disappear. Consequently, in developed countries the efforts have been concentrated on illnesses as cancer. However, at the dawn of the new century, infectious diseases are still causing suffering and mortality in developing countries. Malaria, yellow fever, AIDS, Ebola and other names will have marked the memory of humanity forever.
Among these diseases, dengue fever, especially known in Southeast Asia, is sweeping the world, hitting countries with tropical and warm climates. It is transmitted to the man by the mosquito of the genus Aedes and exists in two forms: the Dengue Fever (DF) or classic dengue and the Dengue Haemorrhagic Fever (DHF) which may evolve toward a severe form known as Dengue Shock Syndrome (DSS). The major problem with dengue is the fact that the disease is caused by four distinct serotypes known as DEN1, DEN2, DEN3 and DEN4. A person infected by one of the four serotypes will never be infected again by the same serotype (homologus immunity), but he looses immunity to the three other serotypes (heterologus immunity) in about 12 weeks and then becomes more susceptible to developing dengue haemorrhagic fever.
The first form (DF) is characterized by a sudden fever without respiratory symptoms, accompanied by intense headaches making its nickname "breakbone fever" well deserved. It lasts between three and seven days but it stays benign. The haemorrhagic form (DHF) is also characterized by a sudden fever, nausea, vomiting and fainting due to low blood pressure caused by fluid leakage. It usually lasts between two and three days and can lead to death. The case of a second infection has therefore a capital importance because of the possibility of evolution toward the haemorrhagic form of the disease.
So far, the strategies of mosquito control by insecticides or similar techniques proved to be inefficient. Moreover, the deterioration of the environment, the climatic changes, the unsanitary habitat, the poverty and the uncontrolled urbanization are as many favorable factors to the infectious illness propagation in general and dengue fever in particular. During the last decades the global prevalence of the dengue progressed dramatically. The disease is now endemic in more than 100 countries of Africa and America. The Southeast of Asia and the Western Pacific are seriously affected by the illness.
Before 1970, only nine countries had known epidemics of dengue haemorrhagic fever, but this number had increased more than fourfold in 1995 and about 2500 million of people are now exposed to the risk of the dengue fever. According to the present evaluations of the World Health Organization (WHO), about 50 million cases of dengue occur in the world every year, with an increasing tendency. In 1998, there were more than 616, 000 cases of dengue in America, of which 11, 000 cases of dengue haemorrhagic fever, that's twice the number of cases recorded in the same region during the year 1995. In 2001 there were 400, 000 cases of haemorrhagic fever in Southeast Asia, whereas, in Rio de Janeiro alone, 500, 000 people were struck by a dengue outbreak in 2002. The epidemic effect reached Florida and southern Texas [1–3].
During epidemics of the dengue fever the rate of the infectious among the susceptibles is often between 40 and 50 % but it may reach 80–90 % in favorable geographic and environmental conditions. Every year more then 500, 000 cases of dengue haemorrhagic require an hospitalization.
By contrast to malaria, caused by a parasite mainly in rural areas and transmitted by mosquito bites only at night, dengue is a viral disease resulting from the interaction of susceptible individuals with any of the four serotypes, mosquitoes of genus Aedes. The two recognized species of the vector transmetting dengue are Aedes aegypti and aedes albopictus. The first is highly anthropophilic, thriving in crowded cities and biting primarily during the day while the later is less anthropophilic and inhabits rural areas. Consequently, the importance of dengue is twofold: (i) Even in the absence of fatal forms, and because of its wide spreading and its multiple serotypes, the disease breeds significant economic and social costs (absenteeism, immobilization debilitation, medication). (ii) The potential risk of evolution towards the haemorrhagic form and the dengue shock syndrome with high economic costs and which may lead to death.
Mathematical modelling became an interesting tool for the understanding of these illnesses and for the proposition of strategies. The formulation of the model and the possibility of a simulation with parameter estimation, allow tests for sensitivity and comparison of conjunctures [4].
In the case of dengue fever, the mathematical models we have found in the literature propose compartmental dynamics with Susceptible, Exposed, Infective and Removed (immunised). In particular, SEIRS models [5] and SIR models [6] with only one virus or two viruses acting simultaneously [7] were considered.
In the present paper, while pointing out that the idea of two viruses coexisting in the same epidemic is controversial [7], a model with two different viruses acting at separated intervals of time is proposed. Our main purpose is to study the dynamics of dengue fever, while concentrating on its progression to the haemorrhagic form, in order to understand the epidemic phenomenon and to suggest strategies for the control of the disease. In search of clarity and simplicity, we assume that the latent period is not crucial for the susceptibleinfective interaction, hence we omit the compartment of exposed and consider a SIR model.
Formulation of the model and stability analysis
Parameters of the model
We suppose that we dispose of a human population (respectively of mosquito population) of size N _{ h }(resp. N _{ v }) formed of Susceptibles S _{ h }, of Infective I_{ h }and of Removed R _{ h }(resp. S _{ v }and I _{ v }).
The model supposes a homogeneous mixing of human and mosquito population so that each bite has an equal probability of being taken from any particular human. While noting b _{ s }the average biting rate of susceptible vectors, p _{ hv }the average transmission probability of an infectious human to a susceptible vector, the rate of exposure for vectors is given by: (p _{ hv } I _{ h } b _{ s } )/N _{ h }. It is admitted [6] that some infections increase the number of bites by the infected mosquitos in relation to the susceptible, therefore, we will assume that the rate of infected mosquito bites b _{ i }is greater than the one of the susceptible mosquitos b _{ s }.
Noting p _{ vh }the average transmission probability of an infectious vector to human and I _{ v }the infectious vector number, the rate of exposure for humans is given by: (p _{ vh } I _{ v } b _{ i })/N _{ h }so:

The adequate contact rate of human to vectors is given by: C _{ hv }= P _{ hv } b _{ s }

The adequate contact rate of vectors to human is given by: C _{ vh }= p _{ vh } b _{ i }.
The man life span is taken equal to 25 000 days (68.5 years), and the one of the vector is of 4 days. Other parameters values are given in the following section.
Basis parameters
definitions and values of basis parameters used in simulations [5]
Name of the parameter  Notation  Base value 

transmission probability of vector to human  p _{ hv }  0.75 
transmission probability of human to vector  p _{ vh }  0.75 
Bites per susceptible mosquito per day  b _{ s }  0.5 
Bites per infectious mosquito per day  b _{ i }  1.0 
Effective contact rate, human to vector  C _{ hv }  0.375 
Effective contact rate, vector to human  C _{ vh }  0.75 
Human life span  25000 days  
Vector life span  4 days  
Host infection duration  3 days 
Equations of the model
Up to now there is no vaccine against dengue viruses but research is going on and the eventuality of an immunization program is not excluded in the medium term. In this study we investigate the effect of such an immunization option and we also discuss the possibility of a partial vaccination against each serotype that will enable the control of the second epidemic and the evolution of dengue to dengue haemorrhagic fever.
In the case of first epidemic, the simplest assumption is that a random fraction, p, of susceptible humans can permanently be immunized against all the four serotypes. While for the second epidemic, a partial immunization is applied to the removed from the first epidemic. The dynamics of this disease in the host and vector populations is given by the following equations:
First epidemic
The model is governed by the following equations:
Human population
Vector population
With the conditions S _{ h }+ I _{ h }+ R _{ h }= N _{ h }and S _{ v }+ I _{ v }= N _{ v }, so:
R _{ h }= N _{ h } S _{ h } I _{ h }and S _{ v }= N _{ v } I _{ v }
then the two previous systems become:
Equilibrium points
Let the set Ω given by:
Ω = {(S _{ h },I _{ h },I _{ v })/0 ≤ I _{ v }≤ N _{ v }; 0 ≤ I _{ h }; 0 ≤ S _{ h }, (1 + p/μ_{ h })S _{ h }+ I _{ h }≤ N _{ h }}
then we have the following theorem:
Theorem 1
The previous system admits two equilibrium points E _{1}(N _{ h }/(1 + p/μ_{ h }), 0, 0) and E _{2} ( , , ) where , and
Proof:
see appendix A.1
Remark 1

if R  1 ≤ p' the trivial state E _{1} is the only equilibrium.

if R  1 >p' then the endemic equilibrium E _{2} will also be in Ω.
stability
Theorem 2
i) For R ≤ 1 + p' the state E _{1}(N _{ h }/(1 + p'), 0, 0) is globally asymptotically stable (ie ).
ii) For R > 1 + p' the state E _{2} ( , , ) is locally asymptotically stable.
Proof
see appendix A.2
Remark 2
The inequality R  1 ≤ p' (ie μ_{ h }(R  1) ≤ p) represents the principle of herd immunity because the susceptible population may be protected from epidemics if enough people are immunized.
Second epidemic
Human population
Vector population
Equilibrium points
Let the set Ω given by:
Theorem 3
The previous system admits two equilibrium points: E _{1}( /(1 + p'), 0, 0) and E _{2}( , , ) with , , and
Remark 1
We have two equilibrium points: the first E _{1} = ( /(1 + p'), 0, 0) is trivial, corresponding to to the state where the whole population is and will remain healthy. The second point is: E _{2} = ( , , ) it corresponds to the endemic state in which the disease persists in the two populations.
stability
Theorem 4
i) For R _{1} ≤ 1 + p' the state E _{1}( , 0, 0) is globally asymptotically stable
ii) For R _{1} > 1 + p' the state E _{2} ( , , ) is locally asymptotically stable.
Results and Discussion
Assuming a vaccination program as an option that would enable a proportion p of susceptible humans to be globally immunized against the four serotypes, the stability analysis shows that the population may be protected from epidemics if p satisfies the principle of herd immunity: μ_{ h }(R  1) ≤ p, where R is a function of the model parameters, defined as the number of new infected by an infective in interaction with susceptibles. Otherwise, there will be an endemic equilibrium and the disease will persist. But the problem is precisely in finding a vaccine against the four serotypes, which makes a strategy based on global immunization irrealistic in the short term. Meanwhile, a search leading to a partial vaccine against each serotype should be more feasible. In this direction, the proposed model shows that the evolution of dengue to dengue haemorrhagic fever can be controlled by a partial vaccine restricted to the people affected by the first epidemic.
In order to illustrate the dynamics of each epidemic and to study different strategies, a simulation was carried out using MATLAB routines with different values of the parameters implied in each model.
Mainly two directions can be envisaged to control the disease. The first may act on the number of mosquitoes and the second may consider the number of susceptible humans.
Conclusion
As mentioned in the introduction, Our main purpose was to study the dynamics of dengue fever and its progression to the dengue haemorrhagic fever in order to understand the epidemic phenomenon and to suggest strategies for the control of the disease in general and the haemorrhagic form in particular. The nature of dengue epidemics is complex since it conjugates human, environmental, biological, geographical and socioeconomic factors. Our model and simulations show that the strategy based on the prevention of dengue epidemic using vector control through environmental management (to eliminate the larval resting places such as containers like bottles, limp of canned foods, tire or other objects susceptible to keep water) or chemical methods (application of insecticides) remains insufficient since it only permits to delay the outbreak of the epidemic (figure 4). This conclusion agrees with the experiences realized by EA. Newton and P. Reiter using insecticides [5].
On the other hand, although the model suggests the reduction of susceptibles via vaccination, such a strategy is unlikely to be applicable in the short term because it faces some hurdles due to the fact that a vaccine must protect against the four serotypes at the same time. However, we consider this option since its eventuality is not excluded in the medium and long term.
In the short term, an intermediate solution would be to combine as much as possible, the environmental prevention and a partial vaccination essentially to avoid the haemorrhagic form of the disease caused by different viruses. This suggestion may help healthcare policy makers to tackle environment causes as preventive measures and researchers to investigate and concentrate on the search for a vaccine against each serotype rather than looking for a vaccine against the four serotypes at the same time.
Appendix
proof of theorem 1
the equilibrium points satisfy the following relations:
From the equation (3) we have: C _{ hv } I _{ h }/N _{ v }(N _{ v } I _{ v })  μ_{ v } I _{ v }= 0
From the equation (1) we have: μ_{ h } N _{ h } (μ_{ h }+ p + C _{ vh } I _{ v }/N _{ h })S _{ h }= 0 then
From the equation (2) we have:
therefore it can be seen that values of I _{ h }are: = 0 or substituting the equilibrium values of in the equations (4) and (5) we obtain two equilibrium points:
i) the first E _{1} = (N _{ h }/(1 + p'), 0, 0) is trivial in the sense that all individual are healthy and stay healthy for all time.
ii) The second point is: E _{2} = ( , , ), (where , and ) that corresponds to the endemic state i.e the case where the disease persists in the two populations. ■
proof of theorem 2
i) For E _{1} the matrix of linearization (Jacobian matrix) is giving by:
so E _{1} is stable for R < 1 + p'.
For the global stability, we consider the following Liapunov function:
thus:
So in Ω and for R ≤ 1 + p' we have: ≤ 0.
⇒ If R < 1 + p' then (N _{ h } (1 + p')S _{ h })I _{ v }= 0, I _{ h }= 0.
and If R = 1 + p' then (N _{ h } (1 + p')S _{ h })I _{ v }= 0, I _{ v } I _{ h }= 0.
Thus the set {E _{1}} is the largest invariant set within the set {(x, y, z)/ (x, y, z) = 0}. So according to the invariant set theorem: every trajectory in Ω tends to E _{1} as time t increases and as E _{1} is locally stable then it is globally asymptotically stable
ii) the point E _{2}
The local stability of E _{2} is governed by the matrix of linearization (Jacobian matrix) of E _{2} is given by:
P(λ) = λ^{3} + A λ^{2} + B λ + C where:
Or we have:
So: AB >C then following RouthHurwitz conditions for the polynomial P, the state E _{2} is locally asymptotically stable for R > 1 + p' ■
Declarations
Acknowledgement
This work was supported by the Royal Society of London under the DWSV scheme and the Moroccan C.N.R under the program PARS MI 23.
Authors’ Affiliations
References
 Derouich M: Modélisation et Simulation de modèles avec et sans structure d'âge : Application au diabète et à la fièvre dengue. Ph.D thesis Faculty of Sciences, Oujda Morocco 2001.Google Scholar
 Sesser S: Plague proportion. The Asian Wall Street Journal 2002 August 30 Septembre 1Google Scholar
 Dengue and Dengue Haemorrhagic Fever [http://www.who.int/inf.fs/en/fact117.html]
 Hethcote HW: The Mathematics of Infectious Diseases. SIAM review 2000, 42: 599–653.MathSciNetView ArticleGoogle Scholar
 Newton EA, Reiter P: A model of the transmission of dengue fever with an evaluation of the impact of ultralow volume (ULV) Insecticide applications on dengue epidemics. Am J Trop Med Hyg 1992, 47: 709–720.Google Scholar
 Esteva L, Vargas C: Analysis of a dengue disease transmission model. Mathematical Biosciences 1998, 150: 131–151. 10.1016/S00255564(98)100032View ArticleGoogle Scholar
 Feng Z, Hernàndez V: Competitive exclusion in a vectorhost model for the dengue fever. Journal of Mathematical Biology 1997, 35: 523–544. 10.1007/s002850050064MathSciNetView ArticleGoogle Scholar
 Luenberger DG: Introduction to Dynamic systems, Models, and Applications. Theory, Models, and Applications, John Wiley & sons 1997.Google Scholar
 Luenberger DG: Introduction to Applied Non linear systems and Chaos. S Wiggins, Springer 1996.Google Scholar
 Dietz K: Transmission and control of arbovirus diseases. In Proceedings of the Society for Industrial and Applied Mathematics, Epidemiology: Philadelphia (Edited by: Ludwing D). Philadelphia 1974, 104–106.Google Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.