<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Microbiol.</journal-id>
<journal-title>Frontiers in Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">1664-302X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmicb.2016.01485</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Microbiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Persistent Infection by <italic>Wolbachia w</italic>AlbB Has No Effect on Composition of the Gut Microbiota in Adult Female <italic>Anopheles stephensi</italic></article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name><surname>Chen</surname> <given-names>Shicheng</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>&#x002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/45445/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhao</surname> <given-names>Jiangchao</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Joshi</surname> <given-names>Deepak</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/353508/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Xi</surname> <given-names>Zhiyong</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Norman</surname> <given-names>Beth</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/368483/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Walker</surname> <given-names>Edward D.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/357617/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Microbiology and Molecular Genetics, Michigan State University, East Lansing</institution> <country>MI, USA</country></aff>
<aff id="aff2"><sup>2</sup><institution>Department of Animal Science, University of Arkansas, Fayetteville</institution> <country>AR, USA</country></aff>
<aff id="aff3"><sup>3</sup><institution>Department of Entomology, Michigan State University, East Lansing</institution> <country>MI, USA</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: <italic>Weiwen Zhang, Tianjin University, China</italic></p></fn>
<fn fn-type="edited-by"><p>Reviewed by: <italic>Michael Strand, University of Georgia, USA; Hao Song, Tianjin University, China</italic></p></fn>
<fn fn-type="corresp" id="fn001"><p>&#x002A;Correspondence: <italic>Shicheng Chen, <email>shicheng@msu.edu</email></italic></p></fn>
<fn fn-type="other" id="fn002"><p>This article was submitted to Microbial Physiology and Metabolism, a section of the journal Frontiers in Microbiology</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>21</day>
<month>09</month>
<year>2016</year>
</pub-date>
<pub-date pub-type="collection">
<year>2016</year>
</pub-date>
<volume>7</volume>
<elocation-id>1485</elocation-id>
<history>
<date date-type="received">
<day>02</day>
<month>06</month>
<year>2016</year>
</date>
<date date-type="accepted">
<day>07</day>
<month>09</month>
<year>2016</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2016 Chen, Zhao, Joshi, Xi, Norman and Walker.</copyright-statement>
<copyright-year>2016</copyright-year>
<copyright-holder>Chen, Zhao, Joshi, Xi, Norman and Walker</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p>The bacteria in the midgut of <italic>Anopheles stephensi</italic> adult females from laboratory colonies were studied by sequencing the V4 region of 16S rRNA genes, with respect to three experimental factors: stable or cured <italic>Wolbachia</italic> infection; sugar or blood diet; and age. <italic>Proteobacteria</italic> and <italic>Bacteroidetes</italic> dominated the community [>90% of operational taxonomic units (OTUs)]; most taxa were in the classes <italic>Flavobacteriia</italic>, <italic>Gammaproteobacteria</italic>, and <italic>Alphaproteobacteria</italic>, and were assigned to <italic>Elizabethkingia</italic> (46.9%), <italic>Asaia</italic> (6.4%) and <italic>Pseudomonas</italic> (6.0%), or unclassified <italic>Enterobacteriaceae</italic> (37.2%). Bacterial communities were similar between <italic>Wolbachia</italic>-cured and <italic>Wolbachia</italic>-infected mosquito lines, indicating that the gut microbiota were not dysregulated in the presence of <italic>Wolbachia</italic>. The proportion of <italic>Enterobacteriaceae</italic> was higher in mosquitoes fed a blood meal compared to those provided a sugar meal. Collectively, the bacterial community had a similar structure in older <italic>Wolbachia</italic>-infected mosquitoes 8 days after the blood meal, as in younger <italic>Wolbachia</italic>-infected mosquitoes before a blood meal, except that older mosquitoes had a higher proportion of <italic>Enterobacteriaceae</italic> and lower proportion of <italic>Elizabethkingia</italic>. Consistent presence of certain predominant bacteria (<italic>Elizabethkingia</italic>, <italic>Asaia</italic>, <italic>Pseudomonas</italic>, and <italic>Enterobacteriaceae</italic>) suggests they would be useful for paratransgenesis to control malaria infection, particularly when coupled to a <italic>Wolbachia</italic>-based intervention strategy.</p>
</abstract>
<kwd-group>
<kwd><italic>Anopheles</italic></kwd>
<kwd><italic>Wolbachia</italic></kwd>
<kwd>malaria vectors</kwd>
<kwd>paratransgenesis</kwd>
</kwd-group>
<counts>
<fig-count count="3"/>
<table-count count="3"/>
<equation-count count="0"/>
<ref-count count="35"/>
<page-count count="10"/>
<word-count count="0"/>
</counts>
</article-meta>
</front>
<body>
<sec><title>Introduction</title>
<p>Strategies for controlling mosquito-vectored pathogens based on infection with the symbiotic bacterium <italic>Wolbachia</italic> include vector population replacement and population suppression (<xref ref-type="bibr" rid="B30">Sinkins, 2004</xref>; <xref ref-type="bibr" rid="B34">Werren et al., 2008</xref>; <xref ref-type="bibr" rid="B19">Iturbe-Ormaetxe et al., 2011</xref>; <xref ref-type="bibr" rid="B6">Bourtzis et al., 2014</xref>). Although many mosquito species are naturally infected with various strains of <italic>Wolbachia</italic>, artificially but stably induced <italic>Wolbachia</italic> infection in mosquito host species not naturally having <italic>Wolbachia</italic> infection was originally established by embryonic injection followed by <italic>trans</italic>-generational spread (<xref ref-type="bibr" rid="B13">Dobson et al., 1999</xref>; <xref ref-type="bibr" rid="B30">Sinkins, 2004</xref>; <xref ref-type="bibr" rid="B6">Bourtzis et al., 2014</xref>). Certain combinations of <italic>Wolbachia</italic> strains and mosquito species suppress propagation and transmission of dengue viruses and malaria parasites in mosquito hosts (<xref ref-type="bibr" rid="B6">Bourtzis et al., 2014</xref>). However, in nature the host range of <italic>Wolbachia</italic> does not include some of the most important vectors of human pathogens, such as <italic>Anopheles stephensi, A. gambiae</italic>, and <italic>Aedes aegypti</italic> (<xref ref-type="bibr" rid="B17">Hilgenboecker et al., 2008</xref>). This host range has been expanded through forced infections by inoculations into mosquito embryos, resulting in stable infection across subsequent generations in <italic>A. aegypti</italic> and <italic>A. stephensi</italic> with little fitness cost (<xref ref-type="bibr" rid="B35">Xi et al., 2005</xref>; <xref ref-type="bibr" rid="B3">Bian et al., 2013</xref>; <xref ref-type="bibr" rid="B20">Joshi et al., 2014</xref>). A recent study found strongly negative interactions between the commensal gut microflora in mosquitoes and propensity of the mosquitoes to harbor <italic>Wolbachia</italic> infection after intrathoracic inoculation, which might explain why natural infections of <italic>Wolbachia</italic> in some mosquito species are not found in nature (<xref ref-type="bibr" rid="B18">Hughes et al., 2014</xref>). <xref ref-type="bibr" rid="B18">Hughes et al. (2014)</xref> reported that <italic>Wolbachia</italic> vertical transmission was achievable in <italic>Anopheles</italic> species only when the gut microbiota (especially, <italic>Asaia</italic> bacteria) were perturbed (i.e., eliminated) with antibiotics. Furthermore, <xref ref-type="bibr" rid="B28">Rossi et al. (2015)</xref> observed that <italic>Wolbachia</italic> failed to colonize the female reproductive organs of anophelines, due to presence of <italic>Asaia</italic> infection in those tissues (<xref ref-type="bibr" rid="B28">Rossi et al., 2015</xref>). The conclusion of both studies is that commensal microbiota impede <italic>Wolbachia</italic> vertical transmission in mosquitoes (<xref ref-type="bibr" rid="B18">Hughes et al., 2014</xref>; <xref ref-type="bibr" rid="B28">Rossi et al., 2015</xref>). However, what remains to be determined is whether the commensal microbiota of mosquitoes that have been stably infected <italic>trans</italic>-generationally with <italic>Wolbachia</italic> by embryonic injection differ from that of <italic>Wolbachia</italic>-uninfected mosquitoes of the same species and held under the same conditions. Based on the results of studies reviewed here (<xref ref-type="bibr" rid="B18">Hughes et al., 2014</xref>; <xref ref-type="bibr" rid="B28">Rossi et al., 2015</xref>), one would predict that stably infected mosquitoes must differ in their microbiota due to the inhibitory effects of <italic>Asaia</italic> and other gut symbionts on <italic>Wolbachia</italic> infection. Because stably infected mosquitoes are those most likely to be used in attempts to control vector populations or pathogen transmission, the interactions between stable <italic>Wolbachia</italic> infection and gut microbiota clearly become important.</p>
<p>After the <italic>w</italic>AlbB <italic>Wolbachia</italic> strain was successfully transferred from <italic>A. albopictus</italic> into <italic>A. stephensi</italic> by embryonic microinjection, this strain of <italic>Wolbachia</italic> spread rapidly into a laboratory population of <italic>A. stephensi</italic> by vertical transmission over the course of several generations without any manipulation of gut microbiota (<xref ref-type="bibr" rid="B3">Bian et al., 2013</xref>). Further, <italic>w</italic>AlbB infection conferred partial resistance to development of the malaria parasite <italic>Plasmodium falciparum</italic> in <italic>A. stephensi</italic> (<xref ref-type="bibr" rid="B3">Bian et al., 2013</xref>), and significantly up-regulated the mosquitoes&#x2019; innate immune response, including generation of reactive oxygen species from midgut epithelial cells (<xref ref-type="bibr" rid="B25">Pan et al., 2012</xref>). One way to increase the inhibitory effect of <italic>Wolbachia</italic> on virus or parasite development in mosquito hosts is to increase the intensity of <italic>Wolbachia</italic> infection (<xref ref-type="bibr" rid="B12">Correa and Ballard, 2014</xref>). Another strategy is to couple paratransgenesis of gut microbiota with a <italic>Wolbachia</italic>-based vector population replacement/suppression strategy. The latter concept is promising because certain symbionts such as <italic>Asaia</italic> and <italic>Pantoea agglomerans</italic>, when engineered to express molecules having anti-malaria parasite properties, result in reduced malaria parasite load in host mosquitoes (<xref ref-type="bibr" rid="B32">Wang et al., 2012</xref>; <xref ref-type="bibr" rid="B5">Bongio and Lampe, 2015</xref>). However, the innate immune response triggered by <italic>Wolbachia</italic> infection is non-specific and it may affect the commensal microflora (<xref ref-type="bibr" rid="B25">Pan et al., 2012</xref>). Thus, it remains unclear if introduced or genetically modified bacteria will persist in <italic>Wolbachia</italic>-infected mosquitoes due to elevated innate immune response and putative incompatibility between <italic>Wolbachia</italic> and gut microbiota. If <italic>Wolbachia w</italic>AlbB infection with its non-specific immune effects in <italic>A. stephensi</italic>, and engineered gut microbiota with specific anti-parasite effects, are found to be compatible; then a synergistic, dual strategy could be developed to suppress malaria parasite or virus development in mosquitoes. In this study, and as a step toward determining the extent of this compatibility, we quantified the diversity of the bacterial community in stably <italic>Wolbachia</italic>-infected and <italic>Wolbachia-</italic>cured <italic>A. stephensi</italic> mosquitoes. Moreover, because of the dynamics of microbiota associated with key life history events in mosquitoes, we further compared bacterial communities when sugar or meals were provided, and when the mosquitoes were young or aged.</p>
</sec>
<sec id="s1" sec-type="materials|methods">
<title>Materials and Methods</title>
<sec><title>Ethics Statement</title>
<p>All procedures involving vertebrate animals were in accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health, and were conducted under protocol 03/14-036-00, approved by the Michigan State University Institutional Animal Care and Use Committee.</p>
</sec>
<sec><title>Insects</title>
<p>Mosquito colonization procedures are described elsewhere (<xref ref-type="bibr" rid="B8">Chen et al., 2015</xref>). Briefly, adult mosquitoes were held in cages at 27&#x00B0;C and 85% humidity under a light/dark 12:12-h photoperiod without dawn/dusk transitions, and were provided a 10% sucrose solution <italic>ad libitum</italic> by cotton wick. To initiate egg development, bovine blood with sodium heparin (Hemostat Lab, Dixon, CA, USA) was fed to adult mosquitoes via an artificial membrane feeder. Two days after the blood meal, oviposition substrates (consisting of filter paper moistened with water in a Petri dish) were placed in the cages. Eggs were transferred to plastic containers with distilled water for hatching. First Bite (Kyorin, Himeji, Japan) and Tetramin tropical fish food flakes (Tetra, Blacksburg, VA, USA) were provided <italic>ad libitum</italic> to first instar larval mosquitoes. Thereafter, pet food (Purina Cat Chow; Nestle&#x00E9;) was given once per day <italic>ad libitum</italic>.</p>
<p>The stable <italic>w</italic>AlbB infection in the <italic>A. stephensi</italic> isofemale line (designated LB1) and the aposymbiotic <italic>A. stephensi</italic> line (<italic>w</italic>AlbB cured, designated LBT) were previously established (<xref ref-type="bibr" rid="B3">Bian et al., 2013</xref>). Aposymbiotic line LBT was derived from the LB1 strain after treatment with tetracycline (<xref ref-type="bibr" rid="B3">Bian et al., 2013</xref>). Approximately 4 years and 21 generations passed since that treatment was done, before this study was conducted. The presence or absence of <italic>Wolbachia</italic> in LB1 or LBT was confirmed by PCR of DNA extractions from abdomens using primers wsp81F (TGGTCC AATAAGTGATGAAGAAAC) and wsp691R (AAAAATTAAACGCTACTCCA). No <italic>Wolbachia</italic> infection was found in LBT individuals but all LB1 individuals were positive across the experimental groups (<bold>Supplementary Figure <xref ref-type="supplementary-material" rid="SM1">S1</xref></bold>). Females of each strain were sampled 1 week after adult emergence, designated here as &#x201C;LB1-SM&#x201D; and &#x201C;LBT-SM&#x201D; given that they had sugar meals but no blood meals yet. Other females were provided a blood meal at 7 days after adult emergence by allowing them to feed on anesthetized BALB/c mice for 20 min. These mosquitoes were sampled 24 h after the blood meal, and were designated &#x201C;LB1-BM&#x201D; and &#x201C;LBT-BM,&#x201D; respectively. Finally, other blood fed females of the LB1 and LBT strains were sampled after an interval of 1 week following the blood meal as &#x201C;aged&#x201D; mosquito samples (15 days after adult emergence and 8 days after the blood meal), during which time they were provisioned 10% sugar solution. These mosquitoes were designated &#x201C;LB1-PBM&#x201D; and &#x201C;LBT-PBM,&#x201D; respectively.</p>
</sec>
<sec><title>DNA Extraction, Library Construction, and 16S rRNA Sequencing</title>
<p>All DNA extractions were performed in a laminar flow biosafety cabinet to avoid contamination. Prior to dissection, mosquitoes were disinfected externally with 70% ethanol (three rinses), followed by a rinse with sterile water. Midguts were dissected with sterile forceps, and transferred to 200 &#x03BC;l of sterile phosphate-buffered saline (PBS) buffer. Midguts from six individuals were homogenized by crushing with a sterile pestle in the same tube as a single sample. The volume was re-suspended in 200 &#x03BC;l of lysis buffer and DNA extracted with the DNeasy Blood &#x0026; Tissue Kit (Qiagen) following the manufacturer&#x2019;s recommendations. Six biological replicates (total 36 midguts) were used for each of the six experimental treatments; a total of 216 mosquitoes were used in this study. The DNA concentration was measured using Qubit<sup>TM</sup> dsDNA HS Assay Kits and the DNA integrity was tested by PCR using primers 63F (CAGGCCTAACACATGCAAGTC) and 1387R (CGGAACATGTGWGGCGGG). Amplicon tagging and sequencing were conducted at the Research Technology Support Facility (RTSF) at Michigan State University. DNA was amplified by using the primers 515f (GTG CCA GCM GCC GCG GTA A) and 806r (NNN NNN GGA CTA CHV GGG TWT CTA AT), which targets the V4 region of bacterial 16S rDNA. The reverse primer contains a 6-bp error-correcting barcode unique to each sample. The purified amplicons were pooled, loaded on an Illumina MiSeq flow cell, and sequenced in a 2 &#x00D7; 250 bp paired end format using a 500 cycle v2 reagent cartridge. Base calling was done by Illumina Real Time Analysis (RTA) v1.18.54 and output of the RTA was demultiplexed and converted to FastQ format by Illumina Bcl2fastq v1.8.4.</p>
</sec>
<sec><title>DNA Sequence Analysis</title>
<p>Sequencing reads from the FastQ format were processed and analyzed using the Mothur package<sup><xref ref-type="fn" rid="fn01">1</xref></sup> v.1.37.0. After denoising procedures by PyroNoise, Uchime, preclustering, good quality sequences (>250 bp) without detectable sequencing errors or chimeras and removing rare operational taxonomic units (OTUs; 5 or fewer) were used for assigning OTUs using an average neighbor algorithm (97% similarity cutoff). OTUs were classified at the genus level using the Bayesian method. When there were no references for the interested bacteria in the database used by Mothur, we retrieved the representative sequences for specific OTUs and submitted the sequences to GenBank (NCBI) for a BLAST search for purposes of classification. Data were further analyzed by first trimming sequence quality with different cutoffs using standard filtering tools (<xref ref-type="bibr" rid="B29">Schloss et al., 2009</xref>). Rarefaction curves were built to estimate sample coverage. Non-parametric estimates of richness were made using the Chao index and Abundance Coverage Estimator (cut-off threshold, 97%). Within-sample diversity (or &#x03B1;-diversity) was calculated using Simpson and Shannon diversity indices. Welch&#x2019;s <italic>t</italic>-test was used to compare diversity indexes between <italic>Wolbachia</italic>-infected and <italic>Wolbachia</italic>-cured mosquitoes with different diets and at different rearing ages. Community composition (or &#x03B2;-diversity) was compared among treatments using Bray&#x2013;Curtis dissimilarities. Differences in community composition among treatments were tested using permutational multivariate analysis of variance (PERMANOVA). Richness and &#x03B1;-diversity analyses were done in mothur. &#x03B2;-diversity analyses were done in R (v. 3.2.0). Bray&#x2013;Curtis dissimilarity matrices were created using the vegdist() function (method = &#x201C;bray&#x201D;) and PERMANOVA tests (999 permutations) were done using the adonis() function in the vegan package [v. 2.3-2 (<xref ref-type="bibr" rid="B24">Oksanen et al., 2015</xref>)]. Principal coordinate analyses (PCoA) ordinations were created using the pcoa() function in the ape package [v. 3.4 (<xref ref-type="bibr" rid="B26">Paradis et al., 2004</xref>)].</p>
</sec>
<sec><title>Nucleotide Sequence Accession Number</title>
<p>The raw sequences of this study have been deposited in the Sequence Read Archive (accession number: SRP072068).</p>
</sec>
</sec>
<sec><title>Results</title>
<sec><title>Characteristics of Midgut Bacteria Community Libraries</title>
<p>Approximately 2,130,826 sequence reads were successfully retained. After removal of short (&#x003C;250 bp), possibly chimeras and rare OTUs (5 or fewer, see <bold>Supplementary Table <xref ref-type="supplementary-material" rid="SM1">S1</xref></bold>), a total of 1,704,538 good quality sequences were available for further analysis (<bold>Supplementary Table <xref ref-type="supplementary-material" rid="SM1">S2</xref></bold>). These sequences were assigned to 394 OTUs (<bold>Supplementary Table <xref ref-type="supplementary-material" rid="SM1">S1</xref></bold>). Sequence number for each sample was normalized to the minimal readings (6,056) by randomly subsampling to minimize the biases generated by sequencing depth (<bold>Supplementary Table <xref ref-type="supplementary-material" rid="SM1">S2</xref></bold>). The average Good&#x2019;s coverage was 99.9% (<bold>Supplementary Table <xref ref-type="supplementary-material" rid="SM1">S2</xref></bold>). The rarefaction curves (<bold>Figure <xref ref-type="fig" rid="F1">1</xref></bold>) showed that the samples from some of these mosquito collections did not reach saturation to an asymptote, indicating that additional rare bacterial taxa likely occur in them, but all deflected as they rose with sampling effort in a characteristic concave function (<bold>Figure <xref ref-type="fig" rid="F1">1</xref></bold>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption><p><bold>Rarefaction analysis of bacterial 16S rRNA gene libraries from <italic>Wolbachia</italic>-infected and <italic>Wolbachia-</italic>cured <italic>A. stephensi</italic> samples that were sugar fed, blood fed, and post-blood fed.</bold> Operational taxonomic units (OTUs) were grouped with a 97% similarity. LB1-SM, <italic>Wolbachia</italic>-infected mosquitoes with sugar meal; LBT-SM, <italic>Wolbachia</italic>-cured mosquitoes with sugar meal; LB1-BM, <italic>Wolbachia</italic>-infected mosquitoes with blood meal; LBT-BM, <italic>Wolbachia-</italic>cured mosquitoes with blood meal; LB1-PBM, post-blood fed <italic>Wolbachia</italic>-infected mosquitoes; LBT-PBM, post-blood fed <italic>Wolbachia-</italic>cured mosquitoes.</p></caption>
<graphic xlink:href="fmicb-07-01485-g001.tif"/>
</fig>
</sec>
<sec><title>Effects of <italic>Wolbachia</italic> Infection, Diet, and Age on Composition of the Gut Microbiota</title>
<p>Operational taxonomic units were assigned to over 16 bacterial phyla and to unclassified bacteria (<bold>Supplementary Table <xref ref-type="supplementary-material" rid="SM1">S3</xref></bold>). Single OTUs assigned to the classes <italic>Flavobacteriia</italic>, <italic>Gammaproteobacteria</italic>, and <italic>Alphaproteobacteria</italic> were most frequent, together accounting for more than 90% of the sequences, and were dominant in experimental treatments (<bold>Figure <xref ref-type="fig" rid="F2">2A</xref></bold>). Sequences of bacteria in the classes <italic>Actinobacteria</italic> and <italic>Betaproteobacteria</italic> were also detected, but they were not uniformly represented in every sample. Some bacteria (2% of the total community) could not be assigned to a certain rank (unclassified), suggesting the existence of previously uncharacterized taxa. However, OTUs assigned to other classes were not abundant (i.e., &#x003C;1%) and/or were absent from some samples or groups. The top 10 OTUs at the genus level comprised 99.2% of the total sequences (<bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>). The predominant OTUs were represented by <italic>Elizabethkingia</italic> (46.9%), unclassified <italic>Enterobacteriaceae</italic> (37.2%), <italic>Asaia</italic> (6.4%) and <italic>Pseudomonas</italic> (6.0%).</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption><p><bold>Gut bacterial composition at class level and the relative abundance of the top 10 genera in <italic>Wolbachia</italic>-infected and <italic>Wolbachia-</italic>cured mosquito samples by dietary treatment.</bold> <bold>(A)</bold> Taxonomic classification of bacterial reads retrieved from the <italic>Wolbachia</italic>-infected and <italic>Wolbachia-</italic>cured mosquitoes. <bold>(B)</bold> Relative abundance of 10 most abundant OTUs (97% similarity). Unclassified <italic>Enterobacteriaceae</italic> or <italic>Micrococcaceae</italic> represent the reads that were assigned to the family <italic>Enterobacteraceae</italic> or <italic>Micrococcaceae</italic>, but could not be assigned to a genus. <italic>Wol</italic>+, <italic>Wolbachia</italic>-infected mosquito; <italic>Wol</italic>-, <italic>Wolbachia</italic>-cured mosquito.</p></caption>
<graphic xlink:href="fmicb-07-01485-g002.tif"/>
</fig>
<p>There was no significant difference in bacterial class abundance between <italic>Wolbachia-</italic>infected (LB1-SM) and <italic>Wolbachia-</italic>cured mosquitoes (LBT-SM) (Welch&#x2019;s <italic>t</italic>-test) (<bold>Figure <xref ref-type="fig" rid="F2">2A</xref></bold>). At the genus level, the most abundant OTU in the midgut of <italic>A. stephensi</italic> LBT (<italic>Wolbachia-</italic>cured) was <italic>Elizabethkingia</italic>, accounting for 65% of the community (<bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>). A similar proportion of <italic>Elizabethkingia</italic> was found in <italic>A. stephensi</italic> LB1 (69%), indicating that <italic>Wolbachia</italic> infection did not significantly alter its relative abundance (<italic>p</italic> > 0.05). Next to <italic>Elizabethkingia</italic>, unclassified <italic>Enterobacteriaceae</italic> were the second most abundant OTU in sugar-fed young <italic>Anopheles</italic> mosquitoes (LB1-SM and LBT-SM), which accounted for between 8 and 15% of the OTUs (<bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>). The number of <italic>Enterobacteriaceae</italic> OTUs was not statistically different between LB1-SM and LBT-SM (<italic>p</italic> > 0.05; <bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>). <italic>Asaia</italic> OTUs ranged in abundance from 5 to 8% of the total bacterial community (<bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>) in LB1-SM and LBT-SM, respectively, indicating that the colonization and persistence of <italic>Asaia</italic> in <italic>A. stephensi</italic> midguts occurred in a stable <italic>Wolbachia</italic> infection mosquito line as well as in <italic>A. stephensi</italic> without <italic>Wolbachia</italic> infection (<italic>p</italic> > 0.05; <bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>). Similar to <italic>Asaia</italic>, <italic>Pseudomonas</italic> OTUs were equally associated with both mosquito lines (<italic>p</italic> > 0.05; <bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>). As expected, the relative abundance of <italic>Wolbachia</italic> was 8.6% of the total bacterial OTUs in LB1-SM while it was nearly absent (&#x003C;0.01% of OTUs) in LBT-BM (<bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>).</p>
<p>At 24 h after a blood meal, the relative abundance of <italic>Elizabethkingia</italic> decreased from 65 to 19% in LBT (<italic>p</italic> &#x003C; 0.05) and decreased from 69 to 34% in LB1 mosquitoes, respectively (<italic>p</italic> &#x003C; 0.05). By contrast, after a blood meal the <italic>Enterobacteriaceae</italic> OTUs increased 7.1 and 4.4-fold in LB1-BM and LBT-BM within 24 h (<italic>p</italic> &#x003C; 0.05), respectively. Compared to sugar meals, the relative abundance of <italic>Asaia</italic> or <italic>Pseudomonas</italic> was not significantly changed by blood feeding (<italic>p</italic> > 0.05; <bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>). The relative abundance of <italic>Wolbachia</italic> in the midgut decreased to 0.18% of the total taxa in 24 h after blood meals (<bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>).</p>
<p>In mosquitoes at age 15 days (that is, 8 days post-blood meal), the frequency of <italic>Elizabethkingia</italic> OTUs represented 42% of the bacterial community in guts of <italic>Wolbachia</italic>-infected mosquitoes (LB1-PBM), showing that relative abundance of these OTUs in 15 days old LB1-PBM was significantly different from 7 days old LBT-SM mosquitoes (<italic>p</italic> &#x003C; 0.05; <bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>). However, there was no statistical difference between LBT-SM and LBT-PBM mosquitoes (<italic>p</italic> > 0.05; <bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>). The relative abundance of <italic>Enterobacteriaceae</italic> OTUs in 15 days old LBT and LB1 mosquitoes was 2.2 and 3.4-fold higher than that in young LBT and LB1 mosquitoes (<italic>p</italic> &#x003C; 0.05), respectively. The abundance of <italic>Asaia</italic> was consistent and not affected by age (<italic>p</italic> > 0.05), and accounted for 8.9 and 6.5% of bacterial community in older <italic>Wolbachia</italic>-infected and <italic>Wolbachia-</italic>cured mosquitoes, respectively. The abundance of <italic>Pseudomonas</italic> in older, <italic>Wolbachia</italic>-infected mosquitoes (LBT-PBM) was 2.6-fold higher than in young LBT mosquitoes (<italic>p</italic> &#x003C; 0.05), while it remained stable between young and old, <italic>Wolbachia-</italic>cured mosquitoes (<italic>p</italic> > 0.05; <bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>). The relative abundance of <italic>Wolbachia</italic> OTUs in LB1-PBM reached up to 7.4% of total bacterial community level, which was comparable to that in LB1-SM (<italic>p</italic> > 0.05; <bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>).</p>
</sec>
<sec><title>Alpha and Beta Bacterial Diversity in <italic>Wolbachia</italic>-Infected and <italic>Wolbachia-</italic>Cured Mosquitoes of Different Diets and Ages</title>
<p>To analyze differences in alpha diversity, we compared diversity (Shannon index) and richness (Chao1). There was no significant difference in Shannon diversity between LB1-SM and LBT-SM or LB1-BM and LBT-BM, showing that intracellular <italic>Wolbachia</italic> infection did not affect microbial diversity in young mosquito midguts (<bold>Table <xref ref-type="table" rid="T1">1</xref></bold>; <bold>Supplementary Figure <xref ref-type="supplementary-material" rid="SM1">S2</xref></bold>). However, the Shannon index was different between LB1-PBM and LBT-PBM samples (<italic>p</italic> &#x003C; 0.05), indicating that <italic>Wolbachia</italic> infection only changed gut microbial diversity in older mosquito. After blood feeding, the Shannon diversity in LB1-BM was not significantly different from that in LB1-SM (<italic>p</italic> > 0.05) and there was no significant difference between LBT-BM and LBT-SM (<italic>p</italic> > 0.05) (<bold>Table <xref ref-type="table" rid="T1">1</xref></bold>; <bold>Supplementary Figure <xref ref-type="supplementary-material" rid="SM1">S2</xref></bold>). Thirdly, microbial diversity was significantly different between older mosquitoes (LBT-PBM) and young mosquitoes (LBT-SM) with <italic>Wolbachia</italic> infection (<italic>p</italic> &#x003C; 0.05). However, in non-<italic>Wolbachia</italic> infection mosquitoes, the microbial diversity was similar between older mosquitoes (LB1-PBM) and young mosquitoes (LB1-SM) (p &#x003C; 0.05; <bold>Table <xref ref-type="table" rid="T1">1</xref></bold>; <bold>Supplementary Figure <xref ref-type="supplementary-material" rid="SM1">S2</xref></bold>). For community richness comparisons (Chao1), no significant differences were observed between <italic>Wolbachia</italic>-infected mosquitoes (LB1-SM, LB1-BM, or LB1-PBM) and <italic>Wolbachia-</italic>cured ones (LBT-SM, LBT-BM, or LBT-PBM) (p > 0.05; <bold>Supplementary Figure <xref ref-type="supplementary-material" rid="SM1">S2</xref></bold>). However, blood-fed mosquitoes with <italic>Wolbachia</italic> infection had a significantly lower number of estimated OTUs (Chao1) than sugar-fed ones with <italic>Wolbachia</italic> infection (i.e., LB1-BM vs. LB1-SM, <italic>p</italic> &#x003C; 0.05) while there was no significant difference between blood-fed and sugar-fed mosquitoes without <italic>Wolbachia</italic> infection (LBT-BM vs. LBT-SM, <italic>p</italic> > 0.05) as shown in <bold>Supplementary Figure <xref ref-type="supplementary-material" rid="SM1">S2</xref></bold>, indicating that the blood meal decreased richness in <italic>Wolbachia</italic> infected mosquitoes. However, the richness (Chao1) in the guts of older mosquitoes (LB1-PBM or LBT-PBM) was not significantly different from that in of younger ones (LB1-BM or LBT-BM) after blooding feeding (<italic>p</italic> > 0.05; <bold>Supplementary Figure <xref ref-type="supplementary-material" rid="SM1">S2</xref></bold>).</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>Richness and diversity estimation of the 16S rRNA gene libraries.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left">Sample</th>
<th valign="top" align="center">Cutoffs</th>
<th valign="top" align="center">Simpson</th>
<th valign="top" align="center">Shannon</th>
<th valign="top" align="center">Chao</th>
<th valign="top" align="center">Ace</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">LBT-SM</td>
<td valign="top" align="center">0.03</td>
<td valign="top" align="center">0.58 &#x00B1; 0.23</td>
<td valign="top" align="center">0.99 &#x00B1; 0.71</td>
<td valign="top" align="center">32.15 &#x00B1; 14.16</td>
<td valign="top" align="left">41.24 &#x00B1; 14.35</td>
</tr>
<tr>
<td valign="top" align="left">LB1-SM</td>
<td valign="top" align="center">0.03</td>
<td valign="top" align="center">0.50 &#x00B1; 0.01</td>
<td valign="top" align="center">1.13 &#x00B1; 0.23</td>
<td valign="top" align="center">39.91 &#x00B1; 10.31</td>
<td valign="top" align="left">41.23 &#x00B1; 10.07</td>
</tr>
<tr>
<td valign="top" align="left">LBT-BM</td>
<td valign="top" align="center">0.03</td>
<td valign="top" align="center">0.57 &#x00B1; 0.12</td>
<td valign="top" align="center">0.78 &#x00B1; 0.19</td>
<td valign="top" align="center">39.81 &#x00B1; 32.14</td>
<td valign="top" align="left">22.07 &#x00B1; 25.94</td>
</tr>
<tr>
<td valign="top" align="left">LB1-BM</td>
<td valign="top" align="center">0.03</td>
<td valign="top" align="center">0.50 &#x00B1; 0.11</td>
<td valign="top" align="center">0.92 &#x00B1; 0.16</td>
<td valign="top" align="center">25.96 &#x00B1; 11.81</td>
<td valign="top" align="left">40.49 &#x00B1; 24.43</td>
</tr>
<tr>
<td valign="top" align="left">LBT-PBM</td>
<td valign="top" align="center">0.03</td>
<td valign="top" align="center">0.41 &#x00B1; 0.04</td>
<td valign="top" align="center">1.12 &#x00B1; 0.06</td>
<td valign="top" align="center">39.23 &#x00B1; 13.51</td>
<td valign="top" align="left">64.59 &#x00B1; 24.09</td>
</tr>
<tr>
<td valign="top" align="left">LB1-PBM</td>
<td valign="top" align="center">0.03</td>
<td valign="top" align="center">0.33 &#x00B1; 0.07</td>
<td valign="top" align="center">1.56 &#x00B1; 0.31</td>
<td valign="top" align="center">36.76 &#x00B1; 10.51</td>
<td valign="top" align="left">37.21 &#x00B1; 9.64</td>
</tr>
<tr>
<td valign="top" align="left"></td>
</tr>
</tbody>
</table>
</table-wrap>
<p>For beta diversity analysis, we examined the relationships in gut microbiota between different diet, age and <italic>Wolbachia</italic> infection status by using principal coordinate analyses (PCoA) based on Bray&#x2013;Curtis distance (<bold>Figure <xref ref-type="fig" rid="F3">3</xref></bold>). The first two components captured 90% of the variance among samples. Community composition differed among all three diet types (PERMANOVA; <italic>F</italic> = 10.93, <italic>p</italic> &#x003C; 0.05, <bold>Table <xref ref-type="table" rid="T2">2</xref></bold>). Bacterial communities from the mosquitoes fed on sugar separated from those fed on blood and those PBM along PC1, which explained 79% of the variation (<bold>Table <xref ref-type="table" rid="T3">3</xref></bold>). However, community composition was similar between <italic>Wolbachia</italic>-infected and <italic>Wolbachia-</italic>cured groups across all diet treatments (PERMANOVA; <italic>F</italic> = 2.86, <italic>p</italic> > 0.05, <bold>Table <xref ref-type="table" rid="T2">2</xref></bold>).</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption><p><bold>Principal coordinate analysis of bacterial the community structure using Bray&#x2013;Curtis distances.</bold> Symbols represent one mosquito collection (six pooled mosquito guts). Distances between symbols on the ordination plot reflect relative dissimilarities in community structures.</p></caption>
<graphic xlink:href="fmicb-07-01485-g003.tif"/>
</fig>
<table-wrap position="float" id="T2">
<label>Table 2</label>
<caption><p>Pseudo <italic>F</italic> table of PERMANOVA analysis based on Bray&#x2013;Curtis dissimilarities.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left">Source of variance</th>
<th valign="top" align="center">Sum of squares</th>
<th valign="top" align="center">Degrees of freedom</th>
<th valign="top" align="center">Mean square</th>
<th valign="top" align="center"><italic>F</italic></th>
<th valign="top" align="center"><italic>R</italic><sup>2</sup></th>
<th valign="top" align="center"><italic>p</italic><sup>&#x2217;</sup></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left"><italic>Wolbachia</italic></td>
<td valign="top" align="center">0.186</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">0.186</td>
<td valign="top" align="center">2.861</td>
<td valign="top" align="center">0.053</td>
<td valign="top" align="center">0.086</td>
</tr>
<tr>
<td valign="top" align="left">Diet</td>
<td valign="top" align="center">1.418</td>
<td valign="top" align="center">2</td>
<td valign="top" align="center">0.709</td>
<td valign="top" align="center">10.931</td>
<td valign="top" align="center">0.404</td>
<td valign="top" align="center"><bold>0.001</bold></td>
</tr>
<tr>
<td valign="top" align="left"><italic>Wolbachia</italic><sup>&#x2217;</sup>Diet</td>
<td valign="top" align="center">0.088</td>
<td valign="top" align="center">2</td>
<td valign="top" align="center">0.044</td>
<td valign="top" align="center">0.681</td>
<td valign="top" align="center">0.025</td>
<td valign="top" align="center">0.549</td>
</tr>
<tr>
<td valign="top" align="left">Residuals</td>
<td valign="top" align="center">1.816</td>
<td valign="top" align="center">28</td>
<td valign="top" align="center">0.065</td>
<td valign="top" align="center"></td>
<td valign="top" align="center">0.518</td>
<td valign="top" align="center"></td>
</tr>
<tr>
<td valign="top" align="left">Total</td>
<td valign="top" align="center">3.509</td>
<td valign="top" align="center">33</td>
<td valign="top" align="center"></td>
<td valign="top" align="center"></td>
<td valign="top" align="center">1</td>
<td valign="top" align="center"></td>
</tr>
<tr>
<td valign="top" align="left"></td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<attrib><italic><sup>&#x2217;</sup>OTUs with 5 or fewer sequences were excluded. Significant <italic>p</italic>-values (&#x2264;0.05) are bolded.</italic></attrib>
</table-wrap-foot>
</table-wrap>
<table-wrap position="float" id="T3">
<label>Table 3</label>
<caption><p>Pseudo <italic>F</italic> table of pairwise comparisons of diet types using PERMANOVA.</p></caption>
<table cellspacing="5" cellpadding="5" frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left">Pairwise comparison</th>
<th valign="top" align="left">Source of variance</th>
<th valign="top" align="center">Sum of squares</th>
<th valign="top" align="center">Degrees of freedom</th>
<th valign="top" align="center">Mean square</th>
<th valign="top" align="center"><italic>F</italic></th>
<th valign="top" align="center"><italic>R</italic><sup>2</sup></th>
<th valign="top" align="center"><italic>p</italic><sup>&#x2217;</sup></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">SM vs. BM</td>
<td valign="top" align="left">Diet</td>
<td valign="top" align="center">1.314</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">1.314</td>
<td valign="top" align="center">19.660</td>
<td valign="top" align="center">0.496</td>
<td valign="top" align="center"><bold>0.001</bold></td>
</tr>
<tr>
<td valign="top" align="left"></td>
<td valign="top" align="left">Residuals</td>
<td valign="top" align="center">1.337</td>
<td valign="top" align="center">20</td>
<td valign="top" align="center">0.067</td>
<td valign="top" align="center"></td>
<td valign="top" align="center">0.504</td>
<td valign="top" align="center"></td>
</tr>
<tr>
<td valign="top" align="left"></td>
<td valign="top" align="left">Total</td>
<td valign="top" align="center">2.652</td>
<td valign="top" align="center">21</td>
<td valign="top" align="center"></td>
<td valign="top" align="center"></td>
<td valign="top" align="center">1</td>
<td valign="top" align="center"></td>
</tr>
<tr>
<td valign="top" align="left">SM vs. PBM</td>
<td valign="top" align="left">Diet</td>
<td valign="top" align="center">0.308</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">0.308</td>
<td valign="top" align="center">4.261</td>
<td valign="top" align="center">0.176</td>
<td valign="top" align="center"><bold>0.021</bold></td>
</tr>
<tr>
<td valign="top" align="left"></td>
<td valign="top" align="left">Residuals</td>
<td valign="top" align="center">1.444</td>
<td valign="top" align="center">20</td>
<td valign="top" align="center">0.072</td>
<td valign="top" align="center"></td>
<td valign="top" align="center">0.824</td>
<td valign="top" align="center"></td>
</tr>
<tr>
<td valign="top" align="left"></td>
<td valign="top" align="left">Total</td>
<td valign="top" align="center">1.752</td>
<td valign="top" align="center">21</td>
<td valign="top" align="center"></td>
<td valign="top" align="center"></td>
<td valign="top" align="center">1</td>
<td valign="top" align="center"></td>
</tr>
<tr>
<td valign="top" align="left">BM vs. PBM</td>
<td valign="top" align="left">Diet</td>
<td valign="top" align="center">0.522</td>
<td valign="top" align="center">1</td>
<td valign="top" align="center">0.522</td>
<td valign="top" align="center">8.208</td>
<td valign="top" align="center">0.272</td>
<td valign="top" align="center"><bold>0.004</bold></td>
</tr>
<tr>
<td valign="top" align="left"></td>
<td valign="top" align="left">Residuals</td>
<td valign="top" align="center">1.399</td>
<td valign="top" align="center">22</td>
<td valign="top" align="center">0.064</td>
<td valign="top" align="center"></td>
<td valign="top" align="center">0.728</td>
<td valign="top" align="center"></td>
</tr>
<tr>
<td valign="top" align="left"></td>
<td valign="top" align="left">Total</td>
<td valign="top" align="center">1.921</td>
<td valign="top" align="center">23</td>
<td valign="top" align="center"></td>
<td valign="top" align="center"></td>
<td valign="top" align="center">1</td>
<td valign="top" align="center"></td>
</tr>
<tr>
<td valign="top" align="left"></td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<attrib><italic><sup>&#x2217;</sup>OTUs with 5 or fewer sequences were excluded. Significant <italic>p</italic>-values (&#x2264;0.05) are bolded.</italic></attrib>
</table-wrap-foot>
</table-wrap>
</sec>
</sec>
<sec><title>Discussion</title>
<p>The bacteria (i.e., microbiota) associated with mosquitoes have significant roles throughout the mosquito lifecycle, in providing energy and nutrients, and regulating development, fecundity, and immune responses (<xref ref-type="bibr" rid="B30">Sinkins, 2004</xref>; <xref ref-type="bibr" rid="B9">Chouaia et al., 2012</xref>; <xref ref-type="bibr" rid="B25">Pan et al., 2012</xref>; <xref ref-type="bibr" rid="B11">Coon et al., 2014</xref>). These symbiotic bacteria find shelter and nutrients within this environment and could properly be referred to as commensals for these beneficial reasons, even if the beneficial associations of the bacteria to mosquitoes are not entirely clear. Stable, intracellular <italic>Wolbachia</italic> infection in <italic>A. stephensi</italic> (LB1) was achieved previously through embryonic inoculation and multi-generational, vertical passage (<xref ref-type="bibr" rid="B3">Bian et al., 2013</xref>). It was done without perturbation of the natural microbiota (as shown here) and without excessive mortality after blood meals (<xref ref-type="bibr" rid="B3">Bian et al., 2013</xref>; <xref ref-type="bibr" rid="B20">Joshi et al., 2014</xref>). More recently and by contrast, vertical persistence of <italic>Wolbachia</italic> infection in <italic>A. stephensi</italic> was achieved only when symbiotic bacteria were perturbed by antibiotics (<xref ref-type="bibr" rid="B18">Hughes et al., 2014</xref>). In that same study, a combination of <italic>Wolbachia</italic> and <italic>Asaia</italic> caused mosquito mortality after a blood meal, a finding suggestive of a physiologic incompatibility between <italic>Wolbachia</italic> infection and presence of a constitutive microflora when the stressors of the blood meal arrive (<xref ref-type="bibr" rid="B18">Hughes et al., 2014</xref>). Indeed, if such an incompatibility exists, it would obviate the possibility of a dual, anti-malaria parasite strategy based on immune-mediating effects of <italic>Wolbachia</italic> infection and paratransgenesis of select members of that microflora. This possibility, in part, motivated the research presented here.</p>
<p>The similarity of the gut bacterial community structure, as measured by composition of taxa as well as by &#x03B1;- and &#x03B2;-diversity, between <italic>Wolbachia</italic>-infected LB1 and <italic>Wolbachia-</italic>cured LBT mosquitoes indicates that there were no significant <italic>Wolbachia</italic>-mediated effects on the structure of the commensal bacterial community; and that the bacterial assembly was resilient to any effects, such as innate immune effects, that the presence of <italic>Wolbachia</italic> might impose. Therefore, coexistence of these gut bacteria with <italic>Wolbachia</italic> was not problematic in the stable <italic>Wolbachia</italic>-infected <italic>Anopheles</italic> progeny. These findings contrast with other studies demonstrating very strongly negative effects of combination of gut bacteria and <italic>Wolbachia</italic> on host fitness (high mortality rate after a blood meal), particularly due to <italic>Asaia</italic> bacteria (<xref ref-type="bibr" rid="B18">Hughes et al., 2014</xref>). One explanation for the difference is in the way the <italic>Wolbachia</italic> infection was originally introduced into the mosquitoes: embryonic microinjection followed by stable, <italic>trans</italic>-generational perpetuation of infection (<xref ref-type="bibr" rid="B3">Bian et al., 2013</xref>); or intrathoracic inoculation without perpetuation (<xref ref-type="bibr" rid="B18">Hughes et al., 2014</xref>). The latter likely led to a larger and more acute dose of the bacteria, than one inherited in the germ line. However, the detailed mechanisms need to be further investigated.</p>
<p>The blood meal profoundly affects bacterial community composition in the midgut of several mosquito species, including <italic>A. gambiae</italic>, <italic>A. stephensi</italic>, and <italic>A. coluzzii</italic> (<xref ref-type="bibr" rid="B33">Wang et al., 2011</xref>; <xref ref-type="bibr" rid="B16">Gimonneau et al., 2014</xref>; <xref ref-type="bibr" rid="B31">Tchioffo et al., 2016</xref>). The main influencing factors included dramatic changes in mosquito gut temperature and appearance of oxidative stress (<xref ref-type="bibr" rid="B22">Luckhart et al., 1998</xref>; <xref ref-type="bibr" rid="B27">Peterson and Luckhart, 2006</xref>; <xref ref-type="bibr" rid="B2">Benoit et al., 2011</xref>). Especially, the ingested red blood cells from animals carry an enormous amount of hemoglobin, which causes a massive release of heme in the midgut, thus leading to a dramatic change in gut conditions (<xref ref-type="bibr" rid="B33">Wang et al., 2011</xref>; <xref ref-type="bibr" rid="B31">Tchioffo et al., 2016</xref>). Blood meals induce midgut epithelia to produce nitric oxide which is also a source for free radicals (<xref ref-type="bibr" rid="B22">Luckhart et al., 1998</xref>). Hematophagous mosquitoes exhibit detoxification mechanisms in response to these free radicals, but the role of commensal bacteria in the gut in that same process and how those bacteria survive them are not known (<xref ref-type="bibr" rid="B27">Peterson and Luckhart, 2006</xref>; <xref ref-type="bibr" rid="B2">Benoit et al., 2011</xref>). After blood meals, <italic>Enterobacteriaceae</italic> dramatically increased, which was consistent with those reported in <italic>A. gambiae</italic> (<xref ref-type="bibr" rid="B33">Wang et al., 2011</xref>). <xref ref-type="bibr" rid="B33">Wang et al. (2011)</xref> found several stress response systems existing in the <italic>Enterobacteriaceae</italic> genomes. Possibly, this group of bacteria interacts with the mosquito host in response to the oxidative stressors present in the blood meal bolus (<xref ref-type="bibr" rid="B33">Wang et al., 2011</xref>). Therefore, blood-induced mortality in transiently <italic>Wolbachia</italic>-infected <italic>Anopheles</italic> mosquito could be caused by a combinative effect, that is, innate immune response elicited by <italic>Wolbachia</italic> infection and by dysregulated microbiota (<xref ref-type="bibr" rid="B22">Luckhart et al., 1998</xref>; <xref ref-type="bibr" rid="B25">Pan et al., 2012</xref>; <xref ref-type="bibr" rid="B18">Hughes et al., 2014</xref>). Both of them caused an unusual high level of active radicals in midgut, which can cause destructive damage to mosquitoes (<xref ref-type="bibr" rid="B22">Luckhart et al., 1998</xref>). Although there was a dramatic change in microbial structure after blood meals (compared to sugar meals), the core microbial compositions and community structure between <italic>Wolbachia</italic>-infected and <italic>Wolbachia-</italic>cured mosquitoes were rather similar (<bold>Figure <xref ref-type="fig" rid="F3">3</xref></bold>).</p>
<p>When switching diet from blood to sugar meals and aging the insects, we found that the microbial community structure was different between the young (before blood meal) and aged mosquitoes (post-blood meal) in both mosquito lines. However, the relative abundance of transiently increased bacteria such as <italic>Enterobacteriaceae</italic> decreased when blood digestion was complete; the dominant <italic>Elizabethkingia</italic> abundance at age 15 days, and after the blood meal, was similar to before the blood meal (<bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>). The dominant bacterial community patterns we observed here were similar to those reported in <italic>A. gambiae</italic> (<xref ref-type="bibr" rid="B33">Wang et al., 2011</xref>). Collectively, <italic>Asaia</italic>, <italic>Elizabethkingia</italic>, and <italic>Enterobacteriaceae</italic> were present and predominant in aged, blood-fed, <italic>Wolbachia</italic>-infected <italic>A. stephensi</italic>, providing the possibility to devise an efficient malaria control reagent based on a combination of effects of <italic>Wolbachia</italic> and effects of paratransgenic gut microbiota (<xref ref-type="bibr" rid="B15">Federici et al., 2003</xref>; <xref ref-type="bibr" rid="B14">Favia et al., 2008</xref>; <xref ref-type="bibr" rid="B10">Cirimotich et al., 2011</xref>; <xref ref-type="bibr" rid="B4">Boissi&#x00E8;re et al., 2012</xref>; <xref ref-type="bibr" rid="B7">Capone et al., 2013</xref>; <xref ref-type="bibr" rid="B5">Bongio and Lampe, 2015</xref>; <xref ref-type="bibr" rid="B8">Chen et al., 2015</xref>).</p>
<p>Previous studies demonstrated that <italic>Wolbachia</italic> infected the midgut in <italic>A. stephensi</italic> and <italic>A. gambiae</italic> mosquitoes though its distribution level was much lower than that in the fat body (<xref ref-type="bibr" rid="B3">Bian et al., 2013</xref>). After a blood meal, the dramatic decrease of <italic>Wolbachia</italic> OTUs in LB1 gut was consistent with a previous observation by <xref ref-type="bibr" rid="B18">Hughes et al. (2014)</xref>. Several factors including heme release, ROS level change, and interactions with gut bacteria contributed to this change (<xref ref-type="bibr" rid="B25">Pan et al., 2012</xref>; <xref ref-type="bibr" rid="B3">Bian et al., 2013</xref>; <xref ref-type="bibr" rid="B18">Hughes et al., 2014</xref>). However, after removing gut bacteria by antibiotic treatment (especially <italic>Asaia</italic>), <italic>Wolbachia</italic> persisted at a level similar to that non-blood fed mosquitoes, indicating gut microbiota suppressed <italic>Wolbachia</italic> level after a blood meal (<xref ref-type="bibr" rid="B18">Hughes et al., 2014</xref>).</p>
<p>Operational taxonomic units of <italic>Asaia</italic> (up to 8.9% of the total taxa) were consistently detected in <italic>Wolbachia</italic>-infected <italic>A. stephensi</italic> guts here. The relative abundance of <italic>Asaia</italic> in LBT or LB1 mosquitoes was similar regardless of <italic>Wolbachia</italic> infection, diet switching and rearing ages. This observation differs from the study conducted by <xref ref-type="bibr" rid="B18">Hughes et al. (2014)</xref> where <italic>Asaia</italic> concentration was much higher in blood fed mosquitoes than in non-blood fed ones. The initial assembly of the gut microbiota might account for this discrepancy: <italic>Asaia</italic>, <italic>Pseudomonas</italic>, and unclassified <italic>Gammaproteobacteria</italic> were the most predominant bacteria (<xref ref-type="bibr" rid="B18">Hughes et al., 2014</xref>). <italic>Asaia</italic> is a common commensal of several mosquito species such as <italic>Aedes</italic>, <italic>Anopheles</italic>, and <italic>Culex</italic> (<xref ref-type="bibr" rid="B14">Favia et al., 2008</xref>; <xref ref-type="bibr" rid="B28">Rossi et al., 2015</xref>). Its physiological roles in mosquitoes have been suggested to be nutrient scavenging and regulation of larval development (<xref ref-type="bibr" rid="B9">Chouaia et al., 2012</xref>). Recent studies showed antagonism between <italic>Asaia</italic> and <italic>Wolbachia</italic> in the reproduction organs (such as ovary; <xref ref-type="bibr" rid="B28">Rossi et al., 2015</xref>). Further, antagonistic interactions between <italic>Wolbachia</italic> and the gut microbiota, particularly <italic>Asaia</italic>, resulted in significant fitness costs in the mosquito host (<xref ref-type="bibr" rid="B18">Hughes et al., 2014</xref>). Our observation here showed that <italic>Asaia</italic> was stably maintained in both LB1 and LBT guts, which was consistent with Rossi et al. (<xref ref-type="bibr" rid="B28">Rossi et al., 2015</xref>).</p>
<p>The relative abundance of <italic>Elizabethkingia</italic> OTUs accounted for more than 60% of the total microbial taxa, highlighting their importance in <italic>Anopheles</italic> mosquitoes. Presumably, they contribute to metabolism of sugars in the gut environment, as evidenced by the SusC/SusD-like polysaccharide transport system(s) and glycosidases in their genome (<xref ref-type="bibr" rid="B21">Kukutla et al., 2014</xref>). Supplementation of <italic>Elizabethkingia</italic> cells in the sugar meal for <italic>A. stephensi</italic> led to approximately 50% more egg production compared to those without supplementation, suggesting that <italic>Elizabethkingia</italic> contributed to animal erythrocyte lysis and thus increased host&#x2019;s fecundity in mosquitoes (S. Chen, unpublished). <italic>Elizabethkingia</italic> was frequently detected in water, sediments and insects including field caught, semi-natural reared, and insectary reared mosquitoes, showing that it adapted to diverse ecological niches in nature, but it is commonly found in <italic>Anopheles</italic> midguts (<xref ref-type="bibr" rid="B8">Chen et al., 2015</xref>). Pyrosequencing analysis showed that <italic>Elizabethkingia</italic> OTUs were more abundant in larval <italic>A. gambiae</italic> mosquitoes than in the aquatic medium in which they were being reared, but were present in the latter suggesting a common environmental source (<xref ref-type="bibr" rid="B33">Wang et al., 2011</xref>). Furthermore, <italic>Elizabethkingia anophelis</italic> and its nearly identical taxon, <italic>Elizabethkingia meningoseptica</italic>, were among the most predominant bacteria in adult <italic>A. gambiae</italic> and <italic>A. stephensi</italic> (<xref ref-type="bibr" rid="B33">Wang et al., 2011</xref>; <xref ref-type="bibr" rid="B1">Akhouayri et al., 2013</xref>; <xref ref-type="bibr" rid="B23">Ngwa et al., 2013</xref>).</p>
</sec>
<sec><title>Author Contributions</title>
<p>SC and EW wrote the manuscript. DJ and SC performed the experiments. SC, JZ, ZX, BN, and EW performed the data analysis.</p>
</sec>
<sec><title>Conflict of Interest Statement</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
<p>The reviewer HS and handling Editor declared their shared affiliation, and the handling Editor states that the process nevertheless met the standards of a fair and objective review.</p>
</sec>
</body>
<back>
<ack>
<p>This project was funded by NIH grant R37 AI21884.</p>
</ack>
<sec sec-type="supplementary material">
<title>Supplementary Material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="http://journal.frontiersin.org/article/10.3389/fmicb.2016.01485">http://journal.frontiersin.org/article/10.3389/fmicb.2016.01485</ext-link></p>
<supplementary-material xlink:href="Data_Sheet_1.PDF" id="SM1" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink">
</supplementary-material>
<supplementary-material xlink:href="Data_Sheet_2.XLSX" id="SM2" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink">
</supplementary-material>
</sec>
<ref-list>
<title>References</title>
<ref id="B1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Akhouayri</surname> <given-names>I. G.</given-names></name> <name><surname>Habtewold</surname> <given-names>T.</given-names></name> <name><surname>Christophides</surname> <given-names>G. K.</given-names></name></person-group> (<year>2013</year>). <article-title>Melanotic pathology and vertical transmission of the gut commensal <italic>Elizabethkingia meningoseptica</italic> in the major malaria vector <italic>Anopheles gambiae</italic>.</article-title> <source><italic>PLoS ONE</italic></source> <volume>8</volume>:<issue>e77619</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0077619</pub-id></citation></ref>
<ref id="B2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Benoit</surname> <given-names>J.</given-names></name> <name><surname>Lopez-Martinez</surname> <given-names>G.</given-names></name> <name><surname>Patrick</surname> <given-names>K.</given-names></name> <name><surname>Phillips</surname> <given-names>Z.</given-names></name> <name><surname>Krause</surname> <given-names>T.</given-names></name> <name><surname>Denlinger</surname> <given-names>D.</given-names></name></person-group> (<year>2011</year>). <article-title>Drinking a hot blood meal elicits a protective heat shock response in mosquitoes.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>108</volume> <fpage>8026</fpage>&#x2013;<lpage>8029</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1105195108</pub-id></citation></ref>
<ref id="B3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bian</surname> <given-names>G.</given-names></name> <name><surname>Joshi</surname> <given-names>D.</given-names></name> <name><surname>Dong</surname> <given-names>Y.</given-names></name> <name><surname>Lu</surname> <given-names>P.</given-names></name> <name><surname>Zhou</surname> <given-names>G.</given-names></name> <name><surname>Pan</surname> <given-names>X.</given-names></name></person-group> (<year>2013</year>). <article-title><italic>Wolbachia</italic> invades <italic>Anopheles stephensi</italic> populations and induces refractoriness to <italic>Plasmodium</italic> infection.</article-title> <source><italic>Science</italic></source> <volume>340</volume> <fpage>748</fpage>&#x2013;<lpage>751</lpage>. <pub-id pub-id-type="doi">10.1126/science.1236192</pub-id></citation></ref>
<ref id="B4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Boissi&#x00E8;re</surname> <given-names>A.</given-names></name> <name><surname>Tchioffo</surname> <given-names>M. T.</given-names></name> <name><surname>Bachar</surname> <given-names>D.</given-names></name> <name><surname>Abate</surname> <given-names>L.</given-names></name> <name><surname>Marie</surname> <given-names>A.</given-names></name> <name><surname>Nsango</surname> <given-names>S. E.</given-names></name><etal/></person-group> (<year>2012</year>). <article-title>Midgut microbiota of the malaria mosquito vector <italic>Anopheles gambiae</italic> and interactions with <italic>Plasmodium falciparum</italic> infection.</article-title> <source><italic>PLoS Pathog.</italic></source> <volume>8</volume>:<issue>e1002742</issue>. <pub-id pub-id-type="doi">10.1371/journal.ppat.1002742</pub-id></citation></ref>
<ref id="B5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bongio</surname> <given-names>N. J.</given-names></name> <name><surname>Lampe</surname> <given-names>D. J.</given-names></name></person-group> (<year>2015</year>). <article-title>Inhibition of <italic>Plasmodium</italic> berghei development in mosquitoes by effector proteins secreted from <italic>Asaia</italic> sp. Bacteria using a novel native secretion signal.</article-title> <source><italic>PLoS ONE</italic></source> <volume>10</volume>:<issue>e0143541</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0143541</pub-id></citation></ref>
<ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bourtzis</surname> <given-names>K.</given-names></name> <name><surname>Dobson</surname> <given-names>S. L.</given-names></name> <name><surname>Xi</surname> <given-names>Z.</given-names></name> <name><surname>Rasgon</surname> <given-names>J. L.</given-names></name> <name><surname>Calvitti</surname> <given-names>M.</given-names></name> <name><surname>Moreira</surname> <given-names>L. A.</given-names></name><etal/></person-group> (<year>2014</year>). <article-title>Harnessing mosquito&#x2013;<italic>Wolbachia</italic> symbiosis for vector and disease control.</article-title> <source><italic>Acta Trop.</italic></source> <volume>132(Suppl.)</volume>, <fpage>S150</fpage>&#x2013;<lpage>S163</lpage>. <pub-id pub-id-type="doi">10.1016/j.actatropica.2013.11.004</pub-id></citation></ref>
<ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Capone</surname> <given-names>A.</given-names></name> <name><surname>Ricci</surname> <given-names>I.</given-names></name> <name><surname>Damiani</surname> <given-names>C.</given-names></name> <name><surname>Mosca</surname> <given-names>M.</given-names></name> <name><surname>Rossi</surname> <given-names>P.</given-names></name> <name><surname>Scuppa</surname> <given-names>P.</given-names></name><etal/></person-group> (<year>2013</year>). <article-title>Interactions between <italic>Asaia</italic>, <italic>Plasmodium</italic> and <italic>Anopheles</italic>: new insights into mosquito symbiosis and implications in malaria symbiotic control.</article-title> <source><italic>Parasit. Vectors</italic></source> <volume>6</volume> <fpage>1</fpage>&#x2013;<lpage>13</lpage>. <pub-id pub-id-type="doi">10.1186/1756-3305-6-182</pub-id></citation></ref>
<ref id="B8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname> <given-names>S.</given-names></name> <name><surname>Bagdasarian</surname> <given-names>M.</given-names></name> <name><surname>Walker</surname> <given-names>E. D.</given-names></name></person-group> (<year>2015</year>). <article-title><italic>Elizabethkingia anophelis</italic>: molecular manipulation and interactions with mosquito hosts.</article-title> <source><italic>Appl. Environ. Microbiol.</italic></source> <volume>81</volume> <fpage>2233</fpage>&#x2013;<lpage>2243</lpage>. <pub-id pub-id-type="doi">10.1128/AEM.03733-14</pub-id></citation></ref>
<ref id="B9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chouaia</surname> <given-names>B.</given-names></name> <name><surname>Rossi</surname> <given-names>P.</given-names></name> <name><surname>Epis</surname> <given-names>S.</given-names></name> <name><surname>Mosca</surname> <given-names>M.</given-names></name> <name><surname>Ricci</surname> <given-names>I.</given-names></name> <name><surname>Damiani</surname> <given-names>C.</given-names></name><etal/></person-group> (<year>2012</year>). <article-title>Delayed larval development in <italic>Anopheles</italic> mosquitoes deprived of <italic>Asaia</italic> bacterial symbionts.</article-title> <source><italic>BMC Microbiol.</italic></source> <volume>12(Suppl. 1)</volume>:<issue>S2</issue>. <pub-id pub-id-type="doi">10.1186/1471-2180-1112-S1181-S1182</pub-id></citation></ref>
<ref id="B10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cirimotich</surname> <given-names>C. M.</given-names></name> <name><surname>Dong</surname> <given-names>Y.</given-names></name> <name><surname>Clayton</surname> <given-names>A. M.</given-names></name> <name><surname>Sandiford</surname> <given-names>S. L.</given-names></name> <name><surname>Souza-Neto</surname> <given-names>J. A.</given-names></name> <name><surname>Mulenga</surname> <given-names>M.</given-names></name><etal/></person-group> (<year>2011</year>). <article-title>Natural microbe-mediated refractoriness to <italic>Plasmodium</italic> infection in <italic>Anopheles gambiae</italic>.</article-title> <source><italic>Science</italic></source> <volume>332</volume> <fpage>855</fpage>&#x2013;<lpage>858</lpage>. <pub-id pub-id-type="doi">10.1126/science.1201618</pub-id></citation></ref>
<ref id="B11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Coon</surname> <given-names>K. L.</given-names></name> <name><surname>Vogel</surname> <given-names>K. J.</given-names></name> <name><surname>Brown</surname> <given-names>M. R.</given-names></name> <name><surname>Strand</surname> <given-names>M. R.</given-names></name></person-group> (<year>2014</year>). <article-title>Mosquitoes rely on their gut microbiota for development.</article-title> <source><italic>Mol. Ecol.</italic></source> <volume>23</volume> <fpage>2727</fpage>&#x2013;<lpage>2739</lpage>. <pub-id pub-id-type="doi">10.1111/mec.12771</pub-id></citation></ref>
<ref id="B12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Correa</surname> <given-names>C. C.</given-names></name> <name><surname>Ballard</surname> <given-names>J. W. O.</given-names></name></person-group> (<year>2014</year>). <article-title>What can symbiont titres tell us about co-evolution of <italic>Wolbachia</italic> and their host?</article-title> <source><italic>J. Invertebr. Pathol.</italic></source> <volume>118</volume> <fpage>20</fpage>&#x2013;<lpage>27</lpage>. <pub-id pub-id-type="doi">10.1016/j.jip.2014.02.009</pub-id></citation></ref>
<ref id="B13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dobson</surname> <given-names>S.</given-names></name> <name><surname>Bourtzis</surname> <given-names>K.</given-names></name> <name><surname>Braig</surname> <given-names>H.</given-names></name> <name><surname>Jones</surname> <given-names>B.</given-names></name> <name><surname>Zhou</surname> <given-names>W.</given-names></name> <name><surname>Rousset</surname> <given-names>F.</given-names></name><etal/></person-group> (<year>1999</year>). <article-title><italic>Wolbachia</italic> infections are distributed throughout insect somatic and germ line tissues.</article-title> <source><italic>Insect Biochem. Mol. Biol.</italic></source> <volume>29</volume> <fpage>153</fpage>&#x2013;<lpage>160</lpage>. <pub-id pub-id-type="doi">10.1016/S0965-1748(98)00119-2</pub-id></citation></ref>
<ref id="B14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Favia</surname> <given-names>G.</given-names></name> <name><surname>Ricci</surname> <given-names>I.</given-names></name> <name><surname>Marzorati</surname> <given-names>M.</given-names></name> <name><surname>Negri</surname> <given-names>I.</given-names></name> <name><surname>Alma</surname> <given-names>A.</given-names></name> <name><surname>Sacchi</surname> <given-names>L.</given-names></name></person-group> (<year>2008</year>). <article-title>Bacteria of the genus <italic>Asaia</italic>: a potential weapon against malaria.</article-title> <source><italic>Adv. Exp. Med. Biol.</italic></source> <volume>627</volume> <fpage>49</fpage>&#x2013;<lpage>59</lpage>. <pub-id pub-id-type="doi">10.1007/978-0-387-78225-6_4</pub-id></citation></ref>
<ref id="B15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Federici</surname> <given-names>B. A.</given-names></name> <name><surname>Park</surname> <given-names>H.-W.</given-names></name> <name><surname>Bideshi</surname> <given-names>D. K.</given-names></name> <name><surname>Wirth</surname> <given-names>M. C.</given-names></name> <name><surname>Johnson</surname> <given-names>J. J.</given-names></name></person-group> (<year>2003</year>). <article-title>Recombinant bacteria for mosquito control.</article-title> <source><italic>J. Exp. Biol.</italic></source> <volume>206</volume> <fpage>3877</fpage>&#x2013;<lpage>3885</lpage>. <pub-id pub-id-type="doi">10.1242/jeb.00643</pub-id></citation></ref>
<ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gimonneau</surname> <given-names>G.</given-names></name> <name><surname>Tchioffo</surname> <given-names>M. T.</given-names></name> <name><surname>Abate</surname> <given-names>L.</given-names></name> <name><surname>Boissi&#x00E8;re</surname> <given-names>A.</given-names></name> <name><surname>Awono-Amb&#x00E9;n&#x00E9;</surname> <given-names>P. H.</given-names></name> <name><surname>Nsango</surname> <given-names>S. E.</given-names></name><etal/></person-group> (<year>2014</year>). <article-title>Composition of <italic>Anopheles coluzzii</italic> and <italic>Anopheles gambiae</italic> microbiota from larval to adult stages.</article-title> <source><italic>Infect. Genet. Evol.</italic></source> <volume>28</volume> <fpage>715</fpage>&#x2013;<lpage>724</lpage>. <pub-id pub-id-type="doi">10.1016/j.meegid.2014.09.029</pub-id></citation></ref>
<ref id="B17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hilgenboecker</surname> <given-names>K.</given-names></name> <name><surname>Hammerstein</surname> <given-names>P.</given-names></name> <name><surname>Schlattmann</surname> <given-names>P.</given-names></name> <name><surname>Telschow</surname> <given-names>A.</given-names></name> <name><surname>Werren</surname> <given-names>J. H.</given-names></name></person-group> (<year>2008</year>). <article-title>How many species are infected with <italic>Wolbachia</italic>?&#x2013;a statistical analysis of current data.</article-title> <source><italic>FEMS Microbiol. Lett.</italic></source> <volume>281</volume> <fpage>215</fpage>&#x2013;<lpage>220</lpage>. <pub-id pub-id-type="doi">10.1111/j.1574-6968.2008.01110.x</pub-id></citation></ref>
<ref id="B18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hughes</surname> <given-names>G. L.</given-names></name> <name><surname>Dodson</surname> <given-names>B. L.</given-names></name> <name><surname>Johnson</surname> <given-names>R. M.</given-names></name> <name><surname>Murdock</surname> <given-names>C. C.</given-names></name> <name><surname>Tsuijimoto</surname> <given-names>H.</given-names></name> <name><surname>Suzuki</surname> <given-names>Y.</given-names></name></person-group> (<year>2014</year>). <article-title>Native microbiome impedes vertical transmission of <italic>Wolbachia</italic> in <italic>Anopheles</italic> mosquitoes.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>111</volume> <fpage>12498</fpage>&#x2013;<lpage>12503</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1408888111</pub-id></citation></ref>
<ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Iturbe-Ormaetxe</surname> <given-names>I.</given-names></name> <name><surname>Walker</surname> <given-names>T.</given-names></name> <name><surname>O&#x2019; Neill</surname> <given-names>S. L.</given-names></name></person-group> (<year>2011</year>). <article-title><italic>Wolbachia</italic> and the biological control of mosquito-borne disease.</article-title> <source><italic>EMBO Rep.</italic></source> <volume>12</volume> <fpage>508</fpage>&#x2013;<lpage>518</lpage>. <pub-id pub-id-type="doi">10.1038/embor.2011.84</pub-id></citation></ref>
<ref id="B20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Joshi</surname> <given-names>D.</given-names></name> <name><surname>Mc Fadden</surname> <given-names>M. J.</given-names></name> <name><surname>Bevins</surname> <given-names>D.</given-names></name> <name><surname>Zhang</surname> <given-names>F.</given-names></name> <name><surname>Xi</surname> <given-names>Z.</given-names></name></person-group> (<year>2014</year>). <article-title><italic>Wolbachia</italic> strain wAlbB confers both fitness costs and benefit on <italic>Anopheles stephensi</italic>.</article-title> <source><italic>Parasit. Vector</italic></source> <volume>7</volume>:<issue>336</issue>. <pub-id pub-id-type="doi">10.1186/1756-3305-7-336</pub-id></citation></ref>
<ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kukutla</surname> <given-names>P.</given-names></name> <name><surname>Lindberg</surname> <given-names>B. G.</given-names></name> <name><surname>Pei</surname> <given-names>D.</given-names></name> <name><surname>Rayl</surname> <given-names>M.</given-names></name> <name><surname>Yu</surname> <given-names>W.</given-names></name> <name><surname>Steritz</surname> <given-names>M.</given-names></name><etal/></person-group> (<year>2014</year>). <article-title>Insights from the genome annotation of <italic>Elizabethkingia anophelis</italic> from the malaria vector <italic>Anopheles gambiae</italic>.</article-title> <source><italic>PLoS ONE</italic></source> <volume>9</volume>:<issue>e97715</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0097715</pub-id></citation></ref>
<ref id="B22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Luckhart</surname> <given-names>S.</given-names></name> <name><surname>Vodovotz</surname> <given-names>Y.</given-names></name> <name><surname>Cui</surname> <given-names>L.</given-names></name> <name><surname>Rosenberg</surname> <given-names>R.</given-names></name></person-group> (<year>1998</year>). <article-title>The mosquito <italic>Anopheles stephensi</italic> limits malaria parasite development with inducible synthesis of nitric oxide.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>95</volume> <fpage>5700</fpage>&#x2013;<lpage>5705</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.95.10.5700</pub-id></citation></ref>
<ref id="B23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ngwa</surname> <given-names>C. J.</given-names></name> <name><surname>Gl&#x00F6;ckner</surname> <given-names>V.</given-names></name> <name><surname>Abdelmohsen</surname> <given-names>U. R.</given-names></name> <name><surname>Scheuermayer</surname> <given-names>M.</given-names></name> <name><surname>Fischer</surname> <given-names>R.</given-names></name> <name><surname>Hentschel</surname> <given-names>U.</given-names></name><etal/></person-group> (<year>2013</year>). <article-title>16S rRNA gene-based identification of <italic>Elizabethkingia meningoseptica</italic> (Flavobacteriales: Flavobacteriaceae) as a dominant midgut bacterium of the Asian malaria vector <italic>Anopheles stephensi</italic> (Diphteria: Culicidae) with antimicrobial activities.</article-title> <source><italic>J. Med. Entomol.</italic></source> <volume>50</volume> <fpage>404</fpage>&#x2013;<lpage>414</lpage>. <pub-id pub-id-type="doi">10.1603/ME12180</pub-id></citation></ref>
<ref id="B24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Oksanen</surname> <given-names>J.</given-names></name> <name><surname>Blanchet</surname> <given-names>F. G.</given-names></name> <name><surname>Kindt</surname> <given-names>R.</given-names></name> <name><surname>Legendre</surname> <given-names>P.</given-names></name> <name><surname>Minchin</surname> <given-names>P. R.</given-names></name> <name><surname>O&#x2019;Hara</surname> <given-names>R. B.</given-names></name></person-group> (<year>2015</year>). <source><italic>Vegan: Community Ecology</italic>.</source> Available at: <ext-link ext-link-type="uri" xlink:href="http://cran.r-project.org/package=vegan">http://cran.r-project.org/package=vegan</ext-link></citation></ref>
<ref id="B25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pan</surname> <given-names>X.</given-names></name> <name><surname>Zhou</surname> <given-names>G.</given-names></name> <name><surname>Wu</surname> <given-names>J.</given-names></name> <name><surname>Bian</surname> <given-names>G.</given-names></name> <name><surname>Lu</surname> <given-names>P.</given-names></name> <name><surname>Raikhel</surname> <given-names>A.</given-names></name><etal/></person-group> (<year>2012</year>). <article-title><italic>Wolbachia</italic> induces reactive oxygen species (ROS)-dependent activation of the Toll pathway to control dengue virus in the mosquito <italic>Aedes aegypti</italic>.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>109</volume> <fpage>E23</fpage>&#x2013;<lpage>E31</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1116932108</pub-id></citation></ref>
<ref id="B26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Paradis</surname> <given-names>E.</given-names></name> <name><surname>Claude</surname> <given-names>J.</given-names></name> <name><surname>Strimmer</surname> <given-names>K.</given-names></name></person-group> (<year>2004</year>). <article-title>APE: analyses of phylogenetics and evolution in R language.</article-title> <source><italic>Bioinformatics</italic></source> <volume>20</volume> <fpage>289</fpage>&#x2013;<lpage>290</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btg412</pub-id></citation></ref>
<ref id="B27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Peterson</surname> <given-names>T. M. L.</given-names></name> <name><surname>Luckhart</surname> <given-names>S.</given-names></name></person-group> (<year>2006</year>). <article-title>A mosquito 2-Cys peroxiredoxin protects against nitrosative and oxidative stresses associated with malaria parasite infection.</article-title> <source><italic>Free Radic. Biol. Med.</italic></source> <volume>40</volume> <fpage>1067</fpage>&#x2013;<lpage>1082</lpage>. <pub-id pub-id-type="doi">10.1016/j.freeradbiomed.2005.10.059</pub-id></citation></ref>
<ref id="B28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rossi</surname> <given-names>P.</given-names></name> <name><surname>Ricci</surname> <given-names>I.</given-names></name> <name><surname>Cappelli</surname> <given-names>A.</given-names></name> <name><surname>Damiani</surname> <given-names>C.</given-names></name> <name><surname>Ulissi</surname> <given-names>U.</given-names></name> <name><surname>Mancini</surname> <given-names>M. V.</given-names></name><etal/></person-group> (<year>2015</year>). <article-title>Mutual exclusion of Asaia and <italic>Wolbachia</italic> in the reproductive organs of mosquito vectors.</article-title> <source><italic>Parasit. Vector</italic></source> <volume>8</volume> <fpage>1</fpage>&#x2013;<lpage>10</lpage>. <pub-id pub-id-type="doi">10.1186/s13071-015-0888-0</pub-id></citation></ref>
<ref id="B29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Schloss</surname> <given-names>P. D.</given-names></name> <name><surname>Westcott</surname> <given-names>S. L.</given-names></name> <name><surname>Ryabin</surname> <given-names>T.</given-names></name> <name><surname>Hall</surname> <given-names>J. R.</given-names></name> <name><surname>Hartmann</surname> <given-names>M.</given-names></name> <name><surname>Hollister</surname> <given-names>E. B.</given-names></name><etal/></person-group> (<year>2009</year>). <article-title>Introducing mothur: open-source, platform-independent, community supported software for describing and comparing microbial communities.</article-title> <source><italic>Appl. Environ. Microbiol.</italic></source> <volume>75</volume> <fpage>7537</fpage>&#x2013;<lpage>7541</lpage>. <pub-id pub-id-type="doi">10.1128/AEM.01541-09</pub-id></citation></ref>
<ref id="B30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sinkins</surname> <given-names>S. P.</given-names></name></person-group> (<year>2004</year>). <article-title><italic>Wolbachia</italic> and cytoplasmic incompatibility in mosquitoes.</article-title> <source><italic>Insect Biochem. Mol. Biol.</italic></source> <volume>34</volume> <fpage>723</fpage>&#x2013;<lpage>729</lpage>. <pub-id pub-id-type="doi">10.1016/j.ibmb.2004.03.025</pub-id></citation></ref>
<ref id="B31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Tchioffo</surname> <given-names>M. T.</given-names></name> <name><surname>Boissi&#x00E8;re</surname> <given-names>A.</given-names></name> <name><surname>Abate</surname> <given-names>L.</given-names></name> <name><surname>Nsango</surname> <given-names>S. E.</given-names></name> <name><surname>Bayib&#x00E9;ki</surname> <given-names>A. N.</given-names></name> <name><surname>Awono-Amb&#x00E9;n&#x00E9;</surname> <given-names>P. H.</given-names></name><etal/></person-group> (<year>2016</year>). <article-title>Dynamics of bacterial community composition in the malaria mosquito&#x2019;s epithelia.</article-title> <source><italic>Front. Microbiol.</italic></source> <volume>6</volume>:<issue>1500</issue>. <pub-id pub-id-type="doi">10.3389/fmicb.2015.01500</pub-id></citation></ref>
<ref id="B32"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>S.</given-names></name> <name><surname>Ghosh</surname> <given-names>A. K.</given-names></name> <name><surname>Bongio</surname> <given-names>N.</given-names></name> <name><surname>Stebbings</surname> <given-names>K. A.</given-names></name> <name><surname>Lampe</surname> <given-names>D. J.</given-names></name> <name><surname>Jacobs-Lorena</surname> <given-names>M.</given-names></name></person-group> (<year>2012</year>). <article-title>Fighting malaria with engineered symbiotic bacteria from vector mosquitoes.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>109</volume> <fpage>12734</fpage>&#x2013;<lpage>12739</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1204158109</pub-id></citation></ref>
<ref id="B33"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>Y.</given-names></name> <name><surname>Gilbreath</surname> <given-names>T. M.</given-names></name> <name><surname>Kukutla</surname> <given-names>P.</given-names></name> <name><surname>Yan</surname> <given-names>G.</given-names></name> <name><surname>Xu</surname> <given-names>J.</given-names></name></person-group> (<year>2011</year>). <article-title>Dynamic gut microbiome across life history of the malaria mosquito <italic>Anopheles gambiae</italic> in Kenya.</article-title> <source><italic>PLoS ONE</italic></source> <volume>6</volume>:<issue>e24767</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0024767</pub-id></citation></ref>
<ref id="B34"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Werren</surname> <given-names>J. H.</given-names></name> <name><surname>Baldo</surname> <given-names>L.</given-names></name> <name><surname>Clark</surname> <given-names>M. E.</given-names></name></person-group> (<year>2008</year>). <article-title><italic>Wolbachia</italic>: master manipulators of invertebrate biology.</article-title> <source><italic>Nat. Rev. Microbiol.</italic></source> <volume>6</volume> <fpage>741</fpage>&#x2013;<lpage>751</lpage>. <pub-id pub-id-type="doi">10.1038/nrmicro1969</pub-id></citation></ref>
<ref id="B35"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xi</surname> <given-names>Z.</given-names></name> <name><surname>Khoo</surname> <given-names>C. C.</given-names></name> <name><surname>Dobson</surname> <given-names>S. L.</given-names></name></person-group> (<year>2005</year>). <article-title>Wolbachia establishment and invasion in an <italic>Aedes aegypti</italic> laboratory population.</article-title> <source><italic>Science</italic></source> <volume>310</volume> <fpage>326</fpage>&#x2013;<lpage>328</lpage>. <pub-id pub-id-type="doi">10.1126/science.1117607</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn id="fn01"><label>1</label><p><ext-link ext-link-type="uri" xlink:href="http://www.mothur.org">http://www.mothur.org</ext-link></p></fn>
</fn-group>
</back>
</article>