
==== Front
Virulence
Virulence
Virulence
2150-5594
2150-5608
Taylor & Francis

39263937
10.1080/21505594.2024.2401985
2401985
Version of Record
Review Article
Review Article
Host-adaptive mutations in Chikungunya virus genome
X. NING ET AL.
VIRULENCE
Ning Xinhang a
Xia Binghui a
Wang Jiaqi a
Gao Rong b
https://orcid.org/0009-0000-7897-4831
Ren Hao a
a Department of Microbiology, Faculty of Naval Medicine, Shanghai Key Laboratory of Medical Biodefense, Naval Medical University , Shanghai, People’s Republic of China
b Department of Respiratory Medicine, The People’s Liberation Army Joint Logistic Support Force 943 Hospital , Wuwei, Gansu, People’s Republic of China
CONTACT Hao Ren hren2020@163.com
12 9 2024
2024
12 9 2024
15 1 2401985Integra12 9 2024
Integra12 9 2024
09 5 2024
08 8 2024
31 8 2024
© 2024 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group.
2024
The Author(s)
https://creativecommons.org/licenses/by/4.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The terms on which this article has been published allow the posting of the Accepted Manuscript in a repository by the author(s) or with their consent.

ABSTRACT

Chikungunya virus (CHIKV) is the causative agent of chikungunya fever (CHIKF), and its primary vectors are the mosquitoes Aedes aegypti and Aedes albopictus. CHIKV was initially endemic to Africa but has spread globally in recent years and affected millions of people. According to a risk assessment by the World Health Organization, CHIKV has the potential seriously impact public health. A growing body of research suggests that mutations in the CHIKV gene that enhance viral fitness in the host are contributing to the expansion of the global CHIKF epidemic. In this article, we review the host-adapted gene mutations in CHIKV under natural evolution and laboratory transmission conditions, which can help improve our understanding of the adaptive evolution of CHIKV and provide a basis for monitoring and early warning of future CHIKV outbreaks.

KEYWORDS

Chikungunya virus
arbovirus
adaptive mutation
evolution
The author(s) reported there is no funding associated with the work featured in this article.
==== Body
pmcChikungunya virus (CHIKV) is an enveloped, single-stranded RNA virus of the genus Alphavirus in the family Togaviridae and is transmitted through the bite of the mosquitoes Aedes albopictus or Aedes aegypti [1]. Infection in humans by CHIKV can lead to chikungunya fever (CHIKF), characterized primarily by fever, rashes, and joint pain [2]. Other symptoms include headache, myalgia, nausea and vomiting, conjunctivitis, and meningitis. Its impacts on the central nervous system are the most harmful [3], and severe acute cases may progress to heart and multi-organ failure [4].

CHIKV was first isolated from a febrile patient in Tanzania in 1952 and has since become endemic only in parts of Africa and Southeast Asia [5,6]. A CHIKF outbreak in eastern Kenya in 2004 was followed by a CHIKV pandemic in the region surrounding the Indian Ocean [7], and a large outbreak on the French island of Réunion in 2005–2006, affecting hundreds of thousands of people [8,9]. The first local transmission of CHIKV occurred in the United States in 2013 and subsequently spread to Central and South America, the Caribbean, and beyond [10]. CHIKF is currently endemic in various regions of Asia, Africa, Europe, and the Americas. According to the World Health Organization, CHIKV has been detected in approximately 60 countries worldwide since 2004, affecting millions of people [11].

Changing host adaptation is a strategy through which viruses expand their host range and endemic area [12]. Genetic mutation of several or even only a single amino acid is one of the primary ways by which the host range of a virus can be expanded or by which the virus undergoes changes in pathogenicity. For example, the alanine-to-valine mutation at residue 188 of a non-structural protein of the Zika virus over the course of its evolution enhanced its infection of mosquitoes and replication in mice and is a cause of the accelerated endemicity of the Zika virus [13]. The S139N mutation in the prM protein, which appeared during a Zika virus epidemic in the Americas, resulted in a significant increase in the viral infectivity of neural progenitor cells in mice and humans, more severe microcephaly in mouse foetuses, and increased mortality in neonatal mice [14]. Other arboviruses such as dengue virus also expand their infected population and endemic area through changes in host adaptation [15]. In addition to vector-borne viruses, other RNA viruses such as human immunodeficiency virus, hepatitis C virus, and poliovirus have developed changes in host adaptation [16–18].

CHIKV genomic structure

CHIKV is a single-stranded positive-sense RNA virus with a diameter of approximately 70 nm and a total genome length of approximately 11.8 kb, which includes untranslated regions (UTR) at the 5' and 3' ends and two open reading frames, ORF1 and ORF2, separated by a non-coding region [19,20]. Stem-loop RNAs formed from conserved sequence elements in the 5’ UTR, 3’ UTR, and sub-genomic promoters in the non-coding region, regulate replication and transcription of the viral genome. ORF1 encodes non-structural proteins (nsP1, nsP2, nsP3, and nsP4), and ORF2 encodes structural proteins, namely, capsid protein (C), envelope glycoproteins (E), and 6K/Transframe (TF) protein (Figure 1 [21]). The four non-structural proteins make up the replication complex. nsP1 has transmethylase activity and is the main enzyme involved in mRNA capping [22]; it also anchors the viral replication enzyme complex to the lipid membranes [23]. nsP2 has RNA helicase activity as well as cysteine protease and nucleotide triphosphatase activity, and can inhibit the immune response of the host cell. nsP3 plays a role in ADP ribose hydrolysis and RNA binding, while it promotes viral infection by interacting with some host cell proteins [24]. nsP4 has RNA polymerase activity and catalyzes RNA synthesis.These proteins assemble into viral replication complexes to synthesize full-length negative-strand RNA intermediates. These intermediates concurrently serve as a template for the synthesis of subgenomic (26S) and genomic (49S) RNAs [25,26]. The capsid protein contains two structural domains that interact with RNA in addition to forming the nucleocapsid. The E1 glycoprotein and the E2glycoproteins are arranged in trimeric spikes, (E1-E2)3, of heterodimers in T = 4 icosahedral symmetry, and the glycoprotein spikes protrude from a host-derived lipid bilayer that surrounds the nucleocapsid core containing the RNAgenome [27]. The glycoproteins play a role in virus entry, egress of mature particles from infected cells, particle stability, and immunogenicity [28,29]. E1 is a class II fusion protein consisting of three structural domains (I, II, and III) that mediate viral membrane fusion via hydrophobic peptides. E2 consists of three structural domains (A, B, and C) that play a key role in binding to cell receptors and are the primary targets of neutralizing antibodies. E3 promotes the folding and activation of p62 (the E2 protein precursor) and stabilizes the spikes formed by E2 and E1. 6K plays roles in virus budding and membrane permeability [27]. Figure 1. Genomic organization of CHIKV.

CHIKV genotypes and vectors

As an RNA virus, CHIKV is prone to errors during replication, resulting in a high mutation rate, which may lead to the emergence of novel endemic variants [30]. With this rapid generation of genetic diversity during replication, viral populations consist not just of a single genotype but form mutant clusters of multiple genotypes. Currently, CHIKV is divided into three majorlineages: the West African lineage, the East/Central/South African (ECSA) lineage, and the Asian lineage [31]. Researchers identified CHIKV collected from Indian Ocean region outbreaks (2004–2006) that emerged independently from the ECSA lineage to form the Indian Ocean lineage(ECSA-IOL) [32]. Phylogenetic analysis of ECSA-IOL indicates that it can be further subdivided into four sub-lineages (sl1, sl2, sl3 and sl4) [33].

The geographic distribution of female Ae. aegypti and Ae. albopictus as CHIKV vectors is strongly associated with the endemic area of CHIKF. Aedes aegypti originated in Africa and was later introduced to the New World; it is now distributed mainly in tropical and subtropical areas such as northern Brazil and Southeast Asia [34]. After the 2004–2006 outbreak in the Indian Ocean, it was discovered that Ae. albopictus could also serve as a vector for CHIKV [35]. Aedes albopictus first appeared in Asia and then spread to islands in the Pacific and Indian Oceans [36]; Aedes albopictus can survive in cooler temperatures than Ae. Aegypti; therefore, it has a wider geographic distribution and greater risk of transmission [37]; currently, its extent reaches temperate regions such as southern Europe, northern China, southern Brazil, the northern United States, and Japan [38]. In addition, more convenient global transportation has further expanded the geographic distribution of both Aedes species and has become a major cause of the global spread of CHIKV.

Impact of adaptive mutations on vector and population susceptibility to CHIKV

The increased frequency of CHIKV outbreaks and the expansion in its transmission since 2004 have been associated with increased host adaptation due to viral mutations. Studies have shown that adaptive mutations in the natural evolution of CHIKV that alter vector susceptibility are responsible for the dissemination of CHIKV to non-endemic areas as well as large-scale local transmission [39]. The ability to distinguish between these variants will provide a better understanding of the evolutionary tendencies of viruses, explain how viruses develop host jumps [40], and provide a theoretical basis for improving the efficiency of virus control.

Adaptive changes arising from the natural evolution of CHIKV

E1-A226V mutation enhances Ae. albopictus susceptibility

A major CHIKF outbreak occurred on Réunion Island in 2005–2006, with 265,000 clinical cases (34% of the population) and 237 deaths reported. As Ae. aegypti was virtually absent from the outbreak area but Ae. albopictus was abundant and CHIKV RNA was detected in Ae. albopictus, it was concluded that Ae. albopictus mediated the spread of the outbreak [35]. Characterization of the CHIKV genome revealed a mutation from alanine to valine at E1 glycoprotein 226 (E1-A226V).Infection of Ae. albopictus with CHIKV/E1-A226V revealed that the mutant virus enhanced midgut infectivity and viral transmission to the salivary glands and to vertebrates in Ae. albopictus, confirming that the E1-A226V mutation resulted in increased viral adaptation in Ae. albopictus and gave it a selective advantage over Ae. aegypti in accelerating its transmission to uninfected populations [41,42]. This mutation also led to the introduction of CHIKV into new regions, with outbreaks in areas such as Italy and France, where Ae. aegypti is absent but Ae. albopictus is abundant, and the gradual development of the IOL.

A second mutation on a E1-A226V background improves the vector adaptation of CHIKV

In a recent study, researchers used two strains of CHIKV with and without the E1-A226V mutation and found no significant difference in their ability to infect Ae. albopictus. This finding suggests that E1-A226V is not the only factor contributing to the enhanced replication of CHIKV in Ae. albopictus and that there may be other mutations affecting the vector adaptation of CHIKV [43].

A comparison of CHIKV sequences from outbreaks at the same location (Kerala, India) in 2007 and 2009 by Tsetsarkin et al. revealed two new mutations in the background of the E1-A226V mutation, namely E2-K252Q and E2-L210Q, which appeared in Southeast Asia and northeastern India in subsequent years.Researchers infected Ae. albopictus with a E2-L210Q mutant strain and found that the mutation increased the infectivity of CHIKV in Ae. albopictus midgut cells 2.3–2.4 fold. This mutation represented a stronger effect than that caused by the E1-A226V mutation, which increased the efficiency of CHIKV transmission via Ae. albopictus 2 fold, but neither had a significant effect on Ae. aegypti adaptation nor led to adaptive changes in primate cells (Vero and 293 cell lines).The mutation did not inhibit virus transmission with Ae. aegypti as the vector but increased transmission via Ae. albopictus, resulting in increased CHIKV endemicity [44]. Further studies revealed that E2-L210Q and E2-K252Q double mutant viruses spread to the brain of Ae. albopictus at higher rates than any of the single mutants, suggesting that the double mutation can lead to further increases in CHIKV host adaptation [33]. However, there is no evidence that the 2009 sub-lineage replaced the 2007 sub-lineage; thus, a more plausible explanation is that the two mutations, E2-K252Q and E2-L210Q, appeared successively and independently in the two outbreaks. Each of these, increased the fitness of Ae. albopictus compared to the parental lineage with only the E1-A226V mutation (The structural visualization of mutant residues is shown in Figure 3d) [33].

Expression of E1:K211E and E2:V264A on a E1:226A background revealed higher fitness for Ae. aegypti

In another study, researchers analysing the sequences of CHIKV isolates endemic to New Delhi, Andhra Pradesh, Tamil Nadu, and West Bengal, India between 2009 and 2010 found two new mutations in the background of E1-226A, namely, E1-K211E and E2-V264A. In Ae. aegypti, introduction of single mutation (E1-K211E, E2-V264A) and double mutation (E1-K211E:E2-V264A)in the 226A virus, respectively, resulted in increased viral titres in the midgut, saliva, legs, and wings, all of which were clearly observed after infection. However, in Ae. albopictus, introduction of a single or double mutation did not lead to any significant differences in midgut, leg, wing, and saliva titres at either day post infection. The results demonstrated that introduction of E1-K211E or E2-V264A mutations in the background of 226A led to increased infectivity, dissemination, and transmission in Ae. aegypti (The structural visualization of mutant residues is shown in Figure 3e) [45,46].

nsP4-R82S enhances viral replication in human hosts

An analysis of the CHIKF strain from an outbreak in Malaysia in 2008–2009 by Jolene et al. revealed an arginine-to-serine substitution at amino acid 82 of nsP4 (nsP4-R82S). Phylogenetic analysis indicated a gradual increase in the frequency of this mutation since the start of the 2008 outbreak. The nsP4-R82S sub-lineage was initially present only in Malaysia and Singapore [47], but later it spread to other Southeast Asian countries including Indonesia, Thailand, China, Cambodia, and Myanmar until it replaced the wild-type strain as the main sub-lineage of CHIKV circulating in Southeast Asia in 2009. Aedes albopictus is the main vector of CHIKV transmission in Malaysia. Studies in Ae. albopictus C6/36 cells have shown that the R82S mutation has no significant effect on viral infection and transmission in Ae. albopictus but is more selective in human cell lines. This finding suggests that the mutant strain nsP4-82S causes widespread CHIKF transmission in Southeast Asia by enhancing the fitness of CHIKV in human hosts rather than Ae. albopictus, leading to more efficient human-mosquito transmission [48].

3’UTR repeats enhance the replication of CHIKV in mosquitoes

The extensive within-species diversity in the 3’UTR of CHIKV (Figure 2) is unique in this genus of mainly mosquito-borne viruses.In conjunction with studies of the evolutionary trajectory of CHIKV and the relative conservation and comparability of sequences across lineages, this diversity of 3’UTR sequences offers a valuable avenue for elucidating their evolutionary and functional significance. Stapleford et al. [49] identified the CHIKV Asian genotype evolved to duplicate a 177 nt repeat element in the 3’UTR of the Caribbean lineage and the Caribbean Lineage 3’UTRreplicates in mammalian cells in a manner similar to another previously reported Asian CHIKV lineage (NC-2011). However, this duplication enhanced replication of the Caribbean lineage in mosquito cells. On this basis, Merwaiss et al. [50] used an engineered variant of the Caribbean strain of CHIKV bearing a deletion of the first 500 nt of the 3’UTR to explore the effect of 3’UTR on the virus in mosquitoes. The positive effect of duplications has also been recapitulated by engineering deletion mutants in the backbone of the Caribbean virus. Experimental infection of mosquitoes resulted in delayed replication in the midgut and impaired dissemination. The results of competition assays indicated that the wild-type virus was able to effectively replace the mutant virus during systemic infection. The above results are further confirmed by another study [51]. A series of mutant viruses was engineered based on two wild-type (wt) CHIKV strains, Mal06 (MY002IMR/06/BP; GenBank Acc. No. EU703759.1), representing the Asian lineage, and SL07 (SL-CK1; GenBank Acc. No. HM045801.1), representing the ECSA lineage.The results showed that deletion of any part of the 3’UTR of the Asian Mal06 strain severely reduced its replication rate in mosquito cells. The findings of a study by Claudia et al. [52] also demonstrated a positive correlation between the number of 3’UTR repeats and the rate of viral replication in mosquito cells.Overall, the results of both studies indicate a negative impact of deletion of 3’UTR repeat regions on replication in mosquito cells in vitro. Figure 2. Lineage-specific structures of the CHIKV 3’UTR.

The 3’UTR structures, based on sequence alignment. Direct repeats are illustrated by different colored blocks, with each of the four colors representing a different homologous sequence region. Sequence gaps in the alignment are indicated by white blocks. The Caribbean strain contains a 177-nucleotide duplication of the 3’ end of 1+2a and complete 1+2b in 3’UTR.

Adaptive changes during the experimental propagation of CHIKV

E1-V80I:A129V double mutation on a E1-A226V background enhances viral adaptation in Ae. albopictus

Stapleford et al. [53] infected Ae. aegypti with a virus containing the E1-A226V mutation isolated from Cambodian patients infected with CHIKV in 2012. Ten days later, the E1-V80I:A129V double mutation was identified in the saliva of all infected mosquitoes. The authors engineered double-mutation V80I:A129V and compared the strains in terms of infection (mosquito midgut), dissemination (legs and wings), and transmission (saliva) in both mosquito hosts. At high-infection doses (106 plaque-forming units [pfu]/ml), infection by and dissemination of the double-mutant virus were significantly higher than those in the control group in Ae. albopictus, whereas they were similar to the control group in Ae. aegypti. The infectivity,dissemination, and transmission of the double-mutant virus in Ae. albopictus were significantly increased at a low infectious dose (103 pfu/ml) compared to controls, but this phenomenon was not observed in Ae. aegypti. Animal experiments in eight-day-old mice inoculated with a sublethal dose of the virus (102 pfu) revealed that the double mutation increased the virus titres in blood and muscle 100-fold. Fusion experiments on baby hamster kidney (BHK) cells and mouse embryonic fibroblasts (MEF) showed that the double-mutant virus had significantly higherfusion at the most relevant pH range (pH 5.4–5.8). The E1-A129V single mutant also exhibited a tendency for increased fusion. Incubation of the mutant viruses in cell-free suspensions at 28 °C (mosquito temperature) and 37 °C (mammalian temperature) for 48 h resulted in increased stability at both temperatures of both double-mutant strains. This suggests that the E1-V80I:A129V mutation has a significant adaptive advantage both in vitro and in vivo in mosquitoes and mammalian hosts, and may be a gene capable of emerging and replacing the parental A226V prevalent strain during natural mosquito-to-mammalian transmission (The structural visualization of mutant residues is shown in Figure 3c) [53].Figure 3. Structural visualization of structural proteins co-mutated residues.

Protein structural data were accessed via the Protein Data Bank (PDB) and analyzed using PyMOL version 2.3.3.(a-i) Crystal structure of the mature envelope glycoprotein complex (spontaneous cleavage) of Chikungunya virus. (PDB entry 3N41) with green for E1; cyan for E2; grey for E3. (a) E1-A226V, E2-R198Q, E3-S18F co-mutation [33]. (b) E1-A226V, E1-N20Y, E1-M88L co-mutation [55]. (c) E1-A226V, E1-V80I, E1-A129V co-mutation [53]. (d) E1-A226V, E2- L210Q, E2- K252Q co-mutation [33]. (e) E1-A226V, E1:K211E, E2-V264A co-mutation [45]. (f) E1-A226V, E1-V156A, E1-K211T co-mutation [56]. (g) E1-A226V, E2-G60D, E2-I211T co-mutation [58]. (h) E1-A226V, E2-K233E/Q co-mutation [33]. (i) Map of residue positions before mutation of structural proteins.

E1-M88L and E1-N20Y mutations on a E1-A226V background enhance viral infectivity

Studies have led to the conclusion that E1-V80 may be a key residue that affects the infectivity and transmission of CHIKV. Further studies found that the mutation E1-V80Q significantly attenuates viral virulence in mammalian and Ae. aegypti cells [54]. However, at 14 dpi of Ae. aegypti with a low dose (104 PFU/mL) of the mutant virus, virulence was found to be restored in viruses isolated from two mosquitoes. Sequencing revealed that the viruses not only retained the E1-V80Q mutation, but also contained two new mutations (E1-M88L and E1-N20Y). To verify the effects of these two mutations on viral infectivity, ammonium chloride was used to block viral dissemination by neutralizing endosomal and lysosomal pH. The E1-M88L mutant virus was found to infect mammalian BHK-21 cells more rapidly and with greater infectivity, whereas the E1-N20Y mutant virus was not significantly altered. Enhanced infectivity of both E1-M88L and E1-N20Y mutant viruses was observed 60 min after infection of Ae. albopictus C6/36 cells. Viral binding assays were performed to confirm whether increased infectivity affects viral binding to cells. The results showed that E1-M88L resulted in increased viral binding compared to the wild-type strain, whereas E1-N20Y had no effect. In addition, compared to E1-N20Y and E1-A226V single mutant viruses, E1-M88L mutant viruses have reduced cholesterol dependence and can enter cells with less cholesterol. In summary, the E1-M88L mutation allows the virus to bind more efficiently to and enter Aedes and mammalian cells, which in turn enhances viral infection of Ae. albopictus and mammals, whereas the E1-N20Y mutation specifically enhances the infectivity of only Ae. albopictus cells, thereby causing increased virulence and enhanced viral transmission (The structural visualization of mutant residues is shown in Figure 3b) [55].

E1-V156A and/or K211T mutations on a E1-A226V background enhance viral pathogenicity

To investigate the evolution of CHIKV during transmission from vectors to vertebrate hosts, laboratory Ae. aegypti mosquitoes were fed bloodmeal containing the virulent E1-A226V strain, single mosquitoes were isolated at 10 dpi, and the mosquitoes were introduced into cages with suckling mice for 5 days. An additional round of vector-borne transmission was carried out, from infected mouse to naïve mouse via mosquito. Saliva and midgut from all stages of mosquitoes and blood from mice were collected for deep sequencing, which revealed a new V156A mutation in the E1 glycoprotein [53]. These experiments were repeated, and the mutation was stable in both the mosquito and mouse populations [56]. Not only was the E1-V156A mutation detected in two samples collected during a chikungunya outbreak in Rio de Janeiro in 2016 but also was the E1-K211T mutation present in all samples [57].

Membrane fusion experiments in BHK-21 cells revealed that the fusion pH thresholds of the two single mutants V156A and K211T were similar to those of the parental viruses; however, the double mutant E1-V156A:K211T had a lower fusion pH threshold, suggesting that the double-mutant viruses can undergo fusion and enter cells more easily. Previous studies suggested that the E1-K211E mutation increases the adaptability of Ae. aegypti [45]. When Ae. aegypti was infected with E1-K211T and E1-V156 single-/double-mutant viruses, there were elevated titres of the E1-K211T single-mutant virus in the legs and wings, which may confirm that the mutation of residue 211 is more dominant for infection and facilitates viral transmission in Ae. aegypti. Animal experiments in mice infected with single- and double-mutant viruses revealed increased swelling at the site of infection and higher muscle and serum viral titres compared to the parental virus, demonstrating that the E1-V156A and/or K211T mutations resulted in increased viral pathogenicity (The structural visualization of mutant residues is shown in Figure 3f) [56].

E2-G60D or E2-I211T on a E1-A226V background enhances ae. Aegypti infectivity

The adaptability of viral populations may not depend only on a single amino acid change; multiple mutations may accumulate during transmission and result in cumulative adaptive effects.Konstantin et al. introduced the E1-A226V mutation into three CHIKV strains (Ag41855, LR2006 OPY1 and 37,997) and observed changes in their infectivity. The researchers found that the Ag41855 strain had weaker infectivity of Ae. albopictus. Sequence analysis revealed the presence of a unique eight-site mutation in strain Ag41855 compared to the other two strains (Table 1). Preliminary experiments revealed that the Ag41855 strain containing the E2-G60D, E2-I211T and E2-V162A mutations did not show reduced virulence in BHK-21 cells.To further investigate the effects of this mutation on the strain, the E2-V162A mutation was introduced into Ag41855. No significant difference in the infectivity of the virus compared to the parental virus was found. In contrast, the introduction of E2-G60D into Ag41855 resulted in enhanced infectivity of the attenuated strain against Ae. albopictus, regardless of whether it was expressed together with E1-A226V. In addition, the E2-G60D mutation alone enhanced the infectivity of CHIKV against Ae. aegypti to a level similar to that of the other two strains.When E2-I211T was introduced into an attenuated strain, the mutation was found to enhance the infectivity of the virus against Ae. albopictus but had no effect on Ae. aegypti.However, this change only occurred when the mutation was co-expressed with the E1-A226V mutation, suggesting that the E2–211 residue acts to specifically regulate the infectivity of CHIKV in Ae. albopictus. Phylogenetic analysis showed that the E2-G60D mutation is primarily limited to laboratory strains, but the E2-I211T substitution may be responsible for the emergence of the 1976 South Africa strain, the 2005–2007 Indian Ocean strain, and the 2007 Gabon strain, thus warranting further research (The structural visualization of mutant residues is shown in Figure 3g) [58].Table 1. Genetic differences between Ag41855, LR2006 OPY1, and 37,997 strains of CHIKV.

 	Amino acid site	37997	LR2006 OPY1	Ag41855	
E1	19	V	V	I	
377	A	A	T	
E2	60	D	D	G	
162	A	A	V	
211	T	T	1	
nsP3	328	Q	Q	P	
358	S	S	P	
461	P	P	L	

E2-K233Q/E on a E1-A226V background enhances adaptation of CHIKV in Ae. albopictus

Introduction of E2-K233Q in the background of E1-226 V increased the efficiency of CHIKV infection of Ae. albopictus 5 fold. Further testing of the effect of the E2-K233Q mutation on viral infection and replication in the midgut of Ae. albopictus revealed that the relative amount of mutant virus RNA in the midgut was increased 3.6 fold at 1 dpi, progressively increasing to 8.9 fold at 3 dpi, indicating that the E2-K233Q mutation plays a role in the initial stages of infection and/or replication in the midgut [33]. In the study by Tsetsarkin et al., the GenBank database was searched for the E2-233Q mutation, and although it was not found, a mutant sequence with E2-233E was found instead. The introduction of E2-K233E into the E1-226 V background resulted in a 5.7-fold increase in the relative amount of viral RNA in the midgut at 1 dpi. Mosquitoes that ingested bloodmeal containing the double-mutant virus exhibited disseminated viral infection not only in the head but also in the mosquito homogenate, suggesting that E2-K233E increases the adaptation of CHIKV in other organs/tissues of the mosquito in addition to its effects on midgut infection (The structural visualization of mutant residues is shown in Figure 3h) [33].

Two double mutations on a E1-A226V background enhance adaptation of CHIKV in Ae. albopictus

The double mutation that emerged in the ECSA-IOL sl1 sub-lineage (nsP2-L539S/E2-K252Q) was introduced into the E1-226 V mutant strain. Measurements of viral infection in mosquito heads showed that the experimental viruses were significantly more efficient at disseminating infections via Ae. albopictus compared to the control group, in which only the E1-A226V mutation was present. In addition, the double-mutant virus resulted in a 7.8-fold average RNA increase in the Ae. albopictus midgut 10 days post infection (dpi). To further validate these results, the sl1 virus strain Thailand-2009_CK11/53 containing the characteristic nsP2-L539S/E2-K252Q substitution isolated from febrile patients was analysed. The results showed that viruses containing the double mutation had 5.8-fold greater transmission efficiency and that the Ae. albopictus viral titres at 10 dpi were higher than that of those control group. These results confirm that nsP2-L539S/E2-K252Q increases viral adaptation in Ae. albopictus infection. Introduction of a non-synonymous mutation in sl3 (E2-R198Q/E3-S18F) into E1-226 V increased the transmission ability of the experimental group of Ae. albopictus11.5 fold compared to the E1-226 V control. In addition, the relative amount of RNA in the midgut of mosquitoes infected with the double-mutant virus was increased 10.2-fold compared to the control group. These results indicate that the sl3 strain achieved further adaptation in Ae. albopictus after acquiring the two synergistic mutations E3-S18F and E2-R198Q (The structural visualization of mutant residues is shown in Figure 3a) [33].

Deletion of the 3’ UTR and E2-K200R enhances CHIKV virulence and transmission

CHIKV was isolated from the serum of mice persistently infected with CHIKV for 28 days and then inoculated into the left hindfoot pad of wild-type mice. These mice developed signs of more severe musculoskeletal lesions, suggesting that this virus is more virulent than its parent virus. Sequencing of this virus revealed three mutations: a substitution at residue 200 of the E2 protein (E2-K200R), a synonymous 10506A→G mutation in the E1 gene, and a 44-nt (11,921 –11,964) deletion in the 3’ UTR (Δ3’UTR). The E2-K200R mutation was also found in a Thai case in 1995, whereas the 44-nt deletion is highly conserved and has only been observed in the Asian lineage [59].

To confirm that these mutations were responsible for increased virulence, a E2-K200R:Δ3’ UTR double mutant virus was constructed and inoculated into the footpad of mice. The increased degree of weight loss and more severe musculoskeletal symptoms in mice infected with the double-mutant virus compared to controls confirmed that the mutations were responsible for the damage to musculoskeletal tissues distant from the inoculation site and aggravated systemic disease. Next, the transmissibility of the mutated virus was assessed, and the mutation was found to enhance the early transmission of CHIKV from the inoculation site to other sites. In-vitro experiments revealed that the amount of the virus detected in the experimental and control groups was similar in mouse fibroblasts, differentiated C2C12 mouse cells, and C6/36 mosquito cells, suggesting that the mutation had negligible effect on viral replication and mainly affected viral transmission in mice. Viruses with the E2-K200R or Δ3’UTR mutations were constructed, and it was found that the E2-K200R single mutation also increased viral transmissibility and pathogenicity, but the 44-nt deletion in the 3’ UTR had to be present in conjunction with the E2-K200R mutation in order to enhance viral pathogenicity [59].

Non-structural protein mutations on a E1-A226V background enhance viral adaptation in Ae. albopictus and vertebrates

CHIKV was sequenced after 19 blind passages in difluoromethylornithine (DFMO)-treated cells, and three mutations were identified, with two in nsP1 (G230R and V326M) and one in which the stop codon between nsP3 and nsP4 was mutated to arginine (*524 R). Mutant viruses were constructed and used to infect Vero-E6 cells, after which plaque formation experiments were performed. Viruses with the nsP1 mutation resulted in smaller plaques, and viruses with the nsP1/*524 R triple mutation resulted in larger plaques. Viral infection experiments in Ae. albopictus C6/36 cells revealed that each of the single mutant viruses (G230R, V326M, and *524 R) exhibited severe defects in replication, whereas the viral titres of the triple mutants were similar to wild-type levels. To investigate their adaptation, each single-mutant virus was mixed with wild-type CHIKV at a 1:1 ratio and used to infect BHK-21 or C6/36 cells for 24 h. Sequencing of the target mutations after five passages revealed that the single mutant was lost and that wild-type and triple-mutant viruses with the *524 R mutation were present in BHK-21 cells, suggesting that the single-mutant virus was more poorly adapted compared to the wild-type virus, whereas the viruses with the *524 R mutation were better adapted. After observing that the mutant viruses exhibited an adaptive advantage in vitro, Ae. albopictus mosquitoes and zebrafish, a vertebrate infection model of CHIKV, were infected with mutant CHIKV. The results showed that single-mutant viruses had decreased titres in mosquitoes compared to the wild-type, whereas the titres of the triple-mutant virus were significantly higher. These results indicate that single-mutant CHIKV exhibits growth defects in mosquitoes, but that triple-mutant viruses have enhanced replication and adaptation in mammals and Ae. albopictus both in vitro and in vivo [60].

The nsP1:R171Q mutation on a background E1-A226V was identified in virus isolates from CHIKF outbreaks in Comoros in 2005, in Sri Lanka in 2006, in India in 2008–2013, and in the Caribbean in 2014. Although the function of this mutation has not been elucidated, there is already evidence that the nsP1:171Q mutant exhibits enhanced pathogenicity in primate cells [61].

Other potentially adaptive mutations arising through evolution

Partial sequencing of the E1 region of CHIKV during the outbreak in the southern Indian state of Karnataka between November 2015 and November 2016 and alignment with other ECSA strains from India revealed six non-synonymous mutations: K211E, M269V, D284E, V322A, I317V, and V220I [62]. The K211E mutation facilitates enhanced host adaptation, as discussed earlier. The M269V and D284E mutations, which were reported in the 2005–2006 outbreaks in India, have been retained in currently prevalent strains, suggesting that they may be associated with host adaptation, but their specific effects have not been experimentally confirmed. The I317V mutation, which was previously found only in north-central India and Bangladesh, was first continuously transmitted in southern India, and deletion of the I317V mutation was observed in the 2020–2021 outbreak in Pune, India. Whether the I317V mutation affects CHIKV transmission is unknown, but absence of the I317V mutation may limit CHIKV transmission and may be a reason for why the number of cases in the 2020–2021 outbreak in Pune, India was lower than the number of cases in 2016 [63].

Comparing the sequences obtained from cases during the 2014–2016 CHIKF outbreak in Mexico with those sampled at the beginning of the outbreak, 14, 7, and 4 non-synonymous mutations were found in nsP3, E2, and E1, respectively, compared to that in the ECSA genotype, and four and one non-synonymous mutations were found in nsP3 and E2, respectively, compared to that in the Asian genotype. Three of these non-synonymous mutations (E2-S248F, nsP3-A437T, and nsP3-L451F) were also found in outbreaks in the Caribbean and the Philippines in 2012. It has been demonstrated that E2-S248L favours viral transmission in Ae. aegypti, whereas further studies are needed to determine whether the E2-S248F mutation affects host adaptation. The nsP3-A437T and nsP3-L451F mutations have not been studied in detail; however, the presence of the A437V mutation, at the same locus as nsP3-A437T, has been reported in previous ECSA samples, suggesting that the amino acid change at this locus may be similar to that of the E2-S248F mutation discussed earlier and affect the adaptation of Ae. aegypti and Ae. albopictus [64].

Concluding remarks

As an RNA virus, the adaptive mutations that occur during CHIKV replication impact the susceptibility of vectors and humans and are a major factor in the increasing global spread of CHIKF (Table 2). To date, the ECSA genotype has dominated, been widely prevalent, and evolving globally [39]. The 2004–2006 outbreak of CHIKF in the Indian Ocean region was identified as being caused by a newly emerged branch of the ECSA lineage, known as ECSA-IOL [65]. ECSA-IOL is prominently characterized by the presence of the E1-A226V mutation, which confers a selective advantage to the virus in Ae. albopictus [40], resulting in the increasing global spread of IOL sub-lineages and leading to over six million infections in subsequent years [35].Table 2. Chikungunya fever outbreaks and viral gene mutations.

Genotype	Year of outbreak	Outbreak area	Main vector	Mutation	
ECSA	1953	Tanzania [66]	Ae. aegypti	 	
1976	South Africa [66]	Ae. aegypti	 	
1982	Central Africa, Uganda[66]	Ae. aegypti	 	
2000	India [67]	Ae. aegypti	nsP1:V326M, L507R	
nsP2:H374Y	
nsP3:Q328P, S358P, I376T, C435R, V438A, M449T, S462N	
nsP4:Q500L	
E1:V19I, M269V, A377T	
E2:A162V, A164T, V318M, I377V	
E3:I23T	
C:K63R	
2002	Equatorial Guinea [68]	Ae. aegypti	 	
2004	Kenyan [7]	Ae. aegypti	 	
2005	La Réunion [69]	Ae. albopictus	nsP1: T301I	
nsP2:Y642N	
nsP3:S358P, 460del	
E1:A226V	
E2:Q146R	
2006	India [67]	Ae. albopictus	nsP1:A101V, T128K, T376M, T376M, Q488R, L507R	
nsP2:S54N, H374Y, T675M, A793V	
nsP3:V175I, Y217H, V331A, T337I, T341M, I376T, L461P, P471S	
nsP4:T75A, T254A, Q500L	
E1:K211N, V213I, M269V, D284E	
E2:A164T, I211T, T312M, S375T, V386A	
E3:I23T	
6K:V8I	
C:P23S, V27I, K63R, N80D	
2007, 2017	Italy [70,71]	Ae. albopictus	E1:A226V	
2007	Gabon [72]	Ae. albopictus	E1:A226V	
2008	China[73]	Ae. albopictus	nsP1:Q120R	
nsP2:G577R, N632S, L539S	
nsP3:M394I, D372N,	
nsP4:P181S	
E1:A306V, A226V	
E2:R178H, K252Q	
6K:V31I	
C:T8A	
2008, 2013, 2018	Thailand [74]	Ae. albopictus	E1:K211E	
E2:V264A	
2008	Malaysia [75]	Ae. albopictus	nsP1:S3P, S34P, T128K, V153I, M253K, T376M, G454S, R473S, A478T, N486D, Q488R, Q491R, H507R	
nsP2:L16P, S54N, S218T, L273Q, M338K, H374Y, V466M, V486I, L539S, I756V, S768N, A793V	
nsP3:T77S, G117R, V175I, I176V, V213M, N283S, V303T, V331A, R332Q, V334A, M336T, T337I, A349V, T353I, del376-382THTLPST, I383T, I413T, Q434L, A437V, I449M, R452Q, I457T, T458A, V459T, L461P, S462N, P471S, D483N, D484E	
nsP4:L42A, T58M, T75A, K85R, A90S, V101T, R235Q, T254A, R271K, D280E, A366T, Q500L, A582V	
C:P23S, V27I, K37Q, V48A, K73R, R78Q, M81T, V93A	
E1:S72N, T98A, A145T, E211K, S225A, A226V, M269V, D284E, S304P, P397L	
E2:I2T, H5N, G118S, R149K, A157V, A164T, S194G, D205G, S207N, S248L, K252Q, V/L255I, T312M, I317V, R318V, S375T, V384M, V386A	
E3:K33E, S44R, R60H, R62Q	
6K:V8I, M45T, T47A, L52M	
2009	India [44]	Ae. albopictus	E2:L210Q, K252Q	
2010, 2014	France [76]	Ae. albopictus	 	
2011	Cambodia[77]	Ae. albopictus	 	
2011	Republic of Congo [78]	Ae. aegypti,	E1:A226V	
Ae. albopictus	
2011	Indonesia [79]	Ae. albopictus	E1:A226V	
2012	Bhutan [80]	Ae. aegypti,	 	
Ae. albopictus	 	
2016	Pakistan[80]	Ae. albopictus	nsP1:K128T	
E1: T98A, S111T, A145T, K157N	
2016	Brazil	Ae. aegypti	E1: V156A, K211T	
2017	China [81]	Ae. albopictus	nsP3:D372E	
E1:M269V, D284E	
WA	1964	Nigeria [82]	Ae. aegypti	 	
1966, 1983	Senegal [82]	Ae. aegypti	 	
Asian	1958, 1988, 1995, 1990	Thailand [66]	Ae. aegypti	 	
1963, 1973	India[67]	Ae. aegypti	 	
1983,2008, 2011	Indonesia [79]	Ae. albopictus,Ae. aegypti	 	
1998, 2006	Malaysia[75]	Ae. albopictus,Ae. aegypti	 	
2011, 2013	New Caledonia [83]	Ae. aegypti	 	
2014	Tonga, American Samoa, the Independent States of Samoa, Tokelau [83]	Ae. aegypti	 	
2014	French Polynesia [84]	Ae. aegypti,	 	
Ae. polynesiensis	 	

To address the threat to human health posed by CHIKF, attention should be paid to host-associated adaptive changes in the virus. Although the mechanisms of adaptive evolution of viruses have not been fully understood, many studies have demonstrated the close relationship between viral mutations and changes in host adaptation, and new genotypes of mutant strains may cause large-scale outbreaks of CHIKF. For example, the E1-A226V and E2-L210Q mutations that appeared during natural evolution enhance the infectivity of Ae. albopictus [41,44]. The amino acid sites in which structural proteins undergo synergistic mutations are shown in the Figure 3i. Of these, E1-A226V is widely believed to have led to the expansion of CHIKV transmission. Subsequently, several secondary mutations have been identified (A226V:E2-L210Q, A226V:E2-K252Q, nsP2:L539S/E2-K252Q, and E2-R198Q/E3-S18F) that further enhance viral adaptation in Ae. albopictus [33], whereas E1-V156A, E1-K211T, E1-K211E and E2-V264A enhance adaptation in Ae. aegypti [45,56], and nsP4:R82S increases the efficiency of human-mosquito transmission. Various adaptive mutations have also been identified during laboratory passage. E2-G60D, E2-I211T, E1-V80I:A129V, and E2-K233Q/E, which appear in the background of the E1-A226V mutation and further enhance viral adaptation in Ae. albopictus [53,58]; E2-K200R:Δ3’UTR enhances viral adaptation, virulence, and transmission in mammals [59]. The nsP1:G230R/V326M/nsP3:*524 R triple mutant has enhanced viral replication and adaptation in mammals and Ae. albopictusin vitro and in vivo [60](Table 3). The fact that secondary mutations can further enhance viral adaptation suggests that adaptive changes to CHIKV may also be found in Ae. albopictus. However, it is worth noting that many ECSA isolates from recent outbreaks do not carry the E1 A226V mutation. In addition to the 3’UTR repeats present in the Caribbean strains mentioned above, Souza et al. sequenced CHIKV endemic to Tocantins, Brazil, between 2021 and 2022 and found that all 27 sequences belonged to the ECSA lineage, with 11 identified highly frequent missense mutations, while none of them had the E1-A226V substitution described earlier [85].Table 3. Amino acid mutations enhance host adaptation by CHIKV.

Gene	Mutation	Genotype	Effect	
E1	A226V	IOL	Increase adaptation and infectivity to Ae. albopictus [41]	
K211E	IOL	Increase adaptation to Ae. aegypti [45]	
V80I/A129V	IOL	Increase fitness for Ae. albopictus on a E1-A226V background [53]	
V80Q/M88L	IOL	Increase adaptation and transmission to Ae. aegypti and mammals on a E1-A226V background [55]	
V80Q/N20Y	IOL	Increase adaptation and transmission to Ae. aegypti on a E1-A226V background [19]	
V156A	ESCA (nature)	Increase pathogenicity of the virus by Aedes aegypti on a E1-A226V background [56]	
K211T	IOL (experiment)	
E2	G60D	ESCA	Increase fitness in Ae. albopictus and Aedes aegypti on a E1-A226V background [58]	
R198Q	IOL	Enhanced infection of Ae. albopictus, synergistic with E3:18F on a E1-A226V background [33]	
K200R	Asian	Enhanced viral transmissibility and pathogenicity in mammals [59]	
L210Q	IOL	Enhanced infection of Ae. albopictus on a E1-A226V background [44]	
I211T	IOL	Increase fitness in Ae. albopictus on a E1-A226V background [58]	
K233E/Q	IOL	Enhanced infection of Ae. albopictus on a E1-A226V background [33]	
K252Q	IOL	Enhanced infection of Ae. albopictus, synergistic with nsP2: L539S on a E1-A226V background [33]	
V264A	IOL	Enhanced fitness in Ae. aegypti on a E1-226A background [45]	
E3	S18F	IOL	Enhanced infection of Ae. albopictus, synergistic with E2:R198Q on a E1-A226V background [33]	
nsP1	G230R	IOL	Increase replication in Ae. albopictus in combination with nsP3:524* on a E1-A226V background [60]	
V326M	
nsP2	L539S	IOL	Enhanced infection of Ae. albopictus, synergistic with E2-K252Q on a E1-A226V background [33]	
nsP3	*524 R	IOL	Enhanced infection of Ae. albopictus, synergistic with nsP1: G230R/V326M on a E1-A226V background [60]	
nsP4	R82S	ESCA	Increase human-mosquito transmission [48]	
3’UTR	Deletion
(11921 –11,964 nt)	Asian	Increase viral pathogenicity for mammals compared to the E2-K200R single mutation [59]	

With the frequent movement of mosquitoes between endemic and non-endemic areas worldwide, it is likely that variants with greater infectivity and transmissivity will emerge in the future. Therefore, improved monitoring of genomic variations in CHIKV and focusing on mutations occurring during its evolution are important for controlling the global spread of CHIKV. Equally valuable is the introduction of adaptive mutations of CHIKV in the laboratory. Identifying new viral subtypes in nature can take years, while the function of newly discovered mutation sites can be validated in the laboratory through cell infection and in vivo mosquito-mouse transmission experiments. This approach would facilitate the targeted surveillance of CHIKV evolution by health agencies such as the CDC, enabling good responses to CHIKF outbreaks that can mitigate the impact of CHIKV mutations on public health.

Disclosure statement

No potential conflict of interest was reported by the author(s).

Author contributions

Ning XH and Xia BH wrote the draft of the paper, Wang JQ and Gao R collected and analysed the data, Ning XH and Ren H reviewed the paper critically, Ren H designed the paper and made the final approval of the version to be published. All authors agree to be accountable for all aspects of the work.

Data availability statement

No primary data is included in this article.
==== Refs
References

[1] Handler MZ, Handler NS, Stephany MP, et al. Chikungunya fever: an emerging viral infection threatening North America and Europe. Int J Dermatol. 2017 Feb;56 (2 ):e19–15. doi: 10.1111/ijd.13439 PubMed PMID: 27667734.27667734
[2] Cai L, Hu X, Liu S, et al. The research progress of Chikungunya fever. Frontiers in Public Health. 2022;10 :1095549. doi: 10.3389/fpubh.2022.1095549 PubMed PMID: 36699921; PubMed Central PMCID: PMCPMC9870324.36699921
[3] Bouquillard E, Fianu A, Bangil M, et al. Rheumatic manifestations associated with Chikungunya virus infection: a study of 307 patients with 32-month follow-up (RHUMATOCHIK study). Joint Bone Spine. 2018 Mar;85 (2 ):207–210. doi: 10.1016/j.jbspin.2017.01.014 PubMed PMID: 28238882.28238882
[4] Deqiao Tian WC. Chikungunya virus and chikungunya fever [J]. J Microbes And Infections. 2016;11 (4 ):194–206.
[5] Burt FJ, Chen W, Miner JJ, et al. Chikungunya virus: an update on the biology and pathogenesis of this emerging pathogen. The Lancet Infect Dis. 2017 Apr;17 (4 ):e107–e117. doi: 10.1016/s1473-3099(16)30385-1 PubMed PMID: 28159534.28159534
[6] Mavalankar D, Shastri P, Raman P. Chikungunya epidemic in India: a major public-health disaster. The Lancet Infectious Diseases. 2007 May;7 (5 ):306–307. doi: 10.1016/s1473-3099(07)70091-9 PubMed PMID: 17448932.17448932
[7] Sergon K, Njuguna C, Kalani R, et al. Seroprevalence of Chikungunya virus (CHIKV) infection on Lamu Island, Kenya, October 2004. Am J Trop Med Hyg. 2008 Feb;78 (2 ):333–337. PubMed PMID: 18256441. doi: 10.4269/ajtmh.2008.78.333 18256441
[8] Renault P, Solet JL, Sissoko D, et al. A major epidemic of chikungunya virus infection on Reunion Island, France, 2005-2006. Am J Trop Med Hyg. 2007 Oct;77 (4 ):727–731. PubMed PMID: 17978079. doi: 10.4269/ajtmh.2007.77.727 17978079
[9] Bessaud M, Peyrefitte CN, Pastorino BA, et al. Chikungunya virus strains, reunion island outbreak. Emerging infectious diseases. Emerg Infect Dis. 2006 Oct;12 (10 ):1604–1606. doi: 10.3201/eid1210.060596 PubMed PMID: 17176585; PubMed Central PMCID: PMCPMC3290959.17176585
[10] Yactayo S, Staples JE, Millot V, et al. Epidemiology of Chikungunya in the Americas. J Infect Dis. [2016 Dec 15];214 (suppl 5 ):S441–s445. doi: 10.1093/infdis/jiw390 PubMed PMID: 27920170; PubMed Central PMCID: PMCPMC5137246.27920170
[11] Jinhua Lu HZ, Zhang Z, Tian Z. Chikungunya fever and chikungunya virus. Chin J Frontier Health And Quarantine. 44 (5 ):375–378. doi: 10.16408/j.1004-9770.2021.05.023
[12] Yu X, Cheng G. Adaptive evolution as a driving force of the emergence and Re-emergence of mosquito-borne viral diseases. Viruses. [2022 Feb 21];14 (2 ). doi: 10.3390/v14020435 PubMed PMID: 35216028; PubMed Central PMCID: PMCPMC8878277.
[13] Liu Y, Liu J, Du S, et al. Evolutionary enhancement of Zika virus infectivity in aedes aegypti mosquitoes. Nature. [2017 May 25];545 (7655 ):482–486. doi: 10.1038/nature22365 PubMed PMID: 28514450; PubMed Central PMCID: PMCPMC5885636.28514450
[14] Yuan L, Huang XY, Liu ZY, et al. A single mutation in the prM protein of Zika virus contributes to fetal microcephaly. Sci (New Y, NY). [2017 Nov 17];358 (6365 ):933–936. doi: 10.1126/science.aam7120 PubMed PMID: 28971967.
[15] Wang E, Ni H, Xu R, et al. Evolutionary relationships of endemic/epidemic and sylvatic dengue viruses. J Virol. 2000 Apr;74 (7 ):3227–3234. doi: 10.1128/jvi.74.7.3227-3234.2000 PubMed PMID: 10708439; PubMed Central PMCID: PMCPMC111823.10708439
[16] Vignuzzi M, Stone JK, Arnold JJ, et al. Quasispecies diversity determines pathogenesis through cooperative interactions in a viral population. Nature. [2006 Jan 19];439 (7074 ):344–348. doi: 10.1038/nature04388 PubMed PMID: 16327776; PubMed Central PMCID: PMCPMC1569948.16327776
[17] Wolinsky SM, Korber BT, Neumann AU, et al. Adaptive evolution of human immunodeficiency virus-type 1 during the natural course of infection. Sci (New Y, NY). [1996 Apr 26];272 (5261 ):537–542. doi: 10.1126/science.272.5261.537 PubMed PMID: 8614801.
[18] Farci P, Shimoda A, Coiana A, et al. The outcome of acute hepatitis C predicted by the evolution of the viral quasispecies. Sci (New Y, NY). [2000 Apr 14];288 (5464 ):339–344. doi: 10.1126/science.288.5464.339 PubMed PMID: 10764648.
[19] Godaert L, Najioullah F, Bartholet S, et al. Atypical clinical presentations of acute phase chikungunya virus infection in older adults. J Am Geriatr Soc. 2017 Nov;65 (11 ):2510–2515. doi: 10.1111/jgs.15004 PubMed PMID: 28940357.28940357
[20] An W, Ge N, Cao Y, et al. Recent progress on chikungunya virus research. Virol Sin. 2017 Dec;32 (6 ):441–453. doi: 10.1007/s12250-017-4072-x PubMed PMID: 29188427; PubMed Central PMCID: PMCPMC6702244.29188427
[21] Yin H, Yin P, Zhao H, et al. Intraviral interactome of Chikungunya virus reveals the homo-oligomerization and palmitoylation of structural protein TF. Biochem Biophys Res Commun. [2019 Jun 11];513 (4 ):919–924. doi: 10.1016/j.bbrc.2019.04.098 PubMed PMID: 31005258; PubMed Central PMCID: PMCPMC7092876.31005258
[22] Singh A, Jain R. Neurological manifestations of Chikungunya in children. Indian Pediatr. [2017 Mar 15];54 (3 ):249. doi: 10.1007/s13312-017-1040-3 PubMed PMID: 28361790.28361790
[23] Peränen J, Laakkonen P, Hyvönen M, et al. The alphavirus replicase protein nsP1 is membrane-associated and has affinity to endocytic organelles. Virology. [1995 Apr 20];208 (2 ):610–620. doi: 10.1006/viro.1995.1192 PubMed PMID: 7747433.7747433
[24] Pradeep P, Sivakumar KC, Sreekumar E, et al. Host factor nucleophosmin 1 (NPM1/B23) exerts antiviral effects against Chikungunya virus by its interaction with viral nonstructural protein 3. Microbiol Spectr. [2023 Aug 17];11 (4 ):e0537122. doi: 10.1128/spectrum.05371-22 PubMed PMID: 37409962; PubMed Central PMCID: PMCPMC10433958.37409962
[25] Li Li TX, Li X. Molecular variation and epidemiology of the Chikungunya virus. China Trop Med. 2021;21 (8 ):809–813. doi: 10.13604/j.cnki.46-1064/r.2021.08.19
[26] Schwartz O, Albert ML. Biology and pathogenesis of chikungunya virus. Nat Rev Microbiol. 2010 Jul;8 (7 ):491–500. doi: 10.1038/nrmicro2368 PubMed PMID: 20551973.20551973
[27] Jose J, Snyder JE, Kuhn RJ. A structural and functional perspective of alphavirus replication and assembly. Future Microbiol. 2009 Sep;4 (7 ):837–856. doi: 10.2217/fmb.09.59 PubMed PMID: 19722838; PubMed Central PMCID: PMCPMC2762864.19722838
[28] Gorchakov R, Wang E, Leal G, et al. Attenuation of Chikungunya virus vaccine strain 181/clone 25 is determined by two amino acid substitutions in the E2 envelope glycoprotein. J Virol. 2012 Jun;86 (11 ):6084–6096. doi: 10.1128/jvi.06449-11 PubMed PMID: 22457519; PubMed Central PMCID: PMCPMC3372191.22457519
[29] Gardner CL, Hritz J, Sun C, et al. Deliberate attenuation of chikungunya virus by adaptation to heparan sulfate-dependent infectivity: a model for rational arboviral vaccine design. PLOS neglected tropical diseases. 2014 Feb;8 (2 ):e2719. doi: 10.1371/journal.pntd.0002719 PubMed PMID: 24587470; PubMed Central PMCID: PMCPMC3930508.
[30] Domingo E. Mechanisms of viral emergence. Vet Res. 2010 Nov–Dec;41 (6 ):38. doi: 10.1051/vetres/2010010 PubMed PMID: 20167200; PubMed Central PMCID: PMCPMC2831534.20167200
[31] Hakim MS, Annisa L, Gazali FM, et al. The origin and continuing adaptive evolution of chikungunya virus. Arch Virol. 2022 Dec;167 (12 ):2443–2455. doi: 10.1007/s00705-022-05570-z PubMed PMID: 35987965.35987965
[32] Chen R, Puri V, Fedorova N, et al. Comprehensive genome scale phylogenetic study provides new insights on the global expansion of Chikungunya virus. J Virol. [2016 Dec 1];90 (23 ):10600–10611. doi: 10.1128/jvi.01166-16 PubMed PMID: 27654297; PubMed Central PMCID: PMCPMC5110187.27654297
[33] Tsetsarkin KA, Chen R, Yun R, et al. Multi-peaked adaptive landscape for chikungunya virus evolution predicts continued fitness optimization in aedes albopictus mosquitoes. Nat Commun. [2014 Jun 16];5 :4084. doi: 10.1038/ncomms5084 PubMed PMID: 24933611; PubMed Central PMCID: PMCPMC7091890.24933611
[34] Brown JE, Evans BR, Zheng W, et al. Human impacts have shaped historical and recent evolution in aedes aegypti, the dengue and yellow fever mosquito. Evolution. 2014 Feb;68 (2 ):514–525. doi: 10.1111/evo.12281 PubMed PMID: 24111703; PubMed Central PMCID: PMCPMC3946797.24111703
[35] Thiberville SD, Moyen N, Dupuis-Maguiraga L, et al. Chikungunya fever: epidemiology, clinical syndrome, pathogenesis and therapy. Antiviral Res. 2013 Sep;99 (3 ):345–370. doi: 10.1016/j.antiviral.2013.06.009 PubMed PMID: 23811281; PubMed Central PMCID: PMCPMC7114207.23811281
[36] Delatte H, Gimonneau G, Triboire A, et al. Influence of temperature on immature development, survival, longevity, fecundity, and gonotrophic cycles of aedes albopictus, vector of chikungunya and dengue in the Indian Ocean. J Med Entomol. 2009 Jan;46 (1 ):33–41. doi: 10.1603/033.046.0105 PubMed PMID: 19198515.19198515
[37] Kraemer MU, Sinka ME, Duda KA, et al. The global distribution of the arbovirus vectors aedes aegypti and ae. albopictus eLife. [2015 Jun 30];4 :e08347. doi: 10.7554/eLife.08347 PubMed PMID: 26126267; PubMed Central PMCID: PMCPMC4493616.26126267
[38] Medlock JM, Hansford KM, Schaffner F, et al. A review of the invasive mosquitoes in Europe: ecology, public health risks, and control options. Vector Borne And Zoonotic Dis (Larchmt, NY). 2012 Jun;12 (6 ):435–447. doi: 10.1089/vbz.2011.0814 PubMed PMID: 22448724; PubMed Central PMCID: PMCPMC3366101.
[39] Deeba F, Haider MSH, Ahmed A, et al. Global transmission and evolutionary dynamics of the Chikungunya virus. Epidemiol Infect. [2020 Feb 19];148 :e63. doi: 10.1017/s0950268820000497 PubMed PMID: 32070451; PubMed Central PMCID: PMCPMC7118414.32070451
[40] Bordería AV, Isakov O, Moratorio G, et al. Group selection and contribution of minority variants during virus adaptation determines virus fitness and phenotype. PloS Pathog. 2015 May;11 (5 ):e1004838. doi: 10.1371/journal.ppat.1004838 PubMed PMID: 25941809; PubMed Central PMCID: PMCPMC4420505.25941809
[41] Tsetsarkin KA, Vanlandingham DL, McGee CE, et al. A single mutation in chikungunya virus affects vector specificity and epidemic potential. PloS Pathog. 2007 Dec;3 (12 ):e201. doi: 10.1371/journal.ppat.0030201 PubMed PMID: 18069894; PubMed Central PMCID: PMCPMC2134949.18069894
[42] Wu D, Zhang Y, Zhouhui Q, et al. Chikungunya virus with E1-A226V mutation causing two outbreaks in 2010, Guangdong, China. Virol J. [2013 Jun 2];10 :174. doi: 10.1186/1743-422x-10-174 PubMed PMID: 23725047; PubMed Central PMCID: PMCPMC3691762.
[43] Fortuna C, Toma L, Remoli ME, et al. Vector competence of aedes albopictus for the Indian Ocean lineage (IOL) chikungunya viruses of the 2007 and 2017 outbreaks in Italy: a comparison between strains with and without the E1: A226V mutation. Euro surveillance: bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin. 2018 May;23 (22 ). doi: 10.2807/1560-7917.Es.2018.23.22.1800246 PubMed PMID: 29871722; PubMed Central PMCID: PMCPMC6152176.
[44] Tsetsarkin KA, Weaver SC, Heise M. Sequential adaptive mutations enhance efficient vector switching by Chikungunya virus and its epidemic emergence. PLOS Pathog. 2011 Dec;7 (12 ):e1002412. doi: 10.1371/journal.ppat.1002412 PubMed PMID: 22174678; PubMed Central PMCID: PMCPMC3234230.22174678
[45] Agarwal A, Sharma AK, Sukumaran D, et al. Two novel epistatic mutations (E1: K211E and E2: V264A) in structural proteins of Chikungunya virus enhance fitness in aedes aegypti. Virology. 2016 Oct;497 :59–68. doi: 10.1016/j.virol.2016.06.025 PubMed PMID: 27423270.27423270
[46] Maljkovic Berry I, Eyase F, Pollett S, et al. Global outbreaks and origins of a Chikungunya virus variant carrying mutations which May increase fitness for aedes aegypti: revelations from the 2016 Mandera, Kenya outbreak. Am J Trop Med Hyg. 2019 May;100 (5 ):1249–1257. doi: 10.4269/ajtmh.18-0980 PubMed PMID: 30860010; PubMed Central PMCID: PMCPMC6493958.30860010
[47] Hapuarachchi HC, Bandara KB, Sumanadasa SD, et al. Re-emergence of Chikungunya virus in South-east Asia: virological evidence from Sri Lanka and Singapore. J Gen Virol. 2010 Apr;91 (Pt 4 ):1067–1076. doi: 10.1099/vir.0.015743-0 PubMed PMID: 19955565.19955565
[48] Fu JYL, Chua CL, Vythilingam I, et al. An amino acid change in nsP4 of chikungunya virus confers fitness advantage in human cell lines rather than in aedes albopictus. J Gen Virol. 2019 Nov;100 (11 ):1541–1553. doi: 10.1099/jgv.0.001338 PubMed PMID: 31613205.31613205
[49] Stapleford KA, Moratorio G, Henningsson R, et al. Whole-genome sequencing analysis from the Chikungunya virus caribbean outbreak reveals novel evolutionary genomic elements. PLOS neglected tropical diseases. 2016 Jan;10 (1 ):e0004402. doi: 10.1371/journal.pntd.0004402 PubMed PMID: 26807575; PubMed Central PMCID: PMCPMC4726740.
[50] Merwaiss F, Filomatori CV, Susuki Y, et al. Chikungunya virus replication rate determines the capacity of crossing tissue barriers in mosquitoes. J Virol. [2021 Jan 13];95 (3 ). doi: 10.1128/jvi.01956-20 PubMed PMID: 33148794; PubMed Central PMCID: PMCPMC7925089.
[51] Chen R, Wang E, Tsetsarkin KA, et al. Chikungunya virus 3’ untranslated region: adaptation to mosquitoes and a population bottleneck as major evolutionary forces. PLOS Pathog. 2013;9 (8 ):e1003591. doi: 10.1371/journal.ppat.1003591 PubMed PMID: 24009512; PubMed Central PMCID: PMCPMC3757053.24009512
[52] Filomatori CV, Bardossy ES, Merwaiss F, et al. RNA recombination at Chikungunya virus 3’UTR as an evolutionary mechanism that provides adaptability. PLOS Pathog. 2019 Apr;15 (4 ):e1007706. doi: 10.1371/journal.ppat.1007706 PubMed PMID: 30986247; PubMed Central PMCID: PMCPMC6502353.30986247
[53] Stapleford KA, Coffey LL, Lay S, et al. Emergence and transmission of arbovirus evolutionary intermediates with epidemic potential. Cell host & microbe. [2014 Jun 11];15 (6 ):706–716. doi: 10.1016/j.chom.2014.05.008 PubMed PMID: 24922573.
[54] Noval MG, Rodriguez-Rodriguez BA, Rangel MV, et al. Evolution-driven attenuation of alphaviruses highlights key glycoprotein determinants regulating viral infectivity and dissemination. Cell Rep. [2019 Jul 9];28 (2 ):460–471.e5. doi: 10.1016/j.celrep.2019.06.022 PubMed PMID: 31291581; PubMed Central PMCID: PMCPMC7141928.31291581
[55] Thannickal SA, Battini L, Spector SN, et al. The chikungunya virus E1 glycoprotein fusion loop and hinge alter glycoprotein dynamics leading to cell and host specific changes in infectivity. bioRxiv: the preprint server for biology. [2023 Nov 5]. doi: 10.1101/2023.11.03.565585 bioRxiv: the preprint server for biology.
[56] Rangel MV, McAllister N, Dancel-Manning K, et al. Emerging chikungunya virus variants at the E1-E1 interglycoprotein spike interface impact virus attachment and inflammation. J Virol. [2022 Feb 23];96 (4 ):e0158621. doi: 10.1128/jvi.01586-21 PubMed PMID: 34935436; PubMed Central PMCID: PMCPMC8865554.34935436
[57] Souza TM, Azeredo EL, Badolato-Corrêa J, et al. First report of the east-central South African genotype of Chikungunya virus in Rio de Janeiro, Brazil. PloS Curr. [2017 Feb 14];9 . doi: 10.1371/currents.outbreaks.4200119978d62ccaa454599cd2735727 PubMed PMID: 28286701; PubMed Central PMCID: PMCPMC5325710.
[58] Tsetsarkin KA, McGee CE, Volk SM, et al. Epistatic roles of E2 glycoprotein mutations in adaption of chikungunya virus to aedes albopictus and ae. aegypti mosquitoes. PLOS ONE. [2009 Aug 31];4 (8 ):e6835. doi: 10.1371/journal.pone.0006835 PubMed PMID: 19718263; PubMed Central PMCID: PMCPMC2729410.19718263
[59] Hawman DW, Carpentier KS, Fox JM, et al. Mutations in the E2 glycoprotein and the 3’ untranslated region enhance chikungunya virus virulence in mice. J Virol. [2017 Oct 15];91 (20 ). doi: 10.1128/jvi.00816-17 PubMed PMID: 28747508; PubMed Central PMCID: PMCPMC5625495.
[60] Mounce BC, Cesaro T, Vlajnić L, et al. Chikungunya virus overcomes polyamine depletion by mutation of nsP1 and the opal stop codon to confer enhanced replication and fitness. J Virol. [2017 Aug 1];91 (15 ). doi: 10.1128/jvi.00344-17 PubMed PMID: 28539441; PubMed Central PMCID: PMCPMC5512238.
[61] Abraham R, Manakkadan A, Mudaliar P, et al. Correlation of phylogenetic clade diversification and in vitro infectivity differences among cosmopolitan genotype strains of Chikungunya virus. Infect Genet Evol. 2016 Jan;37 :174–184. doi: 10.1016/j.meegid.2015.11.019 PubMed PMID: 26611825.26611825
[62] Tandel K, Kumar M, Shergill SPS, et al. Molecular characterization and phylogenetic analysis of Chikungunya virus from Delhi, India. Med J Armed Forces India. 2019 Jul;75 (3 ):266–273. doi: 10.1016/j.mjafi.2018.02.004 PubMed PMID: 31388228; PubMed Central PMCID: PMCPMC6676464.31388228
[63] Harsha PK, Reddy V, Rao D, et al. Continual circulation of ECSA genotype and identification of a novel mutation I317V in the E1 gene of Chikungunya viral strains in southern India during 2015-2016. J Med Virol. 2020 Aug;92 (8 ):1007–1012. doi: 10.1002/jmv.25662 PubMed PMID: 31900943.31900943
[64] Rodríguez-Aguilar ED, Martínez-Barnetche J, González-Bonilla CR, et al. Genetic diversity and spatiotemporal dynamics of chikungunya infections in Mexico during the outbreak of 2014–2016. Viruses. [2021 Dec 31];14 (1 ):70. doi: 10.3390/v14010070 PubMed PMID: 35062275; PubMed Central PMCID: PMCPMC8779743.
[65] Tsetsarkin KA, Chen R, Sherman MB, et al. Chikungunya virus: evolution and genetic determinants of emergence. Curr Opin Virol. 2011 Oct;1 (4 ):310–317. doi: 10.1016/j.coviro.2011.07.004 PubMed PMID: 21966353; PubMed Central PMCID: PMCPMC3182774.21966353
[66] Zeller H, Van Bortel W, Sudre B. Chikungunya: its history in Africa and Asia and its spread to new regions in 2013–2014. J Infect Dis. [2016 Dec 15];214 (suppl 5 ):S436–s440. doi: 10.1093/infdis/jiw391 PubMed PMID: 27920169.
[67] Arankalle VA, Shrivastava S, Cherian S, et al. Genetic divergence of Chikungunya viruses in India (1963–2006) with special reference to the 2005–2006 explosive epidemic. J Gen Virol. 2007 Jul;88 (Pt 7 ):1967–1976. doi: 10.1099/vir.0.82714-0 PubMed PMID: 17554030.17554030
[68] Collao X, Negredo AI, Cano J, et al. Different lineages of Chikungunya virus in Equatorial Guinea in 2002 and 2006. Am J Trop Med Hyg. 2010 Mar;82 (3 ):505–507. doi: 10.4269/ajtmh.2010.09-0435 PubMed PMID: 20207882; PubMed Central PMCID: PMCPMC2829918.20207882
[69] Schuffenecker I, Iteman I, Michault A, et al. Genome microevolution of chikungunya viruses causing the Indian Ocean outbreak. PLOS Med. 2006 Jul;3 (7 ):e263. doi: 10.1371/journal.pmed.0030263 PubMed PMID: 16700631; PubMed Central PMCID: PMCPMC1463904.16700631
[70] Rezza G, Nicoletti L, Angelini R, et al. Infection with chikungunya virus in Italy: an outbreak in a temperate region. Lancet (London, England). 2007 Dec 1;370 (9602 ):1840–1846. doi: 10.1016/s0140-6736(07)61779-6 PubMed PMID: 18061059.18061059
[71] Venturi G, Fabiani M, Amendola A, et al. Lack of evidence of Chikungunya virus infection among blood donors during the Chikungunya outbreak in Lazio Region, Italy, 2017. Viruses. [2022 Mar 16];14 (3 ). doi: 10.3390/v14030619 PubMed PMID: 35337026; PubMed Central PMCID: PMCPMC8950832.
[72] Leroy EM, Nkoghe D, Ollomo B, et al. Concurrent chikungunya and dengue virus infections during simultaneous outbreaks, Gabon, 2007. Emerg Infect Dis. 2009 Apr;15 (4 ):591–593. doi: 10.3201/eid1504.080664 PubMed PMID: 19331740; PubMed Central PMCID: PMCPMC2671412.19331740
[73] Zheng K, Li J, Zhang Q, et al. Genetic analysis of chikungunya viruses imported to mainland China in 2008. Virol J. [2010 Jan 18];7 :8. doi: 10.1186/1743-422x-7-8 PubMed PMID: 20078896; PubMed Central PMCID: PMCPMC2831882.
[74] Khongwichit S, Chansaenroj J, Thongmee T, et al. Large-scale outbreak of Chikungunya virus infection in Thailand, 2018–2019. PLOS ONE. 2021;16 (3 ):e0247314. doi: 10.1371/journal.pone.0247314 PubMed PMID: 33690657; PubMed Central PMCID: PMCPMC7946318.33690657
[75] Sam IC, Loong SK, Michael JC, et al. Genotypic and phenotypic characterization of Chikungunya virus of different genotypes from Malaysia. PLOS ONE. 2012;7 (11 ):e50476. doi: 10.1371/journal.pone.0050476 PubMed PMID: 23209750; PubMed Central PMCID: PMCPMC3507689.23209750
[76] Delisle E, Rousseau C, Broche B, et al. Chikungunya outbreak in Montpellier, France, September to October 2014. Euro Surveillance: Bull Europeen sur les Mal transmissibles = Eur Communicable Disease Bull. [2015 Apr 30];20 (17 ). doi: 10.2807/1560-7917.es2015.20.17.21108 PubMed PMID: 25955774.
[77] Sasayama M, Benjathummarak S, Kawashita N, et al. Chikungunya virus was isolated in Thailand, 2010. Virus Genes. 2014 Dec;49 (3 ):485–489. doi: 10.1007/s11262-014-1105-5 PubMed PMID: 25113745; PubMed Central PMCID: PMCPMC4232745.25113745
[78] Mombouli JV, Bitsindou P, Elion DO, et al. Chikungunya virus infection, Brazzaville, Republic of Congo, 2011. Emerg Infect Dis. 2013;19 (9 ):1542–1543. doi: 10.3201/eid1909.130451 PubMed PMID: 23968609; PubMed Central PMCID: PMCPMC3810930.23968609
[79] Maha MS, Susilarini NK, Hariastuti NI. Chikungunya virus mutation, Indonesia, 2011. Emerg Infect Dis. 2015 Feb;21 (2 ):379–381. doi: 10.3201/eid2102.141121 PubMed PMID: 25625600; PubMed Central PMCID: PMCPMC4313651.25625600
[80] Wangchuk S, Chinnawirotpisan P, Dorji T, et al. Chikungunya fever outbreak, Bhutan, 2012. Emerg Infect Dis. 2013 Oct;19 (10 ):1681–1684. doi: 10.3201/eid1910.130453 PubMed PMID: 24047543; PubMed Central PMCID: PMCPMC3810753.24047543
[81] Pan J, Fang C, Yan J, et al. Chikungunya fever outbreak, Zhejiang Province, China, 2017. Emerg Infect Dis. 2019 Aug;25 (8 ):1589–1591. doi: 10.3201/eid2508.181212 PubMed PMID: 31310205; PubMed Central PMCID: PMCPMC6649353.31310205
[82] Powers AM, Brault AC, Tesh RB, et al. Re-emergence of Chikungunya and O’nyong-nyong viruses: evidence for distinct geographical lineages and distant evolutionary relationships. J Gen Virol. 2000 Feb;81 (Pt 2 ):471–479. doi: 10.1099/0022-1317-81-2-471 PubMed PMID: 10644846.10644846
[83] Cao-Lormeau VM, Musso D. Emerging arboviruses in the Pacific. Lancet (London, England). [2014 Nov 1];384 (9954 ):1571–1572. doi: 10.1016/s0140-6736(14)61977-2 PubMed PMID: 25443481.25443481
[84] Aubry M, Teissier A, Roche C, et al. Chikungunya outbreak, French polynesia, 2014. Emerging Infectious Diseases. 2015 Apr;21 (4 ):724–726. doi: 10.3201/eid2104.141741 PubMed PMID: 25811534; PubMed Central PMCID: PMCPMC4378499.25811534
[85] Souza UJB, Santos RND, Giovanetti M, et al. Genomic epidemiology reveals the circulation of the chikungunya virus East/Central/South African lineage in Tocantins State, North Brazil. Viruses. [2022 Oct 21];14 (10 ). doi: 10.3390/v14102311 PubMed PMID: 36298867; PubMed Central PMCID: PMCPMC9611869.
