The dinucleotide composition of the Zika virus genome is shaped by conflicting evolutionary pressures in mammalian hosts and mosquito vectors

Most vertebrate RNA viruses show pervasive suppression of CpG and UpA dinucleotides, closely resembling the dinucleotide composition of host cell transcriptomes. In contrast, CpG suppression is absent in both invertebrate mRNA and RNA viruses that exclusively infect arthropods. Arthropod-borne (arbo) viruses are transmitted between vertebrate hosts by invertebrate vectors and thus encounter potentially conflicting evolutionary pressures in the different cytoplasmic environments. Using a newly developed Zika virus (ZIKV) model, we have investigated how demands for CpG suppression in vertebrate cells can be reconciled with potentially quite different compositional requirements in invertebrates and how this affects ZIKV replication and transmission. Mutant viruses with synonymously elevated CpG or UpA dinucleotide frequencies showed attenuated replication in vertebrate cell lines, which was rescued by knockout of the zinc-finger antiviral protein (ZAP). Conversely, in mosquito cells, ZIKV mutants with elevated CpG dinucleotide frequencies showed substantially enhanced replication compared to wild type. Host-driven effects on virus replication attenuation and enhancement were even more apparent in mouse and mosquito models. Infections with CpG- or UpA-high ZIKV mutants in mice did not cause typical ZIKV-induced tissue damage and completely protected mice during subsequent challenge with wild-type virus, which demonstrates their potential as live-attenuated vaccines. In contrast, the CpG-high mutants displayed enhanced replication in Aedes aegypti mosquitoes and a larger proportion of mosquitoes carried infectious virus in their saliva. These findings show that mosquito cells are also capable of discriminating RNA based on dinucleotide composition. However, the evolutionary pressure on the CpG dinucleotides of viral genomes in arthropod vectors directly opposes the pressure present in vertebrate host cells, which provides evidence that an adaptive compromise is required for arbovirus transmission. This suggests that the genome composition of arbo flaviviruses is crucial to maintain the balance between high-level replication in the vertebrate host and persistent replication in the mosquito vector.

Introduction Arthropod-borne viruses (arboviruses) are transmitted between vertebrate hosts by invertebrate arthropod vectors such as mosquitoes, ticks, biting midges, and sand flies. When hematophagous vectors feed on infected vertebrate hosts, the virus first needs to enter and pass through the vector's midgut endothelial cells and persistently infect other tissues to disseminate into the vector's salivary glands before the arbovirus is transmitted to the next vertebrate host during subsequent feeding events. Transmission of arboviruses thus requires active virus replication in 2 evolutionary distant hosts that are separated by more than 500 million years of evolution [1]. Hence, arboviruses encounter divergent cellular and molecular environments, exemplified by the different antiviral responses that vertebrate hosts and invertebrate vectors rely on to cope with viral infections. Vertebrate animals possess a plethora of innate and adaptive antiviral immune responses [2], while invertebrates largely rely on RNA interference (RNAi) to control viral infections [3,4]. As intracellular pathogens, all viruses rely on cellular machinery for their replication, which exposes viruses to the intracellular milieu of the host cell. The transmission cycle of arboviruses places them in 2 distinct intracellular environments with dissimilar adaptive pressures for replication and evasion of host antiviral responses (reviewed in [5]). How arboviruses are able to thrive while continuously alternating between these distinct cellular environments remains an intriguing research topic [6].
One of the most striking molecular differences between vertebrate and invertebrate organisms is the dinucleotide usage in their genomic DNA and coding mRNA sequences. In the DNA of vertebrates, the dinucleotide combination of a cytosine followed by a guanine (CpG) occurs less frequently than predicted, based on frequencies of C and G mononucleotides. The cytosine in CpG is prone to methylation, which can cause the cytosine to deaminate and mutate into a thymine [7,8]. In contrast, the genome of most invertebrates contains largely unbiased frequencies of CpG dinucleotides, reflecting greatly reduced DNA methylation activity [9]. Both vertebrates and invertebrates have additionally evolved to reduce UpA dinucleotide frequencies in their mRNA through a range of mechanisms largely related to mRNA stability [10,11]. Most viruses with RNA genomes have evolved to mirror the dinucleotide usage of their host's mRNA, with invertebrate viruses solely suppressing UpA dinucleotides and vertebrate viruses underrepresenting both UpA and CpG dinucleotides, even though CpGs in viral RNA are not subjected to the same methylation and mutation pressure that acts on vertebrate host DNA [12][13][14].
Whether invertebrate cells possess similar non-self RNA detection mechanisms and whether alternating between vertebrate hosts and arthropod vectors puts additional selective pressure on arboviruses is currently unknown. The lack of CpG suppression in invertebrates and their viruses suggests the absence of selection pressure against CpGs in arthropods, while UpA dinucleotides are clearly underrepresented in both arthropods and their viruses [13,14]. We therefore hypothesized that the selection pressures on UpA dinucleotides would be similar in vertebrate hosts and arthropod vectors, while the selection pressure on CpG dinucleotides may be absent in these vectors. To investigate the constraints that an arbovirus encounters during its transmission cycle, we turned to the genus Flavivirus (family Flaviviridae). This genus contains many clinically relevant arboviruses, which include mosquito-borne flaviviruses (MBFV), e.g., Zika virus (ZIKV), dengue virus, West Nile virus, Japanese encephalitis virus, and yellow fever virus, as well as a number of tick-borne flaviviruses (TBFV), e.g., tickborne encephalitis virus. In addition to arboviruses, the genus Flavivirus also contains members with no known vector (NKVF) and insect-specific flaviviruses (ISF). The latter are unable to infect vertebrates and exclusively replicate in insects and cultured insect cells [31]. To investigate whether arbovirus genome composition is subject to distinct evolutionary pressures in vertebrate host and arthropod vectors, we constructed ZIKV mutants with synonymously increased CpG or UpA dinucleotide frequencies and tested their replication kinetics in vertebrate and invertebrate cells, mice, and mosquitoes.

Flavivirus genome composition and mutant Zika virus design
To compare the CpG and UpA dinucleotide usage between different flaviviruses, the ratio of CpG and UpA dinucleotide frequencies from full-length flavivirus genome sequences were plotted as the observed over the expected (O/E) frequencies ( Fig 1A, large data points). Viral sequence data were superimposed onto large randomly selected sets of vertebrate (Homo sapiens) and invertebrate mosquito (Aedes aegypti) mRNAs ( Fig 1A, small data points). Both human host and mosquito vector mRNA suppress UpA dinucleotides (y-axis), while only human mRNA underrepresents CpG dinucleotides (x-axis). All flaviviruses that infect vertebrates (i.e., MBFV, TBFV, and NKVF) have clearly adapted to mimic the genome composition of the vertebrate host, with marked underrepresentation of both CpG and UpA dinucleotides. In contrast, the ISFs and particularly the phylogenetically most distant lineage 1 ISFs contain a much higher ratio of CpG dinucleotides frequencies, more closely resembling the unbiased CpG frequency found in their arthropod hosts. Lineage II ISFs display an insect-specific phenotype and contain more CpG dinucleotides than VIFsAU : PleasedefineVIFsinitsfirstmentioninthesent , while phylogenetically they cluster together with the VIFs (Fig 1A) [14,[31][32][33][34].
To investigate the biological role of CpG and UpA suppression in a typical MBFV, ZIKV mutants were created with increased frequencies of either CpG or UpA dinucleotides by introducing many synonymous changes into a 1654 nucleotide long, coding region of the viral genome ( Fig 1B). The complete genome of wild-type (WT) ZIKV has ratios of 0.42 and 0.52 for the O/E frequencies of CpG and UpA dinucleotides, respectively. The selected region contains 53 CpG and 57 UpA dinucleotides with comparable CpG and UpA suppression (detailed in S1 Table). Before introducing any mutations, we confirmed that the ZIKV open reading frame (ORF), except for the 5 0 -end which contains the cyclisation sequence [35], has a relatively uniform nucleotide variability at synonymous sites and low minimal free energy differences (MFED) (Fig 1C). This suggests that, apart from conserving amino acid sequence, there Small data points denote human mRNAs (gray) and mosquito mRNAs (black). Large data points represent full-length genome sequences of flavivirus species with MBFV (green), TBFV (light green), NKVF (blue), and the lineage 1 and 2 ISF (red and orange), respectively. The O/E ratio of 1.0 (dotted lines) is the mathematically expected frequency of occurrence if the RNA's mononucleotides were randomly distributed. (B) Schematic representation of the ZIKV genome, with the mutated region boxed and displaying the number of CpG and UpA dinucleotides in wild-type ZIKV. Additional CpG and UpA dinucleotides of the CpG_1.0 (+59), CpG_max (+190), and UpA_max (+150) mutant viruses are listed in brackets. (C) MFED and synonymous pairwise distance (black line) were calculated for 19 ZIKV isolates that were >1% dissimilar in nucleotide sequences using a fragment size of 200 with 48 nucleotide increments for MFED and sliding window of 21 codons for synonymous pairwise distance that were >1% dissimilar in nucleotide sequences using a fragment size of 200 with 48 nucleotide increments for MFED and sliding window of 21 codons for synonymous pairwise distance as used in are no other conserved nucleotide sequences or critical RNA structures in the ZIKV ORF. To confirm the absence of critical RNA sequences experimentally, a scrambled control virus (SCR) was designed that contained extensive synonymous changes in the target region but retained WT dinucleotide frequencies. Furthermore, we created mutants that maximally overrepresented CpG or UpA dinucleotides, without altering the number of the other UpA or CpG dinucleotide or the encoded amino acid sequence. This resulted in the CpG_max ZIKV mutant with 190 additional CpG dinucleotides and an O/E ratio of 1.62 and the UpA_max ZIKV mutant (+ 150 UpA dinucleotides, O/E ratio 1.72). As the genome of Ae. aegypti, the main mosquito vector for ZIKV, shows no suppression of CpG dinucleotides (O/E ratio approximately 1.0) (Fig 1A), we additionally created a mutant virus with synonymously increased CpG dinucleotides to the unsuppressed (mosquito-like) ratio of 1.0 (termed CpG_1.0) by addition of 59 CpG dinucleotides. We made use of an infectious cDNA clone to generate virus stocks of these ZIKV mutants in Vero E6 cells. As the fitness between the ZIKV mutants was expected to vary substantially, we used equal numbers of viral RNA copies to infect cells.

ZIKV with elevated CpG or UpA dinucleotide frequencies is attenuated in vertebrate cell cultures
Vertebrate cell lines were infected with the different mutant ZIKV stocks to investigate the replication kinetics of the mutant viruses. Samples were taken daily from cell culture supernatants and checked for the presence of infectious virus (Fig 2). In all cases, the SCR control behaved very similarly to the WT, indicating that the introduced synonymous changes did not affect virus replication. In vertebrate Vero E6 cells and A549 cells, UpA_max ZIKV was clearly attenuated with approximately 10-fold lower titers of progeny virus at any time postinfection. The CpG_max ZIKV mutant displayed 2-to 3-fold lower titers compared to WT, while the replication kinetics of the CpG_1.0 mutant virus were similar to that of the WT virus (Fig 2A and  2B). At 2 days postinfection, CpG_max and UpA_max virus differed significantly from WT and the SCR control ( Fig 2D) (two-way ANOVA, F (4, 30) = 6.315, p < 0.05). Knockout of ZAP rescued the attenuated phenotype of the CpG_max ZIKV mutant, compared to that in wild-type A549 cells (Fig 2C and 2D) (two-way ANOVA, F (1, 30) = 11.49, p < 0.05). Although replication of the UpA_max ZIKV mutant appeared similarly rescued by knockout of ZAP ( Fig 2B-2D), this was not significantly different from replication in wild-type A549 cells (p = 0.06). Similar to what was reported previously for CpG-high E7 mutants [21], knockout of ribonuclease L (RNaseL) or OAS3 resulted in comparable replication rates between CpG_max and WT ZIKV to WT levels (Fig 2 and S1A-S1C Fig).
The cytopathic effects (CPE) caused by ZIKV infections in A549 cells were mild; however, delays in CPE for the CpG_1.0, CpG_max, and UpA_max viruses were observed. To investigate replication phenotypes further, a plaque assay was performed and plaque sizes were determined by immunostaining as a measure for the rate at which infectious virus spread through the culture. Three days postinoculation with WT ZIKV and SCR ZIKV, large plaques with an average diameter just over 1,000 μm were observed in parental A549 cells (Fig 2E and S1D  Fig). Inoculation with CpG_1.0 ZIKV resulted in slightly smaller average plaque diameters (approximately 750 μm), while the CpG_max and UpA_max ZIKV mutants displayed previous analyses of echovirus 7 and murine norovirus [19,59]. MFED values of both the sense (red line) and the antisense (dashed pink line) strand in the selected region are close to the null expectation (dotted black line). Please refer to S1 Data for the numerical values underlying panels A and C. Values were calculated using SSE software v1. 4 significantly smaller average plaque sizes (530 μm and 200 μm, respectively) ( Fig 2E and S1D  Fig) (one-way ANOVA, F (4, 23) = 19.06, p < 0.05). The same virus dilutions were used to infect A549 ZAP knockout cells in parallel. In these cells, plaques observed in response to WT ZIKV infection were smaller compared to WT infections in ZAP expressing A549 cells. However, no further reduction in plaque sizes was observed when cells were infected with CpG-or UpA-high mutant ZIKV (Fig 2DAU : PleasenotethatthefigurecitationsFig2EandgDinthesentenceHowev and 2E) (one-way ANOVA, F (4, 19) = 1.020, p = 0.42). Plaque phenotypes were confirmed in Vero E6 in which also the average plaque diameters of the CpG_1.0 mutant were significantly smaller than those of WT (S1E Fig) (F (4, 25) = 65.29, p < 0.05). The analysis of plaque diameters suggests that ZAP is responsible for the attenuation of both CpG_max and UpA_max mutant viruses. The results further suggest that the CpG_1.0 mutant is mildly attenuated in vertebrate cells that express ZAP, something that was not visible from the virus growth curves (Fig 2). Single asterisks highlight significant differences in WT vertebrate cells from WT ZIKV and 2 asterisks the difference between A549 cells and A549 ZAP knockout cells (two-way ANOVA, FDR adjusted p < 0.05). (E) Average plaque size measured in A549 and A549 ZAP knockout cells 3 dpi with ZIKV or indicated mutant viruses. After removal of the overlay, cells were fixed and stained for ZIKV E protein. Micrographs were taken from plaques from 2 independent experiments. The plaque diameter was measured using the Axiovision software. Data points represent the average plaque size, error bars indicate 1 standard error of the mean, and asterisks highlight significant difference from WT (one-way ANOVA with Dunnett's post hoc test, p < 0.05). Please refer to S1 Data for the numerical values underlying this figure. dpi, days postinfection; FDR, false discovery rate; ko, knockout; SCR, scrambled control virus; WT, wild-type; ZAP, zincfinger antiviral protein; ZIKV, Zika virus. https://doi.org/10.1371/journal.pbio.3001201.g002

Replication of ZIKV mutants is differentially affected by genomic CpG and UpA dinucleotide frequencies in mosquito cells
Cell lines originating from Aedes spp. mosquitoes were infected with WT ZIKV and the mutant viruses to investigate how changes in CpG and UpA dinucleotide frequencies affect ZIKV replication in invertebrate cells. The replication kinetics of WT ZIKV and the SCR control were highly similar. Interestingly, both mutants with increased CpG dinucleotide frequencies ( Fig 3A and 3B and S2A-S2C Fig, circles) displayed an increased fitness compared to WT in all Aedes mosquito cells tested, with a most consistent difference in the exponential growth phase at 2 days postinfection. Two days postinfection, the mean virus titers of the CpG_1.0 and the CpG_max mutants across all 5 cell lines were 6-and 5-fold higher than WT, respectively (two-way ANOVA, F (4, 45) = 4.120, p = < 0.05). In contrast, virus titers of the UpA_ max mutant were lower than WT in all cell lines tested; however, this did not significantly differ from WT ( Fig 3C). To confirm the opposing replicative fitness of the CpG-high mutant viruses in either vertebrate or mosquito cells in a more sensitive assay, a competitive fitness experiment was performed. Human A549 cells, A549 ZAP knockout cells, and AP-61 mosquito cells were infected with 2 ZIKV variants. WT ZIKV was mixed with either the SCR, CpG-high viruses (CpG_1.0 and CpG_max), or UpA_max at equal RNA concentrations. Virus was passaged to new cells twice before total RNA was isolated and part of the mutated region within the ZIKV genome was reverse transcribed and PCR amplified. Amplicon DNA was digested to identify WT and mutant viruses. As expected, WT and SCR viruses were both detected in all 3 cell types ( Fig  4A). In A549 cells, the CpG_1.0 mutant was almost eliminated, while the CpG_max mutant was completely outcompeted by WT (Fig 4B and 4C). In contrast, in AP-61 mosquito cells, more CpG_max ZIKV was detected compared to WT and the CpG_1.0 mutant had almost completely outcompeted WT (Fig 4B and 4C). The UpA_max mutant was outcompeted by WT in A549 cells, while in mosquito cells, it displayed fitness equal to that of WT. Knockout of ZAP strongly reduced any fitness advantages of WT over the mutants in A549 cells as all mutants were readily detected together with WT ( Fig 4A-4D). Together with the virus growth curves, this indicates that while vertebrate cells attenuate the replication of viruses with increased CpG dinucleotide frequencies, these viruses have a replication advantage in mosquito cells.

CpG-high ZIKV is more infectious in Aedes aegypti mosquitoes
To investigate the biological role of dinucleotides frequencies in arbovirus RNA in living invertebrate vectors, Aedes aegypti mosquitoes were orally infected with WT ZIKV and the mutant viruses. Mosquitoes were fed an infectious blood meal containing a fixed number (10 8 / ml) of RNA copies, which for WT ZIKV equalled 5 × 10 5 TCID 50 /ml ( Fig 5). Engorged females were selected and incubated for 14 days before mosquito saliva was obtained. Whole body mosquito homogenates and saliva samples were scored for the presence of ZIKV ( Fig 5A). Strikingly, 70% of the mosquito bodies became infected after ingesting CpG_1.0 ZIKV and 76% after ingesting the CpG_max ZIKV mutant, while this was only between 26% and 37% for mosquitoes that had ingested WT, SCR, or UpA_max viruses (Fig 5B).
Between 13% and 16% of mosquitoes carried infectious WT, SCR, or UpA_max ZIKV in their saliva. Significantly higher percentages of ZIKV-positive saliva were detected among mosquitoes that fed on blood containing the CpG_1.0 (50%) and CpG_max (52%) ZIKV mutants (Fisher exact test, p < 0.01) (Fig 5C). This indicates that dissemination into the saliva also occurred more frequently when mosquitoes were infected with the CpG-high viruses. Interestingly, a notable mortality of mosquitoes that fed on the CpG_max ZIKV mutant was observed during the experiment (log-rank test, p < 0.001) (Fig 5D). Despite the different infection rates, whole body homogenates of WT, SCR, CpG_1.0, and CpG_max ZIKV-infected mosquitoes contained similar viral titers (approximately 10 7 TCID 50 /ml) ( Fig 5E). However, mosquitoes infected with UpA_max ZIKV displayed 10-fold lower viral titers (Mann-Whitney U, p < 0.05) (Fig 5E). In an independent experiment in which mosquitoes were fed blood meals containing equal infectious virus titers (5 × 10 5 TCID 50 /ml) instead of equal RNA copy numbers, mosquitoes were again more often infected with the ZIKV mutants with elevated CpG frequencies, compared to WT and SCR ZIKV (S4 Fig) (Fisher exact test, p < 0.01).

ZIKV mutants with elevated CpG and UpA dinucleotide frequencies show attenuated infection in IFNAR −/− mice
To investigate the behavior of ZIKV mutants in a mammalian host, groups of adult male IFNAR −/− mice (defective in the interferon-α/β receptor) were infected with the WT and ZIKV mutants or were injected with PBS as a negative control (Fig 6A). Weight loss and viremia were determined for each animal. During the first 5 days, mice infected with WT or SCR ZIKV lost 6% to 8% of their total body weight ( Fig 6B). Infections with the CpG_1.0 and CpG_max ZIKV mutants resulted in a 5% to 6% loss of total body weight. HoweverAU : Pleasecheckwh , the mice infected with the CpG_1.0 and CpG_max ZIKV mutants recovered faster than mice infected with WT ZIKV and reached their respective initial body weights around days 23 and 29 postinfection compared to 32 dpi for mice infected with WT ZIKV. Mice that were inoculated with UpA_max ZIKV or PBS showed no signs of disease and maintained a healthy body weight throughout the experiment ( Fig 6B). Next, blood viremia levels in mice infected with both mutants that displayed reduced weight loss were compared to those infected with WT ZIKV. Mice infected with the CpG_max and UpA_max ZIKV mutants still developed substantial viremia, which was detected for 7 days postinfection. However, at peak viremia, titers were on average 5-fold lower than for WT ZIKV (Fig 6C, day 3). All infections resulted in strong anti-ZIKV immunoglobulin G (IgGAU : PleasenotethatIgGhasbeendefinedasimmunoglobulinGinitsfirs ) antibody responses in week 4, though responses in CpG_max and UpA_max infected mice led to mildly lower (approximately 2-fold) end-point ELISA titers compared to infections with WT, SCR, and CpG_1.0 ZIKV (Fig 6D) (one-way ANOVA, F (4, 20) = 6.917, p < 0.05).

Challenge with ZIKV PRVABC59
The mice were subsequently challenged with ZIKV PRVABC59 on day 32 after the immunizing infection, with the ZIKV mutants (see Fig 6A). All the previously infected mice were protected against ZIKV PRVABC59 -induced weight loss ( Fig 7A). Only those mice that were mock infected with PBS showed signs of weight loss and developed ZIKV viremia after the challenge (Fig 7A  and 7B). In this IFNAR −/− mouse model, ZIKV infection results in testicular damage, clearly seen as substantial overt reductions in testis size. The unvaccinated PBS control group displayed a marked reduction in testis size after challenge with ZIKV PRVABC59 . Testis size reductions were also observed in mice immunized with WT ZIKV and the SCR and CpG_1.0 ZIKV mutants, likely as a result of the first infection rather than due to ZIKV PRVABC59 challenge, given that all mice seroconverted after the initial infection ( Fig 6D). Mice that were first infected with the CpG_max and UpA_max ZIKV mutants retained normal testis size after challenge with ZIKV PRVABC59 (Fig 7C) (Kruskal-Wallis test, p < 0.05), illustrating that neither the first infection nor ZIKV PRVABC59 challenge resulted in detectable testes damage. The suitability of using these 2 mutants as a vaccine was further evaluated by infecting additional mice with either CpG_max or UpA_max ZIKV mutants, with the size of their testes determined 49 days postinfection. No reduction in testes size was observed when compared to noninfected (naïve) mice ( Fig 7D). RealAU : PleasenotethatqRT À PCRhasbeendefinedasreal À timequantitativerev -time quantitative reverse transcription PCR (qRT-PCR) of the testes also illustrated that neither CpG_max nor UpA_max ZIKV mutants contained detectable levels of ZIKV RNA (Fig 7E). Together, these results indicate that the CpG_max and UpA_max ZIKV mutants do not cause testes damage and can protect against ZIKV challenge.

Discussion
The motivation for this study was the evident difference in CpG dinucleotide representation in the genomes of RNA viruses infecting vertebrates compared to those infecting insects (Fig 1) [5,13,14,31]. A recent study reported that ZIKV mutants with elevated CpG dinucleotide frequencies were attenuated in vertebrate cells and mice, while no significant effects on virus replication were observed in C6/36 Aedes albopictus cells [36]. This suggests that CpG suppression in arbo flaviviruses is required for sufficient replication in vertebrate hosts and that there is no selection on CpG dinucleotide frequencies in mosquitoes or other arthropod

PLOS BIOLOGY
hosts. However, in our experiments, ZIKV mutants with elevated CpG dinucleotide frequencies (CpG representations of 1.0 and 1.62) reproducibly showed enhanced replication rates in a range of different Aedes spp. mosquito cell lines (Fig 3). Importantly, live vector mosquitoes became more frequently infected with the CpG-high ZIKV mutants and dissemination into the saliva occurred more often compared to infections with WT ZIKV (Fig 4). In marked contrast to these proviral effects of increased CpG frequencies, elevated UpA dinucleotide frequencies resulted in reduced replication rates in mosquito cell cultures and lowered virus titers in infected Ae. aegypti mosquitoes (Figs 3 and 4). Taken together, this study provides the first

PLOS BIOLOGY
experimental evidence that dinucleotide composition of viral genomes is subject to substantial evolutionary pressures in invertebrate mosquito cells.
This new information sheds light on the biological basis for the compositional differences between flaviviruses that are arthropod-borne and those that replicate exclusively in arthropod hosts (ISFs). Both groups of flaviviruses persistently infect their arthropod hosts without any noticeable pathogenicity. The observation that elevated CpG frequencies enhance ZIKV replication in mosquito cells (Figs 3 and 4) suggests that WT ZIKV infections in mosquitoes are normally attenuated. The substantial mortality in mosquitoes infected with the CpG_max ZIKV ( Fig 4D) further illustrates how finely balanced CpG composition has to be maintained for flaviviruses to persistently infect their arthropod hosts, especially given that the mutated region represents only 15% of the ZIKV genome. Lineage 2 ISFs mildly suppress CpG dinucleotides (Fig 1) and are believed to have evolved from arboviruses that lost the ability to infect vertebrate cells [31]. Rather than drifting toward an unbiased CpG dinucleotide frequency due to a lack of selection pressure, the observed evolutionary pressures in mosquito cells suggest that lineage 2 ISFs accumulate CpG dinucleotides by natural selection. However, the accumulation of CpG dinucleotides in lineage 2 ISF genomes is potentially constrained as nonpathogenic persistent replication needs to be maintained.
This delicate balance additionally has to be sustained in the context of antiviral responses. In the arthropod host, the most potent antiviral responses are mediated through RNAi (reviewed in [37]). Suppression of RNAi in mosquitoes infected with Sindbis virus (arbovirus from family Togaviridae, genus Alphavirus) results in increased mosquito mortality [38]. However, replication of the ZIKV mutants with elevated CpG or UpA dinucleotide frequencies was similarly enhanced or suppressed, respectively, in all mosquito cells tested, including C6/ 36 cells (Fig 3), which are defective in the production of antiviral siRNAs [39,40]. This suggests that other pathways or mechanisms are active in discriminating dinucleotide frequencies in mosquitoes.
Alternatively, synonymous changes in dinucleotide frequencies inevitably alter the codon usage of a nucleotide sequence and changes in translation efficiency may affect the biosynthesis rate of the encoded viral polyprotein and as a consequence viral replication rates. However, comparison of the synonymous codon usage of the mutated region in ZIKV and the synonymous codon frequency of human and insect reference sets show no consistent correlation between the elevation of CpGs and similarities in codon usage with either host reference set (S1 Table). Indeed, previous studies have largely discounted effects of codon use or codon pair bias as the explanation for attenuation of CpG-high mutants of vertebrate viruses. For example, elevation of CpG and UpA dinucleotide frequencies in untranslated regions of viral RNA attenuated virus replication to equivalent degrees as similar dinucleotide modifications in coding regions of the viral genome [15,16,41]. Furthermore, it is now established that the vertebrate antiviral protein ZAP mediates the attenuated replication of viruses with high CpG frequencies.
In vertebrates, the N-terminal zinc-fingers of ZAP bind high CpG containing RNA sequences through a number of hydrogen bonds and stacking interactions [22] to initiate a plethora of downstream intracellular interactions with proteins involved in translation, RNA degradation, and antiviral immunity [20,[23][24][25][26][27][28]42]. Specifically, the host proteins RNaseL, activated by OAS3, and both TRIM25 and KHNYN are crucial for the ZAP-mediated attenuation of CpG-high E7 and HIV-1, respectively [16,21,29]. Correspondingly, knockout of ZAP or OAS3/RNaseL abolished the attenuated phenotype of ZIKV mutants with elevated CpG dinucleotide frequencies in vertebrate cells (Fig 2). Additionally, mice vaccinated with CpG_max ZIKV and subsequently challenged with ZIKV PRVABC59 did not show any signs of ZIKV-induced testicular damage as measured by their testis size, whereas vaccination with CpG_1.0 ZIKV resulted in an intermediate reduction in average testis size (Figs 6 and 7). Together, these results indicate that the degree of virus attenuation is correlated to CpG composition and mediated by vertebrate ZAP.
While most recent investigations involving genome composition and virus replication have focussed on CpG dinucleotides, a growing body of evidence indicates a substantial additional effect of UpA dinucleotide frequencies on virus attenuation. UpA dinucleotides are suppressed in vertebrate and invertebrate animals and their viruses (Fig 1) [14]. The UpA_max ZIKV mutant used in this study showed strong attenuation in vertebrate cells and no signs of disease in mice (Figs 2 and 6). Moreover, the attenuated phenotype of UpA_max ZIKV was rescued by knockout of ZAP, suggesting that the interaction between ZAP and UpA dinucleotides is not exclusive for previously published E7 mutant viruses [21]. Whether ZAP directly binds UpA dinucleotides or requires additional host proteins for this interaction is unknown and remains a topic of active investigation [43]. Whether the selection against UpA dinucleotides in insect cells is equally strong as that observed in vertebrate cells remains unclear. Although flaviviruses suppress UpA dinucleotides in their RNA, UpA_max ZIKV replication and dissemination was not consistently reduced compared to WT ZIKV in mosquito cells and live mosquitoes, respectively (Figs 3-5).
The diverse selective pressures that arboviruses encounter have influenced arbovirus disease emergence on multiple levels (reviewed in [44]). The current study identifies CpG frequencies as a key adaptive trait that enables ZIKV to replicate in 2 profoundly different hosts. Furthermore, the finding that increasing CpG frequencies enhanced the replication of ZIKV in mosquitoes provides evidence that CpG frequency differences between vector-borne flaviviruses and ISFs do not simply result from compositional constraints in one host and none in the other. Instead, the data indicate the existence of strongly conflicting evolutionary pressures on ZIKV to enable sufficient replication in both hosts to maintain a vector-borne life cycle. In contrast to flaviviruses, vector-borne viruses from the genus Alphavirus (family Togaviridae) contain relatively high CpG dinucleotide frequencies, which may suggest these viruses have evolved a more optimal genome composition for replication in the mosquito vector. Although most alphaviruses effectively infect vertebrate cells, these viruses are highly sensitive to ZAP [45]. The same constraints may operate in other dual-host viruses. For example, we and others recently discovered that the replication of potato virus Y and plum pox virus (genus Potyvirus) in plants was similarly attenuated by the synonymous introduction of CpG and UpA dinucleotides [46,47]. Plant transcriptomes also display CpG suppression that is also reflected in the majority of plant RNA viruses. This suggests that plant viruses which are maintained in transmission cycles involving persistent replication in both plant host and arthropod vectors [48] may encounter similarly conflicting evolutionary pressures as arboviruses.
The findings can spark investigations into novel strategies to intervene with the transmission cycles of arboviruses. Attenuating vertebrate viruses by increasing CpG and UpA dinucleotide frequencies is a promising and universally applicable approach to generate liveattenuated virus vaccines (Figs 6 and 7) [17,36]. The enhanced replication of CpG-high mutants of ZIKV in mosquito cells may serve to improve production of live-attenuated or inactivated vaccines and can lead to novel immunization strategies [49]. A radical, yet intriguing application is to release CpG-high infected mosquitoes into the environment to serve as vaccine vectors. At the moment, this application remains purely hypothetical as such an approach raises many questions regarding ethics, vaccine safety, and the spread of genetically modified viruses into the environment. However, theoretically mosquitoes could be highly effective (Fig 5) in not only vaccinating humans and livestock but also animal reservoirs by transmitting high doses of flavivirus vaccines with the same attenuated replication in vertebrate animals as observed in the experimentally infected mice (Figs 6 and 7).

In silico analyses
A single full-length flavivirus genome was selected per flavivirus species. For analysis of host gene composition, mRNA sequences of human (GRCh38.p13i) and Ae. aegypti (AaegL5.0) were downloaded from RefSeq (https://www.ncbi.nlm.nih.gov/nuccore). From these datasets, representative selections of 10,000 sequences were extracted at random using Perl with an inhouse script. Sequences were filtered to exclude coding region sequences of length <250 bases. To avoid redundancy, similar or identical sequences (<1% nucleotide divergence) were identified with SSE software (v1.4) [50] and removed from the dataset. This resulted in over 9,000 human and mosquito sequences, which were used for composition analysis (all GenBank accession numbers available in S1 Text). Dinucleotide frequencies were calculated using the program Composition Scan within SSE (v.1.4). Synonymous site variability (SSV) and MFED of the ZIKV genome were calculated as previously described, using a fragment size of 200 and sliding window of 21 codons for SSV and a fragment size 200 with 48 nucleotide increments for MFED using a selection of ZIKV isolates [51]. The codon adaptation index (CAI) of full length ZIKV and the mutated regions was calculated using an online CAI calculator (www. biologicscorp.com) with the codon usage of H. sapiens and Drosophila melanogaster as reference data sets (S1 Table).

Cells and viruses
Ae. albopictus C6/36 (ATCC CRL-1660) and U4.4 were cultured at 28˚C in Leibovitz-L15 medium (Gibco, Carlsbad, California, United States of America) supplemented with 10% heat-inactivated fetal bovine serum (FBS; Gibco), 2% tryptose phosphate broth (Gibco), and 1% nonessential amino acids (Gibco). Ae. albopictus RML-12 cells were cultured in Leibovitz's L-15 medium supplemented with 10% FBS and 10% tryptose phosphate buffer at 28˚C [52]. Ae. aegypti Aag2 cells [53] and Ae. pseudoscutellaris AP-61 cells [54] were cultured at 28˚C in Schneider's Drosophila medium (Gibco) supplemented with 10% FBS. Human A549 cells including ZAP, RNaseL, 2 0 -5 0 -oligoadenylate synthetase 1 (OAS1), and OAS3 knockout cells were described previously [21], and African green monkey kidney Vero E6 (ATCC CRL-1586) were all maintained at 37˚C with 5% CO 2 in Dulbecco's modified Eagle medium (DMEM) (Gibco) supplemented with 10% heat-inactivated FBS, 100 U/ml penicillin, and 100 μg/ml streptomycin (SigmaAU : PleaseprovidemanufacturerlocationdetailsforSigma À AldrichinthesentenceH -Aldrich). Cell lines were routinely checked for mycoplasma using MycoAlert Mycoplasma Detection Kit (LonzaAU : PleaseprovidemanufacturerlocationdetailsforLonzai ). A pCC1-BAC-based and SP6 promoter-driven infectious clone of the ZIKV BeH81901 isolate from Brazil (GenBank KU365778.1) [55] was used for the construction of mutant viruses. A region of 1,653 nucleotides was selected for the described synonymous mutations. Mutant regions with elevated CpG or UpA dinucleotide frequencies were designed using the "mutate sequence" function in SSE, taking care not to alter the coding amino acid sequence or the other UpA or CpG dinucleotide, respectively. The SCR sequence was originally designed for a region of 4,344 nucleotides that included the above mentioned mutated region using the "scramble sequence" function in SSE with the CDLR algorithm. The amino acid sequence, codon usage, and dinucleotide usage of the SCR sequence were kept identical to that of the corresponding WT sequence. Mutated regions were synthesized (Geneart) and incorporated into the viral genome using unique restriction sites RsrII and Bst BI (NewAU : Pleaseprovidemanufacturerlo England Biolabs). Nucleotide sequences of mutated regions are listed in S2 Text. Plasmid DNA was linearized with PmeI (New England Biolabs) and complete digestion was verified on agarose gel. Linearized DNA was purified by phenol/chloroform extraction and RNA was in vitro transcribed in the presence of G(5 0 )ppp(5 0 )G RNA Cap structure analog using SP6 polymerase (New England Biolabs) according to the manufacturer's protocol. RNA was transfected into Vero cells with lipofectamine 2000 (InvitrogenAU : PleaseprovidemanufacturerlocationdetailsforInvitrogeninthesenten ) to launch virus replication. Between 1 and 4 days posttransfection, cytopathic effects were observed and cell culture supernatant was transferred to a fresh 60% confluent monolayer of Vero cells to generate passage 1 (P1). After 2 hours, inoculum was removed and cells were washed with cell culture medium. Three days postinfection, cell culture supernatant was harvested. Newly generated mutant virus stocks were sequence verified by reverse transcriptase PCR (see competition assay below for details) followed by restriction analysis using DdeI, SalI, and PstI restriction enzymes (New England Biolabs) and DNA sequencing (Macrogen EuropeAU : PleaseprovidemanufacturerlocationdetailsforMacrogenEuropeinthe ). Stocks were again verified by restriction digestion prior to the mouse experiments. At the end of the mosquito experiments, total RNA from pools of 8 infected mosquitoes was isolated and ZIKV mutants were sequence verified by RNA-sequencing (BGIAU : PleaseprovidemanufacturerlocationdetailsforBGIGenomicsinthesentenceAttheendofthemo Genomics) as described previously [56].

Virus quantification
Tissue culture infectious dose 50% (TCID 50 ) was determined by end-point dilution assays (EPDA). Serial 10-fold dilutions of virus stocks and experimental samples were prepared in complete DMEM cell culture medium. Vero E6 cells from confluent tissue cultures were detached with Trypsin-EDTA (Gibco) and dissolved in a total volume equal to 1 ml/1.67 cm 2 .
Cell suspension was added to the virus dilutions in a 1:1 ratio, mixed and plated in 6-fold onto microtiter plates (NuncAU : PleaseprovidemanufacturerlocationdetailsforNuncinthesentenceCellsuspen ). Virus titers were determined by CPE at the end-point of infection. For each of the mutant viruses, the viral titers determined by CPE and respective end-points of the EPDA were confirmed by immunofluorescence. Cells were washed with phosphate buffered saline (PBS) and fixed with 4% paraformaldehyde in PBS. ZIKV envelope was stained with 4G2, a pan-flavivirus primary α-E mouse monoclonal antibody diluted in PBS containing 2% FBS (1:50) [57] for 1 hour at 37˚C. Samples were washed 3 times with PBS and stained with a secondary goat α-mouse Alexa Fluor 488 (1:2,000; Invitrogen). RNA copy numbers in virus stocks were determined by comparing RNA from virus stocks with purified SP6 in vitro transcribed RNA. RNA isolation was performed with TRIzol reagent for liquid samples according to the manufacturer's protocol. Reverse transcriptase and a real-time PCR were performed with 2 sets of ZIKV-specific primers F1 5 0 -cctgatgaccatctgtggca-3 0 , R1 5 0 -cccatagagcaccactcctt-3 0 , and F2 5 0 -tcaggaggtggtgttgaagg-3 0 , R2 5 0 -gtgccctttctccatttggt-3 0 using superscript III (Invitrogen) and SYBR green (AppliedAU : PleaseprovidemanufacturerlocationdetailsforAppliedBio Biosystems) according to the respective manufacturer's protocols with an annealing temperature of 58˚C during real-time PCR in a CFX96 realtime system (Bio-Rad) thermocycler.

Virus growth curves
Cell monolayers were seeded in 6-well plates to a 50% to 60% confluency. The cell culture fluid was removed and infections were performed by diluting virus stocks (P2) in standard culture media for each of the respective cell types in a total volume of 1 ml/well. After 2 hours, the inoculum was removed and the monolayers were washed twice with 1 ml of PBS, before 2 ml of fresh culture medium was added. At indicated times postinfection, 50 μl of cell culture supernatant was taken per well and stored at −80˚C until TCID 50 /ml was determined by EPDA.

Competition assay
Equal concentrations of WT and mutant virus (1 RNA/cell) were mixed in DMEM or XXX cell culture medium and added to A549, A549 ZAP knockout, or AP-61 mosquito cells at 50% confluence. After 2 hours, the inoculum was replaced with culture medium. Ten μl of cell culture supernatant was passaged to fresh cells. After 2 passages, total RNA was isolated using TRIzol reagent according to the manufacturer's protocol. RNA was reverse transcribed and part of the mutated region within the ZIKV genome was PCR amplified using SuperScript III One-Step RT-PCR System with Platinum Taq Polymerase (Invitrogen) together with either the 5 0 -aagtggagaaaaagatggg-3 0 (WT and SCR) or 5 0 -gaggaaaaagagtggaagac-3 0 (WT, CpG_1.0, CpG_max, and UpA_max) forward primers in combination with a 5 0 -gtgccctttctccatttggt-3 0 reverse primer. Reverse transcription and PCR amplification were performed according to the manufacturer's protocol with the reverse transcription at 50˚C for 30 minutes and an annealing temperature during PCR amplification of 58˚C in an Applied Biosystems 2720 thermocycler. The resulting DNA amplicons from the comparisons between WT and SRC were digested with BsaI-HFv2 (New England Biolabs), WT and CpG_1.0 or CpG_max with SalI-HF (New England Biolabs) and WT and UpA_max with HaeIII (Invitrogen) and ran on 2% agarose gels.

Infectious blood meal
Virus solutions were made by diluting the virus stocks to equal concentrations in DMEM before mixing with human blood to their final concentrations. Infectious blood meals were offered through Parafilm using the Hemotek PS5 feeder. Mosquitoes from a single bucket were offered a single blood meal containing one of the ZIKV mutants and allowed to feed for 1 hour ad libitum in light conditions, at 24˚C and 70% RH. Mosquitoes were anesthetized with 100% CO 2 and fully engorged females were selected. Engorged mosquitoes were kept for 14 days on glucose solution at 28˚C in a 12:12 light:dark cycle.

Salivation assay
Fourteen days post blood meal, mosquitoes were anesthetized with 100% CO 2 and placed on a CO 2 pad. Mosquitoes that died within the 14 days incubation period were counted and discarded. Mosquitoes were immobilized by removing their legs and wings. The proboscis of each mosquito was inserted into a 200-μl yellow pipet tip (GreinerAU : Pleaseprovidemanufacturerloca Bio-One) containing 5 μl of a 1:1 solution of 50% glucose solution and FBS, for 1 hour. After salivation, the mosquito bodies were added to 1.5 ml Safe-Seal micro tubes (Sarstedt, Nümbrecht, Germany) containing 0.5 mm zirconium beads (Next Advance, Averill Park, New York, USA). Each saliva sample was added to a 1.5-ml micro tube (Sarstedt) with 55 μl DMEM cell culture medium containing FBS, Penicillin and streptomycin supplemented with additional Fungizone (50 μg/ml; Invitrogen), and Gentamycin (50 μg/ml; LifeAU : PleaseprovidemanufacturerlocationdetailsforLifeTechnologi Technologies), hereafter referred to as DMEM-complete. Mosquito bodies and saliva samples were stored at −80˚C until further processing.

Infectivity assays
Frozen mosquito bodies were homogenized for 2 minutes at speed 10 in a Bullet Blender Storm (Next AdvanceAU : PleaseprovidemanufacturerlocationdetailsforNextAdvanceinthesentenceFroz ). Homogenized bodies were centrifuged briefly and resuspended in 100 μl DMEM-complete medium. Homogenization was repeated and mixtures were centrifuged for 1 minute at 14.500 rpm in a table top centrifuge. From the mosquito homogenates or thawed saliva samples, 30 μl was used to inoculate a Vero cell monolayer in a 96-well plate. After 2 to 3 hours, the inoculum was removed and replaced with 100 μl DMEM-complete. Wells were scored for virus induced CPE at 3 and 4 dpi. Bodies and saliva samples of infected mosquitoes were additionally titrated by EPDA.

Immunostained plaque assays
An 80% confluent monolayers of cells were infected with a series of 10-fold dilutions with equal RNA copies of each of the mutant viruses. After 2 hours of incubation, the inoculum was removed and a 1% carboxymethyl cellulose overlay was added. At the end of the experiment, the overlay was removed and cells were fixed and stained for ZIKV envelope protein using 4G2 primary antibodies and green Alexa Fluor 488 secondary fluorophores. Micrographs were taken with an Axio Observer Z1m inverted microscope in combination with an X-Cite 120 series lamp and AxioVision software of all primary plaques at the lowest possible concentration.

Infection of IFNAR −/− male mice with mutant ZIKVs
All viruses were grown in C6/36 cells and titers were determined by immunostained TCID 50 assays using C6/36 cells. Briefly, virus was titrated in 10-fold serial dilutions in 96-well plates with C6/36 cells (2 × 10 4 cell/well). After 5 days, plates were washed and fixed in 80% cold acetone for 1 hour at −20˚C. Acetone was removed. Plates dried and ZIKVs were detected by indirect immunofluorescence using mouse anti-4G2 monoclonal primary antibody and polyclonal goat anti-mouse IgG conjugated IRDye 800CW secondary antibody (LiCor, 926-32210). Staining was visualized using the LI-COR Biosciences Odyssey Infrared Imaging System and application software. Any staining in any well was considered positive and TCID 50 titers determined.
Male type I IFN receptors knockout (IFNAR −/− ) mice (C57BL/6J background) (n = 5 per group, age 2 to 6 months) were infected subcutaneously (s.c.) at the base of the tail with 10 3 TCID 50 (100 ul in Roswell Park Memorial Institute (RPMI) 1640 supplemented with 2% FCS) of each virus per mouse. Mice were weighed at the indicated times and viremias determined using immunostained TCID 50 assays as above. Briefly, mouse serum was titrated in 10-fold serial dilutions into 96-well plates with C6/36 cells and ZIKVs were detected by indirect immunofluorescence as above, except background IgG was blocked prior to addition of the secondary antibody with Goat F(ab) anti-mouse IgG (H&L) (AbcamAU : Pleaseprovidemanufacturerlocationd ).

Antigen-coated ELISA
Anti-ZIKV IgG titers in mouse sera were determined by antigen-coated ELISA. Plates were coated with MR766 antigen in BupH Carbonate-Bicarbonate buffer (ThermoFisherAU : Pleaseprovidem ) and incubated O/N at 4˚C. After incubation, plates were blocked with 5% skim milk in PBS. Next, the mouse sera were added in triplicates diluted 1:30 in PBS-Tween containing 1% skim milk and incubated for 1 hour at room temperature. After incubation, plates were washed with PBS-Tween. Subsequently, secondary antibody HRP (1:2,000 in PBS-Tween with 1% skim milk) was added and incubated for 1 hour at room temperature. Samples were washed with PBS-Tween and scored after 1-hour incubation in the dark with ABTS substrate (1:1,000 H 2 O 2 ).

ZIKV PRVABC59 challenge
An infectious clone ZIKV PRVABC59 (GenBank: MH158237.1) was constructed as described previously [58] and was provided by S. Tajima (Department of Virology I, National Institute of Infectious Diseases, Tokyo, Japan). Infectious virus was recovered by transfection of Vero E6 cells. ZIKV PRVABC59 stocks were grown in subconfluent C6/36 mosquito cells maintained in RPMI1640 supplemented with 2% FBS and stored in −80˚C. Viral stock titers (TCID 50 ) were determined by 10-fold serial titration in C6/36 cells (2 × 10 4 cells/well of a 96-well plate) and incubated for 5 days. After incubation, 25 μl from each well was placed into parallel plates containing Vero E6 (2 × 10 4 cells/well), and the cytopathic effect assessed on day 7.
On day 32 postprimary injections, mice were challenged s.c. base of tail with ZIKV PRVABC59 (10 3 TCID 50 in 100 μl) as described previously [49]. Five mice were mock injected with 100 μl PBS and used as uninfected controls. Mice were weighed at the indicated times. Mouse serum was collected by tail vein bleed and ZIKV PRVABC59 titers determined by TCID 50 assays as described above. Mice were euthanized on day 80 postprimary infection, testes were removed, photographed, and measured (height × width).

Vaccine safety
To illustrate that the CpG_max and UpA_max ZIKV mutants do not cause testes damage and are suitable for use as attenuated vaccines, 2 groups of male IFNAR −/− mice (n = 6 mice per group) were infected with these ZIKV mutants (10 3 TCID 50 ) and testes harvested 49 days postinfection and compared to testes from uninfected mice. qRT-PCR for ZIKV RNA was also performed on the testes as described previously [49].

Statistics
Statistical analyses were done with the software package Prism (v. 9.0.0) and Microsoft Excel 2016 and detailed in the main text and in the relevant figure legend. All two-way analyses of variance (ANOVA) were performed to compare replication of ZIKV variants in multiple cell lines. Outcomes were corrected for multiple comparisons by controlling the false discovery rate (FDR) (Q = 0.05). Reported p-values for two-way ANOVAs are FDR adjusted p-values unless stated otherwise.

Ethics statement for murine experiments
All mouse work was conducted in accordance with the "Australian Code for the care and use of animals for scientific purposes" as defined by the National Health and Medical Research Council of Australia. Mouse experiments and associated statistical treatments were reviewed and approved by the QIMR Berghofer Medical Research Institute animal ethics committee (P2195). Representative images of ZIKV induced immunoplaques in A549 and A549 ZAP knockout cells. Cells were fixed and stained for ZIKV E protein with primary antibody 4G2 and secondary Alexa Fluor 488 (green). Results of an immunoplaque assay on Vero E6 cells (E). Data points represent the average plaque size, error bars indicate 1 standard error of the mean, and asterisks highlight significant differences from WT (p < 0.05, one-way ANOVA with Dunnett post hoc test). Please refer to S1 Data for the numerical values underlying panels A-C and E. dpiAU : AbbreviationlistshavebeencompiledforthoseusedinS1 À S4Figs:Pleaseverifythatallentriesareco , days postinfection; OAS1, 2 0 -5 0 -oligoadenylate synthetase 1; OAS3, 2 0 -5 0 -oligoadenylate synthetase 3; RNaseL, ribonuclease L; SCR, scrambled control virus; WT, wild-type; ZIKV, Zika virus.  Table. Composition of mutated regions. CpG and UpA dinucleotide composition of the mutated region in ZIKV and the codon adaptation index compared to both Homo sapiens and Drosophila melanogaster. AdditionalAU : PleaseconfirmwhethertheedittothesentenceAdditionalnumerica numerical values that underlie the values presented in this table are available in S1 Data. (PDF) S1 Text. Accession numbers used in this study. Accession numbers that refer to RNA sequences from Homo sapiens, Aedes aegypti, and the flavivirus genomes that were used to generate Fig 1A and 1C.  values that have made up Fig 1A, 1C, 2A, 2B, 2C, 2D, 2E, 3A, 3B, 3C, 5B, 5C, 5D, 5E, 6B, 6C,  6D, 7A, 7B, 7C, 7D Table. (XLSX)