arxiv: v1 [q-bio.pe] 5 Aug 2015

Transcription

1 Mathematical Moeling of Dengue Epiemic: Control Methos an Vaccination Strategies Sylvestre Aureliano Carvalho, Stella Olívia a Silva,Iraziet a Cunha Charret arxiv: v1 [q-bio.pe] 5 Aug 2015 Abstract Dengue is one of the most important infectious iseases in the worl, in terms of eath an economic cost. Hence, the moeling of engue is of great importance to help us unerstan the ynamics isease, an interfering with its spreaing mathematical by the proposition of control methos. In this work, control strategies in an attempt to eliminate the Aees aegypti mosquito, as well as proposals for the vaccination campaign are evaluate. In our mathematical moel, the mechanical control is accomplishe through the environmental support capacity affecte by a iscrete function that represents removal of breeing. Chemical control is carrie out from the use of insecticie an larvicial. The efficacious of vaccination is stuie through the removal of a fraction of iniviuals, proportional to the rate of vaccination, from the susceptible compartment an its transfer to the recovere compartment. Our major fin is that the engue fever epiemic is only eraicate with the use of an immunizing vaccine because control measures irecte against the vector are not enough to halt isease spreaing. Even where the infecte mosquitoes are eliminate from the system, susceptible mosquitoes are still present, an infecte humans cause engue fever to reappear in the human population. Keywors: 1. Introuction In recent years the fiel of public health has been benefit from the use of mathematical moels to the stuy of sprea epiemic of infectious isease. One of the first mathematical moel epiemiology has been propose by W. H. Hamer in 1906, who postulate that the evelopment of an epiemic is relate to the rate of contact between infecte an susceptible iniviuals [1] an currently this postulate is known as the principle of mass action [2]. In 1911 Ross has propose one of the first moels use for epiemiological stuy of malaria [3]. In 1957 Maconal brought a greater biological realism to the Ross moel, improving the accuracy an interpretation of the moel parameters. The contributions of Maconal resulte in the creation of the moel known as the Ross-Maconal moel [4], in which two ifferential equations represent the ynamics of the isease regarless of the age structure of populations an the incubation perio. In the cases of infectious iseases, situations may exist in which infecte hosts once recovere again become susceptible to new infections or may become immune. However, there may be situations where the host contract the isease an o not acquire the isease elicite by the same agent again. In the first case, these iseases are usually cause by bacterial or protozoan agents [5, 6], such as malaria for example. The secon case is ue to viral iseases such as engue [7]. Currently, there is great interest in stuying epiemic propagation, as in the case of engue. This isease is consiere the 1 Departamento e Física, Universiae Feeral e e Viçosa, , Viçosa, MG, Brazil 2 Departamento e Ciências Exatas, Universiae Feeral e Lavras, , Lavras, MG, Brazil most significant in terms of eath an economic cost, since it is responsible for more eaths an infections in humans than any other arthropo-borne viruses [8]. This scenario can be even worse, since the exact number of cases of the isease is not known. Several cases go unreporte because they may clinically manifest with few or no symptoms [9]. Due to the large number of people expose, an absence of meicines able to cure infecte people, the increase of engue epiemics have motivate an intensive stuy of mechanical an chemical control mechanisms to use on eliminate the mosquito Aees aegypti. The mechanical control, relying on eucational campaigns for the removal of containers that accumulate water, results in a significant ecrease in population size of the vector when this type of control is performe on any of the seasons [10]. Moreover, the mechanical control is effective only upon the awareness of the human population to ientify an eliminate outbreaks of engue in their own homes [11]. Chemical control, base on the use of insecticie an larvicie, was teste for the reuction of vector population [11, 12]. More recently, in parallel, vaccination strategies have been propose as vaccines. Stuies have been one to evelop a vaccine, evaluate an investigate its efficacy [13, 14, 15]. So, in this paper, we evaluate the set of these three strategies - mechanical control, chemical control, an vaccination - in orer to assess the best way to eraicate engue. Our stuy is base on the engue epiemiological moel propose by Yang an Ferreira [16]. We will moel the effects by the use of mechanical control mechanisms (parameters that ecrease the rate of oviposition vector, known as carrying capacity) an chemical control (insecticie an larvicie). Besies the mortality rates of the vector we have also inclue in our stuy the spreaing of the engue virus in humans an suggest vaccination strategies Preprint submitte to Applie Mathematics an Computation August 6, 2015

2 in a human population. The paper is organize as follows. Section (2) introuces the orinary ifferential equations of the moel. In Section (3), mechanical an chemical control strategies, as well as the factors that justify the use of information base strategies in climate an rainfall moelling are escribe. It is also highlighte the effects of vaccination campaign an reporte a measure of the efficiency of these mechanisms. In section (4), the moel solutions obtaine by numerical integration are presente an the efficiency of each strategy is evaluate. Finally, in Section 5 some conclusions are rawn. 2. Mathematical moel for engue In the stuy of some epiemics, it is useful to separate the vector an host populations in some classes accoring to their stage in the infection cycle. In the case of engue, for instance, it is useful to consier the existence of three classes for the vector an four for the host. Such separation is suggeste, because both the vector an the host contract the virus, go through a perio of latency in which the iniviual has the virus but is not capable of transmitting it, became infective an, for humans, recovere from the isease [17]. The engue mosquito, the Aees aegypti female, can be infecte with the virus when feeing infecte human bloo. Once it have acquire the virus the mosquito, remains infecte throughout its life. So it is feasible to ivie the population of ault miges into susceptible, expose an infecte. Classes in turn, the human population can recover from the illness, an for that reason, in aition to the existing classes for the vector, there is the class of recovere humans [18]. The moel in which the iniviuals are organize in these ifferent classes from which they can entry or exit is known as SEIR(susceptible - expose - infectious - recovere) [19]. In this type of approach, a population is represente by state variables (subpopulation ensities) associate to the ifferent moel compartments, which envolves on time accoring the rules of interaction. Differential equations can be use to moel the rules, an give every subpopulation growth rates [20, 21]. Our mathematical moel, base on the moel evelope by Yang an Ferreira [16], comprises 10 nonlinear first orer orinary ifferential equations. These equations inclue variables escribing two classes of iniviuals: humans an mosquitoes. The phase of the mosquito in life cycle are represente in the system (1) as follows: egg (E), larva (L), pupa (P) an ault (W), the latter being compose of three subpopulations iffering in its stage of infection, namely, susceptible (W 1 ), expose (W 2 ) an infecte (W 3 ). The parameters of the moel are relate to biotic an abiotic factors. They represent averages of values that takes into account all stages of vector in evelopment, effective rates of natural mortality, an contact rates between the stages of infection of the human population an the mosquitoes. The transition rates between the phases of egg, larva, pupa an ault are represente by the parameters σ j, ( j = e, l, p, w). In aition, the system (1) contains other rates: oviposition 2 (φ), mortality of eggs, larvae, pupae an ault miges µ i (i = l, p, w); which are functions of temperature, humiity, an rain precipitation. Thus these rates reflect abiotic factors. However, in the present analysis, these functions will be replace by, constant average values, liste in table (1). t E(t) = [ φ 1 E(t) ] C W(t) (σ E + µ E ) E(t) (1) t L(t) = σ e E(t) (σ l + µ l + µ l ) L(t) t P(t) = σ l L(t) (σ P + µ P + µ P ) P(t) [ ] t W I(t) 1(t) = σ P P(t) β w N + µ w + µ w W 1 (t) t W 2(t) = I(t) β w N W 1(t) (γ w + µ w + µ w) W 2 (t) t W 3(t) = γ w W 2 (t) (µ w + µ w) W 3 (t) Other parameters correspon to biotic factors, i.e. those that have irect interference of the human activity: carrying capacity, represente by the letter C an the chemical control represente by the function µ i. The carrying capacity can be moulate through the elimination of breeing grouns as a result of eucational campaigns to increase the population awareness, concerning the necessities to eliminate the aquatic population reservoirs. The chemical control is another parameter which takes into account the work of health agents responsible for preventing vector proliferation in every stage of its life cycle, from vector aquatic phase to aulthoo. The population ault mosquitoes is compose only by females feeing bloo. The equations for the humans population are couple an separate in susceptible (s), expose (e), infecte (i) an recovere (r) compartiments. In our moel we take into consieration the presence of only one virus serotype. System (2) introuces a set of four first orer, nonlinear orinary ifferential equations for the human population ensities. [ β h M I (t) t s(t) = µ h t e(t) = β h t i(t) = γ h e(t) (σ h + µ h ) i(t) t r(t) = σ h i(t) µ h r(t) 3. Control strategies ] M(t) + µ h s(t) (2) M I (t) M(t) s(t) (γ h + µ h ) e(t) 3.1. The influence of climate an temperature The climate is one of the factors that affects the environment where the Aees aegypti is present [22]. The Aees is foun

3 primarily in tropical an subtropical areas in zones having relatively constant temperatures aroun 24 o C [23]. The temperature has irect an inirect influences on insect s evelopment an feeing [24]. Hence, the effects of temperature on the life cycle of insects rea be moele to unerstan the population ynamics of engue vectors an can be exploite to evelop appropriate vector control strategies. In Brazil, the fall season begins in March, winter in June, the spring starts in September, an the summer in December. However, these four seasons are neatly marke only in the southern region, São Paulo, Mato Grosso o Sul, an in the mountain regions of Minas Gerais an Rio e Janeiro, occupying approximately 15% of the brazilian territory. In such regions summer an winter temperatures are well-efine. In contrast, in the Amazonia, for example, oes not exist significant variation in temperature an rainfall uring the annual perio, so basically there are no seasons. In other regions, there are only two regimes: the rainy an the ry seasons. Taking into consieration the climatic variability in Brazil, we have ivie a year into three perios, efine in terms of the maximum an minimum rainfall a major eterminant of the mosquito population namely, favorable, unfavorable an intermeiate. Inee, the Aees aegypti has a weak resistant at to the low temperature. Temperatures between 22 o C an 30 o C are favorable in its life cycle, an the extremes, below 18 o C an above 34 o C have negative effects on their evelopment [25]. The most favorable perio for the evelopment of the vector are the months of December, January, February an March, ue to high temperatures an high rainfall. The coler months (June, July, August an September), are unfavourable for the evelopment of Aees aegypti. The intermeiate perio inclues the months of April, May, October an November [26]. The values that make up the template have aily values accoring to table (1) an were taken from stuies by Yang an Luz [16, 11] Mechanical Control It is wiely known that a population can not grow inefinitely because the environment has a limite capacity of resources to sustain their iniviuals. This ability is known as environmental supportability [27]. There population size of the engue mosquito increases consierably rainy season ue to the increase in the number of containers with staning water, ieal places for mosquito oviposition [28]. Thus, the environmental support capacity of the engue vector is primarily relate to the amount of available sites for laying eggs, since the resources necessary for the maintenance of aults vectors such, as foo an shelter, are virtually unlimite. In the moel, escribe by equations (1), the environmental support capacity is given by C for the population of eggs, an is interprete as the physical space available for mosquitoes oviposition. A metho to control the growth of mosquito population is to ecrease C through the elimination of breeing grouns. This is usually one mechanically, at ifferent times of the year, accoringly the climatic conitions. In orer to take ino account this control strategy, our moel assumes that the support capacity C is a function of a ranom 3 variable C i. The effectiveness reate of the control performe at given moment. Specifically, C = C i C f ixo (3) where C f ixo is the supportability of the containers in the absence of control mesures an 0 < C i 1 (4) varies accoring to the time of year. When C i receives values close to 0, leaing to a very small C, there are few locations for oviposition. In contrast, when C i 1, the support capacity remains virtually unchange, an there is no reuction of mosquito breeing sites. In turn, the value of C f ixo varies with the annual seasons. Within the favorable perio, characterize by high temperatures an humiity [29], a value C f ixo = 700 was fixe. Within (col an ari) the unfavorable perio [30] we set C f ixo = 300. For the months that comprise the mil perio, a supportability C f ixo = 500 is assume Chemical Control There are two types of chemical that hill the engue vector: insecticies, that affect ault mosquitoes, an larvicies, that kill vector larvae. We assume that the chemical control oes not act on the population of pupae (µ p = 0). The insecticies is the most use instruments for vector control. In turn, larvicies are complementarly imploye to this control. Hence, in our moel, these chemical methos are represente by the variable µ w an µ l that affect winge mosquitoes an larvae, respectively. These mortality rates are exponential ecay functions. They change with the perio of the year as shown in table (2). It is assume that the chemicals are release uring the first weeks of every month, an persists for 15 ays in the environment. Also, at the first ay the chemical release is the maximum, an at most 10% of the total amount of release chemicals remain in the environment Vaccination Campaign Currently the fight against engue is restricte to eraication methos of its vector. The evelopment of new tools to combat engue, as for example, the iscovery of a vaccine requires more knowlege about the biological characteristics of the virus an its interaction with the hosts. Our moel, escribe by equations 1 an 2, oes not take into account both people who coul ie by the isease an oes not escribe the existence of a vaccine that makes humans immune to engue virus. The evelopment of a tetravalent vaccine against the engue virus is a promise of the Worl Health Organization (WHO), but efforts to evelop vaccines have been hampere by the lack of an appropriate animal test subject [31]. In orer to take into account these two features of the epiemics, we extene our moel as follows. The eath of infecte human introuce through the µ = 0.04 aition of a mortality rate in the infecte compartment: t i(t) = γ h e(t) (σ h + µ h + µ ) i(t) (5)

4 Meaning of the Parameters Symbol Value (ay 1 ) Oviposition φ 1 Transformation from egg to larva(f) σ e Transformation from egg to larva (U) σ e Transformation from egg to larva (I) σ e Transformation from larva to pupa (F) σ l Transformation from larva to pupa (U) σ l Transformation from larva to pupa (I) σ l Transformation from pupa to mosquito (F) σ p Transformation from pupa to mosquito (U) σ p Transformation from pupa to mosquito (I) σ p Mortality rates of egg µ e Mortality rates of larva µ l Mortality rates of pupa µ p Mortality rates of mosquito (F) µ w Mortality rates of mosquito (U) µ w Mortality rates of mosquito (I) µ w Force of infection human infectious mosquito susceptible β w Force of infection human susceptible mosquito infectious β h Incubation perio of virus (mosquito) γ w Incubation perio of virus (in humans) γ h Infectious (or recovery) perio (in humans) σ h Rates per-capita mortality human µ h Rates per-capita natality human µ n Table 1: The parameters an their values use in the moel. The transformation rates in the aquatic phase as well as rate mosquitoes mortality are subivie into perios: Favorable (F), Unfavorable (U) an intermeiate (I). The values have units of ay 1. Perio Insecticie function Larvicie function Unfavorable µ w(t) = e t µ l (t) = e t Intermeiate µ w(t) = e t µ l (t) = e t Favorable µ w(t) = e t µ l (t) = e t Table 2: Continuous functions representing the control by insecticie an larvicial, at the favorable, unfavorable an intermeiate perios. To represent a possible vaccination campaign, a fraction proportional to a vaccination rate ν = 0.2 of susceptible iniviuals is remove from an ae irectly to the recovere compartment. Also, the efficiency of the vaccine is takes into account through a variable ξ. If ξ = 0, the vaccine is not effective an oes not change the system, whereas, if ξ = 1, the vaccine is fully effective, an all the vaccine iniviuals became immunize. So, the effect of a vaccination campaing on the susceptible an recovere populations are given by: [ ] t s(t) = µ M I (t) h β h M(t) + µ h s(t) ξ ν s(t) (6) t r(t) = σ h i(t) µ h r(t) + ν ξ s(t) (7) The uration of the campaign is varie to fin out the time neee to eliminate the virus using a vaccination campaign Efficiency For each control methos, its percentual efficiency Σ can be efine as Σ = ( 1 A c A s ) 100% (8) 4 where A C represents the area uner the affecte mosquito curve population as a function of time when the controlmetho is use an A S is the area uner the affecte total population curve without the use of the metho. 4. Results The results refer to the numerical simulation for a time perio of approximately 20 to 50 years. The parameter values use by the table (1) were aopte in orer to escribe coexistence between populations of humans an mosquitoes (aquatic an ault). Initially, we consier that the human population grows exponentially, as suggeste by the Unite Nations (UN). Finally, simulation for a stationary human population is consiere Human population with exponential growth Dengue is a isease recognize in more than 100 countries aroun the worl [32]. Accoring to the UN, over 70% of the worl s population live in unerevelope countries. Asia, with a population of 4.1 billion inhabitants, an Africa, with billion inhabitants exhibit, both a high population growth rate,

5 which can be characterize as exponential [33]. Since in those continents there is the presence of engue viruses, there is much current interest in analyse the engue epiemics in such exponentially growing populations Mechanical control The Aees aegypti female lays its eggs in places containing stagnant water which nurtures the evelopment of immature stages of the vector. Consequently, if breeing grouns are eliminate, a ecrease in the population size of the vector an, consequently, in cases of engue fever will expecte. In orer to test the impact of the mechanical control on the population of ault vectors, a perioical ecrease has been mae in the carrying capacity, following the rule escribe in section (3.2). Figure (1) compares the results when the mechanical control was accomplishe istinct over annual perios (favorable, unfavorable, intermeiary an throughout the year) (eshe line) with those in which no control was performe (full line). The greatest efficiency was obtaine when the mechanical control was carrie out throughout the year, as shown (1 ). Figure (2) escribes the evolution of the infecte mosquito an human populations, when no control is performe (full line) an when the mechanical control is performe only in the intermeiate perio (ashe line). In the time interval between 240 an 360 ays, there was an increase in the number of isease case when compare to simulations in which there were no mechanical control. However, in the next interval (between 360 an 480 ays) this number was smaller. This variation in the number of cases is ue to the interference of mechanical control. Reflexes are observe in the perio after the intervention, since this type of control eliminates the aquatic stages of the mosquito. Figure 1: Comparison between the population sizes of the engue vector with an without mechanical control. The control is hel only in (a) the most favorable perio, (b) unfavorable perio, (c) intermeiate perio an () throughout the year. The case when mechanical control is performe (eshe line) is compare to the case in which no control type was performe (full line). When the control is performe only in one of the annual perios, an increase in the number of humans infecte in perio, even with a reuction in the number of infecte vectors as 5 shown in figure 2. We can realize greater efficiencies through the use of mechanical control throughout the year. However, the Aees aegypti mosquito is not eraicate only using the mechanical control an, consequently, the virus remains in the human population. Hence, the use of insecticie an larvicial was teste in the following section in an attempt to exterminate the isease. The corresponing results are iscusse below. Figure 2: (a)population size of mosquitoes infecte with an without the mechanical control in the intermeiate perio. (b) Population of humans infecte when no control is performe an when the control is performe in the intermeiate perio Chemical Control Christophers foun that the mosquito in larval an ault stages are bille by of chemicals applie uring their evelopment [34]. We use the functions in the table (2), in orer to evaluate the impact of insecticies an larviciaes on the control of engue epiemic. In aition, to iscriminate the impact of ifferent types of chemical control, simulations were performe consiering only one type of control at a time. Figure (3) presents the population size of the Aees aegypti mosquito when an insecticie was applie. As can be seem, the use of the insecticie was not effective in eliminating the mosquito population, an when the application was mae throughout the year (Figure 3-), the impact on the total population was greater. However this control ha only 21% efficiency, as calculate accoring to the equation 8. In (4) are shown the results for the application of larvicial. The approximate efficiency is 5%, an the population practically evolves as if no control mechanism was use. This inicates that such control is not feasible in an attempt to eraicate the Aees aegypti mosquito. Figure (5) shows the effect of chemical control on the population of infecte humans. The result is inepenent on the chemical applications. Even with a ecrease in the number of iniviuals at a given time, this number increases again latter leaing o a control efficiency of less than 1%. This may inicate that the ensity of infecte mosquitoes oes not interfere

6 Figure 3: Mosquito population sizes after applying the insecticie only in (a) the most favorable, (b) unfavorable, (c) intermeiate perio, an () throughout the year. Figure 5: Comparison of the population size of humans infecte when no control is carrie out an when the chemical an mechanical control are performe together throughout the year. cally for agriculture Mechanical an chemical control Since the use of insecticie an larvicial separately was inefficient to eliminate the population of engue ault vectors, we consiere a scenario in which the two chemicals were applie simultaneously with the mechanical control throughout the year. The results are shown in fig. (6) for ifferent control strategies. When chemical control was combine with the mechanical control, the efficiency is 86% (blue curve). In comparison, if only the mechanical control (green curve) was use the efficiency fell to 26%. However, even if both types of control are use, it has not been possible to eliminate engue from the system, still leaving the human population at risk of contracting the isease. Figure 4: Comparison of the population size of mosquitoes after application of larvicial when the control is hel (a) only in the most favorable perio, (b) the unfavorable perio, (c) the intermeiate perio an () throughout the year. with the virus permanence into the system, but simply there is a critical population size that guarantees virus support. Accoringly, in countries in an enemic state, most of their istricts presente a vector inex below the value require for the risk of engue fever epiemics [29]. The incience of infection was high even in areas with low numbers of positive engue mosquito breeing sites [35]. Furthermore, small breeing are as with mosquito outbreaks in their neighborhoos can ensure the presence of engue virus in all other istricts of a city [37]. The long-term control one by chemical, such as those consiere here, shoul be analyse with caution because its continue use can generate resistant vectors populations [37]. In aition, accoring to the technical note n o 109/2010 [38], it is important to consier that there are no insecticies suitable for use in public health, since all these proucts are use specifi- 6 Figure 6: (Color online) Comparison between the ifferent applications of chemical an mechanical controls. The combine use of chemical mechanical controls resulte in an 86% efficiency (blue curve). In turn, the efficiency fell to 26% it only the mechanical control (green curve) is use Death by isease an vaccination campaign In previous simulations the fact that people can ie ue isease complications was not consiere. Dengue fever, especially in its most severe form, is hemorrhagic, an may have a lethality greater than 20% [32]. Hence, in orer to increase the moel s realism, we a to the human infecte compartment,

7 eq. (5), the rate of eath by illness µ. Since, the Epiemiological Bulletin of the Ministry of health in the year of 2005 reporte a 4% lethality level for engue, we aapt µ = 0.04 [? ]. This new factor has a quantitative impact on isease outbreaks, leaing to outbreaks elays at long-term, as can be seen in figure (7). Figure 8: Evolution in time of engue epiemic for ifferent urations of the vaccination campaign. Figure 7: Evolution in time of the human infecte population with an without the isease mortality rate parameter µ. Because, the chemical an mechanical controls were unable to eliminate the engue virus circulation in the human population, we teste the effects of an eventual vaccine against engue fever on its spreaing ynamics. In orer, to simulate a vaccination, a fraction, proportional to the vaccination rate, ν = 0.2, is withrawn from the compartment of susceptibles, equation (6), an irectly ae into the recovere compartment, equation (7). It is wiely known that vaccines license by the health agencies are not 100% effective. Due to this fact, the efficiency of the vaccine was also consiere an represente by a variable ξ. ξ = 0 correspons to an ineffective vaccine that oes not evelop immunity, whereas ξ = 1 means that the vaccine is fully effective, immunizing all the iniviuals receive the vaccine. With a vaccination rate fixe at 0.2, an an efficiency maintaine at ξ = 1 the time uring which the campaign shoul be maintaine in orer to eliminate the virus was estimate. Figure (8) shows that a campaign active for three consecutive months eraicate the virus from the population. Furthermore, figure (9) shows the evolution in time of the population infecte humans for ifferent values of ξ for a campaign lasting three months. A vaccine that is 60% efficient oes not eliminate the virus an with 70% there is a significant ecrease in the number of cases over time. In orer to eliminate the engue fever of the human population, at least 80% efficiency is require. Larger ξ values, such as 90%, only iffer in the time spent to eraicate the virus Human population without exponential growth As variously mentione, engue is present in many countries, incluing Brazil an other South American countries, where the rate of population growth is not enough to reach levels as 7 Figure 9: Comparison of ifferent rates of efficiency for the vaccine against the engue virus: (a) escribes the efficiency of 90%, (b) 80% efficiency, (c) 70% efficiency, an () 60% efficiency. high as that of the African countries. Accoring to a stuy by the Brazilian Institute of geography an statistics (IBGE 3 ), the Brazilian population is ecreasing its growth rate since the 60 s, when this rate was 3.04%. The IBGE estimates that this reuction in the growth rate still shoul last for many years reaching the so-calle zero growth in 2039 [40]. Hence, we also ajuste the birth rate of the human population to a value closer to the mortality rate. So, in these simulation the birth rate is now fixe in Consiering this new birth rate for the human population an the use of mechanical control throughout the year, an efficiency of 80% was reache. However, as shown in figure (10), the reuction of mosquitoes population oes not ecrease the infecte human population. The behavior remains the same when the control was not use. We note that even with a low number of infecte mosquitoes, the engue virus still remains in the human population. This is possible because, as proven through histological stuies, the engue mosquito fees on bloo several times uring a single cycle consequently, a single female 3 IBGE: Brazilian population ageing in rhythm accelerate

8 can contaminate several humans in the same cycle [41]. Figure 10: Comparison of mosquito (total an infecte) an infecte human infecte population when no control is in effect an when it is one using the mechanical control throughout the year with birth rate µ = (a) The total mosquitoes population, (b) the population of infecte mosquitoes, (c) population of infecte humans. Since only the mechanical control was enable to eliminate the engue fever epiemic, a chemical control was applie in orer to check its efficiency taking in to account this new birth rate. The system response to the use of insecticie an larvicial applie in ifferent perios of the year is shown in figure (11). chemical control, which ha a consierable impact to prevent the sprea of engue fever epiemic [12]. But the persistence of Aees aegypti in the system, even if it is not infecte, the engue fever may reappear at any time if any infecte person is introuce, an so the engue epiemic can re-emerge in the human population, since the vector of the virus is still present. It is known that the vector elimination campaigns by means of chemical agents often occurs in the favorable perio of the year. Accoringly, we perme the chemical control in this perio to see how the system respon to the emergence of infecte people. Figure (12) shows the result of simulations in which 10 iniviuals infecte with the engue virus appear at two ifferent perios of the year. If these infecte iniviuals come when chemical control is applie in the favorable perio (continuous line), engue oes not reappear an the system remains free of the isease. On the contrary, if they are introuce in the unfavorable perio, when there was no chemical control (ashe line), we have perioic outbreaks of the isease are observe. This suggests that a control campaign performe only in a perio of the year may not be effective in eliminating engue totally, since the virus can reappear through the emergence of infecte people. Figure 11: Evolution in time of mosquito an humans population with the use of chemical control an birth rate µ = for the favorable, unfavorable an intermeiate perios: (a) is the total population of mosquitoes, (b) is the population of infecte mosquitoes, (c) is population of infecte humans. Now we can notice that the chemical control is not efficient in eraicate the engue viruses regarless of the perio in which it is being applie, eliminating all infecte mosquito in the system, but this oes not reuce the levels significantly below the susceptible vector s population, as Luz [11]. In the work one by Burattini an et. al, it was possible to evaluate the use of 8 Figure 12: Evolution in time of populations of infecte (a) mosquitoes (b) humans for a scenary in with 10 infecte humans appear in a system free of viruses. The analyses were mae in the favorable an unfavorable perios. With the rate µ = the human population oes not grow exponentially, probably because the elimination of infecte persons through the use of chemical control becomes possible. But even with the complete elimination of the population of infecte vectors, it has not been possible to eliminate all the susceptible mosquito population in the system as shown in figure (12). The engue fever may reappear at any time, unless there is a vaccine to protect the population against the virus Vaccination campaign In this simulation the population oes not exhibit a high growth rate an the vaccination campaign will be represente

9 by the withrawal of a fixe number of susceptible iniviuals (ψ = 200) by aing them into the recovere population. Thus, the vaccination rate is no longer proportional to the number of susceptible iniviuals, an is escribe by equations (9) an (10). In this case, a eeme efficiency ξ = 80% for the vaccine is consiere since it is the immunization threshol acceptable by health agencies. s(t) t r(t) t =... ψ ξ (9) =... + ψ ξ (10) Figure (13) shows how the system behaves uner immunization of 200 people per ay for 1 month, with ifferent times of uration of the campaign. It can be observe that, if the campaign remains for only 5 consecutive years, engue virus re-emerge in populations of mosquitoes an humans an the epiemic reappears. However, if the uration of the campaign is 10 years, it is possible to eraicate the engue virus from the human population for at least 30 years, in a total simulations time fixe at 50 years. To fin out the time require for the permanent eraication of the virus from the system, more simulations consiering longer simulation times, as well as the possibility of the emergence of other virus serotypes must be carrie out. Figure 13: The vaccination campaign to immunize 200 iniviuals per ay with ifferent perios of uration. Several stuies show that a vaccine that protects the population from the engue virus is about to be create. An attenuate tetravalent vaccine virus has been teste in orer to characterize the interactions between the ifferent serotypes of engue virus, concluing that it is possible to create an effective vaccine against all four serotypes of the engue virus [42]. However, an important point to be consiere is whether live, attenuate virus vaccines o not offer risks to human health, since it is not clear that the presence of four types of strains into a single vaccine might increase the risk of eveloping severe forms of the isease [43] Conclusions Mathematical moels are quite useful for observing the ynamics of infectious iseases such as the engue fever. The creation of mathematical moels that simulate the interaction between host-vector generates momentum in the system. Through the compartmental moel propose by Yang an Ferreira was use, it was possible to introuce small changes in orer to test new control strategies. Testing the mechanical an chemical control verifie the vaccination strategies for the human population. The mechanical an chemical control were effective in reucing the population of Aees aegypti but assure that the isease has been eliminate from the system. The control that use only the larvicial ha a very low contribution in this reuction, probably because the time interval between the implementation an its impact on the population of aults is very high, making this control to be less appropriate in controlling an epiemic of engue fever. The mechanical control is much more efficient, which inicates that awareness campaigns are more useful than the use of controls through insecticie an larvicial. When consiering a human population birth rate equal to /ay, it was possible to simulate an exponential growth ynamics of this population. However, this consieration oes not affect the effectiveness of mechanical control. In this case the application of mechanical control must be constant throughout the year. Once this control is active only in time after his application, it generates an increase in the number of cases in the other two perios, as there is no control being carrie out. With respect to the use of chemical control, there was an interference mosquito population in a yiel of 21% for the insecticie an larvicie only 5%. When the chemical an mechanical control were use together, there was an efficiency of 26%. In the latter case, even with greater efficiency, the use of these controls i not influence the ynamics of engue in the human population, who remaine infecte. Importantly, the chemicals use in the simulations is fictitious an quite specific. Simulations with insecticie an larvicie with ifferent uration shoul be conucte so that we can guarantee these results. When the use of both controls was compare to situation where no control was accomplishe, it resulte in an efficiency of 86%, compare to the situation where only the mechanical control was carrie out in which the efficiency value fell to 26%, thus showing that the mechanical control is more efficient in the elimination of mosquito population than the chemical control even when applie throughout the year. When the birth rate was change to /ay, it was no longer possible to observe an exponential growth for the human population. In this case the efficiency of mechanical control remaine at 80%. In the same way for the previous simulations, there was no significant reflection on human population infecte. The chemical control can eliminate the infecte population but the outbreaks o not en with a susceptible population. With 10 iniviuals in the system infecte with the engue virus in the perio in which the control is not being applie in the unfavorable time, perioic outbreaks of the isease occur. The only way to contain the engue epiemic, regarless of

10 the form of population growth, is to create a vaccine to immunize the population against the virus. The population grew exponentially uring the immunization vaccination campaign for a fraction proportional to a vaccination rate of iniviuals, ν = 0.2. For this campaign to become really effective it is necessary to have a uration of three months an an efficiency rate of at least 80%. It was possible to note the existence of annual cycles an a strong seasonality of epiemic cycles [22]. A major concern for the safety of live attenuate virus vaccines is the genetic stability. The use of molecular tools is require to ensure quality, safety an authenticity of genetic vaccine [44]. Alreay when the exponential growth is remove, the campaign is base on aily immunization ψ = 200, of people, with an efficiency of 80%. In this case it is necessary to maintain the campaign for at least 10 years to ensure at least 30 years of a system free of the isease. Variable µ was introuce representing the portion of infecte who ie from complications relate to engue fever. The aition of this new factor cause the outbreaks of the isease might arise later, but as the simulation time is fixe at 50 years this variable can have a major stake, because it may elay the outbreak so that it will not reappear uring the time examine. However, fiel surveys are still necessary in orer to better unerstan the mechanisms responsible for generating these ynamics. In aition, more simulations must be carrie out in orer to consier a longer time of observation, while factoring in the possibility of the emergence of other serotypes of the isease an biological ways of eliminating the vector. Further stuies are neee to unerstan the mechanisms responsible for generating the ynamics of engue in orer to learn how to fight it. The nee for a vaccine is essential for the elimination of engue viruses in human populations, thus it is expecte that stuies about immunization are satisfactory an achieve a positive result as soon as possible. References [1] W. H. Hamer, Epiemic isease in Englan - the evience of variability an of persistency of type, Lancet, i, (1906), [2] E. B. Wilson, J. Woecester, The law of mass action in epiemiology, Proc. Natl. Aca. Sci., USA, 31 (1945) (part I) (part II). [3] N. Bacaër, The prevention of malaria - R. Ross, A short History of Mathematical Population Dynamics, Springer Lonon, Lonon, [4] G. Maconal, The Epiemiology an Control of Malaria, Lonon: Oxfor University Press, [5] M. Ghosh, p. Chanra, P. Sinha, J. B. Shukla, Moelling the sprea of bacterial isease: effect of service proviers from an environmentally egrae region, Appl. Math. Comput. 160 (2005) [6] S. Singh, J. B. Shukla, P. Chanra, Moelling an analysis of the sprea of malaria: Environmental an ecological effects, J. Bio. Syst. 13(1) (2005) 111. [7] N. A. Maiana, H. M. Yang, Spatial spreaing of West Nile Virus escribe by traveling waves, Jour. of Theor. Biol., 258 (2009) [8] D. J. Gubler, The global emergence/resurgence of arboviral isease as public health problems, Arc. Me. 33 (4) (2002), [9] I. kurane, T. Takasaki, Dengue fever an engue hemorrhagic fever: challenges of contolling an enemy still at large, Rev. Me. Virol., 11 (5) (2001) [10] M. B. Aruino, Assessment of Aees aegypti pupal prouctivity uring the engue vector control program in a costal urban centre of S o Paulo state, Brazil, Jour. of Insec., 2014 (2014) [11] P. M. Luz, C. T. Coeço, J. Melock, C. J. Struchiner, D. Valle, A. P. Galvani, Profile of insecticie interventions on the abunance an resistence profile ofaees aegypti, Epiemiol. Infect., 137 (8) (2009) [12] M. N. Burattini, M. Chen, A. Chow, F. A. Coutinho, K. T. Goh, L. F. Lopez, S. Ma, E. Massa, Moelling the control strategies against engue in Singapure, Epiemiol Infect., 136 (3) (2008) [13] A. J. johnson, J. T. Roehring, New mouse moel for engue virus vaccine testing, The Jour. of Virol., 73 (1) (1999) [14] M. A. Johansson, J. Hombach, D. A. Cummings, Moels of the impact of engue vaccines: a review of current research an potential approaches, Vaccine, 29 (35) (2011) [15] J. G. Julaner, S. T. Perry, S. Shresta, Important avances in the fiel of anti-engue virus research, Antivir. Chem. Chemother, 21 (3) (2011) [16] H. M. Yang, C. P. Ferreira, Assessing the effects of the vector control on engue transmission, Appl. Math. Comput., 198 (2008) [17] E. A. Newton, P. Reiter, A moel of the transmission of engue fever with an evolution of the impact of ultra-low volume (ULV) insecticie applications on engue epiemics. Am. J. Trop. Me. Hyg. 47 (6) (1992) [18] M. Chan, M. A. Johansson, The incubation perios of engue viruses, Plos One, 7 (11) (2012) pp. E [19] A. Korobeinicov, Global properties of SIR an SEIR epiemic moels with multiple parallel interactions stages, Bull. Math. Bio., 71 (2009) [20] L. Esteva, C. Vargas, Analysis of a engue isease transmission moel, Math. Biosci., 150 (1998) [21] M. Anrau, N. Hens, C. Marais, P. Beutels, Dynamic epiemiological moels for engue transmission: a systematic review of structural approaches, Plos One, 7 (11) (2012) E [22] K. Nakhapakorn, N. K. Tripathi, An information value-base analysis of physical an climatic factors affeting engue fever an engue hemorrhagic fever incience, Inter. Jour. of Hea. Geo., 4 (2005) [23] J. Setzer, Climate an Ecological Atlas of the State of São Paulo, Comissão Interestaual a Bacia o Paraná-Uruguai em Colaboração com as Centrais Elétricas e SP (CESP), São Paulo, Brazil, [24] D. J. Gubler, Dengue an engue hemorrhagic fever: its hystory an resurgence as a global health problem. In: D. J. Gubler, G. Kuno (E.) Dengue an engue hemorrhagic fever. New York: CAB International, [25] W. Köppen, Climatology con un stuio e los climas e la tierra. Mxico: Fono a Cultura Econmica, [26] E. B. Beserra et al., Biologia e exigências térmicas e Aees aegypti (L.) (Diptera: Culiciae) provenientes e quatro regiões bioclimáticas a Paraíba, Neotropical Entomology, 35 (6) (2006) [27] N. J. Gotelli, Ecologia. 4. e. Lonrina: Eitora Planta, [28] Worl Health Organization (WHO), Global strategy for engue prevention an control, 2012, [29] F. P. Câmara et al., Regional an ynamics characteristics of engue in Brazil: a retrospective stuy. Rev. a Soc. Bras. e Me. Trop. 40(2) (2007) [30] J. A. Koopman et al. Determinants an preictors of engue infection in Mexico. American Journal of Epiemiology, 133 (11) (1999) [31] Worl Health Organization (WHO), Dengue hemorrhagic fever: iagnosis, treatment, prevention an control, 2n e. Geneva [32], L. F. Bricks, Vacinas para a engue: perspectivas. Peiatria, 26 (4) (2004) [33] Population Reference Bureau, Worl Population Data Sheet:The Worl at 7 Billion, 2011, [34] S. R. Christophers, Aees aegypti (L.) the yellow fever mosquito. New York: Cambrigee University Press [35] M. G. Teixeira et al., Dynamics of engue virus circulation: a silent epiemic in a complex urban area, Tropical Meicine an International Health, 7 (9) (2002) [36] L. C. e C. Meeiros et a. Moeling the Dynamic Transmission of Dengue Fever: Investigating Disease Persistence. PLoS Neglecte Tropical Diseases, 5 (1) (2011) E942. [37] I. A. Braga, D. Valle, Aees aegypti: inseticias, mecanismos e ação e resistência. Epiemiologia e Serviços a Saúe, 16 (4) (2007) [38] Ministério a Saúe. O uso racional e inseticias no controle o Aees aegypti e sua utilização oportuna em áreas com transmissão e engue.

11 Nota Técnica N o 109/2010 CGPNCD/DEVEP/SVS/MS Brasília, [39] Ministério a Saúe. Guia e Vigilância Epiemiológica. 6n e. Brasília, [40] Instituto Brasileiro e Geografia e Estatstica IBGE. Daos sobre População o Brasil, PNAD (Pesquisa Nacional por Amostra e Domicílios), [41] A. C. Morrison, A. Getis, M. Santiago, Exploratory space-time analysis of reporte engue cases uring an outbreak in Floria, Puerto Rico, Amer. Jour. of Trop. Me. an Hyg., 57 (1997) [42] K. B. Anerson et al. Interference an facilitation between engue serotypes in a tetravalent live engue virus vaccine caniate. Journal of Infectious Diseases, 204 (3) (2011) [43] L. J. a Silva, R. Richtmann, Vaccines uner evelopment: group B streptococcus, herpes-zoster, HIV, malaria an engue. J. Peiatr. 82 (3) (2006) [44] H. S. Rorigues, M. T. Monteiro, D. F. Torres, Vaccination moels an optimal control strategies to engue, Math. Biosci., 247 (2014)