# Bacterial communities associated with host-adapted populations of pea aphids revealed by deep sequencing of 16S ribosomal DNA.

Jean-Pierre Gauthier, Yannick Outreman, Lucie Mieuzet, Jean-Christophe Simon

► To cite this version:

Jean-Pierre Gauthier, Yannick Outreman, Lucie Mieuzet, Jean-Christophe Simon. Bacterial communities associated with host-adapted populations of pea aphids revealed by deep sequencing of 16S ribosomal DNA. PLoS ONE, 2015, 10 (3), pp.e0120664. 10.1371/journal.pone.0120664. hal-01274516v2

HAL Id: hal-01274516

https://hal.science/hal-01274516v2

Submitted on 1 Mar 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Bacterial Communities Associated with Host-Adapted Populations of Pea Aphids Revealed by Deep Sequencing of 16S Ribosomal DNA

Jean-Pierre Gauthier¹, Yannick Outreman², Lucie Mieuzet¹, Jean-Christophe Simon²*

1 INRA, UMR 1349 IGEPP "Institut de Génétique, Environnement et Protection des Plantes", 35653, Le Rheu, France, 2 Agrocampus Ouest, UMR 1349 IGEPP "Institut de Génétique, Environnement et Protection des Plantes", 35042, Rennes, France

* jean-christophe.simon@rennes.inra.fr

## Abstract

Associations between microbes and animals are ubiquitous and hosts may benefit from har-
bouring microbial communities through improved resource exploitation or resistance to en-
vironmental stress. The pea aphid, *Acyrthosiphon pisum*, is the host of heritable bacterial
symbionts, including the obligate endosymbiont *Buchnera aphidicola* and several faculta-
tive symbionts. While obligate symbionts supply aphids with key nutrients, facultative sym-
bionts influence their hosts in many ways such as protection against natural enemies, heat
tolerance, color change and reproduction alteration. The pea aphid also encompasses mul-
tiple plant-specialized biotypes, each adapted to one or a few legume species. Facultative
symbiont communities differ strongly between biotypes, although bacterial involvement in
plant specialization is uncertain. Here, we analyse the diversity of bacterial communities as-
sociated with nine biotypes of the pea aphid complex using amplicon pyrosequencing of
16S rRNA genes. Combined clustering and phylogenetic analyses of 16S sequences al-
lowed identifying 21 bacterial OTUs (Operational Taxonomic Unit). More than 98% of the
sequencing reads were assigned to known pea aphid symbionts. The presence of *Wolba-
chia* was confirmed in *A. pisum* while *Erwinia* and *Pantoea*, two gut associates, were de-
tected in multiple samples. The diversity of bacterial communities harboured by pea aphid
biotypes was very low, ranging from 3 to 11 OTUs across samples. Bacterial communities
differed more between than within biotypes but this difference did not correlate with the ge-
netic divergence between biotypes. Altogether, these results confirm that the aphid micro-
biota is dominated by a few heritable symbionts and that plant specialization is an important
structuring factor of bacterial communities associated with the pea aphid complex. Howev-
er, since we examined the microbiota of aphid samples kept a few generations in controlled
conditions, it may be that bacterial diversity was underestimated due to the possible loss of
environmental or transient taxa.

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Gauthier J-P, Outreman Y, Mieuzet L, Simon J-C (2015) Bacterial Communities Associated with Host-Adapted Populations of Pea Aphids Revealed by Deep Sequencing of 16S Ribosomal DNA. PLoS ONE 10(3): e0120664. doi:10.1371/journal.pone.0120664

**Academic Editor:** Sebastien Duperron, Université Pierre et Marie Curie, FRANCE

**Received:** September 9, 2014

**Accepted:** February 5, 2015

**Published:** March 25, 2015

**Copyright:** © 2015 Gauthier et al. This is an open access article distributed under the terms of the [Creative Commons Attribution License](http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** The dataset including 21 unique bacterial phylotypes have been deposited in Genbank database (accession numbers KP265043-KP265116).

**Funding:** This work was supported by funds from the French Research Agency under the ANR-11-BSV7-005-01 "Speciaphid" to JCS, and from the INRA Département Santé des Plantes et Environnement and Université Rennes 1 to JCS and YO.

**Competing Interests:** The authors have declared that no competing interests exist.
**Introduction**

Sustained associations between microbes and eukaryotic hosts are widespread. This is particu-
larly true for insects which are engaged in multiple forms of interactions with micro-organisms
referred to as symbiosis in its broad sense (i.e. all types of close associations between organisms
of different species). Many arthropods have established intimate and mutualistic partnerships
with obligate symbionts which generally supplement the host's diet in key nutrients [1], [2]. In
addition, researches over the last decades have highlighted the taxonomic diversity and the
many roles of facultative symbionts. Although not essential for insect reproduction and surviv-
al, facultative symbionts may have indeed profound and sometimes spectacular effects on host
ecology, physiology and behavior [3], [4]. However, until recently, the exploration of insect
symbiont diversity was limited by the use of taxon-specific detection techniques or DNA-
limiting, time consuming and low throughput methods. The development of Next Generation
Sequencing techniques (NGS) now permits to assess without *a priori* the diversity and struc-
ture of insect microbial communities to better understand both the influence of microbes on
their host populations and the dynamics of symbiotic associations. Although it is now possible
to characterize and analyze the whole microbial community of a given host, these techniques
have been mainly restricted so far to insect model species or species of medical importance
[5], [6], [7], [8].

Aphids represent a well-studied case of symbiotic associations. In addition to their obligate
symbiont (*Buchnera aphidicola*), these insects may harbour one or several heritable facultative
bacterial symbionts. These symbionts have been mostly studied in the pea aphid, *Acyrthosi-
phon pisum*, a pest of legume crops. So far, eight bacterial facultative symbiont taxa of the pea
aphid are known: the five Gammaproteobacteria *Hamiltonella defensa*, *Regiella insecticola*, *Ser-
ratia symbiotica* [9], *Rickettsiella viridis* [10, 11] and PAXS (Pea Aphid X-type Symbiont), [12];
two Alphaproteobacteria of the genera *Rickettsia* [13] and *Wolbachia* [14]; and finally *Spiro-
plasma* [15] belonging to Mollicutes. In pea aphids, facultative symbionts reside in different
parts of the host including the hemolymph, sheath cells associated with the primary endosym-
biont bacteriome or within bacteriocytes themselves [9]. These facultative symbionts may ben-
efit their host through increased performance on specific plants [16], body colour change [11],
heat tolerance [17] or protection against natural enemies [18], [19]. Although transmitted ver-
tically with high fidelity, these symbionts may occasionally move horizontally within and be-
tween species [3], [20]. Besides heritable symbionts, little is known about gut associates, the
bacteria colonizing the aphid digestive tract. More generally, we lack a deep characterization of
the pea aphid microbiome and a better assessment of changes in bacterial communities accord-
ing to environmental factors and host genotypes. Here, we used deep 454 amplicon pyrose-
quencing of 16S rDNA genes, to analyze the diversity and structure of bacterial communities
associated with pea aphid populations specialized on nine different host plants.

The pea aphid consists of a series of host-adapted biotypes specialized on different host plants and showing a continuum of divergence from host races still exchanging some genes to genetically isolated incipient species. So far, 11 biotypes have been described, each adapted to one or a few legume species [21] and differing in their symbiotic complement [20], [22], [23].

Here, we examined the diversity of bacterial communities associated with nine biotypes of
the pea aphid complex through combined clustering and phylogenetic analyses of 16S rDNA
sequences. In particular, we were looking for unreported symbionts in the pea aphid and other
bacterial taxa such as gut associates present in low abundance. We also assessed the importance
of plant specialization on the structure of bacterial communities by comparing the micro-
biomes of the different pea aphid biotypes. Finally, we tested whether genetic divergence be-
tween aphid biotypes correlated with dissimilarity between their respective bacterial
<table><caption>Table 1. Host plant and geographic origins of pea aphid samples examined for intra-host bacterial diversity with 16S rDNA amplicon sequencing.</caption><thead><tr><th>Host plant</th><th>Common name</th><th>Sample</th><th>Collection date</th><th>Site</th><th>GPS position</th></tr></thead><tbody><tr><td>Cytisus scoparius</td><td>Common broom</td><td>Csco</td><td>August 2011</td><td>Lantenay</td><td>46°03'31" N, 05°32'32" E</td></tr><tr><td>Lathyrus pratensis</td><td>Meadow vetchling</td><td>Lpra</td><td>August 2011</td><td>Lantenay</td><td>46°03'31" N, 05°32'32" E</td></tr><tr><td>Medicago lupulina</td><td>Black medic</td><td>Mlup</td><td>August 2011</td><td>Lantenay</td><td>46°03'31" N, 05°32'32" E</td></tr><tr><td>Medicago sativa</td><td>Alfalfa</td><td>Msat-1</td><td>May 2011</td><td>Domagné</td><td>48°04'15" N, 01°23'34" W</td></tr><tr><td>Medicago sativa</td><td>Alfalfa</td><td>Msat-2</td><td>May 2011</td><td>Domagné</td><td>48°04'15" N, 01°23'34" W</td></tr><tr><td>Melilotus albus</td><td>White melilot</td><td>Malb</td><td>August 2011</td><td>Lantenay</td><td>46°03'31" N, 05°32'32" E</td></tr><tr><td>Ononis spinosa</td><td>Spiny restharrow</td><td>Ospi</td><td>August 2011</td><td>Lantenay</td><td>46°03'31" N, 05°32'32" E</td></tr><tr><td>Pisum sativum</td><td>Pea</td><td>Psat-1</td><td>May 2011</td><td>Domagné</td><td>48°04'15" N, 01°23'34" W</td></tr><tr><td>Pisum sativum</td><td>Pea</td><td>Psat-2</td><td>May 2011</td><td>La Chapelle Thouarault</td><td>48°07'27" N, 01°51'58" W</td></tr><tr><td>Securigera varia</td><td>Crown vetch</td><td>Svar</td><td>August 2011</td><td>Lantenay</td><td>46°03'31" N, 05°32'32" E</td></tr><tr><td>Trifolium pratense</td><td>Red clover</td><td>Tprat-1</td><td>May 2011</td><td>Domagné</td><td>48°04'15" N, 01°23'34" W</td></tr><tr><td>Trifolium pratense</td><td>Red clover</td><td>Tprat-2</td><td>May 2011</td><td>Domagné</td><td>48°04'15" N, 01°23'34" W</td></tr></tbody></table>

doi:10.1371/journal.pone.0120664.t001

communities. In this study, we wanted to avoid the detection of microbes associated with internal parasites such as parasitic wasps (i.e. aphid hymenopteran parasitoids which are frequently encountered in field populations of aphids as “mummies”) and to focus our study on bacteria with prolonged associations with their hosts. Therefore, we analyzed the microbiota of aphid samples that were kept a few generations in controlled conditions, being however aware this could alter microbial diversity associated with pea aphids in the field.

# Material and Methods

## Pea aphid samples

Aphids were sampled from nine different plant species on which *A. pisum* is known to occur as host-specialized populations or biotypes [21]: *Cytisus scoparius* (common broom), *Lathyrus pratensis* (meadow vetchling), *Medicago lupulina* (black medic), *Medicago sativa* (alfalfa), *Melilotus albus* (white melilot), *Ononis spinosa* (spiny restharrow), *Pisum sativum* (pea), *Securigera varia* (crown vetch), and *Trifolium pratense* (red clover). Aphid collections took place in spring and summer 2011 in eastern and western France and consisted exclusively of wingless parthenogenetic females (Table 1). No specific permissions were required for these collection activities and samplings did not involve endangered or protected species. To prevent the detection in our sequencing analysis of bacterial taxa associated with internal parasites, field-collected aphids were grown individually on broad bean (the universal host for all pea aphid biotypes [23]) in controlled conditions ensuring clonal reproduction (16 hours light per day, 18°C). Parasitism rates ranged from 0% to 86% across pea aphid samples and were mainly due to parasitic wasps (e.g. *Aphidius ervi*) and fungal pathogens (e.g. *Pandora neoaphidis*). Healthy parthenogenetic females were let to reproduce for two generations and those that founded clonal lineages were kept alive for subsequent molecular analyses. Each clonal lineage was characterized at several microsatellite markers to assign genotype to biotype using a procedure detailed in [24]. Based on genotyping and clustering analyses, we then discarded all presumable copies of the same clone, the migrant individuals coming from other host races, and the hybrids between biotypes. Table 1 describes the different pea aphid samples considered in the present study.
## DNA extraction and pyrosequencing of 16S rDNA amplicons

Adults of *A. pisum* from the distinct genotypes and biotypes were first surface-sterilized with 70% ethanol for 1 min, 10% bleach for 1 min and three washes of ultrapure water for 1 min. DNA was extracted from each individual (whole body) using the QIAGEN DNeasy kit following the manufacturer's protocol. DNA extractions were then quantified using a Nanodrop spectrophotometer and DNA from 20 individuals (each with distinct genotypes but from the same biotype) were normalized and combined to create DNA pools. Twelve pools were constituted in total: two from alfalfa, clover and pea from each of the two sampled fields and one from the six other plant origins (Table 1). We used universal primers to amplify a V4-V5 variable region of the 16S ribosomal DNA (rDNA) in Eubacteria and Archebacteria that contains sequence variation to distinguish bacterial species. The 16S primers (F: 5' GTGCCAGCM5 GCCGCGGTAATAC 3' and R: 5' CCGTCAATTCCTTTGAGTTT 3' from) are highly conserved in bacteria and amplify a region of about 420 bp. Additional sequences were added to the 5' end of the primers for multiplexing of the samples and for Roche Titanium amplicon sequencing. Each fusion primer consisted of Adaptor A (5' CCATCTCATCCCTGCGTGTCTCC GAC 3') or Adaptor B (5' CCTATCCCCTGTGTGCCTTGGCAGTC 3') followed by a 4-mer key sequence (5' TCAG 3'), a 5-mer to 10-mer Multiplex IDentifier (MID), and finally the 16S rDNA primer. In total, we used 12 different MIDs for both the forward and the reverse primers. The primers were HPLC purified and PCR amplification was performed in a volume of 25 µL with High Fidelity DNA polymerase (Roche). After initial denaturation for 2 min at 94°C, 33 cycles of denaturation for 30 sec at 94°C, annealing for 30 sec at 64°C and elongation at 72°C for 1 min were performed. For the last cycle, the elongation time was extended to 6 min. Polymerase Chain Reactions were duplicated for each pooled sample, each PCR product duplicate being then cleaned using the QIAquick PCR purification kit (Qiagen) following the manufacturer's instructions and sequenced in both directions on half of a Roche 454 FLX Genome Sequencer plate using Titanium Series reagents at the Environmental Genomics platform of BioGenouest, University of Rennes 1, France.

## 16S rDNA sequence analysis

After application of filter criteria through manufacturer's pipeline, sequencing on the GS FLX machine generated 425,410 sequences with a median size of 405 bp. Sequences were then processed using the workflow provided by MOTHUR [25]. During the demultiplexing step, all sequences with more than one mismatch with our MIDs were removed. To be conservative in our approach and to increase reliability in taxonomic assignment, we then kept only sequences with a length between 400 and 430 bp (expected size of the PCR product is about 420 bp). Finally, chimeric sequences that arise during PCR were eliminated using the UCHIME algorithm [26].

After these selection steps, sequences were aligned using the Greengenes Core Set alignment as a template [27]. The aligned sequences were then assigned to Operational Taxonomic Units (OTUs) based on the general approach implemented in MOTHUR which finds the closest Greengenes template for each read. OTUs that were present in only one of the two PCR replicates of each aphid sample were not considered in further analyses. Representative sequences of each OTU found in the different samples were then uploaded in a file that also contained reference 16S rDNA sequences of common bacterial symbionts of arthropods (Table 2). OTU and reference sequences were used to generate a Neighbor-Joining tree to assess the phylogenetic position of bacterial taxa present in pea aphid samples. The NJ tree was constructed using Clustal W version 2.0.11 [28] and visualized with MEGA 5 [29]. Genbank accession numbers of representative sequences of each OTU found in the different samples (74 sequences in total) and used to build the tree are given in S1 Table.
<table><caption>Table 2. Reference sequences of 16S rDNA genes from bacterial taxa associated with insect hosts (except *Spiroplasma* from *Citrus*).</caption><thead><tr><th>Bacterial taxa</th><th>Host</th><th>Accession numbers</th></tr></thead><tbody><tr><td>Arsenophonus sp.</td><td>Stomaphis fagi</td><td>gi|242397733|gb|FJ655540.1|</td></tr><tr><td>Arsenophonus sp.</td><td>Bemisia tabaci</td><td>gi|354508545|gb|JN204477.1|</td></tr><tr><td>Buchnera aphidicola (ref.1)</td><td>Acyrthosiphon pisum</td><td>gi|304030|gb|M27039.1|</td></tr><tr><td>Buchnera aphidicola (ref.2)</td><td>Acyrthosiphon pisum</td><td>gi|9588077|dbj|AB033774.1|</td></tr><tr><td>Cardinium sp.</td><td>Bemisia tabaci</td><td>gi|354508549|gb|JN204481.1|</td></tr><tr><td>Cardinium sp.</td><td>Encarsia hispida</td><td>gi|114054953|gb|DQ854695.1|</td></tr><tr><td>Erwinia aphidicola (ref.1)</td><td>Acyrthosiphon pisum</td><td>gi|340051242|emb|FN547376.1|</td></tr><tr><td>Erwinia aphidicola (ref.2)</td><td>Acyrthosiphon pisum</td><td>gi|359805495|dbj|AB681773.1|</td></tr><tr><td>Erwinia aphidicola (ref.3)</td><td>Acyrthosiphon pisum</td><td>gi|157073730|dbj|AB273744.1|</td></tr><tr><td>Hamiltonella defensa</td><td>Acyrthosiphon pisum</td><td>gi|56785712|gb|AY692361.1|</td></tr><tr><td>Pantoea agglomerans</td><td>Acyrthosiphon pisum</td><td>gi|57864392|dbj|AB004757.2|</td></tr><tr><td>Pantoea agglomerans</td><td>Myzus persicae</td><td>gi|57670121|gb|AY849936.1|</td></tr><tr><td>PAXS</td><td>Acyrthosiphon pisum</td><td>gi|226246999|gb|FJ821502.1|</td></tr><tr><td>Regiella insecticola (ref.1)</td><td>Acyrthosiphon pisum</td><td>gi|10567782|gb|AF293618.1|</td></tr><tr><td>Regiella insecticola (ref.2)</td><td>Acyrthosiphon pisum</td><td>gi|59709604|gb|AY907547.1|</td></tr><tr><td>Rickettsia sp. (ref.1)</td><td>Acyrthosiphon pisum</td><td>gi|70568364|dbj|AB196668.1|</td></tr><tr><td>Rickettsia sp. (ref.2)</td><td>Acyrthosiphon pisum</td><td>gi|1147763|gb|U42084.1|RPU42084</td></tr><tr><td>Rickettsiella viridis (ref.1)</td><td>Acyrthosiphon pisum</td><td>gi|315064938|dbj|AB522703.1|</td></tr><tr><td>Rickettsiella viridis (ref.2)</td><td>Acyrthosiphon pisum</td><td>gi|315064940|dbj|AB522705.1|</td></tr><tr><td>Serratia symbiotica (ref.1)</td><td>Acyrthosiphon pisum</td><td>gi|9588081|dbj|AB033778.1|</td></tr><tr><td>Serratia symbiotica (ref.2)</td><td>Acyrthosiphon pisum</td><td>gi|29569343|gb|AY136139.1|</td></tr><tr><td>Spiroplasma sp.</td><td>Drosophila hydei</td><td>gi|224797645|gb|FJ657213.1|</td></tr><tr><td>Spiroplasma sp.</td><td>Harmonia axyridis</td><td>gi|4582256|emb|AJ132412.1|</td></tr><tr><td>Spiroplasma sp.</td><td>Acyrthosiphon pisum</td><td>gi|13359308|dbj|AB048263.1|</td></tr><tr><td>Spiroplasma citri R8A2HP</td><td>Citrus</td><td>gi|310974985|ref|NR_036849.1|</td></tr><tr><td>Spiroplasma syrphidicola</td><td>Eristalis arbustorum</td><td>gi|219846121|ref|NR_025711.1|</td></tr><tr><td>Wolbachia sp.</td><td>Cinara cedri</td><td>gi|52355655|gb|AY620430.1|</td></tr></tbody></table>

doi:10.1371/journal.pone.0120664.t002

## Bacterial community analyses

To analyze the bacterial community across the pea aphid samples, only presence-absence data were considered as both technological artifacts and innate biological traits bias relative quantification of each OTU abundance by 454 read counts [30]. The dissimilarity between the bacterial communities of pea aphid samples (beta diversity) was quantified by the Bray-Curtis distance, a metric based on the presence/absence of each OTU. Bray-Curtis dissimilarities between all pairwise combinations of pea aphid samples were summarized as a matrix (S2 Table). The Bray-Curtis dissimilarity matrix was ordinated following a non-metric multidimensional scaling (nMDS), a nonparametric ordination-based method where an iterative algorithm projects the multidimensional data of the dissimilarity matrix into a minimal dimensional space. The result of nMDS ordination is visualized on a scatter graph where the position of each pea aphid sample depends on its distance from all other points in the analysis. This method reduces ecological community data complexity and identifies meaningful relationships amongst communities. To analyze the effect of plant-based genetic differentiation in *A. pisum* on the bacterial community composition, we used a Mantel test to calculate the correlation between the matrix of Bray-Curtis dissimilarities between all pairwise combinations of pea aphid biotypes and the matrix of the genetic distance among these host-adapted aphid populations obtained from [31]. Finally, a Student test was performed to test whether the mean of the Bray-Curtis
dissimilarities between biotypes is higher than the mean of Bray-Curtis dissimilarities within biotypes. The 'vegan' package available in R was used to calculate Bray-Curtis dissimilarities and to perform both the nMDS and the Mantel test.

## Results

### Overall bacterial diversity and taxonomic assignments

Following our selection steps on sequence quality, size and redundancy between PCR replicates, 138,539 sequences were obtained, with a mean of 11,545 reads per sample. Overall, these sequence reads were classified into 21 Operational Taxonomic Units (OTUs). However, only 13 OTUs had a frequency exceeding 0.1% (i.e. approx. 140 reads) among which 3 OTUs represented more than 94% of the sequences. The classification from Greengenes database resulted in 21 unique bacterial genera that belonged to 4 phyla, 6 classes, 8 orders and 14 families (S3 Table). The Proteobacteria represented 51% of these classified bacteria, with the Alphaproteobacteria and the Gammaproteobacteria being the commonest taxa with 25% and 26%, respectively. Among the 138,539 sequences of 400–430 bp, 98% were assigned to bacterial taxa already reported as heritable symbionts of aphids. *Spiroplasma* was the most represented taxon in number of sequences (48%) followed by *Rickettsia* (25%) and *Buchnera* (21%). The frequencies of the other pea aphid symbionts were <2% (Fig 1). We confirmed the presence of *Wolbachia*, which was recently reported in North American *A. pisum* by [14], in two pea aphid samples (total frequency of 0.15%). Apart from these symbionts, 454 sequencing detected other bacteria not known to be heritable in arthropods, although most did not exceed 0.1% of sequence reads. Among those with frequencies above 0.1% were *Erwinia*, *Pantoea*, *Propionibacterium* and *Ralstonia*, each being found in multiple populations and for a total of more than 150 reads each. *Erwinia* and *Pantoea* have been described as aphid gut associates [32, 33]. *Ralstonia* encompass several soil borne plant pathogens that could be ingested or vectorized by aphids. *Propionibacterium* are usually found in the skin of humans and other animals and their occurrence in our 454 reads may represent contaminants from human handling. Bacterial taxa with read abundance below 0.1% were either plant associates such as *Rhizobium* and *Xanthomonas* or environmental/contaminant bacteria.

Taxonomic assignment based on the Greengenes database was very well supported by the Neighbor-Joining tree of 16S rDNA sequences from our 454 dataset and from references selected in GenBank (Fig 2). OTUs assigned to reported aphid symbionts clustered with the reference sequence from a pea aphid host. None of our 454 sequences grouped with the PAXS reference, probably because this recently described pea aphid symbiont was not detected due to poor sequence identity with our universal primers (sequence homology with PAXS was 100% for the reverse primer but only 60% for the forward primer). *Erwinia* and *Pantoea* appear as two distinct clades on the phylogenetic tree. Globally, there was a significant range of 16S rDNA variation within bacterial taxa, which could be due to either sequence errors in both 454 and Genbank sequences or variation between strains of bacterial taxa from different populations or biotypes of *A. pisum*.

### Bacterial community composition and host-adapted pea aphid populations

Bacterial diversity in each pea aphid sample (alpha diversity) was relatively low with 3–11 unique OTUs per sample (mean = 6.2). Several bacterial OTU were shared across pea aphid samples but some samples were characterized by specific bacterial communities (Table 3). Not only were the bacterial communities of pea aphid samples dominated by only a few OTUs (low alpha
<figure><img src="image_4.png" /><figcaption>Fig 1. Number of 16S rDNA reads in Log(10) for the 21 bacterial taxa detected in nine plant-specialized biotypes of the pea aphid, *Acyrthosiphon pisum*. *obligatory symbiont and **facultative symbionts.</figcaption></figure>

doi:10.1371/journal.pone.0120664.g001

diversity), but the distribution of the bacterial OTUs was often limited to only a few pea aphid samples (high beta diversity). Bray-Curtis dissimilarity indices calculated on the 21 validated OTUs distributed among our samples revealed that populations from the same pea aphid bio-type harbor similar bacterial communities (Fig 3). The statistical analysis confirmed this pattern since the bacterial community dissimilarity was significantly lower within than among pea aphid biotypes (respectively 0.31±0.07 and 0.55±0.02, t = 3.56, p = 0.028). Finally, the difference of bacterial community composition across the pea aphid biotypes was not related to the phylogenetic distance between these host-adapted aphid populations (Mantel test, r = -0.102, p = 0.634, Fig 4).

Discussion

Analyzing the microbiome of several pea aphid biotypes using deep 454 pyrosequencing of 16S rDNA genes reveals a relatively low bacterial diversity, with a total of 21 valid taxa detected
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Phylogenetic analysis of 16S rDNA based on a Neighbor-Joining tree relating bacterial sequences from references selected in GenBank (REF) and from 454 amplicon sequencing of populations and biotypes of the pea aphid complex. Accession numbers of representative sequences of each OTU found in the different samples (74 sequences in total) are given in Supporting Information (S3 Table).</figcaption></figure>

doi:10.1371/journal.pone.0120664.g002

and about 6 OTUs per sample on average. This study also shows that the pea aphid microbiota is dominated by very few symbionts, with 98% of the 138,539 reads assigned to the eight bacterial taxa already described as symbionts of this species.

Before discussing further these key results, caution should be taken with data generated by sequencing of 16S rDNA gene amplicons as this method suffers from several limitations [34]. First, the choice of universal primers may bias amplification towards taxa with higher sequence affinity, and sequence abundance per OTU may not reflect bacterial abundance [30]. Because of these possible biases, we decided to analyze qualitatively and not quantitatively our sequence data by taking into account only the presence or absence of OTUs among samples. Here, these artifacts could explain why *Spiroplasma*, and not the obligatory symbiont *Buchnera* as
<table><caption>Table 3. Distribution of the 21 bacterial OTU across samples from various populations and biotypes of the pea aphid.</caption><thead><tr><th rowspan="2">OTU type</th><th rowspan="2">OTU name</th><th colspan="12">Pea aphid samples</th><th rowspan="2">OTU occurrence</th></tr><tr><th>Csco</th><th>Lpra</th><th>Mlup</th><th>Msat-1</th><th>Msat-2</th><th>Malb</th><th>Ospi</th><th>Psat-1</th><th>Psat-2</th><th>Svar</th><th>Tpra-1</th><th>Tpra-2</th></tr></thead><tbody><tr><td rowspan="7">Aphid symbionts</td><td>Buchnera</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>12</td></tr><tr><td>Hamiltonella</td><td>X</td><td></td><td></td><td>X</td><td>X</td><td>X</td><td>X</td><td></td><td>X</td><td></td><td></td><td></td><td>6</td></tr><tr><td>Regiella</td><td></td><td>X</td><td></td><td>X</td><td>X</td><td>X</td><td></td><td></td><td></td><td>X</td><td>X</td><td></td><td>6</td></tr><tr><td>Rickettsia</td><td>X</td><td>X</td><td>X</td><td></td><td></td><td></td><td>X</td><td>X</td><td>X</td><td></td><td></td><td></td><td>6</td></tr><tr><td>Rickettsiella</td><td></td><td></td><td></td><td>X</td><td>X</td><td></td><td></td><td>X</td><td>X</td><td></td><td></td><td></td><td>4</td></tr><tr><td>Serratia</td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td>X</td><td>X</td><td></td><td></td><td>4</td></tr><tr><td>Spiroplasma</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td></td><td></td><td>X</td><td>X</td><td></td><td>X</td><td>X</td><td>9</td></tr><tr><td>Wolbachia</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td>X</td><td></td><td>2</td></tr><tr><td rowspan="2">Gut associates</td><td>Erwinia</td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td></td><td>X</td><td></td><td>3</td></tr><tr><td>Pantoea</td><td></td><td></td><td></td><td>X</td><td></td><td></td><td></td><td></td><td>X</td><td></td><td></td><td></td><td>2</td></tr><tr><td rowspan="5">Plant associates</td><td>Bradyrhizobium</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td></td><td>1</td></tr><tr><td>Mesorhizobium</td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>1</td></tr><tr><td>Pseudomonas</td><td></td><td></td><td></td><td></td><td>X</td><td></td><td></td><td></td><td>X</td><td></td><td></td><td></td><td>2</td></tr><tr><td>Ralstonia</td><td>X</td><td></td><td></td><td>X</td><td></td><td></td><td></td><td>X</td><td>X</td><td></td><td>X</td><td>X</td><td>6</td></tr><tr><td>Xanthomonas</td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td>1</td></tr><tr><td rowspan="6">Environmental bacteria</td><td>Cupriavidus</td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>1</td></tr><tr><td>Corynebacterium</td><td></td><td></td><td></td><td></td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td>1</td></tr><tr><td>Lactobacillus</td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>1</td></tr><tr><td>Microbacterium</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td></td><td></td><td>1</td></tr><tr><td>Paenibacillus</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td></td><td></td><td></td><td></td><td>1</td></tr><tr><td>Propionibacterium</td><td>X</td><td></td><td></td><td></td><td></td><td></td><td>X</td><td>X</td><td></td><td>X</td><td></td><td></td><td>4</td></tr><tr><td colspan="2">OTU number</td><td>11</td><td>4</td><td>3</td><td>5</td><td>7</td><td>4</td><td>5</td><td>8</td><td>11</td><td>3</td><td>9</td><td>4</td><td></td></tr></tbody></table>

Names of samples are those indicated in Table 1.

doi:10.1371/journal.pone.0120664.t003

one would expect, showed the highest abundance of 454 reads in our dataset. Second, the sequencing depth may limit the detection of rare taxa, leading to an underestimation of microbial diversity. However, 11,500 reads per sample were obtained on average here, a sequencing effort which is likely to capture most if not all of bacterial diversity and indeed rarefaction curves reached saturation in all samples (data not shown). Third, the size and quality of sequences may also prevent the accurate taxonomic assignment of OTUs formed by the clustering analysis of 16S rDNA sequences. Here, we applied stringent criteria to filter the 454 reads, by removing two third of them and keeping only high quality and long (400–430bp) sequences. Phylogenetic analysis of 16S rDNA sequences generated by 454 amplicon pyrosequencing validated taxonomic assignment since sequences assigned to bacterial taxa associated with aphids clustered with the reference sequences.

With these cautions in mind, our deep sequencing of 16S rDNA amplicons from multiple populations and biotypes of the pea aphid complex reveals the limited diversity of aphid bacterial communities, which are dominated by symbionts among which facultative symbionts vary drastically between host-specialized biotypes and less between localities. This strong link between plant specialization and bacterial communities confirms previous studies in European and American populations of *A. pisum* based on specific PCR-based detection of a range of facultative symbionts [20], [23], [35], [36], [37]. Such a relationship between plant specialization and microbial communities has been described in other insect species, including aphids [6], [38], [39], [40], [41]. However, it is uncertain whether symbiont communities associated with
<figure><img src="image_6.png" /><figcaption>**Fig 3.** Non-metric MDS ordination plot comparing bacterial communities from different pea aphid samples. Each data point represents the bacterial community identified from a single sample (see Table 1 for dot label legends). The Bray-Curtis dissimilarity index was used to rank distances calculated using the presence-absence community data. Stress of the nMDS = 0.152.</figcaption></figure>

doi:10.1371/journal.pone.0120664.g003

herbivores exert a direct influence on plant specialization. In the pea aphid complex, host ecology more than geography seems to shape symbiont composition [20] but we lack firm conclusions about the direct involvement of facultative symbionts in plant adaptation of their aphid hosts. It has been shown that *Regiella*-infected pea aphids present increased performances on white clover compared to those deprived of this facultative symbiont [16]. However, other studies have contradicted this finding [42], [43] and generalization to other symbionts does not support symbiont-driven plant specialization in *A. pisum* [4], [44]. Since the genetic divergence between pea aphid biotypes is not related with their difference in microbial communities (Fig 4), historical factors seem to have no or very little role in shaping symbiont composition in the complex. Therefore, ecological forces other than the direct plant influence (e.g. enemy pressure), as well as horizontal transfers are likely to be key factors in driving symbiont community structure among host-adapted biotypes [20], [44], [45].

Here, we detected *Wolbachia* in two biotypes (*Securigera* and *Trifolium*-adapted biotypes), confirming the recent discovery of this endosymbiont in *A. pisum* in North America and China [14], [46]. However, the prevalence of this symbiont, and therefore its impact on host fitness, must be very limited since it has been searched intensively before its recent detection
<figure><img src="image_7.png" /><figcaption>**Fig 4.** Relationship between the matrix of Bray-Curtis dissimilarities between all pairwise combinations of pea aphid biotypes and the matrix of the genetic distance among these host-adapted aphid populations. Each data point represents a pairwise combination of pea aphid biotype (Mantel test, $p>0.05$).</figcaption></figure>

doi:10.1371/journal.pone.0120664.g004

[47], [48] and was found in only one *A. pisum* individual from North American populations among 318 screened [14]. Larger screenings on more biotypes and individuals are needed to better assess its distribution and prevalence across the pea aphid complex. More work is also required to examine the localization of *Wolbachia* in the insect, its potential effects and their nature on the host's phenotype and fitness [14], [46], [47].

Besides heritable symbionts known for the pea aphid, *Erwinia* and *Pantoea* were detected in four samples with a substantial number of reads (>0.1%, 1068 reads in total). Interestingly these two bacteria are ubiquitous plant pathogens but have been also described as gut associates in *A. pisum* as well as in other aphid species [49], [33], [39], [50]. It is not clear whether these two gut bacteria are aphid pathogens, commensals, mutualists or use aphids as alternative hosts to plants [51], [52], [53]. *Erwinia aphidicola*, which has been isolated from the pea aphid gut and described by [32, 33] shows virtually complete 16S rDNA sequence homology with the three representative sequences of the OTU assigned to *Erwinia* in our dataset. *E. aphidicola* has been reported as pathogen for the pea aphid although negative effects have been detected in laboratory conditions and with concentration that may not be realistic [53]. It may be that
*E. aphidicola* provides some benefit to the aphid under certain ecological conditions as shown for *Erwinia* bacteria which are gut associates of thrips and have diet-dependent effects on their insect hosts [54]. *Pantoea agglomerans*-like bacteria have been identified in phylloxera-inducing species where its prevalence can reach 100% as in some *Daktulosphaira vitifoliae* populations [39], [50]. Sequences of the OTU assigned to *Pantoea* clustered into the same clade as the reference sequences for *P. agglomerans*, suggesting that some pea aphid individuals are infected by this taxon which is related to phytopathogenic bacteria, but can apparently colonize insect gut with totally unknown effects. It is unlikely that these two gut bacteria are transmitted vertically because their chances of transmission from mother to offspring are small [49], [55]. It is more plausible that aphids acquire these bacteria either from the honeydew where they are probably excreted or from the plant surface [56], [57]. Another possibility is that they are ingested from plant sap, providing that these bacteria may circulate into the phloem [58]. Further investigations are required to ascertain the importance of aphid gut associates in both laboratory and field conditions.

The remaining bacterial taxa detected in samples of pea aphid biotypes and populations had low read abundance and were either plant associates such as *Rhizobium* and *Xanthomonas* or environmental/contaminant bacteria. Although their effect on aphids either as hosts or vectors of plant pathogens may not be negligible [52], their presence requires confirmation through a larger survey.

In this study, host-specialized biotypes that were collected from the field were then grown on broad bean, the universal host-plant for all biotypes, during two generations. This step, while reducing the background noise of microbes associated with internal parasites, might have provoked changes in the microbiota associated with the host, in particular in the diversity of non-heritable bacteria. Comparison of microbiota of aphids directly taken from the field with present results is needed to assess quantitative and qualitative differences in microbial communities between field and laboratory conditions. In addition, shift in host-plant (from native plant to universal host) imposed by our experimental design might have altered the microbial community associated with pea aphid biotypes and populations. However, this hypothesis is unlikely, at least for heritable symbionts whose vertical transmission is almost perfect and loss, very uncommon [3].

In conclusion, despite the artifacts inherent to sequencing of 16S rDNA gene amplicons and possible changes in microbial diversity resulting from aphid culture under laboratory conditions, our study provides qualitative results consistent with previous works based on taxon-specific PCR-based detection and more recent ones using without a priori approaches. In line with these studies, we found that bacterial communities associated with pea aphid hosts are largely dominated by a few symbiont taxa whose distribution varies drastically among biotypes. Among herbivore insects examined so far, sap-feeders such as aphids, psyllids and whiteflies show the poorest microbial diversity, with less than 12 OTUs per sample [14], [34], [59], [present study], which has been attributed to several not mutually exclusive factors involving nutrition, immunity and antagonisms [34]. While this diversity may be enriched by transient bacterial taxa (e.g. on the cuticle or in the gut) or by other microbes such as viruses and fungi, the ecological and evolutionary importance of these additional microbial associates for their hosts would probably be limited compared to obligate and facultative symbionts.

Our work leaves unresolved several important issues regarding the impact of microbial communities on their aphid hosts at individual, population and species levels. These include notably the effective importance and role of gut associates as well as the nature and relative influence of forces that shape microbial community structure, with particular attention to ecological factors and transmission patterns. Accessing the range and distribution of microbial diversity with NGS technologies applied to host populations from a large array of environmental conditions
will certainly help in resolving these key issues, completed by more focused biological and
functional analyses.

**Supporting Information**

<caption>S1 Table. Genbank accession numbers of 16S rDNA sequences of each OTU found in the different samples of pea aphids and assigned to bacterial taxa.</caption>

(XLSX)

<caption>S2 Table. Matrix of Bray-Curtis dissimilarities between all pairwise combinations of pea aphid samples.</caption>

(XLSX)

<caption>S3 Table. Taxonomic affiliation of the 21 bacterial OTUs detected in nine plant-specialized biotypes of the pea aphid, *Acyrtosiphon pisum*. Numbers and percentages of 16S rDNA reads are given for each taxon.</caption>

(XLSX)

**Acknowledgments**

We thank Philippe Vandenkoornhuyse, Sophie Coudouel and Alexandra Dheilly from the Environmental genomics platform of Biogenouest. Denis Poinsot is acknowledged for his contribution at various stages of this work and Christophe Mougel for his advice on sequence data analysis.

**Author Contributions**

Conceived and designed the experiments: JCS YO. Performed the experiments: JPG LM. Analyzed the data: JPG YO JCS. Contributed reagents/materials/analysis tools: JCS YO. Wrote the paper: JCS YO JPG.

**References**

1. Douglas AE. Microbial Brokers of Insect-Plant Interactions Revisited. J Chem Ecol. 2013; 39(7):952-61. doi: 10.1007/s10886-013-0308-x PMID: 23793897

2. Moran NA, McCutcheon JP, Nakabachi A. Genomics and Evolution of Heritable Bacterial Symbionts. Annual Review of Genetics. 2008; 42(1):165-90. doi: 10.1146/annurev.genet.41.110306.130119

3. Oliver KM, Degnan PH, Burke GR, Moran NA. Facultative Symbionts in Aphids and the Horizontal Transfer of Ecologically Important Traits. Annual Review of Entomology. Annual Review of Entomology. 2010. p. 247-66.

4. Hansen AK, Moran NA. The impact of microbial symbionts on host plant utilization by herbivorous insects. Molecular Ecology. 2014; 23(6):1473-96. doi: 10.1111/mec.12421 PMID: 23952067

5. Boissiere A, Tchioffo MT, Bachar D, Abate L, Marie A, Nsango SE, et al. Midgut Microbiota of the Malaria Mosquito Vector *Anopheles gambiae* and Interactions with *Plasmodium falciparum* Infection. Plos Pathogens. 2012; 8(5). e1002742 doi: 10.1371/journal.ppat.1002742 PMID: 22693451

6. Chandler JA, Morgan Lang J, Bhatnagar S, Eisen JA, Kopp A. Bacterial Communities of Diverse Drosoptila Species: Ecological Context of a Host-Microbe Model System. PLoS Genet. 2011; 7(9): e1002272. doi: 10.1371/journal.pgen.1002272 PMID: 21966276

7. Engel P, Martinson VG, Moran NA. Functional diversity within the simple gut microbiota of the honey bee. Proceedings of the National Academy of Sciences of the United States of America. 2012; 109(27):11002-7. doi: 10.1073/pnas.1202970109 PMID: 22711827

8. Osei-Poku J, Mbogo CM, Palmer WJ, Jiggins FM. Deep sequencing reveals extensive variation in the gut microbiota of wild mosquitoes from Kenya. Molecular Ecology. 2012; 21(20):5138-50. doi: 10.1111/j.1365-294X.2012.05759.x PMID: 22988916
9. Moran NA, Russell JA, Koga R, Fukatsu T. Evolutionary relationships of three new species of Enterobacteriaceae living as symbionts of aphids and other insects. Appl Environ Microbiol. 2005; 71(6):3302-10. doi: 10.1128/aem.71.6.3302-3310.2005 PMID: 15933033

10. Tsuchida T, Koga R, Fujiwara A, Fukatsu T. Phenotypic Effect of "Candidatus Rickettsiella viridis," a Facultative Symbiont of the Pea Aphid (Acyrthosiphon pisum), and Its Interaction with a Coexisting Symbiont. Appl Environ Microbiol. 2014; 80(2):525-33. doi: 10.1128/aem.03049-13 PMID: 24212575

11. Tsuchida T, Koga R, Horikawa M, Tsunoda T, Maoka T, Matsumoto S, et al. Symbiotic Bacterium Modifies Aphid Body Color. Science. 2010; 330(6007):1102-4. doi: 10.1126/science.1195463 PMID: 21097935

12. Guay JF, Boudreault S, Michaud D, Cloutier C. Impact of environmental stress on aphid clonal resistance to parasitoids: Role of *Hamiltonella defensa* bacterial symbiosis in association with a new facultative symbiont of the pea aphid. Journal of Insect Physiology. 2009; 55(10):919-26. doi: 10.1016/j.jinsphys.2009.06.006 PMID: 19545573

13. Chen DQ, Campbell BC, Purcell AH. A new Rickettsia from a herbivorous insect, the pea aphid *Acyrthosiphon pisum* (Harris). Current Microbiology. 1996; 33(2):123-8. doi: 10.1007/s002849900086 PMID: 8662184

14. Russell JA, Weldon S, Smith AH, Kim KL, Hu Y, Łukasik P, et al. Uncovering symbiont-driven genetic diversity across North American pea aphids. Molecular Ecology. 2013; 22(7):2045-59. doi: 10.1111/mec.12211 PMID: 23379399

15. Fukatsu T, Tsuchida T, Nikoh N, Koga R. Spiroplasma symbiont of the pea aphid, *Acyrthosiphon pisum* (Insecta: Homoptera). Appl Environ Microbiol. 2001; 67(3):1284-91. doi: 10.1128/aem.67.3.1284-1291.2001 PMID: 11229923

16. Tsuchida T, Koga R, Fukatsu T. Host plant specialization governed by facultative symbiont. Science. 2004; 303(5666):1989-. doi: 10.1126/science.1094611 PMID: 15044797

17. Montllor CB, Maxmen A, Purcell AH. Facultative bacterial endosymbionts benefit pea aphids *Acyrthosiphon pisum* under heat stress. Ecological Entomology. 2002; 27(2):189-95. doi: 10.1046/j.1365-2311.2002.00393.x

18. Oliver KM, Russell JA, Moran NA, Hunter MS. Facultative bacterial symbionts in aphids confer resistance to parasitic wasps. Proceedings of the National Academy of Sciences of the United States of America. 2003; 100(4):1803-7. doi: 10.1073/pnas.0335320100 PMID: 12563031

19. Scarborough CL, Ferrari J, Godfray HCJ. Aphid protected from pathogen by endosymbiont. Science. 2005; 310(5755):1781-. doi: 10.1126/science.1120180 PMID: 16357252

20. Henry LM, Peccoud J, Simon JC, Hadfield JD, Maiden MJC, Ferrari J, et al. Horizontally Transmitted Symbionts and Host Colonization of Ecological Niches. Current Biology. 2013; 23(17):1713-7. doi: 10.1016/j.cub.2013.07.029 PMID: 23993843

21. Peccoud J, Ollivier A, Plantegenest M, Simon JC. A continuum of genetic divergence from sympatric host races to species in the pea aphid complex. Proceedings of the National Academy of Sciences of the United States of America. 2009; 106(18):7495-500. doi: 10.1073/pnas.0811117106 PMID: 19380742

22. Frantz A, Calcagno V, Mieuzet L, Plantegenest M, Simon JC. Complex trait differentiation between host-populations of the pea aphid *Acyrthosiphon pisum* (Harris): implications for the evolution of ecological specialisation. Biological Journal of the Linnean Society. 2009; 97(4):718-27.

23. Ferrari J, West JA, Via S, Godfray HCJ. Population genetic structure and secondary symbionts in host-associated populations of the pea aphid complex. Evolution. 2012; 66(2):375-90. doi: 10.1111/j.1558-5646.2011.01436.x PMID: 22276535

24. Jaquiyer J, Stoeckel S, Nouhaud P, Mieuzet L, Maheo F, Legeai F, et al. Genome scans reveal candidate regions involved in the adaptation to host plant in the pea aphid complex. Molecular Ecology. 2012; 21(21):5251-64. doi: 10.1111/mec.12048 PMID: 23017212

25. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: Open-Source, Platform-Independent, Community-Supported Software for Describing and Comparing Microbial Communities. Appl Environ Microbiol. 2009; 75(23):7537-41. doi: 10.1128/aem.01541-09 PMID: 19801464

26. Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimeria detection. Bioinformatics. 2011; 27(16):2194-200. doi: 10.1093/bioinformatics/btr381 PMID: 21700674

27. DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimeracHECKED 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006; 72(7):5069-72. doi: 10.1128/aem.03006-05 PMID: 16820507
28. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and clustal X version 2.0. Bioinformatics. 2007; 23(21):2947-8. doi: 10.1093/bioinformatics/btm404 PMID: 17846036

29. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Molecular Biology and Evolution. 2011; 28(10):2731-9. doi: 10.1093/molbev/msr121 PMID: 21546353

30. Amend AS, Seifert KA, Bruns TD. Quantifying microbial communities with 454 pyrosequencing: does read abundance count? Molecular Ecology. 2010; 19(24):5555-65. doi: 10.1111/j.1365-294X.2010.04898.x PMID: 21050295

31. Nouhaud P, Peccoud J, Mahéo F, Mieuzet L, Jaquiéry J, Simon JC. Genomic regions repeatedly involved in divergence among plant-specialized pea aphid biotypes. Journal of Evolutionary Biology. 2014; 27:2013-20. doi: 10.1111/jeb.12441 PMID: 24953130

32. Harada H, Oyaizu H, Ishikawa H. A consideration about the origin of aphid intracellular symbiont in connection with gut bacterial flora. Journal of General and Applied Microbiology. 1996; 42(1):17-26. doi: 10.2323/jgam.42.17

33. Harada H, Oyaizu H, Kosako Y, Ishikawa H. *Erwinia aphidicola*, a new species isolated from pea aphid, *Acyrthosiphon pisum*. Journal of General and Applied Microbiology. 1997; 43(6):349-54. doi: 10.2323/jgam.43.349 PMID: 12501306

34. Jing X, Wong ACN, Chaston JM, Colvin J, McKenzie CL, Douglas AE. The bacterial communities in plant phloem-sap-feeding insects. Molecular Ecology. 2014; 23(6):1433-44. doi: 10.1111/mec.12637 PMID: 24350573

35. Bilodeau E, Simon JC, Guay JF, Turgeon J, Cloutier C. Does variation in host plant association and symbiont infection of pea aphid populations induce genetic and behaviour differentiation of its main parasitoid, *Aphidius ervi*? Evolutionary Ecology. 2013; 27(1):165-84. doi: 10.1007/s10682-012-9577-z

36. Peccoud J, Figueroa CC, Silva AX, Ramirez CC, Mieuzet L, Bonhomme J, et al. Host range expansion of an introduced insect pest through multiple colonizations of specialized clones. Molecular Ecology. 2008; 17(21):4608-18. doi: 10.1111/j.1365-294X.2008.03949.x PMID: 19140984

37. Simon JC, Carre S, Boutin M, Prunier-Leterme N, Sabater-Munoz B, Latorre A, et al. Host-based divergence in populations of the pea aphid: insights from nuclear markers and the prevalence of facultative symbionts. Proceedings of the Royal Society B-Biological Sciences. 2003; 270(1525):1703-12. doi: 10.1098/rspb.2003.2430 PMID: 12964998

38. Brady CM, White JA. Cowpea aphid (*Aphis craccivora*) associated with different host plants has different facultative endosymbionts. Ecological Entomology. 2013; 38(4):433-7. doi: 10.1111/een.12020

39. Medina RF, Nachappa P, Tamborindeguy C. Differences in bacterial diversity of host-associated populations of *Phylloxera notabilis* Pergande (Hemiptera: Phylloxeridae) in pecan and water hickory. Journal of Evolutionary Biology. 2011; 24(4):761-71. doi: 10.1111/j.1420-9101.2010.02215.x PMID: 21261774

40. Shi WB, Xie SX, Chen XY, Sun S, Zhou X, Liu LT, et al. Comparative Genomic Analysis of the Endosymbionts of Herbivorous Insects Reveals Eco-Environmental Adaptations: Biotechnology Applications. Plos Genetics. 2013; 9(1). e1003131 doi: 10.1371/journal.pgen.1003131 PMID: 23326236

41. Toju H, Fukatsu T. Diversity and infection prevalence of endosymbionts in natural populations of the chestnut weevil: relevance of local climate and host plants. Molecular Ecology. 2011; 20(4):853-68. doi: 10.1111/j.1365-294X.2010.04980.x PMID: 21199036

42. Leonardo TE. Removal of a specialization-associated symbiont does not affect aphid fitness. Ecol Lett. 2004; 7(6):461-8. doi: 10.1111/j.1461-0248.2004.00602.x

43. Ferrari J, Scarborough CL, Godfray HCJ. Genetic variation in the effect of a facultative symbiont on host-plant use by pea aphids. Oecologia. 2007; 153(2):323-9. doi: 10.1007/s00442-007-0730-2 PMID: 17415589

44. McLean AHC, van Asch M, Ferrari J, Godfray HCJ. Effects of bacterial secondary symbionts on host plant use in pea aphids. Proceedings of the Royal Society B-Biological Sciences. 2011; 278(1706):760-6. doi: 10.1098/rspb.2010.1654 PMID: 20843842

45. Frago E, Dicke M, Godfray HCJ. Insect symbionts as hidden players in insect-plant interactions. Trends in Ecology & Evolution. 2012; 27(12):705-11. doi: 10.1016/j.tree.2012.08.013

46. Wang Z, Su XM, Wen J, Jiang LY, Qiao GX. Widespread infection and diverse infection patterns of Wolbachia in Chinese aphids. Insect Sci. 2014; 21(3):313-25. doi: 10.1111/1744-7917.12102 PMID: 24395812

47. Augustinos AA, Santos-Garcia D, Dionyssopoulou E, Moreira M, Papapanagiotou A, Scarvelakis M, et al. Detection and Characterization of Wolbachia Infections in Natural Populations of Aphids: Is the Hidden Diversity Fully Unraveled? PLoS ONE. 2011; 6(12). e28695 doi: 10.1371/journal.pone.0028695 PMID: 22174869
48. Jeyaprakash A, Hoy MA. Long PCR improves Wolbachia DNA amplification: wsp sequences found in 76% of sixty-three arthropod species. Insect Molecular Biology. 2000; 9(4):393-405. doi: 10.1046/j.1365-2583.2000.00203.x PMID: 10971717

49. Clark EL, Daniell TJ, Wishart J, Hubbard SF, Karley AJ. How Conserved Are the Bacterial Communities Associated With Aphids? A Detailed Assessment of the *Brevicoryne brassicae* (Hemiptera: Aphididae) Using 16S rDNA. Environmental Entomology. 2012; 41(6):1386-97. doi: 10.1603/en12152 PMID: 23321084

50. Vorwerk S, Martinez-Torres D, Forneck A. *Pantoea agglomerans*-associated bacteria in grape phylloxera (*Daktulosphaira vitifoliae*, Fitch). Agricultural and Forest Entomology. 2007; 9(1):57-64. doi: 10.1111/j.1461-9563.2006.000319.x

51. Costechareyre D, Balmand S, Condemine G, Rahbe Y. *Dickeya dadantii*, a Plant Pathogenic Bacterium Producing Cyt-Like Entomotoxins, Causes Septicemia in the Pea Aphid *Acyrthosiphon pisum*. PLoS ONE. 2012; 7(1). e30702 doi: 10.1371/journal.pone.0030702 PMID: 22292023

52. Stavrinides J, No A, Ochman H. A single genetic locus in the phytopathogen *Pantoea stewartii* enables gut colonization and pathogenicity in an insect host. Environmental Microbiology. 2010; 12(1):147-55. doi: 10.1111/j.1462-2920.2009.02056.x PMID: 19788413

53. Harada H, Ishikawa H. Experimental pathogenicity of *Erwinia aphidicola* to pea aphid, *Acyrthosiphon pisum*. Journal of General and Applied Microbiology. 1997; 43(6):363-7. doi: 10.2323/jgam.43.363 PMID: 12501308

54. de Vries EJ, Jacobs G, Sabelis MW, Menken SBJ, Breeuwer JAJ. Diet-dependent effects of gut bacteria on their insect host: the symbiosis of *Erwinia* sp and western flower thrips. Proceedings of the Royal Society B-Biological Sciences. 2004; 271(1553):2171-8. doi: 10.1098/rspb.2004.2817 PMID: 15475338

55. Engel P, Moran NA. The gut microbiota of insects—diversity in structure and function. Fems Microbiol-gy Reviews. 2013; 37(5):699-735. doi: 10.1111/1574-6976.12025 PMID: 23692388

56. Stavrinides J, McCloskey JK, Ochman H. Pea Aphid as both Host and Vector for the Phytopathogenic Bacterium *Pseudomonas syringae*. Appl Environ Microbiol. 2009; 75(7):2230-5. doi: 10.1128/aem.02860-08 PMID: 19201955

57. Sabri A, Vandermoten S, Leroy PD, Haubruge E, Hance T, Thonart P, et al. Proteomic Investigation of Aphid Honeydew Reveals an Unexpected Diversity of Proteins. PLoS ONE. 2013; 8(9). e74656 doi: 10.1371/journal.pone.0074656 PMID: 24086359

58. Caspi-Fluger A, Inbar M, Mozes-Daube N, Katzir N, Portnoy V, Belausov E, et al. Horizontal transmission of the insect symbiont Rickettsia is plant-mediated. Proceedings of the Royal Society B-Biological Sciences. 2012; 279(1734):1791-6. doi: 10.1098/rspb.2011.2095 PMID: 22113034

59. Colman DR, Toolson EC, Takacs-Vesbach CD. Do diet and taxonomy influence insect gut bacterial communities? Molecular Ecology. 2012; 21(20):5124-37. doi: 10.1111/j.1365-294X.2012.05752.x PMID: 22978555