Detection and Characterization of Wolbachia Infections in Natural Populations of Aphids: Is the Hidden Diversity Fully Unraveled?

Abstract

Aphids are a serious threat to agriculture, despite being a rather small group of insects. The about 4,000 species worldwide engage in highly interesting and complex relationships with their microbial fauna. One of the key symbionts in arthropods is Wolbachia, an α-Proteobacterium implicated in many important biological processes and believed to be a potential tool for biological control. Aphids were thought not to harbour Wolbachia; however, current data suggest that its presence in aphids has been missed, probably due to the low titre of the infection and/or to the high divergence of the Wolbachia strains of aphids. The goal of the present study is to map the Wolbachia infection status of natural aphids populations, along with the characterization of the detected Wolbachia strains. Out of 425 samples from Spain, Portugal, Greece, Israel and Iran, 37 were found to be infected. Our results, based mainly on 16S rRNA gene sequencing, indicate the presence of two new Wolbachia supergroups prevailing in aphids, along with some strains belonging either to supergroup B or to supergroup A.

About the Authors:

Antonis A. Augustinos

Affiliations Department of Environmental and Natural Resources Management, University of Ioannina, Agrinio, Greece, Department of Environmental and Natural Resources Management, University of Western Greece, Agrinio, Greece

Diego Santos-Garcia

Affiliation: Institut Cavanilles de Biodiversitat i Biologia Evolutiva, Universitat de València, Valencia, Spain

Eva Dionyssopoulou

Affiliations Department of Environmental and Natural Resources Management, University of Ioannina, Agrinio, Greece, Department of Environmental and Natural Resources Management, University of Western Greece, Agrinio, Greece

Marta Moreira

Affiliation: ISOPlexis Gene Bank, Universidade da Madeira, Funchal, Portugal

Aristeidis Papapanagiotou

Affiliation: Department of Greenhouse Crops and Floriculture, Technological Educational Institute of Messolonghi, Messolonghi, Greece

Marios Scarvelakis

Affiliation: Department of Environmental and Natural Resources Management, University of Ioannina, Agrinio, Greece

Vangelis Doudoumis

Affiliations Department of Environmental and Natural Resources Management, University of Ioannina, Agrinio, Greece, Department of Environmental and Natural Resources Management, University of Western Greece, Agrinio, Greece

Silvia Ramos

Affiliation: Institut Cavanilles de Biodiversitat i Biologia Evolutiva, Universitat de València, Valencia, Spain

Antonio F. Aguiar

Affiliation: Laboratório de Qualidade Agrícola, Núcleo de Fitopatologia,, Madeira, Portugal

Paulo A. V. Borges

Affiliation: Departamento de Ciências Agrárias CITA-A (Azorean Biodiversity Group), Universidade dos Açores, Angra do Heroísmo, Terceira – Azores

Manhaz Khadem

Affiliation: ISOPlexis Gene Bank, Universidade da Madeira, Funchal, Portugal

Amparo Latorre

Affiliations Institut Cavanilles de Biodiversitat i Biologia Evolutiva, Universitat de València, Valencia, Spain, Área de Genómica y Salud, Centro Superior de Investigación en Salud Pública (CSISP), Valencia, Spain

George Tsiamis

Affiliations Department of Environmental and Natural Resources Management, University of Ioannina, Agrinio, Greece, Department of Environmental and Natural Resources Management, University of Western Greece, Agrinio, Greece

Kostas Bourtzis

* E-mail: [email protected]

Affiliations Department of Environmental and Natural Resources Management, University of Ioannina, Agrinio, Greece, Biomedical Sciences Research Center Al. Fleming, Vari, Greece, Department of Environmental and Natural Resources Management, University of Western Greece, Agrinio, Greece

Introduction

Wolbachia is a diverse group of obligatory intracellular and maternally transmitted α-Proteobacteria [1]–[3]. Several studies suggest that these bacteria are present in at least 65% of arthropod species as well as in filarial nematodes and in some plant parasitic nematodes [4]–[8]. Wolbachia strains infecting arthropod and nematode hosts are represented by a single species, Wolbachia pipientis [9]; however, there is extensive diversity which has resulted in the assignment of the bacterial strains into at least eleven Wolbachia supergroups, named A to F and H to L (supergroup G is considered a recombinant between A and B) [4], [10]–[19].Wolbachia diversity was initially characterized using the genes wsp, 16S rRNA, ftsZ, gltA and groEL as molecular markers, while strain genotyping is based on multi locus sequence typing systems (MLST), as well as on the amino acid sequences of the four hypervariable regions (HVRs) of the WSP protein [20], [21].

Wolbachia have been reported in the somatic tissues of arthropod hosts; however, they mainly reside in the reproductive tissues and organs [2]. This tissue localization pattern has been associated with the induction of different reproductive alterations such as feminization, parthenogenesis, male killing and cytoplasmic incompatibility [2], [22], which aid the spread of Wolbachia infections in host populations [23]. The widespread distribution of Wolbachia and their ability to manipulate the reproductive properties of arthropod hosts has attracted interest in its role in host biology, ecology and evolution, as well as in the development of novel, symbiont-based and environment friendly Wolbachia-based methods for pest and disease management [2], [3], [24]–[26]. It has been suggested that Wolbachia-induced cytoplasmic incompatibility can be used either for the control of agricultural pests and disease vectors through the Incompatible Insect Technique (IIT), or by spreading a desirable genotype through populations, such as the inability of a vector species to transmit a pathogen [27]–[33]. The introduction of life-shortening Wolbachia strains could modify the population age structure of insect vector species, thus reducing pathogen transmission [34], [35]. Furthermore, recent studies provide evidence that the presence of Wolbachia in some insect species may provide anti-viral protection as well as inhibit the infection with and transmission of certain pathogens such as Dengue, Chikungunya and Plasmodium [35]–[40].

Aphids are a rather small group of insects but their threat to agricultural ecosystems is enormous. Currently, there are about 4,000 recognized species worldwide [41]. Aphids do great damage to their host plants in several ways [42]. They feed on plant sap and inject saliva (which can be phytotoxic) during feeding. Their honeydew is used by saprophytic ascomycetes that grow on plants. More importantly, aphids have been shown to be vectors of numerous plant viruses. Due to their feeding behavior, they are by far the most important virus vectors, transmitting ∼30% of all plant virus species [43].

Aphids exhibit many interesting biological traits. They have a complicated life cycle, being able to reproduce both sexually and asexually. They are specialized in probing and using phloem sap as sole food source, which leads to a tight association with their host plants. They are also important for the feeding of other insects; they modify phloem sap, which has a high ratio of non-essential to essential amino acids and elevated sugar content, and produce substances more suitable for other species [44].

Aphids have established sophisticated symbiotic relationships and many of their unique properties can be attributed to their symbiotic bacteria [45]. They have established an obligate mutualistic symbiosis with Buchnera aphidicola, whichprovides them with essential amino acids lacking from their phloem diet [46]–[49]. Occasionally, aphids harbour secondary or facultative symbionts that coexist with Buchnera, and can have positive effects on the aphid host [45]. It has been reported that ‘CandidatusHamiltonella defensa’ and ‘CandidatusRegiella insecticola’can protect aphids against parasitoids [50], [51], whereas CandidatusSerratia symbioticais implicated in heat tolerance [52]. Finally, studies showing lateral gene transfer from secondary symbionts to their aphid host and the fact that these genes are expressed in some cases [53], [54], along with a reported case of metabolic complementation between B. aphidicola and “Ca S. symbiotica” in the aphid Cinara cedri [55], [56] illustrate the very complex relationship between aphids and their symbionts. All the above suggest that aphids, together with their host plants and their microbial fauna, not only constitute an interesting biological model worth investigating, but that it is furthermore crucial to study and understand these relationships in order to devise appropriate control methods for these species and the plant diseases they transmit.

A small number of studies has investigated the presence of Wolbachia in aphids [6], [57]–[60]. Most of them failed to detect Wolbachia[57]–[59]. The first report of aphids (Toxoptera citricidusand Aphis craccivora) harboring Wolbachia was based on Long-PCR and the sequencing of the wsp gene [6]. Stronger evidence for the presence of Wolbachia in aphid species was based on 16S rDNA sequencing, electron microscopy and in situ localization of this endosymbiont in C. cedri [60]. It was recently reported that Chinese natural populations of the wheat aphid, Sitobion miscanthi, harbour single and/or double Wolbachia infections belonging to the A and B supergroup [61].

We undertook extensive screening and report here on the presence of Wolbachia infections in natural populations of aphid species. The characterization of these Wolbachia strains is based on the use of gene markers 16S rRNA, ftsZ, gltA, groEL, wsp and MLST. Our study suggests that neither the detection nor the unraveling of Wolbachia diversity in the aphid fauna is an easy task; they demand the development of novel tools.

Results

Screening for Wolbachia infections in natural populations of aphids

A total of 425 natural samples of aphids were screened for the presence of Wolbachia with a 16S rRNA-based PCR approach using the wspecF/wspecR set of primers (Figure S1). The samples were collected in five countries (Greece, Spain, Portugal, Israel and Iran) and on a variety of host plants (at least 165 different species). Collections were in some cases diachronic. This collection represents 144 different aphid species within 69 genera ofnine subfamilies of the family Aphididae (Aphidinae, Chaitophorinae, Pterocommatinae, Myzocallidinae, Drepanosiphinae, Thelaxinae, Lachninae, Mindarinae and Eriosomatinae) (Table S1). The majority of samples screened belong to the subfamily Aphidinae (tribes Aphidini and Macrosiphini)followed by the subfamily Lachninae (mainly from the Eulachnini tribe), considered by recent studies as the most basal lineage among the aphid subfamilies [62],. and from all three tribes of Eriosomatinae subfamily (Pemphigini, Eriosomatini and Fordini).

The results of the screen, which are presented in Table 1, show that the prevalence of Wolbachia infections varied significantly between different aphid populations and can be summarized as follows: (a) Wolbachia infection was detected in only 37 out of 425 aphid populations tested; (b) Wolbachia was detected in aphid species of the subfamilies Lachninae, Aphidinae, Chaitophorinae Eriosomatinae, and Drepanosiphinae, while no infection was found in the rest of subfamilies; (c) at least eight species of the Lachninae subfamily were found infected including seven Cinara species (C. fresai, C. maritimae, C. juniperi, C. pinea,C. tujafilina, C. cedri and Cinara sp. from the Eulachnini tribe, Tuberolachnussalignus and Maculolachnus submacula from the Lachninitribe; (d) at least eleven species of the Aphidinae subfamily were found to be infected; nine of them belong to the Aphidini tribe, including three Aphis species (A. fabae, A. nerii and A. hederae ), three samples assigned as Aphis sp., and two Toxoptera species (T. auranti and T. citricidus). The remaining belong to the Macrosiphini tribe, two samples assigned as Cavariella sp., Macrosiphum euphorbiae, Metopolophium dirhodum and Aulacorthum solani; (e) a single infected species belongs to the Chaitophorinae subfamily (Sipha maydis); (f) a single sample of the Eriosomatinae subfamily, Baizongia pistaciae (tribe Fordini), was found to harbour Wolbachia; and (g) a single sample of the Drepanosiphinae subfamily, Neophyllaphis podocarpi, was found to harbour Wolbachia.

[Figure omitted. See PDF.]

Table 1. Aphid populations positive for Wolbachia and PCR amplification results for 16S rDNA, MLST, wsp, gltA and groEL genes.

https://doi.org/10.1371/journal.pone.0028695.t001

It should be noted that at least four different individuals of Cinara pinea (Madeira), Metopolophium dirhodum, Aphis fabae, Aphis hederae, Toxoptera citricidus (Madeira), Sipha maydis and Baizongia pistaciae populations were tested. All individuals were Wolbachia positive. For the rest of the populations, the screening was performed on a pool of four individuals.

Taken together, these results suggest that Wolbachia may be more abundant in aphids than previously thought, and that new universal primers coupled with new sequencing technologies will enable a better detection and investigation of the Wolbachia diversity.

Genotyping aphid Wolbachia strains

The current genotyping of Wolbachia strains is based on MLST approaches [20], [21]. Efforts were made to amplify the MLST genes for the Wolbachia-infected aphid samples; however, the majority of PCRs failed. Only for a few of the samples, some of the genes were successfully amplified (Table 1). Due to these difficulties, attempts were undertaken to characterize the bacterial strains present in each of the thirty-seven Wolbachia-infected aphid populations using a near-full length sequence of the 16S rRNA gene. Additional markers were also used, such as groEL,gltA, wsp and/or other individual MLST gene markers(gatB, coxA, ftsZ, hcpA and fbpA), which could be amplified from the Wolbachia-infected aphid samples.

The results of these efforts can be summarized as follows: (a) a near-full length sequence of the 16S rRNA gene was amplified and analyzed for 35 out of the 37 Wolbachia-infected aphid samples, using PCR-sequencing approaches and the primers as presented in Table 1 and Figure S1. For two samples, GRA4 and GRA40, the amplification of 16S rRNA gene was not possible, and the characterization was based on other genes (see Table 1); (b) genes gatB, coxA, ftsZ, hcpA, fbpA, wsp, gltA and groEL were amplified only from ten, eight, two, two, four, seven, three and five Wolbachia-infected aphid samples, respectively (Table 1); (c) the sequence analysis of gatB, coxA, ftsZ, hcpA, fbpA and wsp revealed the presence of eight, three, two, two, four and two alleles respectively (Tables 2 and 3); (d) the sequence analysis also indicated the presence of novel alleles: seven for gatB, one for ftsZ, one for hcpA and three for fbpA (Tables 2 and 3); (e) gltA and groEL gene fragments were amplified only in three and five Wolbachia-infected aphid samples, respectively (Table 1).

[Figure omitted. See PDF.]

Table 2. Wolbachia MLST allele profiles for positive aphid populations.

https://doi.org/10.1371/journal.pone.0028695.t002

[Figure omitted. See PDF.]

Table 3. Wolbachia WSP HVR profiles for aphid populations.

https://doi.org/10.1371/journal.pone.0028695.t003

These results indicate that there are differences in the Wolbachia infection status among different aphid species and populations and, more importantly, that the currently available genotyping tools of Wolbachia [10], [20], [21], [63], [64]cannot be universally applied for Wolbachia of aphids.

Phylogenetic analysis

Failure to amplify the majority of the MLST and/or other protein coding genes meant that the phylogenetic analysis had to be based mainly on partial 16S rRNA gene sequences (at least 1100 bp). Our Wolbachia sequences appear to cluster in four different groups: few belong either to supergroup A or to supergroup B, while the majority of the sequences form two new clades M and N, distinct from the so far reported supergroups, as depicted in Figure 1. Genetic distances of all the samples of the new groups are more than 2% from the so far reported A to L supergroups (Table S2). The 2% distance is a value necessary for the establishment of a new supergroup [65], [66]. Supergroup M includes 30 new Wolbachia sequences and shows the smallest genetic distance to supergroup B (0.021) and the largest distance to supergroup J (0.059) (Table S3). Supergroup N includes 3 new Wolbachia sequences and shows the smallest genetic distance to supergroup K (0.022) and the largest to supergroup I (0.054) (Table S3). The - within the group - genetic distances of these new groups are only 0.013 and 0.002 for supergroup M and N, respectively, supporting the classification of the strains in new clades. Given the tree topology, presence of recombination events was also examined between M, N and B supergroups, along with Wolbachia strains fromaphids placed into supergroups A and B. No indication of recombination events were detected using the RDP3 package.

[Figure omitted. See PDF.]

Figure 1. Unrooted phylogenetic tree of Wolbachia 16S rRNA gene sequences.

Strains are designated with the names of their host species, followed by the collection site and the sample name. Bayesian posterior probabilities (bottom number) and ML bootstrap values based on 1000 replicates are given. Wolbachia supergroups are shown to the right of the host species names. New supergroups are shaded while aphid Wolbachia strains that belong to supergroup A or B are boxed.

https://doi.org/10.1371/journal.pone.0028695.g001

A thorough phylogenetic analysis based on protein coding genes could not be completed, due to PCR amplification failure in most of the cases (Table 1); however the phylogenetic analysis that was based on the available aphid WolbachiagltA, gatB, fbpA and groEL gene sequences provided several important findings. gltA-based data indicate that the three amplified Wolbachia sequences belong to Supergroup B, while the corresponding 16S rRNA sequences group with the new supergroup M, which has, as stated above, the smallest genetic distance to supergroup B. Interestingly, all three sequences were amplified from C. cedri populations (Figure S2). gatB sequence analysis indicates that almost all amplified sequences group together in a new phylogenetic lineage close to that of supergroup B, except one that groups with supergroup A (Figure S3). A similar picture can be seen with the fbpA-based data with three Wolbachia sequences forming a new phylogenetic cluster and one grouping with supergroup A sequences (Figure S4). For the groEL-based data, one sequence makes a new phylogenetic cluster with supergroup L while the other four group with supergroup B sequences (Figure S5).

Taken together, these results suggest that the aphid fauna may contain an unprecedented range of highly diverse Wolbachia strains, which requires the development of new tools for their detection. In addition, these data clearly indicate the need for the development of new (MLST) tools for the genotyping of Wolbachia strains belonging to new and/or less characterized supergroups.

Discussion

Extending our knowledge on Wolbachia infection of aphids

The presence of Wolbachia was investigated for 425 samples belonging to 153 different species and 70 genera, using a Wolbachia specific 16S rRNA-based PCR assay. The screen included aphid subfamilies with no previous reports of Wolbachia infection and included aphid species from different geographic locations and a variety of hosts (at least 165 different species) [6], [57]–[60]. Despite difficulties with PCR amplification (see below), Wolbachia infections were detected, adding important information to previous studies on aphids which had detected Wolbachia in only four species: three of these species belong to the Aphidinae subfamily and one to the Lachninae. The present analysis showed that the prevalence of Wolbachia infections varied significantly between different aphid populations (Table S1). Wolbachia were detectedin eighteen new aphid species, belonging to the subfamilies Chaitophorinae, Eriosomatinae and Drepanosiphidae, while they were not found in 146 species tested belonging to the seven aphid subfamilies: Aphidinae, Chaitophorinae, Pterocommatinae, Drepanosiphinae, Lachninae, Mindarinae and Eriosomatinae.

A direct comparison with previous screening efforts is difficult since: (a) aphid hosts, sample origin and even screening approaches differed and (b) not many aphid species were common in these studies. Our study confirmed previous results regarding the absence of Wolbachia in members of the subfamily Aphidinae: (i) Acyrthosiphon pisum [58], [59]; (ii) different species of Uroleucon genus [59]; (iii) A. craccivora, Myzuspersicae, Rhopalosiphumpadi, Rhopalosiphummaidis and Schizaphis graminum [59] and (iv) Aphis jacobaeae, Capitophorus carduinis and Sitobium fragariae [57]. It should be noted that Wolbachia was not detected in any species tested of the generaUroleucon, Capitophorus, Myzus and Sitobion although Wolbachia infection was reported in a previous study [61].

Our study also confirmed previous results regarding the absence of Wolbachia in A. craccivora [6]and the presence [60] in all but one C. cedri samples tested (originating from different geographic locations: Spain, Portugal, Iran and Israel). Wolbachia were also detected in five more Cinara species (C. pinea, C. fresai, C. juniperi, C. tujafilina and C. maritimae), suggesting that the genus Cinara has a well-established symbiotic association with Wolbachia. However, it is difficult to speculate about a possible role of Wolbachia in this genus because in 20 out of the 37 samples screened, Wolbachia was not detected. In any case, most members of the Lachninae subfamily harbor S. symbiotica as a second symbiont [67], [68] and, thus the possibility that these species are more prone to accept other infections cannot be ruled out. Finally, the possibility of a co-evolution with the host can be discarded. First, samples from the same species and the same or different location are found in different supergroups (i.e. C. cedri from Israel and Valencia, Spain are found in M and B supergroups; samples from C. pinea are found in M and A supergroup). Second, due to the fact that several of the Cinara species were previously studied in a work analyzing the presence of Serratia in the subfamily Lachninae (94), we can compare the phylogenetic tree obtained in the present work, with those of Buchnera and Serratia previously obtained. The topology obtained regarding the samples from Cinara sp is non-congruent either with Buchnera or with Serratia. A very interesting result is the identification of multiple infections in C. cedri samples. PCR-sequencing analysis of 16S rRNA clones from Israeli populations of C. cedri indicates the presence of two Wolbachia strains: one from supergroup B and a second from the new supergroup M (see below; Figure 1). The fact that DNA was extracted from a mix of four aphids leaves the possibility open that these two strains derive from different individuals.

There are two limitations in our study, regarding the detection of superinfections: the first is the low body mass of many aphid species, which did not allow isolation of high quality and quantity single-aphid DNA for multiple PCRs. The second is the small number of individuals analyzed per population, since we focused on the screening of as many populations as possible, which, in association with the differential abundance of strains and the non-optimized PCR protocols can lead to under-estimation of multiple infections. It should be noted that Wolbachia superinfections have repeatedly been reported in different insect taxa, including Chinese populations of the wheat aphidSitobion miscanthi [61], [69]–[75].

Extending our knowledge on Wolbachia diversity - Two new supergroups

The 16S rRNA gene sequence analysis strongly supports the existence of two new Wolbachia supergroups in aphids and raises questions about the robustness of supergroups E, F and H (Figure 1). Thirty-three of the aphid Wolbachia-specific 16S rRNA gene sequences cluster in two new clades, which are at least 2% genetically distant from all previously described supergroups and from each other (Table S3). However, the analysis also shows that the distance of supergroup A 16S rRNA gene sequences is less than 2% from the sequences present in supergroups E, F and H, suggesting that the overall classification of Wolbachia strains in supergroups (A to N) should be re-evaluated (see Table S3, figures in bold).

Our analysis indicates that the within-supergroup diversity of M and N is 1.3% and 0.2%, respectively (Table S3). Similar phylogenetic analysis with the rest of the genes that are currently being used for the designation of supergroups, could not be completed due to failure of most PCR amplifications. However, the analysis performed with the limited available protein encoding gene sequences (groEL and some MLST genes) also support the presence of new supergroups (Figures S3, S4 and S5).

Earlier efforts to characterize Wolbachia infections were based on the 16S rDNA and ftsZ genes, and later groEL and gltA were included [10], [64], [76]. In 2006, MLST-based systems were proposed for systematic genotyping and strain classification of Wolbachia infections [20], [21]. However, the bacterial strains present in 37 Wolbachia-infected aphid populations, representing 25 aphid species, could not be genotyped using MLST analysis due to failure of PCR amplification despite great effort (Table 1). We managed to obtain sequences from ten samples for gatB, eight for coxA, two for ftsZ, two for hcpA, four forfbpA, seven for wsp, three for gltA and five for groEL (Table 1). Although the sequence analysis in the MLST and wsp databases indicated the presence of new alleles ( Tables 2 and 3), it also clearly shows that the currently available tools cannot be applied universally for the genotyping of the highly diverse aphid Wolbachia strains, and a new MLST system may need to be developed.

The challenge of detection and strain classification of Wolbachia infections in aphids

A major crossroad will be the choice of genes for a new MLST system, given that in the present study there were two cases [see Table 1: Aphis fabae (GRA4) and Metopolophium dirhodum (GRA40)] where Wolbachia-specific amplicons were obtained and confirmed by sequencing analysis, also for some MLST genes, but not for the 16S rRNA gene, which is considered one of the most conserved genes. Our data are in agreement with recent efforts on the assessment of PCR protocols for the detection of Wolbachia,which suggested that the current tools are far from optimal [77].

The development of robust and efficient Wolbachia detection and classification protocols is certainly hindered by the presence of low titre infections and multiple infections [78]–[80]. It has been reported that Wolbachia density may be affected and/or regulated by co-infection with other Wolbachia strains or other vertically transmitted symbionts, as well as by host genotype [81]–[83].

Another important factor is horizontal transfer of Wolbachia genes to host genomes, which further complicates both Wolbachia detection and strain classification. Horizontal transfer events of Wolbachia genome fragments have been reported for several invertebrate species [84]–[88].It is evident that such phenomena can complicate phylogenetic analysis, since nuclear gene copies would evolve in a different way than cytoplasmic copies of Wolbachia genes. Also, Wolbachia detection is compromised in populations that carry nuclear copies of Wolbachia genes but lost the cytoplasmic Wolbachia [87]. The draft genome sequence of the pea aphid Acyrthosiphon pisum revealed the existence of 12 genes of bacterial origin [89], nine of which were intact and closely related to genes of α-Proteobacteria. There is, however, no evidence for horizontal transfer of Wolbachia genes in aphids, and Wolbachia was neither detected in the pea aphid in the present nor in previous studies [58], [59]. The PCR detection of some MLST genes, but not of the 16S rRNA gene in two aphid samples of the present study (Aphis fabae (GRA4) and Metopolophium dirhodum (GRA40)) could be explained by the integration of genomic sequences of a former Wolbachia symbiont into the host genome although alternative causes can not be excluded.

Possible role of Wolbachia in aphids

Aphids feed on phloem sap, which has an unbalanced nitrogen/carbon content and is deficient in a number of nutrients, mainly amino acids, which insects cannot synthesize and are provided by Buchnera aphidicola, their primary endosymbiont. The relationship is mutualistic, since aphids need B. aphidicola for normal growth and reproduction, whereas the bacteriumcannot live outside the aphid [46], [47], [90]. In addition to B. aphidicola, some aphid populations harbor other heritable bacterial symbionts that are not required for host growth and reproduction, referred to as facultative or secondary symbionts [50], [91]. The most common facultative symbionts found in aphids are ‘Ca.Regiella insecticola’, ‘Ca.Hamiltonella defensa’ and ‘Ca.Seratia symbiotica’ [45], [91]. Several studies, mainly in A. pisum, a member of the Aphidinae subfamily, have shown that these symbionts can provide some benefits to the host; however, as mentioned above, no Wolbachia has so far been detected in A. pisum. The genome sequence of these endosymbionts shows that they have lost the ability to synthesize some amino acids and are thus dependent on Buchnera [92]–[94].

C. cedri, a member of the subfamily Lachninae that possess the B. aphidicola with the smallest genome reported so far, and has established a permanent association with the co-primary endosymbiont Serratia symbotica, deserves special attention. Both bacteria are needed for the survival of the whole consortium. When Wolbachia was found in C. cedri, it was postulated that its presence could increase the prevalence of asexual lineages, (C. cedri has a cyclic parthenogenetic life cycle) (see below). In the present study, Wolbachiahas been found in all analyzed C. cedri populations, corroborating their tight association with this species.

Facultative endosymbionts are a common feature of the Lachninae subfamily, to which C. cedri belongs [67], [68]. These symbionts are somehow compensating the drastic metabolic losses that have occurred in B. aphidicola as it has been recently shown for C. tujafilina [95]. The presented data indicate that the members of the Lachninae subfamily tend to be infected with Wolbachia. The possibility that Wolbachia may have a nutritional function in these cases cannot be discarded, as it has been recently proven in thebedbug, Cimex lectularius [96].

Wolbachia is well known for its ability to induce reproductive alterations, such as parthenogenesis, feminization, male-killing and, most commonly, cytoplasmic incompatibility, in its hosts [2], [3]. Aphids are known to have complicated life cycles, which include sexually and asexually reproducing species, as well as species with both sexual and asexual phases [97]. Whether Wolbachia is somehow involved in these phenomena remains to be investigated. Specifically, it would be interesting to check the life cycle of Wolbachia-infected versus non-infected aphids, as its presence could increase the prevalence of asexual lineages, as previously reported for the Hymenopteran group [22].

Conclusion

We report the largest screening effort so far for Wolbachia in aphids. Our results indicate the presence of two new supergroups prevailing in aphids, previously well hidden, probably due to low titer, genetic variability and lack of optimized identification and classification tools. Although Wolbachia was unambiguously identified only in a fraction of the samples analyzed, we believe that its presence is underestimated, and the development of more universal Wolbachia-screening tools is needed. Clarifying the Wolbachia status of aphids can help in the development of novel and environment-friendly methods for the efficient control of aphids, major pests and disease vectors.

Materials and Methods

Sample collection and DNA extraction

Aphid taxa examined in this study, information about their taxonomy, collection locations and the host plants they have been isolated from are listed in Table S1. Natural aphid populations were sampled in different years in Greece (2006, 2007, 2009), Iran (2009, 2010), Israel (2005), Portugal (2009, 2010, 2011) and Spain (2003, 2005, 2009) from a variety of host plants. Aphid species were identified based on morphological criteria [98]–[101] and were stored in 100% ethanol at −20°C. Total DNA of the Greek aphid populations was extracted from single aphids (at least three individuals per sample) while for the Spanish, Portuguese, Israeli and Iranian samples, extractions were done from a pool of four adults. DNA extraction was performed as described previously [102]or by using a modified CTAB protocol [103].

PCR screen

A total of 425 specimens from five subfamilies of the 148 different aphid species were screened for the presence of Wolbachia strains. Detection was based on the amplification of a 16S rRNA gene fragment (438 base pairs) with the Wolbachia specific primers wspecF and wspecR (Figure S1) [5]. For those samples that appeared negative for Wolbachia infection, the quality of DNA was further examined by amplifying part of the mitochondrial 12S rRNA gene (420 base pairs) using primers 12SCFR 5′-GAGAGTGACGGGCGATATGT-3′ and 12SCRR 5′-AAACCAGGATTAGATACCCTATTAT-3′ [104]. PCR amplifications were performed in 20 µl reactions containing 1 µl of DNA, 4 µl 5× reaction buffer (Promega), 1.6 µl MgCl2 (25 mM), 0.1 µl deoxynucleotide triphosphate mixture (25 mM each), 0.5 µl of each primer (25 µM), 0.1 µl of Taq polymerase (Promega, 1 U/µl) and 12.2 µl water. Amplification was performed in a PTC-200 Thermal Cycler (MJ Research), using the following cycling conditions: 95°C for 5 min, followed by 34 cycles of 30 s at 94°C, 30 s at 54°C, 1 min at 72°C and a final extension of 10 min at 72°C. PCR reactions were electrophoresed on a 1.5% agarose gel. Positive samples were further analysed.

PCR, cloning and sequencing of 16S rRNA, groEl, gltA, wsp and MLST gene fragments

Amplification of near full size 16S rRNA sequences proved to be a rather difficult task and required the deployment of a number of approaches (see Figure S1). These involved the use of (a) a new Wolbachia specific primer, W169F, designed for the purposes of this study and the universal eubacterial primer 1513R, followed by a nested PCR using the same forward primer (W169F) and wspecR and (b) the newly designed primer W169F and the new Wolbachia specific primer 16S_woR1 as reverse primer (Figure S1). For some of the populations, a direct PCR with 16S_169F/16S_woR1 was used. PCR amplifications were performed in 20 µl reactions containing 1 µl of DNA, 4 µl 5× reaction buffer (Promega), 1.6 µl MgCl2 (25 mM), 0.1 µl deoxynucleotide triphosphate mixture (25 mM each), 0.5 µl of each primer (25 µM), 0.1 µl of Taq polymerase (Promega 1 U/µl) and 12.2 µl water. Amplification was performed in a PTC-200 Thermal Cycler (MJ Research), using the following cycling conditions: 95°C for 5 min, followed by 34 cycles of 30 s at 94°C, 30 s at 51°C for W169F/1513R and 53°C for W169F/16SwolR1, 1 min at 72°C and a final extension of 10 min at 72°C. The annealing temperature for the nested PCR was 53°C.

The Wolbachia strains of infected aphid populations were genotyped by MLST, wsp, groEL and gtlA based approaches. Gene fragments of the groEL, gtlA,wsp and the MLST genes (gatB, coxA, hcpA, fbpA and ftsZ) were amplified using the respective primers reported previously [17], [20], [64].

Cloning and sequencing

To determine the sequence of 16S rRNA, wsp, groEL, gtlA and MLST gene fragments, PCR fragments were cloned in cases of poor sequencing quality or multiple chromatographic peaks in direct sequencing of PCR products. PCR products from 18 out of the 37 populations harboring Wolbachia were ligated into a T-vector (pGEM-T Easy) and then transformed into DH5α competent cells according to the manufacturer's instructions. Four to six clones were directly subjected to PCR using the primers T7 and SP6. The colony PCR products were purified using the PEG-NaCl method [105] or using NucleoFast® 96 PCR Plates (Macherey-Nagel) according to the manufacturer's instructions. Inserts were fully sequenced with the same primers and with the internal 16S rRNA primer 960R [106]. A dye terminator-labelled cycle sequencing reaction was conducted with the BigDye Terminator v3.1 Cycle Sequencing Kit (PE Applied Biosystems). Reaction products were analysed using an ABI PRISM 310 or an ABI 3730 Genetic Analyzer (PE Applied Biosystems).All Wolbachia gene sequences generated in this study were assembled and manually edited with SeqManII by DNAStar. For each sample, a majority-rule consensus sequence was created.

Nucleotide sequence accession numbers

All 16S rRNA, wsp, groEL, gtlA and MLST gene sequences generated in this study have been deposited in the GenBank database under accession numbers JN384025–JN384106.

Phylogenetic analysis

All Wolbachia 16S rRNA, gatB, fbpA, hcpA, ftsZ, coxA, groEL and gltA gene sequences generated in this study were aligned using MUSCLE [107] and ClustalW [108]. Sequences obtained from GenBank representing all currently known supergroups of Wolbachia were included in the analysis (Table S2). Phylogenetic analyses were performed using maximum-likelihood (ML) and Bayesian methods. PAUP version 4.0b10 was used to select the optimal evolution model by critically evaluating the selected parameters using the Akaike Information Criterion [109]. For the 16S rRNA and gltA gene sequence data the submodel GTR+I+G was selected. For the groEL, gatB and fbpA sequence data, the submodel GTR+G was selected. ML analysis was performed in PAUP using a heuristic search with a random addition of sequences with ten replicates and TBR swapping. The robustness was assessed with 1,000 bootstrap replicates. Bayesian analyses were performed as implemented in MrBayes 3.1 [110]. Analyses were initiated from random starting trees. Four separate runs, each composed of four chains were run for 6,000,000 generations. The cold chain was sampled every 100 generations, and the first 20,000 generations were discarded. Posterior probabilities were computed for the remaining trees.

Recombination events were examined with the default options of the RDP3 software package (Heath et al. 2006).To test for recombination events, we used the RDP3 software package, with all available softwares implemented in it [111]. We used the default options for all analyses.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Position of the primers used in this study, relative to the 16S rRNA gene from wMel.

https://doi.org/10.1371/journal.pone.0028695.s001

(TIF)

Figure S2.

Bayesian inference phylogeny based on gltA data. The three new Wolbachia strains are indicated with bold letters, and the other strains represent supergroups A, B, C, D, F, H, I, and K. Strains are designated with the names of their host species, followed by the collection site and the sample name. Bayesian posterior probabilities (top numbers) and ML bootstrap values based on 100 replicates (bottom numbers) are given.

https://doi.org/10.1371/journal.pone.0028695.s002

(TIF)

Figure S3.

Bayesian inference phylogeny based on gatB data. The 10 new Wolbachia strains are indicated with bold letters, and the other strains represent supergroups A, B, D, F, and H. Strains are designated with the names of their host species, followed by the collection site and the sample name. Bayesian posterior probabilities (top numbers) and ML bootstrap values based on 100 replicates (bottom numbers) are given.

https://doi.org/10.1371/journal.pone.0028695.s003

(TIF)

Figure S4.

Bayesian inference phylogeny based on fbpA data. The four new Wolbachia strains are indicated with bold letters, and the other strains represent supergroups A, B, D, and F. Strains are designated with the names of their host species, followed by the collection site and the sample name. Bayesian posterior probabilities (top numbers) and ML bootstrap values based on 100 replicates (bottom numbers) are given.

https://doi.org/10.1371/journal.pone.0028695.s004

(TIF)

Figure S5.

Bayesian inference phylogeny based on groEL data. The five new Wolbachia strains are indicated with bold letters and the other strains represent supergroups A, B, C, D, F, H, I, K, and L. Strains are designated with the names of their host species, followed by the collection site and the sample name. Bayesian posterior probabilities (top numbers) and ML bootstrap values based on 100 replicates (bottom numbers) are given.

https://doi.org/10.1371/journal.pone.0028695.s005

(TIF)

Table S1.

Wolbachia detection of all aphid populations examined in this study, based on 16S rDNA gene sequencing.

https://doi.org/10.1371/journal.pone.0028695.s006

(DOC)

Table S2.

Taxonomic details of Wolbachia hosts and accession numbers of analyzed sequences.

https://doi.org/10.1371/journal.pone.0028695.s007

(DOC)

Table S3.

Estimates of Evolutionary Divergence (average) over Sequence Pairs between and within Wolbachia Supergroups.

https://doi.org/10.1371/journal.pone.0028695.s008

(DOC)

Acknowledgments

We would like to thank: Prof. J.M. Michelena for his help in collecting and identifying aphid species.

Author Contributions

Conceived and designed the experiments: KB AL MK GT. Performed the experiments: AAA DSG ED MM AP MS VD SR AFA GT. Analyzed the data: AAA DSG MM AL MK GT KB. Contributed reagents/materials/analysis tools: KB AL MK AP AFA PAVB GT. Wrote the paper: KB AAA AL MK GT.

Citation: Augustinos AA, Santos-Garcia D, Dionyssopoulou E, Moreira M, Papapanagiotou A, Scarvelakis M, et al. (2011) Detection and Characterization of Wolbachia Infections in Natural Populations of Aphids: Is the Hidden Diversity Fully Unraveled? PLoS ONE6(12): e28695. https://doi.org/10.1371/journal.pone.0028695

References

1. Werren JH (1997) Biology of Wolbachia. Annu Rev Entomol 42: 587–609.JH Werren1997Biology of Wolbachia.Annu Rev Entomol42587609

2. Aguiar FM, Ilharco FA (2001) Aphids (Homoptera: Aphidoidea) from Madeira Island – New records and corrections. Boletin sanidade Vegetal, Plagas 27: 323–336.FM AguiarFA Ilharco2001Aphids (Homoptera: Aphidoidea) from Madeira Island – New records and corrections.Boletin sanidade Vegetal, Plagas27323336

3. Saridaki A, Bourtzis K (2010) Wolbachia: more than just a bug in insects genitals. Curr Opin Microbiol 13: 67–72.A. SaridakiK. Bourtzis2010Wolbachia: more than just a bug in insects genitals.Curr Opin Microbiol136772

4. Bandi C, Anderson TJ, Genchi C, Blaxter ML (1998) Phylogeny of Wolbachia in filarial nematodes. Proc Biol Sci 265: 2407–2413.C. BandiTJ AndersonC. GenchiML Blaxter1998Phylogeny of Wolbachia in filarial nematodes.Proc Biol Sci26524072413

5. Werren JH, Windsor DM (2000) Wolbachia infection frequencies in insects: evidence of a global equilibrium? Proc Biol Sci 267: 1277–1285.JH WerrenDM Windsor2000Wolbachia infection frequencies in insects: evidence of a global equilibrium?Proc Biol Sci26712771285

6. Jeyaprakash A, Hoy MA (2000) Long PCR improves Wolbachia DNA amplification: wsp sequences found in 76% of sixty-three arthropod species. Insect Mol Biol 9: 393–405.A. JeyaprakashMA Hoy2000Long PCR improves Wolbachia DNA amplification: wsp sequences found in 76% of sixty-three arthropod species.Insect Mol Biol9393405

7. Hilgenboecker K, Hammerstein P, Schlattmann P, Telschow A, Werren JH (2008) How many species are infected with Wolbachia? statistical analysis of current data. FEMS Microbiol Lett 281: 215–220.K. HilgenboeckerP. HammersteinP. SchlattmannA. TelschowJH Werren2008How many species are infected with Wolbachia? statistical analysis of current data.FEMS Microbiol Lett281215220

8. Ferri E, Bain O, Barbuto M, Martin C, Lo N, et al. (2011) New insights into the evolution of Wolbachia infections in filarial nematodes inferred from a large range of screened species. PLoS One 6: e20843.E. FerriO. BainM. BarbutoC. MartinN. Lo2011New insights into the evolution of Wolbachia infections in filarial nematodes inferred from a large range of screened species.PLoS One6e20843

9. Lo N, Paraskevopoulos C, Bourtzis K, O'Neill SL, Werren JH, et al. (2007) Taxonomic status of the intracellular bacterium Wolbachia pipientis. Int J Syst Evol Microbiol 57: 654–657.N. LoC. ParaskevopoulosK. BourtzisSL O'NeillJH Werren2007Taxonomic status of the intracellular bacterium Wolbachia pipientis.Int J Syst Evol Microbiol57654657

10. Werren JH, Zhang W, Guo LR (1995) Evolution and phylogeny of Wolbachia: reproductive parasites of arthropods. Proc Biol Sci 261: 55–63.JH WerrenW. ZhangLR Guo1995Evolution and phylogeny of Wolbachia: reproductive parasites of arthropods.Proc Biol Sci2615563

11. Vandekerckhove TT, Watteyne S, Willems A, Swings JG, Mertens J, et al. (1999) Phylogenetic analysis of the 16S rDNA of the cytoplasmic bacterium Wolbachia from the novel host Folsomia candida (Hexapoda, Collembola) and its implications for wolbachial taxonomy. FEMS Microbiol Lett 180: 279–286.TT VandekerckhoveS. WatteyneA. WillemsJG SwingsJ. Mertens1999Phylogenetic analysis of the 16S rDNA of the cytoplasmic bacterium Wolbachia from the novel host Folsomia candida (Hexapoda, Collembola) and its implications for wolbachial taxonomy.FEMS Microbiol Lett180279286

12. Lo N, Casiraghi M, Salati E, Bazzocchi C, Bandi C (2002) How many wolbachia supergroups exist? Mol Biol Evol 19: 341–346.N. LoM. CasiraghiE. SalatiC. BazzocchiC. Bandi2002How many wolbachia supergroups exist?Mol Biol Evol19341346

13. Vaishampayan PA, Dhotre DP, Gupta RP, Lalwani P, Ghate H, et al. (2007) Molecular evidence and phylogenetic affiliations of Wolbachia in cockroaches. Mol Phylogenet Evol 44: 1346–1351.PA VaishampayanDP DhotreRP GuptaP. LalwaniH. Ghate2007Molecular evidence and phylogenetic affiliations of Wolbachia in cockroaches.Mol Phylogenet Evol4413461351

14. Dunn AK, Stabb EV (2005) Culture-independent characterization of the microbiota of the ant lion Myrmeleon mobilis (Neuroptera: Myrmeleontidae). Appl Environ Microbiol 71: 8784–8794.AK DunnEV Stabb2005Culture-independent characterization of the microbiota of the ant lion Myrmeleon mobilis (Neuroptera: Myrmeleontidae).Appl Environ Microbiol7187848794

15. Bordenstein S, Rosengaus RB (2005) Discovery of a novel Wolbachia super group in Isoptera. Curr Microbiol 51: 393–398.S. BordensteinRB Rosengaus2005Discovery of a novel Wolbachia super group in Isoptera.Curr Microbiol51393398

16. Haegeman A, Vanholme B, Jacob J, Vandekerckhove TT, Claeys M, et al. (2009) An endosymbiotic bacterium in a plant-parasitic nematode: member of a new Wolbachia supergroup. Int J Parasitol 39: 1045–1054.A. HaegemanB. VanholmeJ. JacobTT VandekerckhoveM. Claeys2009An endosymbiotic bacterium in a plant-parasitic nematode: member of a new Wolbachia supergroup.Int J Parasitol3910451054

17. Ros VI, Fleming VM, Feil EJ, Breeuwer JA (2009) How diverse is the genus Wolbachia? Multiple-gene sequencing reveals a putatively new Wolbachiasupergroup recovered from spider mites (Acari: Tetranychidae). Appl Environ Microbiol 75: 1036–1043.VI RosVM FlemingEJ FeilJA Breeuwer2009How diverse is the genus Wolbachia? Multiple-gene sequencing reveals a putatively new Wolbachiasupergroup recovered from spider mites (Acari: Tetranychidae).Appl Environ Microbiol7510361043

18. Rowley SM, Raven RJ, McGraw EA (2004) Wolbachia pipientis in Australian spiders. Curr Microbiol 49: 208–214.SM RowleyRJ RavenEA McGraw2004Wolbachia pipientis in Australian spiders.Curr Microbiol49208214

19. Baldo L, Prendini L, Corthals A, Werren JH (2007) Wolbachia are present in southern african scorpions and cluster with supergroup F. Curr Microbiol 55: 367–373.L. BaldoL. PrendiniA. CorthalsJH Werren2007Wolbachia are present in southern african scorpions and cluster with supergroup F.Curr Microbiol55367373

20. Baldo L, Dunning Hotopp JC, Jolley KA, Bordenstein SR, Biber SA, et al. (2006) Multilocus sequence typing system for the endosymbiont Wolbachia pipientis. Appl Environ Microbiol 72: 7098–7110.L. BaldoJC Dunning HotoppKA JolleySR BordensteinSA Biber2006Multilocus sequence typing system for the endosymbiont Wolbachia pipientis.Appl Environ Microbiol7270987110

21. Paraskevopoulos C, Bordenstein SR, Wernegreen JJ, Werren JH, Bourtzis K (2006) Toward a Wolbachia multilocus sequence typing system: discrimination of Wolbachia strains present in Drosophila species. Curr Microbiol 53: 388–395.C. ParaskevopoulosSR BordensteinJJ WernegreenJH WerrenK. Bourtzis2006Toward a Wolbachia multilocus sequence typing system: discrimination of Wolbachia strains present in Drosophila species.Curr Microbiol53388395

22. Stouthamer R, Breeuwer JA, Hurst GD (1999) Wolbachia pipientis: microbial manipulator of arthropod reproduction. Annu Rev Microbiol 53: 71–102.R. StouthamerJA BreeuwerGD Hurst1999Wolbachia pipientis: microbial manipulator of arthropod reproduction.Annu Rev Microbiol5371102

23. Hancock PA, Sinkins SP, Godfray HC (2011) Population dynamic models of the spread of Wolbachia. Am Nat 177: 323–333.PA HancockSP SinkinsHC Godfray2011Population dynamic models of the spread of Wolbachia.Am Nat177323333

24. Bourtzis K, Dobson SL, Braig HR, O'Neill SL (1998) Rescuing Wolbachia have been overlooked. Nature 391: 852–853.K. BourtzisSL DobsonHR BraigSL O'Neill1998Rescuing Wolbachia have been overlooked.Nature391852853

25. Bourtzis K, Robinson AS (2006) Insect pest control using Wolbachia and/or radiation. In: Bourtzis K, Miller T, editors. Insect Symbiosis 2. CRC Press, Talylor and Francis Group, LLC, Florida, USA. pp. 225–246.K. BourtzisAS Robinson2006Insect pest control using Wolbachia and/or radiation.K. BourtzisT. MillerInsect Symbiosis 2CRC Press, Talylor and Francis Group, LLC, Florida, USA225246

26. Bourtzis K (2008) Wolbachia-based technologies for insect pest population control. Adv Exp Med Biol 627: 104–113.K. Bourtzis2008Wolbachia-based technologies for insect pest population control.Adv Exp Med Biol627104113

27. Laven H (1967) Eradication of Culex pipiens fatigans through cytoplasmic incompatibility. Nature 216: 383–384.H. Laven1967Eradication of Culex pipiens fatigans through cytoplasmic incompatibility.Nature216383384

28. Zabalou S, Riegler M, Theodorakopoulou M, Stauffer C, Savakis C, et al. (2004) Wolbachia-induced cytoplasmic incompatibility as a means for insect pest population control. Proc Natl Acad Sci U S A 101: 15042–15045.S. ZabalouM. RieglerM. TheodorakopoulouC. StaufferC. Savakis2004Wolbachia-induced cytoplasmic incompatibility as a means for insect pest population control.Proc Natl Acad Sci U S A1011504215045

29. Xi Z, Dean JL, Khoo C, Dobson SL (2005) Generation of a novel Wolbachia infection in Aedes albopictus (Asian tiger mosquito) via embryonic microinjection. Insect Biochem Mol Biol 35: 903–910.Z. XiJL DeanC. KhooSL Dobson2005Generation of a novel Wolbachia infection in Aedes albopictus (Asian tiger mosquito) via embryonic microinjection.Insect Biochem Mol Biol35903910

30. Kambris Z, Cook PE, Phuc HK, Sinkins SP (2009) Immune activation by life-shortening Wolbachia and reduced filarial competence in mosquitoes. Science 326: 134–136.Z. KambrisPE CookHK PhucSP Sinkins2009Immune activation by life-shortening Wolbachia and reduced filarial competence in mosquitoes.Science326134136

31. Turley AP, Moreira LA, O'Neill SL, McGraw EA (2009) Wolbachia infection reduces blood-feeding success in the dengue fever mosquito, Aedes aegypti. PLoS Negl Trop Dis 3: e516.AP TurleyLA MoreiraSL O'NeillEA McGraw2009Wolbachia infection reduces blood-feeding success in the dengue fever mosquito, Aedes aegypti.PLoS Negl Trop Dis3e516

32. Zabalou S, Apostolaki A, Livadaras I, Franz G, Robinson AS, et al. (2009) Incompatible insect technique: incompatible males from a Ceratitis capitata genetic sexing strain. Entomologia Experimentalis et Applicata 132: 232–240.S. ZabalouA. ApostolakiI. LivadarasG. FranzAS Robinson2009Incompatible insect technique: incompatible males from a Ceratitis capitata genetic sexing strain.Entomologia Experimentalis et Applicata132232240

33. Apostolaki A, Livadaras I, Saridaki A, Chrysargyris A, Savakis C, et al. (2011) Transinfection of the olive fruit fly Bactrocera oleae with Wolbachia: towards a symbiont-based population control strategy. Journal of Applied Entomology 135: 546–553.A. ApostolakiI. LivadarasA. SaridakiA. ChrysargyrisC. Savakis2011Transinfection of the olive fruit fly Bactrocera oleae with Wolbachia: towards a symbiont-based population control strategy.Journal of Applied Entomology135546553

34. McMeniman CJ, Lane RV, Cass BN, Fong AW, Sidhu M, et al. (2009) Stable introduction of a life-shortening Wolbachia infection into the mosquito Aedes aegypti. Science 323: 141–144.CJ McMenimanRV LaneBN CassAW FongM. Sidhu2009Stable introduction of a life-shortening Wolbachia infection into the mosquito Aedes aegypti.Science323141144

35. Moreira LA, Iturbe-Ormaetxe I, Jeffery JA, Lu G, Pyke AT, et al. (2009) A Wolbachia symbiont in Aedes aegypti limits infection with dengue, Chikungunya and Plasmodium. Cell 139: 1268–1278.LA MoreiraI. Iturbe-OrmaetxeJA JefferyG. LuAT Pyke2009A Wolbachia symbiont in Aedes aegypti limits infection with dengue, Chikungunya and Plasmodium.Cell13912681278

36. Hedges LM, Brownlie JC, O'Neill SL, Johnson KN (2008) Wolbachia and virus protection in insects. Science 322: 702.LM HedgesJC BrownlieSL O'NeillKN Johnson2008Wolbachia and virus protection in insects.Science322702

37. Teixeira L, Ferreira A, Ashburner M (2008) The bacterial symbiont Wolbachia induces resistance to RNA viral infections in Drosophila melanogaster. PLoS Biol 6: e2.L. TeixeiraA. FerreiraM. Ashburner2008The bacterial symbiont Wolbachia induces resistance to RNA viral infections in Drosophila melanogaster.PLoS Biol6e2

38. Bian G, Xu Y, Lu P, Xie Y, Xi Z (2010) The endosymbiotic bacterium Wolbachia induces resistance to dengue virus in Aedes aegypti. PLoS Pathog 6: e1000833.G. BianY. XuP. LuY. XieZ. Xi2010The endosymbiotic bacterium Wolbachia induces resistance to dengue virus in Aedes aegypti.PLoS Pathog6e1000833

39. Cook PE, McGraw EA (2010) Wolbachia pipientis: an expanding bag of tricks to explore for disease control. Trends Parasitol 26: 373–375.PE CookEA McGraw2010Wolbachia pipientis: an expanding bag of tricks to explore for disease control.Trends Parasitol26373375

40. Glaser RL, Meola MA (2010) The native Wolbachia endosymbionts of Drosophila melanogaster and Culex quinquefasciatus increase host resistance to West Nile virus infection. PLoS One 5: e11977.RL GlaserMA Meola2010The native Wolbachia endosymbionts of Drosophila melanogaster and Culex quinquefasciatus increase host resistance to West Nile virus infection.PLoS One5e11977

41. Dixon AFG, Kindlmann P, Leps J, Holman J (1987) Why There are So Few Species of Aphids, Especially in the Tropics. The American Naturalist 129: 580–592.AFG DixonP. KindlmannJ. LepsJ. Holman1987Why There are So Few Species of Aphids, Especially in the Tropics.The American Naturalist129580592

42. Schadler M, Brandl R, Haase J (2007) Antagonistic interactions between plant competition and insect herbivory. Ecology 88: 1490–1498.M. SchadlerR. BrandlJ. Haase2007Antagonistic interactions between plant competition and insect herbivory.Ecology8814901498

43. Brault V, Tanguy S, Reinbold C, Le Trionnaire G, Arneodo J, et al. (2010) Transcriptomic analysis of intestinal genes following acquisition of pea enation mosaic virus by the pea aphid Acyrthosiphon pisum. J Gen Virol 91: 802–808.V. BraultS. TanguyC. ReinboldG. Le TrionnaireJ. Arneodo2010Transcriptomic analysis of intestinal genes following acquisition of pea enation mosaic virus by the pea aphid Acyrthosiphon pisum.J Gen Virol91802808

44. Douglas AE (2006) Phloem-sap feeding by animals: Problems and solutions. Journal of Experimental Botany 57: 747–754.AE Douglas2006Phloem-sap feeding by animals: Problems and solutions.Journal of Experimental Botany57747754

45. Oliver KM, Degnan PH, Burke GR, Moran NA (2010) Facultative Symbionts in Aphids and the Horizontal Transfer of Ecologically Important Traits. Annual Review of Entomology 55: 247–266.KM OliverPH DegnanGR BurkeNA Moran2010Facultative Symbionts in Aphids and the Horizontal Transfer of Ecologically Important Traits.Annual Review of Entomology55247266

46. Douglas AE (1997) Parallels and contrasts between symbiotic bacteria and bacterial-derived organelles: evidence from Buchnera, the bacterial symbiont of aphids. FEMS MicrobiologyEcology 24: 1–9.AE Douglas1997Parallels and contrasts between symbiotic bacteria and bacterial-derived organelles: evidence from Buchnera, the bacterial symbiont of aphids.FEMS MicrobiologyEcology2419

47. Baumann P (2005) Biology bacteriocyte-associated endosymbionts of plant sap-sucking insects. Annu Rev Microbiol 59: 155–189.P. Baumann2005Biology bacteriocyte-associated endosymbionts of plant sap-sucking insects.Annu Rev Microbiol59155189

48. Moya A, Pereto J, Gil R, Latorre A (2008) Learning how to live together: genomic insights into prokaryote-animal symbioses. Nat Rev Genet 9: 218–229.A. MoyaJ. PeretoR. GilA. Latorre2008Learning how to live together: genomic insights into prokaryote-animal symbioses.Nat Rev Genet9218229

49. Brinza L, Vinuelas J, Cottret L, Calevro F, Rahbe Y, et al. (2009) Systemic analysis of the symbiotic function of Buchnera aphidicola, the primary endosymbiont of the pea aphid Acyrthosiphon pisum. C R Biol 332: 1034–1049.L. BrinzaJ. VinuelasL. CottretF. CalevroY. Rahbe2009Systemic analysis of the symbiotic function of Buchnera aphidicola, the primary endosymbiont of the pea aphid Acyrthosiphon pisum.C R Biol33210341049

50. Oliver KM, Moran NA, Hunter MS (2005) Variation in resistance to parasitism in aphids is due to symbionts not host genotype. Proc Natl Acad Sci U S A 102: 12795–12800.KM OliverNA MoranMS Hunter2005Variation in resistance to parasitism in aphids is due to symbionts not host genotype.Proc Natl Acad Sci U S A1021279512800

51. Vorburger C, Gehrer L, Rodriguez P (2010) A strain of the bacterial symbiont Regiella insecticola protects aphids against parasitoids. Biol Lett 6: 109–111.C. VorburgerL. GehrerP. Rodriguez2010A strain of the bacterial symbiont Regiella insecticola protects aphids against parasitoids.Biol Lett6109111

52. Montllor CB, Maxmen A, Purcell AH (2002) Facultative bacterial endosymbionts benefit pea aphids Acyrthosiphon pisum under heat stress. Ecological Entomology 27: 189–195.CB MontllorA. MaxmenAH Purcell2002Facultative bacterial endosymbionts benefit pea aphids Acyrthosiphon pisum under heat stress.Ecological Entomology27189195

53. Moran NA, Jarvik T (2010) Lateral transfer of genes from fungi underlies carotenoid production in aphids. Science 328: 624–627.NA MoranT. Jarvik2010Lateral transfer of genes from fungi underlies carotenoid production in aphids.Science328624627

54. Nikoh N, McCutcheon JP, Kudo T, Miyagishima SY, Moran NA, et al. (2010) Bacterial genes in the aphid genome: absence of functional gene transfer from Buchnera to its host. PLoS Genet 6: e1000827.N. NikohJP McCutcheonT. KudoSY MiyagishimaNA Moran2010Bacterial genes in the aphid genome: absence of functional gene transfer from Buchnera to its host.PLoS Genet6e1000827

55. Perez-Brocal V, Gil R, Ramos S, Lamelas A, Postigo M, et al. (2006) A small microbial genome: the end of a long symbiotic relationship? Science 314: 312–313.V. Perez-BrocalR. GilS. RamosA. LamelasM. Postigo2006A small microbial genome: the end of a long symbiotic relationship?Science314312313

56. Gosalbes MJ, Lamelas A, Moya A, Latorre A (2008) The striking case of tryptophan provision in the cedar aphid Cinara cedri. J Bacteriol 190: 6026–6029.MJ GosalbesA. LamelasA. MoyaA. Latorre2008The striking case of tryptophan provision in the cedar aphid Cinara cedri.J Bacteriol19060266029

57. West SA, Cook JM, Werren JH, Godfray HC (1998) Wolbachia in two insect host-parasitoid communities. Mol Ecol 7: 1457–1465.SA WestJM CookJH WerrenHC Godfray1998Wolbachia in two insect host-parasitoid communities.Mol Ecol714571465

58. Tsuchida T, Koga R, Shibao H, Matsumoto T, Fukatsu T (2002) Diversity and geographic distribution of secondary endosymbiotic bacteria in natural populations of the pea aphid, Acyrthosiphon pisum. Mol Ecol 11: 2123–2135.T. TsuchidaR. KogaH. ShibaoT. MatsumotoT. Fukatsu2002Diversity and geographic distribution of secondary endosymbiotic bacteria in natural populations of the pea aphid, Acyrthosiphon pisum.Mol Ecol1121232135

59. Nirgianaki A, Banks GK, Frohlich DR, Veneti Z, Braig HR, et al. (2003) Wolbachia infections of the whitefly Bemisia tabaci. Curr Microbiol 47: 93–101.A. NirgianakiGK BanksDR FrohlichZ. VenetiHR Braig2003Wolbachia infections of the whitefly Bemisia tabaci.Curr Microbiol4793101

60. Gomez-Valero L, Soriano-Navarro M, Perez-Brocal V, Heddi A, Moya A, et al. (2004) Coexistence of Wolbachia with Buchnera aphidicola and a secondary symbiont in the aphid Cinara cedri. J Bacteriol 186: 6626–6633.L. Gomez-ValeroM. Soriano-NavarroV. Perez-BrocalA. HeddiA. Moya2004Coexistence of Wolbachia with Buchnera aphidicola and a secondary symbiont in the aphid Cinara cedri.J Bacteriol18666266633

61. Wang Z, Shen ZR, Song Y, Liu HY, Li ZX (2009) Distribution and diversity of Wolbachia in different populations of the wheat aphid Sitobion miscanthi (Hemiptera: Aphididae) in China. European Journal of Entomology 106: 49–55.Z. WangZR ShenY. SongHY LiuZX Li2009Distribution and diversity of Wolbachia in different populations of the wheat aphid Sitobion miscanthi (Hemiptera: Aphididae) in China.European Journal of Entomology1064955

62. Ortiz-Rivas B, Martinez-Torres D (2010) Combination of molecular data support the existence of three main lineages in the phylogeny of aphids (Hemiptera: Aphididae) and the basal position of the subfamily Lachninae. Mol Phylogenet Evol 55: 305–317.B. Ortiz-RivasD. Martinez-Torres2010Combination of molecular data support the existence of three main lineages in the phylogeny of aphids (Hemiptera: Aphididae) and the basal position of the subfamily Lachninae.Mol Phylogenet Evol55305317

63. Zhou W, Rousset F, O'Neil S (1998) Phylogeny and PCR-based classification of Wolbachia strains using wsp gene sequences. Proc Biol Sci 265: 509–515.W. ZhouF. RoussetS. O'Neil1998Phylogeny and PCR-based classification of Wolbachia strains using wsp gene sequences.Proc Biol Sci265509515

64. Casiraghi M, Bordenstein SR, Baldo L, Lo N, Beninati T, et al. (2005) Phylogeny of Wolbachia pipientis based on gltA, groEL and ftsZ gene sequences: clustering of arthropod and nematode symbionts in the F supergroup, and evidence for further diversity in the Wolbachia tree. Microbiology 151: 4015–4022.M. CasiraghiSR BordensteinL. BaldoN. LoT. Beninati2005Phylogeny of Wolbachia pipientis based on gltA, groEL and ftsZ gene sequences: clustering of arthropod and nematode symbionts in the F supergroup, and evidence for further diversity in the Wolbachia tree.Microbiology15140154022

65. Breeuwer JA, Stouthamer R, Barns SM, Pelletier DA, Weisburg WG, et al. (1992) Phylogeny of cytoplasmic incompatibility micro-organisms in the parasitoid wasp genus Nasonia (Hymenoptera: Pteromalidae) based on 16S ribosomal DNA sequences. Insect Mol Biol 1: 25–36.JA BreeuwerR. StouthamerSM BarnsDA PelletierWG Weisburg1992Phylogeny of cytoplasmic incompatibility micro-organisms in the parasitoid wasp genus Nasonia (Hymenoptera: Pteromalidae) based on 16S ribosomal DNA sequences.Insect Mol Biol12536

66. Stouthamer R, Breeuwert JA, Luck RF, Werren JH (1993) Molecular identification of microorganisms associated with parthenogenesis. Nature 361: 66–68.R. StouthamerJA BreeuwertRF LuckJH Werren1993Molecular identification of microorganisms associated with parthenogenesis.Nature3616668

67. Burke GR, Normark BB, Favret C, Moran NA (2009) Evolution and diversity of facultative symbionts from the aphid subfamily Lachninae. Appl Environ Microbiol 75: 5328–5335.GR BurkeBB NormarkC. FavretNA Moran2009Evolution and diversity of facultative symbionts from the aphid subfamily Lachninae.Appl Environ Microbiol7553285335

68. Lamelas A, Gosalbes MJ, Moya A, Latorre A (2011) New Clues about the Evolutionary History of Metabolic Losses in Bacterial Endosymbionts, Provided by the Genome of Buchnera aphidicola from the Aphid Cinara tujafilina. Appl Environ Microbiol 77: 4446–4454.A. LamelasMJ GosalbesA. MoyaA. Latorre2011New Clues about the Evolutionary History of Metabolic Losses in Bacterial Endosymbionts, Provided by the Genome of Buchnera aphidicola from the Aphid Cinara tujafilina.Appl Environ Microbiol7744464454

69. Sinkins SP, Braig HR, O'Neill SL (1995) Wolbachia superinfections and the expression of cytoplasmic incompatibility. Proc Biol Sci 261: 325–330.SP SinkinsHR BraigSL O'Neill1995Wolbachia superinfections and the expression of cytoplasmic incompatibility.Proc Biol Sci261325330

70. Dobson SL, Marsland EJ, Veneti Z, Bourtzis K, O'Neill SL (2002) Characterization of Wolbachia host cell range via the in vitro establishment of infections. Appl Environ Microbiol 68: 656–660.SL DobsonEJ MarslandZ. VenetiK. BourtzisSL O'Neill2002Characterization of Wolbachia host cell range via the in vitro establishment of infections.Appl Environ Microbiol68656660

71. Reuter M, Keller L (2003) High levels of multiple Wolbachia infection and recombination in the ant Formica exsecta. Mol Biol Evol 20: 748–753.M. ReuterL. Keller2003High levels of multiple Wolbachia infection and recombination in the ant Formica exsecta.Mol Biol Evol20748753

72. Kikuchi Y, Fukatsu T (2003) Diversity of Wolbachia endosymbionts in heteropteran bugs. Appl Environ Microbiol 69: 6082–6090.Y. KikuchiT. Fukatsu2003Diversity of Wolbachia endosymbionts in heteropteran bugs.Appl Environ Microbiol6960826090

73. Duron O, Lagnel J, Raymond M, Bourtzis K, Fort P, et al. (2005) Transposable element polymorphism of Wolbachia in the mosquito Culex pipiens: evidence of genetic diversity, superinfection and recombination. Mol Ecol 14: 1561–1573.O. DuronJ. LagnelM. RaymondK. BourtzisP. Fort2005Transposable element polymorphism of Wolbachia in the mosquito Culex pipiens: evidence of genetic diversity, superinfection and recombination.Mol Ecol1415611573

74. Li ZX, Hu DX, Song Y, Shen ZR (2005) Molecular differentiation of the B biotype from other biotypes of Bemisia tabaci (Hemiptera: Aleyrodidae), based on internally transcribed spacer 1 sequences. European Journal of Entomology 102: 293–297.ZX LiDX HuY. SongZR Shen2005Molecular differentiation of the B biotype from other biotypes of Bemisia tabaci (Hemiptera: Aleyrodidae), based on internally transcribed spacer 1 sequences.European Journal of Entomology102293297

75. Li ZX, Lin HZ, Guo XP (2007) Prevalence of Wolbachia infection in Bemisia tabaci. Curr Microbiol 54: 467–471.ZX LiHZ LinXP Guo2007Prevalence of Wolbachia infection in Bemisia tabaci.Curr Microbiol54467471

76. O'Neill SL, Giordano R, Colbert AM, Karr TL, Robertson HM (1992) 16S rRNA phylogenetic analysis of the bacterial endosymbionts associated with cytoplasmic incompatibility in insects. Proc Natl Acad Sci U S A 89: 2699–2702.SL O'NeillR. GiordanoAM ColbertTL KarrHM Robertson199216S rRNA phylogenetic analysis of the bacterial endosymbionts associated with cytoplasmic incompatibility in insects.Proc Natl Acad Sci U S A8926992702

77. Simoes PM, Mialdea G, Reiss D, Sagot MF, Charlat S (2011) Wolbachia detection: an assessment of standard PCR protocols. Mol Ecol Resour 11: 567–572.PM SimoesG. MialdeaD. ReissMF SagotS. Charlat2011Wolbachia detection: an assessment of standard PCR protocols.Mol Ecol Resour11567572

78. Arthofer W, Riegler M, Schneider D, Krammer M, Miller WJ, et al. (2009) Hidden Wolbachia diversity in field populations of the European cherry fruit fly, Rhagoletis cerasi (Diptera, Tephritidae). Mol Ecol 18: 3816–3830.W. ArthoferM. RieglerD. SchneiderM. KrammerWJ Miller2009Hidden Wolbachia diversity in field populations of the European cherry fruit fly, Rhagoletis cerasi (Diptera, Tephritidae).Mol Ecol1838163830

79. Miller WJ, Ehrman L, Schneider D (2010) Infectious speciation revisited: impact of symbiont-depletion on female fitness and mating behaviour of Drosophila paulistorum. PLoS Pathog e1001214.WJ MillerL. EhrmanD. Schneider2010Infectious speciation revisited: impact of symbiont-depletion on female fitness and mating behaviour of Drosophila paulistorum.PLoS Pathoge1001214

80. Arthofer W, Riegler M, Schuler H, Schneider D, Moder K, et al. (2011) Allele intersection analysis: a novel tool for multi locus sequence assignment in multiply infected hosts. PLoS One 6: e22198.W. ArthoferM. RieglerH. SchulerD. SchneiderK. Moder2011Allele intersection analysis: a novel tool for multi locus sequence assignment in multiply infected hosts.PLoS One6e22198Epub 2011 Jul 15. Epub 2011 Jul 15.

81. Mouton L, Henri H, Bouletreau M, Vavre F (2003) Strain-specific regulation of intracellular Wolbachia density in multiply infected insects. Mol Ecol 12: 3459–3465.L. MoutonH. HenriM. BouletreauF. Vavre2003Strain-specific regulation of intracellular Wolbachia density in multiply infected insects.Mol Ecol1234593465

82. Kondo N, Shimada M, Fukatsu T (2005) Infection density of Wolbachia endosymbiont affected by co-infection and host genotype. Biol Lett 1: 488–491.N. KondoM. ShimadaT. Fukatsu2005Infection density of Wolbachia endosymbiont affected by co-infection and host genotype.Biol Lett1488491

83. Vautrin E, Vavre F (2009) Interactions between vertically transmitted symbionts: cooperation or conflict? Trends Microbiol 17: 95–99.E. VautrinF. Vavre2009Interactions between vertically transmitted symbionts: cooperation or conflict?Trends Microbiol179599

84. Kondo N, Nikoh N, Ijichi N, Shimada M, Fukatsu T (2002) Genome fragment of Wolbachia endosymbiont transferred to X chromosome of host insect. Proc Natl Acad Sci U S A 99: 14280–14285.N. KondoN. NikohN. IjichiM. ShimadaT. Fukatsu2002Genome fragment of Wolbachia endosymbiont transferred to X chromosome of host insect.Proc Natl Acad Sci U S A991428014285

85. Fenn K, Conlon C, Jones M, Quail MA, Holroyd NE, et al. (2006) Phylogenetic relationships of the Wolbachia of nematodes and arthropods. PLoS Pathog 2: e94.K. FennC. ConlonM. JonesMA QuailNE Holroyd2006Phylogenetic relationships of the Wolbachia of nematodes and arthropods.PLoS Pathog2e94

86. Dunning Hotopp JC, Clark ME, Oliveira DC, Foster JM, Fischer P, et al. (2007) Widespread lateral gene transfer from intracellular bacteria to multicellular eukaryotes. Science 317: 1753–1756.JC Dunning HotoppME ClarkDC OliveiraJM FosterP. Fischer2007Widespread lateral gene transfer from intracellular bacteria to multicellular eukaryotes.Science31717531756

87. McNulty SN, Foster JM, Mitreva M, Dunning Hotopp JC, Martin J, et al. (2010) Endosymbiont DNA in endobacteria-free filarial nematodes indicates ancient horizontal genetic transfer. PLoS One 5: e11029.SN McNultyJM FosterM. MitrevaJC Dunning HotoppJ. Martin2010Endosymbiont DNA in endobacteria-free filarial nematodes indicates ancient horizontal genetic transfer.PLoS One5e11029

88. Doudoumis V, Tsiamis G, Wamwiri F, Brelsfoard C, Alam U, et al. (2011) Detection and characterization of Wolbachia infections in laboratory and natural populations of different species of tse-tse flies (genus Glossina). BMC Microbiology. V. DoudoumisG. TsiamisF. WamwiriC. BrelsfoardU. Alam2011Detection and characterization of Wolbachia infections in laboratory and natural populations of different species of tse-tse flies (genus Glossina).BMC MicrobiologyIn press. In press.

89. IAGC (2010) Genome sequence of the pea aphid Acyrthosiphon pisum. PLoS Biol 8: e1000313.IAGC2010Genome sequence of the pea aphid Acyrthosiphon pisum.PLoS Biol8e1000313

90. Buchner P (1965) Endosymbiosis of Animals with Plant Microorganisms. New York: Interscience. P. Buchner1965Endosymbiosis of Animals with Plant MicroorganismsNew YorkInterscience

91. Moran NA, Degnan PH, Santos SR, Dunbar HE, Ochman H (2005) The players in a mutualistic symbiosis: insects, bacteria, viruses, and virulence genes. Proc Natl Acad Sci U S A 102: 16919–16926.NA MoranPH DegnanSR SantosHE DunbarH. Ochman2005The players in a mutualistic symbiosis: insects, bacteria, viruses, and virulence genes.Proc Natl Acad Sci U S A1021691916926

92. Degnan PH, Yu Y, Sisneros N, Wing RA, Moran NA (2009a) Hamiltonella defensa, genome evolution of protective bacterial endosymbiont from pathogenic ancestors. Proc Natl Acad Sci U S A 106: 9063–9068.PH DegnanY. YuN. SisnerosRA WingNA Moran2009aHamiltonella defensa, genome evolution of protective bacterial endosymbiont from pathogenic ancestors.Proc Natl Acad Sci U S A10690639068

93. Degnan PH, Leonardo TE, Cass BN, Hurwitz B, Stern D, et al. (2009b) Dynamics of genome evolution in facultative symbionts of aphids. Environ Microbiol 12: 2060–2069.PH DegnanTE LeonardoBN CassB. HurwitzD. Stern2009bDynamics of genome evolution in facultative symbionts of aphids.Environ Microbiol1220602069

94. Burke GR, Normark BB, Favret C, Moran NA (2009) Evolution and diversity of facultative symbionts from the aphid subfamily Lachninae. Appl Environ Microbiol 75: 5328–5335.GR BurkeBB NormarkC. FavretNA Moran2009Evolution and diversity of facultative symbionts from the aphid subfamily Lachninae.Appl Environ Microbiol7553285335

95. Lamelas A, Gosalbes MJ, Moya A, Latorre A (2011) New Clues about the Evolutionary History of Metabolic Losses in Bacterial Endosymbionts, Provided by the Genome of Buchnera aphidicola from the Aphid Cinara tujafilina. Appl Environ Microbiol 77: 4446–4454.A. LamelasMJ GosalbesA. MoyaA. Latorre2011New Clues about the Evolutionary History of Metabolic Losses in Bacterial Endosymbionts, Provided by the Genome of Buchnera aphidicola from the Aphid Cinara tujafilina.Appl Environ Microbiol7744464454

96. Hosokawa T, Koga R, Kikuchi Y, Meng XY, Fukatsu T (2010) Wolbachia as a bacteriocyte-associated nutritional mutualist. Proc Natl Acad Sci U S A 107: 769–774.T. HosokawaR. KogaY. KikuchiXY MengT. Fukatsu2010Wolbachia as a bacteriocyte-associated nutritional mutualist.Proc Natl Acad Sci U S A107769774

97. Heie OE (2004) The history of the studies on aphid palaeontology and their bearing on the evolutionary history of aphids. Paris: INRA Editions. OE Heie2004The history of the studies on aphid palaeontology and their bearing on the evolutionary history of aphidsParisINRA Editions

98. Aguiar FM, Ilharco FA (2001) Aphids (Homoptera: Aphidoidea) from Madeira Island – New records and corrections. Boletin sanidade Vegetal, Plagas 27: 323–336.FM AguiarFA Ilharco2001Aphids (Homoptera: Aphidoidea) from Madeira Island – New records and corrections.Boletin sanidade Vegetal, Plagas27323336

99. Blackman RL, Eastop VF (1984) Aphids on the world's crops: an identification and information guide. Chichester West Sussex; New York: Wiley. RL BlackmanVF Eastop1984Aphids on the world's crops: an identification and information guideChichester West Sussex; New YorkWiley

100. Blackman RL, Eastop VF (2000) Aphids on the world's crops: an identification and information guide. Chichester, West Sussex, England; New York: Wiley. RL BlackmanVF Eastop2000Aphids on the world's crops: an identification and information guideChichester, West Sussex, England; New YorkWiley

101. Remaudiere G, Remaudiere M (1997) Catalogue des Aphididae du Monde. G. RemaudiereM. Remaudiere1997Catalogue des Aphididae du Monde.Institut National de la Recherche Agronomique: Paris. Institut National de la Recherche Agronomique: Paris.

102. Latorre A, Moya A, Ayala FJ (1986) Evolution of mitochondrial DNA in Drosophila subobscura. Proc Natl Acad Sci U S A 83: 8649–8653.A. LatorreA. MoyaFJ Ayala1986Evolution of mitochondrial DNA in Drosophila subobscura.Proc Natl Acad Sci U S A8386498653

103. Doyle JJ, Doyle JL (1990) Isolation of plant DNA from fresh tissue. Focus 12: 13–15.JJ DoyleJL Doyle1990Isolation of plant DNA from fresh tissue.Focus121315

104. Hanner R, Fugate M (1997) Branchiopod phylogenetic reconstruction from 12S rDNA sequence data. Journal of Crustacean Biology 17: 174–183.R. HannerM. Fugate1997Branchiopod phylogenetic reconstruction from 12S rDNA sequence data.Journal of Crustacean Biology17174183

105. Hartley JL, Bowen H (2003) PEG precipitation for selective removal of small DNA fragments. Focus 25: 18.JL HartleyH. Bowen2003PEG precipitation for selective removal of small DNA fragments.Focus2518

106. Reed DW, Fujita Y, Delwiche ME, Blackwelder DB, Sheridan PP, et al. (2002) Microbial communities from methane hydrate-bearing deep marine sediments in a forearc basin. Appl Environ Microbiol 68: 3759–3770.DW ReedY. FujitaME DelwicheDB BlackwelderPP Sheridan2002Microbial communities from methane hydrate-bearing deep marine sediments in a forearc basin.Appl Environ Microbiol6837593770

107. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32: 1792–1797.RC Edgar2004MUSCLE: multiple sequence alignment with high accuracy and high throughput.Nucleic Acids Res3217921797

108. Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22: 4673–4680.JD ThompsonDG HigginsTJ Gibson1994CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice.Nucleic Acids Res2246734680

109. Swofford DL (2000) PAUP: phylogenetic analysis using parsimony, 4.0, beta version 4a ed. Sunderland, Md.: Sinauer Associates. DL Swofford2000PAUP: phylogenetic analysis using parsimony, 4.0, beta version 4a ed.Sunderland, Md.Sinauer Associates

110. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574.F. RonquistJP Huelsenbeck2003MrBayes 3: Bayesian phylogenetic inference under mixed models.Bioinformatics1915721574

111. Heath L, van der Walt E, Varsani A, Martin DP (2006) Recombination patterns in aphthoviruses mirror those found in other picornaviruses. J Virol 80: 11827–11832.L. HeathE. van der WaltA. VarsaniDP Martin2006Recombination patterns in aphthoviruses mirror those found in other picornaviruses.J Virol801182711832