key: cord-0052047-6kr6cvi8
authors: Juma, Elijah O.; Kim, Chang-Hyun; Dunlap, Christopher; Allan, Brian F.; Stone, Chris M.
title: Culex pipiens and Culex restuans egg rafts harbor diverse bacterial communities compared to their midgut tissues
date: 2020-10-27
journal: Parasit Vectors
DOI: 10.1186/s13071-020-04408-4
sha: 96962e22c4ed051b3480005ad3f2e87075a02ca3
doc_id: 52047
cord_uid: 6kr6cvi8

BACKGROUND: The bacterial communities associated with mosquito eggs are an essential component of the mosquito microbiota, yet there are few studies characterizing and comparing the microbiota of mosquito eggs to other host tissues. METHODS: We sampled gravid female Culex pipiens L. and Culex restuans Theobald from the field, allowed them to oviposit in the laboratory, and characterized the bacterial communities associated with their egg rafts and midguts for comparison through MiSeq sequencing of the 16S rRNA gene. RESULTS: Bacterial richness was higher in egg rafts than in midguts for both species, and higher in Cx pipiens than Cx. restuans. The midgut samples of Cx. pipiens and Cx. restuans were dominated by Providencia. Culex pipiens and Cx. restuans egg rafts samples were dominated by Ralstonia and Novosphingobium, respectively. NMDS ordination based on Bray-Curtis distance matrix revealed that egg-raft samples, or midgut tissues harbored similar bacterial communities regardless of the mosquito species. Within each mosquito species, there was a distinct clustering of bacterial communities between egg raft and midgut tissues. CONCLUSION: These findings expand the list of described bacterial communities associated with Cx. pipiens and Cx. restuans and the additional characterization of the egg raft bacterial communities facilitates comparative analysis of mosquito host tissues, providing a basis for future studies seeking to understand any functional role of the bacterial communities in mosquito biology. [Image: see text]

Studies applying high throughput, culture-independent sequencing of the bacterial 16S rRNA gene have advanced understanding of the association between mosquitoes and their bacterial communities [1] [2] [3] [4] . The bulk of studies characterizing mosquito-associated bacterial communities have focused on the mosquito gut.

However, other mosquito organs or tissues, including the ovaries, the male reproductive system, the salivary glands, and eggs, are also known to harbor bacterial communities that may play essential roles in mosquito biology [5] [6] [7] [8] [9] [10] [11] . The research focus on the mosquito gut, especially in the adult stage, is underpinned by the understanding that the mosquito midgut environment typically is the first barrier that mosquito-borne pathogens must overcome to develop successfully within the mosquito host and be transmitted to the next susceptible host [12] [13] [14] [15] .

Mosquito egg bacterial communities likely also play important roles in mosquito ecology. Studies with Aedes aegypti and Ae. triseriatus have shown that bacterial density on the egg surface mediates egg hatch rates and timeto-hatch, while mosquito eggs subjected to heavy larval grazing in high larval density habitats can exhibit delayed time-to-hatch [16] [17] [18] . The presence of aerobic microorganisms on the egg surface and in the larval environment has been associated with reduced oxygen tension, providing the stimulus for mosquito egg hatching [19, 20] . To what extent the bacterial communities on the egg surface drive the hatching effect relative to the microbes in the water column remains unclear.

Previous studies characterizing and comparing the bacterial communities of mosquito eggs with those of their other host tissues (e.g. midguts) have been conducted with well-known Afro-tropical or Asian vectors, including Anopheles or Aedes species, but not important North American vector species such as Culex pipiens L. or Culex restuans Theobald [5] [6] [7] . Culex pipiens is an introduced European species that arrived in North America in the early 16th century through trade and has been naturalized in the United States north of 39° latitude [21] [22] [23] . It serves as both an amplifying and bridge vector for West Nile virus (WNV) and St Louis encephalitis due to its preference for feeding on birds [24] . Culex restuans, native to North America, is also an important vector of WNV and distributed in the northeast and Great Lakes regions of the USA [22, 24] . The two species are common in urban, residential neighborhoods and woodlots and are ecologically similar and spatially-overlapping throughout much of their ranges, primarily utilizing birds as blood-meal sources and artificial container habitats for juvenile development [25] [26] [27] [28] . Comparative studies with Cx. pipiens and Cx. restuans are necessary to expand the known library of bacterial communities associated with mosquito host tissues and to facilitate further studies on the role of these bacterial communities in mosquito vector biology and ecology and potentially benefit mosquitoborne disease control.

We used Illumina MiSeq sequencing of the V3-V4 hypervariable regions of the 16S rRNA gene to characterize the bacterial communities associated with egg rafts and midguts of Cx. pipiens and Cx. restuans, to gain insights into bacterial community structure and diversity, and to observe how they compare between the two host tissues. We tested the hypothesis that, within each mosquito species, egg-raft-and the midgut samples will harbor distinct bacterial communities given that the external egg surface and the mosquito midgut represent physiologically distinct environments, and thus are likely to support distinct bacterial communities. We also hypothesized that similar tissues across species (e.g. eggs or midguts for both Cx. pipiens and Cx. restuans) will harbor similar bacterial communities because similar tissues represent similar physiological environments. This study expands upon understanding of the bacterial communities associated with mosquito host tissues other than the mosquito gut and provides a basis for further studies focused on the role of mosquito egg-associated bacterial communities in mosquito biology and mosquito-borne disease control. 14 ha "virgin" deciduous upland forest. It is primarily composed of mature oak (Quercus spp.), ash (Fraxinus spp.), and maple (Acer saccharinum L.) forest with a high, closed canopy and fairly open understory. Sugar maple is becoming the dominant tree species. South Farms is a 8.15-ha woodland composed of low canopy trees consisting mainly of sugar maple, sycamore (Platanus occidentalis L.), and pine (Pinus spp.), with oak and patchy grass undergrowth in some sections [27] . Approximately two years prior to this study, the invasive Amur honeysuckle (Lonicera maackii), that was the dominant shrub in this woodlot was removed. Trelease Woods is a 28.80ha deciduous forest consisting mainly of mature oak, ash, hackberry (Celtis occidentalis) and maple species, with a high, closed canopy and moderately dense understory. The site includes two small seasonal ponds which provide suitable habitat for many aquatic and semi-aquatic invertebrates, including mosquitoes. At each sampling site, two CDC gravid traps baited with 3.8 l each of grass infusion [29] were deployed beginning on June 18, 2018 and sampling was conducted three times weekly up to July 20, 2018. Traps were placed in the evening just before dusk and the collection bags collected approximately 12 h later the next morning [30] . Individual gravid female Culex species mosquitoes were transferred separately to individual 270 ml paper cups to facilitate oviposition; each consisted of a 30 ml inner plastic oviposition cup filled to half capacity with distilled water. The mosquitoes were maintained in a walk-in environmental chamber at 27 ± 1°C and ~75 ± 5% relative humidity with a 16:8 (L:D) photoperiod. The paper cups were provisioned with cotton balls soaked in distilled water to provide additional humidity and a source of water for the gravid females. Monitoring for egg rafts was conducted every 12 h. Following oviposition, the egg raft and the parous female were separately preserved at -80 °C for future bacterial DNA extraction.

Egg rafts and mosquito samples were thawed and adult mosquito samples surface sterilized in 70% ethanol for 5 min, transferred to 3% bleach solution for 3 min, transferred again to 70% ethanol for 5 min, and then rinsed 3 times in sterile water and 4 times in Dulbeccoʼs phosphate-buffered saline (DPBS) (Thermo Fisher Scientific, Waltham, MA, USA) [4] . Each sample was dissected in a small drop of sterile DPBS using a dissecting stereomicroscope and the midguts were transferred to PowerSoil bead tubes. Similarly, egg rafts were individually transferred to PowerSoil bead tubes. Due to their hydrophobicity and ease of disintegration during handling, egg rafts were processed for bacterial DNA without surface sterilization. Samples were homogenized using Retsch MM 300 TissueLyser (Retsch, Haan, Germany) and genomic DNA was extracted using MoBio PowerSoil DNA Isolation Kit (MoBio Laboratories, Inc., CA, USA) according to the manufacturer's instructions. DNA was quantified using the Nanodrop 1000 (Thermo Fisher Scientific, Pittsburgh, PA, USA).

Sequencing was performed at the National Center for Agricultural Utilization Research, Peoria, IL. The V3-V4 hypervariable region of bacterial 16S rRNA gene was PCR-amplified using previously published universal primers 341f and 806r [31, 32] . The V3-V4 hypervariable region has been shown to have higher sensitivity in bacterial phylogenetic analysis compared to the rest of the hypervariable regions of the 16S rRNA gene [1] . The following primer set specific for the V3-V4 region of the 16S rRNA gene was used: forward (5'-CCT ACG GGN GGC WGC AG-3'); reverse (5'-GAC TAC HVG GGT ATC TAA TCC-3'). The primers were incorporated into fusion primers for dual indexing and incorporation of adapters prior to genome sequencing using Illumina MiSeq (Illumina Inc., San Diego, CA, USA) [33] . The V3-V4 hypervariable region of the bacterial 16S rRNA gene was PCR-amplified using the following primer set: forward (5'-CCT ACG GGN GGC WGC AG-3'); reverse (5'-GAC TAC HVG GGT ATC TAA TCC-3'). PCR was conducted in 25 µl reactions containing 12 A negative control sample made up of DNA extracted from molecular biology grade water was sequenced with the same protocol to allow detection of the contamination background. PCR amplicons were cleaned and normalized using a SequalPrep normalization plate (Thermo Fisher Scientific, Waltham, MA, USA). The pooled library was mixed with Phix control spike-in of 5% as a sequencing control. The samples were sequenced on Illumina MiSeq system with a MiSeq V3 2 × 300 bp sequencing kit. The demultiplexed reads were quality-trimmed to Q30 using CLC genomics workbench v12.0 (Qiagen Inc., Valencia, CA, USA). Read pairing, fixed-length trimming and OTU clustering were done using CLC Bio Microbial Genomics module (Qiagen Inc., Valencia, CA, USA) utilizing the reference sequences from the Greengenes ribosomal RNA gene database [34] . The operational taxonomic unit (OTU) assignment was done at 97% sequence similarity, which is considered adequate for bacterial identification to the genus level [35] .

A duplex real-time TaqMan PCR assay [36] was used for molecular identification of Cx. pipiens and Cx. restuans using primers and probes targeting the acetylcholinesterase gene (Ace2) adopted from [36] . Primers and probes for Cx. pipiens consisted of: CxPip-F1 (5'-GGT GGA AAC GCA TGA CCA GAT A-3'); CxPip-R1 (5'-TGC AAT AAA GAG GTG GCC ACG-3'); and probe (FAM/ AGC CAC GAA CAA CTA AAT CAT CAC AAG CAC AGC/3BHQ). Culex restuans primers and probes were as follows: CxRest-F1 (5'-ATC GGT CTG GCT TCC TTT CAG AT-3'); CxRest-R1 (5'-TTA GTC AAG TTA ACT GGC CTA CAT CCT A-3'), and the probe (JOE/ AGC AAA CTG GCC GTC GTC CAC CGA TAT AA AT/3BHQ_1). The target DNA used as a template was taken from DNA samples extracted from midgut samples for bacterial DNA analysis. Separate DNA samples were extracted from Cx. pipiens and Cx. restuans adults initially identified from larval stages and used as positive controls. Additionally, a reaction mixture consisting of Ae. albopictus DNA template minus reverse transcriptase was used as a negative control to rule out any chance of contamination in the PCR reaction. Each PCR sample was assayed in 25 µl reaction mixture consisting of 5 µl of the target DNA, 12.5 µl SensiFAST ™ Probe Hi-ROX Kit master mix (Bioline, Tauton, MA, USA), 1.25 µl each of Cx. pipiens forward and reverse primer; 0.625 µl Cx. pipiens probe (Fam); 1.25 µl each of Cx. restuans forward and reverse primer; 0.625 µl Cx. restuans probe (Fam) and 1.25 µl of nuclease-free water. Thermocycling was performed on an ABI 7300 HT sequence detection system (Applied Biosystems, Foster City, CA, USA) using the following reaction conditions: 95 °C for 5 min followed by 40 cycles of 95 °C for 30 s, 58 °C for 30 s, and 72 °C for 60 s [36] .

All analyses were conducted using R version 3.6.1 [37] within the Rstudio environment version 1.2.1335 [38] and PAST version 3.15 [39] . OTUs accounting for < 0.005% of the total number of sequences were removed prior to analysis to eliminate spurious OTUs [40] . The OTU sequence numbers varied markedly within samples (mean ± SE = 7621.28±547.17 per sample). Bacterial sequences were rarefied to an even depth of 1007 reads per sample to standardize the sampling coverage [41, 42] . From an initial sample size of 188 samples, 44 samples had < 1007 reads and were excluded from further analysis. To estimate sample coverage, rarefaction curves were fitted on the unrarefied data using the phyloseq package version 1.24.0 in R [41] [42] [43] . Alpha diversity metrics, including Shannon diversity index, observed species, and Chao1, were generated in QIIME 2 [44] . The means and 95% confidence intervals were calculated in R to test for significant differences in the alpha diversity indices between treatments. The Kruskal-Wallis test was used to test for differences in means between egg-raft samples and midgut samples and the Wilcoxon rank sum test with Bonferroni correction was performed pairwise to separate significant treatments. Beta-diversity measures were estimated using the Bray-Curtis dissimilarity index using the phyloseq package and non-metric multidimensional scaling (NMDS) ordination plots were generated to visualize the results. The non-parametric Analysis of Similarity (ANOSIM) test with Bonferroni adjustments was performed in PAST version 3.26 [39] to determine degree of dissimilarity in bacterial composition between treatment groups. Similarity percentage (SIMPER) analysis was performed in PAST to identify the bacterial species characterizing each treatment group. Venn diagrams were generated using the R package limma [45] version 3.40.2 to visualize OTUs that were shared between egg rafts and midgut samples of Culex pipiens and Culex restuans. Using the function chisq.test from the R package RVAidemMemoire version 0.9-74, differences in bacterial OTU per sample were tested and pairwise multiple comparison test with Bonferroni correction was applied.

Sequencing of the V3-V4 regions of the 16S rRNA gene from 188 samples (66 Cx. pipiens egg rafts; 66 Cx. pipiens midguts; 28 Cx. restuans egg rafts; and 28 Cx. restuans midguts) generated 1,432,800 raw sequences (mean ± SE = 7621.28 ± 547.17 per sample). After quality-filtering to remove chimeric sequences, other non-bacterial sequences, and bacterial OTUs constituting < 0.005% of the total sequences and rarefying the reads to an even depth of 1007 sequences per sample to standardize sampling effort, a total of 144 samples were retained ( Rarefaction analysis of the bacterial OTU samples revealed that the sequencing depth coverage sufficiently recovered most of the bacterial OTUs. Chao1 estimator revealed that up to 86.3 ± 0.07% (mean ± SE) of the bacterial OTUs were recovered. The highest bacterial OTU richness was reported in Cx. pipiens eggs, while Cx. restuans midgut had the lowest bacterial OTU richness (Additional file 1: Figure S1 ). Culex pipiens egg-raft samples had significantly higher observed and expected (Chao1) bacterial OTU richness compared to Cx. pipiens midgut samples, or Cx. restuans egg-raft and midgut samples. Culex restuans egg samples had significantly higher observed and expected (Chao1) bacterial OTU richness compared to Cx. restuans midgut samples (Observed OTUs: Kruskal-Wallis χ 2 = 92.72, df = 3, P < 0.0001; Chao1: Kruskal-Wallis χ 2 = 89.02, df = 3, P < 0.0001; Shannon index: Kruskal-Wallis χ 2 = 85.82, df = 3, P < 0.0001) ( Table 1) .

The 153 bacterial OTUs were classified into 7 phyla, 13 classes, 27 orders, 40 families, and 54 genera. The most dominant phyla were Proteobacteria (85.5%), consisting of Alphaproteobacteria (34.0%), Betaproteobacteria (18.4%), Gammaproteobacteria (33.0%) and Deltaproteobacteria (0.01%). Other phyla included Spirochaetes (9.3%), Bacteroidetes (3.3%) and Firmicutes (1.1%), and the rest were < 1% cumulatively (Fig. 1a) . Alpha-and Betaproteobacteria were dominant in egg-raft samples of either species, while Gammaproteobacteria was dominant in the midgut samples of both species. Spirochaetes was abundant in Cx. restuans midgut samples (Fig. 1a) . The top five most abundant families accounted for 68.7% of all sequences. They included Enterobacteriaceae (25.8%), Sphingomonadaceae (15.9%), Oxalobacteraceae (10.6%), Borreliaceae (9.3%), and Rickettsiaceae (7.1%). Enterobacteriaceae was dominant in the midgut samples of both species, Sphingomonadaceae in Cx. restuans egg samples, Oxalobacteraceae in Cx. pipiens egg samples, and Borreliaceae in Cx. restuans midgut samples (Fig. 1b) . At the genus level, the top five most abundant OTUs accounted for 58.1% of all sequences. They included Providencia (17.8%), Novosphingobium (13.6%), Ralstonia (10.3%), Spironema (9.3%) and Wolbachia (7.1%) (Fig. 1c) . Providencia was the dominant taxon in the midgut samples of both species, Novosphingobium in Cx. restuans egg samples, Ralstonia in Cx. pipiens egg samples, and Spironema in Cx. restuans (Fig. 1c) . Overall, 64 (41.8%) bacterial OTUs were shared between all sample type combinations of mosquito species and host tissue (Additional file 1: Figure S2a ). Seventy-two bacterial OTUs (61.5%) were shared between Cx. pipiens midgut samples and Cx. restuans midgut samples and 128 of the bacterial OTUs (91%) were shared between the egg samples of the two mosquito species. Ninety-three (62.8%) bacterial OTUs were shared between egg and midgut samples of Cx. pipiens, while 73 bacterial OTUs (51.4%) were shared between egg and midgut samples of Cx. restuans mosquitoes. Overall, there was higher bacterial OTU richness and diversity in egg rafts compared to midgut samples for both Cx. pipiens and Cx. restuans. One hundred and five bacterial OTUs were detected in Table 2 ). SIMPER analysis identified 9 bacterial OTUs that were responsible for 70% of the observed differences between groups, with Providencia (19.14%), Ralstonia (10.57%), Novosphingobium (10.34%), and Spironema (8.39%) constituting the largest variation (Additional file 1: Table S1 ). Providencia was the most dominant bacterial OTU in midgut samples of both Cx. pipiens and Cx. restuans. Ralstonia was the dominant bacterial OTU in Cx. pipiens egg-raft samples but was also present in high proportions in Cx. restuans egg-raft samples. Novosphingobium was dominant in egg-raft samples from both species. Spironema was the dominant bacterial taxon in Cx. restuans midgut samples (Fig. 1c ).

In this study, we characterized the composition and diversity of the bacterial communities associated with egg rafts and midgut samples of Cx. pipiens and Cx. restuans. The egg-raft samples in both species were more diverse compared to midgut samples, with Cx. restuans midgut samples supporting the lowest bacterial diversity. Bacterial communities clustered by mosquito host tissue, such that the egg rafts from Cx. pipiens and Cx. restuans shared substantially similar bacterial communities and so did their midgut samples. However, both species had significantly different bacterial communities between their egg-raft and midgut tissues.

The bacterial communities associated with mosquito eggs are mostly localized on the external surfaces of the eggs. Previous studies have not been able to isolate bacterial communities from within the egg cytoplasm [4, 6, 7, 17, 46] . The egg-raft samples from both species in our study had highly overlapping bacterial communities with up to 91% of the bacterial OTUs shared between them. We presume, based on existing evidence, that these bacterial communities were mostly localized on the exterior of the egg rafts, thus representing the natural bacterial communities supported by the exterior egg-raft surface environment. In this study, female mosquito oviposition took place in deionized water, which is deficient in microorganisms [47] , and the egg rafts were preserved at -80 °C within hours of oviposition. These measures limit the possibility that any significant level of bacterial colonization of the egg rafts may have occurred immediately following oviposition. We suspect that most of the bacterial communities colonizing the egg rafts were inherited maternally from the ovaries through egg-smearing or through a yet-to-be-described form of transovarial transmission method [6, 7] . However, this study did not characterize the bacterial communities associated with the ovarian tissues from either species to validate this possibility. Additional future studies characterizing the bacterial communities of mosquito ovaries in addition to the egg rafts and midgut tissues would shed more light on this question. This study did not assess the potential bacterial composition of the deionized water used as oviposition substrate by the gravid females. While many studies assume deionized water to be microbe-free, this requires further investigation in future studies of the microbial composition of mosquito tissues. The midgut samples of Cx. pipiens and Cx. restuans mostly shared similar bacterial community composition with over 61% of the bacterial OTUs shared between them. Previous studies with the adults of these two mosquito species have generated variable results, with one study showing that they harbor distinctly different bacterial communities [48] , while another study did not report unique clustering in the bacterial communities of the two species [49] . Both studies were conducted with non-blood-fed adult mosquito samples whose parity status was not assessed, a factor that might have contributed to the observed differences. We can partly attribute the overlap in the bacterial OTUs of the adult midgut samples of the two species to the shared environment resulting in colonization by similar bacterial communities. Mosquitoes from different species sampled from common habitats have been shown to share more similar bacterial communities, a possible indication of horizontal acquisition from the surrounding environment [8, 49, 50] . The two species also overlap in their larval habitat use, which potentially predisposes them to acquisition of similar bacterial communities that are then propagated transstadially all the way to adulthood [28] . Similarly, the blood-fed and gravid status of the females at the time of sampling may have contributed to similar internal gut oxidative environment subsequently supporting comparable bacterial communities [51, 52] . It may also be possible that the two mosquito species largely blood-fed on similar blood-meal sources, most likely birds, given their ecological and behavioral similarity, as well as spatial and seasonal overlap, although the two species also display some degree of seasonal separation [22, 25, 53] . Future experimental designs could incorporate blood-meal analysis and location data to refine other correlates of mosquito tissue bacterial composition and diversity.

The significant separation of the bacterial community composition between egg rafts and midguts within each mosquito species, is not surprising since the physiological environment in the mosquito gut is expected to differ substantially from that of the egg rafts, and thus is likely to facilitate colonization by very different consortia of bacterial communities. The disproportionate dominance of a few distinct bacterial taxa in the egg rafts compared to midgut samples for both species further indicates that egg and midgut environments were substantially different. However, studies comparing the bacterial communities of mosquito midguts and eggs are rare. One such study focusing on Aedes aegypti shows that mosquito guts share a significant proportion of their bacterial community composition with those of the eggs. Most of the taxa that have been reported to be shared between mosquito egg and midgut life stage are widespread in mosquito species and have been described in many other mosquito microbial studies [5] .

There was high bacterial richness in the egg-raft samples compared to midgut samples for both Cx. pipiens and Cx. restuans. The midgut samples were from gravid females that possibly experienced a sudden and sharp decline in bacterial diversity due to their prior bloodmeal diet. Wang et al. (2011) previously showed that mosquito gut bacterial diversity falls markedly in bloodfed compared to sugar-fed adults [51] . They hypothesized that this was mediated by the breakdown of the heme proteins from the blood-meal diet releasing reactive oxygen species that alter the physiological configuration of the gut environment, thus limiting the bacterial communities to those that can tolerate the high oxidative stress [51, 54] . However, more recent study by Muturi et al. [54] reported increased bacterial diversity in mosquitoes exposed to different blood-meal sources relative to sugar meal, indicating that blood-meal diet may produce variable effects on mosquito gut bacterial diversity, and that bacterial diversity is additionally determined by mosquito species, as well as the source of the blood meal [54] . The digestive process involving the movement of the blood bolus along the midgut endoperitrophic space may also have contributed to physical propulsion and excretion of a significant proportion of the midgut bacterial communities further reducing the midgut bacterial diversity [55] . The lower bacterial diversity in the guts of field-collected mosquitoes relative to their egg tissues may also be influenced by exposure to insecticides in the wild, or the insecticide resistance status of the female mosquitoes. The guts of insecticide-resistant mosquitoes may be enriched with bacterial communities with selective advantage over the pesticide constituents thus reducing their bacterial diversity [56, 57] . The gravid female adults used in this study were all field-collected. We did not determine their chronological (calendar ages) or biological age (number of gonotrophic cycles) as this was beyond the scope of this study. Future studies can explore how chronological age, biological age, and possibly insecticide resistance status of field-collected gravid mosquitoes may influence the bacterial communities of mosquito tissues such as the midguts and the egg stages.

The dominance of Ralstonia and Novosphingobium in egg-raft samples of both species point to their possible adaptation to colonizing the mosquito egg raft stages. It also could be related to their potential role in mosquito ecology, such as inducing egg hatch, but this requires further research. Literature is scanty on the isolation and characterization of Ralstonia or Novosphingobium from mosquito egg stages, but this may be attributed to the dearth of studies characterizing mosquito egg bacterial communities. However, Ralstonia and Novosphingobium have been isolated in mosquito midgut samples from several mosquito species including Aedes japonicus, Ae. aegypti, Ae. albopictus and Anopheles coluzzii. They have also been described from the natural environment, including soil and aquatic sources, indicating that mosquitoes may acquire them horizontally [58] [59] [60] [61] [62] [63] . Their isolation from egg rafts in this study shows that these taxa may be part of the common mosquito bacterial commensals shared between different mosquito life stages including the eggs. Additional studies characterizing the bacterial communities of Culex ovaries, in addition to egg rafts, would help shed light on microbial presence in the ovary tissues and whether they are passed maternally to the egg rafts through egg-smearing or other forms of transovarial transmission.

The midgut samples of Cx. pipiens and Cx. restuans were dominated by Providencia, while Wolbachia and Spironema were the second most abundant bacterial taxa in Cx. pipiens and Cx. restuans midgut samples, respectively. We attribute the dominance of Providencia in both species potentially to the blood-meal diet prior to oviposition, and potentially from an avian blood-meal source, given its previous isolation in birds [65] . This may be an indication that Providencia could tolerate the high oxidative stress and the enzymatic conditions in the midgut environment provided by the breakdown of heme proteins from the blood-meal diet [51, 54] . The genus Providencia has been isolated from several mosquito species, including Anopheles albumanus [56] , Aedes aegypti [59, 64] and Aedes vexans [49] . However, no existing study has isolated the bacterial genus Providencia in Cx. pipiens or Cx. restuans. Providencia is a genus consisting of gram-negative rods with peritrichous flagella belonging to the family Enterobacteriaceae. It is an enteric bacterial pathogen commonly isolated from human intestines. It has also been isolated from other organisms, including birds and pigs, pointing to its ubiquity in the natural environment. However, its role in mosquito biology has not been described. The ease of culturing Providencia in bacteriological media, its ready availability in the natural environment, and its abundance in Cx. pipiens in this study points to its potential suitability as a candidate for manipulation for paratransgenesis for mosquito vector management [65] . The bacterial genus Spironema has been characterized in Culiseta melanura [66] , Ae. aegypti un-infected with Wolbachia [67] , Cx. pipiens [68] , and was also dominant in Culex nigripalpus [69] . This taxon has been isolated from soil as well as river water samples, providing evidence of its potential horizontal acquisition in the Cx. restuans in this study [70] . High abundance in Cx. restuans has not been reported previously and opens an avenue to conduct further studies on its potential role in mosquito biology and disease transmission. The presence of Wolbachia in Cx. pipiens was expected as Cx. pipiens naturally harbor Wolbachia, a maternally inherited endosymbiont common in many arthropods, where they mediate several reproductive manipulations in their hosts [71] . Our study did not report any trace of Wolbachia in Cx. restuans egg-raft samples. However, it was reported in the Cx. restuans midgut samples at < 0.1%, providing further evidence of its recent detection in Cx. restuans elsewhere [72] . The presence of Wolbachia in Cx. restuans and Cx. pipiens has potential to alter the epidemiology of WNV infection especially in the regions of the US where WNV has become endemic. Transient infection of Cx. tarsalis with the wAlbB strain has been associated with enhanced replication of WNV in the host [73] . In other studies, transient infection of Ae. aegypti with WNV was associated with enhanced viral replication [14] . However, it is likely that the interaction between Wolbachia and WNV in the naturally infected Culex populations could manifest differently. This area requires further investigation to disentangle the role of Wolbachia in WNV incubation and subsequently the epidemiology of the disease in the endemic regions.

In conclusion, our study has shown that Cx. pipiens and Cx. restuans egg-raft samples are more diverse in their bacterial communities and the bacterial communities differ significantly between egg-raft and midgut tissues within each mosquito species. However, the bacterial communities of the egg-raft versus the midgut tissues of the two species are mostly similar in their community composition. Whereas previous studies with Cx. pipiens or Cx. restuans have prioritized characterizing the bacterial communities from adult midguts [48, 49, 74] , the additional characterization of the egg raft bacterial communities in this study fills an important gap in our understanding of the bacterial communities associated with the egg rafts and how they compare with those of the midguts. These findings open the way for further studies on their role in mosquito biology and ecology and their potential to be exploited for mosquito vector management such as strategies to interfere with egg hatching in the environment.

Supplementary information accompanies this paper at https ://doi. org/10.1186/s1307 1-020-04408 -4.

Additional file 1: Table S1 . SIMPER analysis of the major bacterial OTUs driving differences between sample treatments (mosquito species and life stage). The Bray-Curtis average dissimilarity between sample treatments was >1% for 18 bacterial taxa. Overall average dissimilarity. CXP.EG -Cx. pipiens egg raft samples; CXP.MG -Cx. pipiens midgut samples; CXR.EG -Cx. restuans egg raft samples; CXR.MG -Cx. restuans midgut samples. 

Sensitivity and correlation of hypervariable regions in 16S rRNA genes in phylogenetic analysis

The mosquito holobiont: fresh insight into mosquito-microbiota interactions. Microbiome

Mosquitoes host communities of bacteria that are essential for development but vary greatly between local habitats

Mosquitoes rely on their gut microbiota for development

Culture-dependent and culture-independent characterization of microorganisms associated with Aedes aegypti (Diptera: Culicidae) (L.) and dynamics of bacterial colonization in the midgut

Bacteria of the genus Asaia stably associate with Anopheles stephensi, an Asian malarial mosquito vector

Mosquito-bacteria symbiosis: The case of Anopheles gambiae and Asaia

Composition of Anopheles coluzzii and Anopheles gambiae microbiota from larval to adult stages

Dynamics of bacterial community composition in the malaria mosquito's epithelia

The reproductive tracts of two malaria vectors are populated by a core microbiome and by gender-and swarm-enriched microbial biomarkers

Salivary glands harbor more diverse microbial communities than gut in Anopheles culicifacies

Implication of the mosquito midgut microbiota in the defense against malaria parasites

Natural microbe-mediated refractoriness to Plasmodium infection in Anopheles gambiae

Effect of Wolbachia on replication of West Nile virus in a mosquito cell line and adult mosquitoes

Chromobacterium Csp_P reduces malaria and dengue infection in vector mosquitoes and has entomopathogenic and in vitro anti-pathogen activities

To hatch or not to hatch? Egg hatch response to larval density and to larval contact in a treehole mosquito

Erratic hatching in Aedes eggs: a new interpretation

The complex hatching response of Aedes eggs to larval density

Eggs of Aedes triseriatus and Ae hendersoni: a method to stimulate optimal hatch

The hatching stimulus for eggs of Aedes punctor (Diptera: Culicidae)

Population genetic structure of the Culex pipiens (Diptera: Culicidae) complex, vectors of West Nile virus, in five habitats

Interspecific competition between larval Culex restuans Theobald and Culex pipiens L. (Diptera: Culicidae) in Michigan

The colonization of temperate North America by mosquitoes and man

Modeling the impact of variable climatic factors on the crossover of Culex restauns and Culex pipiens (Diptera: Culicidae), vectors of West Nile virus in Illinois

Population dynamics of Culex restuans and Culex pipiens (Diptera: Culicidae) related to climatic factors in northwest Ohio

A standardized procedure for the quantitative surveillance of certain Culex mosquitoes by egg raft collection

Oviposition preferences of Culex pipiens and Culex restuans for infusion-baited traps

Habitat characteristics of larval mosquitoes in zoos of South Carolina, USA

Oviposition preferences of Culex restuans and Culex pipiens (Diptera: Culicidae) for selected infusions in oviposition traps and gravid traps

Comparative analysis of the circadian rhythm genes period and timeless in Culex pipiens Linnaeus, 1758 (Diptera, Culicidae)

Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample

Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA

Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms

Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB

Estimating bacterial diversity for ecological studies: methods, metrics, and assumptions

A real-time TaqMan polymerase chain reaction for the identification of Culex vectors of West Nile and Saint Louis encephalitis viruses in North America

R: A language and environment for statistical computing

RStudio: Integrated development for R

PAST: Paleontological statistics software package for education and data analysis

Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing

Counting the uncountable: statistical approaches to estimating microbial diversity

A unified mathematical framework for the measurement of richness and evenness within and among multiple communities

Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data

QIIME allows analysis of high-throughput community sequencing data

limma powers differential expression analyses for RNA-sequencing and microarray studies

Failure to isolate microorganisms from within mosquito eggs

Discovery and exploitation of a natural ecological trap for a mosquito disease vector

Culex pipiens and Culex restuans mosquitoes harbor distinct microbiota dominated by few bacterial taxa

Comparative analysis of gut microbiota of mosquito communities in central Illinois

Bacterial communities and midgut microbiota associated with mosquito populations from waste tires in east-central Illinois

Dynamic gut microbiome across life history of the malaria mosquito Anopheles gambiae in Kenya

Bacterial diversity of field-caught mosquitoes, Aedes albopictus and Aedes aegypti, from different geographic regions of Madagascar

Host-feeding habits of Culex and other mosquitoes (Diptera: Culicidae) in the borough of queens in New York City, with characters and techniques for identification of Culex mosquitoes

Host blood-meal source has a strong impact on gut microbiota of Aedes aegypti

First isolation of microorganisms from the gut diverticulum of Aedes aegypti (Diptera: Culicidae): new perspectives for an insectbacteria association

Whole metagenome sequencing reveals links between mosquito microbiota and insecticide resistance in malaria vectors

Effect of life stage and pesticide exposure on the gut microbiota of Aedes albopictus and Culex pipiens

Midgut fungal and bacterial microbiota of Aedes triseriatus and Aedes japonicus shift in response to La Crosse virus infection

Serratia odorifera a midgut inhabitant of Aedes aegypti mosquito enhances its susceptibility to dengue-2 virus

Association between fertilizer-mediated changes in microbial communities and Aedes albopictus growth and survival

Investigation of the seasonal microbiome of Anopheles coluzzii mosquitoes in Mali

French invasive Asian tiger mosquito populations harbor reduced bacterial microbiota and genetic diversity compared to Vietnamese autochthonous relatives

The Anopheles mosquito microbiota and their impact on pathogen transmission

Mosquito vector-associated microbiota: metabarcoding bacteria and eukaryotic symbionts across habitat types in Thailand endemic for dengue and other arthropod-borne diseases

The Prokaryotes: A Handbook on the Biology of Bacteria: Proteobacteria: Gamma subclass

Microbial communities within fieldcollected Culiseta melanura and Coquillettidia perturbans

Wolbachia infection alters the relative abundance of resident bacteria in adult Aedes aegypti mosquitoes, but not larvae

Characterization of spirochetal isolates from arthropods collected in South Moravia, Czech Republic, using fatty acid methyl esters analysis

Microbiota variations in Culex nigripalpus disease vector mosquito of west nile virus and saint louis encephalitis from different geographic origins

Revision of the family Spironemidae Doflein (Protista, Hemimastigophora), with description of two new species, Spironema terricola n. sp. and Stereonema geiseri n. g., n. sp

Wolbachia: master manipulators of invertebrate biology

Investigating the presence of Wolbachia pipientis in various mosquito species

Wolbachia enhances west Nile virus (WNV) infection in the mosquito Culex tarsalis

Comparative analysis of gut microbiota of Culex restuans (Diptera: Culicidae) females from different parents

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

We thank Millon Blackshear for his technical assistance and all members of the Medical Entomology Laboratory (INHS), Yani Kaldis, Rosemary Ogbonna, and Eddie Lewis, for assisting with fieldwork and running the experiments, and members of Allan Laboratory for providing valuable reviews to the manuscript. The authors would like to thank Heather Walker of the National Center for Agricultural Utilization Research (U.S. Department of Agriculture) for helpful technical assistance.