Genetic evidence from mitochondrial, nuclear, and endosymbiont markers for the evolution of host plant associated species in the aphid genus Hyalopterus (Hemiptera: Aphididae)

Characteristics and Comparative Analysis of Mitochondrial Genomes of the Aphid Genus Hyalopterus Koch (Hemiptera: Aphididae: Aphidinae)

Using Illumina sequencing technology, we generated complete mitochondrial genomes (mitogenomes) of three constituent species of the aphid genus Hyalopterus Koch, Hyalopterus amygdali (Blanchard), Hyalopterus arundiniformis Ghulamullah, and Hyalopterus pruni (Geoffroy). The sizes of the Hyalopterus mitogenomes range from 15,306 to 15,410 bp, primarily due to variations in the length of non-coding regions. The Hyalopterus mitogenomes consist of 37 coding genes arranged in the order of the ancestral insect mitogenome, a control region, and a repeat region between trnE and trnF. According to the COI-based analysis, one previously reported mitogenome of H. pruni should be assigned to H. arundiniformis. The gene order, nucleotide composition, and codon usage in the Hyalopterus mitogenomes are highly conserved and similar to those of other species of Aphidinae. The tandem repeat units differ in nucleotide composition, length, and copy number across three Hyalopterus species. Within the widespread Eurasian species H. arundiniformis, variation in repeat units among different geographic populations is observed, indicating that the repeat region may provide valuable insights for studying the intraspecific diversification of aphids. Phylogenetic analyses based on 28 complete mitogenomes of Aphidinae supported the monophyly of Aphidinae, Aphidini, Macrosiphini, and two subtribes of Aphidini. Hyalopterus was monophyletic. H. amygdali and H. pruni formed a sister group, while H. arundiniformis was placed basally. Characterization of the mitogenomes of Hyalopterus provides valuable resources for further comparative studies and for advancing our understanding of the aphid mitogenome architecture.

Complex evolution in Aphis gossypii group (Hemiptera: Aphididae), evidence of primary host shift and hybridization between sympatric species

Aphids provide a good model system to understand the ecological speciation concept, since the majority of the species are host-specific, and they spend their entire lifecycle on certain groups of host plants. Aphid species that apparently have wide host plant ranges have often turned out to be complexes of host-specialized biotypes. Here we investigated the various host-associated populations of the two recently diverged species, Aphis gossypii and A. rhamnicola, having multiple primary hosts, to understand the complex evolution with host-associated speciation. Using mitochondrial DNA marker and nine microsatellite loci, we reconstructed the haplotype network, and analyzed the genetic structure and relationships. Approximate Bayesian computation was also used to infer the ancestral primary host and host-associated divergence, which resulted in Rhamnus being the most ancestral host for A. gossypii and A. rhamnicola. As a result, Aphis gossypii and A. rhamnicola do not randomly use their primary and secondary host plants; rather, certain biotypes use only some secondary and specific primary hosts. Some biotypes are possibly in a diverging state through specialization to specific primary hosts. Our results also indicate that a new heteroecious race can commonly be derived from the heteroecious ancestor, showing strong evidence of ecological specialization through a primary host shift in both A. gossypii and A. rhamnicola. Interestingly, A. gossypii and A. rhamnicola shared COI haplotypes with each other, thus there is a possibility of introgression by hybridization between them by cross-sharing same primary hosts. Our results contribute to a new perspective in the study of aphid evolution by identifying complex evolutionary trends in the gossypii sensu lato complex.

The complete mitochondrial genome of the mealy plum aphid, Hyalopterus pruni (Hemiptera: Aphididae)

This study completes the sequencing and annotation of the mitochondrial genome (mitogenome) of Hyalopterus pruni (Hemiptera: Aphididae) by using the high-throughput sequencing. The mitogenome is a typical circular DNA of 15,410 bp with 86.2% A + T content, and consists of 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, a repeat region between tRNA-Glu and tRNA-Phe, and a control region. The gene order follows the putative ancestral arrangements of insect mitogenome. All 13 protein-coding genes start with codon ATN and terminate with TAA or a single T. All tRNA genes have typical clover-leaf structure except for tRNA-Ser^AGN. The control region is 638 bp in length with 86.0% A + T content. The phylogenetic tree supports the monophyly of Aphidini and Macrosiphini in Aphidinae and the sister relationship between Hyalopterus pruni and Schizaphis graminum.

Natural Occurrence of Secondary Bacterial Symbionts in Aphids from Tunisia, with a Focus on Genus Hyalopterus

Aphids (Hemiptera: Aphididae) can harbor two types of bacterial symbionts. In addition to the obligate endosymbiont Buchnera aphidicola Munson, Baumann and Kinsey 1991 (Enterobacteriales: Enterobacteriaceae), several facultative symbiotic bacteria, called secondary (S) symbionts, have been identified among many important pest aphid species. To determine interpopulational diversity of S-symbionts, we carried out a survey in a total of 18 populations of six aphid species collected from six localities in Tunisia, by performing a diagnostic polymerase chain reaction analysis of partial 16S-23S rRNA operon sequences. While 61.7% of individuals contained only Buchnera, three S-symbionts were found at different frequencies. Arsenophonus sp. Gherna et al. 1991 (Enterobacteriales: Enterobacteriaceae) was found in all species under study except for Acyrtosiphon pisum (Harris 1776) (Aphidinae: Macrosiphini); Serratia symbiotica Moran et al. 2005 (Enterobacteriales: Enterobacteriaceae) was present in all analyzed individuals of A. pisum but only sporadically in Aphis spiraecola (Patch 1914) (Aphidinae: Aphidini) and Hyalopterus amygdali (Blanchard 1840) (Aphidinae: Aphidini), while Hamiltonella defensa Moran et al. 2005 (Enterobacteriales: Enterobacteriaceae) was found in all analyzed individuals of one population of Aphis gossypii (Glover 1877) (Aphidinae: Aphidini) and sporadically in two populations of Hyalopterus. The lysogenic bacteriophage APSE-1 (A. pisum secondary endosymbiont, type 1) was detected in the three populations infected with H. defensa. This bacteriophage has been associated with moderate protection against braconid parasitoids in pea aphids. The high prevalence of Arsenophonus sp. in our samples is in accordance with previous studies indicating that, among gammaproteobacteria, this genus is one of the most widespread insect facultative symbionts.

Comparative genetics of invasive populations of walnut aphid, Chromaphis juglandicola, and its introduced parasitoid, Trioxys pallidus, in California

Coevolution may be an important component of the sustainability of importation biological control, but how frequently introduced natural enemies coevolve with their target pests is unclear. Here we explore whether comparative population genetics of the invasive walnut aphid, Chromaphis juglandicola, and its introduced parasitoid, Trioxys pallidus, provide insights into the localized breakdown of biological control services in walnut orchards in California. We found that sampled populations of C. juglandicola exhibited higher estimates of genetic differentiation (FST) than co-occurring populations of T. pallidus. In contrast, estimates of both the inbreeding coefficient (GIS) and contemporary gene flow were higher for T. pallidus than for C. juglandicola. We also found evidence of reciprocal outlier loci in some locations, but none showed significant signatures of selection. Synthesis and applications. Understanding the importance of coevolutionary interactions for the sustainability of biological control remains an important and understudied component of biological control research. Given the observed differences in gene flow and genetic differentiation among populations of T. pallidus and C. juglandicola, we suspect that temporary local disruption of biological control services may occur more frequently than expected while remaining stable at broader regional scales. Further research that combines genomewide single nucleotide polymorphism genotyping with measurements of phenotypic traits is needed to provide more conclusive evidence of whether the occurrence of outlier loci that display significant signatures of selection can be interpreted as evidence of the presence of a geographic mosaic of coevolution in this system.

Partitioning of herbivore hosts across time and food plants promotes diversification in the Megastigmus dorsalis oak gall parasitoid complex

Communities of insect herbivores and their natural enemies are rich and ecologically crucial components of terrestrial biodiversity. Understanding the processes that promote their origin and maintenance is thus of considerable interest. One major proposed mechanism is ecological speciation through host-associated differentiation (HAD), the divergence of a polyphagous species first into ecological host races and eventually into more specialized daughter species. The rich chalcid parasitoid communities attacking cynipid oak gall wasp hosts are structured by multiple host traits, including food plant taxon, host gall phenology, and gall structure. Here, we ask whether the same traits structure genetic diversity within supposedly generalist parasitoid morphospecies. We use mitochondrial DNA sequences and microsatellite genotypes to quantify HAD for Megastigmus (Bootanomyia) dorsalis, a complex of two apparently generalist cryptic parasitoid species attacking oak galls. Ancient Balkan refugial populations showed phenological separation between the cryptic species, one primarily attacking spring galls, and the other mainly attacking autumn galls. The spring species also contained host races specializing on galls developing on different host-plant lineages (sections Cerris vs. Quercus) within the oak genus Quercus. These results indicate more significant host-associated structuring within oak gall parasitoid communities than previously thought and support ecological theory predicting the evolution of specialist lineages within generalist parasitoids. In contrast, UK populations of the autumn cryptic species associated with both native and recently invading oak gall wasps showed no evidence of population differentiation, implying rapid recruitment of native parasitoid populations onto invading hosts, and hence potential for natural biological control. This is of significance given recent rapid range expansion of the economically damaging chestnut gall wasp, Dryocosmus kuriphilus, in Europe.

Hoverfly preference for high honeydew amounts creates enemy-free space for aphids colonizing novel host plants

The existence of an enemy-free space can play an important role in aphid host race formation processes, but little is known about the mechanisms that create an area of low predation pressure on particular host plants. In this paper, we identify a mechanism generating lower predation pressure that promotes the maintenance of the different host races of the pea aphid (Acyrthosiphon pisum) complex, a well-studied model for ecological speciation. The pea aphid consists of at least 15 genetically distinct host races which are native to specific host plants of the legume family, but can all develop on the universal host plant Vicia faba. Previous work showed that hoverfly (Episyrphus balteatus) oviposition preferences contribute to the enemy-free space that helps to maintain the different pea aphid host races, and that higher amounts of honeydew are more attractive to ovipositing hoverflies. Here we demonstrated that aphid honeydew is produced in large amounts when aphid reproduction rate was highest, and is an important oviposition cue for hoverflies under field conditions. However, on less suitable host plants, where honeydew production is reduced, pea aphids enjoy lower predation rates. A reduction in enemy pressure can mitigate the performance disadvantages of aphids colonizing a novel host and probably plays an important role in pea aphid host race formation.

Molecular and quantitative genetic differentiation in Sitobion avenae populations from both sides of the Qinling Mountains

Quantitative trait differences are often assumed to be correlated with molecular variation, but the relationship is not certain, and empirical evidence is still scarce. To address this issue, we sampled six populations of the cereal aphid Sitobion avenae from areas north and south of the Qinling Mountains, and characterized their molecular variation at seven microsatellite loci and quantitative variation at nine life-history traits. Our results demonstrated that southern populations had slightly longer developmental times of nymphs but much higher lifetime fecundity, compared to northern populations. Of the nine tested quantitative characters, eight differed significantly among populations within regions, as well as between northern and southern regions. Genetic differentiation in neutral markers was likely to have been caused by founder events and drift. Increased subdivision for quantitative characters was found in northern populations, but reduced in southern populations. This phenomenon was not found for molecular characters, suggesting the decoupling between molecular and quantitative variation. The pattern of relationships between FST and QST indicated divergent selection and suggested that local adaptation play a role in the differentiation of life-history traits in tested S. avenae populations, particularly in those traits closely related to reproduction. The main role of natural selection over genetic drift was also supported by strong structural differences in G-matrices among S. avenae populations. However, cluster analyses did not result in two groups corresponding to northern and southern regions. Genetic differentiation between northern and southern populations in neutral markers was low, indicating considerable gene flow between them. The relationship between molecular and quantitative variation, as well as its implications for differentiation and evolution of S. avenae populations, was discussed.

Armored scale insect endosymbiont diversity at the species level: genealogical patterns of Uzinura diasipipdicola in the Chionaspis pinifoliae-Chionaspis heterophyllae species complex (Hemiptera: Coccoidea: Diaspididae)

Armored scale insects and their primary bacterial endosymbionts show nearly identical patterns of co-diversification when viewed at the family level, though the persistence of these patterns at the species level has not been explored in this group. Therefore we investigated genealogical patterns of co-diversification near the species level between the primary endosymbiont Uzinura diaspidicola and its hosts in the Chionaspis pinifoliae-Chionaspis heterophyllae species complex. To do this we generated DNA sequence data from three endosymbiont loci (rspB, GroEL, and 16S) and analyzed each locus independently using statistical parsimony network analyses and as a concatenated dataset using Bayesian phylogenetic reconstructions. We found that for two endosymbiont loci, 16S and GroEL, sequences from U. diaspidicola were broadly associated with host species designations, while for rspB this pattern was less clear as C. heterophyllae (species S1) shared haplotypes with several other Chionaspis species. We then compared the topological congruence of the phylogenetic reconstructions generated from a concatenated dataset of endosymbiont loci (including all three loci, above) to that from a concatenated dataset of armored scale hosts, using published data from two nuclear loci (28S and EF1α) and one mitochondrial locus (COI-COII) from the armored scale hosts. We calculated whether the two topologies were congruent using the Shimodaira-Hasegawa test. We found no significant differences (P = 0.4892) between the topologies suggesting that, at least at this level of resolution, co-diversification of U. diaspidicola with its armored scale hosts also occurs near the species level. This is the first such study of co-speciation at the species level between U. diaspidicola and a group of armored scale insects.

Is ecological speciation a major trend in aphids? Insights from a molecular phylogeny of the conifer-feeding genus Cinara

Introduction

In the past decade ecological speciation has been recognized as having an important role in the diversification of plant-feeding insects. Aphids are host-specialised phytophagous insects that mate on their host plants and, as such, they are prone to experience reproductive isolation linked with host plant association that could ultimately lead to species formation. The generality of such a scenario remains to be tested through macroevolutionary studies. To explore the prevalence of host-driven speciation in the diversification of the aphid genus Cinara and to investigate alternative modes of speciation, we reconstructed a phylogeny of this genus based on mitochondrial, nuclear and Buchnera aphidicola DNA sequence fragments and applied a DNA-based method of species delimitation. Using a recent software (PhyloType), we explored evolutionary transitions in host-plant genera, feeding sites and geographic distributions in the diversification of Cinara and investigated how transitions in these characters have accompanied speciation events.

Results

The diversification of Cinara has been constrained by host fidelity to conifer genera sometimes followed by sequential colonization onto different host species and by feeding-site specialisation. Nevertheless, our analyses suggest that, at the most, only half of the speciation events were accompanied by ecological niche shifts. The contribution of geographical isolation in the speciation process is clearly apparent in the occurrence of species from two continents in the same clades in relatively terminal positions in our phylogeny. Furthermore, in agreement with predictions from scenarios in which geographic isolation accounts for speciation events, geographic overlap between species increased significantly with time elapsed since their separation.

Conclusions

The history of Cinara offers a different perspective on the mode of speciation of aphids than that provided by classic models such as the pea aphid. In this genus of aphids, the role of climate and landscape history has probably been as important as host-plant specialisation in having shaped present-day diversity.

Mitochondrial COI and morphological specificity of the mealy aphids (Hyalopterus ssp.) collected from different hosts in Europe (Hemiptera, Aphididae)

Forty three European population samples of mealy aphids from various winter and summer host plants were attributed to respective species of Hyalopterus by means of their partial sequences of mitochondrial COI gene. Used Hyalopterus samples emerged as monophyletic relative to outgroup and formed three major clades representing three host specific mealy aphid species in the Neighbor joining, Maximum parsimony, Maximum likelihood and Bayesian inference trees. Hyalopterus pruni and Hyalopterus persikonus emerged as a sister species, whilst Hyalopterus amygdali was located basally. Samples representing different clades in the molecular trees were used for canonical discrimination analysis based on twenty two morphological characters. Length of the median dorsal head hair enabled a 97.3 % separation of Hyalopterus amygdali from the remaining two species. No single character enabled satisfactory discrimination between apterous viviparous females of Hyalopterus pruni and Hyalopterus persikonus. A modified key for the morphological identification of Hyalopterus species is suggested and their taxonomic status discussed.

Regular or covert sex defines two lineages and worldwide superclones within the leaf-curl plum aphid (Brachycaudus helichrysi, Kaltenbach)

Asexual reproduction occurs widely in plants and animals, particularly in insects. Aphid species usually reproduce by cyclic parthenogenesis, but many species include obligate asexual lineages. We recently showed that the leaf-curl plum aphid, Brachycaudus helichrysi, actually encompasses two lineages, B. helichrysi H1 and H2. Ecological data suggest that these lineages have different life cycles. We conducted a large population genetics study, based on 14 microsatellite loci, to infer their respective life cycles and investigate their population structure and geographical distribution. Brachycaudus helichrysi H1 displayed the genetic signature of cyclical parthenogenesis, using plum trees as primary hosts for sexual reproduction, as classically described for B. helichrysi. This global survey showed that the Central Asian population of H1 was clearly differentiated from American-European populations. By contrast, B. helichrysi H2 displayed the typical signature of obligate asexual reproduction. H2 encompassed at least eight highly successful genotypes or superclones. This lack of ability to undergo sexual reproduction was confirmed for one of the superclones by sex induction experiments. We found only one B. helichrysi H2 population that underwent sexual reproduction, which was collected from peach trees, in Northern India. Our results confirm that H1 and H2 have different life cycles. Brachycaudus helichrysi H1 is clearly heteroecious using plum trees as primary hosts, while B. helichrysi H2 encompasses several anholocyclic lineages, and some heteroecious populations that until now have only been found associated with peach trees as primary hosts. We discuss implications of these findings for the pest status of B. helichrysi lineages.

Inconsistent genetic structure among members of a multitrophic system: did bruchid parasitoids (Horismenus spp.) escape the effects of bean domestication?

Anthropogenic range expansion and cultural practices have modified the distribution, abundance and genetic diversity of domesticated organisms, thereby altering multitrophic assemblages through space and time. The putative Mesoamerican domestication centre of the common bean, Phaseolus vulgaris L., in Mexico allows investigating the effects of plant domestication on the genetic structure of members of a multitrophic system. The aim of this study was to compare the evolutionary history of Horismenus parasitoids (Hymenoptera: Eulophidae) to those of their bruchid beetle hosts (Coleoptera: Bruchidae) and their domesticated host plant (P. vulgaris), in the context of traditional agriculture in Mexico. We analyzed the population genetic structure of four Horismenus species in Mexico using mitochondrial COI haplotype data. The two most abundant parasitoid species were Horismenus depressus and Horismenus missouriensis. Horismenus missouriensis were infected by Wolbachia endosymbionts and had little to no population differentiation (F(ST) = 0.06). We suspect the mitochondrial history of H. missouriensis to be blurred by Wolbachia, because differentiation among infected vs. non-infected individuals exists (F(ST) = 0.11). Populations of H. depressus were found to be highly differentiated (F(ST) = 0.34), but the genetic structuring could not be explained by tested spatial components. We then compared the genetic structure observed in this parasitoid species to previously published studies on bruchid beetles and their host plants. Despite extensive human-mediated migration and likely population homogenization of its two Acanthoscelides bruchid beetle hosts, H. depressus populations are structured like its host plant, by a recent dispersal from a diverse ancestral gene pool. Distinct evolutionary dynamics may explain inconsistent patterns among trophic levels. Parasitoids likely migrate from wild bean populations and are poorly adapted to bean storage conditions similar to their bruchid beetle hosts. Integrating several trophic levels to the study of evolutionary history has proven to be fruitful in detecting different ecological responses to human-mediated disturbances and host parasite interactions.

Speciation by symbiosis

In the Origin of Species, Darwin struggled with how continuous changes within a species lead to the emergence of discrete species. Molecular analyses have since identified nuclear genes and organelles that underpin speciation. In this review, we explore the microbiota as a third genetic component that spurs species formation. We first recall Ivan Wallin's original conception from the early 20th century on the role that bacteria play in speciation. We then describe three fundamental observations that justify a prominent role for microbes in eukaryotic speciation, consolidate exemplar studies of microbe-assisted speciation and incorporate the microbiota into classic models of speciation.

A molecular phylogeny of the Cephinae (Hymenoptera, Cephidae) based on mtDNA COI gene: a test of traditional classification

Cephinae is traditionally divided into three tribes and about 24 genera based on morphology and host utilization. There has been no study testing the monophyly of taxa under a strict phylogenetic criterion. A molecular phylogeny of Cephinae based on a total of 68 sequences of mtDNA COI gene, representing seven genera of Cephinae, is reconstructed to test the traditional limits and relationships of taxa. Monophyly of the traditional tribes is not supported. Monophyly of the genera are largely supported except for Pachycephus. A few host shift events are suggested based on phylogenetic relationships among taxa. These results indicate that a more robust phylogeny is required for a more plausible conclusion. We also report two species of Cephus for the first time from Turkey.

Morphometric relationship, phylogenetic correlation, and character evolution in the species-rich genus Aphis (Hemiptera: Aphididae)

Background

The species-rich genus Aphis consists of more than 500 species, many of them host-specific on a wide range of plants, yet very similar in general appearance due to convergence toward particular morphological types. Most species have been historically clustered into four main phenotypic groups (gossypii, craccivora, fabae, and spiraecola groups). To confirm the morphological hypotheses between these groups and to examine the characteristics that determine them, multivariate morphometric analyses were performed using 28 characters measured/counted from 40 species. To infer whether the morphological relationships are correlated with the genetic relationships, we compared the morphometric dataset with a phylogeny reconstructed from the combined dataset of three mtDNA and one nuclear DNA regions.

Principal Findings

Based on a comparison of morphological and molecular datasets, we confirmed morphological reduction or regression in the gossypii group unlike in related groups. Most morphological characteristics of the gossypii group were less variable than for the other groups. Due to these, the gossypii group could be morphologically well separated from the craccivora, fabae, and spiraecola groups. In addition, the correlation of the rates of evolution between morphological and DNA datasets was highly significant in their diversification.

Conclusions

The morphological separation between the gossypii group and the other species-groups are congruent with their phylogenetic relationships. Analysis of trait evolution revealed that the morphological traits found to be significant based on the morphometric analyses were confidently correlated with the phylogeny. The dominant patterns of trait evolution resulting in increased rates of short branches and temporally later evolution are likely suitable for the modality of Aphis speciation because they have adapted species-specifically, rapidly, and more recently on many different host plants.

A phylogenetic analysis of armored scale insects (Hemiptera: Diaspididae), based upon nuclear, mitochondrial, and endosymbiont gene sequences

Armored scale insects (Hemiptera: Diaspididae) are among the most invasive insects in the world. They have unusual genetic systems, including diverse types of paternal genome elimination (PGE) and parthenogenesis. Intimate relationships with their host plants and bacterial endosymbionts make them potentially important subjects for the study of co-evolution. Here, we expand upon recent phylogenetic work (Morse and Normark, 2006) by analyzing armored scale and endosymbiont DNA sequences from 125 species of armored scale insect, represented by 253 samples and eight outgroup species. We used fragments of four different gene regions: the nuclear protein-coding gene Elongation Factor 1α (EF1α), the large ribosomal subunit (28S) rDNA, a mitochondrial region spanning parts of cytochrome oxidase I (COI) and cytochrome oxidase II (COII), and the small ribosomal subunit (16S) rDNA from the primary bacterial endosymbiont Uzinura diaspidicola. Maximum likelihood, and Bayesian analyses were performed producing highly congruent topological results. A comparison of two datasets, one with and one without missing data, found that missing data had little effect on topology. Our results broadly corroborate several major features of the existing classification, although we do not find any of the subfamilies, tribes or subtribes to be monophyletic as currently constituted. Using ancestral state reconstruction we estimate that the ancestral armored scale had the late PGE sex system, and it may as well have been pupillarial, though results differed between reconstruction methods. These results highlight the need for a complete revision of this family, and provide the groundwork for future taxonomic work in armored scale insects.

Post-Pleistocene radiation of the pea aphid complex revealed by rapidly evolving endosymbionts

Adaptation to different resources has the potential to cause rapid species diversification, but few studies have been able to quantify the time scale of recent adaptive radiations. The pea aphid, Acyrthosiphon pisum, a model of speciation for host-specialized parasites, consists of several biotypes (races or species) living on distinct legume hosts. To document this radiation, we used rapidly evolving sequences from Buchnera, the maternally transmitted bacterial endosymbiont of aphids. Analyses of Buchnera pseudogene sequences revealed that 11 host-associated biotypes sort mostly into distinct matrilines despite low sequence divergence. A calibration based on divergence times of 7 sequenced genomes of Buchnera allowed us to date the last maternal ancestor of these biotypes between 8,000 and 16,000 years, with a burst of diversification at an estimated 3,600-9,500 years. The recency of this diversification, which is supported by microsatellite data, implies that the pea aphid complex ranks among the most rapid adaptive radiations yet documented. This diversification coincides with post-Pleistocene warming and with the domestication and anthropogenic range expansion of several of the legume hosts of pea aphids. Thus, we hypothesize that the new availability or abundance of resources triggered a cascade of divergence events in this newly formed complex.

Ecological niche models and coalescent analysis of gene flow support recent allopatric isolation of parasitoid wasp populations in the Mediterranean

Background

The integration of multiple complementary approaches is a powerful way to understand the processes of diversification and speciation. The parasitoid wasp Aphidius transcaspicus Telenga (Hymenoptera: Braconidae) is a parasitoid of Hyalopterus aphids across a wide geographic range. This species shows a remarkable degree of genetic structure among western, central, and eastern Mediterranean population clusters. In this paper we attempt to better characterize this genetic structure.

Methodology/Principal Findings

We use a Bayesian coalescent analysis of gene flow under the Isolation with Migration model using mitochondrial and microsatellite markers together with climate-based ecological niche models to better understand the genetic structure of A. transcaspicus in the Mediterranean. The coalescent analysis revealed low levels of migration among western and eastern Mediterranean populations (Nm<1) that were not statistically distinguishable from zero. Niche models showed that localities within population clusters each occupy areas of continuously high environmental suitability, but are separated from each other by large regions of completely unsuitable habitat that could limit dispersal. Overall, environmental characteristics were similar among the population clusters, though significant differences did emerge.

Conclusions/Significance

These results support contemporary allopatric isolation of Mediterranean populations of A. transcaspicus, which together with previous analyses indicating partial behaviorally mediated reproductive isolation, suggest that the early stages of cryptic speciation may be in progress.

Fine-scale cospeciation between Brachycaudus and Buchnera aphidicola: bacterial genome helps define species and evolutionary relationships in aphids

Aphids harbour an obligatory symbiont, Buchnera aphidicola, providing essential amino acids not supplied by their diet. These bacteria are transmitted vertically and phylogenic analyses suggest that they have 'cospeciated' with their hosts. We investigated this cospeciation phenomenon at a fine taxonomic level, within the aphid genus Brachycaudus. We used DNA-based methods of species delimitation in both organisms, to avoid biases in the definition of aphid and Buchnera species and to infer association patterns without the presumption of a specific interaction. Our results call into question certain 'taxonomic' species of Brachycaudus and suggest that B. aphidicola has diversified into independently evolving entities, each specific to a 'phylogenetic' Brachycaudus species. We also found that Buchnera and their hosts simultaneously diversified, in parallel. These results validate the use of Buchnera DNA data for inferring the evolutionary history of their host. The Buchnera genome evolves rapidly, making it the perfect tool for resolving ambiguities in aphid taxonomy. This study also highlights the usefulness of species delimitation methods in cospeciation studies involving species difficult to conceptualize--as is the case for bacteria--and in cases in which the taxonomy of the interacting organisms has not been determined independently and species definition depends on host association.

Molecular markers reveal strong geographic, but not host associated, genetic differentiation in Aphidius transcaspicus, a parasitoid of the aphid genus Hyalopterus

Host plant associated genetic differentiation is a common phenomenon in phytophagous insects, but the degree to which such associations sequentially drive diversification at higher trophic levels is not as well understood. A recent study examining neutral molecular markers in Hyalopterus aphids revealed that genetic structure in this genus is strongly determined by primary host plant use (Prunus spp.). In this paper, we take a similar approach to determine whether this host plant specificity has affected genetic structure in the parasitoid Aphidius transcaspicus, an important natural enemy of Hyalopterus spp. in the Mediterranean. Mitochondrial DNA (428 bp) and seven microsatellite loci were examined in parasitoids collected from aphid populations on almond, apricot, peach and plum trees from Spain and Greece. In contrast to the previous findings for Hyalopterus from the same regions, here we find no evidence for host associated diversification in A. transcaspicus at the species level or below, though geographic structure between regional populations is exceptionally high. These findings have several implications for our understanding of the ecology and evolution of A. transcaspicus as well as for its use as a biological control agent for Hyalopterus, suggesting that a consideration of host plant specificity may be less critical than factors such as climatic suitability or geographic origins of invasive populations.