|Year : 2018 | Volume
| Issue : 8 | Page : 460-466
Dengue outbreaks in Taiwan, 1998-2017: Importation, serotype and temporal pattern
Department of Public Health, China Medical University, Taichung Taiwan 40402
|Date of Submission||10-Jun-2018|
|Date of Decision||15-Jul-2018|
|Date of Acceptance||20-Jul-2018|
|Date of Web Publication||31-Aug-2018|
Department of Public Health, China Medical University Taichung, 91 Hsueh-Shih Road, Taichung, Taiwan, 40402
Source of Support: None, Conflict of Interest: None
Objective: To ascertain the role of imported cases and serotypes on dengue outbreaks in Taiwan which have been sporadic yet highly volatile during the past two decades, exhibiting record-breaking magnitude in recent years. Methods: Confirmed case and serotype data from Taiwan Centers for Disease Control during 1998-2017 were fully examined, with fitting of weekly and daily case data of each city/county to a mathematical model to pinpoint the waves of cases and their locations. Moreover, we quantify the timing of turning point and transmission potential of each wave and determine its circulating serotype, to ascertain any pattern or connection between the variations in circulating serotypes and the magnitude/transmissibility of outbreak. Results: While the number of imported case increased steadily during past two decades, the yearly number of indigenous cases fluctuated wildly. Moreover, while yearly percentages of serotypes for imported cases remains steady, that of indigenous cases does not exhibit any clear pattern. There was at least one wave of reported cases somewhere in Taiwan every year from 1998 to 2015, except in 2016-2017. The effective reproduction number R for all waves in all locations ranged from 1.14 to 2.87, with the exception of two Tainan waves, in 2010 (3.95) and 2015 (6.84). Four major outbreaks of over 2000 cases reveal circulation of one dominant serotype. Conclusions: Correlation between imported cases and indigenous outbreak prove to be difficult to ascertain, even with the availability of serotype data. However, although there had been occasional co-circulation of serotypes in one location, and for some years with different serotypes circulating in different locations, all major outbreaks of over 2 000 cases during the past two decades are due to circulation of mainly a single serotype, perhaps indicating greater transmission potential with one dominating serotype.
Keywords: DENV, Serotype, Taiwan, Imported cases, Mathematical model, Reproduction number
|How to cite this article:|
Hsieh YH. Dengue outbreaks in Taiwan, 1998-2017: Importation, serotype and temporal pattern. Asian Pac J Trop Med 2018;11:460-6
|How to cite this URL:|
Hsieh YH. Dengue outbreaks in Taiwan, 1998-2017: Importation, serotype and temporal pattern. Asian Pac J Trop Med [serial online] 2018 [cited 2018 Sep 21];11:460-6. Available from: http://www.apjtm.org/text.asp?2018/11/8/460/240081
Foundation project: YHH is supported by funding from grants from Taiwan Ministry of Science and Technology (MOST 103-2314-B-039-010-MY3) and China Medical University-Taiwan (CMU106-S-02)
| 1. Introduction|| |
Incidence of dengue has grown dramatically around the world in recent decades, when dengue is ranked among the top re-emerging diseases, posing serious public health threat. In the past, the actual numbers of dengue cases were often severely underreported with many cases misclassified, mainly due to large portions of asymptomatic infections and difficulty in accurate diagnosis in early decades. Diagnosis of dengue can be challenging, highly dependent on timing of sampling, phase of infection, serotype and immune response which varies depending on whether the individual has a primary (i.e., first dengue or other flavivirus infection) or a secondary (i.e., had dengue or other flavivirus infection in past) infection. A recent study on prevalence of dengue estimates that 3.9 billion people in 128 countries are at risk for dengue infection.
There are four distinct but closely related serotypes of the virus that cause dengue, namely DENV1, DENV2, DENV3 and DENV4. Recovery from infection by one provides lifelong immunity against that particular serotype. However, cross-immunity to the other serotypes after recovery is not only partial but also temporary. Antibody-mediated enhancement (ADE) of dengue virus infection has been known to further complicate disease severity. Subsequent infections by other serotypes could actually increase the risk of developing severe dengue. Hyper-endemicity of multiple dengue virus serotypes in many countries has also contributed to the difficult challenge to predict and control dengue outbreaks.
In Taiwan, located in the tropical-subtropical region of the Northern Hemisphere, there had been many dengue fever/ dengue haemorrhagic fever outbreaks in the first half of 20th century,,, including a widespread outbreak island-wide in 1915-1916, when merchant ship from Southeast Asia brought the disease to Kaohsiung and subsequently spread via railroad to other cities on the island such as Keelung, Taipei and Taichung. It has been reported that several million people were infected during this outbreak. Another island-wide outbreak also occurred in 1942-1943,, during World War II, when Taiwan as a main transport hub for Japanese invasion in Southeast Asia which led to numerous imported dengue cases from affected regions.
However, no outbreak had occurred in the main island of Taiwan from 1944 until 1987,. From 1998, mainly initiated by imported cases, indigenous dengue outbreak has occurred in Taiwan mainly in the southern cities of Kaohsiung and Tainan every summer, except in 2004 and 2013 when most of cases were reported in Pingtung County in the southern tip of Taiwan. In particular, recording-breaking outbreaks occurred for two consecutive years in 20142015. The first of which in 2014 occurred mainly in Kaohsiung, while the 2015 outbreak, totaling more than 40 000 cases, struck both Tainan and Kaohsiung. Typically, Aedes albopictus mosquito is distributed throughout Taiwan while Aedes aegypti appears only in the tropical southern Taiwan, divided by the Tropic of Cancer which cuts across central Taiwan.
It has been proposed that the July 31, 2014 gas explosion in Kaohsiung had contributed significantly to the recording-breaking outbreak in Kaohsiung that year,,. Moreover, the Kaohsiung Rapid Transit System also played a role in spreading the disease within the metropolitan area. However, no certain explanation can be given for the 2015 explosive outbreak in Taiwan. Moreover, since the end of that outbreak there has been an extremely drastic drop in subsequent reported indigenous case number in 2016-2017 from a few hundred cases in 2016, down to only 10 reported indigenous cases in all of Taiwan in 2017. The sudden explosion of indigenous cases in 2014-2015 and its subsequent disappearance in 2016-2017 in Taiwan mirrors the unpredictability of dengue outbreaks in many countries in Asia and Latin America, albeit sometimes in a less dramatic fashion, which has left scientists with many theories yet no clear explanation for this unusual phenomenon. Dengue outbreaks in Taiwan has been typically initiated by imported cases in early summer, continuing through summer and spreading locally being boosted by summer international travels, peaking in late summer and early fall, and finally ending in winter as the temperature drops. The pattern repeats almost yearly, albeit with decidedly different magnitude from year to year.
In this work, we attempt to explore the relationship between dengue outbreak and local/imported serotype evolution, by using a mathematical model to pinpoint each wave of infections in Tainan and Kaohsiung during 1998-2017. We further quantify the transmission potential of each wave via its effective reproduction number R, to ascertain its temporal changes in relation to the circulating serotype(s) during this wave as well as the circulating serotypes(s) and transmission potential of the preceding waves. The focus is whether our results can offer some clues to the perplexing puzzle of emergence and magnitude of dengue outbreaks, for the purpose of predicting future outbreaks.
| 2. Methods and materials|| |
We make use of the Taiwan daily dengue confirmed case statistics and dengue serotyping data on the city/county level during the period of 1998-2017, respectively available from Taiwan Centers for Disease Control and Government Open Data Platform websites to generate weekly data, and to fit the Richards model, in order to determine the exact number of waves that had occurred for each year in each city/county in Taiwan.
2.2. Mathematical model
The Richards model, first proposed by ecologists to study biological growth, has been found to be useful in modeling the increase in cumulative case number in infectious diseases in recent years and has become a useful choice of modeling among modelers in the world,,,. It has also been found to be useful in modeling dengue outbreaks,,,,.
The Richards model, with C(t) denoting the cumulative number of reported dengue cases at time t, is given by the analytic formula:
where K is the total cumulative case number of a wave of cases, r is the per capita growth rate of cumulative case number, a is the exponent of deviation of cumulative case curve, and t¡ is time at which a turning point (or the peak) occurs, which signifies the exact moment of an upturn or downturn in the rate of increase for the cumulative case number of a wave of cases.
The Richards model is a phenomenological model which models the growth of cumulative case number. Since dengue outbreak is known to occur in waves, the model is most suitable to study temporal progression of dengue infections. Three model parameters of epidemiological importance are K, r, and the turning point t¡ of one single wave of the outbreak, which can be estimated by fitting the Richards model to the cumulative case curve of the outbreak, using any standard software with nonlinear least-squares (NLS) approximation subroutine, e.g., SAS (which is used in this work) or MATLAB.
We can subsequently compute the well-known basic reproduction number R0, the mean number of secondary human infections produced by an infective individual in a totally susceptible population in the absence of intervention measures. It is given by R0 = exp(rT), where r is the per capita growth rate estimated for a given wave and T is the serial interval, or the mean time interval from onset of one infected individual to the onset of his/her infectees. It has been shown mathematically that, given r, the expression exp(rT) provides an upper bound for basic reproduction number over any estimates of T obtained from all assumed distributions. However, in countries with sustained dengue outbreaks such as Taiwan, some level of herd immunity is present in the community, hence the estimate we obtain in this study is not the basic reproduction number, but rather the effective reproduction number R. See Hsieh and Chen for related discussions.
| 3. Results|| |
A summary of dengue cases and serotype data during 1998-2017 are given in [Table 1] & [Table 2], and [Figure 1] and [Figure 2]. Here we also provide the circulating serotype(s) in every city/county for each year between 1998-2017, using the definition that a serotype is “circulating” in a location in a given year if: (1) there are 10 or more indigenous cases with isolates of this serotype in this location during the year; and (2) this particular serotype also constitutes 35% or more of all indigenous serotype isolates obtained in this location during this year. Furthermore, the serotype is “dominant” in that year if it constitutes more than 50% of all isolates. There is at least one serotype circulating somewhere in Taiwan every year except in 1999-2000 and 2016-2017, when the respective yearly number of serotype isolates for indigenous cases was less than 10 [Table 1].
|Figure 1: Yearly percentages of serotype isolate results in Taiwan among indigenous dengue cases during 1998-2016.|
Click here to view
|Figure 2: Percentages of serotype isolate results in Taiwan among imported dengue cases during 1998-2016.|
Click here to view
The results of data fitting to the Richards model are provided in [Table 3], with time intervals for each wave occurring in each city during 1998-2017, along with the corresponding estimated effective reproduction number R with 95% confidence interval (CI). Through model fitting, there were at least one wave of dengue cases detected somewhere in Taiwan every year during 1998 to 2015. Interestingly, there is no wave pinpointed anywhere in Taiwan in 2016 or 2017. In some years, namely 1998 and 2013-2015, there were multiple waves in some locations. In 2011, waves were detected in all three southern cities/county: Kaohsiung, Tainan and Pingtung, plus the nearby island county of Penghu.
|Table 2: 1998-2017 dengue serotyping results and circulating serotypes (defined by the serotype with 10 serotyping results and 35% of all results in that year) in each location for indigenous and imported cases (n,%).|
Click here to view
|Table 3: Summary table for waves of weekly dengue cases with the Richards model fit during 1998-2008 in Taiwan.|
Click here to view
To further illustrate the possible connection between circulating serotype and transmission potential (via the effective reproduction number R), we provide in [Figure 3] and [Figure 4], for each year during 19982015 in Tainan and Kaohsiung, the percentages of serotypes and the largest estimate of R if there is more than one wave in a city that year. In some years, there were too few serotyping results in a city [Table 2]. We also note that in some years, there are no waves in either Tainan or Kaohsiung and hence no estimate for the effective reproduction number R.
|Figure 3: Yearly percentages of indigenous serotype isolate results in Taiwan during 1998-2016, with black diamond denoting largest effective reproduction number R of that year.|
Click here to view
|Figure 4: Yearly percentages of indigenous serotype isolate results in Kaohsiung during 1998-2016, with black diamond denoting largest effective reproduction number R of that year.|
Click here to view
| 4. Discussion|| |
4.1. Role of importation
While the number of indigenous cases in Taiwan has fluctuated greatly with large swings since it was first recorded in 1998, the number of imported cases exhibits a slow but steady increasing trend [Table 1], partly attributable to a similar pattern of increases in the number of international tourists to and from Taiwan during the past two decades, and also in the number of migrant workers in Taiwan from the Association of Southeast Asian Nations (ASEAN) countries with historically frequent dengue outbreaks.
It has been long conjectured that dengue outbreaks in Taiwan is related to outbreaks in neighboring Association of Southeast Asian Nations countries, especially after 1987 with the end of martial law which commenced an era of unrestricted tourism. King et al. show that three major dengue fever/dengue haemorrhagic fever outbreaks in Taiwan between 1981 and 1998 had statistically significant association with the increasing numbers of dengue cases in several Asian countries before or during these outbreaks in Taiwan, suggesting imported cases played an important role in indigenous outbreaks in Taiwan. Shang et al. conclude that imported dengue cases could initiate indigenous outbreaks in Taiwan, albeit only under suitable climate conditions. However, our results indicate that, while such relationship is likely to exist, it is by no means simply quantifiable, nor can it be easily introduced in an early warning system with a set of quantities that are widely different, such as vector indices and climatologic factors.
4.2. Role of serotype
Serotyping results of indigenous cases reveal that all four serotypes had been reported in Taiwan since 1998, with each serotype taking the role of a circulating serotype during at least some years, although DENV2 and DENV1 are clearly the most frequently circulating serotypes. In every year except 2010, the circulating serotype is also the dominant serotype (>50% of the total results) of that year. In 2010, Tainan had an outbreak of DENV 4 while Kaohsiung had an outbreak of both DENV2 and DENV3. Subsequently as all three serotypes coexisted more or less evenly with no one dominant serotype.
Except in 1999-2010, 2016-2017 (mainly due to a scarcity of positive isolates), as well as in 2010 (when both DENV2 and DENV3 were circulating in Kaohsiung), there was typically one highly dominant (>95%) serotype every year, including all four major outbreaks with more than 2000 cases - 2002 (DENV2 in Tainan), 2014 (DENV1 In Kaohsiung) and 2015 (DENV2 in Tainan and Kaohsiung), perhaps indicating the greater transmission potential of one dominating serotype.
Coincidentally, in 2002 the last serotype to emerge in Taiwan, DENV4, was first found in one imported case in Pingtung County. Furthermore, [Figure 2] reveals that since 2002 a clear and consistent pattern in the serotypes of imported cases in Taiwan emerges, with frequency of serotypes in the order of mostly DENV1 and DENV2, followed by DENV3, with DENV4 the least but being present every year. Subsequently, there seems to be no noticeable correlation between the serotypes of imported cases and circulating serotypes of indigenous cases in Taiwan. However, how the four serotypes interact in antibody-mediated immunity and enhancement remains a mystery. One might speculate that this lack of knowledge has played a significant role in our difficulty in predicting dengue outbreak.
4.3. Transmission potential
Data fitting with the Richards model results in at least one wave somewhere in Taiwan each year from 1998 to 2015 [Table 3], but there had been very few indigenous cases reported in Taiwan after January of 2016. The mean estimate of the effective reproduction number R for all waves ranges consistently between 1.21 and 2.87, except in Tainan in 2010 (R=3.95) and in 2015 (R=6.84). We note the exceptionally high estimate in 2015, when a historically large number of cases were reported, and the wave of cases fitted was unusual in both its early start in May (week 21) and its length (over 30 weeks). We speculate that the early starting point might contribute to a high initial growth rate and subsequently a high estimate for R.
Many studies have shown a significant correlation between dengue outbreak and many factors, including climate, serotype, imported cases, timing, geographic location, human mobility, etc. Subsequently, there are ample studies endeavoring to predict future dengue outbreak. However, a scientifically reliable prediction or early warning system of dengue outbreak still very much eludes us. The reason is that dengue outbreak often occurs due a combination of the factors mentioned above, not the least of which is the interaction of serotypes in hosts. As technology advances that enable the collection of data pertaining to these factors, the primary challenge becomes that of consolidation of these different datasets in a multi-layered model constructed in a suitable platform to make the result explicitly understandable to policy makers. Multi-layered information on infectious diseases pertaining to its epidemiology, etiology, immunology, and related climatology/geography/sociology, must be consolidated into one single model incorporating different types of data, in order to truly ascertain its potential threat to humans.
Conflict of interest statement
We declare that we have no conflict of interest.
Yu-Han Li contributed to the data analysis of 2013-2016 dengue outbreaks. YHH is supported by funding from grants from Taiwan Ministry of Science and Technology (MOST 103-2314-B-039-010-MY3) and China Medical University-Taiwan (CMU106-S-02).
| References|| |
Braga IA., Valle D. Aedes aegypti:
histórico do controle no Brasil. [Aedes aegypti:
history of control in Brazil]. Epidemiologia e Serviços de Saúde
Brady OJ, Gething PW, Bhatt S, Messina JP, Brownstein JS, Hoen AG, et al. Refining the global spatial limits of dengue virus transmission by evidence-based consensus. PLoS Negl Trop Dis
Flipse J, Diosa-Toro MA, Hoornweg TE, van de Pol DPI, Urcuqui-Inchima S, Smit JM. Antibody-dependent enhancement of dengue virus infection in primary human macrophages; balancing higher fusion against antiviral responses. Sci Rep
Akashi K. Tainan dengue epidemic study. Taiwan J Med
(7): 767777. (In Japanese).
Clarke EJ Jr, Suitor EC Jr, Jenkin HM. A serologic survey for arbovirus antibodies in inhabitants of a southwest Taiwan village. Trop Geogr Med
Kuno G. Research on dengue and dengue-like illness in East Asia and the Western Pacific during the first half of the 20th century. Rev Med Virol
Ooi T. dengue epidemic in Taiwan: Epidemiology and clinical symptoms in Kaohsiung. Taiwan J Med
734-741. (In Japanese).
Ko YC. Epidemiology of dengue fever in Taiwan. Gaoxiong Yi Xue Ke Xue Za Zhi
1-11. (In Chinese).
Tung TTH, Tsai KH, King CP, Huang SYJ, King CC. The relationships among socio-political/environmental changes, prevention and control strategies, and epidemics of dengue fever in Taiwan—future prospects. Taiwan Med J
Hsieh, YH. Ascertaining the impact of catastrophic events on dengue outbreak: The 2014 gas explosions in Kaohsiung, Taiwan. PLoS one
Wang SF, Chang K, Lu RW, Wang WH, Chen YH, Chen M, et al. Large dengue virus type 1 outbreak in Taiwan. Emerg Microbes & Infect
Yang CF, Hou JN, Chen TH, Chen WJ. Discriminable roles of Aedes aegypti
and Aedes albopictus
in establishment of dengue outbreaks in Taiwan. Acta Tropica
Wang SF, Wang WH, Chang K, Chen YH, Tseng SP, Yen CH, et al. Severe dengue fever outbreak in Taiwan. Am J Trop Med Hyg
Sanna M, Hsieh YH. Ascertaining the impact of public rapid transit system on spread of dengue in urban settings. Sci Total Environ
Taiwan Centers for Prevention and Control. TCDC: Dengue surveillance data; Taiwan CDC Notifiable Diseases Surveillance System (NDSS) [Online]. Available from: http://nidss.cdc.gov.tw/ch
[Accessed on 1st July 2018].
Lin CC, Huang YH, Shu PY, Wu HS, Lin YS, Yeh TM, et al. Characteristic of dengue disease in Taiwan: 2002-2007. Am J Trop Med Hyg
Taiwan Open Data. Daily confirmed dengue case data; Taiwan Government Open Database[Online]. Available from: https://data.gov.tw/ dataset/21025
[Accessed on 1st July 2018].
Richards FJ. A flexible growth function for empirical use. J Exp Bot
Hsieh YH, Lee JY, Chang HL. SARS epidemiology modeling. Emerg Infect Dis
Hsieh YH, Cheng YS. Real-time forecast of multiphase outbreak. Emerg Infect Dis
Mostaço-Guidolin LC, Greer A, Sander B, Wu J, Moghadas SM. Variability in transmissibility of the 2009 H1N1 pandemic in Canadian communities. BMC Res Notes
Chan CH, Tuite AR, Fisman DN. Historical epidemiology of the second cholera pandemic: relevance to present day disease dynamics. PloS One
Hsieh YH, Ma S. Intervention measures, turning point, and reproduction number for dengue, Singapore, 2005. Am J Trop Med Hyg
Hsieh YH, Chen CWS. Turning points, reproduction number, and impact of climatological events for multi-wave dengue outbreaks. Trop Med & Int Health
Hsieh YH, de Arazoza H, Lounes R. Temporal trends and regional variability of 2001-2002 Multi-wave DENV-3 epidemic in Havana City: Did hurricane michelle contribute to its severity? Trop Med Internat Health
Sanna M, Wu J, Zhu Y, Yang Z, Lu J, Hsieh YH. Spatial and temporal characteristics of 2014 dengue outbreak in Guangdong, China. Sci Rep
Wallinga J, Lipsitch M. How generation intervals shape the relationship between growth rates and reproductive numbers. Proc Royal Socie London
B: Biol Sci
King CC, Wu YC, Chao DY, Lin TH, Chow L, Wang H T, et. al. Major epidemics of dengue in Taiwan in 1981-2000: related to intensive virus activities in Asia. Dengue Bulletin
Shang CS, Fang CT, Liu CM, Wen TH, Tsai KH, King CC. The role of imported cases and favorable meteorological conditions in the onset of dengue epidemics. PLoS Negl Trop Dis
Chang FS, Tseng YT, Hsu PS, Chen CD, Lian IB, Chao DY. Re-assess vector indices threshold as an early warning tool for predicting dengue epidemic in a dengue non-endemic country. PLoS Negl Trop Dis
[Figure 1], [Figure 2], [Figure 3], [Figure 4]
[Table 1], [Table 2], [Table 3]