http://orcid.org/0000-0003-3402-1670
Correction
2 May 2018: Prayitno A, Taurel AF, Nealon J, Satari HI, Karyanti MR, et al. (2018) Correction: Dengue seroprevalence and force of primary infection in a representative population of urban dwelling Indonesian children. PLOS Neglected Tropical Diseases 12(5): e0006467. https://doi.org/10.1371/journal.pntd.0006467 View correction
Figures
Abstract
Background
Indonesia reports the second highest dengue disease burden in the world; these data are from passive surveillance reports and are likely to be significant underestimates. Age-stratified seroprevalence data are relatively unbiased indicators of past exposure and allow understanding of transmission dynamics.
Methodology/Principal Findings
To better understand dengue infection history and associated risk factors in Indonesia, a representative population-based cross-sectional dengue seroprevalence study was conducted in 1–18-year-old urban children. From October to November 2014, 3,210 children were enrolled from 30 geographically dispersed clusters. Serum samples were tested for anti-dengue IgG antibodies by indirect ELISA. A questionnaire investigated associations between dengue serologic status and household socio-demographic and behavioural factors. Overall, 3,194 samples were tested, giving an adjusted national seroprevalence in this urban population of 69.4% [95% CI: 64.4–74.3] (33.8% [95% CI: 26.4–41.2] in the 1–4-year-olds, 65.4% [95% CI: 69.1–71.7] in the 5–9-year-olds, 83.1% [95% CI: 77.1–89.0] in the 10–14-year-olds, and 89.0% [95% CI: 83.9–94.1] in the 15–18-year–olds). The median age of seroconversion estimated through a linear model was 4.8 years. Using a catalytic model and considering a constant force of infection we estimated 13.1% of children experience a primary infection per year. Through a hierarchical logistic multivariate model, the subject’s age group (1–4 vs 5–9 OR = 4.25; 1–4 vs. 10–14 OR = 12.60; and 1–4 vs 15–18 OR = 21.87; p<0.0001) and the number of cases diagnosed in the household since the subject was born (p = 0.0004) remained associated with dengue serological status.
Conclusions/Significance
This is the first dengue seroprevalence study in Indonesia that is targeting a representative sample of the urban paediatric population. This study revealed that more than 80% of children aged 10 years or over have experienced dengue infection at least once. Prospective incidence studies would likely reveal dengue burdens far in excess of reported incidence rates.
Author summary
Indonesia reported to the WHO the world’s second highest average number of dengue cases and the highest in Asia from 2004 to 2010. These passive surveillance reports vary widely within the country and are likely to be a severe under-estimation of the full disease burden as frequently only dengue haemorrhagic fever is captured. Understanding the intensity of dengue virus transmission and associated risk factors nationwide is necessary to guide and prioritize appropriate prevention and control measures against dengue disease, especially considering the availability of the first dengue vaccine and recent recommendations for its use in areas of high endemicity, as measured by seroprevalence and other indicators. Age-stratified seroprevalence data provide robust estimates of past exposure and can inform on transmission intensity. Therefore, we conducted a seroprevalence study of anti-dengue IgG antibodies in a representative sample of urban-dwelling Indonesian children. We found an overall dengue seroprevalence of 69.4% with half of the children having been infected at least once by the age of 5 years. Age of the subject and the number of dengue cases diagnosed in the household were associated with serostatus. These results confirm the high dengue disease burden in Indonesia and the urgency of implementation of effective prevention and control measures.
Citation: Prayitno A, Taurel A-F, Nealon J, Satari HI, Karyanti MR, Sekartini R, et al. (2017) Dengue seroprevalence and force of primary infection in a representative population of urban dwelling Indonesian children. PLoS Negl Trop Dis 11(6): e0005621. https://doi.org/10.1371/journal.pntd.0005621
Editor: Ricardo E. Gürtler, Universidad de Buenos Aires, ARGENTINA
Received: January 2, 2017; Accepted: May 3, 2017; Published: June 15, 2017
Copyright: © 2017 PRAYITNO et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This work was supported by Sanofi Pasteur. The sponsor participated in the study design, data collection and analysis, decision to publish, and preparation of the manuscript. The interpretation of the data and the opinions expressed in the manuscript are those of the authors. Manuscript editing, assistance and support for the submission provided by inScience Communications, Springer Healthcare, funded by Sanofi Pasteur.
Competing interests: I have read the journal's policy and the authors of this manuscript have the following competing interests: AP, HIS, MRK, RS, SS, HG, BEM, RTS, SRH have been clinical trial and/or study investigators for, and received associated payments from, Sanofi Pasteur; AFT, JN, MS and AB are employees of Sanofi Pasteur
Introduction
Dengue is an arbovirus transmitted to humans via the bites of infected Aedes mosquitoes. It is the most rapidly spreading mosquito-borne viral disease with a global incidence that has increased 30-fold over the last 50 years [1]. While reliable burden estimates remain elusive, two studies have estimated the global symptomatic disease burden to be 96 million and 58.4 million cases/year, with 70–80% of cases occurring in the Asia-Pacific region [2, 3]. Traditionally an urban disease, dengue disease is increasingly reported in rural areas and its geographic range has expanded to more than 125 tropical countries [1]. There is no specific antiviral treatment; clinical management is focused on careful fluid management and detection of early warning signs of severe disease. Historically, prevention measures have focused on vector control, education and behavioural changes to reduce interactions between humans and vector mosquitoes [4, 5]. Improved clinical management and public awareness have contributed to declining case fatality rates to below 1% in most countries [1]. While this represents important progress, overall dengue incidence continues to rise and fatalities remain unacceptably high, suggesting that traditional control approaches are not sufficient. Vector control measures are important yet operationally challenging, of variable effectiveness and costly to sustain [6]. Routine vaccination is becoming a reality: several dengue vaccines are at different stages of clinical development [7] and a chimeric tetravalent vaccine from Sanofi Pasteur is being licensed in an increasing number of countries in Latin America and Asia [7, 8]. In this new era of dengue as a vaccine-preventable disease, an accurate understanding of disease burden and transmission patterns will be essential to inform vaccine policy decisions.
Dengue is hyper-endemic with frequent epidemic cycles in Indonesia. The disease is most common in urban areas and in recent years has reportedly spread to smaller, more rural villages. Reported incidence remains highest in children 1–15 years of age, but since the 1980s incidence in persons over 15 years of age has gradually increased [9, 10]. Reporting of dengue haemorrhagic fever (DHF) is mandatory in Indonesia and the country typically reports the highest number of cases in the WHO Southeast Asia Region [1]. Between 2001 and 2011, there was an average of 94,564 reported cases and between 472 and 1,446 reported deaths per year [1, 11]. Dengue disease reporting is acknowledged by Indonesian experts to be incomplete and to vary widely between provinces, with reported incidence rates ranging from 2.2 to 168.5 cases per 100,000 inhabitants in 2013 [12].
An improved understanding of dengue epidemiology, burden and its dynamic characteristics are important for public health planning. Seroprevalence studies in healthy volunteers provide information on infection history in the population, from which inferences about disease burden may be drawn. Since age reflects duration of exposure, age-stratified data provide insights into transmission dynamics [13–17]. There is a lack of dengue seroepidemiological data from Indonesia and no previous study has used a population representative sample of urban Indonesian children [18–20]. This is a particularly important gap as it will provide information on whether the variations in reported incidence from different Indonesian provinces are reflective of underlying transmission dynamics or to the result of the reporting or surveillance practices employed. We conducted a seroprevalence study in urban-dwelling Indonesian children to improve understanding of dengue epidemiology and infection risk factors and inform future dengue vaccine policy decisions.
Methods
The present study is reported according to STrengthening the Reporting of OBservational studies in Epidemiology (STROBE) recommendations (supporting information file).
Ethic statement
The protocol was reviewed and ethical approval was obtained from the Health Research Ethics Committee of Faculty of Medicine of University of Indonesia.
Study area
Indonesia is the largest country in Southeast Asia, with an area of 1.91 million km2. The country has a population of 252.2 million living on five main islands and four archipelagos (>17,000 islands) administratively divided into 34 provinces [21]. In 2014/2015, approximately 60% of Indonesians were living on the island of Java and 53.3% lived in urban areas [21, 22]. Indonesia is divided into five administrative levels: provinces (n = 34), regencies (n = 416), cities (n = 98), subdistricts (n = 7,024), and villages (n = 81,626). Villages are considered either as rural (desa) or urban (kelurahan) based on population density, percentage of agricultural household and number of urban facilities such as schools and hospitals [21, 23].
Sampling design
A population-based cross-sectional study design was adapted from the World Health Organization (WHO) Expanded Program on Immunization (EPI) cluster survey method. This approach considers 30 clusters as an adequate number for their means to be normally distributed, thus permitting statistical theory based on the normal distribution to be used to analyse the data [24, 25]. Based on the probability proportional to population size, 30 urban subdistricts were selected using demographic data from 2009 or 2010, provided by the Sub-Directorate of Statistical Services and Promotion, Statistics Indonesia.
The geographical coordinates of Indonesian administrative units were retrieved from the Global Rural-Urban Mapping Project, maintained by the Socioeconomic Data and Applications Center [26]. Provinces were listed based on their mean geographical coordinates from West to East (Fig 1) and the cumulative urban population of their subdistricts was calculated using 2010 population data. To ensure the population of clusters was sufficient to enrol the desired sample, a minimum population of 1,000 persons per subdistrict was defined and any smaller subdistricts were removed from the list. The first cluster was selected by generating a random number between 1 and 1/30th of the total urban population, using Epi Info Version 7, and selecting the first subdistrict for which the cumulative population was superior or equal to this random number. Subsequent clusters were selected by adding 1/30th of the urban population to the random number and selecting the first corresponding subdistrict for which cumulative population was higher or equal so that:
Provinces with at least one site are coloured in grey. Developed with QGIS 1.8.0 using data from Global Rural-Urban Mapping Project18.
The 30 subdistricts selected by this method are listed in Appendix 1. Each subdistrict in Indonesia contains one main health centre (puskesmas kecamatan) whose catchment area was the site of the study. Households in the five neighbourhood associations located closest to the health centre (each comprising 30–50 households, giving a total of 150–250 households) were eligible to participate in the study. Household visits were conducted, inviting one child from each household to participate, until the sample size was reached. A table indicating the required number of children from each of four age groups was provided to the health centre study teams. If a household had only one eligible child, the child was invited. When a household had several eligible children, a child in the age group with the fewest children already participating was selected. Towards the end of the survey, survey teams were allocated a specific number of subjects in each age group to recruit to avoid over-sampling. If the parents refused the participation of the selected child, the household was not included. This process was continued until the desired sample size was achieved in each of the 30 clusters.
Sample size
The sample size was calculated using EpiInfo Version 7 to estimate seroprevalence in each of four age groups (1–4, 5–9, 10–14 and 15–18 years old) with 95% confidence, a margin error of 5% and accounting for clustering with a design effect of 2. The expected national seroprevalence, based on Indonesian expert opinion and published regional data [14, 19, 27, 28], was 25% in the 1–4-year-old group, 45% in the 5–9-year-old group, 55% in the 10–14-year-old group and 65% in the 15–18-year-old group. To account for incomplete data, a 10% contingency was applied. The total sample size was 3,210 children, 660 from the 1–4-year-old group (22 per cluster), 870 from the 5–9-year-old group (29 per cluster), 870 from the 10–14-year-old group (29 per cluster) and 810 from the 15–18-year-old group (27 per cluster). In total, 107 children were enrolled in each cluster.
Enrolment
The study was presented to families during monthly neighbourhood association meetings. After household visits, eligible subjects were invited to the healthcare centre for enrolment and blood sampling if they were healthy, 1–18 years of age on inclusion day, and had lived in the location for at least 1 year. An informed consent form was signed by a parent or legal guardian, and by the subject if aged 13–18 years. Subjects aged 8–12 years provided signed assent.
A questionnaire was administered to collect information on demographics, knowledge of dengue symptoms and transmission, vector control practice, and medical history in the household.
Blood sampling and laboratory analysis
For each subject, 2mL of venous blood was drawn into plain vacutainer tubes. After centrifugation, serum aliquots were frozen at -20°C before refrigerated transport by courier to a central laboratory for analysis. Each specimen was tested for dengue IgG antibodies by ELISA using the commercial Panbio Dengue IgG Indirect ELISA kit (sensitivity = 96.3%; specificity = 91.4–100% according to manufacturer’s instructions; Panbio, Alere, Australia) [29]. Samples were considered positive for previous dengue infection according to the standard protocols of the manufacturer (Panbio units <9 is negative; 9–11 is equivocal; and >11 is positive).
Data analysis and statistics
All analyses were run using SAS 9.4.
Dengue antibody seroprevalence and associations between serologic status and socio-demographic and behavioural factors.
The statistical unit was the individual subject.
Seroprevalence and the 95% confidence interval (95% CI) were calculated taking account of the cluster effect. Univariate logistic regression was used to identify variables significantly associated with serologic status. As the data structure was hierarchical with subjects included in clusters, hierarchical logistic regression models were used to consider subject intra-cluster correlation. The clusters account for the random effect and the covariates were taken as fixed effects. As these analyses were considered exploratory, a level of significance (p-value) of <0.15 was applied at univariate level.
The multivariate hierarchical model was reduced by applying a backward descending selection of the non-significant variables at p-value >0.05.
The final model was:
Where
was the probability for a j subject from a i cluster to be seropositive, the βs were the fixed effect describing the subject variables associated with socio-demographic and behavioural factors, μ the cluster random effect and ɛ the error term.
Median age of conversion.
The median age of seroconversion was estimated by fitting a weighted linear regression model to age-specific seroprevalence data. Seroprevalence data were transformed into probits and age values were log transformed to fit the model [30, 31]. However, goodness of fit parameters were not respected. Therefore, a simple linear regression was used.
Force of infection.
Catalytic models use seroprevalence data as cumulative markers of past infections that result in life-long immunity from which force of primary infection estimates can be derived. [32, 33], Two force of infection models were developed to describe the rate of infection over the last 18 years and to examine its variability over time. The first model assumed a constant force of infection (model 1) and the second one assumed a force of infection that varied with age (model 2) [13].
The probability of a person living in the area being infected in one year, the force of infection, is estimated by [34]:
Where μ is the mean number of infections per year.
The variable force of infection model can be estimated by allowing a separate risk of infection for each age group, were pi is the mean number of infections per year for the ith age group and A is the age midpoint of the ith age group [34]. By fitting a binomial model with a complementary log-log link function and by using X = log(A) as an offset term, α = log(μ) can be estimated as an intercept parameter [34]. The probability of being infected for the ith group at midpoint age A is pi = 1- exp(-μi Ai), so that:
Results
Site selection and baseline demographics
From a total of 6,299 Indonesian subdistricts, 2,823 with urban population were identified, 2,756 of which had an urban population >1,000 and were thus used for sampling. A map of the 30 selected clusters is presented in Fig 1. From 30 October 2014 to 27 November 2014, a total of 3,210 subjects were enrolled in the study; 39 subjects (1.2%) were excluded due to at least one criteria of eligibility not being fulfilled and four subjects (0.1%) due to missing or incomplete data (demographic or serologic status result). A total of 3,194 subjects (98.7%) were included in the analyses (Fig 2); there were 107 subjects per site with the exception of four sites with 106 subjects, three sites with 105 subjects and one site with 101 subjects.
There were 672 subjects in the 1–4-year-old age group, 861 subjects in the 5–9-year-old age group, 886 in the 10–14-year-old age group and 775 in the 15–18-year-old age group. Among them, 47.8% were male and the mean age was 9.7 years.
Dengue antibody seroprevalence and association between serologic status and socio-demographic and behavioural factors
The age-specific seroprevalence ranged from 26.4% (95% CI: 15.8–37.1) in those aged 1-year-old to 95.3% (95% CI: 89.8–100) in the 18-year-old subjects (Fig 3). The median age at seroconversion was 4.8 years. The overall nationwide seroprevalence was 69.4%, with a minimum of 34.6% and a maximum of 87.9% observed per site, and the seroprevalence per age group was 33.8% in the 1–4year-old group, 65.4% in the 5–9-year-old group, 83.1% in the 10–14-year-old group and 89.0% in the 15–18-year-old group (Table 1).
In the final data set, the level of non-response (“no data”) varied from 0.4 to 14.0% (Table 1). Subjects were familiar with dengue disease, with 92% having heard about dengue and 91.4% able to cite at least one symptom. Control practices reported included use of repellent cream or mosquito spray (43.8%), elimination of mosquito breeding sites by covering water containers (59.0%) and eliminating stagnant water around the home (85.1%). Most subjects (75.3%) reported they had never been diagnosed with dengue.
Age and gender were associated with dengue serological status, with seroprevalences increasing with age (p<0.0001) and values of 71.1% (95% CI: 65.9–76.3) in females versus 67.4% (95% CI: 62.4–72.5) in males (p = 0.018) (Table 1). After univariate analysis, the type of household (p = 0.08), the level of education of the parents/guardians (p<0.0001), the number of persons living in the household (p<0.0001), knowledge about dengue symptoms (p = 0.14), sleeping under an untreated bed net (p = 0.10), the number of dengue cases identified since the subject was born (p<0.0001), and a previous clinical diagnosis of dengue for the subject (p<0.0001) were also associated with dengue serological status. In the multivariate model (Table 2), two variables remained associated with the dengue serologic status, the subject age group (1–4 vs 5–9 OR = 4.25; 1–4 vs. 10–14 OR = 12.60; and 1–4 vs 15–18 OR = 21.87; p<0.0001) and the number of cases diagnosed in the household since the subject was born (p = 0.0004).
Force of infection
The constant force of infection model was valid and estimated a force of primary infection of 13.1% per year in dengue-naïve children. As a result of the goodness of fit statistic being close to 0.05, a model of varying force of infection (age groups of one year) was run to examine the homogeneity of the force of primary infection estimates per age group. As suggested by the first model, there was no clear trend in changes in force of infection with age; the estimates were overlapping, ranging from 10.2% to 18.5% per year. The highest force of primary infection was observed in the 1-year-old age group (Table 3).
Discussion
This is the first dengue antibody seroprevalence study conducted in a representative population of urban dwelling Indonesian children. The findings benefit from a cluster sampling design with probability proportional to size method, and sensitive and specific dengue diagnostic assays performed in the same laboratory.
This study found that 69.4% of children had been previously infected with dengue virus, more than 80% of children aged 10 years or over, indicating that the disease burden is extremely high. A seroprevalence study conducted in 1995 in healthy children in Yogyakarta, Indonesia, using the plaque reduction neutralization test to determine previous exposure, reported the presence of neutralizing antibodies in 56.2% of 4–9-year-old children, ranging from 37.2% in 4-year-old subjects to 69.7% in those 9 years of age. These are slightly lower than the rates observed in our study (Fig 3 and Table 1) and may be reflective of increasing dengue endemicity in the intervening decades, or geographic variability [19]. Our results also show higher levels of dengue virus exposure than those reported in other dengue endemic countries such as Sri Lanka (Colombo, 2008, 52.0% in those <12 years of age, and median age of seroconversion of 4.7 years) [13, 35], and Vietnam (Binh Thuan, 2003, 65.7% in 7–13 year olds) [14]. This elevated dengue exposure risk was also observed during a 2011 dengue vaccine trial in 5 Asian countries, where baseline dengue seroprevalence was highest in Indonesian children [36].
Our constant force of infection model estimated a 13.1% annual rate of primary infection among 1–18-year-old children, while the variable model estimated a force of infection that varied from 10.2% to 18.5%. These estimates are similar to those reported in Sri Lanka in 2008 (14.1% in those aged <12 years) and Southern Vietnam in 2003 (11.7% in 7–13-year-old children) [13, 14]. Despite these similarities between Vietnam and Indonesia in terms of transmission dynamics, the reported incidence of disease in Vietnam is more than twice that in Indonesia. [37]. A number of hypotheses could explain this difference in findings: most likely, it is reflective of Indonesia’s specific case definition for reported dengue disease (only DHF is reported), but underlying virological, genetic or epidemiological differences could play a role. From the constant force of primary infection model, it can be assumed that the average rate of primary infection was not highly variable over the past 18 years. Additional analysis may be needed to better understand infection risk over time. The recently observed increase in age distribution of reported cases may have been driven by more variable virologic, demographic, reporting or other determinants of disease [10]. A similar phenomenon was illustrated by a study conducted in Thailand showing that the upward shift in dengue case age was associated with demographic changes [38].
It can be assumed that dengue awareness, through social mobilization and education campaigns, begun in the 1970s, and the increasing public health importance associated with high media coverage, has steadily increased [39]. Knowledge of dengue transmission and symptoms was high within the study subjects; 92% of households had heard about dengue before our study and were able to cite at least one of the disease symptoms, and more than 80% knew that dengue virus is transmitted by diurnal mosquito bites. In term of exposure, household practices were focused on destroying mosquito breeding sites rather than personal protection. The level of exposure to the virus, however, is strong evidence that these reported behaviours are inadequate to protect against infection and additional prevention and control measures are urgently required.
In the multivariate model, only subject age group and the number of dengue cases that occurred in the household were associated with seropositive status. Some of the parameters significantly associated with dengue seropositivity in univariate models were also implicated in other dengue studies conducted in Latin America and Asia. For example, parental level of education and dengue illness history in the household have been associated with dengue seropositivity [17]. Other parameters, such as household size, exhibit an association inverse to that previously reported in the literature [40]. This is most likely explained by confounding effects from known risk factors such as age or unknown, socio-demographic drivers of exposure risk. The lack of significant associations between socio-demographic and behavioural factors with serological status provides evidence that essentially everyone is at risk of infection; that knowledge of prevention and control at the individual/household level is not protective against infection; and that additional measures to prevent transmission are required. The retrospective nature of our questionnaire limits the robustness of our results; recall bias may have been an issue.
A recent expansion in dengue virus transmission from urban to peri-urban and rural areas has been described [15] and the identification of provinces or areas of high transmission risk is a focus of prevention and control planning. This study showed a high level of exposure across urban Indonesia and, while we excluded rural areas from this study for operational reasons, it is likely that nearby peri-urban populations may have experienced similar high levels of exposure [40]. Another possible limitation is that cross-reaction between flaviviruses has been documented and the risk of false positives cannot be excluded. We consider this risk as low, because reports of other viruses such as Japanese encephalitis and Zika, in Indonesia, are rare. This study was not designed to make national-level infection or disease burden estimates but the observation that 13.1% of children suffer a primary infection per year translates into many millions of infections per year. Adults are presumably infected with a similar frequency. A proportion of these infections will be secondary, predisposing to symptomatic and severe disease. While a modelling approach would be required to quantify this burden, these data are strongly suggestive that dengue infections result in a significant burden of symptomatic and severe disease in urban Indonesia.
Supporting information
S1 Appendix. List of the 30 clusters selected in Indonesia.
https://doi.org/10.1371/journal.pntd.0005621.s002
(DOCX)
Acknowledgments
The authors would like to thank Annick Moureau from Sanofi Pasteur for her technical advice on statistical methods, Benediktus Yohan and Rahma F. Hayati from Eijkman Institute for Molecular Biology for their technical assistance in laboratory analysis, the health centre employees who conducted the study in the field, the subjects and their family who agreed to take part to this study, the employee of the Indonesian Sub-Directorate of Statistical Services and Promotion Statistics (BPS) for providing the data for the sampling of urban areas. Juliette Gray of inScience Communications, Springer Healthcare provided manuscript editing assistance and provided support for the submission.
Author Contributions
- Conceptualization: AFT JN AP JMS AB.
- Formal analysis: AFT.
- Funding acquisition: AFT JN JMS AB.
- Investigation: AP SRH HIS RMK RS SS HG RTS BEM.
- Methodology: AFT JN JMS.
- Project administration: AP SRH HIS MRK RS SS HG RTS BEM.
- Supervision: SRH AP.
- Validation: AFT JN AP JMS AB AP SRH HIS MRK RS SS HG RTS BEM.
- Visualization: AFT.
- Writing – original draft: AFT JN JMS.
- Writing – review & editing: AFT JN AP JMS AB AP SRH HIS MRK RS SS HG RTS BEM.
References
- 1.
WHO, editor. Global Strategy For Dengue Prevention And Control, 2012–2020. Geneva: WHO Press; 2012.
- 2. Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013.
- 3. Stanaway JD, Shepard DS, Undurraga EA, Halasa YA, Coffeng LE, Brady OJ, et al. The global burden of dengue: an analysis from the Global Burden of Disease Study 2013. Lancet Infect Dis. 2016. Epub 2016/02/15. pmid:26874619.
- 4.
WHO. Dengue Haemorrhagic Fever: Diagnosis, treatment, prevention and control: World Health Organization; 1997.
- 5. Guzman MG, Halstead SB, Artsob H, Buchy P, Farrar J, Gubler DJ, et al. Dengue: a continuing global threat. Nat Rev Microbiol. 2010;8(12 Suppl):S7–16. pmid:21079655.
- 6. Horstick O, Runge-Ranzinger S, Nathan MB, Kroeger A. Dengue vector-control services: how do they work? A systematic literature review and country case studies. Transactions of the Royal Society of Tropical Medicine and Hygiene. 2010;104(6):379–86. Epub 2010/04/20. pmid:20400169.
- 7. Schwartz LM, Halloran ME, Durbin AP, Longini IM. The dengue vaccine pipeline: Implications for the future of dengue control. Vaccine. 2015.
- 8. Capeding MR, Tran NH, Hadinegoro SR, Ismail HI, Chotpitayasunondh T, Chua MN, et al. Clinical efficacy and safety of a novel tetravalent dengue vaccine in healthy children in Asia: a phase 3, randomised, observer-masked, placebo-controlled trial. Lancet. 2014;384(9951):1358–65. Epub 2014/07/16. pmid:25018116.
- 9.
Indonesia S. Population of Indonesia by Village—Result of the 2010 population census 2010. http://bps.go.id/eng/download_file/Population_of_Indonesia_by_Village_2010.pdf.
- 10. Karyanti M, Uiterwaal CSP, Kusriastuti R, Hadinegoro S, Rovers M, Heesterbeek H, et al. The changing incidence of Dengue Haemorrhagic Fever in Indonesia: a 45-year registry-based analysis. BMC infectious diseases. 2014;14(1):412. pmid:25064368
- 11.
SEARO. Reported cases of dengue and deaths from the SEA countries 2000–2011 2013 [cited 2013 November 2013]. http://www.searo.who.int/entity/vector_borne_tropical_diseases/data/ReportedCasesDeaths.pdf.
- 12.
Indonesia MoH. Situasi Demam Berdarah Dengue Di Indonesia. [Situation of dengue fever in Indonesia]. April 2014 ed2014.
- 13. Tam CC, Tissera H, de Silva AM, De Silva AD, Margolis HS, Amarasinge A. Estimates of dengue force of infection in children in Colombo, Sri Lanka. PLoS Negl Trop Dis. 2013;7(6):e2259. pmid:23755315.
- 14. Thai KT, Binh TQ, Giao PT, Phuong HL, Hung le Q, Van Nam N, et al. Seroprevalence of dengue antibodies, annual incidence and risk factors among children in southern Vietnam. Tropical medicine & international health: TM & IH. 2005;10(4):379–86. pmid:15807802.
- 15. Muhammad Azami NA, Salleh SA, Neoh HM, Syed Zakaria SZ, Jamal R. Dengue epidemic in Malaysia: Not a predominantly urban disease anymore. BMC Res Notes. 2011;4:216. Epub 2011/07/01. pmid:21714858.
- 16. Wilder-Smith A, Foo W, Earnest A, Sremulanathan S, Paton NI. Seroepidemiology of dengue in the adult population of Singapore. Tropical Medicine & International Health. 2004;9(2):305–8.
- 17. Siqueira JB, Martelli CM, Maciel IJ, Oliveira RM, Ribeiro MG, Amorim FP, et al. Household survey of dengue infection in central Brazil: spatial point pattern analysis and risk factors assessment. The American journal of tropical medicine and hygiene. 2004;71(5):646–51. pmid:15569799
- 18. Chairulfatah A, Setiabudi D, Agoes R, van Sprundel M, Colebunders R. Hospital based clinical surveillance for dengue haemorrhagic fever in Bandung, Indonesia 1994–1995. Acta tropica. 2001;80(2):111–5. pmid:11600087
- 19. Graham R, Juffrie M, Tan R, Hayes C, Laksono I, Ma'roef C, et al. A prospective seroepidemiologic study on dengue in children four to nine years of age in Yogyakarta, Indonesia I. studies in 1995–1996. The American journal of tropical medicine and hygiene. 1999;61(3):412–9. pmid:10497982
- 20. Porter KR, Beckett CG, Kosasih H, Tan RI, Alisjahbana B, RUDIMAN PIF, et al. Epidemiology of dengue and dengue hemorrhagic fever in a cohort of adults living in Bandung, West Java, Indonesia. The American journal of tropical medicine and hygiene. 2005;72(1):60–6. pmid:15728868
- 21.
Indonesia B-S. Statistical Yearbook of Indonesia 2015. In: Publication S-doSCa, editor. Indonesia: BPS-Statistics Indonesia; 2015. p. 670.
- 22.
Indonesia B-S. Percentage of Urban Population by Province, 2010–2035 2015 [cited 2015 01/09/2015]. http://www.bps.go.id/linkTabelStatis/view/id/1276.
- 23.
Mulyana W. Rural-Urban Linkages: Indonesia Case Study. In: Urban and Regional Development Institute (URDI) J, Indonesia, editor. Santiago, Chile: RIMISP; 2014.
- 24. Henderson RH, Sundaresan T. Cluster sampling to assess immunization coverage: a review of experience with a simplified sampling method. Bulletin of the World Health Organization. 1982;60(2):253. pmid:6980735
- 25. Turner AG, Magnani RJ, Shuaib M. A not quite as quick but much cleaner alternative to the Expanded Programme on Immunization (EPI) Cluster Survey design. International journal of epidemiology. 1996;25(1):198–203. pmid:8666490
- 26.
Center for International Earth Science Information Network—CIESIN—Columbia University, International Food Policy Research Institute—IFPRI, The World Bank, Centro Internacional de Agricultura Tropical—CIAT. Global Rural-Urban Mapping Project, Version 1 (GRUMPv1): Urban Extents Grid. Palisades, NY: NASA Socioeconomic Data and Applications Center (SEDAC); 2011.
- 27. Malavige G, Fernando S, Aaskov J, Sivayogan S, Dissanayaka T, Peelawattage M K, et al. Seroprevalence of Anti-dengue Virus Antibodies in Children in Colombo District, Sri Lanka. Dengue Bulletin. 2006;30(Dec.).
- 28. Pengsaa K, Limkittikul K, Luxemburger C, Yoksan S, Chambonneau L, Ariyasriwatana C, et al. Age-specific prevalence of dengue antibodies in Bangkok infants and children. The Pediatric infectious disease journal. 2008;27(5):461–3. Epub 2008/03/25. pmid:18360303.
- 29.
matters AKw. Panbio® Dengue IgG Indirect ELISA 2010–2015. http://www.alere.com/ww/en/product-details/panbio-dengue-igg-indirect-elisa.html.
- 30. Ferreira MU, Hiramoto RM, Aureliano DP, da Silva-Nunes M, da Silva NS, Malafronte RS, et al. A Community-based Survey of Human Toxoplasmosis in Rural Amazonia: Seroprevalence, Seroconversion Rate, and Associated Risk Factors. The American journal of tropical medicine and hygiene. 2009;81(1):171–6. pmid:19556584
- 31. Ferreira MU, Camargo LMA, Ferreira CS. The median age of primary malaria infection in an Amazonian community: Probit analysis of cross-sectional data. Acta Tropica. 1994;58(1):73–8. http://dx.doi.org/10.1016/0001-706X(94)90122-8. pmid:7863855
- 32. Ferrari MJ, Djibo A, Grais RF, Grenfell BT, Bjornstad ON. Episodic outbreaks bias estimates of age-specific force of infection: a corrected method using measles as an example. Epidemiology and infection. 2010;138(1):108–16. pmid:19538818.
- 33. Whitaker HJ, Farrington CP. Estimation of infectious disease parameters from serological survey data: the impact of regular epidemics. Statistics in medicine. 2004;23(15):2429–43. Epub 2004/07/27. pmid:15273957.
- 34.
Inc. SI. SAS/STAT® 13.1 User’s Guide. Inc. SI, editor. Cary, NC, USA: SAS Institute Inc.; 2013.
- 35. Tissera HA, De Silva AD, Abeysinghe MRN, de Silva AM, Palihawadana P, Gunasena S, et al. Dengue Surveillance in Colombo, Sri Lanka: Baseline seroprevalence among children. Procedia in Vaccinology. 2010;2(1):109–12. http://dx.doi.org/10.1016/j.provac.2010.03.020.
- 36. L’Azou M, Moureau A, Sarti E, Nealon J, Zambrano B, Wartel TA, et al. Symptomatic Dengue in Children in 10 Asian and Latin American Countries. New England Journal of Medicine. 2016;374(12):1155–66. pmid:27007959.
- 37.
SEARO. Comprehensive Guidelines for Prevention and Control of Dengue and Dengue Haemorrhagic Fever. WHO SEARO, 2011.
- 38. Rodriguez-Barraquer I, Buathong R, Iamsirithaworn S, Nisalak A, Lessler J, Jarman RG, et al. Revisiting Rayong: shifting seroprofiles of dengue in Thailand and their implications for transmission and control. American journal of epidemiology. 2014;179(3):353–60. Epub 2013/11/08. pmid:24197388.
- 39. Kusriastuti R, Sutomo S. Evolution of dengue prevention and control programme in Indonesia. Dengue bulletin. 2005;29:1.
- 40. Tipayamongkholgul M, Lisakulruk S. Socio-geographical factors in vulnerability to dengue in Thai villages: a spatial regression analysis. Geospatial health. 2011;5(2):191–8. pmid:21590669