Figures
Abstract
The first documented chikungunya virus (CHIKV) outbreak in Lao People’s Democratic Republic (Lao PDR) occurred in 2012–2013. Since then, several imported and a few autochthonous cases were identified by the national arbovirus surveillance network. The present study aimed to summarize the main genetic features of the CHIKV strains detected in Lao PDR between 2014 and 2020. Samples from Lao patients presenting symptoms compatible with a CHIKV infection were centralized in Vientiane Capital city for real-time RT-PCR screening. Molecular epidemiology was performed by sequencing the E2-6K-E1 region. From 2014 to 2020, two Asian lineage isolates (e.g. French Polynesia; Indonesia), one ECSA-IOL lineage isolate (e.g. Thailand) and one unclassified (e.g. Myanmar) were imported in Vientiane Capital city. Sequences from the autochthonous cases recorded in the Central and Southern parts of the country between July and September 2020 belonged to the ECSA-IOL lineage and clustered with CHIKV strains recently detected in neighboring countries. These results demonstrate the multiple CHIKV introductions in Lao PDR since 2014 and provide evidence for sporadic and time-limited circulation of CHIKV in the country. Even if the circulation of CHIKV seems to be geographically and temporally limited in Lao PDR, the development of international tourism and trade may cause future outbreaks of CHIKV in the country and at the regional level.
Citation: Calvez E, Bounmany P, Somlor S, Xaybounsou T, Viengphouthong S, Keosenhom S, et al. (2022) Multiple chikungunya virus introductions in Lao PDR from 2014 to 2020. PLoS ONE 17(7): e0271439. https://doi.org/10.1371/journal.pone.0271439
Editor: Pierre Roques, CEA, FRANCE
Received: October 21, 2021; Accepted: June 30, 2022; Published: July 15, 2022
Copyright: © 2022 Calvez 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 article and its Supporting Information files.
Funding: MG received funding from the Agence Française de Développement grant n° CZZ 2146 01 (Ecomore2 project). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Chikungunya virus (CHIKV) is transmitted to humans by Aedes mosquito species [1]. With a single-stranded RNA positive sense genome, CHIKV belongs to the Togaviridae family, Alphavirus genus. During the acute phase of human infection, CHIKV infection is clinically characterized by sudden onset of high fever often accompanied by severe joint pain and other symptoms such as headache, maculopapular rash evolving to chronic poly arthralgia in up to one third of the patients [2]. Since 1952, CHIKV was recorded in more than 60 countries in Africa, Asia, Europe, the Southern Pacific region and the Americas [3–6]. The genetic polymorphism of the CHIKV genome allows to discriminate four lineages: the West African (WA), East/Central/South African (ECSA), ECSA-derived Indian Ocean (ECSA-IOL) and Asian [7]. In Southeast Asia, the circulation of the Asian lineage was described since the 1950s and an active co-circulation of the Asian, ECSA and ECSA-IOL lineages was highlighted since the beginning of the 2000’s [3, 5, 6].
Lao People’s Democratic Republic (Lao PDR) is a land-locked, least developed country located in the middle of the Indochinese peninsula that progressively opens its borders to international exchanges for more than a decade now. This central location increases the risk of arbovirus introductions as previously reported for dengue viruses (DENV) [8–13] and CHIKV [14, 15]. In 2012–2013, the first chikungunya outbreak was described in the Southern provinces of the country and phylogenetic studies revealed that the Lao strains belonged to the ECSA-IOL lineage with a Cambodian origin [14, 16]. Since then and until 2020, no outbreak of CHIKV has been recorded in the country. Here, we studied the phylogenic relationships of the few imported and autochthonous CHIKV cases identified between 2014 and 2020 in Lao PDR.
Methods
Ethic statement
Ethics approval was obtained from the Lao National Ethics Committee for Health Research of the Ministry of Health of Lao PDR (N°114/NECHR).
For the patient samples included in this study, participant or a parent or a legal guardian provided a written informed consent.
Human samples collection and screening
Samples of clinically-suspected CHIKV cases were collected by hospitals belonging to the arbovirus surveillance network implemented in 2012 in Lao PDR and coordinated by the Institut Pasteur du Laos [8, 9, 14]. Viral genomic RNA was extracted from human plasmas or culture supernatants using a Nucleospin DX or Nucleospin 96 core kit purification kit according to the manufacturer’s instructions (Macherey-Nagel) and samples were screened using real-time RT-PCR assays [17].
Samples from autochthonous cases from Bolikhamxay have been tested for the presence of anti-CHIKV antibodies. Detection of IgM and IgG was performed using SD Biosensor STANDARD Q Arbo Panel I (Z/D/C/Y) kit (SD Biosensor) and an in-house ELISA method as described elsewhere, respectively [18, 19].
Virus amplification
Our previous experience revealed that direct sequencing failed on samples (either plasmas or culture supernatants) displaying Ct values above 28 [14]. Therefore, virus isolation was attempted for samples with a low viral load, and for which a leftover of at least 150 μL of original plasma was available. Part of the remaining volume (between 100 and 200μL) was used to perform one or two passages on Vero E6 cells maintained in culture at 37°C with 5% CO2 in Medium 199 (Gibco™, Thermo Fisher Scientific, Waltham, MA, USA) supplemented with 2% fetal bovine serum (FBS; Gibco™, Thermo Fisher Scientific, Waltham, MA, USA). After five days of incubation, supernatants were collected and stored at -80°C before analysis.
E2-6K-E1 gene amplification
A DNA fragment encompassing the E2-6K-E1 region (2,772 bp) was prepared from RNA purified from plasma or culture supernatant by RT-PCR using primers set described elsewhere [14, 20] (S1 Data). Amplification products were sequenced using primers used for the RT-PCR amplification as well as internal primers (S1 Data) [9]. Briefly, after generating the first stand cDNA using a Maxima H Minus First Stand cDNA Synthesis kit (Thermo Scientific, Waltham, MA, USA), Phusion Flash High-Fidelity PCR Master Kit (New England Biolabs® Inc, Waltham, MA, USA) was used for the overlapping PCRs. Amplification products were then purified by ExoSAP-ITTM PCR Product Cleanup Reagent (Thermo Fisher Scientific, Waltham, MA, USA). Sequencing was performed on both strands using BigDye Terminator v3.1 Cycle sequencing kit (Applied Biosystem, Waltham, MA, USA) on a Genetic Analyzer 3500xL (Applied Biosystem, Waltham, MA, USA).
Phylogenetic analysis
Raw sequences were analyzed and edited in BioEdit version 7.0.5.3 software (Manchester, United Kingdom). Then, a multiple-sequence alignment was built using a panel of reference sequences retrieved from GenBank using the ClustalW program and the alignment was checked manually (S2 Data). The sequences were translated into amino acids, and both nucleotide and amino acid sequences were checked for irregularities. Pairwise sequence identity (at the nucleotide and amino acid levels) of the E2-6K-E1 coding sequences was calculated with MEGA version 7 (www.megasoftware.net) using uncorrected p-distances. A phylogenetic tree was inferred from the aligned nucleotide sequences by using a maximum likelihood phylogenetic approach [9, 14]. The best-fitted model of nucleotide substitution was determined using MEGA v7 under corrected Akaike information criteria (AICc). The Kimura-2 model was identified as the optimal model of nucleotide evolution. One thousand bootstrap replicates were generated.
Results
Epidemiological analysis of the patients infected with CHIKV
Suspected CHIKV cases, matching with a case definition of CHIKV infection, were identified in Lao PDR between 2014 and 2020, corresponding to either imported cases or to autochthonous cases in Bolikhamsay and Savannakhet provinces.
A total of 18 cases were confirmed in Lao PDR between 2014 and 2020 of which four were imported cases from French Polynesia, Indonesia, Myanmar and Thailand, respectively (S3 Data). All autochthonous cases were recorded between July and September 2020 in one district of the central province of Bolikhamxay and in another one of the Southern province of Savannakhet (Fig 1). Among the samples collected in Bolikhamxay (N = 32), four were positive for anti-CHIKV IgM only, ten for IgM/IgG and 11 for IgG. These 25 anti-CHIKV antibody positive samples include four of the eight confirmed cases by RT-PCR from Bolikhamxay. Given the low volume (less than 250μL) of the samples from Savannakhet and from the imported cases, molecular tests have been prioritized.
Data based on the analysis of the samples collected by the Institut Pasteur du Laos arbovirus surveillance network between 2012 and 2020. The provinces indicated in yellow correspond to those where at least one sample was found positive for chikungunya virus during the study period (2014–2020). The red triangle indicates the province where the chikungunya outbreak occurred in 2012–2013.
In Bolikhamxay, eight CHIKV positive samples were detected. The mean age of the positive cases was 23.25 years (range 10–39 years). Among them, seven samples were from men and one from a woman. Of the six confirmed cases from Savannakhet, three were men and the other three women. The mean age for the Savannakhet’ patients was 60.33 years (range 30–68 years).
Genetic diversity of sequences
In this series, E2-6K-E1 sequences were obtained from six samples from either plasmas or culture supernatants. These six samples had a low Ct value ≤ 25 (Table 1). We were unable to obtain the sequence from the imported case from Myanmar due to the low viral load (Ct = 38) and virus isolation attempts were unsuccessful.
The six strains exhibited between themselves from 92.6 to 99.9% nucleotide identity each other and from 96.1 to 99.9% amino acid identity over the 2,772 bp / 923 aa sequences. The three sequences from the 2020 autochthonous cases as well as the one from the imported case from Thailand displayed from 99.7 to 99.9% nucleotide identity and from 99.6 to 99.9% amino acid identity. They exhibited from 99.6 to 100% nucleotide identity and from 99.7 to 100% amino acid identity with CHIKV strains identified in different Asian countries (Thailand, Myanmar, China, Bangladesh and Malaysia) since 2017. They also showed respectively from 98.4 to 98.7% and from 98.8 to 99.1% nucleotide and amino acid identity with the 2012–2013 Lao strains from Champassak. In addition, the sequence from the imported case from French Polynesia was identical to the one (GenBank accession number KR559473) described in 2015. At last, the sequence from the case imported from Indonesia exhibited from 99.2 to 99.7% nucleotide identity and from 99.6 to 99.8% amino acid identity with sequences from Indonesia described from 2013 to 2018.
Looking at specific amino acid signatures, the three autochthonous sequences here described from Bolikhamxay and Savannakhet and the four ones identified in Champassak in 2013 showed a cumulative number of 30 amino acid substitutions over 923 positions as compared to the reference S27 strain (NP_690589.2): 20 common to all strains, seven specific to either the 2013 or the 2020 strains, one (F84S in E2) shared between the two sequences from Bolikhamxay with one of the three sequences from Champassak (LN901352), and two (S60N in 6K and A359G in E1) restricted to a single strain (LN901359 and MZ292731, respectively). Twenty polymorphic positions were located in E2 (G57K, I74M, G79E, F84S, T110S, N160T, A164T, L181M, S194G, G205S, I211T, K252Q, V264A, M267R, N273I, S299N, T312M, A344T, S375T and V386A), three in 6K (V8I, I54V and S60N) and eight in E1 (K211E, A226V, M269V, D284E, I317V, V322A and A359G). Focusing on the main genetic signatures for the E1 and E2 envelope glycoproteins, K211E and I317V in E1 were specific to the 2020 strains while A226V was specific to the 2013 strains. Regarding E2, T110S and K252Q were restricted to the 2013 strains while G205S and V264A were specific to the 2020 strains.
Phylogenetic analysis of CHIKV
On a phylogenetic viewpoint, the imported strains from French Polynesia and Indonesia belonged to the Asian lineage whereas the one imported from Thailand belonged to the ECSA-IOL lineage (Fig 2). Within the Asian lineage, the strain from Indonesia grouped into a cluster of sequences from Indonesia identified between 2013 and 2018 while the one imported from French Polynesia was associated with another strain from French Polynesia (KR559473) within a cluster of sequences from the Americas. In addition, the 2020 autochthonous Lao strains as well as the one imported from Thailand belonged to the ECSA-IOL lineage. They all fell within a cluster of Asian sequences from Thailand, China, Myanmar, Bangladesh and Malaysia identified between 2017 and 2020. Within the ECSA-IOL lineage, these strains grouped in a cluster independent from the former one gathering the Champassak 2013 isolates (Fig 2).
Only bootstrap values >70 are shown. The Lao strains sequenced in this study are indicated in red, the previously identified ones in blue.
Discussion
The circulation of CHIKV has been poorly documented in Lao PDR since the 2012–2013 outbreak [14, 15]. As CHIKV actively circulates around the world and in neighboring countries, such as in China, Thailand and Myanmar [21–24], this creates a threat for CHIKV re-introduction and spread in Lao PDR. Since 2014, the geographic diversity of the imported cases (French Polynesia; Indonesia; Thailand; Myanmar) and the occurrence of local transmissions support this assumption.
As for the 2012–2013 outbreak [14–16], the demographic profile of the patients seemed to demonstrate an active circulation of CHIKV among the general population in Lao PDR in 2020. However, the virus circulation appeared to be geographically and temporally limited. After the 2012–2013 outbreak in Champassak province, seroprevalence study revealed a medium to high IgM/IgG seroprevalence levels ranging from 33% to 94% in selected villages in the province [14]. The population in Lao PDR is young (half of the population is under 24 years old) (www.lsb.gov.la). Since the 1960s, only few studies investigated the presence of anti-CHIKV antibodies among the Lao population but in restricted areas, Champassak province and Vientiane capital [14, 25, 26]. In absence of a seroprevalence study at the country level this suggests a low immunity against CHIKV among the Lao population that could promote the introduction and emergence of CHIKV in the country.
The CHIKV sequences identified in 2020 in Bolikhamxay and Savannakhet provinces belonged to the ECSA-IOL lineage. Interestingly, they grouped into a different cluster independent from the ones identified during the 2012–2013 outbreak (98.6% of nucleotide identity) [14]. In addition, they clustered with recent sequences isolated in Thailand between 2019 and 2020 as well as in neighboring countries such as Myanmar, where CHIKV outbreaks were recorded, or China where imported cases from Myanmar and Thailand were detected [22, 27–29].
At the amino acid level, CHIKV sequences from autochthonous Lao cases identified in 2020 harbored common substitutions with those of 2012–13. However, others substitutions are specific to the 2012–13 strains from Champassak or to the 2020 strains from Bolikhamxay and Savannakhet. Three positions, one in E1 (A226V) and two in E2 (T110S and K252Q) were restricted to the 2012–13 strains while four others, two in E1 (K211E and I317V) and the other two (G205S and V264A) in E2, were specific to the 2020 strains. These specific substitutions, as well as the phylogenetic closeness of the 2020 autochthonous Lao strains with others identified in neighboring countries at the same period are in support of an independent reintroduction of the virus in 2020 rather than their evolution from the 2012–13 strains. These genetic links are thus in favor of an active circulation of the virus at a regional level. It is noteworthy that substitutions within the CHIKV genome can have an impact on its transmissibility by the vectors [20, 30, 31]. Indeed, as it has previously been demonstrated, the presence of a single mutation in the E1 and E2 region could influence the virus transmission by Aedes [32–37]. Further investigations on the impact of the mutations here described on virus transmission should be carried out with the main Aedes vectors in Lao PDR, Ae. aegypti and Ae. albopictus [9, 14, 38, 39].
Overall, the results obtained in this study highlight the re-introduction and a limited circulation of novel CHIKV strains in Lao PDR. Indeed, even if CHIKV sequences obtained in 2012–2013 and 2020 are closely related (98.4–98.7% of nucleotide identity) between each other, the presence of specific mutations within E2-6K-E1 region is correlated with the circulation of new variants in the country as already observed in several countries in Southeast Asia especially in Thailand [23, 28] and in China [40]. Furthermore, as described in this study, phylogenetic studies on CHIKV and DENV in Lao PDR demonstrated that introduction and spread of arboviruses in this country are correlated with their circulation at the regional level, particularly in countries neighboring such as Thailand, Cambodia, Myanmar and China [8–10, 12, 14]. The central geographic position of Lao PDR in Southeast Asia could promote this risk of arbovirus introduction from neighboring countries due to human mobility [41–43] and should be further investigated due to development of the international trade through, for example, the Lao-China rail connectivity.
Up until now, only the ECSA-IOL lineage has been involved in the autochthonous circulation of CHIKV in Lao PDR, which is intriguingly limited to the Central and Southern regions of the country while the imported cases were detected in Vientiane Capital. However, several factors associated with the virus (genetic, viremia), its vector(s) (density, lifespan, vector competence), the host (immunity), and the environment (temperature, humidity) can influence the transmission of arboviruses. These factors and their interactions have to be investigated to prevent and control future CHIKV outbreaks in Lao PDR.
Supporting information
S1 Data. List and position of primers used for RT-PCR and sequencing of the E2-6K-E1 region.
https://doi.org/10.1371/journal.pone.0271439.s001
(DOCX)
S2 Data. References of CHIKV E2-6K-E1 sequences from GenBank used in this study.
https://doi.org/10.1371/journal.pone.0271439.s002
(DOCX)
S3 Data. Information relative to the Lao chikungunya samples collected between 2014 and 2020 and found positive by RT-PCR.
https://doi.org/10.1371/journal.pone.0271439.s003
(DOCX)
Acknowledgments
The authors thank all the clinicians implicated in the arbovirus surveillance network for their assistance in samples’ collection.
References
- 1. Higgs S, Vanlandingham D. Chikungunya Virus and Its Mosquito Vectors. Vector-Borne and Zoonotic Diseases. 2015;15: 231–240. pmid:25674945
- 2. Schwartz O, Albert ML. Biology and pathogenesis of chikungunya virus. Nature Reviews Microbiology. 2010;8: 491–500. pmid:20551973
- 3. World Health Organization. Chikungunya update 15 September 2020. https://www.who.int/news-room/fact-sheets/detail/chikungunya (accessed 21 December 2020). 2020. Available: https://www.who.int/news-room/fact-sheets/detail/chikungunya
- 4. Robinson MC. An epidemic of virus disease in Southern Province, Tanganyika territory, in 1952–1953. Transactions of the Royal Society of Tropical Medicine and Hygiene. 1955;49: 28–32. pmid:14373834
- 5. Sharif N, Sarkar MK, Ferdous RN, Ahmed SN, Billah MdB, Talukder AA, et al. Molecular Epidemiology, Evolution and Reemergence of Chikungunya Virus in South Asia. Frontiers in Microbiology. 2021;12: 1411. pmid:34163459
- 6. Wimalasiri-Yapa BMCR, Stassen L, Huang X, Hafner LM, Hu W, Devine GJ, et al. Chikungunya virus in Asia–Pacific: a systematic review. Emerging Microbes & Infections. 2019;8: 70–79. pmid:30866761
- 7. Volk SM, Chen R, Tsetsarkin KA, Adams AP, Garcia TI, Sall AA, et al. Genome-scale phylogenetic analyses of chikungunya virus reveal independent emergences of recent epidemics and various evolutionary rates. Journal of Virology. 2010;84: 6497–504. pmid:20410280
- 8. Lao M, Caro V, Thiberge J-M, Bounmany P, Vongpayloth K, Buchy P, et al. Co-Circulation of Dengue Virus Type 3 Genotypes in Vientiane Capital, Lao PDR. Coffey LL, editor. PLoS ONE. 2014;9: e115569. pmid:25551768
- 9. Calvez E, Pommelet V, Somlor S, Pompon J, Viengphouthong S, Bounmany P, et al. Trends of the Dengue Serotype-4 Circulation with Epidemiological, Phylogenetic, and Entomological Insights in Lao PDR between 2015 and 2019. Pathogens. 2020;9: 728. pmid:32899416
- 10. Calvez E, Somlor S, Viengphouthong S, Balière C, Bounmany P, Keosenhom S, et al. Rapid genotyping protocol to improve dengue virus serotype 2 survey in Lao PDR. Roques P, editor. PLOS ONE. 2020;15: e0237384. pmid:32764809
- 11. Dubot-Pérès A, Vongphrachanh P, Denny J, Phetsouvanh R, Linthavong S, Sengkeopraseuth B, et al. An Epidemic of Dengue-1 in a Remote Village in Rural Laos. Morrison AC, editor. PLoS Neglected Tropical Diseases. 2013;7: e2360. pmid:23951379
- 12. Castonguay-Vanier J, Klitting R, Sengvilaipaseuth O, Piorkowski G, Baronti C, Sibounheuang B, et al. Molecular epidemiology of dengue viruses in three provinces of Lao PDR, 2006–2010. Aguilar PV, editor. PLOS Neglected Tropical Diseases. 2018;12: e0006203. pmid:29377886
- 13. Calvez E, Bounmany P, Balière C, Somlor S, Viengphouthong S, Xaybounsou T, et al. Using Background Sequencing Data to Anticipate DENV-1 Circulation in the Lao PDR. Microorganisms. 2021;9: 2263. pmid:34835389
- 14. Somlor S, Vongpayloth K, Diancourt L, Buchy P, Duong V, Phonekeo D, et al. Chikungunya virus emergence in the Lao PDR, 2012–2013. Ng LFP, editor. PLOS ONE. 2017;12: e0189879. pmid:29284012
- 15. Soulaphy C, Souliphone P, Phanthavong K, Phonekeo D, Phimmasine S, Bouaphan K, et al. Emergence of chikungunya in Moonlapamok and Khong Districts, Champassak Province, the Lao People’s Democratic Republic, May to September 2012. Western Pacific Surveillance and Response Journal. 2013;4: 46–50. pmid:23908956
- 16. Phommanivong V, Kanda S, Shimono T, Lamaningao P, Darcy AW, Mishima N, et al. Co-circulation of the dengue with chikungunya virus during the 2013 outbreak in the southern part of Lao PDR. Tropical Medicine and Health. 2016;44. pmid:27524929
- 17. Pastorino B, Bessaud M, Grandadam M, Murri S, Tolou HJ, Peyrefitte CN. Development of a TaqMan RT-PCR assay without RNA extraction step for the detection and quantification of African Chikungunya viruses. J Virol Methods. 2005;124: 65–71. pmid:15664052
- 18. Grandadam M, Caro V, Plumet S, Thiberge J-M, Souarès Y, Failloux A-B, et al. Chikungunya Virus, Southeastern France. Emerg Infect Dis. 2011;17: 910–913. pmid:21529410
- 19. Prat CM, Flusin O, Panella A, Tenebray B, Lanciotti R, Leparc-Goffart I. Evaluation of Commercially Available Serologic Diagnostic Tests for Chikungunya Virus. Emerg Infect Dis. 2014;20: 2129–2132. pmid:25418184
- 20. Schuffenecker I, Iteman I, Michault A, Murri S, Frangeul L, Vaney M-C, et al. Genome Microevolution of Chikungunya Viruses Causing the Indian Ocean Outbreak. PLOS Medicine. 2006;3: e263. pmid:16700631
- 21. Chinnawirotpisan P, Chusri S, Manasatienkij W, Joonlasak K, Huang AT, Poolpanichupatam Y, et al. Complete Coding Sequences of 22 East/Central/South African Genotype Chikungunya Virus Isolates from Thailand (2018 to 2019). Matthijnssens J, editor. Microbiology Resource Announcements. 2020;9. pmid:33060262
- 22. Kyaw AK, Tun MMN, Nabeshima T, Soe AM, Thida T, Aung TH, et al. Chikungunya Virus Infection in Blood Donors and Patients During Outbreak, Mandalay, Myanmar, 2019. Emerging Infectious Diseases. 2020;26: 2741–2745. pmid:33079056
- 23. Khongwichit S, Chansaenroj J, Thongmee T, Benjamanukul S, Wanlapakorn N, Chirathaworn C, et al. Large-scale outbreak of Chikungunya virus infection in Thailand, 2018–2019. Roques P, editor. PLoS ONE. 2021;16: e0247314. pmid:33690657
- 24. Puntasecca CJ, King CH, LaBeaud AD. Measuring the global burden of chikungunya and Zika viruses: A systematic review. Mostafa A, editor. PLoS Negl Trop Dis. 2021;15: e0009055. pmid:33661908
- 25. Kanamitsu M, Taniguchi K, Urasawa S, Ogata T, Wada Y, Wada Y, et al. Geographic distribution of arbovirus antibodies in indigenous human populations in the Indo-Australian Archipelago. The American Journal of Tropical Medicine and Hygiene. 1979;28: 351–363. pmid:453438
- 26. Halstead SB. Mosquito-borne haemorrhagic fevers of South and South-East Asia. Bull World Health Organ. 1966;35: 3–15. pmid:5297536
- 27. Díaz-Menéndez M, Esteban ET, Ujiie M, Calleri G, Rothe C, Malvy D, et al. Travel-associated chikungunya acquired in Myanmar in 2019. Eurosurveillance. 2020;25. pmid:31937394
- 28. Phadungsombat J, Imad H, Rahman M, Nakayama EE, Kludkleeb S, Ponam T, et al. A Novel Sub-Lineage of Chikungunya Virus East/Central/South African Genotype Indian Ocean Lineage Caused Sequential Outbreaks in Bangladesh and Thailand. Viruses. 2020;12: 1319. pmid:33213040
- 29. Feng Y, Xiao H, Li X, Lu Y, Li J, Zheng M, et al. Nonindigenous East/Central/South African genotype of chikungunya virus identified in febrile returning travellers in Yunnan, China. Journal of Infection. 2020;80: 469–496. pmid:31891726
- 30. Powers A. Chikungunya virus outbreak expansion and microevolutionary events affecting epidemiology and epidemic potential. RRTM. 2015; 11.
- 31. Sahadeo NSD, Allicock OM, De Salazar PM, Auguste AJ, Widen S, Olowokure B, et al. Understanding the evolution and spread of chikungunya virus in the Americas using complete genome sequences. Virus Evolution. 2017;3. pmid:28480053
- 32. Viglietta M, Bellone R, Blisnick AA, Failloux A-B. Vector Specificity of Arbovirus Transmission. Front Microbiol. 2021;12: 773211. pmid:34956136
- 33. Lambrechts L. Quantitative genetics of Aedes aegypti vector competence for dengue viruses: towards a new paradigm? Trends Parasitol. 2011;27: 111–4. pmid:21215699
- 34. Vega-Rúa A, Marconcini M, Madec Y, Manni M, Carraretto D, Gomulski LM, et al. Vector competence of Aedes albopictus populations for chikungunya virus is shaped by their demographic history. Commun Biol. 2020;3: 326. pmid:32581265
- 35. Coffey L, Failloux A-B, Weaver S. Chikungunya Virus–Vector Interactions. Viruses. 2014;6: 4628–4663. pmid:25421891
- 36. Chen R, Plante JA, Plante KS, Yun R, Shinde D, Liu J, et al. Lineage Divergence and Vector-Specific Adaptation Have Driven Chikungunya Virus onto Multiple Adaptive Landscapes. Duggal NK, editor. mBio. 2021;12: e02738–21. pmid:34749526
- 37. Tsetsarkin KA, Vanlandingham DL, McGee CE, Higgs S. A single mutation in chikungunya virus affects vector specificity and epidemic potential. PLOS Pathogens. 2007;3: e201. pmid:18069894
- 38. Marcombe S, Fustec B, Cattel J, Chonephetsarath S, Thammavong P, Phommavanh N, et al. Distribution of insecticide resistance and mechanisms involved in the arbovirus vector Aedes aegypti in Laos and implication for vector control. Barker CM, editor. PLOS Neglected Tropical Diseases. 2019;13: e0007852. pmid:31830027
- 39. Tangena J-AA, Marcombe S, Thammavong P, Chonephetsarath S, Somphong B, Sayteng K, et al. Bionomics and insecticide resistance of the arboviral vector Aedes albopictus in northern Lao PDR. Moreira LA, editor. PLOS ONE. 2018;13: e0206387. pmid:30359425
- 40. Qiu S, Guo J, Li P, Li P, Du X, Hao R, et al. Source-tracking of the Chinese Chikungunya viruses suggests that Indian subcontinent and Southeast Asia act as major hubs for the recent global spread of Chikungunya virus. Virol J. 2021;18: 203. pmid:34635129
- 41. Kiang MV, Santillana M, Chen JT, Onnela J-P, Krieger N, Engø-Monsen K, et al. Incorporating human mobility data improves forecasts of Dengue fever in Thailand. Sci Rep. 2021;11: 923. pmid:33441598
- 42. Lai S, Johansson MA, Yin W, Wardrop NA, van Panhuis WG, Wesolowski A, et al. Seasonal and interannual risks of dengue introduction from South-East Asia into China, 2005–2015. PLOS Neglected Tropical Diseases. 2018;12: e0006743. pmid:30412575
- 43. Wesolowski A, Qureshi T, Boni MF, Sundsøy PR, Johansson MA, Rasheed SB, et al. Impact of human mobility on the emergence of dengue epidemics in Pakistan. Proc Natl Acad Sci USA. 2015;112: 11887–11892. pmid:26351662