- Open Access
Zika virus and temperature modulate Elizabethkingia anophelis in Aedes albopictus
Parasites & Vectors volume 14, Article number: 573 (2021)
Vector-borne pathogens must survive and replicate in the hostile environment of an insect’s midgut before successful dissemination. Midgut microbiota interfere with pathogen infection by activating the basal immunity of the mosquito and by synthesizing pathogen-inhibitory metabolites.
The goal of this study was to assess the influence of Zika virus (ZIKV) infection and increased temperature on Aedes albopictus midgut microbiota. Aedes albopictus were reared at diurnal temperatures of day 28 °C/night 24 °C (L) or day 30 °C/night 26 °C (M). The mosquitoes were given infectious blood meals with 2.0 × 108 PFU/ml ZIKV, and 16S rRNA sequencing was performed on midguts at 7 days post-infectious blood meal exposure.
Our findings demonstrate that Elizabethkingia anophelis albopictus was associated with Ae. albopictus midguts exposed to ZIKV infectious blood meal. We observed a negative correlation between ZIKV and E. anophelis albopictus in the midguts of Ae. albopictus. Supplemental feeding of Ae. albopictus with E. anophelis aegypti and ZIKV resulted in reduced ZIKV infection rates. Reduced viral loads were detected in Vero cells that were sequentially infected with E. anophelis aegypti and ZIKV, dengue virus (DENV), or chikungunya virus (CHIKV).
Our findings demonstrate the influence of ZIKV infection and temperature on the Ae. albopictus microbiome along with a negative correlation between ZIKV and E. anophelis albopictus. Our results have important implications for controlling vector-borne pathogens.
Zika virus (ZIKV), a Flavivirus of the family Flaviviridae with an ~ 11-kb positive-sense, single-stranded RNA genome , is classified as a member of the Spondweni group . There are two main ZIKV lineages: one lineage with two genotypes from Africa and another from Asia [3, 4].
ZIKV is mainly transmitted by the bite of infected Aedes mosquitoes, although sexual, transplacental, and blood transfusion transmissions have also been documented . There is no Food and Drug Administration (FDA)-approved medication or vaccine to treat or prevent ZIKV infection, despite the World Health Organization (WHO) declaring ZIKV a public health emergency of international concern . Medications for symptomatic relief are the only source of relief for infected individuals . Control measures based on preventing sexual transmission and controlling the mosquito vector are the main prevention strategies against ZIKV . An upsurge in insecticide resistance, as well as climate change, has expanded the geographical range of vector territories , highlighting the urgent need for alternative mosquito control approaches.
The interaction of gut microbial taxa and arboviruses of medical importance has been studied extensively [9,10,11,12,13,14,15,16,17] due to the potential for microbes to act as candidates for the development of novel therapeutic and transmission-blocking agents.
Flavobacteria have been shown to dominate the Anopheles mosquito midgut . Elizabethkingia (EBK)—a gram-negative, rod-like, aerobic, non-fermenting, non-motile, and non-spore-forming flavobacterium—is widely distributed in Aedes [18,19,20] and Anopheles [21,22,23] mosquitoes. EBK has been isolated from both field-caught and laboratory-reared mosquitoes in Africa, Europe, and North America [20, 23, 24] and can be transmitted among mosquitoes through vertical, transstadial, or horizontal transmission [18, 22]. In addition, Elizabethkingia meningoseptica metabolites obtained from Anopheles stephensi midguts possess antiparasitic effects, as well as gametocyte toxicity . Similarly, Elizabethkingia anophelis has a significant effect on Plasmodium parasite development, reducing the oocyst load when introduced at a low bacterial dose .
A previous metagenomics study aimed at understanding the effect of ZIKV on the dynamic bacterial community harbored by Aedes aegypti  demonstrated that midgut microbes belonging to the family Flavobacteriaceae were the most abundant bacterial taxa among ZIKV-infected mosquitoes. This was particularly the case among the ZIKV-exposed, blood-fed individuals. Flavobacteriaceae abundance was high during the first days of infection but decreased at 7 days post-infection (dpi) among ZIKV-exposed gravid individuals . The aims of this study were to identify the bacterial species within the Aedes albopictus midgut that significantly correlate with ZIKV infectivity and replication, assess whether these microbes possessed antiviral properties, and characterize the influence of temperature on these interactions. Once identified, the symbionts associated with altered vector competence can potentially inform novel vector-borne pathogen control measures.
The ZIKV strain used in this study was HND 2016–19563 (GenBank accession no. KX906952), which was isolated in 2016 in New York state from an infected patient who had traveled to Honduras . To create stock virus for this study, the virus was isolated and passaged three times on Vero cell culture and once on C6/36 cell culture. Stocks were frozen at −70 °C prior to use. Chikungunya virus (CHIKV) IDR 140025961 was isolated from a patient in the Dominican Republic in 2014, while the dengue virus serotype 2 (DENV-2) NI2B strain 306 was isolated from a patient in Nicaragua in 2007. All virus samples were treated similarly.
Mosquito sampling and screening for ZIKV infection
The Intergovernmental Panel on Climate Change (IPCC) predicts a 2–4 °C mean global temperature rise over the next century due to global warming . We modeled a 2 °C increase (M temperature regimen) over our baseline temperatures (L temperature regimen), representing the mean temperatures during peak transmission in Suffolk County (NY), where the Ae. albopictus were collected in 2015 (kindly provided by I. Rochlin). Ae. albopictus were colonized at the New York State Department of Health (NYSDOH) Arbovirus Laboratory.
The 15th generation following laboratory colonization of Ae. albopictus were vacuum-hatched and maintained at 28 °C under standard rearing conditions.
To initiate the experiment aimed at exposing them to ZIKV, the mosquitoes were hatched and reared under two thermal conditions: 14 h light at 30 °C/10 h dark at 26 °C (M temperature regimen), or 14 h light at 28 °C/10 h dark at 24 °C (L temperature regimen).
At 8 days post-eclosion, female adults were orally exposed to blood meals consisting of either 1:1 dilution of defibrinated sheep blood plus 2.5% sucrose, sodium bicarbonate (to adjust pH to 8.0) and the virus, or a non-infectious blood meal (NIBM) containing a final concentration of 2.5% sucrose solution. Infectious blood meals contained 2.00E+08 PFU/ml ZIKV HND . The female mosquitoes were exposed to blood meals in a 37 °C preheated Hemotek membrane feeding system (Discovery Workshops, Accrington, UK) with a porcine sausage casing membrane. After 1 h, the mosquitoes were anesthetized with CO2 and immobilized on a pre-chilled tray connected to CO2. Engorged females were separated and placed in three separate 0.6 L cardboard cartons (30 individuals/carton). Blood-fed females were maintained on 10% sucrose solution provided ad libitum.
At 7 dpi, the female mosquitoes were immobilized using triethylamine (Sigma-Aldrich, St. Louis, MO, USA). To examine for ZIKV dissemination in the mosquito, the legs were removed before dissecting the gut. To assess transmission, saliva was collected under sterile conditions by inserting the proboscis of each female mosquito into a capillary tube containing ~ 20 µl fetal bovine serum (FBS) plus 50% sucrose 1:1 for 30 min. The mixture was then ejected into 125 µl of mosquito diluent (MD) (20% heat-inactivated FBS in Dulbecco phosphate-buffered saline plus 50 µg/ml penicillin/streptomycin, 50 µg/ml gentamicin, and 2 µg/ml Fungizone [Sigma-Aldrich, St. Louis, MO, USA]). The midgut was dissected from surface-sterilized individual mosquitoes and rinsed twice in sterile phosphate-buffered saline (PBS) before transfer to sterile microcentrifuge tubes and storage at −80 °C until testing. Mosquito carcass and legs were then placed in individual tubes containing 500 µl MD and a bead. All samples were held at −80 °C until assaying. An absolute quantity of ZIKV in the legs, carcass, and saliva was obtained by isolating RNA using TRIzol (Thermo Fisher Scientific, Waltham, MA, USA). A ZIKV-specific quantitative polymerase chain reaction (qPCR) assay targeting the NS1 region was used for detection (Table 1). Serially diluted ZIKV of known quantities was used to generate the standard curve and extrapolate viral load (Additional file 1: Figure S1). Midguts obtained from Ae. albopictus adults maintained on 10% sucrose solution were used as a control.
Sequencing and analysis of microbiome of Aedes mosquito midguts
A total of 29 samples were included in this study: 10 midguts exposed to infected blood meal (IBM), 10 midguts exposed to NIBM, and nine midguts exposed to 10% sucrose only (non-blood-fed [NBF]), along with one negative control (NC), one non-template control (NTC), and one positive spike (PC). Aedes albopictus midgut, legs, carcass, and saliva samples were screened for virus at 7 dpi. The infection and dissemination status of each individual was determined prior to downstream processing. The complementary DNA (cDNA) from individual midguts was prepared from the RNA samples using the iScript cDNA kit (Bio-Rad, Hercules, CA, USA). Thereafter, 16S rRNA was amplified using the 16S primer set and conditions from Table 1.
The barcoded high-throughput amplicon sequencing of the bacterial 16S V3-V4 hypervariable rRNA was completed at Wadsworth Center’s Advanced Genomic Technologies Core (WCAGT). PCR reactions were carried out in a total volume of 50 µl: 5 µM of 16S rRNA V3-V4 hypervariable region primers consisting of an Illumina barcode (Table 1), 2 µl cDNA, 8 µl deionized filter sterilized water, and 36 µl AccuStart II PCR SuperMix (Quanta Biosciences, Beverly, MA, USA). To ensure identification of any potential contamination at either the PCR stage or the sequencing step, two negative controls (water and non-template control) were included at the PCR step and sequenced. A positive spike was included during the sequencing run. A fragment size of ~ 460 base pairs (bp) of each sample, as well as the negative controls, was submitted to the WCAGT for sequencing. Automated cluster generation and paired-end sequencing (250-bp reads) were performed using the Illumina MiSeq 500-cycle reagent kit.
Data analysis was carried out on QIITA (https://qiita.ucsd.edu/) and MicrobiomeAnalyst . In summary, to convert the de-multiplexed FASTQ files to the format used by QIITA, the command split libraries FASTQ was used. The clustering of the sequences into operational taxonomic units (OTUs) was done via closed-reference OTU picking based on the Greengenes 16S reference database, which matches the sequences based on a 97% sequence identity threshold before assigning taxonomies. A BIOM-formatted OTU table was subsequently generated, and in order to reduce the potential for alpha and beta diversity biases, the data were rarefied to a depth of 30,000 reads per sample. Finally, to visualize the taxonomic profiles of each sample, taxa bar plots were generated using the rarefied data artifact. Following the generation of the relative quantity of each taxon, we generated taxa which were deposited in the NCBI GenBank Short Read Archives with accession numbers PRJNA608507, PRJNA608505, PRJNA608503, PRJNA608501, PRJNA608500, PRJNA608499, PRJNA608472, PRJNA608496, PRJNA608495, PRJNA608493, PRJNA608492, PRJNA608491, PRJNA608489, PRJNA608483, PRJNA608471, PRJNA608482, PRJNA608470, PRJNA608456, PRJNA608450, PRJNA608481, PRJNA608449, PRJNA608448, PRJNA608480, PRJNA608479, PRJNA608477, PRJNA608437, PRJNA608434, PRJNA608427, and PRJNA608428.
The alpha diversity profiling (number of species of taxa in a single sample) was calculated using the Chao1 diversity measure, and significance testing was done using analysis of variance (ANOVA) in MicrobiomeAnalyst. The Chao1 index estimates the richness of taxa in a sample by extrapolating the number of rare organisms that may have been omitted due to under-sampling. The index assumes that the number of observations for a taxon has a Poisson distribution and thus corrects for variance . Beta diversity, a comparison of microbial community composition between two samples or two conditions, was analyzed by utilizing the Bray–Curtis dissimilarity statistic . The distance matrix generated after comparing each sample to the other was visualized for dissimilarity between samples using the principal coordinates analysis method.
Phylogenetic analysis of the Elizabethkingia (EBK) isolated from Ae. albopictus midguts
In order to validate the identity of EBK taxa identified in the Ae. albopictus midguts, the rpoB gene (accession number MT096047) was amplified from the midgut cDNA using the rpoB primer pairs and PCR conditions as described in Table 1. The 558-bp PCR amplicon was excised from an agarose gel and cleaned using the QIAquick gel extraction kit (Qiagen, UK). The amplicons were Sanger-sequenced at the WCAGT. To construct the EBK rpoB phylogenetic tree, Molecular Evolutionary Genetic Analysis (MEGA) software  was utilized.
A total of nine gene sequences published with GenBank accession numbers CP023010.2, KY587659.1, CP016377.1, KY587657.1, CP035811.1, CP035809.1, CP011059.1, CP039929.1, and MN327643.1 were added to the phylogenetic analysis.
Clustering and biomarker analysis
Patterns of relative abundance of taxa species in response to different experimental factors were analyzed in MicrobiomeAnalyst. To provide an estimate of the most important taxa for classification of data resulting from different experimental factors, biomarker analysis was carried out using the random forest algorithm within MicrobiomeAnalyst. The default setting of the number of trees to grow and number of predictors to try (500 and 7, respectively) was applied with the randomness setting left on. The random forest algorithm is a supervised classification algorithm of trees created by using bootstrap samples with training data and random feature selection in tree induction. It is an ensemble of unpruned classifications or regression trees trained with the bagging method .
Analysis of the impact of supplemental feeding of Ae. albopictus with E. anophelis aegypti on their survivability and vector competence for ZIKV
The Elizabethkingia strain associated with the Ae. albopictus in this study was identified as E. anophelis albopictus (Fig. 3). Elizabethkingia anophelis aegypti cultured from Ae. aegypti midguts was used for supplemental feeding of Ae. albopictus, which were reared to adult stage as described previously. Septic adult Ae. albopictus were fed a supplemental diet of E. anophelis aegypti grown to optical density at 600 nm (OD600) = 1.0 [(3.16E+08 colony-forming units per milliliter (CFU/ml)] (N = 30) and OD600 = 0.5 (1.0E+08 CFU/ml) (N = 30), after which the survivability was assessed. Sterile PBS was used as a negative control (N = 30), and the differences in survivability were calculated using the Log-rank test.
To assess the impact of E. anophelis aegypti supplemental feeding on Ae. albopictus vector competence for ZIKV, Ae. albopictus mosquitoes were reared to the adult stage, and upon emergence, Ae. albopictus adults were antibiotic-treated with 20% penicillin, 20% streptomycin, and 75ug/ml gentamicin in a 10% sucrose solution for 8 days (Fig. 4b). The cotton strips that were used to deliver the antibiotic treatment were changed daily. Both groups were maintained on filter-sterilized 10% sucrose solution. The effectiveness of the sterility rearing was confirmed by CFU assays on randomly selected adult female midguts prior to blood-feeding or bacterial challenge (Fig. 4c). Elizabethkingia anophelis aegypti was introduced to a total of n = 75 antibiotic-treated Ae. albopictus through a sugar meal (cotton strips moistened with 1.5% sucrose solution containing E. anophelis aegypti to a final concentration of 1.0E+08 CFU/ml) . The control consisted of septic Ae. albopictus (n = 75) without exposure to E. anophelis aegypti. The cotton strips were withdrawn after 24 h and replaced with sterile cotton strips moistened with sterile 10% sucrose solution. Both the experimental (antibiotic-treated Ae. albopictus infected with E. anophelis aegypti) and control groups (septic Ae. albopictus) were deprived of sucrose solution for 18–24 h and offered 2.00E+08 PFU/ml ZIKV blood meal as described previously. The midgut, carcass, saliva, and legs were obtained from each individual of each population at 7 dpi. RNA was isolated from the carcass, saliva, and legs, and screened for ZIKV. The infection, dissemination, and transmission percentage rates were calculated and compared between the two groups using Chi-square tests. The assay was carried out in two independent experiments.
Assessment of the broad antiviral effect of E. anophelis in vitro
The antiviral effect of E. anophelis aegypti was assessed in the Vero and C636 cell lines. A total of 3 ml/well of cells (2.0 × 105 cells/ml) was used to seed six-well plates that were incubated at 37 °C (Vero) and 28 °C (C636) and 5% CO2 until they attained 90% confluence. Upon attaining 90% confluence, the experimental wells were inoculated with 75 μl of 1.0E+08 CFU/ml of E. anophelis aegypti or Escherichia coli, while the control wells were inoculated with either E. anophelis aegypti or E. coli that were heat-inactivated for 6 h at 90 °C. The experiment was performed on two separate days, and each assay was performed in triplicate.
At 24 h post-infection (hpi) with bacteria (E. anophelis aegypti or E. coli), Vero cell monolayers were infected at a multiplicity of infection (MOI) of 0.1 infectious particles/cell with ZIKV, CHIKV, or DENV. C636 cells were inoculated with ZIKV at 24 hpi with E. anophelis aegypti or E. coli. The supernatant was harvested at 24 hpi with the virus. The viral load was compared between the experimental and control supernatants. RNA was extracted from the supernatant using the TRIzol procedure (Thermo Fisher Scientific, Waltham, MA, USA) and the virus quantified by qPCR as described previously, using the primers in Table 1. To assess the effect of timing on virus replication, 75 μl of 8.0 log10 CFU/ml of E. anophelis aegypti or E. coli was co-inoculated with 0.1 MOI of ZIKV. Differences in virus quantities were measured using ANOVA in GraphPad Prism.
In order to assess the cell-independent effect of E. anophelis aegypti antiviral activity, U4.4 cells (Ae. albopictus) were co-infected with E. anophelis aegypti and ZIKV.
The E. anophelis aegypti bacterial strain was streaked on Luria–Bertani (LB) plates and incubated in a 10% CO2 incubator at 28 °C overnight. A single colony was picked using a sterile plastic loop and used to inoculate 500 ml of sterile LB broth. The inoculated broth medium was incubated overnight at 28 °C in an orbital shaker shaking at 225 rpm. The E. anophelis aegypti overnight culture was subsequently diluted to OD600 = 1.0 (8.5 log10 CFU/ml), OD600 = 0.5 (8.0 log10 CFU/ml), and OD600 = 0.2 (7.7 log10 CFU/ml).
Individual wells were seeded with 3 ml of 2 × 105 U4.4 cells/ml suspended in minimum essential medium medium without FBS and incubated at 28 °C for 1 h before adding FBS to each well at a final concentration of 20%. The U4.4 cells were then incubated for 5 days in a 10% CO2 incubator at 28 °C. After 5 days, the supernatant was carefully removed, and the plates were infected with E. anophelis aegypti as described previously at OD600 = 1.0 (8.5 log10 CFU/ml), OD600 = 0.5 (8.0 log10 CFU/ml), and OD600 = 0.2 (7.7 log10 CFU/ml) for three biological replicates. After 1 h of incubation at 28 °C, 3 ml of M&M media containing 20% FBS was added to each well and incubated at 28 °C.
To determine whether the U4.4 cells were successfully infected with E. anophelis aegypti, a single well from each replicate of the different dilutions of E. anophelis aegypti infection was harvested at 2 dpi.
The cells were spun at 1000 rpm and washed three times with sterile PBS. After the third wash, 100 µl of the washed cells was used to inoculate 3 ml of LB media that was incubated overnight at 28 °C in an orbital shaker at 225 rpm. The experimental wells were infected with 0.1 MOI ZIKV as described previously, as were the ZIKV quantities. Supernatants were harvested from the two treatments at both 2 dpi and 4 dpi. RNA was extracted from the cells using TRIzol (Thermo Fisher Scientific, Waltham, MA, USA). The ZIKV quantities were measured as described previously.
Influence of temperature, blood feeding, and ZIKV infection on microbial profile of Ae. albopictus
A total of 12 out of 30 Ae. albopictus reared in the L temperature regimen disseminated the virus, with a virus titer that ranged from 1 to 2 log10 PFU/ml, while 26 out of 30 individuals reared in the M temperature regimen disseminated the virus, with a virus titer that ranged from 1 to 5 log10 PFU/ml (Additional file 5). None of the individuals reared in the L temperature regimen transmitted the virus in the saliva, while three individuals reared in the M temperature regime had positive saliva. The viral titer range in saliva was 1–2 log10 PFU/ml (Additional file 6).
Illumina MiSeq 16S rRNA sequencing resulted in 6,089,629 reads. The number of reads for individual samples ranged from 531,220 to 34,833, with a mean of 152,448. Both the positive spike-in (PhiX) and negative controls (water and non-template control) had a limited number of sequence reads, ranging from 1620 to 2823, and hence were discarded in downstream analysis. In order to control for alpha and beta diversity biases, each sample was rarefied to a total frequency of 34,000 reads. A total of 1878 OTUs were identified in this study. The core phyla identified were Bacteroidetes, Proteobacteria, Actinobacteria, Cyanobacteria, and Firmicutes.
Overall, midgut taxa diversity was impacted by feeding [permutational multivariate analysis of variance (Beta-diversity analysis; PERMANOVA test F-value = 13.22; R-squared = 0.5042; P-value < 0.001; Fig. 1a)]. Increasing the rearing temperature did not result in significant changes in taxa diversity (Beta-diversity analysis; PERMANOVA test; F-value: 2.1707; R-squared: 0.074415; P-value < 0.1; Alpha-diversity index: Chao1: p-value: 0.85872; t-test = 0.17997; P = 0.9; Fig. 1b and Additional file 2: Figure S2). Despite this, at the individual taxa level, Elizabethkingia was unique in that it was identified at higher proportions at the lower rearing temperature regime (Fig. 2b, Additional file 3: Figure S3). Blood meal intake, irrespective of the infection status, decreased the taxon diversity (PERMANOVA; F-value: 20.007; R-squared: 0.42562; P-value < 0.001; Fig. 1c). Additionally, ZIKV infectious blood meal uptake significantly reduced the individual midgut bacterial diversity (Alpha-diversity index: Chao1; p-value: ANOVA test = 24.2116; p < 0.05; Fig. 2a; Additional file 4: Figure S4). Elizabethkingia was significantly associated with the midgut exposure to an infectious blood meal (Fig. 2c).
Phylogenetic analysis of Elizabethkingia
To taxonomically classify the Elizabethkingia strain that we identified in the midgut of the Ae. albopictus, we sequenced the Elizabethkingia rpoB gene amplified from the midgut cDNA of Ae. albopictus. The phylogenetic analysis demonstrated that the Elizabethkingia strain identified in our study has 100% sequence homology to E. anophelis, which was originally isolated from the midgut of Anopheles gambiae . Our strain has 99% homology to E. endophytica, originally isolated from sweet corn , and was closely related to the E. meningoseptica genomospecies 4 , E. miricola (isolated from condensation water in the space station) , and E. bruuniana . Elizabethkingia anophelis has 92% homology to E. meningoseptica, originally isolated from human cerebrospinal fluid  (Fig. 3).
Elizabethkingia anophelis colonization of Ae. albopictus midguts is associated with reduced ZIKV infection rates and viral loads
We assessed the impact of E. anophelis aegypti colonization of Ae. albopictus midguts on survivability and vector competence for ZIKV. Colonization of Ae. albopictus midguts with 3.16E+08 CFU/ml and 1.0E+08 CFU/ml did not impact their survivability [log-rank (Mantel-Cox) test [Chi-square = 2.345; df = 2; P = 0.31] (Fig. 4a). To study the impact of E. anophelis aegypti on the vector competence of Ae. albopictus for ZIKV, conventionally reared Ae. albopictus was utilized as a positive control; antibiotic-treated Ae. albopictus exposed to 1.0E+08 CFU/ml EBK constituted the experimental group. The E. anophelis aegypti-treated group demonstrated increased blood-feeding rates (40/75) relative to the conventionally reared individuals (19/75) (Chi-square test 12.3207; P = 0.0004) (Additional file 7). In spite of the differences in blood meal intake, E. anophelis aegypti-treated individuals demonstrated a 14% lower average ZIKV infection rate than the conventionally reared Ae. albopictus, which was 37% (Fisher’s exact test, P = 0.0003; Fig. 4d). There was no difference observed in the dissemination rates [conventionally reared (5%); E. anophelis aegypti-treated (7%) (Fisher’s exact test, P = 0.7673)], and no mosquito transmitted the virus at the measured time point (Additional file 8). Despite the fact that feeding on a ZIKV infectious blood meal was associated with a higher percentage of E. anophelis aegypti relative to other taxa (Fig. 2a), the fold change in expression of the Elizabethkingia rpoB gene from ZIKV-infected Ae. albopictus (N = 22) correlated negatively with ZIKV absolute quantities (Spearman r = −0.3985; P < 0.05); (Fig. 4e). ZIKV viral loads ranged from 1.00E+03 to 1.00E+11 ZIKV RNA copies/ml, while the log fold change of E. anophelis aegypti ranged from 0.016 to 158.49 CFU/ml (Additional file 9).
Infection of Vero cells with E. anophelis results in a broad-spectrum reduction of viral load
To clarify the potential for a broad-spectrum antiviral effect of E. anophelis aegypti, we performed in vitro assays on both Ae. albopictus and mammalian cells. A broad-spectrum E. anophelis aegypti antiviral effect was measured by significantly lower viral loads associated with Vero cells sequentially inoculated with E. anophelis aegypti/ZIKV (ANOVA test, P < 0.0001; Fig. 5a; Additional file 10), E. anophelis aegypti/CHIKV (ANOVA test, P < 0.0001) (Fig. 5a); (Additional file 10), and E. anophelis aegypti/DENV (ANOVA test, P < 0.0001; Fig. 5a; Additional file 10).
This antiviral property was specific to E. anophelis aegypti. We demonstrated no difference in viral loads when E. anophelis aegypti was replaced by E. coli, E. coli/ZIKV (ANOVA test, P = 0.9867), E. coli/CHIKV (ANOVA test, P > 0.9999), or E. coli/DENV (ANOVA test, P > 0.9999; Fig. 5b; Additional file 11). It is noteworthy that the antiviral effect was only measured when E. anophelis and the virus were inoculated sequentially, but not when co-inoculated with E. anophelis aegypti/ZIKV (ANOVA test, P = 0.6433) or E. anophelis aegypti/ZIKV (ANOVA test, P = 0.9987; Fig. 5c; Additional file 12). When both E. anophelis aegypti and E. coli were used to inoculate C636 cells, we observed excessive bacteria growth (E. anophelis aegypti and E. coli) after 24 hpi, and ZIKV replication was attenuated in C636 cells co-inoculated with E. anophelis aegypti/ZIKV (ANOVA test, P < 0.0001) as well as E. coli/ZIKV (ANOVA test, P < 0.001; Fig. 6a; Additional file 13). Viability assays demonstrated 90% viability of E. anophelis aegypti-infected cells at 2 dpi, suggesting direct cell killing was not responsible for decreased viral replication (Additional file 14). We performed an in vitro assay of E. anophelis aegypti/ZIKV on the Ae. albopictus-derived U4.4 cell line and measured a significant reduction in ZIKV titers in all E. anophelis aegypti/ZIKV-infected U4.4 cells relative to the control (ANOVA; Kruskal–Wallis test, P < 0.01; Fig. 6b; Additional file 15).
Recently, ZIKV has become a focus of intense research due to its unprecedented global spread, [4, 40,41,42,43,44] leading to outbreaks in different parts of the world. An upsurge in insecticide resistance, as well as climate change, has resulted in expansion of the geographical range of vector territories , which suggests an urgent need for alternative approaches to control mosquitoes and the pathogens they transmit.
Due to their potential as candidates for development of therapeutics and transmission-blocking agents, investigators have studied the interaction of gut microbial taxa and arboviruses of medical importance [9,10,11,12,13,14,15,16,17]. Aedes midgut bacterial species associated with ZIKV are poorly understood. The aim of this study was to assess the impact of increased temperature on the Ae. albopictus midgut microbiome and identify bacterial species that are associated with altered ZIKV infectivity and Ae. albopictus competence.
Reduction in abundance of midgut E. anophelis due to temperature increase
Given that the global increase in temperature is projected to accelerate [45, 46], a comprehensive assessment of the potential effects of climatic change on vectorial capacity requires an assessment of the influence of temperature on the mosquito microbiome. Studies have shown that temperature is an important component in shaping the microbial communities of organisms [47,48,49,50,51]. In an effort to understand how climate change might ultimately alter microbial taxa that influence vectorial capacity for ZIKV, we measured the effect of increases in temperature on the midgut microbial community. Specifically, we observed a reduction in levels of E. anophelis albopictus, which our data suggest is associated with decreased levels of ZIKV and competence. These data suggest that future increases in temperatures could alter microbial profiles and therefore further exacerbate documented increases in transmission resulting directly resulting from increases in viral replication at higher temperatures.
Impact of ZIKV infection and blood meal on Ae. albopictus gut microbiota diversity
A significant reduction in Ae. albopictus midgut microbial diversity resulting from ZIKV infection and/or blood meal digestion was evident, yet the effect on individual genera was highly variable. Our results are consistent with the findings of a previous study assessing the total gut microbiota diversity of Aedes japonicus and Aedes triseriatus following infection with La Crosse orthobunyavirus, which revealed higher midgut bacterial richness among the midguts obtained from newly emerged individuals versus lower richness among midguts exposed to viral infection or a non-infectious blood meal . In addition, our study corroborates the finding of a study aimed at assessing the impact of West Nile virus infection on the bacterial diversity of Culex pipiens, which revealed higher bacterial diversity at the genus level among the unexposed individuals relative to both exposed infected and the exposed uninfected individuals .
Mosquitoes exposed to non-infectious blood meals demonstrated substantial increases in the relative abundance of E. anophelis albopictus as compared to unfed individuals, clearly indicating that the gut environment is more amenable to this species’ growth following blood meal digestion. This alteration in microbial composition could be beneficial to Ae. albopictus. Specifically, E. anophelis has been shown to produce OxyR regulon and antioxidants that, in turn, could protect the mosquitoes from the oxidative stress associated with blood digestion. Elizabethkingia anophelis has also been shown to possess hemolytic activity and may therefore help facilitate erythrocyte digestion in the mosquito midgut .
Elizabethkingia anophelis reduction of viral loads
Our findings demonstrate that the E. anophelis aegypti isolated from Aedes midguts is associated with broad-spectrum reduction of viral load of in vitro ZIKV, DENV, and CHIKV, as well as a reduction in ZIKV infection rates, when fed as a supplemental diet to Ae. albopictus. We also measured a negative correlation between ZIKV and native E. anophelis albopictus in the midguts of Ae. albopictus. Further, our findings demonstrate a cell-independent effect such that E. anophelis aegypti maintains its antiviral activity with both intact and deficient RNA interference responses. Past studies show that E. anophelis attenuates Plasmodium parasite development . Taken together, these results suggest that E. anophelis has the capacity to directly or indirectly elicit a broad antipathogen phenotype. Further studies delineating the potential antiviral mechanism of E. anophelis are needed, as this could lead to the development of novel control measures.
We observed higher infection rates in the conventionally reared Ae. albopictus relative to the E. anophelis aegypti-infected Ae. albopictus group. However, this did not translate to increased virus dissemination, which corroborates the results of other studies [55,56,57,58,59]. While these data suggest that E. anophelis may primarily act to block viral infection, both the in vitro data and negative correlation of native E. anophelis with viral load suggest that increasing doses could further inhibit viral replication and dissemination. An antiviral effect was only measured when E. anophelis aegypti and the virus were inoculated sequentially. While the mechanistic basis is unknown, this suggests that establishment and/or replication of E. anophelis is required to effectively perturb viral infection and replication.
To the best of our knowledge, our study provides the first description of broad-spectrum antiviral association of E. anophelis. We demonstrate both a negative correlation of E. anophelis and ZIKV in Ae. albopictus and an association with reduced viral loads of ZIKV, DENV, or CHIKV in vitro. Further, E. anophelis was found to be temperature-sensitive, with decreased levels identified at higher temperatures. Together, these findings have significant implications for our understanding of the interactions of bacterial and viral agents in mosquitoes, not to mention how climate change could influence these interactions. Identifying the mechanistic basis for these associations could improve our understanding of disease epidemiology and ultimately inform the development of novel vector control measures.
Availability of data and materials
Datasets generated during the current study are available in GenBank, as well as included in this manuscript.
- ZIKV :
- EBK :
- CHIKV :
New York state
New York State Department of Health
- MD :
Infectious blood meal
- NIBM :
Non-infectious blood meal
Operational taxonomic units
Multiplicity of infection
Molecular Evolutionary Genetic Analysis
Faye O, Freire CCM, Iamarino A, Faye O, de Oliveira JVC, Diallo M, Zanotto PMA, Sall AA. Molecular evolution of Zika Virus during Its emergence in the 20th Century. PLoS Negl Trop Dis. 2014;8:1–10.
Cook S, Holmes EC. A multigene analysis of the phylogenetic relationships among the flaviviruses (Family: Flaviviridae) and the evolution of vector transmission. Arch Virol. 2006;151:309–25.
Chambers T. Flavivirus genome organization, expression, and replication. Annu Rev Microbiol. 1990;44:649–88.
Lanciotti RS, Kosoy OL, Laven JJ, Velez JO, Lambert AJ, Johnson AJ, Stanfield SM, Duffy MR. Genetic and serologic properties of Zika virus associated with an epidemic, Yap State, Micronesia, 2007. Emerg Infect Dis. 2008;14:1232–9.
Petersen LR, Jamieson DJ, Powers AM, Honein MA. Zika Virus. N Engl J Med. 2016;374:1552–63.
Gulland A. Zika virus is a global public health emergency, declares WHO. BMJ. 2016;352:657.
Kazmi SS, Ali W, Bibi N, Nouroz F. A review on Zika virus outbreak, epidemiology, transmission and infection dynamics. J Biol Res. 2020;27:1–11.
Wells MB, Andrew J. Anopheles salivary gland architecture shapes Plasmodium sporozoite availability for transmission. Am Soc Microbiol. 2019;10:1–17.
Ramirez JL, Short SM, Bahia AC, Saraiva RG, Dong Y, Kang S, Tripathi A, Mlambo G, Dimopoulos G. Chromobacterium Csp_P reduces malaria and dengue infection in vector mosquitoes and has entomopathogenic and in vitro anti-pathogen activities. PLoS Pathog. 2014;10:e1004398.
Ramirez JL, Souza-Neto J, Cosme RT, Rovira J, Ortiz A, Pascale JM, Dimopoulos G. Reciprocal tripartite interactions between the Aedes aegypti midgut microbiota, innate immune system and dengue virus influences vector competence. PLoS Negl Trop Dis. 2012;6:1–11.
Anglero-Rodriguez Y, Talyuli A, Blumberg B, Kang S, Demby C, Shields A, Carlson J, Jupatanakul N, Dimopoulos G. An Aedes aegypti -associated fungus increases susceptibility to dengue virus by modulating gut trypsin activity. Elife. 2017;6:1–20.
Wu P, Sun P, Nie K, Zhu Y, Shi M, Xiao C, Liu H, Liu Q, Zhao T, Chen X, Zhou H, Wang P, Cheng G. A gut commensal bacterium promotes mosquito permissiveness to arboviruses. Cell Host Microbe. 2019;25:101-112.e5.
Lindsey ARI, Bhattacharya T, Newton ILG, Hardy RW. Conflict in the intracellular lives of endosymbionts and viruses: a mechanistic look at Wolbachia-mediated pathogen-blocking. Viruses. 2018;10:141.
Asad S, Hussain M, Hugo L, Osei-Amo S, Zhang G, Watterson D, Asgari S. Suppression of the pelo protein by Wolbachia and its effect on dengue virus in Aedes aegypti. PLoS Negl Trop Dis. 2018;12:1–20.
Moretti R, Yen PS, Houé V, Lampazzi E, Desiderio A, Failloux AB, Calvitti M. Combining Wolbachia-induced sterility and virus protection to fight Aedes albopictus-borne viruses. PLoS Negl Trop Dis. 2018;12:1–20.
Mohanty I, Rath A, Swain SP, Pradhan N, Hazra RK. Wolbachia population in vectors and non-vectors: a sustainable approach towards dengue control. Curr Microbiol. 2019;76:133–43.
McLean BJ, Dainty KR, Flores HA, O’Neill SL. Differential suppression of persistent insect specific viruses in trans-infected wMel and wMelPop-CLA Aedes-derived mosquito lines. Virology. 2018;2019(527):141–5.
Chen S, Bagdasarian M, Walker ED. Elizabethkingia anophelis: molecular manipulation and interactions with mosquito hosts. Appl Environ Microbiol. 2015;81:2233–43.
Perrin A, Larsonneur E, Nicholson AC, Edwards DJ, Gundlach KM, Whitney AM, Gulvik CA, Bell ME, Rendueles O, Cury J, Hugon P, Clermont D, Enouf V, Loparev V, Juieng P, Monson T, Warshauer D, Elbadawi LI, Walters MS, Crist MB, Noble-wang J, Borlaug G, Rocha EPC, Criscuolo A, Touchon M, Davis JP, Holt KE, Mcquiston JR, Brisse S. Evolutionary dynamics and genomic features of the Elizabethkingia anophelis 2015 to 2016 Winconsin outbreak strain. Nat Commun. 2017;8(May):1–12.
Terenius O, Lindh JM, Eriksson-Gonzales K, Bussière L, Laugen AT, Bergquist H, Titanji K, Faye I. Midgut bacterial dynamics in Aedes aegypti. FEMS Microbiol Ecol. 2012;80:556–65.
Bahia AC, Dong Y, Blumberg BJ, Mlambo G, Tripathi A, Chandra R, Dimopoulos G. Exploring Anopheles gut bacteria for Plasmodium blocking activity. Environ Microbiol. 2015;16:2980–94.
Akhouayri IG, Habtewold T, Christophides GK. Melanotic pathology and vertical transmission of the gut commensal Elizabethkingia meningoseptica in the major malaria vector Anopheles gambiae. PLoS ONE. 2013;8:1–12.
Boissière A, Tchioffo MT, Bachar D, Abate L, Marie A, Nsango SE, Shahbazkia HR, Awono-Ambene PH, Levashina EA, Christen R, Morlais I. Midgut microbiota of the malaria mosquito vector Anopheles gambiae and interactions with Plasmodium falciparum infection. PLoS Pathog. 2012;8:1–12.
Wang Y, Gilbreath TM, Kukutla P, Yan G, Xu J. Dynamic gut microbiome across life history of the malaria mosquito Anopheles gambiae in Kenya. PLoS ONE. 2011;6:1–9.
Ngwa C, Glockner V, Abdelmohsen U, Scheuermayer M, Fischer R, Hentschel U, Pradel G. 16S rRNA gene-based identification of Elizabethkingia meningoseptica (Flavobacteriales:Flavobacteriaceae) as a dominant midgut bacterium of the Asian malaria vector Anopheles stephensi (Diptera: Culicidae) with antimicrobial activities. J Med Entomol. 2013;50:404–14.
Villegas LEM, Campolina TB, Barnabe NR, Orfano AS, Chaves BA, Norris DE, Pimenta PFP, Secundino NFC. Zika virus infection modulates the bacterial diversity associated with Aedes aegypti as revealed by metagenomic analysis. PLoS One. 2018;13:e0190352.
Ciota A, Bialosuknia S, Zink S, Brecher M, Ehrbar D, Morrissette M, Kramer LD. Effects of Zika virus strain and Aedes mosquito species on vector competence. Emerg Infect Dis. 2017;23:1110–7.
IPCC: Third Assessment Report. Cambridge, United Kingdom; 2007.
Dhariwal A, Chong J, Habib S, King IL, Agellon LB, Xia J. MicrobiomeAnalyst: a web-based tool for comprehensive statistical, visual and meta-analysis of microbiome data. Nucleic Acids Res. 2017;45:W180–8.
Chao A. Nonparametric estimation the number of classes in a population. Ann Math Stat. 1984;20:265–70.
Bray RJ, Curtis JT. An ordination of the upland forest communities of Southern Wisconsin. Ecol Monogr. 1957;27:325–49.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.
Liaw A, Wiener M. Classification and regression by randomForest. R news. 2002;2(December):18–22.
Kämpfer P, Matthews H, Glaeser SP, Martin K, Lodders N, Faye I. Elizabethkingia anophelis sp. nov., isolated from the midgut of the mosquito Anopheles gambiae. Int J Syst Evol Microbiol. 2011;61:2670–5.
Kämpfer P, Busse HJ, McInroy JA, Glaeser SP. Elizabethkingia endophytica sp. nov. isolated from Zea mays and emended description of Elizabethkingia anophelis Kämpfer et al. 2011. Int J Syst Evol Microbiol. 2015;65:2187–93.
Holmes B, Steigerwalt AG, Nicholson AC. DNA-DNA hybridization study of strains of Chryseobacterium, Elizabethkingia and Empedobacter and of other usually indole-producing non-fermenters of CDC groups IIc, IIe, IIh and IIi, mostly from human clinical sources, and proposals of Chryseobacterium be. Int J Syst Evol Microbiol. 2013;63(12):4639–62.
Li Y, Kawamura Y, Fujiwara N, Naka T, Hongsheng L, Huang X, Kobayashi K, Ezaki T. Chryseobacterium miricola sp. nov., A novel species isolated from condensation water of space station Mir. Syst Appl Microbiol. 2003;26:523–8.
Nicholson AC, Gulvik CA, Whitney AM, Humrighouse BW, Graziano J, Emery B, Bell M, Loparev V, Juieng P, Gartin J, Bizet C, Clermont D, Criscuolo A, Brisse S, McQuiston JR. Revisiting the taxonomy of the genus Elizabethkingia using whole-genome sequencing, optical mapping, and MALDI-TOF, along with proposal of three novel Elizabethkingia species: Elizabethkingia bruuniana sp. nov. Elizabethkingia ursingii sp. nov., and Eliz. Antonie van Leeuwenhoek. Int J Gen Mol Microbiol. 2018;111:55–72.
King E. Studies on a group of previously unclassified bacteria associated with meningitis in infants. Am J Clin Pathol. 1959;31:241–7.
Hayes EB. Zika virus outside Africa. Emerg Infect Dis. 2009;15:1347–50.
Hennessey M, Fischer M, Staples JE. Zika virus spreads to new areas—region of the Americas, May 2015–January 2016. Morb Mortal Wkly Rep. 2016;65:55–8.
Zanluca C, Campos V, De MA, Luiza A, Mosimann P, Igor G, Nunes C, Luz K. First report of autochthonous transmission of Zika virus in Brazil. Mem Inst Oswaldo Cruz. 2015;110(June):569–72.
Musso D, Nilles E, Cao-Lormeau V. Rapid spread of emerging Zika virus in the Pacific area. Clin Microbiol Infect. 2014;20:1–2.
Duffy MR, Chen T-H, Hancock WT, Powers AM, Kool JL, Lanciotti RS, Pretrick M, Marfel M, Holzbauer S, Dubray C, Guillaumot L, Griggs A, Bel M, Lambert AJ, Laven J, Kosoy O, Panella A, Biggerstaff BJ, Fischer M, Hayes EB. Zika Virus outbreak on Yap Island, Federated States of Micronesia. N Engl J Med. 2009;360:2536–43.
Vose R, Easterling K, Kunkel A, LeGrande, Wehner M. Temperature changes in the United States. In Clim Sci Spec Rep Fourth Natl Clim Assesment. Edited by Wuebbles D, Fahey D, Hibbard K, Dokken D, Stewart B, Maycock TUS. Global Change Research Program, Washington; 2017:185–206.
Romero-Lankao P, Smith J, Davidson D, Diffenbaugh N, Kinney P, Kirshen P, Kovacs P, Ruiz L. Climate change 2014: impacts, adaptation and vulnerability. In: Barros V, Field D, Dokken M, Mastrandrea K, Mach T, Bilir M, Chatterjee M, Ebi K, Estrada Y, Genova R, Girma E, Kissel A, Levy S, MacCracken S, Mastrandrea P, White L, editors. Reg Asp Contrib Work Gr II to Fifth Assesment Rep Intergov Panel Clim Chang. New York: Cambridge University Press; 2014. p. 1439–98.
Fuhrman JA, Steele JA, Hewson I, Schwalbach MS, Brown MV, Green JL, Brown JH. A latitudinal diversity gradient in planktonic marine bacteria. Proc Natl Acad Sci. 2008;105:7774–8.
Lokmer A, Mathias Wegner K. Hemolymph microbiome of Pacific oysters in response to temperature, temperature stress and infection. ISME J. 2015;9:670–82.
Ross PA, Wiwatanaratanabutr I, Axford JK, White VL, Endersby-harshman NM, Hoffmann AA. Wolbachia infections in Aedes aegypti differ markedly in their response to cyclical heat stress. PLoS Pathog. 2017;13:1–17.
Shan H, Deng W, Luan J, Zhang M-J, Zhen Z, Liu S-S, Liu Y-Q. Thermal sensitivity of bacteriocytes constrains the persistence of intracellular bacteria in whitefly symbiosis under heat stress. Environ Microbiol. 2017;9:706–16.
Hussain M, Akutse KS, Ravindran K, Lin Y, Bamisile BS, Qasim M, Dash CK, Wang L. Effects of different temperature regimes on survival of Diaphorina citri and its endosymbiotic bacterial communities. Environ Microbiol. 2017;19:3439–49.
Muturi EJ, Bara JJ, Rooney AP, Hansen AK. Midgut fungal and bacterial microbiota of Aedes triseriatus and Aedes japonicus shift in response to La Crosse virus infection. Mol Ecol. 2016;25:4075–90.
Zink SD, van Slyke GA, Palumbo MJ, Kramer LD, Ciota AT. Exposure to west nile virus increases bacterial diversity and immune gene expression in culex pipiens. Viruses. 2015;7:5619–31.
Kukutla P, Lindberg BG, Pei D, Rayl M, Wanqin Y, Steritz M, Faye I, Xu J. Insights from the genome annotation of Elizabethkingia anophelis from the malaria vector Anopheles gambiae. PLoS ONE. 2014;9:e97715.
Diagne CT, Diallo D, Faye O, Ba Y, Faye O, Gaye A, Dia I, Faye O, Weaver SC, Sall AA, Diallo M. Potential of selected Senegalese Aedes spp mosquitoes (Diptera: Culicidae) to transmit Zika virus. BMC Infect Dis. 2015;15:1–7.
Chouin-Carneiro T, Vega-Rua A, Vazeille M, Yebakima A, Girod R, Goindin D, Myrielle DR, Lourenco de Oliveira R, Failloux A. Differential susceptibilities of Aedes aegypti and Aedes albopictus from the Americas to Zika Virus. PLoS Negl Trop Dis. 2016;10:1–11.
Diagne CT, Faye O, Guerbois M, Knight R, Diallo D, Faye O, Ba Y, Dia I, Faye O, Weaver SC, Sall AA, Diallo M. Vector competence of Aedes aegypti and Aedes vittatus (Diptera:Culicidae) from Senegal and Cape Verde Archipelago for West African lineages of Chikungunya Virus. Am J Trop Med Hyg. 2014;91:635–41.
Calvez E, Guillaumot L, Girault D, Richard V, O’Connor O, Paoaafaite T, Teurlai M, Pocquet N, Cao-Lormeau V-M, Dupont-Rouzeyrol M. Dengue-1 virus and vector competence of Aedes aegypti (Diptera:Culicidae) populations from New Caledonia. Parasit Vectors. 2017;10:1–8.
Gaye A, Wang E, Vasilakis N, Guzman H, Diallo D, Talla C, Ba Y, Dia I, Weaver SC, Diallo M. Potential for sylvatic and urban Aedes mosquitoes from Senegal to transmit the new emerging dengue serotypes 1, 3 and 4 in West Africa. PLoS Negl Trop Dis. 2019;13(2):e0007043.
Naze F, Le Roux K, Schuffenecker I, Zeller H, Staikowsky F, Grivard P, Michault A, Laurent P. Simultaneous detection and quantitation of Chikungunya, Dengue and West Nile viruses by multiplex RT-PCR assays and Dengue virus typing using High Resolution Melting. J Virol Methods. 2009;162:1–7.
Laurent P, Le Roux K, Grivard P, Bertil G, Naze F, Picard M, Staikowsky F, Barau G, Schuffenecker I, Michault A. Development of a sensitive real-time reverse transcriptase PCR assay with an internal control to detect and quantify chikungunya virus. Clin Chem. 2007;53:1408–14.
Callahan JD, Wu SL, Dion-schultz A, Mangold BE, Peruski LF, Watts DM, Porter KR, Murphy GR, Suharyono W, King C, Hayes CG. Development and evaluation of serotype- and group-specific fluorogenic reverse transcriptase PCR (TaqMan) assays for dengue virus. J Clin Microbiol. 2001;39:4119–24.
We would like to extend our gratitude to Illia Rochlin of the Suffolk County (NY) Health Department for providing us with Ae. albopictus mosquitoes and the Arbovirus Laboratory insectary crew for mosquito maintenance. We also thank the Wadsworth Center Applied Genomics Core for sequencing and the Wadsworth Center Tissue and Media Facility for providing cells and media. We would like to extend our unreserved gratitude to Geoffrey Pollock Texas Tech University for editing the manuscript. Finally, we would like to express our gratitude to Steven Blaire at the Connecticut Agricultural Experiment Station for providing us with the E. anophelis strain used in this study.
This publication was supported by the Cooperative Agreement Number U01CK000509 funded by the Centers for Disease Control and Prevention. Its content is solely the responsibility of the authors and do not necessarily represent the official views of the Centers for Disease Control and Prevention or the Department of Health and Human Services. These studies were also partially funded by the National Institutes of Health awards R21AI131683 and R21 AI138344.
Ethics approval and consent to participate
Consent for publication
The authors declare they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
A standard curve showing R2, PCR efficiency and slope used to extrapolate viral load.
A box plot file showing the impact of temperature increase on individual midgut diversity.
A biomarker analysis cluster file bearing patterns of relative abundance of taxa species in response to temperature variations.
A box plot file showing the impact of blood meal on the midgut taxa.
Excel file showing the dissemination of ZIKV at 7 dpi by Ae. albopictus reared in the L and M temperature regimens.
Excel file showing the transmission of ZIKV at 7 dpi by Ae. albopictus reared in the L and M temperature regimens.
Excel file of the blood feeding rates of Ae. albopictus supplementarily fed with E. anophelis aegypti.
Excel file listing the ZIKV infection rates of Ae. albopictus supplementarily fed with E. anophelis aegypti.
Excel file of the Ae. albopictus ZIKV/E. anophelis aegypti correlation analysis.
Excel file of virus titer of sequential inoculation of Vero cells with ZIKV, CHIKV or DENV and E. anophelis aegypti.
Excel file of virus titer of sequential inoculation of Vero cells with ZIKV, CHIKV or DENV and E. coli.
Excel file of virus titer of co-inoculation of Vero cells with ZIKV and E. anophelis aegypti or E. coli.
An excel sheet of ZIKV titer of sequential inoculation of C636 cells with ZIKV and E. anophelis aegypti or E. coli.
Excel file of Vero cell viability at 24 hpi with E. anophelis aegypti.
Excel file of ZIKV titer following in-vitro assay of E. anophelis aegypti and ZIKV in U4.4 cells.
About this article
Cite this article
Onyango, M.G., Lange, R., Bialosuknia, S. et al. Zika virus and temperature modulate Elizabethkingia anophelis in Aedes albopictus. Parasites Vectors 14, 573 (2021). https://doi.org/10.1186/s13071-021-05069-7
- Zika virus
- Aedes albopictus