Infection, dissemination, and transmission efficiencies of Zika virus in Aedes aegypti after serial passage in mosquito or mammalian cell lines or alternating passage in both cell types
Parasites & Vectors volume 14, Article number: 261 (2021)
Zika virus (ZIKV) is an arthropod-borne virus (arbovirus) with an urban transmission cycle that primarily involves humans and Aedes aegypti. Evidence suggests that the evolution of some arboviruses is constrained by their dependency on alternating between disparate (vertebrate and invertebrate) hosts. The goals of this study are to compare the genetic changes that occur in ZIKV after serial passaging in mosquito or vertebrate cell lines or alternate passaging in both cell types and to compare the replication, dissemination, and transmission efficiencies of the cell culture-derived viruses in Ae. aegypti.
An isolate of ZIKV originally acquired from a febrile patient in Yucatan, Mexico, was serially passaged six times in African green monkey kidney (Vero) cells or Aedes albopictus (C6/36) cells or both cell types by alternating passage. A colony of Ae. aegypti from Yucatan was established, and mosquitoes were challenged with the cell-adapted viruses. Midguts, Malpighian tubules, ovaries, salivary glands, wings/legs and saliva were collected at various times after challenge and tested for evidence of virus infection.
Genome sequencing revealed the presence of two non-synonymous substitutions in the premembrane and NS1 regions of the mosquito cell-adapted virus and two non-synonymous substitutions in the capsid and NS2A regions of both the vertebrate cell-adapted and alternate-passaged viruses. Additional genetic changes were identified by intrahost variant frequency analysis. Virus maintained by continuous C6/36 cell passage was significantly more infectious in Ae. aegypti than viruses maintained by alternating passage and consecutive Vero cell passage.
Mosquito cell-adapted ZIKV displayed greater in vivo fitness in Ae. aegypti compared to the other viruses, indicating that obligate cycling between disparate hosts carries a fitness cost. These data increase our understanding of the factors that drive ZIKV adaptation and evolution and underscore the important need to consider the in vivo passage histories of flaviviruses to be evaluated in vector competence studies.
Zika virus (ZIKV) is an arthropod-borne virus (arbovirus) that belongs to the genus Flavivirus (family Flaviviridae). ZIKV was first isolated from a febrile rhesus monkey in Zika Forest, Uganda, in 1947, with subsequent isolations made from Aedes africanus in the same region the following year . The first known infections in humans occurred in Nigeria in 1954 [2, 3]. Afterwards, the virus was detected elsewhere in Africa and Asia, but cases of human disease were sporadic . ZIKV emerged in the Pacific Islands in 2007, causing an explosive outbreak of febrile illness on Yap Island, with infections occurring in most inhabitants . Additional outbreaks of ZIKV have occurred elsewhere in the Pacific Islands [6, 7]. The virus underwent a dramatic expansion of its geographic distribution in 2015, when it was reported for the first time in the Western Hemisphere . ZIKV now occurs throughout the tropical Americas.
ZIKV is primarily maintained in transmission cycles involving Aedes spp. mosquitoes and human or nonhuman primates, although other modes of transmission (i.e. sexual transmission between humans) have been documented. The principal urban transmission vector of ZIKV in the Americas is Aedes aegypti, but other Aedes spp. mosquitoes also transmit the virus [8,9,10]. In humans, ZIKV infections are usually asymptomatic or self-limiting but serious manifestations, such as congenital microcephaly and Guillain-Barré syndrome in adults, can occur [11,12,13].
ZIKV has a positive-sense, single-stranded RNA genome. The genomes of ZIKV and most other RNA viruses encode an RNA-dependent RNA polymerase, an enzyme that incorporates a relatively high number of mutations, ranging from 10–4 to 10–6 mutations per site per round of replication, due to the lack of a proofreading exonuclease activity [14, 15]. Most mutations have no effect on viral fitness or are quickly removed from the virus population by natural selection. Some mutations, however, are advantageous and can result in enhanced virus transmission or increased disease severity, resulting in epidemics of increased severity, as documented for several arboviruses, including West Nile virus (WNV), Venezuelan equine encephalitis virus (VEEV), chikungunya virus (CHIKV) and ZIKV [16,17,18,19,20]. An amino acid substitution (V159A) in domain I of the envelope protein of WNV appeared in 2002 and resulted in a significantly decreased extrinsic incubation period in Culex mosquitoes, thereby enhancing vectorial capacity . Additionally, an amino acid substitution (T249P) in the NS3 protein of a relatively attenuated strain of WNV produced a virus that was highly lethal to American crows . A mutation (S218N) in the E2 protein of VEEV enhanced viral infection, replication and transmission in Aedes taeniorhynchus. The mutation was linked to an epizootic outbreak of VEEV in horses from Mexico in 1993–1996 . In 2007, a substitution (A226V) in the E1 protein of CHIKV was responsible for significantly increased infectivity and more efficient dissemination and transmission of the virus by Aedes albopictus, facilitating its rapid spread . Genetic changes in the viral RNA genome can therefore have a major effect on arbovirus transmission, epidemiology and disease presentation.
RNA viruses that cycle between disparate (vertebrate and invertebrate) hosts typically experience lower evolutionary rates than RNA viruses that replicate in a single host [21, 22]. Alternate host replication has been hypothesized to impose a fitness trade-off that constrains genome evolution and host-specific adaptation. Many studies have been performed to test this hypothesis and to evaluate the genotypic and phenotypic consequences of sequential and alternating host cell passaging of arboviruses [23,24,25,26,27,28,29,30,31,32,33]. One of the earliest studies was performed by Novella and colleagues, who reported a lack of evolutionary stasis during alternating replication of vesicular stomatitis virus in baby hamster kidney (BHK-21) and sandfly (LL-5) cell lines . In another study, VEEV serially passaged in Ae. aegypti exhibited increased mosquito infectivity when compared to parental virus while sequential passage in mice yielded a virus that produce higher viremias . Alternately passaged VEEV exhibited no apparent fitness gains in either host. The goal of this study was to determine the effects of in vitro host specialization on the genomic composition of ZIKV and in vivo replicative ability of the virus in geographically matched mosquitoes.
Cell lines and virus
African Green Monkey kidney (Vero) and Aedes albopictus (C6/36) cells were obtained from the American Type Culture Collection (Manassas, VA). Vero cells were cultured in Dulbecco’s modified Eagle medium (Invitrogen, Carlsbad, CA), and C6/36 cells were cultured in Leibovitz L15 medium (Invitrogen). All media were supplemented with 10% fetal bovine serum (FBS), 100 units/ml penicillin and 100 μg/ml streptomycin. Mammalian cells were cultured at 37 °C with 5% CO2 whereas C6/36 cells were cultured at 28 °C. ZIKV (strain ZIKV-Mer17) was originally isolated from a febrile patient in Merida, Yucatan, in 2017 and had been passaged four times in C6/36 cells prior to this study. The virus was titrated by plaque assay following published protocols , and titers were expressed as plaque-forming units per milliliter (PFU/ml).
Serial passaging of ZIKV in cell culture
ZIKV was sequentially passaged six times in one cell type (Vero or C6/36 cells) or alternately between both cell types (Vero then C6/36 cells) in triplicate 75 cm2 flasks. Original cultures were inoculated with ZIKV at a multiplicity of infection (m.o.i) of 0.01, and subsequent cultures were inoculated with 100 µl of cell culture supernatant collected from the prior passage. Supernatants were harvested from Vero cells when 70–90% of the cell monolayer exhibited cytopathic effect and from C6/36 cells at 6 days post-infection. Aliquots of the triplicate supernatants harvested after the final passage were consolidated before use in the vector competence experiments.
Immature Ae. aegypti were collected at field sites in 2018 in the city of Merida, Yucatan. Mosquitoes were transported to the insectaries at the Laboratorio de Arbovirología at the Universidad Autónoma de Yucatán and reared to adults. After mating, eggs were collected, dehydrated and sent to Colorado State University (CSU). Eggs were hatched and reared to adults in the insectaries at the CSU Arthropod-Borne and Infectious Diseases Laboratory. Adults were maintained with water and sucrose ad libitum under standard conditions (28 °C, 70–80% humidity, light-dark photoperiod of 12:12 h). To facilitate egg production, females were provided defibrinated calf blood at 37 °C using an artificial membrane and water-jacketed glass feeder. All experimental infections were performed using F3 adults.
Infectious blood meals were prepared by mixing virus (1 × 106 PFU/ml) and defibrinated calf blood at a ratio of 1:1. Blood meals were offered for 30 min to 5–9 day-old Ae. aegypti using a water-jacketed feeder and hog gut membrane. The mosquitoes were deprived of sucrose and water for 24 h prior to viral infection. Mosquitoes were cold-anesthetized, and engorged females were selected and held in the insectary. At 3, 5, 7, 14 and 21 days post-infection (dpi), 27–30 females were selected from each group and cold-anesthetized. Wings and legs were removed and placed into 2-ml tubes containing 200 μl of diluent and one stainless steel bead. The diluent consisted of phosphate buffer saline (PBS), pH 7.4, supplemented with 20% FBS, antibiotics and fungizone. Once the wings and legs were removed, saliva was obtained from each female as previously described by Smith et al. in 2005 . Briefly, the proboscis was placed into a capillary tube containing immersion oil, allowing females to expectorate in the tube. After 30 min, the capillary tube was broken into a conical tube containing 100 μl of diluent. Once salivation was complete, midguts, Malpighian tubules, ovaries, and salivary glands were harvested. Tissues were rinsed individually in a drop of PBS and placed in separate tubes containing 200 μl of diluent and a stainless steel bead. Tissue homogenization was carried out at 24 Hz for 5 min. Tubes containing homogenates and saliva were centrifuged at 20,000× g for 5 min and stored at − 80 °C. These experiments were performed in duplicate in the Biosafety Level 3 facilities at CSU.
Midgut infection rates (MIRs), dissemination infection rates (DIRs), salivary gland infection rates (SGIRs), and transmission rates (TRs) were calculated and compared between viruses. MIR was defined as the number of mosquitoes with infectious ZIKV in their midguts as determined by plaque assay divided by the total number of mosquitoes examined multiplied by 100. DIR was defined as the number of mosquitoes with ZIKV in non-midgut tissues (e.g. wings/legs, salivary glands, Malpighian tubules or ovaries) divided by the number of mosquitoes with ZIKV in their midguts multiplied by 100. SGIR was defined as the number of mosquitoes with ZIKV in their salivary glands divided by the number of mosquitoes with ZIKV-infected midguts multiplied by 100. TR was defined as the number of mosquitoes with ZIKV in their saliva divided by the total number of mosquitoes with disseminated infections multiplied by 100. DIRs were often not determined if no more than five of the mosquitoes in the group had midgut infections. Likewise, SGIRs and TRs were often not determined if no more than five of the mosquitoes in the group contained virus in non-midgut tissues and their salivary glands, respectively. Due to the large number of tissues tested and, more importantly, the small size of select tissues (particularly the salivary glands), it was not feasible to perform duplicate serial dilutions of each tissue homogenate for the plaque assay analysis in order to calculate virus titers. Instead, plaque assay data were recorded as positive or negative.
Total RNA was extracted from cell culture supernatants using the Mag-Bind® Viral DNA/RNA 96 kit (Omega Bio-Tek) on the KingFisher Flex Magnetic Particle processor (Thermo Fisher Scientific, Waltham WA) and assayed for ZIKV RNA by quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR). Full-length ZIKV RNA generated from a plasmid using T7 polymerase and diluted tenfold (108 to 102 ng/µl) to create a standard curve, which was used to calculate viral RNA copy numbers in select mosquito organs. Reactions were performed using SsoAdvanced Universal SYBR Green supermix (Bio-Rad Laboratories, Hercules, CA) and ZIKV-reactive primers that amplify a region of the NS5 gene and 3′ untranslated region . The QuantStudio 3 Real-Time PCR system (Applied Biosystems, Foster City CA) was used, with amplifying conditions of 50 °C for 15 min, 95 °C for 60 s, followed by 40 cycles of 95 °C for 15 s and 60 °C for 30 s.
Viral RNA was prepared for next-generation sequencing using the Nextera XT library preparation kit (New England BioLabs, Ipswich, MA) following the manufacturer’s instructions. Briefly, viral RNA was quantified, with amounts expressed as genome equivalents (GE/ml), treated with RQ1 DNase (Promega) and used as templates in reverse transcription reactions performed using Superscript IV (ThermoFisher) and random 15-mer primers. Reaction products were treated with RNase H and amplification using the NEBNext® Ultra™ II Q5® Master Mix (New England Biolabs) and random 15-mer primers. Amplified products were fragmented and adapters and indexes added using Nextera XT (Ilumina) with in house Nextera XT dual indexes per manufacturer protocol. The index libraries were then amplified with a KAPA Real-time Library Amplification Kit (KAPA BioSystems, Indianapolis, IN) and purified using Ampure XP beads (Beckman). Final libraries were pooled by Qubit (Invitrogen) concentration and analyzed for size distribution using the Agilent High Sensitivity D1000 Screen Tape on an Agilent Tapestation 2200. Final quantification was performed using the NEBNext® Library Quant Kit for Illumina® (New England Biolabs), and 300 nt pair-end reads were generated using the Illumina MiSeq at the CSU Next-Generation Sequencing Facility.
Next-generation sequencing data were analyzed using a custom Snakemake  workflow, designated as “RNA virus Population Genetics (RPG) Workflow.” Briefly, Read 1 and Read 2.fastq files from paired-end Illumina MiSeq data were trimmed to remove adaptor sequences and quality trimming of phred scores < 30 from the 3′ and 5′ read ends using cutadapt . Reads were then mapped to the ZIKV-I reference sequence using MOSAIK as previously described [40, 41]. Picard (BroadInstitute), Genome Analysis Toolkit (GATK) (BroadInstitute), and SAMtools  were used for variant calling preprocessing. Single-nucleotide variants (SNVs), insertions and deletions (INDELS) were called using LoFreq with the –call-indels command; otherwise, default settings were used . Consensus sequences were generated using the.vcf files generated above and VCFtools .
Markers of genetic diversity were assessed using the RPG Workflow analysis scripts. Richness was calculated by the sum of the unique (not found in the virus input ZIKV-I population) intrahost single-nucleotide variant (iSNV) sites detected in the coding sequence (CDS) in each population. FST was used to estimate genetic divergence between two viral populations as described by Fumagalli et al. :
where pi,s, pj,s and ps are the frequencies of the input ZIKV consensus nucleotide at sites from populations i, j, and combined, respectively. Only LoFreq iSNVs were used, all other sites p = 1. The estimate of FST for the CDS of m sites is:
MIRs, DIRs, SGIRs and TRs were converted to the square root of the arcsine and analyzed by Kruskal-Wallis test and Mann-Whitney U-test. Differences were considered significant when P ≤ 0.05. Comparisons of RNA copy numbers among viruses, organs and days post-infection were performed by Friedman test. The Kruskal-Wallis test and Mann-Whitney U-test were used to compare RNA copy numbers in different organs collected on the same day post-infection. Data were transformed to log10, and results were considered significant when significance was < 5%. The statistical analyses were performed using SPSS version 22 (IBM Corporation Armonk, NY). ZIKV-6V and ZIKV-6A were not evaluated because a small number of mosquitoes tested positive for these viruses, prohibiting a robust statistical analysis.
Genetic characterization of the host-cell adapted viruses
ZIKV was sequentially passaged six times in Vero or C6/36 cells or alternately between both cell types, and the viruses harvested from the final cultures were designated as ZIKV-6V, ZIKV-6C and ZIKV-6A, respectively. The genomes of the input virus (designated ZIKV-I) and viruses collected after the final passage were sequenced (Genbank Accession Nos. MT507047 to MT507050). ZIKV-6C contained two non-synonymous substitutions when compared to ZIKV-I, and these were located in the prM and NS1 protein genes (Table 1). The genome sequences of ZIKV-6V and ZIKV-6A were identical to each other and both contained two non-synonymous substitutions when compared to ZIKV-I. These substitutions were located in the capsid (C) and NS2A protein genes. The estimated reads per sample was between 600,000 and 2,500,000 with an average coding of 5621 (ZIKV-0), 6926 (ZIKV-1), 16110 (ZIKV-2), and 1276 (ZIKV-3) for each nucleotide position.
Intrahost genetic diversity
An intrahost frequency analysis was performed on ZIKV-I and the three final passage viruses (Fig. 1a). Standing variation was detected in ZIKV-I, with 12 iSNVs identified, mostly in the 5′ half of the genome. Four iSNVs identified in ZIKV-6C were absent from viruses passaged in Vero cells and were located at prM-R116G, prM-K132N, NS1-E146K and NS1-K265E (Fig. 1c, d, f, g). prM-K132N and NS1-E146K became fixed. Two sSNVs unique to viruses passaged in Vero cells (ZIKV-6V and ZIKV-6A) also became fixed and were located at capsid-I50T and NS2A-A117V (Figs. 1b and 1h). The richness of ZIKV-6C, ZIKV-6V and ZIKV-6A (Fig. 1i) and their divergence from ZIKV-I (Fig. 1j) was also assessed. Novel variants accumulated more often when passaged exclusively on Vero cells (n = 9) or by alternating between both cell types (n = 11). The fixation index (FST) of ZIKV-6C, ZIKV-6V and ZIKV-6A were estimated and compared to the FST of ZIKV-I (Fig. 1j). ZIKV-6V and ZIKV-6A diverged from ZIKV-I to a similar degree, with FST values of 0.62 and 0.58 respectively, while ZIKV-6C diverged less (FST = 0.36). No iSNVs were identified in the 5′ and 3′ UTRs.
Comparison of infection, dissemination and transmission rates
Select organs and saliva were collected from 1148 Ae. aegypti females that had been orally challenged with ZIKV and were tested for virus by plaque assay. These data were used to calculate midgut infection, salivary gland infection, dissemination and transmission rates (Table 2). MIRs reached as high as 0.97 and 0.83 for mosquitoes infected with ZIKV-I and ZIKV-6C, respectively. In contrast, the MIRs for mosquitoes infected with ZIKV-6V and ZIKV-6A were considerably lower and never exceeded 0.27 and 0.10, respectively. The results were analyzed by Kruskal-Wallis test, revealing that the MIRs for ZIKV-I and ZIKV-6C were significantly higher than ZIKV-6V at the end of the 21 dpi (X2 = 13.574, df = 3, P ≤ 0.05; Fig. 2). The MIR for ZIKV-6A at 21 dpi was not calculated because the MIRs for the earlier times were low.
Disseminated infections were detected in all mosquitoes with midguts infected with ZIKV-I, ZIKV-6C and ZIKV-6V at 21 dpi; therefore, the DIRs were calculated as 1.0 (Table 1). Successful dissemination of all viruses occurred at 14 dpi, with ZIKV-I, ZIKV-6C, ZIKV-6V and ZIKV-6A yielding DIRs of 0.89, 1.0, 1.0 and 0.67, respectively. Kruskal-Wallis test revealed that there were no statistical differences in DIRs between the different groups of mosquitoes at 14 and 21 dpi (X2 = 3.280, df = 2, P ≥ 0.05).
SGIRs as high as 0.95, 0.86, 1.0 and 0.66 were reported for mosquitoes infected with ZIKV-I, ZIKV-6C, ZIKV-6V and ZIKV-6A, respectively (Table 1). The SGIRs for the different groups of mosquitoes were not significantly different at 14 and 21 dpi (X2 = 1.129, df = 3, P ≥ 0.05 and X2 = 1.858, df = 2, P ≥ 0.05, respectively; Fig. 3).
Most saliva samples did not contain detectable virus, with TRs not exceeding 0.17 and 0.05 for mosquitoes challenged with ZIKV-I and ZIKV-6C, respectively (Table 1). Mann-Whitney U-test revealed that there were no significance differences in the TRs for these two groups of mosquitoes at 14 and 21 dpi (Z = 1.449, P ≥ 0.05 and Z = − 0.062, P ≥ 0.05, respectively; Fig. 3). Saliva from mosquitoes infected with the other two viruses were not tested because of the low numbers of mosquitoes (no more than five) that had salivary gland infections, preventing a reliable statistical analysis from being performed.
Viral RNA copy numbers in Ae. aegypti organs
Total RNA was extracted from select organs and other body parts of Ae. aegypti that had been infected with ZIKV-I or ZIKV-6C an then tested by quantitative RT-PCR to measure the amounts of genomic RNA. Evidence of infection was detected in all organs and other body parts (midguts, ovaries, wings/legs, salivary glands, and Malpighian tubules) tested (Fig. 4). Mosquitoes infected with ZIKV-6V and ZIKV-6A were not tested because the viruses rarely infected organs other than the midgut. Likewise, no saliva samples were tested by quantitative RT-PCR because most did not contain detectable amounts of infectious virus.
For mosquitoes infected with ZIKV-I, there was a significant statistical difference in viral RNA copy numbers in midguts at 3 to 21 dpi, with copies numbers highest at the first and last time points (See Additional file 1: Table S1). There was also a significant difference in viral RNA copy numbers in Malpighian tubules, ovaries, and wings/legs at 7 to 21 dpi, with copy numbers highest at the final time point. In salivary glands, there was a significant difference in viral RNA copy numbers at 14 and 21 dpi, with copy numbers highest at the final time point. Similar findings were observed in mosquitoes infected with ZIKV-6C, with a significant statistical difference in viral RNA copy numbers in the midguts and all other body parts at 3 to 21 and 7 to 21 dpi, respectively (See Additional file 1: Table S2). On each occasion, copy numbers were highest at 21 dpi.
Viral RNA copy numbers were compared according to body part and time in mosquitoes infected with ZIKV-I and ZIKV-6C (Fig. 4). Significantly higher amounts of genomic RNA were detected in the ZIKV-I group on several occasions: midguts at 3 dpi (Z = − 5.723, P = 0.000), ovaries at 14 dpi (Z = − 2.103, P = 0.035) and Malpighian tubules at 21 dpi (Z = − 3.070, P = 0.002). Significantly higher amounts of genomic RNA were also detected in the ZIKV-6C group on several occasions: salivary glands at 14 dpi (Z = − 2.639, P = 0.008), salivary glands at 21 dpi (Z = − 3.085, P = 0.002) and wings/legs at 21 dpi (Z = − 2.103, P = 0.035).
RNA viruses have faster evolutionary rates than DNA viruses because they have shorter replication times and larger population sizes and their RNA-dependent RNA polymerases are error-prone and lack proofreading activity . However, RNA viruses that alternate between disparate (vertebrate and invertebrate) hosts typically experience lower rates of evolution than RNA viruses that replicate in a single host [21, 22]. It has been hypothesized that alternate host replication imposes a fitness trade-off that constrains genome evolution and host-specific adaptation. Here, we provide evidence that ZIKV adopts a compromised fitness level for alternate replication in vertebrate and mosquito cells. When ZIKV was released from alternate host cell replication by sequential in vitro passaging, the mosquito cell-adapted virus exhibited similar in vivo fitness in its principal mosquito vector compared to the input virus and greater in vivo fitness than the alternate-passaged virus. ZIKV passaged only in Vero cells were also released from alternate host cell replication. The vertebrate cell-adapted virus exhibited lower in vivo fitness in Ae. aegypti compared to both the input virus as the virus sequentially passaged in mosquito cell cultures. Others have also shown diminished arbovirus replication in mosquito cells or mosquitoes after sequential vertebrate cell passage [28, 34, 45]. For example, sequential passage of VEEV in BHK-21 cells yielded virus that exhibited significantly reduced in vivo fitness in both mosquitoes and vertebrates compared to both the parental (input) virus and the virus maintained by alternating cell culture passages . In contrast, VEEV sequentially passaged in mice produced significantly higher viremias compared to the parental virus as well as virus maintained by alternating passages in mice and Ae. aegypti, underscoring the limitations of in vitro systems. Our study is the first to investigate how sequential and alternating host passaging effects on the in vivo fitness of ZIKV in mosquitoes.
Several other studies have examined the fitness implications of releasing flaviviruses from alternate host passaging [23,24,25,26,27,28]. Serial passaging of WNV and St. Louis encephalitis virus (SLEV) in C6/36 cells increased the relative fitness and replicative ability of each virus in the mosquito cell line, but not an avian cell line . Subsequent in vivo experiments revealed that serial passaging of WNV and SLEV in Culex pipiens increased the replicative ability of WNV, but not SLEV, in this mosquito species [25, 26]. The consequences of releasing arboviruses from alternate host cycling has also been studied using viruses in the families Phenuiviridae, Rhabdoviridae and Togaviridae [29,30,31,32,33]. For example, eastern equine encephalitis virus passaged exclusively in BHK or C6/36 cells exhibited reduced fitness in the bypassed cell line but replicated no more efficiently in the non-bypassed cell line compared to the alternate-passaged virus . Although other studies have compared the replicative abilities of host-cell adapted and alternate-passaged arboviruses, there are several unique aspects to our study. We are the first to investigate the in vivo fitness implications of releasing an Aedes-borne flavivirus from alternate host cycling. Intrahost frequency analysis has rarely been used to determine the genetic composition of a flavivirus released from alternate host passaging [41, 46]. Additionally, geographically matched virus and Aedes mosquitoes were used, with the isolate having a low-passage history and the mosquitoes having been maintained in the insectary for only a few generations to minimize the effects of laboratory adaptation.
Full-genome sequencing identified two nonsynonymous substitutions in ZIKV after serial passage in C6/36 cells when compared to the input virus. One substitution (NS1-E146K) was conservative, and the other (prM-K132N) was non-conservative. Infection, dissemination and transmission rates were not significantly different in Ae. aegypti challenged with ZIKV-6C when compared to ZIKV-I. Our findings indicate that these substitutions do not have a favorable or detrimental effect on viral fitness in mosquitoes. Another, albeit less likely, explanation is that one substitution could be deleterious while the other is advantageous, but together they function as compensatory mutations with a net neutral fitness compared to the initial strain . Compensatory mutations have previously been identified in other flaviviruses, including dengue, Japanese encephalitis and yellow fever viruses [48,49,50]. Two additional mutations were detected in ZIKV-6C (prM-R116G and NS1-K265E) by intrahost analysis but neither had become fixed, although the frequency of each increased after passages.
Serial passaging of ZIKV in Vero cells yielded two fixed, nonsynonymous substitutions. One substitution (NS2A-A117V) was conservative, and the other (C-I50T) was non-conservative. Both substitutions were also detected in the virus generated by six alternating passages. These mutations significantly decreased ZIKV fitness in Ae. aegypti. A non-fixed variant (E-S148G) that increased in frequency after cell culture passages was also identified in both the vertebrate cell-adapted and alternate-passaged viruses. The NS2A-A117V substitution has previously been detected in ZIKV . Site-directed mutagenesis and in vivo experimental infection studies revealed that this single change enhances ZIKV virulence in mice. This change also reduces the host innate immune responses and delays the induction of apoptosis in ZIKV-infected human lung adenocarcinoma (A549) cells. The effects, if any, of the substitutions that occurred in the capsid and envelope proteins are not known. The capsid protein is required for the assembly of the nucleocapsid that surrounds and protects the flavivirus genome, but also has many other roles including apoptosis induction, lipid droplet binding, functioning as an RNA chaperone and suppression of anti-viral RNA silencing [52,53,54]. The envelope protein is the main structural component of flavivirus particle and participates in key steps during the viral life cycle, such as virus-mediated membrane fusion, receptor recognition and entry in susceptible host cells . Site-directed mutagenesis and in vivo experimental infection studies would be required to determine the potential effect of these mutations on the vector competence of Ae. aegypti for ZIKV.
Infection and dissemination rates of ZIKV-6A in Ae. aegypti were low, even though the virus was maintained by alternate passaging in vertebrate and mosquito cells. This finding was not unexpected because several studies have shown large variations in the vector competence of Ae. aegypti from Mexico for ZIKV [56,57,58]. For example, Garcia-Luna and colleagues reported transmission rates of 8 to 51% using Ae. aegypti collected from different geographic regions in Mexico . Flavivirus susceptibility in Ae. aegypti is dependent on genetic differences and gene flow between mosquito populations and environmental factors [59,60,61]. Ae. aegypti from Mexico are more susceptible to infection with ZIKV isolates from the African genotype than the Asian/American genotype . A notable genetic difference between isolates from these genotypes is that all from the Asian/American group encode an N-linked glycosylation site in the envelope protein at E-154N while most from the African group do not [62, 63]. This glycosylation site has been linked to increased ZIKV pathogenesis in mice and implicated in the occurrence of complications in patients infected during recent outbreaks [63,64,65].
The vector competence experiments, which were performed using temporally and spatially matched mosquitoes and virus, revealed that our mosquitoes were inefficient vectors of ZIKV, with TRs of no more than 17%. Other vector competence studies performed using Ae. aegypti and ZIKV strains from the Western Hemisphere have reported TRs exceeding 50% [56, 66,67,68,69]. Although vector competence may vary according to virus dose and the route of inoculation, some of the aforementioned studies used similar doses and the same inoculation route as used in our study. Low TRs have also been documented. For example, in a study where Ae. aegypti from Brazil, the USA and the Dominican Republic were orally challenged with a ZIKV isolate from Mexico, TRs of 0%, 0% and 20% were reported . Garcia-Luna and colleagues revealed considerable variation in the vector competence of Ae. aegypti from different geographic locations within Mexico, with TRs of 14 to 52% at 14 days post-blood meal . Variations in vector competence could also be attributed to differences in the microbiomes and viromes of mosquito populations [71, 72]. Another consideration is that our tissue homogenates were tested by plaque assay, which has a lower limit of detection than PCR-based assays but has the advantage of specifically detecting infectious virus.
There are several limitations of the study. First, the in vitro experiments were performed using C6/36 cells when it would have been preferable to use a cell line derived from Ae. aegypti because it is the principal urban vector and, more importantly, it would have provided consistency in the mosquito species used for the cell culture and vector competence work. Nevertheless, we provide compelling evidence that cell culture history can affect the outcome of vector competence experiments. It is also important to note that Ae. albopictus is also a competent vector of ZIKV in urban settings and that many other Aedes spp. mosquitoes collected in sylvatic environments have tested positive for the virus [9, 73]. Another consideration is that the mosquito and vertebrate cell lines were maintained at different temperatures (28 °C and 37 °C, respectively), which could have selected for genetic changes in the virus. The inclusion of in vitro experiments where both cell types are maintained at the same temperature (i.e. 32 °C) as done by others would have identified temperature-specific mutations [28, 30]. Viruses collected after each cell culture passage were not titrated at each passage were not considered, and founding populations were not taken into consideration, although we estimate that the cells were inoculated with virus at MOIs well below 0.1 and that there were one or two generations per passage. The study could also have been strengthened by evaluating the vector competence of Ae. aegypti for the viruses collected after every cell culture passage but, due to logistical constraints, it was not feasible to challenge and process the number of mosquitoes needed to complete this work. The absence of in vitro fitness assays is another limitation, and a more robust analysis could have been performed if more than six cell culture passages were used. Viruses were not sequenced after each cell culture passage. Likewise, viruses in mosquito tissues were not sequenced, and the inclusion of these experiments would have identified reversions as well as other mutations that could have affected vector competence. Additionally, the virus had been passaged four times in C6/36 cells prior to the sequential and alternating cell culture passaging experiments. Although we purposefully used a local isolate with a low passage history, we cannot dismiss the possibility that the four initial passages in C6/36 cells had a differential effect on viral infectivity in the diverse cell types used in the studies described in this report.
Despite the limitations of our experimental design, we provide evidence that the passage history of ZIKV in cell culture can influence the outcome of vector competence experiments, with mosquito cell-adapted ZIKV exhibiting greater fitness in its natural vector compared to its vertebrate cell-adapted and alternate-passaged counterparts. No statistical differences were observed among the infection, dissemination and transmission rates of the input and mosquito cell-adapted viruses in Ae. aegypti suggesting that infectivity was at best maintained by sequential passage in mosquito cells while the vertebrate cell-adapted and alternate-passaged viruses apparently lost infectivity. Studies on the adaptation of arboviruses to local mosquito populations are important because they provide insight into their transmissibility, vector and host range and epidemic potential.
Availability of data and materials
The nucleotide sequences obtained during the current study are available in GenBank under the accession numbers MT507047–MT507050.
Dick GW, Kitchen SF, Haddow AJ. Zika Virus. (I) Isolation and serological specificity. Trans R Soc Trop Med Hyg. 1952;46(5):509–20.
MacNamara FN. Zika virus: a report on three cases of human infection during an epidemic of jaundice in Nigeria. Trans R Soc Trop Med Hyg. 1954;48(2):139–45.
Boorman J, Draper C. Isolations of arboviruses in the Lagos area of Nigeria, and survey of antibodies to them in man and animals. Trans R Soc Trop Med Hyg. 1968;62(2):269–77.
Simpson DIH. Zika virus infection in man. Trans R Soc Trop Med Hyg. 1964;58(4):335–7.
Duffy MR, Chen TH, Hancock WT, Powers AM, Kool JL, Lanciotti RS, et al. Zika virus outbreak on Yap Island, Federated States of Micronesia. N Engl J Med. 2009;360(24):2536–43.
Musso D, Gubler DJ. Zika virus. Clin Microbiol Rev. 2016;29(3):487–524.
Cao-Lormeau VM, Roche C, Teissier A, Robin E, Berry A-L, Mallet H-P, et al. Zika virus, French Polynesia, South Pacific, 2013. Emerg Infect Dis. 2014;20(6):1085–6.
Epelboin Y, Talaga S, Epelboin L, Dusfour I. Zika virus: an updated review of competent or naturally infected mosquitoes. PLoS Negl Trop Dis. 2017;11(11):1–22.
Boyer S, Calvez E, Chouin-Carneiro T, Diallo D, Failloux AB. An overview of mosquito vectors of Zika virus. Microbes Infect. 2018;20(11–12):646–60.
Souza-Neto JA, Powell JR, Bonizzoni M. Aedes aegypti vector competence studies: a review. Infect Genet Evol. 2019;67(September):191–209. https://doi.org/10.1016/j.meegid.2018.11.009.
Petersen LR, Jamieson DJ, Powers AM, Honein MA. Zika virus. New Engl J Med. 2016;374:1552–63.
Pardy RD, Richer MJ. Zika virus pathogenesis: from early case reports to epidemics. Viruses. 2019;11(10):1–9.
Weaver SC, Costa F, Garcia-Blanco MA, Ko AI, Ribeiro GS, Saade G, et al. Zika virus: history, emergence, biology, and prospects for control. Antivir Res. 2016;130:69–80. https://doi.org/10.1016/j.antiviral.2016.03.010.
Venkataraman S, Prasad BVLS, Selvarajan R. RNA dependent RNA polymerases: insights from structure, function and evolution. Viruses. 2018;10(2):1–23.
Drake JW, Holland JJ. Mutation rates among RNA viruses. Proc Natl Acad Sci USA. 1999;96(24):13910–3.
Ebel GD, Carricaburu J, Young D, Bernard KA, Kramer LD. Genetic and phenotypic variation of West Nile virus in New York, 2000–2003. Am J Trop Med Hyg. 2004;71(4):493–500.
Brault AC, Huang CYH, Langevin SA, Kinney RM, Bowen RA, Ramey WN, et al. A single positively selected West Nile viral mutation confers increased virogenesis in American crows. Nat Genet. 2007;39(9):1162–6.
Brault AC, Powers AM, Ortiz D, Estrada-Franco JG, Navarro-Lopez R, Weaver SC. Venezuelan equine encephalitis emergence: enhanced vector infection from a single amino acid substitution in the envelope glycoprotein. Proc Natl Acad Sci USA. 2004;101(31):11344–9.
Greene IP, Paessler S, Austgen L, Anishchenko M, Brault AC, Bowen RA, et al. Envelope glycoprotein mutations mediate equine amplification and virulence of epizootic Venezuelan equine encephalitis virus. J Virol. 2005;79(14):9128–33.
Tsetsarkin KA, Vanlandingham DL, McGee CE, Higgs S. A single mutation in Chikungunya virus affects vector specificity and epidemic potential. PLoS Pathog. 2007;3(12):1895–906.
Weaver SC, Rico-Hesse R. Genetic diversity and slow rates of evolution in new world alphaviruses. In: Holland JJ, editor. Genetic diversity of RNA viruses current topics in microbiology and immunology. Berlin: Springer; 1992.
Jenkins GM, Rambaut A, Pybus OG, Holmes EC. Rates of molecular evolution in RNA viruses: a quantitative phylogenetic analysis. J Mol Evol. 2002;54(2):156–65.
Chen WJ, Wu HR, Chiou SS. E/NS1 modifications of dengue 2 virus after serial passages in mammalian and/or mosquito cells. Intervirology. 2003;46(5):289–95.
Ciota AT, Lovelace AO, Ngo KA, Le AN, Maffei JG, Franke MA, et al. Cell-specific adaptation of two flaviviruses following serial passage in mosquito cell culture. Virology. 2007;357(2):165–74.
Ciota AT, Lovelace AO, Jia Y, Davis LJ, Young DS, Kramer LD. Characterization of mosquito-adapted West Nile virus. J Gen Virol. 2008;89(7):1633–42.
Ciota AT, Jia Y, Payne AF, Jerzak G, Davis LJ, Young DS, et al. Experimental passage of St. Louis encephalitis virus in vivo in mosquitoes and chickens reveals evolutionarily significant virus characteristics. PLoS ONE. 2009;4(11):e7876.
Deardorff ER, Fitzpatrick KA, Jerzak GVS, Shi PY, Kramer LD, Ebel GD. West nile virus experimental evolution in vivo and the trade-off hypothesis. PLoS Pathog. 2011;7(11):e1002335.
Vasilakis N, Deardorff ER, Kenney JL, Rossi SL, Hanley KA, Weaver SC. Mosquitoes put the brake on arbovirus evolution: experimental evolution reveals slower mutation accumulation in mosquito than vertebrate cells. PLoS Pathog. 2009;5(6):e1000467.
Novella IS, Hershey CL, Escarmis C, Domingo E, Holland JJ. Lack of evolutionary stasis during alternating replication of an arbovirus in insect and mammalian cells. J Mol Biol. 1999;287(3):459–65.
Weaver SC, Brault AC, Kang W, Holland JJ. Genetic and fitness changes accompanying adaptation of an arbovirus to vertebrate and invertebrate cells. J Virol. 1999;73(5):4316–26.
Greene IP, Wang E, Deardorff ER, Milleron R, Domingo E, Weaver SC. Effect of alternating passage on adaptation of sindbis virus to vertebrate and invertebrate cells. J Virol. 2005;79(22):14253–60.
Arias-Goeta C, Moutailler S, Mousson L, Zouache K, Thiberge JM, Caro V, et al. Chikungunya virus adaptation to a mosquito vector correlates with only few point mutations in the viral envelope glycoprotein. Infect Genet Evol. 2014;24:116–26. https://doi.org/10.1016/j.meegid.2014.03.015.
Moutailler S, Roche B, Thiberge JM, Caro V, Rougeon F, Failloux AB. Host alternation is necessary to maintain the genome stability of Rift Valley fever virus. PLoS Negl Trop Dis. 2011;5(5):e1156.
Coffey LL, Vasilakis N, Brault AC, Powers AM, Tripet FR, Weaver SC. Arbovirus evolution in vivo is constrained by host alternation. Proc Natl Acad Sci USA. 2008;105(19):6970–5.
Beaty BJ, Calisher C, Shope R. Arboviruses. In: Diagnostic procedures for viral, rickettsial, and chlamydial infections. 6th ed. Washington, D.C: American Public Health Association; 1989. p. 797–855.
Smith DR, Carrara AS, Aguilar PV, Weaver SC. Evaluation of methods to assess transmission potential of Venezuelan equine encephalitis virus by mosquitoes and estimation of mosquito saliva titers. Am J Trop Med Hyg. 2005;73(1):33–9.
Xu MY, Liu SQ, Deng CL, Zhang QY, Zhang B. Detection of Zika virus by SYBR green one-step real-time RT-PCR. J Virol Methods. 2016;236:93–7. https://doi.org/10.1016/j.jviromet.2016.07.014.
Köster J, Rahmann S. Corrigendum: Snakemake—a scalable bioinformatics workflow engine (Bioinformatics (2012) 28:19 (2520–2522) https://doi.org/10.1093/bioinformatics/bts480). Bioinformatics. 2018;34(20):3600.
Wilm A, Aw PPK, Bertrand D, Yeo GHT, Ong SH, Wong CH, et al. LoFreq: a sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets. Nucleic Acids Res. 2012;40(22):11189–201.
Lee WP, Stromberg MP, Ward A, Stewart C, Garrison EP, Marth GT. MOSAIK: a hash-based algorithm for accurate next-generation sequencing short-read mapping. PLoS ONE. 2014;9(3):e90581.
Grubaugh ND, Fauver JR, Rückert C, Weger-Lucarelli J, Garcia-Luna S, Murrieta RA, et al. Mosquitoes transmit unique West Nile virus populations during each feeding episode. Cell Rep. 2017;19(4):709–18.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, et al. The variant call format and VCFtools. Bioinformatics. 2011;27(15):2156–8.
Fumagalli M. Assessing the effect of sequencing depth and sample size in population genetics inferences. PLoS ONE. 2013;8(11):14–7.
Moser LA, Boylan BT, Moreira FR, Myers LJ, Svenson EL, Fedorova NB, et al. Growth and adaptation of Zika virus in mammalian and mosquito cells. PLoS Negl Trop Dis. 2018;12(11):1–19.
Grubaugh ND, Smith DR, Brackney DE, Bosco-Lauth AM, Fauver JR, Campbell CL, et al. Experimental evolution of an RNA virus in wild birds: evidence for host-dependent impacts on population structure and competitive fitness. PLoS Pathog. 2015;11(5):1–19.
Storz JF. Compensatory mutations and epistasis for protein function. Curr Opin Struct Biol. 2018;50:18–25. https://doi.org/10.1016/j.sbi.2017.10.009.
Li X-D, Deng C-L, Ye H-Q, Zhang H-L, Zhang Q-Y, Chen D-D, et al. Transmembrane domains of NS2B contribute to both viral RNA replication and particle formation in Japanese encephalitis virus. J Virol. 2016;90(12):5735–49.
Kummerer BM, Rice CM. Mutations in the yellow fever virus nonstructural protein NS2A selectively block production of infectious particles. J Virol. 2002;76(10):4773–84.
Naik NG, Wu H-N. Mutation of putative N-glycosylation sites on dengue virus NS4B decreases RNA replication. J Virol. 2015;89(13):6746–60.
Ávila-Pérez G, Nogales A, Park JG, Márquez-Jurado S, Iborra FJ, Almazan F, et al. A natural polymorphism in Zika virus NS2A protein responsible of virulence in mice. Sci Rep. 2019;9(1):1–17.
Oliveira ERA, Mohana-Borges R, de Alencastro RB, Horta BAC. The flavivirus capsid protein: structure, function and perspectives towards drug design. Virus Res. 2017;227:115–23. https://doi.org/10.1016/j.virusres.2016.10.005.
Freire JM, Santos NC, Veiga AS, Da Poian AT, Castanho MARB. Rethinking the capsid proteins of enveloped viruses: multifunctionality from genome packaging to genome transfection. FEBS J. 2015;282(12):2267–78.
Sotcheff S, Routh A. Understanding flavivirus capsid protein functions: the tip of the iceberg. Pathogens. 2020;9(1):42.
Zhang X, Jia R, Shen H, Wang M, Yin Z, Cheng A. Structures and functions of the envelope glycoprotein in flavivirus infections. Viruses. 2017;9(11):1–14.
Garcia-Luna SM, Weger-Lucarelli J, Rückert C, Murrieta RA, Young MC, Byas AD, et al. Variation in competence for ZIKV transmission by Aedes aegypti and Aedes albopictus in Mexico. PLoS Negl Trop Dis. 2018;12(7):1–21.
Weger-Lucarelli J, Rückert C, Chotiwan N, Nguyen C, Garcia Luna SM, Fauver JR, et al. Vector competence of American mosquitoes for three strains of Zika virus. PLoS Negl Trop Dis. 2016;10(10):1–16.
Lozano-Fuentes S, Kenney JL, Varnado W, Byrd BD, Burkhalter KL, Savage HM. Susceptibility and vectorial capacity of American Aedes albopictus and Aedes aegypti (Diptera: Culicidae) to American Zika virus strains. J Med Entomol. 2019;56(1):233–40.
Gorrochotegui-Escalante N, De Lourdes MM, Fernandez-Salas I, Beaty BJ, Black WC IV. Genetic isolation by distance among Aedes aegypti populations along the northeastern coast of Mexico. Am J Trop Med Hyg. 2000;62(2):200–9.
Gorrochotegui-Escalante N, Gomez-Machorro C, Lozano-Fuentes S, Fernandez-Salas I, Munoz MDL, Farfan-Ale JA, et al. Breeding structure of Aedes aegypti populations in Mexico varies by region. Am J Trop Med Hyg. 2002;66(2):213–22.
Black WC IV, Bennett KE, Gorrochótegui-Escalante N, Barillas-Mury CV, Fernández-Salas I, Muñoz MDL, et al. Flavivirus susceptibility in Aedes aegypti. Arch Med Res. 2002;33(4):379–88.
Faye O, Freire CCM, Iamarino A, Faye O, de Oliveira JVC, Diallo M, et al. Molecular evolution of Zika virus during its emergence in the 20th century. PLoS Negl Trop Dis. 2014;8(1):36.
Haddow AD, Schuh AJ, Yasuda CY, Kasper MR, Heang V, Huy R, et al. Genetic characterization of Zika virus strains: geographic expansion of the Asian lineage. PLoS Negl Trop Dis. 2012;6(2):e1477.
Lanciotti RS, Kosoy OL, Laven JJ, Velez JO, Lambert AJ, Johnson AJ, et al. Genetic and serologic properties of Zika virus associated with an epidemic, Yap State, Micronesia, 2007. Emerg Infect Dis. 2008;14(8):1232–9.
Carbaugh DL, Baric RS, Lazear HM. Envelope protein glycosylation mediates Zika virus pathogenesis. J Virol. 2019;93(12):1–16.
Dutra HLC, Rocha MN, Dias FBS, Mansur SB, Caragata EP, Moreira LA. Wolbachia blocks currently circulating Zika virus isolates in Brazilian Aedes aegypti mosquitoes. Cell Host Microbe. 2016;19(6):771–4. https://doi.org/10.1016/j.chom.2016.04.021.
Kenney JL, Romo H, Duggal NK, Tzeng WP, Burkhalter KL, Brault AC, et al. Transmission incompetence of Culex quinquefasciatus and Culex pipiens pipiens from North America for Zika virus. Am J Trop Med Hyg. 2017;96(5):1235–40.
Main BJ, Nicholson J, Winokur OC, Steiner C, Riemersma KK, Stuart J, et al. Vector competence of Aedes aegypti, Culex tarsalis, and Culex quinquefasciatus from California for Zika virus. PLoS Negl Trop Dis. 2018;12(6):1–13.
Fernandes RS, O’connor O, Bersot MIL, Girault D, Dokunengo MR, Pocquet N, et al. Vector competence of Aedes aegypti, Aedes albopictus and Culex quinquefasciatus from Brazil and New Caledonia for three Zika virus lineages. Pathogens. 2020;9(7):1–17.
Roundy CM, Azar SR, Rossi SL, Huang JH, Leal G, Yun R, et al. Variation in Aedes aegypti mosquito competence for Zika virus transmission. Emerg Infect Dis. 2017;23(4):625–32.
Dennison NJ, Jupatanakul N, Dimopoulos G. The mosquito microbiota influences vector competence for human pathogens. Curr Opin Insect Sci. 2014;3:6–13.
Azar SR, Weaver SC. Vector competence: what has Zika virus taught us? Viruses. 2019;11(9):1–24.
McKenzie BA, Wilson AE, Zohdy S. Aedes albopictus is a competent vector of Zika virus: a meta-analysis. PLoS ONE. 2019;14(5):1–16.
The authors would like to thank the staff of Arbovirology Laboratory from Universidad Autónoma de Yucatán and the staff of Arthropod-Borne and Infectious Diseases Laboratory, Department of Microbiology, Immunology and Pathology, Colorado State University, for the support and assistance during this study. Additionally, Talavera-Aguilar would like to thank Consejo Nacional de Ciencia y Tecnología (CONACyT) for the doctoral fellowship awarded.
This research was supported by Bio & Medical Technology Development Program of the National Research Foundation (NRF) funded by the Korean government (MSIT) (no. 2016M3A9B8942289) (CMW, SMK) and partially by a grant from the US National Institute of Allergy and Infectious Disease, National Institutes of Health, under grant AI067380 (GDE).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Comparison of genomic RNA copy numbers in organs and other body parts of Ae. aegypti challenged with input Zika virus (ZIKV-I). Table S2. Comparison of genomic RNA copy numbers in organs and other body parts of Ae. aegypti challenged with C6/36 cell-adapted Zika virus (ZIKV-6C).
About this article
Cite this article
Talavera-Aguilar, L.G., Murrieta, R.A., Kiem, S. et al. Infection, dissemination, and transmission efficiencies of Zika virus in Aedes aegypti after serial passage in mosquito or mammalian cell lines or alternating passage in both cell types. Parasites Vectors 14, 261 (2021). https://doi.org/10.1186/s13071-021-04726-1