- Research article
- Open Access
- Open Peer Review
Chikungunya virus infection in Indonesia: a systematic review and evolutionary analysis
BMC Infectious Diseases volume 19, Article number: 243 (2019)
Despite the high number of chikungunya cases in Indonesia in recent years, comprehensive epidemiological data are lacking. The systematic review was undertaken to provide data on incidence, the seroprevalence of anti-Chikungunya virus (CHIKV) IgM and IgG antibodies, mortality, the genotypes of circulating CHIKV and travel-related cases of chikungunya in the country. In addition, a phylogenetic and evolutionary analysis of Indonesian CHIKV was conducted.
A systematic review was conducted to identify eligible studies from EMBASE, MEDLINE, PubMed and Web of Science as of October 16th 2017. Studies describing the incidence, seroprevalence of IgM and IgG, mortality, genotypes and travel-associated chikungunya were systematically reviewed. The maximum likelihood phylogenetic and evolutionary rate was estimated using Randomized Axelerated Maximum Likelihood (RAxML), and the Bayesian Markov chain Monte Carlo (MCMC) method identified the Time to Most Recent Common Ancestors (TMRCA) of Indonesian CHIKV. The systematic review was registered in the PROSPERO database (CRD42017078205).
Chikungunya incidence ranged between 0.16-36.2 cases per 100,000 person-year. Overall, the median seroprevalence of anti-CHIKV IgM antibodies in both outbreak and non-outbreak scenarios was 13.3% (17.7 and 7.3% for outbreak and non-outbreak events, respectively). The median seroprevalence of IgG antibodies in both outbreak and non-outbreak settings was 18.5% (range 0.0–73.1%). There were 130 Indonesian CHIKV sequences available, of which 120 (92.3%) were of the Asian genotype and 10 (7.7%) belonged to the East/Central/South African (ECSA) genotype. The ECSA genotype was first isolated in Indonesia in 2008 and was continually sampled until 2011. All ECSA viruses sampled in Indonesia appear to be closely related to viruses that caused massive outbreaks in Southeast Asia countries during the same period. Massive nationwide chikungunya outbreaks in Indonesia were reported during 2009–2010 with a total of 137,655 cases. Our spatio-temporal, phylogenetic and evolutionary data suggest that these outbreaks were likely associated with the introduction of the ECSA genotype of CHIKV to Indonesia.
Although no deaths have been recorded, the seroprevalence of anti-CHIKV IgM and IgG in the Indonesian population have been relatively high in recent years following re-emergence in early 2001. There is sufficient evidence to suggest that the introduction of ECSA into Indonesia was likely associated with massive chikungunya outbreaks during 2009–2010.
Arboviruses (arthropod-borne viruses) are a group of viruses that exist in a transmission cycle between blood-feeding arthropod vectors and amplifying, vertebrate hosts. With most arboviruses, human involvement in this transmission cycle is incidental . In terms of public health significance, the mosquito is the most important vector of arbovirus transmission. It is estimated that approximately 3.9 billion people, living in more than 120 different countries, are at risk of becoming infected with any of the three major arboviruses: Chikungunya virus (CHIKV), Dengue virus (DENV) and Zika virus (ZIKV) . Chikungunya virus, which is primarily transmitted by Aedes aegypti and Ae. albopictus mosquitoes, is a positive-sense single-stranded RNA virus and a member of the family Togaviridae [3, 4]. The first well-characterised chikungunya outbreak was reported in Southern province, Tanganyika territory of Tanzania in 1952 [5, 6]. Sporadic chikungunya outbreaks were subsequently identified in parts of Africa and Asia during the 1950s and 1960s, followed by an apparent re-emergence in the 2000s . Since 2005, large-scale outbreaks of chikungunya sweeping across south-western Indian Ocean and Southeast Asia [8,9,10,11,12,13,14,15,16,17,18,19,20]. In La Réunion, the outbreak affected about a third of the population [9, 21] and in India, the viruses infected more than 1.3 million persons during 2005–2006 . In Sri Lanka, the viruses infected more than 100,000 people  and CHIKV subsequently spread to Southeast Asia including Indonesia.
The illnesses caused by CHIKV and DENV are clinically indistinguishable and the accurate diagnosis of these infections on clinical grounds alone is usually problematic [22, 23]. Although prior literature has stated that a higher proportion of people infected with CHIKV are symptomatic than those infected with DENV , a recent systematic review revealed that asymptomatic chikungunya had a very high variability in percentages ranging from 3.2% in La Réunion (2005–2006) to 82.1% in the Philippines (2012–2013) . The highest percentage of asymptomatic chikungunya have been recorded with the Asian genotype in the Philippines with 82.1% . The common symptoms of chikungunya include rash, high fever, severe joint and muscle pain, headache and photophobia [7, 27]. Severe symptoms, involving vital organs, may develop during CHIKV infection such as encephalitis [28, 29], encephalopathy [29,30,31], optic neuropathy [29, 32], neuroretinitis , myelopathy and myelitis , Guillain-Barré syndrome [29, 32], myocarditis , hepatitis , acute interstitial nephritis , severe sepsis , septic shock  and multi-organ failure [31,32,33, 36, 37]. In rare cases, infections may be fatal [28, 33,34,35, 37]. Perinatal CHIKV infection can cause sequelae such as microcephaly and cerebral palsy . In adults, persisting arthralgia/arthritis, alopecia and depression are the most commonly recorded sequelae [39,40,41,42]. A meta-analysis found that approximately 25% of chikungunya cases develop chronic inflammatory rheumatism and 14% develop chronic arthritis , creating a major burden on society in terms of morbidity and economic productivity [41, 42, 44, 45].
CHIKV has an approximately 12 kb genome that encodes four non-structural proteins (NSP1–4) and five structural proteins (C, E3, E2, 6 K, and E1) . Genetic analysis has identified three genotypes of CHIKV: the West African, East/Central/South African (ECSA), and Asian genotype . The ECSA genotype consists of three lineages: Central African, East/South African, and the Indian Ocean lineage. These genotypes are spreading sporadically worldwide, with ECSA and Asian genotypes being the predominately isolated .
In Indonesia, based on official documents from the Ministry of Health (MoH), chikungunya cases were reported for the first time in Samarinda (Kalimantan island) in 1973 . However, there is evidence to suggest that CHIKV infections have occurred in Indonesia prior to 1973. Serum samples collected between 1969 and 1972 demonstrated significant titres of anti-CHIKV antibodies when tested using haemagglutination inhibition assay (HI) and the plaque reduction neutralization test (PRNT) in most of the Indonesian archipelago, except Java [49, 50]. In addition, evidence from historical reports suggest that the first circulation of CHIKV was back in 1779 in Jakarta, when this infection called as kidinga pepo [51, 52]. This is widely acknowledged by experts in the field of arboviruses as the first report of chikungunya in Indonesia, although it is impossible to demonstrate by molecular clock analysis . The first virologically confirmed chikungunya outbreak was reported in June 1982 in Jambi province of Sumatra island, followed by multiple outbreaks between 1983 and 1984 . Chikungunya cases were not recorded in Indonesia for approximately 20 years, before the infection appeared to re-emerge and cause multiple outbreaks in South Sumatera, Aceh and West Java in early 2001 . Since then multiple outbreaks have been reported [48, 55,56,57,58,59,60,61,62,63,64,65]. Despite the high number of chikungunya cases in Indonesia in recent years, comprehensive epidemiological data are lacking.
The overarching aim of the study was to provide a comprehensive overview of chikungunya epidemiology in Indonesia including an evolutionary analysis of CHIKV that have circulated in the country between 1983 and 2016. Through a systematic review of available literature and available sequence data on GenBank, the primary objectives of this study were: a) to estimate the incidence of chikungunya in Indonesia; b). to characterize the seroprevalence of anti-CHIKV IgM and IgG antibodies in Indonesia; c) to describe the mortality of chikungunya in Indonesia and; d) to genetically characterize CHIKV circulating in Indonesia that were isolated locally or from travelers returning from Indonesia. The secondary objective was to provide data of travel-associated chikungunya originating from Indonesia.
Protocol and registration
The systematic review was conducted as recommended by the Preferred Reporting Items for Systematic Reviews and Meta-analyses (PRISMA) guidelines . The protocol of the systematic review was registered at PROSPERO, an international database of prospectively registered systematic reviews at the University of York (CRD42017078205).
Role of the funding source
The funders of the study had no role in study design, data collection, data analysis, data interpretation, writing of the report and decision to publish.
In this study, all full articles or abstracts published in any year and language were included while editorials, reviews, commentaries, and qualitative studies were excluded for all outcomes except for mortality. This is because study reporting chikungunya mortality is limited. Studies in any setting (outbreak or non-outbreak) with any study design (cross-sectional or prospective) were considered eligible.
Any publications that reported the number of acute chikungunya infections or seroconversions over any time interval were eligible, as proposed previously . Studies reporting chikungunya outbreaks with an attack rate were also considerate eligible if information about the population at risk was given.
Seroprevalence of anti-CHIKV IgM and IgG antibodies
Any publications that provided primary data on the seroprevalence of anti-CHIKV IgM and IgG antibodies in Indonesia were considered eligible. Eligible studies included those reporting anti-DENV IgM and/or IgG seroprevalence, measured by the presence of IgM and IgG antibodies, among individuals with suspected acute chikungunya infection, with suspected other arbovirus infection, with undifferentiated acute febrile illness or who were healthy. Studies reporting anti-DENV IgM and/or IgG during outbreak and non-outbreak settings were considered eligible.
Studies reporting chikungunya-associated deaths either during an outbreak or in a non-outbreak setting were considered eligible. For mortality, all type of articles (full papers, case studies, letters, editorials, and commentaries) with adequate information were considered eligible. Studies reporting chikungunya-associated death among travellers returning from Indonesia were excluded.
Chikungunya virus genotype
All CHIKV sequences originating from Indonesia, either isolated locally or in travellers returning from Indonesia, were included. CHIKV sequences originating from Indonesian mosquito samples were also included. Single and partial gene sequences, of either non-structural or structural genes, as well as complete genome sequences were included. Studies reporting CHIKV genotype information, that failed to deposit sequence data into GenBank, were excluded.
All chikungunya cases reported in travellers returning from Indonesia were included. Cases were confirmed by anti-CHIKV IgM, reverse transcription polymerase chain reaction (RT-PCR) or viral isolation. Cases without laboratory confirmation were excluded.
Information sources and search strategy
A systematic search was conducted using four bibliographical databases (EMBASE, MEDLINE, PubMed and Web of Science as of October 16th 2017) to identify potential articles using Medical Subject Headings (MeSH) encompassing the terms “Chikungunya” AND “Indonesia”. No limit was set for the publication year or language. In addition, grey literature and official reports from the MoH, published in either Indonesian or English, were included. Reference lists of the identified articles and references of deposited sequences in GenBank were also searched manually to find additional potential studies. Some authors were contacted to obtain or clarify information.
All titles and abstracts of identified articles were imported into the local library of EndNote X7 (Thompson Reuters, Philadelphia, PA, USA) and duplicate records between multiple databases were removed. Retrieved articles were initially screened based on title and abstract to identify possible eligible studies. The full texts of potential eligible articles were then reviewed. The screening and review process was conducted by two authors (HH and AM). After reviewing full texts, the eligibility of each study was decided. If there was disagreement between two authors, a third made the final decision.
Data collection process and synthesis
The number of new chikungunya cases and population at risk or overall attack rate during an outbreak were recorded. Location, date, setting of study, diagnosis method, number of confirmed or suspected cases and other related information were also extracted. The incidence of chikungunya was expressed in terms of person-year.
Besides systematically reviewing the available literature, we also accessed and extracted data from the National Disease Surveillance run by the Directorate General of Disease Prevention and Control of the Indonesian MoH. Chikungunya is a notifiable disease in Indonesia. The chikungunya case definition and case ascertainment used in this surveillance system have previously been published at the National Guideline for Prevention and Control of Chikungunya from MoH . This case definition follow the World Health Organization (WHO) criteria . In brief, chikungunya cases were classified into three categories: a) possible case, diagnosed based on clinical criteria alone as acute onset of fever > 38.5°C and severe arthralgia/arthritis not explained by other medical conditions; b) probable case, diagnosed based on the clinical criteria as mentioned and epidemiological criteria (residing or having visited epidemic areas) and; c) confirmed case, diagnosed based on laboratory criteria which show a positive result for virus isolation, RT-PCR, IgM antibodies or a four-fold increase in IgG antibodies . This surveillance system includes all types of chikungunya cases. The number of chikungunya cases were retrieved from the database. This primary dataset was analysed separately from other studies from systematic review. The annual chikungunya incidence rate was determined by dividing the number of new chikungunya cases, by the size of the population at risk based on MoH data. The incidence rate was expressed as per 100,000 person-years. To obtain a more comprehensive distribution of the chikungunya cases in Indonesia, geographical maps of provincial incidence rates from 2008 to 2016 were created using ArcGIS software . The provincial incidence rate was expressed as the number of chikungunya cases per 100,000 civilians.
Seroprevalence of anti-CHIKV IgM and IgG antibodies
Because each study had different settings, designs and tested sera (collected from febrile or non-febrile subjects), which may potentially affect the percentage of prevalence, the seroprevalence in this study was stratified. The seroprevalence of anti-CHIKV IgM and IgG antibodies were divided into: a) setting of study (outbreak and non-outbreak) and; b) status of tested sera (collected from healthy individuals and patients with undifferentiated acute febrile illness or suspected arbovirus infection). Therefore, the seroprevalence of anti-CHIKV IgM and IgG antibodies were stratified into four categories: a) among febrile patients during a chikungunya outbreak or post outbreak; b) among non-febrile patients during a chikungunya outbreak or post outbreak; c) among febrile patients not in an outbreak; and d) among non-febrile patients not in an outbreak. If one reference consisted of two or more studies with different settings or statuses of participants, then the data were divided into the relevant categories. The median value and the ranges of the seroprevalence of IgM and IgG antibodies were calculated for each setting.
For mortality, the number of deaths associated with CHIKV infections either during an outbreak or in a non-outbreak setting was included. Data such as location, date, setting of study, sample size, diagnosis method, number of cases and other related information were extracted. The case fatality rate was calculated as the number of deaths associated with chikungunya, confirmed by doctor resume on death certificate during passive surveillance or confirmation by outbreak investigation staff run by MoH, divided by the number of cases, expressed as a percentage (%). All deaths associated with chikungunya are reported continuously in the National Disease Surveillance managed by MoH.
Chikungunya virus genotype
To determine the CHIKV genotype, CHIKV sequences derived from locals and travellers returning from Indonesia were included. In addition, a search was conducted in GenBank to confirm the findings from systematic review and to obtain additional sequence data. To avoid missing sequences from our search on GenBank especially from travel-related cases, for example unlabelled origin of sequences, intensive manual search was conducted on references from the systematic review. Then manual searches were conducted in GenBank. Authors were contacted if there was disagreement between the information in the article and with the GenBank data. To provide more comprehensive data of circulating genotype over time, data from a study that consisted of genotype data over multiple years was divided and reported per year. Sequences were extracted and imported to Geneious v.10.1.3  for further phylogenetic analysis.
To fully assess epidemiology data of chikungunya in Indonesia, all confirmed CHIKV infections (as defined by serology or RT-PCR) reported from other countries that are known to have originated from Indonesia were also included. Essential data such as date, country, and number of cases, diagnosis method, and the genotype of CHIKV were collected. Data from a study that consisted of genotype data over multiple years were also divided and reported per year, if possible.
The primary outcomes of the systematic review, together with MoH database analysis, were: a) the incidence of chikungunya in Indonesia; b) the seroprevalence of anti-CHIKV IgM and IgG antibodies in Indonesia; c) chikungunya-related mortality in Indonesia and; d) a description of the CHIKV genotypes circulating in Indonesia. The secondary outcome was chikungunya cases exported from Indonesia to other countries.
Risk of bias assessment
The quality of eligible studies of incidence and seroprevalence IgM and IgG was evaluated by evaluating the risk of bias and the precision of the reported measures in accordance with a modified tool from previous study . The risk of bias was evaluated by critically appraising diagnostic methods, sampling, response rate, and study setting. Studies were considered to have a low, high, or unclear risk of bias based on those four domains. Risk of bias was considered low if: a) seroprevalence of CHIKV-specific antibodies was measured using a standard method such as enzyme-linked immunosorbent assay (ELISA), HI, and PRNT-based assays and chikungunya was diagnosed with laboratory confirmation such as IgM, RT-PCR or viral isolation; b) samples was selected using probability-based strategy - for incidence studies that enrolled individuals presenting to a health facility with acute infection this domain was not assessed; c) response rate was ≥80% and; d) the study setting was clearly stated (during outbreak or not). Studies were considered to have high precision if the number of individuals recruited or tested was ≥100 . For genotyping studies, the quality of study was judged based on: a) whether the origin of sequences is clearly mentioned (human or mosquito) and; b) whether the location and date of isolation are clearly provided. The quality of studies about exported chikungunya was assessed based on two criteria: a) clear information of travel history and; b) chikungunya cases were diagnosed using standard methods such as IgM, RT-PCR or viral isolation. If at least one criterion was missing in studies on genotyping or exported chikungunya, then the study was considered to have a high risk of bias. If a reference consists more than one studies, the risk of bias was assessed for each study.
Phylogenetic and evolutionary analysis
To provide solid data on CHIKV circulating in Indonesia, phylogenetic and evolutionary analyses were conducted for all Indonesian CHIKV sequences. To achieve this, Indonesian CHIKV sequences were downloaded from the GenBank database using the search keywords: “Chikungunya” AND “Indonesia” (last searched: October 20th, 2017). For comparison, all CHIKV sequences from other countries were also downloaded and those full genome sequences were trimmed leaving only the E1 gene. All Indonesian and comparison sequences were aligned using a Multiple Alignment using Fast Fourier Transform (MAFFT) v.7.309  as implemented in Geneious v.10.1.3 . To provide the relationship between Indonesian viruses and viruses isolated from other parts of world, an estimation of the maximum likelihood-based phylogenetic tree was performed using Randomized Axelerated Maximum Likelihood (RAxML) v.7.2.8 [73, 74]. General Time Reversible with gamma substitution model (GTR + Γ) and a rapid bootstrap procedure were employed with 1589 comparison sequences.
To provide evolutionary information of Indonesian viruses, their evolutionary rate and the Time to Most Recent Common Ancestors (TMRCA) was estimated using the Bayesian Markov chain Monte Carlo (MCMC) method as implemented in Bayesian Evolutionary Analysis Sampling Trees (BEAST) v.2.4.6 . The input file for BEAST was prepared using Bayesian Evolutionary Analysis Utility (BEAUti) and runs were performed using General Time Reversible with four gamma category and invariant sites (GTR + Γ4 + I). The runs used the relaxed lognormal molecular clock with an initial estimated evolutionary rate of 4.33 × 10− 4 substitutions per site per year  and with a Coalescent Bayesian skyline model. A hundred million chains were run and sampled every thousand cycles with 53 representative reference sequences. The output of MCMC runs was analyzed using Tracer v.1.6 to assess effective sampling size. The Maximum Clade Credibility (MCC) of phylogenetic trees was created using TreeAnnotator v.2.4.6 with 10% burn-in and visualized using FigTree v.1.4.3.
Study eligibility results
The searches yielded 182 references of which 84 were excluded as duplicates. After conducting a title and abstract screen of 98 references, an additional 55 references were excluded (Fig. 1). The full text references were retrieved for 43 titles; 19 had been excluded for one of four reasons: a) the full-text did not include data of any outcome of interest (n = 9); b) full-texts were not available and the abstract did not include adequate data for any outcome of interest (n = 8); c) duplicated data in another study (n = 1); and d) did not meet eligibility criteria (n = 1). This screen resulted in 24 eligible studies [13, 49, 50, 54, 55, 77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95]. An additional 16 references were identified from reference lists, reviews and GenBank references [20, 23, 76, 96,97,98,99,100,101,102,103,104,105,106,107,108] in which some of the references reported more than one of the outcomes of interest. A total of 40 references were included in the systematic review [13, 20, 23, 49, 50, 54, 55, 76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108]. Of these references, two studies assessed the incidence of chikungunya [94, 95], five studies assessed the seroprevalence of anti-CHIKV IgM antibodies [54, 55, 77, 78, 95] and seven assessed the seroprevalence of anti-CHIKV IgG antibodies [49, 50, 54, 55, 77,78,79]. Five studies reported mortality rates of chikungunya [13, 54, 55, 78, 80]. Genotype analyses of CHIKV circulating in Indonesia were reported in sixteen studies [13, 20, 23, 76, 78, 81,82,83,84,85, 94, 98, 100, 102, 103, 108], and exported chikungunya or travel-related chikungunya were reported in 24 studies [20, 23, 82, 84, 86,87,88,89,90,91,92,93, 96,97,98,99,100,101,102,103,104,105,106,107]. Detailed characteristics of each study is provided in Additional file 1. In addition, eleven official reports from the MoH of Indonesia, all reporting the incidence and mortality of chikungunya, were identified [48, 56,57,58,59,60,61,62,63,64,65].
Precision and risk of bias assessment
Quality assessment for each study included on chikungunya incidence and seroprevalence anti-CHIKV IgM and IgG is provided in Table 1. Both chikungunya incidence studies, four out five (75%) and five out seven (71.4%) of references of seroprevalence anti-CHIKV IgM and IgG, respectively, contained high precision as defined by a sample size of > 100 subjects. For the risk of bias assessment, response rate was not reported in all studies and all studies used valid essay. Among seroprevalence studies, 16% (2/12) studies used probability sampling (i.e. low risk of bias) and all studies had clear setting of the study. There was no variability of risk of bias among studies on CHIKV genotype in Indonesia and exported chikungunya cases. All studies on CHIKV genotypes provided clear information of sequences (origin of the sequences, location and date of isolation) and all exported chikungunya studies had robust methods in diagnosing chikungunya and clear information of travel history (i.e. low risk of bias). See Additional file 1 for detailed risk of bias assessment for each study.
Incidence of chikungunya
Two studies assessed the incidence of chikungunya in Indonesia [94, 95]. One study took place in Bandung (West Java province) between 2000 and 2004 and between 2006 and 2008 . This study found that the incidence of chikungunya was 10.1/1000 persons-year . Another study in three localities in Indonesia found similar incidence, 8.8/1000 persons-year during 2010–2011 . In addition, eleven annual reports from the MoH of Indonesia were identified [48, 56,57,58,59,60,61,62,63,64,65], providing the incidence rate of chikungunya between 2004 and 2015. This data revealed that the lowest chikungunya incidence was in 2005 with 0.16/100,000 person-year  while the highest incidence rate was recorded in 2009 with 36.2 cases per 100,000 persons-year . In 2009, more than 83 thousand cases were reported and CHIKV been reported circulate in seventeen out of 34 provinces in Indonesia .
To provide a more comprehensive picture of incidence rate of chikungunya, we accessed and extracted the National Disease Surveillance database at the MoH. The database covered chikungunya incidence between 2001 and 2016. In this 16-year period, the highest number of chikungunya cases was recorded in 2009 with 83,756 cases, equivalent to an incidence rate of 36.2 cases per 100,000 person-year (Fig. 2). The map of provinvial indicence rates indicated that chikungunya cases did not distribute equally among Indonesian regions (Fig. 3). The highest incidence of chikungunya occurred in Sumatera, Kalimantan and Java. Chikungunya was not reported in Papua and West Papua province of Indonesia between 2008 and 2016.
Seroprevalence of anti-CHIKV IgM antibodies
Five studies were identified that assessed the seroprevalence of anti-CHIKV IgM antibodies (with or without the combination of other tests), during both outbreak and non-outbreak settings [54, 55, 77, 78, 95]. These studies analyzed a total of 1183 serum samples between 1998 and 2010. The combination of all studies that were conducted in outbreak and non-outbreak settings yielded a median IgM seroprevalence of 13.3% (range 0.0–60.4%) (Table 2).
Our data indicate that the seroprevalence of anti-CHIKV IgM antibodies during an outbreak is much higher compared to in a non-outbreak setting (median 17.7%, range 3.2–60.4% vs. median 7.3%, range 0.0–26.3%, respectively). In some outbreaks the seroprevalence of IgM antibodies was more than 50%, such as the outbreak in Bogor 2001 . However, an investigation conducted 3 months after a chikungunya outbreak in Yogyakarta during 1998–1999 found that IgM antibodies were only present in 3.2% of suspected patients with CHIKV infection while IgG seroprevalence was at 70.4%  (Table 2).
Seroprevalence of anti-CHIKV IgG antibodies
Seven studies provided information on the seroprevalence of anti-CHIKV IgG antibodies in humans, in both outbreak and non-outbreak settings [49, 50, 54, 55, 77,78,79]. These studies covered a 47-year period (between 1969 and 2015) with a total of 3717 serum samples tested (Table 3). Overall, in both outbreak and non-outbreak settings, the median seroprevalence of IgG antibodies was 18.5% (ranging from 0.0% in non-febrile inhabitants in non-outbreak area to 73.1% among patients who were febrile during an outbreak). In a non-outbreak setting, for both febrile and non-febrile patients, the median seroprevalence of anti-CHIKV IgG was 14.1% (range 0.0–43.9%). In contrast, the seroprevalence of IgG antibodies was much higher among patients during an outbreak (both febrile and non-febrile individuals) (median 47.4%, range 0.0–73.1%). The median seroprevalence of anti-CHIKV IgG among non-febrile subjects in a non-outbreak setting was 12.2% (range 0.0 to 26.8%) in Indonesia.
Mortality of chikungunya in Indonesia
Five outbreak investigations [13, 54, 55, 78, 80] were identified and no deaths due to chikungunya were reported. In addition, in eleven annual reports from the MoH of Indonesia, covering 44 years (1973 to 2016), there was no deaths related to CHIKV infection [48, 56,57,58,59,60,61,62,63,64,65].
Genotypes of CHIKV circulating in Indonesia
Sixteen studies reporting the genotype of circulating Indonesian CHIKV were identified with a total of 130 viruses sequences [13, 20, 23, 76, 78, 81,82,83,84,85, 94, 98, 100, 102, 103, 108] (Table 4). Of these studies, seven were conducted among local inhabitants [13, 78, 81, 83, 85, 94, 108], and eight were conducted with viruses isolated from travellers returning from Indonesia [20, 23, 82, 84, 98, 100, 102, 103]. One study did not mention whether the virus was isolated in local inhabitant or in a traveller . Of the seven local studies, four were conducted in a non-outbreak setting [81, 83, 85, 94], two investigations were carried out during a chikungunya outbreak [13, 78] and one study did not specify the setting . Most of the viruses that were isolated from travellers were from Taiwan [23, 82, 100]. Other viruses were collected from travellers returning to Singapore , France , the Netherlands , Russia , and Germany . A vast majority of the 130 CHIKV isolated from Indonesia, belonged to Asian genotype (92.3%, 120/130). The remaining isolated viruses were identified as being viruses of the ECSA genotype (7.7%, 10/130). Of these ECSA viruses, two were isolated locally in 2011  and the remaining eight were isolated from travelers returning from Indonesia between 2008 and 2010 [20, 23, 102]. Viruses of the West African genotype have not yet been identified in or from Indonesia.
To confirm the results of the search from systematic review for CHIKV genotype literature and the database, a search was conducted in GenBank which yielded 119 sequences, rather than the expected 130 sequences. Discrepancy between these findings is due to three reasons: a) eight sequences were from an in press study and had not yet been deposited into GenBank when the search was conducted ; b) two sequences of CHIKV, which were isolated from Singaporean travelers returning from Indonesia, were deposited into GenBank as Singaporean viruses  and; c) one sequence had no information of country of origin in GenBank.
For phylogenetic and evolutionary analyses based on CHIKV E1 gene sequences, 127 sequences were included. One paper, reporting non-structural gene sequence of three Asian genotype viruses, was excluded . The relationship of Indonesian CHIKV to other CHIKV isolated worldwide is shown in Fig. 4. The Asian genotype was first identified in Indonesia in 1983 , while viruses of the ECSA genotype were isolated for the first time in 2008, following the isolation of this genotype from returning Taiwanese  and Singaporean travellers .
MCMC analysis using 53 representative reference sequences reveals that the ECSA viruses sampled in Indonesia between 2008 and 2011 belong to the Indian Ocean Lineage (IOL) and are closely related to viruses that have circulated in Southeast Asian countries such as Sri Lanka, Malaysia, Singapore, Thailand, and Myanmar, as well as China and South Korea during the same period (Fig. 5). Interestingly, viruses of the ECSA genotype have not been sampled in Indonesia, either in local or from travelers studies, from 2011 onwards.
Exported chikungunya cases
Twenty four studies, reporting exported chikungunya cases that originated from Indonesia covering a 26 year-period from 1989 to 2014, were identified [20, 23, 82, 84, 86,87,88,89,90,91,92,93, 96,97,98,99,100,101,102,103,104,105,106,107] (Table 5). During this period, a total of 195 cases of chikungunya cases were reported from travellers returning from Australia (128 cases) [86,87,88, 99, 104,105,106,107], Taiwan (47 cases) [23, 82, 100], Japan (4 cases) [89, 90] and other countries in Asia, Europe and the Pacific region (16 cases) [76, 84, 91,92,93, 96, 98, 101,102,103, 109]. Most of these cases were diagnosed with a combination of serology and RT-PCR, and sequence data was generated from only 27.6% (54/195) of cases [20, 23, 82, 84, 98, 100,101,102,103].
Our systematic review and analysis of MoH database examine the epidemiology of chikungunya in Indonesia by summarizing published studies on chikungunya incidence, the seroprevalence of IgM and IgG, mortality, CHIKV genotypes and exported cases. Existing literature on incidence rate of chikungunya Indonesia is lacking. Nevertheless, using MoH databases, we are able to generate comprehensive incidence rates of chikungunya in Indonesia over a 16 year period, from 2001 to 2016. Incidence data indicate there is no typical pattern of chikungunya in the country. The most notable finding is the massive nationwide outbreak during 2009–2010 with 137,655 cases followed by a smaller nationwide outbreak in 2013 with 15,324 cases (Fig. 2). Prior to 2008 the incidence of chikungunya was less than 10,000 cases/year and this increased significantly in 2009 and 2010 with 83,756 and 53,899 cases, respectively [61, 62] (Fig. 2) indicating that chikungunya incidence increased more than 20 times in 2009. There is no scientific evidence to explain what occurred during this outbreak. The MoH suggests that this increase was due to the fact that prior to 2009, many regions did not include CHIKV case reports in their annual reports to the MoH.
Despite the multiple outbreaks that have occurred in Indonesia since its re-emergence, there is also a lack of data regarding the magnitude of CHIKV exposure in the Indonesian population. Our systematic review, covering the period of 1998–2010, indicates that CHIKV-specific IgM antibodies were detected approximately 13% of sera collected from healthy and acute febrile illness patients (range 0.0–60.4%). This estimate can be interpreted as evidence of recent CHIKV infection. Although IgM may persist for months to years following the resolution of infection, IgM titers reduce to undetectable levels after weeks for most chikungunya patients and are replaced with long lasting IgG . The high seroprevalence of anti-CHIKV IgM indicates that chikungunya is endemic in Indonesia with on-going transmission.
Interestingly, 10.8% of sera from febrile patients (presenting with dengue-symptoms or undifferentiated acute febrile illness) in a non-outbreak setting had anti-CHIKV IgM [77, 95]. This finding suggests that many patients presenting with dengue-like symptoms had a CHIKV infection. This proposed consideration is supported by the evidence from other studies: a) in a non-outbreak setting, approximately 10% of febrile patients with dengue-like symptoms demonstrated titers of anti-CHIKV IgM ; and b) in most Asian countries, DENV and CHIKV co-circulate in the same location [111, 112]. The confusion between dengue and chikungunya infection could explain underreporting of CHIKV in Indonesia . Therefore, further tests is required to exclude CHIKV infection in patients with dengue-like symptoms. However, many regions in Indonesia lack a specific diagnostic test for chikungunya. In clinical settings in Indonesia, however, DENV infection must be considered first as the main diagnosis to minimize the incidence of severe forms of dengue and therefore to reduce dengue mortality. This rationale is based on the fact that dengue more frequently presents with severe forms compared to CHIKV, and has a high mortality rate in Indonesia .
We found that in both outbreak and non-outbreak localities, the median seroprevalence of asymptomatic residents with current CHIKV infection (the present of IgM) was 9.1% (range 0.0–48.0%) (Table 2). This finding highlights the important fact that not all CHIKV infections present with specific symptoms that are noticed by the infected individuals. A recent systematic review revealed that the proportion of CHIKV infections that were asymptomatic was highly variable, ranging from 3.2 to 82.1% . Interestingly, current studies indicate there is an increasing trend of asymptomatic chikungunya over time [7, 21, 27]. Except for one study, all studies published from 2013 onwards reported that the percentage of asymptomatic chikungunya was more than 30% .
Our systematic review indicates that the median seroprevalence of CHIKV-specific IgG antibodies was 18.5% (range 0.0–73.1%) among the total population in all settings. However, from an epidemiological perspective, data from non-febrile subjects collected in a non-outbreak setting is the best representation of the true seroprevalence of CHIKV infection. The median seroprevalence of anti-CHIKV IgG in this population was 12.7% (range 0.0–26.8%) in Indonesia. This figure is higher than the estimated seroprevalence of anti-CHIKV IgG antibodies in Malaysian adults, which is at 5.9% . As predicted, the median seroprevalence of IgG antibodies among healthy residents living in post-outbreak areas is much higher (approximately three times higher) compared to healthy residents in non-outbreak areas (median 44.8%, range 55.5–50.5% vs. median 12.7%, range 0.0–26.8%). This finding is comparable with other outbreak investigations, including 20.0% in Brazil (2016) , 26.8% in Thailand (2014) , 37.2% in Mayotte, Indian Ocean (2006)  and 67.9% in India (2007) . The variability between these estimates is influenced by several factors such as serum collection time relative to outbreak, criteria of sample selection (whether asymptomatic residents during the outbreak are included or excluded) and the location of the serosurvey relative to the epicentrum of outbreak.
Surprisingly, no chikungunya-associated deaths have been recorded in Indonesia since this disease was officially recognized by MoH of Indonesia in 1973 . Since 2004, chikungunya has been included in the National Diseases Surveillance run by the MoH and has been reported in the Annual Report of Indonesia Health Profile, the formal annual report from MoH of Indonesia, since 2005 . Although a substantial incidence rate has been reported, no deaths have been recorded since then. In the National Guideline of Prevention and Control of Chikungunya from MoH , there is no clear criteria for assessing or reporting chikungunya-related deaths that should be included in the surveillance system. Unsurprisingly, this lack of guidance has led to no reported deaths from chikungunya in Indonesia; and this number is probably due to underreporting of fatal chikungunya cases within the current passive surveillance system. In fact several studies in India [120, 121], Mauritius  and Brazil [122, 123] during ECSA genotype outbreaks and in Dominican Republic  during an Asian genotype outbreak revealed an increase in the mortality rate, with deaths not adequately identified by passive surveillance systems. In addition, those studies revealed an increase of excess mortality [12, 120, 122,123,124] indicating that the actual proportion of chikungunya-associated deaths is underestimated in many settings. In Indonesia, during the 2008–2009 chikungunya outbreaks, there were 137,655 chikungunya cases officially reported to MoH, with no reported death. In addition, given the fact that the illnesses caused by CHIKV and DENV are clinically indistinguishable and not all reported dengue cases in Indonesia are confirmed by laboratory testing, leads to the possibility that deaths occur due to chikungunya which may be wrongly attributed to dengue as a consequence of the prioritization of dengue surveillance. Chikungunya has been included in the Annual Report of Indonesia Health Profile since 2005 [48, 56,57,58,59,60,61,62,63,64,65], but was excluded in the 2016 MoH Annual Report . In the 2017 Annual Report, only 126 chikungunya cases were reported to MoH, which dropped from 1702 cases that were reported in 2016 . These points reflect that chikungunya is not a priority disease in Indonesia and that there is poor chikungunya surveillance and an insufficient reporting system in Indonesia. If fact, in some studies, higher fatality rate for chikungunya was reported compared to dengue [122, 127, 128]. Therefore, the MoH should provide clear criteria for chikungunya-associated death in the national guidelines and institute a clear reporting system. In addition, in outbreaks, active surveillance might be required to reduce underreporting of chikungunya deaths in the country.
To date, two genotypes of CHIKV have been isolated in Indonesia (Asian and ECSA). Viruses of the West African genotype have not yet been isolated from Indonesia. Viruses of the Asian genotype were first isolated in Indonesia 35 years ago . The ECSA genotype however was identified for the first time in 2008 [20, 100] and at the same time this genotype caused several major outbreaks in Southeast Asia countries [16, 100, 129]. This suggests that chikungunya outbreaks that occurred prior to 2008 in Indonesia were associated with viruses of the Asian genotype.
A phylogenetic analysis revealed that the ECSA viruses sampled from Indonesia during 2008–2011 are closely related with viruses that caused contemporaneous outbreaks in Southeast and East Asian countries, such as Malaysia (2008–2009) , Singapore (2008) [19, 20], Thailand (2008–2009) [16,17,18], and China (2010) . This suggests that the viruses in Indonesia may have been introduced from another country in Southeast Asia, most likely Malaysia . In 2008, a nationwide chikungunya outbreak occurred in Malaysia, which was found to be caused by viruses of the ECSA genotype . In the same year ECSA viruses were also reported in Indonesia [20, 100] and continued to be sampled until 2011 . Furthermore, our phylogenetic analysis revealed that the ECSA viruses that circulated in Southeast Asia, including Indonesia, were introduced from India and Sri Lanka, where CHIKV circulated in 2005–2007 and infected more than 1.3 million persons during 2005–2006 in India alone .
Our data, together with other evidence, suggest that the introduction of the ECSA genotype to Indonesia may be the reason for the large chikungunya outbreak reported between 2009 and 2010 in the country (Fig. 2). First, the Asian genotype had been the dominant genotype of Indonesia and had circulated for several decades. However, the annual number of cases never exceeded more than 10,000 cases prior to 2009. Second, the ECSA was first reported in Indonesia in 2008, and had not been isolated in the multiple studies of Indonesian CHIKV (conducted both locally and in neighboring countries) prior to this time [81, 82, 94, 100]. This evidence indicates a time relationship between the introduction of the ECSA genotype into Indonesia and a sharp increase of reported chikungunya cases. Third, within the same time-frame, nationwide outbreaks of chikungunya associated with viruses of the ECSA genotype occurred in Malaysia , Singapore  and Thailand . The ECSA viruses isolated from these outbreaks are all closely related, indicating that these outbreaks were probably caused by the same type of virus that then dispersed throughout the region and likely into Indonesia as well (Fig. 5). Other scientific evidence also supports this finding of chikungunya outbreaks in Southeast Asian countries during that time-frame being driven by the emergence of the same strain of ECSA . In addition, MoH data reveal that the highest number of CHIKV cases in 2009 were reported in Bangka Belitung and provinces in Sumatera Island (Jambi, South Sumatera, Lampung, Bengkulu, North Sumatera and Aceh) . These provinces are in close geographical proximity to Malaysia and Singapore (Fig. 3). Finally, evolutionally analysis reveals that ECSA viruses isolated from Indonesia are closely related to those viruses associated with massive outbreaks in Malaysia, Singapore and Thailand (Fig. 5). This possibility is also supported by the fact that ECSA is more adaptive to Ae. albopictus due to a mutation in the envelope protein gene (E1-A226V) . This mutation associated with a significant increase in viral infectivity for Ae. albopictus, and led to more efficient viral dissemination . In the area where both Ae. aegypti and Ae. albopictus are endemic, like Indonesia, introduction of ECSA viruses could cause a massive outbreak. This evidence, in aggregate, suggests that the massive outbreaks in Indonesia during 2009 and 2010 were caused by the introduction of the ECSA genotype to the country.
There are some limitations of this study that should be discussed. First, full-texts of some potential references that assess the seroprevalence of the anti-CHIKV antibodies were not available, especially those published prior to 1985 [133,134,135,136,137,138,139,140,141,142]. Second, in the seroprevalence section of this paper, we included a study that utilized the HI test as a diagnostic method . This test is quite sensitive for detecting alphavirus antibodies; however, cross-reactions often occur among the viruses of the same group. Nevertheless, a great majority of sera from this study reacted with only one of the alphaviruses by the HI test and therefore the results from that study were valid. Third, the collection time of the sera, relative to time of infection, among studies were not the same. This might influence the percentage of samples that were positive for IgM and IgG among studies. For example, an outbreak investigation in Yogyakarta reported a significantly lower of seroprevalence of anti-CHIKV IgM and a higher IgG seroprevalence compared to other studies, as the investigation was conducted after the resolution of the outbreak . Fourth, in Indonesia the absence of serological surveillance means that the precise numbers of chikungunya cases reported to the MoH during outbreaks are often diagnosed clinically only. Fifth, some of the known CHIKV of the Asian genotype, that have been isolated previously, were not available in GenBank . Finally, some CHIKV sequences originating from Indonesia may have been reported as deriving from elsewhere, if a returning traveller had an incomplete or inaccurate travel history. Studies may have failed to report Indonesia as the country of origin of isolated viruses, and instead, reported the country in which the virus was isolated as the origin . Despite these limitations, this study is, to the best of our knowledge, the first systematic review on CHIKV infection in Indonesia. Strengths also include a search strategy that avoided missing potential references and Indonesian CHIKV sequences. In addition, this study is also able to provide the incidence rate of chikungunya for a 16-year period using the National Diseases Surveillance database from the MoH.
Evidence suggests that CHIKV has circulated in most of the Indonesian archipelago for at least 50 years, despite the first virologically confirmed outbreak being recorded as recently as 1982. Although anti-CHIKV IgM and IgG seroprevalence is high and two massive nationwide outbreaks have been recorded in Indonesia, no deaths have been reported. The lack of reported deaths possibly reflects under-reporting of fatal cases in the country using the current surveillance system. Data from both locals and travelers revealed that the Asian genotype has been identified in Indonesia for more than 30 years, while viruses of the ECSA genotype were sampled only during the 2008–2011 period. A phylogenetic analysis reveals that the ECSA viruses sampled in Indonesia during this period are closely related to viruses that have circulated in neighboring Southeast Asian countries, as well as China and South Korea, within the same timeframe. It is likely that ECSA viruses were introduced to Indonesia from a neighboring country within the region, in 2008. Spatio-temporal, phylogenetic and evolutionary data suggest that this genotype was likely associated with the large 2009–2010 chikungunya outbreak in Indonesia.
Bayesian Evolutionary Analysis Sampling Trees
Enzyme-linked immunosorbent assay
- GTR + Γ:
General Time Reversible, gamma substitution model
haemagglutination inhibition assay
Multiple Alignment using Fast Fourier Transform
Maximum Clade Credibility
Bayesian Markov chain Monte Carlo
Ministry of Health, Republic of Indonesia
Preferred reporting items for systematic reviews and meta-analyses
Plaque reduction neutralization test
Randomized Axelerated Maximum Likelihood
Reverse transcription polymerase chain reaction
Time to most recent common ancestors
World Health Organization
Weaver SC, Reisen WK. Present and future arboviral threats. Antivir Res. 2010;85(2):328–45.
Shragai T, Tesla B, Murdock C, Harrington LC. Zika and chikungunya: mosquito-borne viruses in a changing world. Ann N Y Acad Sci. 2017;1399(1):61–77.
Schwartz O, Albert ML. Biology and pathogenesis of chikungunya virus. Nat Rev Microbiol. 2010;8(7):491–500.
Burt FJ, Chen W, Miner JJ, Lenschow DJ, Merits A, Schnettler E, Kohl A, Rudd PA, Taylor A, Herrero LJ, et al. Chikungunya virus: An update on the biology and pathogenesis of this emerging pathogen. Lancet Infect Dis. 2017;17(4):e107–17.
Robinson MC. An epidemic of virus disease in Southern Province, Tanganyika territory, in 1952-53. I. Clinical features. Trans R Soc Trop Med Hyg. 1955;49(1):28–32.
Lumsden WH. An epidemic of virus disease in Southern Province, Tanganyika territory, in 1952-53. II. General description and epidemiology. Trans R Soc Trop Med Hyg. 1955;49(1):33–57.
Weaver SC, Lecuit M. Chikungunya virus and the global spread of a mosquito-borne disease. N Engl J Med. 2015;372(13):1231–9.
WHO. Outbreak and spread of chikungunya. Wkly Epidemiol Rec. 2007;82(47):409–15.
Renault P, Solet JL, Sissoko D, Balleydier E, Larrieu S, Filleul L, Lassalle C, Thiria J, Rachou E, de Valk H, et al. A major epidemic of chikungunya virus infection on Reunion Island, France, 2005-2006. Am J Trop Med Hyg. 2007;77(4):727–31.
Sissoko D, Malvy D, Giry C, Delmas G, Paquet C, Gabrie P, Pettinelli F, Sanquer MA, Pierre V. Outbreak of chikungunya fever in Mayotte, Comoros archipelago, 2005-2006. Trans R Soc Trop Med Hyg. 2008;102(8):780–6.
Jain M, Rai S, Chakravarti A. Chikungunya: a review. Trop Dr. 2008;38(2):70–2.
Beesoon S, Funkhouser E, Kotea N, Spielman A, Robich RM. Chikungunya fever, Mauritius, 2006. Emerg Infect Dis. 2008;14(2):337–8.
Maha MS, Susilarini NK, Hariastuti NI. Chikungunya virus mutation, Indonesia, 2011. Emerg Infect Dis. 2015;21(2):379–81.
Dash PK, Parida MM, Santhosh SR, Verma SK, Tripathi NK, Ambuj S, Saxena P, Gupta N, Chaudhary M, Babu JP, et al. East Central South African genotype as the causative agent in reemergence of chikungunya outbreak in India. Vector-Borne Zoonot. 2007;7(4):519–27.
Razmy A. Clinical features of chikungunya infection in Sri Lanka. Asian Pac J Trop Dis. 2014;4(2):131–4.
Rianthavorn P, Prianantathavorn K, Wuttirattanakowit N, Theamboonlers A, Poovorawan Y. An outbreak of chikungunya in southern Thailand from 2008 to 2009 caused by African strains with A226V mutation. Int J Infect Dis. 2010;14:E161–5.
Theamboonlers A, Rianthavorn P, Praianantathavorn K, Wuttirattanakowit N, Poovorawan Y. Clinical and molecular characterization of chikungunya virus in South Thailand. Jpn J Infect Dis. 2009;62(4):303–5.
Sasayama M, Benjathummarak S, Kawashita N, Rukmanee P, Sangmukdanun S, Masrinoul P, Pitaksajjakul P, Puiprom O, Wuthisen P, Kurosu T, et al. Chikungunya virus was isolated in Thailand, 2010. Virus Genes. 2014;49(3):485–9.
Leo YS, Chow ALP, Tan LK, Lye DC, Lin L, Ng LC. Chikungunya outbreak, Singapore, 2008. Emerg Infect Dis. 2009;15(5):836–7.
Ng LC, Tan LK, Tan CH, Tan SS, Hapuarachchi HC, Pok KY, Lai YL, Lam-Phua SG, Bucht G, Lin RT, et al. Entomologic and virologic investigation of chikungunya, Singapore. Emerg Infect Dis. 2009;15(8):1243–9.
Gerardin P, Guernier V, Perrau J, Fianu A, Le Roux K, Grivard P, Michault A, de Lamballerie X, Flahault A, Favier F. Estimating chikungunya prevalence in La Reunion Island outbreak by serosurveys: two methods for two critical times of the epidemic. BMC Infect Dis. 2008;8:99.
Carey DE. Chikungunya and dengue: a case of mistaken identity? J Hist Med Allied Sci. 1971;26(3):243–62.
Yang CF, Su CL, Hsu TC, Chang SF, Lin CC, Huang JC, Shu PY. Imported chikungunya virus strains, Taiwan, 2006-2014. Emerg Infect Dis. 2016;22(11):1981–4.
Endy TP, Chunsuttiwat S, Nisalak A, Libraty DH, Green S, Rothman AL, Vaughn DW, Ennis FA. Epidemiology of inapparent and symptomatic acute dengue virus infection: a prospective study of primary school children in Kamphaeng Phet, Thailand. Am J Epidemiol. 2002;156(1):40–51.
Carrillo FB, Gordon A, Harris E. Reply to Gérardin et al. Clin Infect Dis. 2018;68(1):172–4.
Yoon IK, Alera MT, Lago CB, Tac-An IA, Villa D, Fernandez S, Thaisomboonsuk B, Klungthong C, Levy JW, Velasco JM, et al. High rate of subclinical chikungunya virus infection and association of neutralizing antibody with protection in a prospective cohort in the Philippines. PLoS Negl Trop Dis. 2015;9(5):e0003764.
Mourya DT, Mishra AC. Chikungunya fever. Lancet. 2006;368(9531):186–7.
Gerardin P, Couderc T, Bintner M, Tournebize P, Renouil M, Lemant J, Boisson V, Borgherini G, Staikowsky F, Schramm F, et al. Chikungunya virus-associated encephalitis: a cohort study on La Reunion Island, 2005-2009. Neurology. 2016;86(1):94–102.
Mehta R, Gerardin P, de Brito CAA, Soares CN, Ferreira MLB, Solomon T. The neurological complications of chikungunya virus: a systematic review. Rev Med Virol. 2018;28(3):e1978.
Lemant J, Boisson V, Winer A, Thibault L, Andre H, Tixier F, Lemercier M, Antok E, Cresta MP, Grivard P, et al. Serious acute chikungunya virus infection requiring intensive care during the Reunion Island outbreak in 2005-2006. Crit Care Med. 2008;36(9):2536–41.
Rajapakse S, Rodrigo C, Rajapakse A. Atypical manifestations of chikungunya infection. Trans R Soc Trop Med Hyg. 2010;104(2):89–96.
Cerny T, Schwarz M, Schwarz U, Lemant J, Gerardin P, Keller E. The range of neurological complications in chikungunya fever. Neurocrit Care. 2017;27(3):447–57.
Chua HH, Abdul Rashid K, Law WC, Hamizah A, Chem YK, Khairul AH, Chua KB. A fatal case of chikungunya virus infection with liver involvement. Med J Malaysia. 2010;65(1):83–4.
Mercado M, Acosta-Reyes J, Parra E, Guzman L, Beltran M, Gasque P, Mejia-Garcia C, Viasus D. Renal involvement in fatal cases of chikungunya virus infection. J Clin Virol. 2018;103:16–8.
Rolle A, Schepers K, Cassadou S, Curlier E, Madeux B, Hermann-Storck C, Fabre I, Lamaury I, Tressieres B, Thiery G, et al. Severe Sepsis and septic shock associated with chikungunya virus infection, Guadeloupe, 2014. Emerg Infect Dis. 2016;22(5):891–4.
Gauri LA, Ranwa BL, Nagar K, Vyas A, Fatima Q. Post chikungunya brain stem encephalitis. J Assoc Physicians India. 2012;60:68–70.
Hoz JM, Bayona B, Viloria S, Accini JL, Juan-Vergara HS, Viasus D. Fatal cases of chikungunya virus infection in Colombia: diagnostic and treatment challenges. J Clin Virol. 2015;69:27–9.
Gerardin P, Samperiz S, Ramful D, Boumahni B, Bintner M, Alessandri JL, Carbonnier M, Tiran-Rajaoefera I, Beullier G, Boya I, et al. Neurocognitive outcome of children exposed to perinatal mother-to-child chikungunya virus infection: the CHIMERE cohort study on Reunion Island. PLoS Negl Trop Dis. 2014;8(7):e2996.
van Aalst M, Nelen CM, Goorhuis A, Stijnis C, Grobusch MP. Long-term sequelae of chikungunya virus disease: a systematic review. Travel Med Infect Dis. 2017;15:8–22.
Zaid A, Gerardin P, Taylor A, Mostafavi H, Malvy D, Mahalingam S. Chikungunya arthritis: implications of acute and chronic inflammation mechanisms on disease management. Arthritis Rheumatol. 2018;70(4):484–95.
Paixao ES, Rodrigues LC, Costa M, Itaparica M, Barreto F, Gerardin P, Teixeira MG. Chikungunya chronic disease: a systematic review and meta-analysis. Trans R Soc Trop Med Hyg. 2018;112(7):301–16.
Duvignaud A, Fianu A, Bertolotti A, Jaubert J, Michault A, Poubeau P, Fred A, Mechain M, Gauzere BA, Favier F, et al. Rheumatism and chronic fatigue, the two facets of post-chikungunya disease: the TELECHIK cohort study on Reunion island. Epidemiol Infect. 2018;146(5):633–41.
Rodriguez-Morales AJ, Cardona-Ospina JA, Fernanda Urbano-Garzon S, Sebastian Hurtado-Zapata J. Prevalence of post-chikungunya infection chronic inflammatory arthritis: a systematic review and meta-analysis. Arthritis Care Res (Hoboken). 2016;68(12):1849–58.
Marimoutou C, Vivier E, Oliver M, Boutin JP, Simon F. Morbidity and impaired quality of life 30 months after chikungunya infection: comparative cohort of infected and uninfected French military policemen in Reunion Island. Medicine. 2012;91(4):212–9.
Gerardin P, Fianu A, Malvy D, Mussard C, Boussaid K, Rollot O, Michault A, Gauzere BA, Breart G, Favier F. Perceived morbidity and community burden after a chikungunya outbreak: the TELECHIK survey, a population-based cohort study. BMC Med. 2011;9:5.
Caglioti C, Lalle E, Castilletti C, Carletti F, Capobianchi MR, Bordi L. Chikungunya virus infection: an overview. New Microbiol. 2013;36:211–27
Wahid B, Ali A, Rafique S, Idrees M. Global expansion of chikungunya virus: mapping the 64-year history. Int J Infect Dis. 2017;58:69–76.
Kemenkes RI. Profil kesehatan Indonesia tahun 2005. Jakarta: Kementerian Kesehatan RI; 2007.
Tesh RB, Gajdusek DC, Garruto RM, Cross JH, Rosen L. The distribution and prevalence of group a arbovirus neutralizing antibodies among human populations in Southeast Asia and the Pacific islands. Am J Trop Med Hyg. 1975;24(4):664–75.
Kanamitsu M, Taniguchi K, Urasawa S, Ogata T, Wada Y, Saroso JS. Geographic distribution of arbovirus antibodies in indigenous human populations in the indo-Australian archipelago. Am J Trop Med Hyg. 1979;28(2):351–63.
Christie J. On epidemics of dengue fever: their diffusion and etiology. Glasgow Med J. 1881;3:161–76.
Halstead SB. Reappearance of chikungunya, formerly called dengue, in the Americas. Emerg Infect Dis. 2015;21(4):557–61.
Weaver SC, Forrester NL. Chikungunya: evolutionary history and recent epidemic spread. Antivir Res. 2015;120:32–9.
Porter KR, Tan R, Istary Y, Suharyono W, Sutaryo WS, Ma'Roef C, Listiyaningsih E, Kosasih H, Hueston L, et al. A serological study of chikungunya virus transmission in Yogyakarta, Indonesia: evidence for the first outbreak since 1982. Southeast Asian J Trop Med Public Health. 2004;35(2):408–15.
Laras K, Sukri NC, Larasati RP, Bangs MJ, Kosim R, Djauzi WT, Master J, Kosasih H, Hartati S, et al. Tracking the re-emergence of epidemic chikungunya virus in Indonesia. Trans R Soc Trop Med Hyg. 2005;99(2):128–41.
Kemenkes RI. Profil kesehatan Indonesia tahun 2014. Jakarta: Kementerian Kesehatan RI; 2015.
Kemenkes RI. Profil kesehatan Indonesia tahun 2015. Jakarta: Kementerian Kesehatan RI; 2016.
Kemenkes RI. Profil kesehatan Indonesia tahun 2013. Jakarta: Kementerian Kesehatan RI; 2014.
Kemenkes RI. Profil kesehatan Indonesia tahun 2012. Jakarta: Kementerian Kesehatan RI; 2013.
Kemenkes RI. Profil kesehatan Indonesia tahun 2011. Jakarta: Kementerian Kesehatan RI; 2012.
Kemenkes RI. Profil kesehatan Indonesia tahun 2010. Jakarta: Kementerian Kesehatan RI; 2011.
Kemenkes RI. Profil kesehatan Indonesia tahun 2009. Jakarta: Kementerian Kesehatan RI; 2010.
Kemenkes RI. Profil kesehatan Indonesia tahun 2008. Jakarta: Kementerian Kesehatan RI; 2009.
Kemenkes RI. Profil kesehatan Indonesia tahun 2007. Jakarta: Kementerian Kesehatan RI; 2008.
Kemenkes RI. Profil kesehatan Indonesia tahun 2006. Jakarta: Kementerian Kesehatan RI; 2007.
Moher D, Liberati A, Tetzlaff J, Altman DG, Grp P. Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. PLoS Med. 2009;6(7):e1000097.
Humphrey JM, Cleton NB, Reusken C, Glesby MJ, Koopmans MPG, Abu-Raddad LJ. Urban chikungunya in the Middle East and North Africa: a systematic review. PLoS Negl Trop Dis. 2017;11(6):e0005707.
Kemenkes RI. Pedoman pencegahan dan pengendalian demam chikungunya di Indonesia. In. Jakarta: Kementerian Kesehatan RI; 2017.
WHO. Guidelines for prevention & control of chikungunya fever. India: World Health Organization, Regional Office for South-East Asia; 2009.
Environmental Systems Research Institute (ESRI). ArcGIS Release 10.1. Redlands, CA.: ESRI; 2012.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, et al. Geneious basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28(12):1647–9.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80.
Stamatakis A, Hoover P, Rougemont J. A rapid bootstrap algorithm for the RAxML web servers. Syst Biol. 2008;57(5):758–71.
Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22(21):2688–90.
Bouckaert R, Heled J, Kuhnert D, Vaughan T, Wu CH, Xie D, Suchard MA, Rambaut A, Drummond AJ. BEAST 2: A software platform for Bayesian evolutionary analysis. Plos Comput Biol. 2014;10(4):e1003537.
Volk SM, Chen R, Tsetsarkin KA, Adams AP, Garcia TI, Sall AA, Nasar F, Schuh AJ, Holmes EC, Higgs S, et al. Genome-scale phylogenetic analyses of chikungunya virus reveal independent emergences of recent epidemics and various evolutionary rates. J Virol. 2010;84(13):6497–504.
Ngwe T, Inoue S, Thant KZ, Talemaitoga N, Aryati A, Dimaano EM, Matias RR, Buerano CC, Natividad FF, Abeyewickreme W, et al. Retrospective seroepidemiological study of chikungunya infection in South Asia, Southeast Asia and the Pacific region. Epidemiol Infect. 2016;144(11):2268–75.
Sari K, Myint K, Andayani A, Adi P, Dhenni R, Perkasa D, Maroef C, Witaria N, Megawati D, Powers A, et al. Chikungunya fever outbreak identified in North Bali, Indonesia. Trans R Soc Trop Med Hyg. 2017; [In press].
Suharti C, van Gorp EC, Dolmans WM, Groen J, Hadisaputro S, Djokomoeljanto RJ, MEO D, van der Meer JW. Hanta virus infection during dengue virus infection outbreak in Indonesia. Acta medica Indonesiana. 2009;41(2):75–80.
Sitepu Y, Suprayogi A, Pramono D. Risk factors for an outbreak of chikungunya fever in West Borneo, Indonesia, 2010. Am J Epidemiol. 2011;173:S328.
Riswari SF, Ma'roef CN, Djauhari H, Kosasih H, Perkasa A, Yudhaputri FA, Artika IM, Williams M, van der Ven A, Myint KS, et al. Study of viremic profile in febrile specimens of chikungunya in Bandung, Indonesia. J Clin Virol. 2016;74:61–5.
Shu PY, Yang CF, Su CL, Chen CY, Chang SF, Tsai KH, Cheng CH, Huang JH. Two imported chikungunya cases, Taiwan. Emerg Infect Dis. 2008;14(8):1326–7.
Mulyatno KC, Susilowati H, Yamanaka A, Soegijanto S, Konishi E. Primary isolation and phylogenetic studies of chikungunya virus from Surabaya, Indonesia. Jpn J Infect Dis. 2012;65(1):92–4.
Shchelkanov M, L'Vov DK, Kolobukhina LV, Al’khovskii SV, Shchetinin AM, Saifullin MA, Kruzhkova IS, Aristova VA, Morozova TV, Samokhvalov EI, et al. Isolation of the chikungunya virus in Moscow from the Indonesian visitor (September, 2013). Vopr Virusol. 2014;59(3):28–34.
Sasmono R, Perkasa A, Yohan B, Haryanto S, Yudhaputri FA, Hayati R, Maaroef C, Ledermann J, Myint K, Powers A. Chikungunya deterection during dengue outbreak in Sumatera, Indonesia: clinical manifestatation and virological profile. Am J Trop Med Hyg. 2017;97(5):1393–8.
Knope KE, Doggett SL, Kurucz N, Johansen CA, Nicholson J, Feldman R, Sly A, Hobby M, El Saadi D, Muller M, et al. Arboviral diseases and malaria in Australia, 2011-12: annual report of the National Arbovirus and malaria advisory committee. Commun Dis Intell Q Rep. 2014;38(2):E122–42.
Knope KE, Kurucz N, Doggett SL, Muller M, Johansen CA, Feldman R, Hobby M, Bennett S, Sly A, Lynch S, et al. Arboviral diseases and malaria in Australia, 2012-13: annual report of the National Arbovirus and malaria advisory committee. Commun Dis Intell Q Rep. 2016;40(1):E17–47.
Knope KE, Muller M, Kurucz N, Doggett SL, Feldman R, Johansen CA, Hobby M, Bennett S, Lynch S, Sly A, et al. Arboviral diseases and malaria in Australia, 2013-14: annual report of the National Arbovirus and malaria advisory committee. Commun Dis Intell Q Rep. 2016;40(3):E400–36.
Kobashi K, Kobayashi T, Nakamura-Uchiyama F, Ohnishi K. A Japanese patient with chikungunya fever returning from Flores Island, Indonesia. Kansenshogaku Zasshi. 2010;84(4):457–9.
Mizuno Y, Kato Y, Takeshita N, Ujiie M, Kobayashi T, Kanagawa S, Kudo K, Lim CK, Takasaki T. Clinical and radiological features of imported chikungunya fever in Japan: a study of six cases at the National Center for Global Health and medicine. J Infect Chemother. 2011;17(3):419–23.
Cha GW, Cho JE, Lee EJ, Ju YR, Han MG, Park C, Jeong YE. Travel-associated chikungunya cases in South Korea during 2009-2010. Osong Public Health Res Perspectives. 2013;4(3):170–5.
Panning M, Grywna K, van Esbroeck M, Emmerich P, Drosten C. Chikungunya fever in travelers returning to Europe from the Indian Ocean region, 2006. Emerg Infect Dis. 2008;14(3):416–22.
Chaves T, Pellini A, Mascheretti M, Jahnel M, Ribeiro A, Rodrigues S, Vasconcelos P, Boulos M. Travelers as sentinels for chikungunya fever, Brazil. Emerg Infect Dis. 2012;18(3):529–30.
Kosasih H, de Mast Q, Widjaja S, Sudjana P, Antonjaya U, Ma'roef C, Riswari SF, Porter KR, Burgess TH, Alisjahbana B, et al. Evidence for endemic chikungunya virus infections in Bandung, Indonesia. PLoS Negl Trop Dis. 2013;7(10):e2483.
Capeding MR, Chua MN, Hadinegoro SR, Hussain II, Nallusamy R, Pitisuttithum P, Rusmil K, Thisyakorn U, Thomas SJ, Huu Tran N, et al. Dengue and other common causes of acute febrile illness in Asia: an active surveillance study in children. PLoS Negl Trop Dis. 2013;7(7):e2331.
Albuquerque IG, Marandino R, Mendonca AP, Nogueira RM, Vasconcelos PF, Guerra LR, Brandao BC, Aguiar GR, Bacco PA. Chikungunya virus infection: report of the first case diagnosed in Rio de Janeiro, Brazil. Rev Soc Bras Med Trop. 2012;45(1):128–9.
Dupont-Rouzeyrol M, Caro V, Guillaumot L, Vazeille M, D'Ortenzio E, Thiberge JM, Baroux N, Gourinat AC, Grandadam M, Failloux AB. Chikungunya virus and the mosquito vector Aedes aegypti in New Caledonia (South Pacific Region). Vector Borne Zoonotic Dis. 2012;12(12):1036–41.
Grandadam M, Caro V, Plumet S, Thiberge JM, Souares Y, Failloux AB, Tolou HJ, Budelot M, Cosserat D, Leparc-Goffart I, et al. Chikungunya virus, southeastern France. Emerg Infect Dis. 2011;17(5):910–3.
Harnett GB, Bucens MR. Isolation of chikungunya virus in Australia. Med J Aust. 1990;152(6):328–9.
Huang JH, Yang CF, Su CL, Chang SF, Cheng CH, Yu SK, Lin CC, Shu PY. Imported chikungunya virus strains, Taiwan, 2006-2009. Emerg Infect Dis. 2009;15(11):1854–6.
Roth A, Hoy D, Horwood PF, Ropa B, Hancock T, Guillaumot L, Rickart K, Frison P, Pavlin B, Souares Y. Preparedness for threat of chikungunya in the pacific. Emerg Infect Dis. 2014;20(8):e130696.
Scholte FE, Tas A, Martina BE, Cordioli P, Narayanan K, Makino S, Snijder EJ, van Hemert MJ. Characterization of synthetic chikungunya viruses based on the consensus sequence of recent E1-226V isolates. PLoS One. 2013;8(8):e71047.
Wolfel S, Vollmar P, Poluda D, Zange S, Antwerpen MH, Loscher T, Dobler G. Complete genome sequence of a chikungunya virus imported from Bali to Germany. Genome announcements. 2015;3(2):e00164–15.
Wright P, Fitzsimmons GJ, Johansen CA, Whelan PI, Co NAMA. Arboviral diseases and malaria in Australia, 2009-10: annual report of the National Arbovirus and malaria advisory committee. Commun Dis Intell Q Rep. 2012;36(1):70–81.
Knope K, Whelan P, Smith D, Johansen C, Moran R, Doggett S, Sly A, Hobby M, Kurucz N, Wright P, et al. Arboviral diseases and malaria in Australia, 2010-11: annual report of the National Arbovirus and malaria advisory committee. Commun Dis Intell Q Rep. 2013;37(1):E1–E20.
Fitzsimmons GJ, Wright P, Johansen CA, Whelan PI, Ad NAM. Arboviral diseases and malaria in Australia, 2007/08: annual report of the National Arbovirus and malaria advisory committee. Commun Dis Intell Q Rep. 2009;33(2):155–69.
Fitzsimmons GJ, Wright P, Johansen CA, Whelan PI, Arbovirus N. Arboviral diseases and malaria in Australia, 2008-09: annual report of the National Arbovirus and malaria advisory committee. Commun Dis Intell Q Rep. 2010;34(3):225–40.
Powers AM, Brault AC, Tesh RB, Weaver SC. Re-emergence of chikungunya and O'nyong-nyong viruses: evidence for distinct geographical lineages and distant evolutionary relationships. J Gen Virol. 2000;81(Pt 2):471–9.
Alibert A, Pfannstiel A, Grangeon J. Chikungunya outbreak in New Caledonia in 2011, Status report as at 22 August 2011. Inform’Action. 2011;34:3–9.
Velasco JM, Valderama MT, Lopez MN, Chua D, Latog R, Roque V, Corpuz J, Klungthong C, Rodpradit P, Hussem K, et al. Chikungunya virus infections among patients with dengue-like illness at a tertiary care hospital in the Philippines, 2012-2013. Am J Trop Med Hyg. 2015;93(6):1318–24.
Furuya-Kanamori L, Liang SH, Milinovich G, Magalhaes RJS, Clements ACA, Hu WB, Brasil P, Frentiu FD, Dunning R, Yakob L. Co-distribution and co-infection of chikungunya and dengue viruses. BMC Infect Dis. 2016;16:84.
Kaur M, Singh K, Sidhu SK, Devi P, Kaur M, Soneja S, Singh N. Coinfection of chikungunya and dengue viruses: A serological study from North Western region of Punjab, India. J Lab Physicians. 2018;10(4):443–7.
de Jong W, Rusli M, Bhoelan S, Rohde S, Rantam FA, Noeryoto PA, Hadi U, Gorp E, Goeijenbier M. Endemic and emerging acute virus infections in Indonesia: an overview of the past decade and implications for the future. Crit Rev Microbiol. 2018;44(4):487–503.
Karyanti MR, Uiterwaal CSPM, Kusriastuti R, Hadinegoro SR, Rovers MM, Heesterbeek H, Hoes AW, Bruijning-Verhagen P. The changing incidence of dengue haemorrhagic fever in Indonesia: a 45-year registry-based analysis. BMC Infect Dis. 2014;14.
Azami NA, Salleh SA, Shah SA, Neoh HM, Othman Z, Zakaria SZ, Jamal R. Emergence of chikungunya seropositivity in healthy Malaysian adults residing in outbreak-free locations: chikungunya seroprevalence results from the Malaysian cohort. BMC Infect Dis. 2013;13:67.
Cunha RV, Trinta KS, Montalbano CA, Sucupira MV, de Lima MM, Marques E, Romanholi IH, Croda J. Seroprevalence of chikungunya virus in a rural community in Brazil. PLoS Negl Trop Dis. 2017;11(1):e0005319.
Vongpunsawad S, Intharasongkroh D, Thongmee T, Poovorawan Y. Seroprevalence of antibodies to dengue and chikungunya viruses in Thailand. PLoS One. 2017;12(6):e0180560.
Sissoko D, Moendandze A, Malvy D, Giry C, Ezzedine K, Solet JL, Pierre V. Seroprevalence and risk factors of chikungunya virus infection in Mayotte, Indian Ocean, 2005-2006: a population-based survey. PLoS One. 2008;3(8):e3066.
Kumar NP, Suresh A, Vanamail P, Sabesan S, Krishnamoorthy KG, Mathew J, Jose VT, Jambulingam P. Chikungunya virus outbreak in Kerala, India, 2007: a seroprevalence study. Memorias do Instituto Oswaldo Cruz. 2011;106(8):912–6.
Mavalankar D, Shastri P, Bandyopadhyay T, Parmar J, Ramani KV. Increased mortality rate associated with chikungunya epidemic, Ahmedabad, India. Emerg Infect Dis. 2008;14(3):412–5.
Manimunda SP, Mavalankar D, Bandyopadhyay T, Sugunan AP. Chikungunya epidemic-related mortality. Epidemiol Infect. 2011;139(9):1410–2.
Brito CAA, Teixeira MG. Increased number of deaths during a chikungunya epidemic in Pernambuco, Brazil. Memorias do Instituto Oswaldo Cruz. 2017;112(9):650–1.
Freitas ARR, Cavalcanti L, Von Zuben AP, Donalisio MR. Excess mortality related to chikungunya epidemics in the context of Co-circulation of other arboviruses in Brazil. PLoS Curr. 2017;9.
Freitas ARR, Alarcon-Elbal PM, Paulino-Ramirez R, Donalisio MR. Excess mortality profile during the Asian genotype chikungunya epidemic in the Dominican Republic, 2014. Trans R Soc Trop Med Hyg. 2018;112(10):443–9.
Kemenkes RI. Profil kesehatan Indonesia tahun 2016. Jakarta: Kementerian Kesehatan RI; 2017.
Kemenkes RI. Profil kesehatan Indonesia tahun 2017. Jakarta: Kementerian Kesehatan RI; 2018.
Cavalcanti LPG, Freitas ARR, Brasil P, Cunha RVD. Surveillance of deaths caused by arboviruses in Brazil: from dengue to chikungunya. Memorias do Instituto Oswaldo Cruz. 2017;112(8):583–5.
Freitas ARR, Alarcon-Elbal PM, Donalisio MR. Excess mortality in Guadeloupe and Martinique, islands of the French West Indies, during the chikungunya epidemic of 2014. Epidemiol Infect. 2018;146(16):2059–65.
Sam IC, Chan YF, Chan SY, Loong SK, Chin HK, Hooi PS, Ganeswrie R, Abubakar S. Chikungunya virus of Asian and Central/East African genotypes in Malaysia. J Clin Virol. 2009;46(2):180–3.
Wu D, Wu J, Zhang QL, Zhong HJ, Ke CW, Deng XL, Guan DW, Li H, Zhang YH, Zhou HQ, et al. Chikungunya outbreak in Guangdong Province, China, 2010. Emerg Infect Dis. 2012;18(3):493–5.
Pulmanausahakul R, Roytrakul S, Auewarakul P, Smith DR. Chikungunya in Southeast Asia: understanding the emergence and finding solutions. Int J Infect Dis. 2011;15(10):E671–6.
Tsetsarkin KA, Vanlandingham DL, McGee CE, Higgs S. A single mutation in chikungunya virus affects vector specificity and epidemic potential. PLoS Pathog. 2007;3(12):e201.
Hotta S, Aoki H, Samoto S, Yasui T, Kawabe M. Virologic-epidemiological studies on Indonesia. II. Measurement of anti-arbovirus antibodies in sera from residents of Lombok, south Sumatera and West Java, in comparison with results concerning sera from residents of Japanese main islands. Kobe J Med Sci. 1970;16(4):215–34.
Olson JG, Ksiazek TG, Gubler DJ, Lubis SI, Simanjuntak G, Lee VH, Nalim S, Juslis K, See R. A survey for arboviral antibodies in sera of humans and animals in Lombok, Republic of Indonesia. Ann Trop Med Parasitol. 1983;77(2):131–7.
Hotta S. Sero-epidemiology of arboviruses in Indonesia. Kobe J Med Sci. 1972;18(4):211–4.
Hotta S, Aoki H, Samoto S, Yasui T, Noerjasin B. Virologic-epidemiological studies on Indonesia. 3. HI antibodies against selected arboviruses (groups A and B) in human and animal sera collected in Surabaja, East Java, in 1968. Kobe J Med Sci. 1970;16(4):235–50.
Iswara A. Soeharto, Biroum: anti-arbovirus HI antibodies in sera from suspected hemorrhagic fever patients found in Surabaja, Indonesia, in 1968-69. Kobe J Med Sci. 1970;16(4):211–4.
Miura Y, Inaba Y, Tsuda T, Tokuhisa S, Sato K, Akashi H, Matumoto M. A survey of antibodies to arthropod-borne viruses in Indonesian cattle. Jpn J Vet Sci. 1982;44(6):857–63.
Verlinde JD, Molron JH, Wyler R. Antibodies that neutralize neurotropic arthropod-borne viruses in residents of curacao and Indonesia. Doc Med Geogr Trop. 1955;7(1):94–6.
Vantongeren HAE, Vankammen A. Arthropod-borne virus-infections in Western new-Guinea (Irian-Jaya) - serological retrospection. Trop Geogr Med. 1978;30(4):413–27.
Smithburn KC. Neutralizing antibodies against arthropod-borne viruses in the sera of long time residents of Malaya and Borneo. Am J Hyg. 1954;59(2):157–63.
Hotta S, Aoki H, Yasui T, Samoto S. Virologic-epidemiological studies on Indonesia survey of anti-arboviral antibodies in sera from residents of Lombok, Sumatera and Djawa. Kobe J Med Sci. 1967;13(3):221–34.
We would like to thank David Harley (Mater Research Institute, The University of Queensland) and Patrick Gérardin (Centre Hospitalier Universitaire de la Réunion) for suggestions that improved the manuscript.
HH is a recipient of an Australia Awards Scholarship (AAS) from the Department of Foreign Affairs and Trade (ST000DMX2), Scholarship for International Research Fees (SIRF) from The University of Western Australia (F51416) and NHMRC Centre of Research Excellence in Emerging Infectious Diseases (CREID) Scholarship. AM is supported by West Australian Mosquito-Borne Diseases Scholarship from the Department of Health of Western Australia, Australia. The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.