Ecological divergence despite common mating sites: Genotypes and symbiotypes shed light on cryptic diversity in the black bean aphid species complex

Different host plants represent ecologically dissimilar environments for phytophagous insects. The resulting divergent selection can promote the evolution of specialized host races, provided that gene flow is reduced between populations feeding on different plants. In black bean aphids belonging to the Aphis fabae complex, several morphologically cryptic taxa have been described based on their distinct host plant preferences. However, host choice and mate choice are largely decoupled in these insects: they are host-alternating and migrate between specific summer host plants and shared winter hosts, with mating occurring on the shared hosts. This provides a yearly opportunity for gene flow among aphids using different summer hosts, and raises the question if and to what extent the ecologically defined taxa are reproductively isolated. Here, we analyzed a geographically and temporally structured dataset of microsatellite genotypes from A. fabae that were mostly collected from their main winter host Euonymus europaeus, and additionally from another winter host and fourteen summer hosts. The data reveals multiple, strongly differentiated genetic clusters, which differ in their association with different summer and winter hosts. The clusters also differ in the frequency of infection with two heritable, facultative endosymbionts, separately hinting at reproductive isolation and divergent ecological selection. Furthermore, we found evidence for occasional hybridization among genetic clusters, with putative hybrids collected more frequently in spring than in autumn. This suggests that similar to host races in other phytophagous insects, both prezygotic and postzygotic barriers including selection against hybrids maintain genetic differentiation among A. fabae taxa, despite a common mating habitat.

Similar cost of Hamiltonella defensa in experimental and natural aphid-endosymbiont associations

Endosymbiont-conferred resistance to parasitoids is common in aphids, but comes at a cost to the host in the absence of parasitoids. In black bean aphids (Aphis fabae), costs in terms of reduced lifespan and lifetime reproduction were demonstrated by introducing 11 isolates of the protective symbiont Hamiltonella defensa into previously uninfected aphid clones. Transfection of H. defensa isolates into a common genetic background allows to compare the costs of different endosymbiont isolates unconfounded by host genetic variation, but has been suggested to overestimate the realized costs of the endosymbiont in natural populations, because transfection creates new and potentially maladapted host-symbiont combinations that would be eliminated by natural selection in the field. In this experiment, we show that removing H. defensa isolates from their natural host clones with antibiotics results in a fitness gain that is comparable to the fitness loss from their introduction into two new clones. This suggests that estimating cost by transfecting endosymbiont isolates into a shared host genotype does not lead to gross overestimates of their realized costs, at least not in the two recipient genotypes used here. By comparing our data with data reported in previous publications using the same lines, we show that symbiont-induced costs may fluctuate over time. Thus, costs estimated after extended culture in the laboratory may not always be representative of the costs at the time of collection in the field. Finally, we report the accidental observation that two isolates from a distinct haplotype of H. defensa could not be removed by cefotaxime treatment, while all isolates from two other haplotypes were readily eliminated, which is suggestive of variation in susceptibility to this antibiotic in H. defensa.

Strong genotype-by-genotype interactions between aphid-defensive symbionts and parasitoids persist across different biotic environments

The dynamics of coevolution between hosts and parasites are influenced by their genetic interactions. Highly specific interactions, where the outcome of an infection depends on the precise combination of host and parasite genotypes (G × G interactions), have the potential to maintain genetic variation by inducing negative frequency‐dependent selection. The importance of this effect also rests on whether such interactions are consistent across different environments or modified by environmental variation (G × G × E interaction). In the black bean aphid, Aphis fabae, resistance to its parasitoid Lysiphlebus fabarum is largely determined by the possession of a heritable bacterial endosymbiont, Hamiltonella defensa, with strong G × G interactions between H. defensa and L. fabarum. A key environmental factor in this system is the host plant on which the aphid feeds. Here, we exposed genetically identical aphids harbouring three different strains of H. defensa to three asexual genotypes of L. fabarum and measured parasitism success on three common host plants of A. fabae, namely Vicia faba, Chenopodium album and Beta vulgaris. As expected, we observed the pervasive G × G interaction between H. defensa and L. fabarum, but despite strong main effects of the host plants on average rates of parasitism, this interaction was not altered significantly by the host plant environment (no G × G × E interaction). The symbiont‐conferred specificity of resistance is thus likely to mediate the coevolution of A. fabae and L. fabarum, even when played out across diverse host plants of the aphid.

Genetic origin and dispersal of the invasive soybean aphid inferred from population genetic analysis and approximate Bayesian computation

Biological invasion is considered among the most important global environmental problems. Knowledge of the source and dispersal routes of invasion could facilitate the eradication and control of invasive species. Soybean aphid, Aphis glycines, is among the most destructive soybean pests. For effective management of this pest, we conducted genetic analyses and approximate Bayesian computation analysis to determine the origins and dispersal of the aphid species, as well as the source of its invasion in the USA, using 8 microsatellite loci and the mitochondrial cytochrome c oxidase subunit I gene. We were able to identify a significant isolation by distance pattern and 3 genetic lineages in the microsatellite data but not in the mtDNA dataset. The genetic structure showed that the USA population had the closest relationship with those from Korea and Japan, indicating that the 2 latter populations might be the sources of the invasion to the USA. Both population genetic analyses and approximate Bayesian computation showed that the northeastern populations in China were the possible sources of the further spread of A. glycines to Indonesia. The dispersal history of this aphid can provide useful information for pest management strategies and can further help predict areas at risk of invasion.

The effect of host plants on genotype variability in fitness and honeydew composition of Aphis fabae

Aphid species can be polyphagous, feeding on multiple host plants across genera. As host plant species can have large variation in their phloem composition, this can affect aphid fitness and honeydew composition. Previous research showed significant intraspecific genotype variation in the composition of the honeydew carbohydrates of the black bean aphid Aphis fabae, with the ant attractant trisaccharide melezitose showing especially large variation across different genotypes. In this study, we test if variation in melezitose and carbohydrate composition of aphid honeydew could be linked to the adaptation of specific aphid genotypes to particular host plants. To this end, 4 high and 5 low melezitose secreting genotypes of the black bean aphid Aphis fabae were reared on 4 common host plants: broad bean, goosefoot, beet, and poppy. The carbohydrate composition, and in particular melezitose secretion, showed important aphid genotype and host plant interactions, with some genotypes being high melezitose secreting on 1 host plant but not on another. However, the interaction effects were not paralleled in the fitness measurements, even though there were significant differences in the average fitness across the different host plants. On the whole, this study demonstrates that aphid honeydew composition is influenced by complex herbivore-plant interactions. We discuss the relevance of these findings in the context of ant-aphid mutualisms and adaptive specialization in aphids.

Rapid evolution of symbiont-mediated resistance compromises biological control of aphids by parasitoids

There is growing interest in biological control as a sustainable and environmentally friendly way to control pest insects. Aphids are among the most detrimental agricultural pests worldwide, and parasitoid wasps are frequently employed for their control. The use of asexual parasitoids may improve the effectiveness of biological control because only females kill hosts and because asexual populations have a higher growth rate than sexuals. However, asexuals may have a reduced capacity to track evolutionary change in their host populations. We used a factorial experiment to compare the ability of sexual and asexual populations of the parasitoid Lysiphlebus fabarum to control caged populations of black bean aphids (Aphis fabae) of high and low clonal diversity. The aphids came from a natural population, and one‐third of the aphid clones harbored Hamiltonella defensa, a heritable bacterial endosymbiont that increases resistance to parasitoids. We followed aphid and parasitoid population dynamics for 3 months but found no evidence that the reproductive mode of parasitoids affected their effectiveness as biocontrol agents, independent of host clonal diversity. Parasitoids failed to control aphids in most cases, because their introduction resulted in strong selection for clones protected by H. defensa. The increasingly resistant aphid populations escaped control by parasitoids, and we even observed parasitoid extinctions in many cages. The rapid evolution of symbiont‐conferred resistance in turn imposed selection on parasitoids. In cages where asexual parasitoids persisted until the end of the experiment, they became dominated by a single genotype able to overcome the protection provided by H. defensa. Thus, there was evidence for parasitoid counteradaptation, but it was generally too slow for parasitoids to regain control over aphid populations. It appears that when pest aphids possess defensive symbionts, the presence of parasitoid genotypes able to overcome symbiont‐conferred resistance is more important for biocontrol success than their reproductive mode.

Aphid specialization on different summer hosts is associated with strong genetic differentiation and unequal symbiont communities despite a common mating habitat

Specialization on different host plants can promote evolutionary diversification of herbivorous insects. Work on pea aphids (Acyrthosiphon pisum) has contributed significantly to the understanding of this process, demonstrating that populations associated with different host plants exhibit performance trade-offs across hosts, show adaptive host choice and genetic differentiation and possess different communities of bacterial endosymbionts. Populations specialized on different secondary host plants during the parthenogenetic summer generations are also described for the black bean aphid (Aphis fabae complex) and are usually treated as different (morphologically cryptic) subspecies. In contrast to pea aphids, however, host choice and mate choice are decoupled in black bean aphids, because populations from different summer hosts return to the same primary host plant to mate and lay overwintering eggs. This could counteract evolutionary divergence, and it is currently unknown to what extent black bean aphids using different summer hosts are indeed differentiated. We addressed this question by microsatellite genotyping and endosymbiont screening of black bean aphids collected in summer from the goosefoot Chenopodium album (subspecies A. f. fabae) and from thistles of the genus Cirsium (subspecies A. f. cirsiiacanthoides) across numerous sites in Switzerland and France. Our results show clearly that aphids from Cirsium and Chenopodium exhibit strong and geographically consistent genetic differentiation and that they differ in their frequencies of infection with particular endosymbionts. The dependence on a joint winter host has thus not prevented the evolutionary divergence into summer host-adapted populations that appear to have evolved mechanisms of reproductive isolation within a common mating habitat.

Strong specificity in the interaction between parasitoids and symbiont-protected hosts

Coevolution between hosts and parasites may promote the maintenance of genetic variation in both antagonists by negative frequency-dependence if the host-parasite interaction is genotype-specific. Here we tested for specificity in the interaction between parasitoids (Lysiphlebus fabarum) and aphid hosts (Aphis fabae) that are protected by a heritable defensive endosymbiont, the γ-proteobacterium Hamiltonella defensa. Previous studies reported a lack of genotype specificity between unprotected aphids and parasitoids, but suggested that symbiont-conferred resistance might exhibit a higher degree of specificity. Indeed, in addition to ample variation in host resistance as well as parasitoid infectivity, we found a strong aphid clone-by-parasitoid line interaction on the rates of successful parasitism. This genotype specificity appears to be mediated by H. defensa, highlighting the important role that endosymbionts can play in host-parasite coevolution.

Maintaining genetic diversity and population panmixia through dispersal and not gene flow in a holocyclic heteroecious aphid species

Heteroecious holocyclic aphids exhibit both sexual and asexual reproduction and alternate among primary and secondary hosts. Most of these aphids can feed on several related hosts, and invasions to new habitats may limit the number of suitable hosts. For example, the aphid specialist Aphis glycines survives only on the primary host buckthorn (Rhamnus spp.) and the secondary host soybean (Glycine max) in North America where it is invasive. Owing to this specialization and sparse primary host distribution, host colonization events could be localized and involve founder effects, impacting genetic diversity, population structure and adaptation. We characterized changes in the genetic diversity and structure across time among A. glycines populations. Populations were sampled from secondary hosts twice in the same geographical location: once after secondary colonization (early season), and again immediately before primary host colonization (late season). We tested for evidence of founder effects and genetic isolation in early season populations, and whether or not late-season dispersal restored genetic diversity and reduced fragmentation. A total of 24 single-nucleotide polymorphisms and 6 microsatellites were used for population genetic statistics. We found significantly lower levels of genotypic diversity and more genetic isolation among early season collections, indicating secondary host colonization occurred locally and involved founder effects. Pairwise F(ST) decreased from 0.046 to 0.017 in early and late collections, respectively, and while genetic relatedness significantly decreased with geographical distance in early season collections, no spatial structure was observed in late-season collections. Thus, late-season dispersal counteracts the secondary host colonization through homogenization and increases genetic diversity before primary host colonization.

Climate effects on life cycle variation and population genetic architecture of the black bean aphid, Aphis fabae

Aphid species may exhibit different reproductive modes ranging from cyclical to obligate parthenogenesis. The distribution of life cycle variation in aphids is generally determined by ecological forces, mainly climate, because only sexually produced diapausing eggs can survive harsh winters or periods of absence of suitable host plants. Aphids are thus interesting models to investigate intrinsic and environmental factors shaping the competition among sexual and asexual lineages. We conducted a Europe-wide sampling of black bean aphids, Aphis fabae, and combined population genetic analyses based on microsatellite data with an experimental determination of life cycle strategies. Aphids were collected from broad beans (Vicia faba) as well as some Chenopodiaceae, but we detected no genetic differentiation between aphids from different host plants. Consistent with model predictions, life cycle variation was related to climate, with aphids from areas with cold winters investing more in sexual reproduction than aphids from areas with mild winters. Accordingly, only populations from mild areas exhibited a clear genetic signature of clonal reproduction. These differences arise despite substantial gene flow over large distances, which was evident from a very low geographic population structure and a lack of isolation-by-distance among 18 sites across distances of more than 1000 km. There was virtually no genetic differentiation between aphids with different reproductive modes, suggesting that new asexual lineages are formed continuously. Indeed, a surprising number of A. fabae genotypes even from colder climates produced some parthenogenetic offspring under simulated winter conditions. From this we predict that a shift to predominantly asexual reproduction could take place rapidly under climate warming.

Levels of clonal mixing in the black bean aphid Aphis fabae, a facultative ant mutualist

Aphids are a worldwide pest and an important model in ecology and evolution. Little is known, however, about the genetic structure of their colonies at a microgeographic level. For example, it remains largely unknown whether most species form monoclonal or polyclonal colonies. Here, we present the first detailed study on levels of clonal mixing in a nonsocial facultative ant mutualist, the black bean aphid Aphis fabae. In contrast to the earlier suggestion that colonies of this species are generally monoclonal, we found that across two subspecies of the black bean aphid, A. fabae cirsiiacanthoidis and A. fabae fabae, 32% and 67% of the aphid colonies were in fact polyclonal, consisting of a mix of up to four different clones, which resulted in an overall average relatedness within colonies of 0.90 and 0.79 in the two subspecies. Data further show that the average relatedness in A. f. cirsiiacanthoidis remained relatively constant throughout the season, which means that clonal erosion due to clonal selection more or less balanced with the influx of new clones from elsewhere. Nevertheless, relatedness tended to decrease over the lifetime of a given colony, implying that clonal mixing primarily resulted from the joining of pre-existing colonies as opposed to via simultaneous host colonisation by several foundresses. Widespread clonal mixing is argued to affect the ecology and evolution of the aphids in various important ways, for example with respect to the costs and benefits of group living, the evolution of dispersal and the interaction with predators as well as with the ant mutualists.

Genetic diversity and differentiation among laboratory and field populations of the soybean aphid, Aphis glycines

The soybean aphid, Aphis glycines Matsumura, is a recent invasive pest of soybean in North America. Currently, much research is focused on developing and characterizing soybean cultivars expressing host-plant resistance. During the initial phases of host-plant resistance screening, many of these studies use soybean aphid laboratory populations. Previous studies in other systems have documented substantial differences among laboratory and field populations. Whether or not this pattern exists in A. glycines is unknown, but it is extremely important when estimating the level of selection and virulence to host-plant resistant soybeans. In this study, we used seven microsatellite markers to estimate and compare genetic diversity and differentiation among five laboratory and 12 field populations. Our results indicate that soybean aphid laboratory populations are severely lacking in genotypic diversity and show extreme genetic differentiation among each other and to field populations. Continued use of laboratory populations for initial soybean aphid resistance screening could lead to erroneous estimations of the potential success for host-plant resistance.

Ample genetic variation but no evidence for genotype specificity in an all-parthenogenetic host-parasitoid interaction

Antagonistic coevolution between hosts and parasites can result in negative frequency-dependent selection and may thus be an important mechanism maintaining genetic variation in populations. Negative frequency-dependence emerges readily if interactions between hosts and parasites are genotype-specific such that no host genotype is most resistant to all parasite genotypes, and no parasite genotype is most infective on all hosts. Although there is increasing evidence for genotype specificity in interactions between hosts and pathogens or microparasites, the picture is less clear for insect host-parasitoid interactions. Here, we addressed this question in the black bean aphid (Aphis fabae) and its most important parasitoid Lysiphlebus fabarum. Because both antagonists are capable of parthenogenetic reproduction, this system allows for powerful tests of genotype x genotype interactions. Our test consisted of exposing multiple host clones to different parthenogenetic lines of parasitoids in all combinations, and this experiment was repeated with animals from four different sites. All aphids were free of endosymbiotic bacteria known to increase resistance to parasitoids. We observed ample genetic variation for host resistance and parasitoid infectivity, but there was no significant host clone x parasitoid line interaction, and this result was consistent across the four sites. Thus, there is no evidence for genotype specificity in the interaction between A. fabae and L. fabarum, suggesting that the observed variation is based on rather general mechanisms of defence and attack.

Population genetic structure of Aphis glycines

The soybean aphid (Aphis glycines Matsumura) is an invasive pest of cultivated soybean (Glycine max L.) in North America. After the initial invasion in 2000, the aphid has quickly spread across most of the United States and Canada, suggesting large-scale dispersal and rapid adaptation to new environments. Using microsatellite markers from closely related species, we compared the genetic diversity and the amount of genetic differentiation within and among 2 South Korean and 10 North American populations. Overall allelic polymorphism was low, never exceeding four alleles per locus. However, differences in genetic diversity were seen among South Korean and North American populations in terms of heterozygote excesses and genotypic richness. Within North America, two populations (Michigan and Ontario), had lower genetic diversities and exhibited high genetic differentiation compared with the remaining eight populations. The earlier collection time of Michigan and Ontario samples explained the genetic differences better than geographic subdivisions. These data indicate a pattern of small colonizing populations on soybeans, followed by rapid clonal amplification and subsequent large-scale dispersal across North America.

Genotypic variation and the role of defensive endosymbionts in an all-parthenogenetic host-parasitoid interaction

Models of host-parasite coevolution predict pronounced genetic dynamics if resistance and infectivity are genotype-specific or associated with costs, and if selection is fueled by sufficient genetic variation. We addressed these assumptions in the black bean aphid, Aphis fabae, and its parasitoid Lysiphlebus fabarum. Parasitoid genotypes differed in infectivity and host clones exhibited huge variation for susceptibility. This variation occurred at two levels. Clones harboring Hamiltonella defensa, a bacterial endosymbiont known to protect pea aphids against parasitoids, enjoyed greatly reduced susceptibility, yet clones without H. defensa also exhibited significant variation. Although there was no evidence for genotype-specificity in the H. defensa-free clones' interaction with parasitoids, we found such evidence in clones containing the bacterium. This suggests that parasitoid genotypes differ in their ability to overcome H. defensa, resulting in an apparent host x parasitoid genotype interaction that may in fact be due to an underlying symbiont x parasitoid genotype interaction. Aphid susceptibility to parasitoids correlated negatively with fecundity and rate of increase, due to H. defensa-bearing clones being more fecund on average. Hence, possessing symbionts may also be favorable in the absence of parasitoids, which raises the question why H. defensa does not go to fixation and highlights the need to develop new models to understand the dynamics of endosymbiont-mediated coevolution.

Microsatellite flanking region similarities among different loci within insect species

Although microsatellites are ubiquitous in eukaryota, the number of available markers varies strongly among taxa. This meta-analysis was conducted on 32 insect species. Sequences were obtained from two assembled whole genomes, whole genome shotgun (WGS) sequences from 10 species and screening partial genomic libraries for microsatellites from 23 species. We have demonstrated: (1) strong differences in the abundance of microsatellites among species; (2) that microsatellites within species are often grouped into families based on similarities in their flanking sequences; (3) that the proportion of microsatellites grouped into families varies strongly among taxa; and (4) that microsatellite families were significantly more often associated with transposable elements - or their remnants - than unique microsatellite sequences.

Phylogeny of the genus Aphis Linnaeus, 1758 (Homoptera: Aphididae) inferred from mitochondrial DNA sequences

Aphis is the largest aphid genus in the world and contains several of the most injurious aphid pests. It is also the most reluctant aphid genus to any comprehensive taxonomic treatment: while most species are easily classified into "species groups" that form well defined entities, numerous species within these groups are difficult to tell apart morphologically and identification keys remain ambiguous and mostly rely on host plant affiliation. In this paper, we used partial sequences of COI/COII and CytB genes to reconstruct the first phylogeny of Aphis and discuss the present systematics. The monophyly of the subgenus Bursaphis and of the tree major species groups, Black aphid, Black backed aphid and frangulae-like species was recovered by all phylogenetic analyses. However our data suggested that the nominal subgenus was not monophyletic. Relationships between major species groups were often ambiguous but "Black" and "Black backed" species groups appeared as sister clades. The most striking result of this study was that our molecular data met the same limits as the morphological characters used in classifications: mitochondrial DNA did not allow the differentiation of species that are difficult to identify. Further, interspecies relationships within groups of species for which taxonomic treatment is difficult stayed unresolved. This suggests that species delineation in the genus Aphis is often ambiguous and that diversification might have been a rapid process.