key: cord-0023034-zydm23pf
authors: Chen, Rubing; Plante, Jessica A.; Plante, Kenneth S.; Yun, Ruimei; Shinde, Divya; Liu, Jianying; Haller, Sherry; Mukhopadhyay, Suchetana; Weaver, Scott C.
title: Lineage Divergence and Vector-Specific Adaptation Have Driven Chikungunya Virus onto Multiple Adaptive Landscapes
date: 2021-11-09
journal: nan
DOI: 10.1128/mbio.02738-21
sha: 5b16027d430b7783876394f6c4e992017a044ce7
doc_id: 23034
cord_uid: zydm23pf

Previous studies have shown that the adaptation of Indian Ocean lineage (IOL) chikungunya virus (CHIKV) strains for Aedes albopictus transmission was mediated by an E1-A226V substitution, followed by either a single substitution in E2 or synergistic substitutions in the E2 and E3 envelope glycoproteins. Here, we examined whether Asian lineage strains, including those that descended from the 2014 Caribbean introduction, are likely to acquire these A. albopictus-adaptive E2 substitutions. Because Asian lineage strains cannot adapt through the E1-A226V substitution due to an epistatic constraint, we first determined that the beneficial effect of these E2 mutations in IOL strains is independent of E1-A226V. We then introduced each of these E2 adaptive mutations into the Asian lineage backbone to determine if they improve infectivity for A. albopictus. Surprisingly, our results indicated that in the Asian lineage backbone, these E2 mutations significantly decreased CHIKV fitness in A. albopictus. Furthermore, we tested the effects of these mutations in Aedes aegypti and observed different results from those in A. albopictus, suggesting that mosquito species-specific factors that interact with the envelope proteins are involved in vector infection efficiency. Overall, our results indicate that the divergence between Asian lineage and IOL CHIKVs has led them onto different adaptive landscapes with differing potentials to expand their vector host range.

. Aedes aegypti is broadly distributed in tropical and subtropical regions throughout the world, whereas A. albopictus can also survive colder climates and thus is found in temperate regions of China, the Americas, Europe, and eastern Australia (3) . Additionally, whereas A. aegypti is highly endophilic and prefers taking blood meals from humans, A. albopictus is more exophilic and in some locations bites a variety of domestic and wild animals in addition to humans. Vector competence is therefore an important factor dictating the distribution of mosquito-borne viruses, and understanding the mechanism(s) of adaptation to a given vector may allow us to anticipate patterns of spread by mosquito-borne viruses and redirect public health efforts accordingly.

CHIKV is an alphavirus in the family Togaviridae and is arguably the most medically important of all alphaviruses due to its widespread distribution, history of explosive outbreaks, and frequency of chronic disease symptoms. CHIKV infection leads to chikungunya fever, accompanied by severe, debilitating, and often chronic arthralgia that has major economic as well as public health impacts (4) . Rooted in sub-Saharan African sylvatic transmission cycles, CHIKV strains have emerged at least twice during the past century from the East, Central, South Africa (ECSA) enzootic lineage into Asia, forming independent urban lineages (5) . The first lineage identified during the modern scientific era, termed the Asian lineage, was introduced from Africa between 1879 and 1956 and subsequently generated epidemics in Southeast Asia and India. Although not detected in India since 1973, the Asian lineage is still causing outbreaks in Southeast Asia and Oceania. In addition, it was the Asian lineage of CHIKV that first appeared in the Americas via the Caribbean in late 2013 (6) , followed by its rapid spread throughout much of Latin America (7-9) as well as Florida (10) and Texas (11) in the United States. The only mosquito species directly incriminated as an urban vector during these American epidemics involving approximately 2.4 million persons has been A. aegypti (12) (13) (14) , although A. albopictus is also present throughout most of the affected region.

The second modern urban CHIKV emergence occurred in 2004 when an ECSA lineage strain initiated an outbreak in coastal Kenya and then independently spread to the Indian Ocean islands and India in 2005 and 2006, respectively (5) , forming the epidemic Indian Ocean lineage (IOL). Initial IOL circulation in East Africa relied principally on A. aegypti (15) . However, soon after this CHIKV strain reached La Réunion Island, A. albopictus began to play a major role in transmission. This vector switch was associated with an A226V substitution in the E1 envelope glycoprotein (16) that was later shown to enhance infection and dissemination in A. albopictus (17, 18) . Similarly, whereas the initial outbreaks in India were caused by strains with E1-226A (19) , E1-226V emerged beginning in 2007 in the state of Kerala, where A. albopictus is abundant, and the subsequent distribution of E1-226A versus E1-226V has reflected the abundance of these two vectors (20) (21) (22) . From India, IOL strains rapidly spread not only through the entire Indian Subcontinent and Southeast Asia but also to Europe and Oceania, where A. aegypti is not abundant (5, 23) .

The continued, rapid adaptation of the CHIKV IOL to A. albopictus has been attributed to a series of point mutations on all three of the envelope proteins, E1, E2, and E3 (17, 24) . Although infection of midgut epithelial cells following an infectious blood meal initiates vector infection (25) , it is not clear which stage(s) of the viral replication cycle is directly targeted or altered to determine the expanded vector range mediated by these envelope protein substitutions. The initial adaptation of IOL CHIKV strains to A. albopictus occurred via the convergent E1-A226V substitution. This residue is in the E1ij loop that is crucial for the pH-sensitive E2/E1 conformational change that exposes the nearby fusion peptide during endosomal entry (26) . Although the E1-A226V substitution improves the infectivity of IOL strains for A. albopictus by about 40-fold, it has little or no effect on fitness in A. aegypti, suggesting that a mosquito species-specific factor(s) determines the importance of this substitution (17) .

Following the spread of different CHIKV IOL sublineages, some of which acquired E1-226V, secondary adaptive mutations in E2 (L210Q and K252Q) or an E2-R198Q/E3-S18F synergistic pair further enhanced infection of A. albopictus (24) . In addition, an E2-K233E substitution that was not naturally associated with E1-226V also enhances infection of A. albopictus in the presence of E1-226V (24) . Based on the atomic resolution structure of the CHIKV spike, these E2 mutations are located in or near the acid-sensitive region (ASR) that is critical for triggering E1/E2 disassociation in response to low pH, which initiates E1 trimerization and endosomal membrane fusion (24) . Moreover, combining E2-210Q and E2-252Q, or introducing the substitution E2-233E, further enhances infection of A. albopictus, indicating that IOL strains can readily explore their complex fitness landscape for additional adaptive pathways to improve infectivity for this widespread mosquito vector. However, as with E1-226V, these E2 mutations do not significantly influence the infectivity of IOL strains in A. aegypti (24) , indicating that their beneficial effects in A. albopictus are mediated by species-specific factors.

In contrast to the IOL, there is no evidence that the Asian CHIKV lineage has adapted to A. albopictus as a primary vector despite circulating in its native territory for at least 60 years. Furthermore, the A. albopictus-adaptive E1, E2, and E3 mutations observed in IOL strains have never been detected in an Asian lineage strain. A mutation unique to the Asian lineage, E1-A98T, presumably acquired via a founder effect during its introduction from Africa, imposes an epistatic constraint on the E1-226 residue, thus precluding at least one A. albopictus-adaptive pathway of the Asian lineage (27, 48) . However, it is unclear why the Asian lineage has not adapted to A. albopictus via other mutations such as those in the E2 gene observed in some IOL strains. One possibility is that these E2 mutations are dependent on E1-226V since most of them (except E2-K233E) appeared only after the E1-A226V mutation was fixed in certain IOL sublineages. In addition, an epistatic interaction also occurs between E1-226V and E2-211T, an amino acid located in the same E2 ASR (28) , suggesting that epistasis could affect the entire spectrum of E2 vector-adaptive mutations. Alternatively, if the appearance of the A. albopictus-adaptive E2 mutations was independent of E1-226V, their later appearance (after E1-A226V) may reflect only their smaller fitness benefit and, hence, lower selection coefficient (;10-fold versus 40-fold [17, 24] ). In that case, these IOL E2 adaptive mutations, namely, E2-198Q, E2-210Q, E2-233E, and E2-252Q, might eventually also be selected by A. albopictus in Asian lineage CHIKV strains to enhance transmission efficiency, especially in temperate climates and rural locations.

To test these hypotheses, we engineered CHIKV strains from the IOL and Asian lineage with these specific E2 mutations and evaluated their effects on infection and dissemination in A. albopictus and A. aegypti mosquitoes. Our main goals were to (i) determine if the A. albopictus-adaptive E2 mutations in IOL strains are dependent on E1-226V, (ii) evaluate the potential of these E2 mutations to emerge in Asian lineage CHIKV as an adaptation to A. albopictus, and (iii) determine if these E2 mutations affect infection of A. aegypti.

Beneficial effects of IOL E2 A. albopictus-adaptive mutations do not depend on E1-226V. Four A. albopictus-adaptive E2 substitutions (E2-R198Q, E2-L210Q, E2-K233E, and E2-K252Q) previously shown to increase fitness in an IOL (SL07) backbone in the context of E1-226V were introduced into the SL07 IOL backbone encoding E1-226A, referred to here as the SL07 wild type (wt). Approximately equal amounts of the wt and E2 mutant viruses as measured by genome ratios, one with an introduced synonymous ApaI/PspOMI restriction site marker and one without, were mixed, and A. albopictus mosquitoes were infected via an artificial blood meal. To account for potential bias caused by the marker, reciprocal experiments were conducted where the marker was swapped into the competing backbone. Mosquito heads were collected on day 10, and the differences in wt/mutant ratios in the heads compared to the blood meals were used to assess the mutation's impact on fitness for infection and dissemination.

The majority of mosquito heads had a relatively pure strain composition, with 248 of 265 containing $95% of a single competitor strain. This may reflect, in part, the population bottleneck that occurs upon initial midgut infection and the resultant introduction of a stochastic component to the competition outcome. However, the E2 mutations consistently conferred a fitness advantage to SL07 E1-226A in A. albopictus (Fig. 1 ). This advantage was measured regardless of which strain contained the restriction marker, which did not significantly affect the calculated relative replicative fitness values (two-tailed paired t test, t = 0.4683; df = 3; P = 0.6715). The fitness advantage was highly significant for E2-198Q, E2-233E, and E2-252Q (P , 0.0001): E2-198Q resulted in an 84-or 22-fold advantage, E2-233E resulted in a 30-or 50-fold advantage, and E2-252Q resulted in a 17-or 20-fold advantage depending on whether the marker was located in the mutant or wt strains, respectively. The E2-210Q mutation, on the other hand, conferred a more modest 4-or 8-fold advantage. These results indicate that the beneficial effects of the E2 A. albopictus-adaptive mutations are not contingent upon E1-226V in an IOL backbone.

IOL E2 A. albopictus-adaptive mutations in the Asian lineage are neutral or deleterious in A. albopictus. The E2-198Q, E2-233E, and E2-252Q mutations were introduced into the Asian lineage strain Mal06, which naturally contains E1-226A. Because two point mutations are required to convert L to Q at E2-210 in the Asian lineage, this substitution is FIG 1 Impact of E2 mutations on the fitness of IOL CHIKV in Aedes albopictus in the absence of E1-226V. Aedes albopictus mosquitoes were fed artificial blood meals with approximately equal genomic ratios of wt and E2 mutant CHIKV SL07, one with an introduced ApaI marker and one without. At 10 dpi, heads were collected, and their virus mixture was passaged for 2 days in Vero cells. The resulting cell culture supernatant was amplified by RT-PCR, and the ratio of ApaI-marked to unmarked CHIKV was determined by restriction digestion. Reciprocal experiments were conducted such that the ApaI marker was present in either the mutant (A) or the wt (B) strain. Relative replicative fitness values were calculated as described in Materials and Methods. The Holm-Sidak correction for multiple P values was applied. *, P , 0.05; **, P , 0.01; ***, P , 0.001; ****, P , 0.0001.

® considered unlikely in nature and was therefore excluded from the investigation. As was the case in the SL07 backbone, the strain containing the introduced restriction marker did not significantly impact the resulting changes in fitness (two-tailed paired t test, t = 0.8188; df = 2; P = 0.4990), and the vast majority of heads (134 of 146) contained $95% of a single strain.

In contrast to the IOL lineage SL07 backbone, the Asian lineage Mal06 backbone was either unaffected or negatively impacted by the E2 mutations (Fig. 2) . E2-198Q resulted in no significant fitness change, with a 2-fold decrease in replicative fitness regardless of which strain contained the restriction marker. E2-252Q also resulted in a FIG 2 Impact of E2 mutations on the fitness of Asian lineage CHIKV in Aedes albopictus. Aedes albopictus mosquitoes were fed artificial blood meals with approximately equal genomic ratios of wt and E2 mutant CHIKV Mal06, one with an introduced ApaI marker and one without. At 10 dpi, heads were collected, and their virus mixture was passaged for 2 days in Vero cells. The resulting cell culture supernatant was amplified by RT-PCR, and the ratio of ApaI-marked to unmarked CHIKV was determined by restriction digestion. Reciprocal experiments were conducted such that the ApaI marker was present in either the mutant (A) or the wt (B) strain. Relative replicative fitness values were calculated as described in Materials and Methods. The Holm-Sidak correction for multiple P values was applied. ns, not significant (P $ 0.05); **, P , 0.01; ****, P , 0.0001.

Lineage-Specific Vector Adaptation of CHIKV ® minimal but negative impact on relative replicative fitness, decreasing fitness by 7-fold or 13-fold depending on whether the marker was in the mutant or the wt strain, respectively. The only mutation to result in a substantial impact in the Mal06 Asian lineage backbone was E2-233E, which decreased fitness by 100-or 50-fold depending on whether the marker was in the mutant or the wt strain, respectively. These data demonstrate that the suite of E2 mutations is not absent from Asian lineages solely because of the inability of the Asian lineage to select E1-226V but also because the E2 mutations themselves are neutral at best and frequently deleterious.

IOL E2 A. albopictus-adaptive mutations in the Asian lineage are neutral or deleterious in A. aegypti. While the IOL E2 mutations were not advantageous to the Mal06 Asian lineage backbone in A. albopictus, their effect in A. aegypti was unknown. Although the IOL E2 mutations conferred no advantage to strain SL07 in A. aegypti (24) , the dominance of A. aegypti in Asian lineage CHIKV circulation merits an investigation of their impact in the context of Mal06 in this vector. As described above, a mixture of wt and E2mutated viruses was fed to mosquitoes in an artificial blood meal. The heads of engorged mosquitoes were collected on day 10, and the differences in wt-to-mutant ratios in the heads compared to the blood meals were measured by Sanger sequencing and used to assess each mutation's impact on fitness for infection and transmission.

Consistent with the results in A. albopictus, the IOL E2 mutations were either neutral or deleterious in the Mal06 backbone in A. aegypti (Fig. 3) . The E2-233E mutation was nearly neutral, resulting in a 2-fold decrease in fitness that did not reach the level of statistical significance. The E2-198Q and E2-252Q mutations, on the other hand, were significantly deleterious, resulting in 9-fold and 50-fold decreases in fitness, respectively. As with the A. albopictus experiments assessed by the restriction digest method to estimate strain ratios, the A. aegypti mosquitoes assessed by the Sanger method overwhelmingly (91 of 108) contained $95% of a single strain.

Genetic differences in IOL and Asian lineage CHIKV strains result in contrasting fitness effects of E2 mutations in A. albopictus. The structures of the glycoproteins (E3, E2, and E1) were examined for potential explanations for the dramatically different fitness effects of the adaptive E2 mutations in IOL versus Asian lineage CHIKV strains. The atomic structure of the IOL glycoproteins, but not the Asian lineage glycoproteins, has been solved (26) . Overall, there are 35 amino acids (aa) that differ in the structural proteins between the two lineages, including 4 in E3 (64 aa total), 17 in E2 (423 aa total), 4 in 6K/TF (61 aa total), and 10 in E1 (439 aa total) ( Fig. 4 ; see also Table S1 in the supplemental material). In silico, we mutated residues from the IOL strain to those in the Asian lineage strain. Based on the predicted structures, all amino acids in E2 that differ between the two strains lie along the top of the E2 protein in domains A and B, suggesting that they likely affect host cell binding. Some of the E1 amino acids that differ between these strains are clustered in the hinge region between domains I and II, which is predicted to move during the fusion process (29, 30) .

We specifically investigated amino acids in the two strains that might interact with adaptive residues in E2 positions 198, 233, and 252. All residues within 4 Å of these E2 positions, and therefore most likely to be involved in direct interactions, were conserved between the two CHIKV strains. When the potential interaction radius was expanded to 20 Å to identify residues that might be interacting during a conformational change during the assembly or fusion process, 8 residues spanning all three envelope glycoproteins were identified (Table S2) .

Only one residue, E3-33, which exists as E3-33E in the IOL (SL07) strain and as E3-33K in the Asian lineage (Mal06) strain, was identified as a site of potential long-range interactions with all four E2 adaptive mutations. In the native IOL strain, E3-33E and E2-233K are 10 Å apart and consist of oppositely charged amino acids, representing a likely attractive interaction. In the presence of the IOL adaptive substitution pair (E3-33E/E2-233E), the presence of two like charges may weaken this interaction in a way that promotes infectivity or assembly. In Asian lineage strains, the wt E3-33K/E2-233K residues may be repelled based on like charges, in a manner similar to that of the adaptive E3-33E/E2-233E IOL strain; in the presence of the introduced adaptive E2 mutation, the Asian pair becomes E3-33K/E2-233E, a potentially attractive interaction similar to the wt IOL strain. Thus, in both the IOL and Asian lineages, repelling charges between E3-33 and E2-233 appear to confer a fitness advantage over attractive charges.

To test this hypothesis, E3-33E was introduced into the Asian Mal06 strains with and without each E2 mutation, and their fitness effects in A. albopictus were compared using the competition assay with Sanger sequencing (Fig. 5) . The E2-198Q mutation, which was neutral in the context of its native E3-33K backbone, was deleterious in the context of E3-33E and resulted in an 8-fold decrease in replicative fitness. The E2-233E mutation, which was highly deleterious in its native E3-33K backbone, was also highly disadvantaged in the context of E3-33E, resulting in a 33-fold decrease in replicative fitness. Finally, E2-252Q, which caused a modest decrease in fitness in its native E3-33K backbone, became neutral in the E3-33E backbone, causing only a 2-fold increase in fitness that failed to attain statistical significance. As with the other competition assays, most mosquito heads (69 of 83) contained $95% of a single strain.

CHIKV has expanded its distribution in recent years through the emergence of the IOL from the ECSA lineage from 2005 to 2006, the expansion of the Asian lineage into (31) . Whereas the evolution of the IOL has been characterized in part by sequential adaptation to a widespread and invasive mosquito vector, A. albopictus, there is no evidence that the Asian lineage has undergone comparable adaptation to any mosquito vector. Given the largely sympatric distributions of A. aegypti and A. albopictus (32) as well as the increasing cocirculation of different CHIKV lineages (33) , understanding the vector usage and the potential for adaptation to expand vector host range, as well as their underlying mechanisms and constraints, is essential for the prevention and control of future epidemics.

We sought to understand why the Asian CHIKV lineage has not expanded its urban vector usage to A. albopictus. Specifically, can the Asian lineage CHIKV strains adopt the same evolutionary pathways as the IOL strains to adapt for A. albopictus transmission? To address this question in a feasible manner, capitalizing on previous work, we initially focused on envelope glycoprotein substitutions previously shown in IOL strains to enhance CHIKV infection of A. albopictus.

First, our study showed that the secondary adaptive E2 mutations observed in different IOL sublineages do not require the presence of the primary E1-A226V substitution for their fitness advantages in A. albopictus in an IOL backbone. Therefore, the adaptation of IOL strains to A. albopictus may be mediated by multiple mechanisms during viral infection and replication. In vitro experiments have shown that the E1-A226V substitution leads to an increased requirement for membrane cholesterol and a lower endocytic pH for fusion in A. albopictus C6/36 cells (34), although it is unclear whether and how these two effects differ between A. albopictus and A. aegypti. The IOL A. albopictus-adaptive E2 substitutions, located in domain B (residues 198 and 210) and the ASR (residues 233 and 252), likely affect the pH-dependent conformational change that exposes the fusion loop on E1. Although the two sets (domain B and ASR) of mutations are both hypothesized to affect the fusion and cell entry of IOL CHIKV strains in A. albopictus, they may have different mechanisms.

Second, our results demonstrate that the IOL and Asian CHIKV lineages have diverged onto different adaptive landscapes for A. albopictus. Unlike the IOL lineage, when we introduced the E2 adaptive mutations (E2-R198Q, E2-K233E, and E2-K252Q) into the Asian lineage strain Mal06, fitness was either unchanged or reduced in A. albopictus. The functions of E2-K233E and E2-K252Q, both located in the E2 ASR, could include interactions with cellular factors specific to A. albopictus, which may regulate the structural change of the ASR under low endosomal pH. In IOL strains, these interactions are advantageous, but in Asian strains, these interactions may be abolished or constrained.

These adaptive E2 mutations in the IOL strains were previously shown to have little or no effect on viral fitness in A. aegypti (24) . Our data revealed that these mutations are either neutral or detrimental in the context of Asian lineage CHIKV infection of A. aegypti. Specifically, in the Mal06 backbone, E2-198Q was neutral in A. albopictus and moderately deleterious in A. aegypti, E2-233E was strongly deleterious in A. albopictus and neutral in A. aegypti, and E2-252Q was moderately deleterious in A. albopictus and strongly deleterious in A. aegypti. These results confirm that residues E2-198, 2233, and 2252 interact with Aedes mosquitoes in a species-specific manner. Research to date on CHIKV entry into mosquito cells is largely based on ECSA/IOL strains in cultured A. albopictus C6/36 cells (34) (35) (36) (37) . Comparative studies between the multiple lineages and in both A. albopictus and A. aegypti cell lines could enhance the understanding of CHIKV-vector interactions.

Our studies of the E3-33 substitution suggest that, despite its position as the only residue to both differ between SL07 and Mal06 and to be located within 20 Å of all four E2 adaptive residues, it is not responsible for the contrasting fitness effects of these E2 mutations on IOL or Asian lineage CHIKV. Although E3-33 was considered the most likely candidate based on proximity to all E2 adaptive mutations in the mature glycoprotein spike, it is possible that the change(s) responsible for the divergent SL07 and Mal06 phenotypes is more pertinent to viral assembly or fusion and thus is not reflected in the conformation of a mature virion. It is also possible that multiple mutations jointly shape the conformation of the protein(s) during different stages of replication or influence its interaction with cellular factors in A. albopictus. None of the E2 residues that we studied contact the Mxra8 receptor identified for several alphaviruses (38) . However, this receptor occurs only in vertebrates, so CHIKV interactions with another receptor could be involved in mosquito infection. Regardless, further work is required to fully define the differentiating factors that confine Asian lineage and IOL CHIKV to different adaptive landscapes.

There are several limitations of our studies. We focused only on mutations identified previously in IOL strains as they spread through Asia following the 2005 emergence. However, it is possible that other genome regions have also contributed to CHIKV fitness for transmission by A. albopictus, so further studies starting with chimeric viruses mixing other genes among early-and late-stage IOL strains as well as Asian strains may yield further insights.

In conclusion, our results suggest that the IOL and Asian lineage CHIKVs have diverged onto different evolutionary trajectories and fitness landscapes and no longer share the same adaptability for the two urban vectors. Given that CHIKV had no known exposure to A. albopictus in Africa prior to 2006 (39), it is unclear whether its ability to readily adapt to this mosquito is a consequence of its ancestral state. It is also unclear why the Asian lineage has lost or never gained the ability to use A. albopictus as a major vector in Southeast Asia, where this mosquito is native and abundant. While we cannot exclude the possibility that Asian lineage CHIKV strains could adapt to A. albopictus through other mechanisms yet to be realized, our results strongly indicate that any such adaptation must arise through a different mutation or series of mutations than those found in IOL strains.

Cell culture. Vero African green monkey kidney cells (ATCC CCL-81; American Type Culture Collection, Manassas, VA) were maintained in Dulbecco's modified Eagle's essential medium (DMEM) supplemented with 5% fetal bovine serum (FBS) and penicillin/streptomycin at 37°C under 5% CO 2 . C7/ 10 A. albopictus cells were maintained in DMEM supplemented with 10% FBS, 5% tryptose phosphate broth, and penicillin/streptomycin at 28°C under 5% CO 2 .

Construction of infectious cDNA clones. Plasmids encoding two wild-type CHIKV strains, SL07 and Mal06, from the IOL and Asian lineage, respectively, were described previously (27) . Both wt plasmids contain their ancestral E1-226A codons. Single mutations representing all four of the E2 substitutions (E2-R198Q, E2-L210Q, E2-K233E, and E2-K252Q) were introduced into the SL07 plasmid, and three of the Lineage-Specific Vector Adaptation of CHIKV ® November/December 2021 Volume 12 Issue 6 e02738-21 mbio.asm.org 9

E2 substitutions (E2-R198Q, E2-K233E, and E2-K252Q) in either the presence or absence of an additional E3-K33E substitution were introduced into the Mal06 plasmid using conventional PCR-based methods (40) . In addition, to compare relative fitness levels using a restriction digest-based competition test, a synonymous point mutation was introduced into the nsP4 gene in each plasmid, including all of the wt strains and mutants, to form restriction sites cleavable by the endonucleases ApaI and PspOMI (27) . All PCR-generated genome regions used for cloning were Sanger sequenced to verify their integrity. Plasmids were purified by centrifugation in cesium chloride gradients.

Rescue of viruses from infectious clones. To generate infectious RNA, plasmids were linearized with the NotI restriction endonuclease, followed by in vitro transcription from the SP6 promoter as described previously (40) . Approximately 10 mg of RNA was electroporated into C7/10 cells. Cell culture supernatants were harvested at 24 to 48 h postelectroporation and stored at 280°C. Infectious viral titers were determined by plaque assays on Vero cells (titers ranged from 10 7 to 10 8 PFU/ml). Viruses recovered from electroporated cells were used directly without additional passages.

Mosquito feeds and harvests. To determine the relative fitness levels of mutant versus wt CHIKV strains, competition tests were conducted in mosquitoes as described previously (27) . To achieve the approximately 1:1 ratio desired for the competition assay, wt and mutant viruses were initially combined at a 1:1 ratio as determined by Vero cell PFU. This ratio, as well as a number of ratios closely flanking it, was further refined to achieve a 1:1 genome copy ratio because the final output assays, whether restriction digest based or PCR based, are all genome-based assays. In the case of experiments with restriction digests as their final output, this adjustment was based on the brightness of reverse transcription-PCR (RT-PCR) amplicon bands. In the case of experiments with Sanger sequencing as their final output, the ratio of the two viruses was estimated based on the single nucleotide polymorphism (SNP) heights in the electropherogram generated by Sanger sequencing. Initial ratios measured genetically varied from 1:1.00 to 1:3.04 (mean = 1:1.40; median = 1:1.17).

The competing viruses were mixed at the determined 1:1 genomic ratio with defibrinated sheep blood to a final concentration of 1 Â 10 6 to 1 Â 10 7 PFU/ml. The blood meal was offered for approximately 1 h to A. aegypti or A. albopictus mosquitoes from colonies established from mosquito eggs collected in Bangkok, Thailand. This colony was chosen to correspond to a location where selection for adaptation to both mosquito vectors has presumably occurred for decades. Fully engorged mosquitoes were incubated at 28°C (a typical mean tropical temperature in locations where CHIKV is endemic) with 10% sucrose and 80% relative humidity under a 16-h-light/8-h-dark photoperiod.

Mosquitoes were harvested at 10 days postinfection (dpi). The heads were separated from the bodies and homogenized in 500 ml of DMEM. Heads were selected rather than salivary glands to permit the efficient harvest of large numbers of samples. As salivary gland infection and escape barriers have not been described for CHIKV and none of the mutants studied appear to affect salivary gland replication (41) , the CHIKV mixture present in the head is assumed to accurately reflect the CHIKV mixture present for transmission in the saliva.

To determine if these mosquito heads were infected, homogenized samples were passaged in Vero cells prepared in 96-well plates. Supernatants from Vero cells showing cytopathic effect (CPE) were collected at 2 dpi and used for determining the ratio of the competing viruses. Analysis of the Vero cell-passaged ratios was undertaken to screen for infected heads more efficiently than using the more complex and less sensitive RT-PCR. Also, amplification on Vero cells minimized the risk of direct RT-PCR amplicons from saliva being near the limit of detection and thus representing a very small number of CHIKV genomes sampled, compared to the greater sensitivity of Vero cells to amplify more viral genomes. Passaging known mixtures of CHIKV wt and mutant strains in Vero cells in duplicate produced minimal differences in input versus output ratios, and none of the differences were statistically significant whether considering individual strain combinations (P value range of 0.617811 to 0.904460 by repeatedmeasures one-way analysis of variance [ANOVA] comparing the input and duplicate passage outputs at three ratios per strain with a Holm-Sidak correction for multiple comparisons) or the combined data set (P = 0.300570) (see Fig. S1 in the supplemental material).

Determination of competition outcomes. Viral RNA was purified from each passaged mosquito head as well as from the initial blood meal using the QIAamp viral RNA kit (Qiagen). Ratios of wild-type and mutant viruses were determined by one of two methods: restriction digestion or Sanger sequencing of RT-PCR amplicons.

For samples evaluated by restriction digestion, the genome region from nucleotides (nt) 6106 to 6794, which includes in its center the introduced ApaI/PspOMI restriction site marker, was amplified by RT-PCR using the Qiagen OneStep RT-PCR kit with the primers 41855-nsF5 (59-ATATCTAGACATGGTGGA-39) and 41855-nsR1 (59-TATCAAAGGAGGCTATGTC-39). The amplicons were digested with the restriction endonucleases ApaI and PspOMI (which are isoschizomers sharing the recognition site GGGCCC) for 30 min at 27°C followed by 4 h at 37°C. Specificity and complete digestion were confirmed using controls that contained only one of the competing viruses (with or without the marker). Band ratios were determined using GelQuant.NET software (Biochem Lab Solutions). Gel images, along with the corresponding band ratios, are available in Fig. S2 . Validation of the restriction digest method is available in Fig. S3 .

Because control experiments demonstrated that the restriction marker mutation described above unexpectedly affected fitness in A. aegypti (it had been shown to have no phenotype in A. albopictus), we switched to using a Sanger sequencing method to estimate mutant/wild-type ratios. This method avoids potential effects on the fitness of the engineered restriction site marker and has been shown in other studies to be highly accurate and reproducible (42) (43) (44) 49) . The regions containing the mutations of interest were amplified by RT-PCR using the primers Ch18(1) 8730 (59-TTGGACYAAGCTGCG-39) and Ch19(2) 9419 (59-GTCGGATGGTCAGGATACAG-39). The amplicons were subsequently purified by a Qiagen PCR purification kit and sequenced using a BigDye Terminator v3.1 cycle sequencing kit and an ABI Prism model 3100 genetic analyzer (Applied Biosystems, Foster City, CA). The ratio of the two viruses was estimated based on the SNP heights in the sequencing electropherogram using QSVanalyser software (45) . This method has been shown to be highly accurate and reproducible for RNA viruses (42) (43) (44) . Validation of the Sanger method is available in Fig. S4 .

Using these assays, the frequencies of minority genotypes below an expected frequency of 10% may be observed with variation as high as 20%. However, the variability of accurately measuring low-frequency genotypes does not impact our determination of which genotype is dominant in a sample.

Molecular modeling. The structure of the immature envelope glycoprotein complex of a CHIKV IOL strain (PDB accession number 3N40) was used in PyMOL (PyMOL molecular graphics system version 1.8; Schrödinger, LLC, New York, NY) to determine amino acid distances. Figures were generated in UCSF ChimeraX version 1.2.5 (46, 47) .

Statistical analysis. Relative replicative fitness was calculated as previously described (42, 43) . Briefly, replicative fitness values were modeled as w = f 0 /i 0 , where i 0 is the initial ratio of a given strain in the blood meal and f 0 is the final ratio of that strain in the passaged mosquito head. The model was fitted to log 10 (mosquito head output/blood meal input) ; strain, with log 10 conversion performed to improve the normality of the data set. Relative replicative fitness values were determined by the model's strain coefficient, which was transformed to the original scale as 10 coefficient . When the band ratios following restriction digestion were 100 to 0, the ratios were considered 99.9 to 0.1 for the calculation of relative replicative fitness to avoid dividing by 0. Similarly, if one of the peak heights following Sanger sequencing was 0, it was changed to 0.1 to avoid dividing by 0 (in this instance, the peak height of the other strain was not changed). The Holm-Sidak method was used to correct for multiple comparisons when considering the P values generated by the relative replicative fitness calculations.

Supplemental material is available online only. 

We thank Jing Huang for excellent technical assistance. This research was supported by NIH grants AI121452 and R24 AI120942.

Present and future arboviral threats

Zika virus

Global temperature constraints on Aedes aegypti and Ae. albopictus persistence and competence for dengue virus transmission

Chikungunya virus and the global spread of a mosquito-borne disease

Genome-scale phylogenetic analyses of chikungunya virus reveal independent emergences of recent epidemics and various evolutionary rates

Chikungunya in the Americas

Emergence and potential for spread of chikungunya virus in Brazil

Chikungunya virus sequences across the first epidemic in Nicaragua

Chikungunya virus as cause of febrile illness outbreak

Chikungunya virus transmission found in the United States: US health authorities brace for wider spread

DSHS announces first Texas-acquired chikungunya case

Reemergence of chikungunya virus

First report of Aedes aegypti transmission of chikungunya virus in the Americas

Arbovirus surveillance and first report of chikungunya virus in wild populations of Aedes aegypti from Guerrero, Mexico

Lineage-Specific Vector Adaptation of CHIKV ®

Droughtassociated chikungunya emergence along coastal East Africa

Genome microevolution of chikungunya viruses causing the Indian Ocean outbreak

A single mutation in chikungunya virus affects vector specificity and epidemic potential

Two chikungunya isolates from the outbreak of La Reunion (Indian Ocean) exhibit different patterns of infection in the mosquito, Aedes albopictus

Genetic divergence of chikungunya viruses in India (1963-2006) with special reference to the 2005-2006 explosive epidemic

Genetic characterization of 2006-2008 isolates of chikungunya virus from Kerala, South India, by whole genome sequence analysis

Appearance of E1: A226V mutant chikungunya virus in coastal Karnataka, India during 2008 outbreak

evolutionary dynamics and serotype analyses

Comprehensive genome scale phylogenetic study provides new insights on the global expansion of chikungunya virus

Multi-peaked adaptive landscape for chikungunya virus evolution predicts continued fitness optimization in Aedes albopictus mosquitoes

Sequential adaptive mutations enhance efficient vector switching by chikungunya virus and its epidemic emergence

Glycoprotein organization of chikungunya virus particles revealed by X-ray crystallography

Chikungunya virus emergence is constrained in Asia by lineagespecific adaptive landscapes

Epistatic roles of E2 glycoprotein mutations in adaption of chikungunya virus to Aedes albopictus and Ae. aegypti mosquitoes

A ligand-binding pocket in the dengue virus envelope glycoprotein

East/Central/ South African genotype chikungunya virus, Brazil

The global distribution of the arbovirus vectors

Zika, chikungunya, and other emerging vector-borne viral diseases

pH-dependent entry of chikungunya virus into Aedes albopictus cells

Mosquito cellular factors and functions in mediating the infectious entry of chikungunya virus

Dynamics of chikungunya virus cell entry unraveled by single-virus tracking in living cells

Early events in chikungunya virus infection-from virus cell binding to membrane fusion

Cryo-EM structure of chikungunya virus in complex with the Mxra8 receptor

The invasive mosquito species Aedes albopictus: current knowledge and future perspectives

Infectious clones of chikungunya virus (La Réunion isolate) for vector competence studies

Tissue barriers to arbovirus infection in mosquitoes

Role of mutational reversions and fitness restoration in Zika virus spread to the Americas

Spike mutation D614G alters SARS-CoV-2 fitness

Submergence" of Western equine encephalitis virus: evidence of positive selection argues against genetic drift and fitness reductions

Inferring relative proportions of DNA variants from sequencing electropherograms

UCSF ChimeraX: meeting modern challenges in visualization and analysis

UCSF ChimeraX: structure visualization for researchers, educators, and developers

Population bottlenecks and founder effects: implications for mosquito-borne arboviral emergence

The N501Y spike substitution enhances SARS-CoV-2 infection and transmission. Nature