- Research
- Open Access
- Published:

# Assessing potential countermeasures against the dengue epidemic in non-tropical urban cities

*Theoretical Biology and Medical Modelling*
**volume 13**, Article number: 12 (2016)

## Abstract

### Background

Dengue is a common mosquito-borne viral disease epidemic especially in tropical and sub-tropical regions where water sanitation is not substantially controlled. However, dengue epidemics sometimes occur in non-tropical urban cities with substantial water sanitary control. Using a mathematical model, we investigate what conditions can be important for a dengue epidemic to occur in an urban city such as Tokyo, where vectors are active only in summer and there are little number of vectors around hosts.

### Methods

The model, which is a modified Ross-Macdonald model, consists of two sets of host-vector compartments. The two sets correspond to high-risk and low-risk areas, and only hosts can move between them. Assuming that mosquitoes have constant activity for only 90 days, we assess five potential countermeasures: (1) restricted movement between the two areas, (2) insecticide application, (3) use of repellents, (4) vector control, and (5) isolation of the infected.

### Results

The basic reproduction number *R*
_{0} and the cumulative number of infected hosts for 90 days are evaluated for each of the five countermeasures. In the cases of Measures 2–5, the cumulative number of the infected for 90 days can be reduced substantially for small *R*
_{0} even if *R*
_{0}>1. Although *R*
_{0} for Measure 1 monotonically decreases with the mobility rates, the cumulative number of the infected for 90 days has a maximum at a moderate mobility rate. If the mobility rate is sufficiently small, the restricted movement effectively increases the number density of vectors in the high-risk area, and the epidemic starts earlier in the high-risk area than in the low-risk one, while the growth of infections is slow.

### Conclusions

Measures 2–5 are more or less effective. However, Measure 1 can have the opposite effect, depending on the mobility rates. The restricted movement results in the formation of a kind of core population, which can promote the epidemic in the entire population.

## Background

Dengue is a mosquito-borne viral infection. The number of the global incidence of dengue has grown dramatically in recent years, and 3900 million people, in 128 countries are under the risk of infection [1, 2]. The main vector for transmission of the dengue virus is *Aedes aegypti* (*A. aegypti*), and second less effective vector is *Aedes albopictus* (*A. albopictus*) [2–5]. Before 1970, only 9 countries suffered from dengue epidemic. However, now dengue is endemic in more than 100 countries [2, 6]. The vaccine for dengue fever has been under development, and therefore, in order to prevent the spread of dengue virus, it is important to focus on vector control [2, 7, 8].

According to Gubler, dengue epidemic in tropical and sub-tropical countries has been caused by dramatic population increase, urbanization and globalization [6]. The end of World War II era brought a rapid economic growth in many south tropic countries, which made urban area very congested. Moreover, the situation of inadequate housing and the few or no basic service such as water, sewer and waste management have created the ideal reproductive environment for mosquitoes. The crowded human communities and a large number of vectors have increased the risk of epidemic in these countries. Furthermore, globalization in recent decades enabled easy invasion of a disease. For example, airplanes coming from endemic regions have brought careers of infection to non-endemic countries.

On the other hand, cases of dengue in the non-tropical countries has been reported in past few years. In France and Croatia, the first cases occurred in 2010, and in Florida (United States) was reported in 2013 [2, 7]. Also dengue was reported in Japan after a lapse of 70 years [2]. As features of developed cities such as Tokyo, water sanitation system is well-ordered, and there are many urban green areas, which means the population of mosquitoes varies from place to place. In green area such as a park and a forest, the population of mosquitoes is large, which makes the risk of transmission of dengue high. However, in business districts and residential area, the population of mosquitoes is small, and the risk of infection is low.

In this study, we assess potential countermeasures against the dengue epidemic in a non-tropical urban city with substantial water sanitary control, using a modified Ross-Macdonald model. Here, we suppose that the non-tropical city is the region where mosquitoes are active only for a few months in a year. The Ross-Macdonald model is a mathematical model for mosquito-transmitted diseases [8], and we apply the model to the two sets of host-vector compartments that correspond to high-risk and low-risk areas. This model is a kind of two-patch metapopulation model. Metapopulation models of mosquito-transmitted diseases have been studied by several groups not only in two-patch systems [9, 10] but also in multi-patch systems [11–15], although most of those models are for malaria [10–15]. While both humans and mosquitoes are suppose to move between different patches in some models [14, 15], we assume in our model that only humans move between the two different areas. We investigate what conditions can be important for the dengue epidemic under the considered situation.

## Methods

### Model

We use a modified Ross-Macdonald model which consists of two sets of host-vector compartments (Fig. 1). Each set has susceptible-infected-recovered compartments for host and susceptible-infected compartments for vector. The two sets correspond to high-risk and low-risk areas, and only hosts can move between them. Dengue virus is classified into four serotypes, however, we here consider the situation where only one serotype is imported. The equations for the change in host and vector populations in the high-risk area are

where *S*
_{h}, *I*
_{h} and *R*
_{h} are the populations of susceptible, infected and recovered hosts, and *S*
_{v} and *I*
_{v} are those of susceptible and infected vectors, respectively. The superscripts H and L represent high-risk and low-risk areas, respectively. *γ* is the host recovery rate. *w*
_{LH} and *w*
_{HL} are mobility rates from the high-risk area to the low-risk one and that of the opposite direction, respectively. *μ*
_{b} and *μ*
_{m} are the birth and death rates of vector, respectively. The forces of infection are given by

where *a* is the biting rate, *b* and *c* are the transmission rates from vector to host and from host to vector, respectively, and *q* is the quarantine population of host. The total numbers of hosts and vectors in the high-risk area are \(N_{\mathrm {h}}^{\mathrm {H}}=S_{\mathrm {h}}^{\mathrm {H}}+I_{\mathrm {h}}^{\mathrm {H}}+R_{\mathrm {h}}^{\mathrm {H}}\) and \(N_{\mathrm {v}}^{\mathrm {H}}=S_{\mathrm {v}}^{\mathrm {H}}+I_{\mathrm {v}}^{\mathrm {H}}\), respectively. The equations for the change in host and vector populations in the low-risk area are

Here, the forces of infection are given by

where \(N_{\mathrm {h}}^{\mathrm {L}}=S_{\mathrm {h}}^{\mathrm {L}}+I_{\mathrm {h}}^{\mathrm {L}}+R_{\mathrm {h}}^{\mathrm {L}}\) and \(N_{\mathrm {v}}^{\mathrm {L}}=S_{\mathrm {v}}^{\mathrm {L}}+I_{\mathrm {v}}^{\mathrm {L}}\). In order to keep the host population constant in each area, we impose a condition about the mobility rates,

### Parameters

The parameters we use is summarized in Table 1. We assume that the primary case occurs in the central part of Tokyo. The parameters about population are extracted from the data of Meguro City (Meguro-ku), which is a typical urban city in Tokyo. The population of Meguro City is about *N*
_{h}=2.7×10^{5} (in 2015) [16]. The high-risk area is supposed to be the area where the number density of mosquitoes is high, e.g., parks. The low-risk area is supposed to be the area where it is low, e.g., business districts and residential areas. We allocate the populations in the high-risk and low-risk areas as \(N_{\mathrm {h}}^{\mathrm {H}}=N_{\mathrm {h}}\times 1/31\) and \(N_{\mathrm {h}}^{\mathrm {L}}=N_{\mathrm {h}}\times 30/31\), respectively, according to the proportion of the area of parks in Meguro City [17].

The host mobility rate from the low-risk area to high-risk one *w*
_{HL} for the base case is estimated from the number of visitors to parks [16]. That of the opposite direction *w*
_{LH} is given by Eq. (1).

The number of vectors per host in the high-risk area, \(r^{H}=N_{\mathrm {v}}^{\mathrm {H}}/N_{\mathrm {h}}^{\mathrm {H}}\), for the base case is evaluated from the biting density in a park in urban Tokyo and the attraction rage of a human bite [18, 19]. That in the low-risk area, \(r^{L}=N_{\mathrm {v}}^{\mathrm {L}}/N_{\mathrm {h}}^{\mathrm {L}}\), is given so that the basic reproduction number for an isolated low-risk area

is about 0.5.

The biting rate *a* for the base case is estimated by fitting to the data of the dengue epidemic in Tokyo 2014 [20]. Since the closure of a park was started after 19 days from the index case of the epidemic [21], the curve of the cumulative number of infected hosts is fitted to the cumulative cases during 19 days by a least square method. The result of the fitting is illustrated in Fig. 2.

### Countermeasures

We assess five potential countermeasures (Table 2): (1) restriction of the travel between the high-risk and low-risk areas, (2) insecticide application after the outbreak, (3) call on citizens to use mosquito repellents, (4) reducing carrying capacity of mosquitoes and (5) isolation of infected people. Measures 1 and 2 are implemented after 7 days from the primary case. For Measure 1, the mobility rates *w*
_{HL} and *w*
_{LH} change to take another value from the 7th day to the 90th day. The vector mortality rate is the sum of the natural mortality rate \(\mu _{\mathrm {m}}^{\text {nt}}\) and the excess mortality rate \(\mu _{\mathrm {m}}^{\text {ex}}\). For Measure 2, the excess mortality rate \(\mu _{\mathrm {m}}^{\text {ex}}\) takes a nonzero value during a certain period from the 7th day. Measure 3 is relevant to the biting rate *a*. When people use mosquito repellents, the biting rate is expected to decrease. Measure 4 is relevant to the numbers of vectors per host, *r*
^{H} and *r*
^{L}. Measure 5 is relevant to the quarantine proportion *q*. Here, isolation means to keep an infected person from mosquitoes.

Here, we note that each of populations in the high-risk and low-risk areas does not change with time. In other words, some people live in the high-risk area, and the others live in the low-risk area. Measure 1 is to reduce the traffic between the high-risk and low-risk areas. Although the parameters for the high-risk area are taken as those for parks, the high-risk area is not limited to parks. The high-risk area includes the area where the sanitary control of water is not substantial.

### Seasonality

In Tokyo today, mosquitoes are active only in summer and cannot be alive in winter. In this paper, we assume that mosquitoes are active for only 90 days in summer and that their activity is constant during the period as illustrated in Fig. 3.

## Results and discussion

### Basic reproduction number

The basic reproduction number is given by

where

The sensitivity analysis of *R*
_{0} is shown in Fig. 4. In Fig. 4(a), *R*
_{0} is plotted as the function of the ratio of each parameter to each value of the base case for Measures 1 (restricted movement), 3 (use of repellents) and 4 (vector control). Reduction of the biting rate and the number of mosquitoes per host reduce *R*
_{0} as illustrated by the green and blue curves. Since *R*
_{0} is linearly proportional to the biting rate, Measure 3 is especially effective.

The red curve in Fig. 4(a) implies that *R*
_{0} increases as the mobility rate becomes small. Actually, Eq. (2) is nearly equal to the basic reproduction number for an isolated high-risk area when *w*
_{HL}=*w*
_{LH}=0 since *r*
^{H}≫*r*
^{L} and *r*
^{L} is negligible. In other words, the restricted movement effectively increases the number density of vectors in the high-risk area.

In Figs. 4(b) and (c), *R*
_{0} is plotted as the functions of the vector mortality rate *μ*
_{m} and the quarantine proportion *q*, which are relevant to Measures 2 and 5, respectively. *R*
_{0} decreases when *μ*
_{m} or *q* increases.

### Cumulative number of infected hosts

The sensitivity analysis of the cumulative number of infected hosts for 90 days is shown in Fig. 5. In Fig. 5(a), the cumulative number of the infected is plotted as the function of the ratio of each parameter to each value of the base case for Measures 1 (restricted movement), 3 (use of repellents) and 4 (vector control). As illustrated by the green and blue curves, the reduction of the biting rate and the number of vectors per host substantially reduce the cumulative number of the infected for 90 days even when *R*
_{0}>1.

Although *R*
_{0} decreases monotonically with the mobility rates as shown in Fig. 4(a), the cumulative number of the infected for 90 days, i.e. the red curve in Fig. 5(a), has a maximum. When the mobility rates are large, the growth of infection is slow because the effective number density of vectors becomes small. Thus, a relatively small number of people are infected for 90 days. When the mobility rates are small enough, the growth of infection becomes fast in the high-risk area, however, the epidemic in the low-risk area does not immediately follow that in the high-risk area, and the cumulative number of the infected diminishes. Thus, the maximum of the cumulative number of the infected for 90 days appears at a moderate mobility rate. We discuss the details of effects of Measure 1 (restricted movement) later in the next subsection.

In Fig. 5(b), the cumulative number of the infected is plotted as the function of the implementation period of Measure 2 (insecticide application). When the excess mortality rate \(\mu _{\mathrm {m}}^{\text {ex}}=1.0\), only one day is effective enough to decrease the number of the infected substantially. When the excess mortality rate is half, the period to obtain the same result doubles.

In Fig. 5(c), the cumulative number of the infected is plotted as the function of the quarantine proportion *q*, which is relevant to Measure 5 (isolation). Although *R*
_{0}<1 for *q*>0.9 as shown in Fig. 4(c), the cumulative number of the infected for 90 days is small enough for *q*>0.6.

### Effects of restricted movement

The cumulative number of infected hosts for 90 days decreases when the mobility rates are sufficiently small, however, the epidemic starts earlier than the base case. That is illustrated in Fig. 6. The time dependence of the cumulative number of the infected in the base case is shown as the blue curve in Fig. 6(a). It starts to increase noticeably after 30 days. However, the blue curve in Fig. 6(b), where the mobility rates are 1/10 of the base case, starts to increase noticeably after 20 days.

The red curves in Fig. 6 show that the early start of the epidemic is induced by the epidemic in the high-risk area. The red and green curves are the rates of new infections of hosts in the high-risk and low-risk areas, respectively. In the base case, Fig. 6(a), the difference in day on which the peak of new infections appears is rather small between the two areas. However, in Fig. 6(b), where the mobility rates are 1/10 of the base case, the peak of the new infections in the high-risk area is more than 40 days earlier than that in the low-risk area.

Effects of restricted movement are summarized as follows. First, the effectively high number density of vectors accelerates the epidemic in the high-risk area. Then, a certain amount of infected hosts moves to the low-risk area although the mobility rates are small. Consequently, the epidemic starts early, while the growth of infections is slow.

### Effects of a core population

The nontrivial results for Measure 1 (restricted movement) is caused by the heterogeneous populations which are divided into the high-risk and low-risk areas. When the mobility rates between the two areas are small, the population in the high-risk area is a kind of core population. It is often suggested that targeting a core population is the key to prevent the spread of disease. For example, several works on sexually transmitted diseases suggested control or treatment targeting core members, who are highly active in sexual behavior, should be better than targeting the whole population [22, 23].

The earlier peak of new infections in the high-risk area than in the low-risk area, which is shown in Fig. 6(b), is also one of the effects caused by the core population. The cascade of infection from core to other populations is observed also in other diseases. For example, in the case of pandemic influenza, the burden of disease shifts from children to adults [24]. This means that new infections in the core population, which has many contacts, occur earlier than in the other population.

Heterogeneity is essential to form a core population. Although the effect of heterogeneity is often discussed in studies of metapopulation models, formation of a core population has not attracted attention in metapopulation models for mosquito-transmitted diseases. Lee and Castillo-Chavez assessed the dengue transmission dynamics in heterogeneous environment using a two-patch model [9]. In their study, however, the local basic reproduction number of each patch is taken to be larger than 1 in the base case, which is different from our situation. In studies of two-patch malaria models [10, 15], it is indicated that human movement can cause the persistence of malaria in a patch with a local basic reproduction number below 1. Although the situation is similar to ours, their focus is on the endemic condition caused by human movement. In contrast, our results for Measure 1 (restricted movement) give an insight that the restriction of human movement can sometimes cause a larger-size epidemic than that of a no-restriction case.

### Limitations

We consider only one serotype of dengue virus in the above. It is unlikely that several serotypes are imported at the same time into a non-tropical region. However, if epidemics of several serotypes were took into account, the model should be changed into a more complicated model.

In this paper, particular seasonality is assumed: Mosquitoes keep constant activity only for 90 days in a year. Since the parameter values are taken for the worst case, the results may be overestimated. This limitation is the compensation for the simplicity of the model. The length of the period when mosquitoes are active affects the cumulative number of infected hosts but not the basic reproduction number.

## Conclusions

We have assessed the five potential countermeasures against the dengue epidemic in a non-tropical urban city such as Tokyo, using a mathematical model in which compartments are divided into the high-risk and low-risk areas. Measures 2 (insecticide application), 3 (use of repellents), 4 (vector control) and 5 (isolation) are more or less effective, and substantial reduction of the cumulative number of infected hosts for 90 days can be expected for a small *R*
_{0} even if *R*
_{0}>1. However, Measure 1 (restricted movement) can have the opposite effect, depending on the mobility rates. Reduction of the mobility rates, namely, the restriction of travel between the high-risk and low-risk areas, results in the formation of a kind of core population. The epidemic in the core population (in the high-risk area) can promote that in the entire population when there are moderate contacts between them.

## References

Brady OJ, Gething PW, Bhatt S, Messina JP, Brownstein JS, Hoen AG, Moyes CL, Farlow AW, Scott TW, Hay SI. Refining the global spatial limits of dengue virus transmission by evidence-based consensus. PLoS Negl Trop Dis. 2012; 6(8):e1760.

(World Health Organization, editor.)2015. Dengue and severe dengue. No. 117 in WHO Fact sheet, World Health Organization. http://www.who.int/mediacentre/factsheets/fs117/en. Accessed 31 Dec 2015.

Gubler DJ. The global emergence/resurgence of arboviral diseases as public health problems. Arch of Med Res. 2002; 33:330–42.

WHO, et al.Global strategy for dengue prevention and control 2012–2020. World Health Organization. 2012.

Murray NEA, Quam MB, Wilder-Smith A. Epidemiology of dengue: past, present and future prospects. Clin Epidemiol. 2013; 5:299–309.

Gubler DJ. Dengue, urbanization and globalization: the unholy trinity of the 21st century. Trop Med Health. 2011; 39(4 Suppl):3–11.

Rodriguez-Roche R, Gould EA. Understanding the Dengue Viruses and Progress towards Their Control. BioMed Res Int. 2013; 2013:e690835. ISSN 2314-6133.

Smith DL, Battle KE, Hay SI, Barker CM, Scott TW, McKenzie FE. Ross, Macdonald, and a Theory for the Dynamics and Control of Mosquito-Transmitted Pathogens. PLoS Pathog. 2012; 8:e1002588.

Lee S, Castillo-Chavez C. The role of residence times in two-patch dengue transmission dynamics and optimal strategies. J Theor Biol. 2015; 374:152–64. ISSN 0022-5193.

Prosper O, Ruktanonchai N, Martcheva M. Assessing the role of spatial heterogeneity and human movement in malaria dynamics and control. J Theor Biol. 2012; 303:1–14.

Arino J, Ducrot A, Zongo P. A metapopulation model for malaria with transmission-blocking partial immunity in hosts. J Math Biol. 2012; 64(3):423–48.

Agusto F. Malaria drug resistance: The impact of human movement and spatial heterogeneity. Bull Math Biol. 2014; 76(7):1607–41. doi:10.1007/s11538-014-9970-6.

Auger P, Kouokam E, Sallet G, Tchuente M, Tsanou B. The Ross-Macdonald model in a patchy environment. Math Biosci. 2008; 216(2):123–31.

Rodríguez DJ, Torres-Sorando L. Models of infectious diseases in spatially heterogeneous environments. Bull Math Biol. 2001; 63(3):547–571.

Cosner C, Beier J, Cantrell R, Impoinvil D, Kapitanski L, Potts M, Troyo A, Ruan S. The effects of human movement on the persistence of vector-borne diseases. J Theor Biol. 2009; 258(4):550–60.

Statistics of Meguro City in Tokyo. http://www.city.meguro.tokyo.jp/gyosei/tokei/. Accessed 31 Dec 2015.

Proportions of the area of parks in tokyo. http://www.kensetsu.metro.tokyo.jp/kouen/kouenannai/kouen_menseki.html. Accessed 31 Dec 2015.

Tsuda Y, Kim KS. Ecology of mosquitoes inhabiting a park in urban Tokyo, Japan: Density of biting

*Aedes albopictus*and laboratory estimations of the residual longevity. Med Entomol Zool. 2012; 63:223–230. doi:10.1186/s12976-016-0039-0.Mogi M, Yamamura N. Estimation of the attraction range of a human bait for Aedes albopictus (diptera, culicidae) adults and its absolute density by a new removal method applicable to populations with immigrants. Res Pop Ecol. 1981; 23:328–43. ISSN 0034-5466. doi:10.1007/BF02515635.

Report on the cases of dengue fever in Japan (October 2014). http://www.mhlw.go.jp/stf/houdou/0000063557.html. Accessed 31 Dec 2015.

Seki N. Outbreak of dengue fever in Tokyo’s Yoyogi park and its vicinity, 2014. Infectious Agents Surveillance Report Vol.36 No.3, Infectious Disease Surveillance Center, National Institute of Infectious Diseases. 2015.

Yorke JA, Hethcote HW, Nold A. Dynamics and control of the transmission of gonorrhea. Sex Transm Dis. 1978; 5:51–56.

Newman MEJ. Spread of epidemic disease on networks. Phys Rev E. 2002; 66:016128. doi:10.1103/PhysRevE.66.016128.

Bansal S, Pourbohloul B, Hupert N, Grenfell B, Meyers LA. The shifting demographic landscape of pandemic influenza. PLoS ONE. 2010; 5:e9360. doi:10.1371/journal.pone.0009360.

Oki M, Yamamoto T. Simulation of the probable vector density that caused the Nagasaki dengue outbreak vectored by Aedes albopictus in 1942. Epidemiol Infect. 2013; 141:2612–2622. ISSN 1469–4409, doi:10.1017/S0950268813000447.

World Health Organization. 2011. Comprehensive guidelines for prevention and control of dengue and dengue haemorrhagic fever. Rev. and expanded. ed., no. 60 in SEARO Technical publication series. New Delhi: World Health Organization Regional Office for South-East Asia. ISBN 9789290223948 9290223944 9789290223870.

## Acknowledgments

The authors thank Yoshihiro Fujiya, Ayaka Kobayashi, Shinichi Sakakibara, Kentaro Koide and Bok Ryun Chang for sharing knowledge and fruitful discussions in the group work of the Summer School of Infectious Disease Modeling in 2015; and Hiroshi Nishiura for valuable suggestions and discussions.

## Author information

### Authors and Affiliations

### Corresponding author

## Additional information

### Competing interests

The authors declare that they have no competing interests.

### Authors’ contributions

HM derived the equation of the basic reproduction number; HM, KH and KK performed simulations and analyzed the data; IK, SU and KK wrote the manuscript; MS provided contribution for modeling and for discussion. All authors read and approved the final manuscript.

## Rights and permissions

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## About this article

### Cite this article

Masui, H., Kakitani, I., Ujiyama, S. *et al.* Assessing potential countermeasures against the dengue epidemic in non-tropical urban cities.
*Theor Biol Med Model* **13**, 12 (2016). https://doi.org/10.1186/s12976-016-0039-0

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/s12976-016-0039-0

### Keywords

- Dengue mathematical model
- Sensitivity analysis
- Core population