Volume 22, Number 5—May 2016
Molecular Characterization of Chikungunya Virus, Philippines, 2011–2013
Chikungunya fever is a mosquitoborne infection that causes large outbreaks mainly in tropical and subtropical countries. The causative agent is chikungunya virus (CHIKV), an enveloped, single-stranded positive-sense RNA virus (family Togaviridae, genus Alphavirus). Phylogenetic analysis of the E1 gene of CHIKV revealed 3 genotypes: West African, East/Central/South African (ECSA), and Asian (1). In 2005, large outbreaks occurred in the islands in the Indian Ocean and India that were caused by the Indian Ocean lineage (IOL) virus, which newly emerged from the ECSA genotype (2). More recently, the emergence and potential spread of ECSA and Asian genotypes in the Americas have been major public health concerns in the region. (3).
In the Philippines, CHIKV was first isolated in 1965 (4). Since then, sporadic cases have been reported, including those among US Peace Corps volunteers stationed in the islands of Mindanao, Cebu, and Masbate in 1986 (5), and a local community outbreak in Cavite, Luzon Island, was reported in 1996 (6). However, the first nationwide CHIKV outbreak was identified starting in 2011 (Philippines Department of Health, unpub. data). Previous studies of chikungunya fever in the Philippines focused on clinical and serologic analyses (7), and recent reports on the molecular surveillance of CHIKV in the country were limited and analyzed only samples collected in 2012 (8,9). We conducted genetic analysis to characterize recent CHIKV infections that caused a large nationwide outbreak in the Philippines during 2011–2013.
Serum samples were collected through the chikungunya fever surveillance under the Philippine Integrated Disease Surveillance and Response of the Department of Health Epidemiology Bureau from different provinces in the Philippines during 2011–2013. Samples were collected from patients suspected to have chikungunya fever manifesting with symptoms such as fever, rash, and arthralgia. Samples were sent to the Research Institute for Tropical Medicine, which serves as the National Reference Laboratory for Dengue and Other Arboviruses. We screened serum samples were screened by using CHIKV IgM-capture ELISA (NovaLisa, NovaTec Immundiagnostica GmbH, Dietzenbach, Germany). We extracted viral RNA using the QIAamp Viral RNA kit (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions and amplified the partial E1 gene using 1-step reverse transcription PCR followed by direct Sanger sequencing with primers as previously described (10,11). We conducted phylogenetic analysis using the maximum-likelihood method as implemented in MEGA 6 software (http://www.megasoftware.net/). Molecular clock analysis and Bayes factor calculation were performed using BEAST software 1.8.0 (http://beast.bio.ed.ac.uk/) to select the best migration event model of CHIKV among countries. Bayes factor analysis was used to test phylogeographic hypothesis whether posterior migration rate between locations in whole evolutionary history was significantly higher than the expected prior migration rate, assuming truncated Poisson probability (12).
A total of 5,729 serum samples were collected from persons suspected to have chikungunya within 5 days after symptom onset. Fever, rash, and arthralgia were the most common symptoms (53%, 47%, and 34% of patients, respectively). Of the 5,729 serum samples, 2,891 were IgM positive by ELISA. We conducted reverse transcription PCR on 382 representative samples among IgM-negative patients in accordance with the CHIKV outbreak surveillance strategy of the Philippines, of which 131 samples tested positive. Partial E1 gene sequence (733 nt) was obtained from 31 samples. Sequences were submitted to GenBank (accession nos. LC064714–LC064744).
Phylogenetic analysis identified 28 Asian genotype viruses and 3 ECSA genotype viruses (Figure 1). Sequence analysis of E1 gene showed >99% nt similarity among 28 Asian genotype viruses (data not shown). Bayes factor analysis for the migration events also showed that migration might have occurred between Indonesia and Philippines with high probability (Table). CHIKV-positive patients were found only in Mindanao Island in 2011 but were detected in other parts of the country in 2012 and 2013. All CHIKV detected in parts of the country other than Mindanao were identified as Asian genotype.
The 3 ECSA genotype viruses were detected in the Davao Del Sur and Davao Oriental in Mindanao Island in 2012 and 2013, and all were clustered into IOL (Figure 1,2). This finding indicates that the 2 genotypes co-circulated in the island during the outbreak. ECSA genotype viruses in the Philippines were closely related to the Indonesian viruses (GenBank accession no. KC862329) detected in 2010. Sequencing analysis showed that all Philippine ECSA viruses possess the alanine to valine substitution (A226V) in the E1 gene (Technical Appendix[PDF - 280 KB - 6 pages] Tables 1, 2).
The recent large CHIKV outbreaks in other Asian countries were caused mainly by the IOL of ECSA genotype (2). However, the 2011–2013 outbreak in the Philippines was caused mainly by the Asian genotype; the reason that this large outbreak was caused by this genotype is unknown. Previous reports have confirmed the reemergence of Asian genotype viruses in the Philippines (8,9). In this study, phylogenetic analysis and Bayes factor calculation showed that the Philippines viruses were closely related to Indonesian viruses, which might explain why the outbreak started in southern Mindanao, near Indonesia. Although sequence data of CHIKV in the database are not enough to identify the exact origin of the virus, we tried to estimate when CHIKV was introduced into the Philippines with molecular clock analysis using the dataset of Asian genotype viruses (Technical Appendix[PDF - 280 KB - 6 pages] Figures 1, 2). The results suggested that circulating Asian genotype viruses in the Philippines were introduced from Indonesia before 2010.
We also detected IOL of ECSA genotype, which possesses the E1-A226V substitution, and its co-circulation with the Asian genotype in the Philippines. However, the geographic distribution of ECSA genotype in 2012 and 2013 was limited to southern provinces in Mindanao. After the large outbreak in the Indian Ocean region in 2005, IOL of ECSA genotype rapidly spread to Asian countries and then co-circulated with the endemic Asian viruses, then eventually became the predominant genotype (13). However, this circulation pattern differs from what we observed in the Philippines. Until the early 2010s, most of the viruses circulating in the Philippines and Indonesia were still Asian genotype, and ECSA genotype viruses had been reported in West Kalimantan, Indonesia, in 2011 (14) and in the Philippines in 2012 several years after the large outbreak in the Indian Ocean region (2). The movement of persons near the border of these countries might play a key role in CHIKV transmission.
As part of the national vector control program, surveillance of Aedes aegypti and Ae. albopictus mosquitoes has been conducted in several areas in the Philippines. A previous report showed that the proportion of these 2 mosquito species was almost the same in Metro Manila (15). If the proportion of Ae. albopictus mosquitoes increases, the ESCA genotype virus with A226V mutation could spread more rapidly in the country. Thus, monitoring the spread of ECSA genotype viruses and the proportion of the Aedes mosquitoes in the Philippines is important.
We have demonstrated that the Asian genotype CHIKV, which is closely related to the Indonesian viruses, was identified in Mindanao in 2011 and spread to other regions in 2012 and 2013. Like the Asian genotype, ECSA genotype virus was first detected in Mindanao in 2012. Mindanao might play a key role for the introduction of the CHIKV into the Philippines. Further monitoring is necessary to define the transmission pattern of CHIKV, including cross-border transmission.
Ms. Sy is a senior science research specialist and technical supervisor of the National Reference Laboratory for Dengue and Other Arboviruses, Virology Department, Research Institute for Tropical Medicine, Philippines. Her research interests include epidemiology and diagnosis of arboviral infections.
We thank the staff of the Epidemiology Bureau, Department of Health, Philippines, for the epidemiologic information about the samples. We also thank the staff of the Department of Virology and the Molecular Biology Laboratory, Research Institute for Tropical Medicine, Department of Health, Philippines, and staff of Tohoku–Research Institute for Tropical Medicine Collaborating Research Center on Emerging and Re-emerging Infectious Diseases for their excellent laboratory work.
This study was funded by the Epidemiology Bureau and the Disease Prevention and Control Bureau of the Department of Health, Philippines, and the Japan Initiative for Global Research Network for Infectious Diseases (J-GRID).
- 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:471–9.
- Schuffenecker I, Iteman I, Michault A, Murri S, Frangeul L, Vaney MC, Genome microevolution of chikungunya viruses causing the Indian Ocean outbreak. PLoS Med. 2006;3:e263.
- Nunes MR, Faria NR, de Vasconcelos JM, Golding N, Kraemer MU, de Oliveira LF, Emergence and potential for spread of Chikungunya virus in Brazil.BMC Med. 2015;13:102.
- Campo LE, San Juan A, Cenabre LC, Almagro EF. Isolation of chikungunya virus in the Philippines. Acta Med Philipp. 1969;5:152–5.
- Centers for Disease Control (CDC). Chikungunya fever among U.S. Peace Corps volunteers—Republic of the Philippines. MMWR Morb Mortal Wkly Rep. 1986;35:573–4 .
- Retuya TJA, Ting DL, Dacula BD, Lanada JM, Roque VG, Hugo CT, Chikungunya fever outbreak in an agricultural village in Indang, Cavite, Philippines.Philipp J Microbiol Infect Dis. 1998;27:93–6.
- Capeding MR, Chua MN, Hadinegoro SR, Hussain II, Nallusamy R, Pitisuttithum P, Dengue and other common causes of acute febrile illness in Asia: an active surveillance study in children. PLoS Negl Trop Dis. 2013;7:e2331.
- Yoon IK, Alera MT, Lago CB, Tac-An IA, Villa D, Fernandez S, 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:e0003764.
- Tan KK, Sy AKD, Tandoc AO, Khoo JJ, Sulaiman S, Chang LY, Independent emergence of the cosmopolitan Asian chikungunya virus, Philippines 2012. Sci Rep. 2015;5:12279.
- Hasebe F, Parquet MC, Pandey BD, Mathenge EG, Morita K, Balasubramaniam V, Combined detection and genotyping of Chikungunya virus by a specific reverse transcription-polymerase chain reaction. J Med Virol. 2002;67:370–4.
- Arias-Goeta C, Mousson L, Rougeon F, Failloux AB. Dissemination and transmission of the E1-226V variant of chikungunya virus in Aedes albopictusare controlled at the midgut barrier level. PLoS ONE. 2013;8:e57548.
- Lemey P, Rambaut A, Drummond AJ, Suchard MA. Bayesian phylogeography finds its roots. PLOS Comput Biol. 2009;5:e1000520.
- Tsetsarkin KA, Chen R, Sherman MB, Weaver SC. Chikungunya virus: evolution and genetic determinants of emergence. Curr Opin Virol.2011;1:310–7 .
- Maha MS, Susilarini NK, Hariastuti NI. Subangkit. Chikungunya virus mutation, Indonesia, 2011. Emerg Infect Dis. 2015;21:379–81.
- Duncombe J, Espino F, Marollano K, Velazco A, Ritchie SA, Hu WB, Characterising the spatial dynamics of sympatric Aedes aegypti and Aedes albopictus populations in the Philippines. Geospat Health. 2013;8:255–65.
Suggested citation for this article: Sy AK, Saito-Obata M, Medado IA, Tohma K, Dapat C, Segubre-Mercado E, et al. Molecular characterization of chikungunya virus, Philippines, 2011–2013. Emerg Infect Dis. 2016 May [date cited]. http://dx.doi.org/10.3201/eid2205.151268
1These authors contributed equally to this article.