A simple periodic-forced model for dengue fitted to incidence data in Singapore

https://doi.org/10.1016/j.mbs.2013.04.001Get rights and content

Highlights

  • Fit of a periodic model to dengue incidence data in Singapore from 2003 to 2007.

  • Estimation of the parameters governing the vector population dynamics.

  • Estimation of the impact of vector-control.

  • Close relationship between vector density and temperature.

  • Hypothesis explaining partially the 2007 dengue outbreak in Singapore.

Abstract

Dengue is the world’s major arbovirosis and therefore an important public health concern in endemic areas. The availability of weekly reports of dengue cases in Singapore offers the opportunity to analyze the transmission dynamics and the impact of vector control strategies. Based on a previous model studying the impact of vector control strategies in Singapore during the 2005 outbreak, a simple vector-host model accounting for seasonal fluctuation in vector density was developed to estimate the parameters governing the vector population dynamics using dengue fever incidence data from August 2003 to December 2007. The impact of vector control, which consisted principally of a systematic removal of actual and potential breeding sites during a six-week period in 2005, was also investigated. Although our approach does not account for the complex life cycle of the vector, the good fit between data and model outputs showed that the impact of seasonality on the transmission dynamics is highly important. Moreover, the periodic fluctuations of the vector population were found in phase with temperature variations, suggesting a strong climate effect on the vector density and, in turn, on the transmission dynamics.

Introduction

Dengue is the major arbovirosis in the world with more than 50 million dengue fever cases per year [1], [2]. Although most dengue cases remain asymptomatic or present mild symptoms, Dengue Hemorrhagic Fever (DHF) and Dengue Shock Syndrome (DSS) are the most severe expressions with case-fatality rates (CFR) generally close to 1% [3], [4]. Four immunologically distant serotypes were identified and coexist in many hyperendemic areas [5], [6]. The first dengue cases in Singapore were described in the 1960s and dengue became rapidly a major cause of childhood death. The Singapore government initiated a vector control programme in 1968 resulting in a reduction of the premises index, i.e. the percentage of inspected premises found to have containers with Aedes aegypti larvae or pupae, from 16% to 2% in 1973. Since then, the premises index has remained around 2% and Singapore is now recognized as having one of the most effective vector control policies in the world [7], [8]. Despite relatively low dengue activity over more than 15 years, Singapore experienced several dengue outbreaks since the 1990s which might be due to different factors [8]. The largest outbreak occurred in 2005, with more than 14,000 identified symptomatic cases. During this epidemic, emergency vector-control measures, mainly based on removal of breeding sites, were set up, allowing for a rapid decrease of dengue case notifications. However, another large outbreak occurred in 2007, after one year with low dengue activity.

A wide range of mathematical models was developed to study the transmission process of the dengue virus [9]. Most of these models were derived from the simple vector host model originally proposed by Bailey [10], with different adaptations (e.g. age-structure [11], [12], detailed vector lifecycle [13], [14], [15], [16]; multi-serotypes [17], [18], [19], [20]; etc). However, following the multi-strain model of Ferguson et al. [21], several models considered direct transmission between hosts, assuming the time scale for transmission to be sufficiently short and the mosquito population sufficiently dense to ignore the vector component [22], [23], [24], [25]. Although dengue outbreaks often exhibit 3- to 5-year inter-epidemic periods, which could be reproduced with multi-serotype host-to-host transmission models, intra-annual fluctuations of dengue incidence were shown to be linked with seasonality, acting on the vector population [14], [17], [20], [26]. Moreover, Yang and Ferreira [16] showed that seasonality can also influence the effectiveness of vector strategies by identifying the optimal period to initiate larvicide or adulticide application.

Burattini et al. [14] proposed a transmission model assuming seasonal variations in the development rate from eggs to adult stages. Moreover, in addition to a periodic forcing governing the seasonal fluctuation of the vector population, the “egg carrying capacity” was assumed to increase linearly with time in relation with the observed increase in temperature. Using this framework, Burattini et al. tested the impact of different control strategies and reproduced qualitatively the dengue incidence data in Singapore from 2003 to 2006. This model was further used by Massad et al. who assumed that the presence of haze in Singapore could impact on mosquito behavior, explaining the reduction in dengue activity in 2006 [27]. Although this model provided a good qualitative fit to data, the comparison of dengue incidence time series with the pollution index showed no evidence of any relationship between dengue activity and haze [28], [29].

The aim of the present study was to evaluate the vector population dynamics from dengue fever incidence data. As suggested by Burattini et al. [14], the incidence data in Singapore offer the opportunity to quantify the impact of emergency control measures. For that purpose, we focused our analysis on a restricted time-window running from August 2003 to December 2007. To avoid overparametrization we developed a simplified version of Burattini’s model, accounting exclusively for the adult vector population with a periodic recruitment rate. Although the emergency vector control strategy carried out in September 2005 focused principally on breeding site removal, we evaluated the impact on the adult mosquito population. Finally, the analysis of data post-intervention permitted to formulate a hypothesis to explain the dynamics of infection observed after the period of control.

Section snippets

Model formulation

Based on the model of Burattini et al. [14], the host population, denoted by the index h, is divided into three compartments according to their status : susceptible (Sh), infectious (Ih) and (Rh) recovered individuals. The intrinsic incubation period (IIP) was not represented in the model because of its relatively low impact on the transmission dynamics, as shown by Newton and Reiter [30]. Conversely, the extrinsic incubation period (EIP) lasts between 10 and 12 days for an average mosquito

Results

In the period prior to emergency vector control, the peak in dengue incidence was found to occur in October each year from 2003 to 2005, showing a lag of 16.8 weeks with the maximal vector density (Fig. 1). Studies on the impact of climate variability on dengue incidence in Singapore also reported a correlation between temperature and dengue incidence with a possible time lag up to 20 weeks [35], [47]. Recently, Hii et al. analyzed dengue incidence in Singapore using weather data highlighting an

Conclusion

The aim of the present study was to estimate the parameters governing the adult vector population dynamics from dengue fever case notifications in Singapore. For that purpose, a simple vector-host transmission model was adapted from a previous model, which reproduced Singaporean data between 2003 and 2006 qualitatively well. Burattini et al. developed a vector-host transmission model including eggs and adult stages for the mosquito population [14]. Seasonality was incorporated into the model

Acknowledgements

The authors want to thank the two referees and the associate editor for their comments, which greatly improved the manuscript. NH also acknowledges support from the UA Scientific Chair in evidence based vaccinology.

References (53)

  • M. Derouich et al.

    Dengue fever: mathematical modelling and computer simulation

    Appl. Math. Comput.

    (2006)
  • L. Billings et al.

    Instabilities in multiserotype disease models with antibody-dependent enhancement

    J. Theor. Biol.

    (2007)
  • T.W. Scott et al.

    Feeding strategies of anthropophilic mosquitoes result in increased risk of pathogen transmission

    Trends Parasitol.

    (2012)
  • E.E. Ooi et al.

    Dengue seroepidemiology in Singapore

    The Lancet

    (2001)
  • M. Aguiar et al.

    The role of seasonality and import in a minimalistic multi-strain dengue model capturing differences between primary and secondary infections: complex dynamics and its implications for data analysis

    J. Theor. Biol.

    (2011)
  • D.J. Gubler

    Dengue and dengue hemorrhagic fever

    Clin. Microbiol. Rev.

    (1998)
  • J.L. Kyle et al.

    Global spread and persistence of dengue

    Annu. Rev. Microbiol.

    (2008)
  • E.E. Ooi et al.

    Dengue prevention and 35 years of vector control in Singapore

    Emerg. Infect. Dis.

    (2006)
  • M. Andraud et al.

    Dynamic epidemiological models for dengue transmission: a systematic review of structural approaches

    PLoS ONE

    (2012)
  • N.T.J. Bailey

    The Mathematical Theory of Infectious Diseases and its Applications

    (1975)
  • M.N. Burattini et al.

    Modelling the control strategies against dengue in Singapore

    Epidemiol. Infect.

    (2008)
  • L.M. Bartley et al.

    The seasonal pattern of dengue in endemic areas: mathematical models of mechanisms

    Trans. R. Soc. Trop. Med. Hyg.

    (2002)
  • E. Chikaki, H. Ishikawa, A dengue transmission model in Thailand considering sequential infections with all four...
  • H.J. Wearing et al.

    Ecological and immunological determinants of dengue epidemics

    Proc. Natl. Acad. Sci. USA.

    (2006)
  • N. Ferguson et al.

    The effect of antibody-dependent enhancement on the transmission dynamics and persistence of multiple-strain pathogens

    Proc. Natl. Acad. Sci. USA

    (1999)
  • S. Bianco et al.

    Asymmetry in the presence of migration stabilizes multistrain disease outbreaks

    Bull. Math. Biol.

    (2011)
  • Cited by (41)

    • Spatially weak syncronization of spreading pattern between Aedes Albopictus and dengue fever

      2022, Ecological Modelling
      Citation Excerpt :

      For example, statistical models have widely been developed by defining an empirical relationship between the actual occurrence and a set of hypothesized underlying covariates (e.g., Simon et al., 2002; Aström et al., 2012), and used to predict the distribution or abundance of dengue vectors (Caminade et al., 2012; Monaghan et al., 2018; Liu-Helmersson et al., 2019a, 2019b; Kraemer et al., 2019; Messina et al., 2019). Furthermore, mechanistic models have also been proposed to describe the dynamic behavior of mosquito-borne pathogen transmission (e.g., Macdonald, 1956; Anderson et al., 1991; Massad et al., 2012), and used to predict the pattern and distribution of dengue under climate change in terms of vectorial capacity or basic reproduction number (e.g., Jetten et al., 1997; Patz et al., 1998; Medlock et al., 2006; Andraud et al., 2013; Liu-Helmersson et al., 2014, J. 2016; Mordecai et al., 2017; Lee et al., 2018; Ryan et al., 2019; Paul et al., 2021). However, among all these studies, the impact of climate change on the duration of the dengue transmission process has not been considered when assessing DF's future pattern.

    • A time-periodic and reaction–diffusion Dengue fever model with extrinsic incubation period and crowding effects

      2020, Nonlinear Analysis: Real World Applications
      Citation Excerpt :

      Not only spatial heterogeneity, the environmental factors influencing Dengue fever also manifest in temporal heterogeneity. It has been shown that climate change, such as temperature and rainfall, intervene sensitively the distribution, reproduction and development of the Aedes mosquito vector [17,18], which facilitates relative research on the seasonal dynamics of Dengue fever [1,9,13,19]. It is worthy mentioning that the above factors actually also affect other vector-borne diseases like malaria [20–22].

    • Learning the seasonality of disease incidences from empirical data

      2019, Ecological Complexity
      Citation Excerpt :

      It addresses some conditions under which a periodic solution is locally asymptotically stable, which can bear a relation to the basic reproductive number from the autonomous counterpart or bear a definition of the true basic reproductive number. Some epidemic models have previously been introduced together with periodic parameters, with Tian and Wang (2014), Rocha et al. (2015) or without the analysis of periodic solutions (Aguiar et al., 2011; 2013; Andraud et al., 2013). One apparent issue from those with the analysis is that their parameters are never calibrated using real field data.

    View all citing articles on Scopus
    1

    Address: Centre for Health Economics Research and Modelling Infectious Diseases (CHERMID), Vaccine & Infectious Disease Institute, University of Antwerp, Campus Drie Eiken, Universiteitsplein 1, 2610 Antwerpen, Belgium.

    2

    Address: Interuniversity Institute for Biostatistics and statistical Bioinformatics, Hasselt University, Agoralaan 1, Building D, 3590 Diepenbeek, Belgium

    3

    Address: Faculty of Medicine, The University of New South Wales, Sydney, NSW 2052, Australia

    View full text