Genome scans reveal candidate regions involved in the adaptation to host plant in the pea aphid complex

A major goal in evolutionary biology is to uncover the genetic basis of adaptation. Divergent selection exerted on ecological traits may result in adaptive population differentiation and reproductive isolation and affect differentially the level of genetic divergence along the genome. Genome‐wide scan of large sets of individuals from multiple populations is a powerful approach to identify loci or genomic regions under ecologically divergent selection. Here, we focused on the pea aphid, a species complex of divergent host races, to explore the organization of the genomic divergence associated with host plant adaptation and ecological speciation. We analysed 390 microsatellite markers located at variable distances from predicted genes in replicate samples of sympatric populations of the pea aphid collected on alfalfa, red clover and pea, which correspond to three common host‐adapted races reported in this species complex. Using a method that accounts for the hierarchical structure of our data set, we found a set of 11 outlier loci that show higher genetic differentiation between host races than expected under the null hypothesis of neutral evolution. Two of the outliers are close to olfactory receptor genes and three other nearby genes encoding salivary proteins. The remaining outliers are located in regions with genes of unknown functions, or which functions are unlikely to be involved in interactions with the host plant. This study reveals genetic signatures of divergent selection across the genome and provides an inventory of candidate genes responsible for plant specialization in the pea aphid, thereby setting the stage for future functional studies.     

Identifying genomic hotspots of differentiation and candidate genes involved in the adaptive divergence of pea aphid host races

Identifying the genomic bases of adaptation to novel environments is a long‐term objective in evolutionary biology. Because genetic differentiation is expected to increase between locally adapted populations at the genes targeted by selection, scanning the genome for elevated levels of differentiation is a first step towards deciphering the genomic architecture underlying adaptive divergence. The pea aphid Acyrthosiphon pisum is a model of choice to address this question, as it forms a large complex of plant‐specialized races and cryptic species, resulting from recent adaptive radiation. Here, we characterized genomewide polymorphisms in three pea aphid races specialized on alfalfa, clover and pea crops, respectively, which we sequenced in pools (poolseq). Using a model‐based approach that explicitly accounts for selection, we identified 392 genomic hotspots of differentiation spanning 47.3 Mb and 2,484 genes (respectively, 9.12% of the genome size and 8.10% of its genes). Most of these highly differentiated regions were located on the autosomes, and overall differentiation was weaker on the X chromosome. Within these hotspots, high levels of absolute divergence between races suggest that these regions experienced less gene flow than the rest of the genome, most likely by contributing to reproductive isolation. Moreover, population‐specific analyses showed evidence of selection in every host race, depending on the hotspot considered. These hotspots were significantly enriched for candidate gene categories that control host–plant selection and use. These genes encode 48 salivary proteins, 14 gustatory receptors, 10 odorant receptors, five P450 cytochromes and one chemosensory protein, which represent promising candidates for the genetic basis of host–plant specialization and ecological isolation in the pea aphid complex. Altogether, our findings open new research directions towards functional studies, for validating the role of these genes on adaptive phenotypes.     

LARGE-SCALE CANDIDATE GENE SCAN REVEALS THE ROLE OF CHEMORECEPTOR GENES IN HOST PLANT SPECIALIZATION AND SPECIATION IN THE PEA APHID

Understanding the drivers of speciation is critical to interpreting patterns of biodiversity. The identification of the genetic changes underlying adaptation and reproductive isolation is necessary to link barriers to gene flow to the causal origins of divergence. Here, we present a novel approach to the genetics of speciation, which should complement the commonly used approaches of quantitative trait locus mapping and genome-wide scans for selection. We present a large-scale candidate gene approach by means of sequence capture, applied to identifying the genetic changes underlying reproductive isolation in the pea aphid, a model system for the study of ecological speciation. Targeted resequencing enabled us to scale up the candidate gene approach, specifically testing for the role of chemosensory gene families in host plant specialization. Screening for the signature of divergence under selection at 172 candidate and noncandidate loci, we revealed a handful of loci that show high levels of differentiation among host races, which almost all correspond to odorant and gustatory receptor genes. This study offers the first indication that some chemoreceptor genes, often tightly linked together in the genome, could play a key role in local adaptation and reproductive isolation in the pea aphid and potentially other phytophagous insects. Our approach opens a new route toward the functional genomics of ecological speciation.     

Ecological specialization correlates with genotypic differentiation in sympatric host-populations of the pea aphid

The pea aphid, Acyrthosiphon pisum, encompasses distinct host races specialized on various Fabaceae species, but the extent of genetic divergence associated with ecological specialization varies greatly depending on plant and geographic origins of aphid populations. Here, we studied the genetic structure of French sympatric pea aphid populations collected on perennial (pea and faba bean) and annual (alfalfa and red clover) hosts using 14 microsatellite loci. Classical and Bayesian population genetics analyses consistently identified genetic clusters mostly related to plant origin: the pea/faba bean cluster was highly divergent from the red clover and the alfalfa ones, indicating they represent different stages along the continuum of genetic differentiation. Some genotypes were assigned to a cluster differing from the one expected from their plant origin while others exhibited intermediate genetic characteristics. These results suggest incomplete barriers to gene flow. However, this limited gene flow seems insufficient to prevent ecological specialization and genetic differentiation in sympatry.     

REPRODUCTIVE ISOLATION BETWEEN SYMPATRIC RACES OF PEA APHIDS. I. GENE FLOW RESTRICTION AND HABITAT CHOICE

Determining the extent and causes of barriers to gene flow between genetically divergent populations or races of single species is an important complement to post facto analyses of the causes of reproductive isolation between recognized species. Sympatric populations of pea aphids (Acyrthosiphon pisum Harris, Homoptera: Aphididae) on alfalfa and red clover are highly genetically divergent and locally adapted. Here, hierarchical estimates of population structure based on F_st suggest that gene exchange between closely adjacent aphid populations on the two hosts is highly restricted relative to that among fields of the same host plant. Although these host-associated races are presently considered to be the same subspecies, they appear to be significantly reproductively isolated, suggesting incipient speciation. Habitat (host) choice was investigated as the first in a temporal series of factors that could reduce gene exchange between these sympatric populations. Field studies of winged colonists to newly planted fields of each host suggest pronounced habitat fidelity. This result was verified using replicated observations of the host choice behavior of different aphid genotypes for which the relative demographic performance on each host was known. These laboratory observations of behavior revealed a strong genetic correlation between habitat choice (or acceptance) and the relative performance in each habitat. Because mating occurs on the host plant, habitat choice in this system leads to assortative mating and is therefore a major cause of reproductive isolation between the sympatric pea aphid populations on alfalfa and clover. However, the extent of dispersal between hosts estimated from the field study of winged colonists (9–11%) is too great to be consistent with the genetic divergence estimated between the races. This suggests that barriers to gene flow other than host choice also exist, such as selection against migrants or hybrids in the parental environments, hybrid sterility, or hybrid breakdown.     

Patterns,causes and consequences of defensive microbiome dynamics across multiple scales

The microbiome can significantly impact host phenotypes and serve as an additional source of heritable genetic variation. While patterns across eukaryotes are consistent with a role for symbiotic microbes in host macroevolution, few studies have examined symbiont‐driven host evolution or the ecological implications of a dynamic microbiome across temporal, spatial or ecological scales. The pea aphid, Acyrthosiphon pisum, and its eight heritable bacterial endosymbionts have served as a model for studies on symbiosis and its potential contributions to host ecology and evolution. But we know little about the natural dynamics or ecological impacts of the heritable microbiome of this cosmopolitan insect pest. Here we report seasonal shifts in the frequencies of heritable defensive bacteria from natural pea aphid populations across two host races and geographic regions. Microbiome dynamics were consistent with symbiont responses to host‐level selection and findings from one population suggested symbiont‐driven adaptation to seasonally changing parasitoid pressures. Conversely, symbiont levels were negatively correlated with enemy‐driven mortality when measured across host races, suggesting important ecological impacts of host race microbiome divergence. Rapid drops in symbiont frequencies following seasonal peaks suggest microbiome instability in several populations, with potentially large costs of ‘superinfection’ under certain environmental conditions. In summary, the realization of several laboratory‐derived, a priori expectations suggests important natural impacts of defensive symbionts in host‐enemy eco‐evolutionary feedbacks. Yet negative findings and unanticipated correlations suggest complexities within this system may limit or obscure symbiont‐driven contemporary evolution, a finding of broad significance given the widespread nature of defensive microbes across plants and animals.     

Fusion-fission experiments in Aphidius: evolutionary split without isolation in response to environmental bimodality

Studying host-based divergence naturally maintained by a balance between selection and gene flow can provide valuable insights into genetic underpinnings of host adaptation and ecological speciation in parasites. Selection-gene flow balance is often postulated in sympatric host races, but direct experimental evidence is scarce. In this study, we present such evidence obtained in host races of Aphidius ervi, an important hymenopteran agent of biological control of aphids in agriculture, using a novel fusion-fission method of gene flow perturbation. In our study, between-race genetic divergence was obliterated by means of advanced hybridisation, followed by a multi-generation exposure of the resulting genetically uniform hybrid swarm to a two-host environment. This fusion-fission procedure was implemented under two contrasting regimes of between-host gene flow in two replicated experiments involving different racial pairs. Host-based genetic fission in response to environmental bimodality occurred in both experiments in as little as six generations of divergent adaptation despite continuous gene flow. We demonstrate that fission recovery of host-based divergence evolved faster and hybridisation-induced linkage disequilibrium decayed slower under restricted (6.7%) compared with unrestricted gene flow, directly pointing at a balance between gene flow and divergent selection. We also show, in four separate tests, that random drift had no or little role in the observed genetic split. Rates and patterns of fission divergence differed between racial pairs. Comparative linkage analysis of these differences is currently under way to test for the role of genomic architecture of adaptation in ecology-driven divergent evolution.     

Host-based divergence in populations of the pea aphid: insights from nuclear markers and the prevalence of facultative symbionts

In North America, the pea aphid Acyrthosiphon pisum encompasses ecologically and genetically distinct host races that offer an ideal biological system for studies on sympatric speciation. In addition to its obligate symbiont Buchnera, pea aphids harbour several facultative and phylogenetically distant symbionts. We explored the relationships between host races of A. pisum and their symbiotic microbiota to gain insights into the historical process of ecological specialization and symbiotic acquisition in this aphid. We used allozyme and microsatellite markers to analyse the extent of genetic differentiation between populations of A. pisum on pea, alfalfa and clover in France. In parallel, we examined: (i) the distribution of four facultative symbionts; and (ii) the genetic variation in the Buchnera genome across host-associated populations of A. pisum. Our study clearly demonstrates that populations of A. pisum on pea, clover and alfalfa in France are genetically divergent, which indicates that they constitute distinct host races. We also found a very strong association between host races of A. pisum and their symbiotic microbiota. We stress the need for phylogeographic studies to shed light on the process of host-race formation and acquisition of facultative symbionts in A. pisum. We also question the effects of these symbionts on aphid host fitness, including their role in adaptation to a host plant.     

Attachment forces of pea aphids (Acyrthosiphon pisum) on different legume species

1. Sympatric populations of insects adapted to different host plants are good model systems not only to study how they adapt to the chemistry of their food plant, but also to investigate whether morphological modifications evolved enabling them to live successfully on a certain plant species. 2. The pea aphid, Acyrthosiphon pisum (Harris) encompasses at least 11 genetically distinct sympatric host races, each showing a preference for a certain legume species. The leaflet surfaces of these legumes differ considerably in their wax coverage. 3. It was investigated whether the attachment structures of three pea aphid genotypes from different host races are adapted to the different surface properties of their host plants and whether they show differences in their attachment ability on the respective host and non‐host plants. 4. The surface morphology of plants and aphid tarsi was examined using SEM (scanning electron microscopy). The ability of the aphids to walk on specific surfaces was tested using traction force measurements. 5. The presence of wax blooms on the leaflets lowers the aphids' attachment ability considerably and diminishes their subsequent attachment on ‘neutral’ surfaces like glass. The pea aphid host races differ in their ability to walk on certain surfaces. However, the genotype from the adapted aphid host race was not necessarily the one with the best walking performance on their host plant. All aphids, regardless of the original host plant, were most efficient on the neutral control surface glass. The general host plant Vicia faba was the plant with the most favourable surface for all aphid host races.     

Recognition of host-specific chemical stimulants in two sympatric host races of the pea aphid Acyrthosiphon pisum

Abstract.  1. In ecological speciation , adaptation to variation in the external environment provides the crucial push that starts the process of genetic divergence and eventually leads to speciation. This emphasis on the role of ecological specialisation in speciation events has brought with it a renewed interest in its proximate mechanisms in recently diverged groups such as host races. Here, the proximate mechanisms of feeding specialisation are investigated in two host races of the pea aphid Acyrthosiphon pisum .
2. Using alfalfa and clover extracts, enclosed in diet chambers or applied on whole plants, it is shown that feeding specialisation depends on recognition of stimulants specific to the host plant, not on deterrents or toxins specific to the non-host plants.
3. Because pea aphids mate on their host plant, feeding specialisation leads to de facto assortative mating. This study suggests that behavioural recognition of host-specific chemicals, rather than avoidance of deterrents or/and plant toxins, contributes to gene flow restriction between the alfalfa and clover host races.     

Localizing FST outliers on a QTL map reveals evidence for large genomic regions of reduced gene exchange during speciation‐with‐gene‐flow

Populations that maintain phenotypic divergence in sympatry typically show a mosaic pattern of genomic divergence, requiring a corresponding mosaic of genomic isolation (reduced gene flow). However, mechanisms that could produce the genomic isolation required for divergence‐with‐gene‐flow have barely been explored, apart from the traditional localized effects of selection and reduced recombination near centromeres or inversions. By localizing F_ST outliers from a genome scan of wild pea aphid host races on a Quantitative Trait Locus (QTL) map of key traits, we test the hypothesis that between‐population recombination and gene exchange are reduced over large ‘divergence hitchhiking’ (DH) regions. As expected under divergence hitchhiking, our map confirms that QTL and divergent markers cluster together in multiple large genomic regions. Under divergence hitchhiking, the nonoutlier markers within these regions should show signs of reduced gene exchange relative to nonoutlier markers in genomic regions where ongoing gene flow is expected. We use this predicted difference among nonoutliers to perform a critical test of divergence hitchhiking. Results show that nonoutlier markers within clusters of F_ST outliers and QTL resolve the genetic population structure of the two host races nearly as well as the outliers themselves, while nonoutliers outside DH regions reveal no population structure, as expected if they experience more gene flow. These results provide clear evidence for divergence hitchhiking, a mechanism that may dramatically facilitate the process of speciation‐with‐gene‐flow. They also show the power of integrating genome scans with genetic analyses of the phenotypic traits involved in local adaptation and population divergence.     

Differential introgression causes genealogical discordance in host races of Acrocercops transecta (Insecta: Lepidoptera)

Recently diverged populations often exhibit incomplete reproductive isolation, with a low level of gene flow continuing between populations. Previous studies have shown that, even under a low level of gene flow, genetic divergence between populations can proceed at the loci governing local adaptation and reproductive isolation but not at other neutral loci. A leaf‐mining moth, Acrocercops transecta, consists of Juglans‐ and Lyonia‐associated host races. The two host races differ in host preferences of ovipositing females and in larval adaptation to host plants but mate readily in the laboratory, producing fertile hybrids. The Juglans and Lyonia races are often sympatric in the wild, implying that gene introgression could occur in nature between the two host races. We tested this hypothesis by combining phylogenetic analyses with coalescent simulations, focusing on mitochondrial genes (COI and ND5) and the nuclear Tpi, Per and Ldh genes located on the Z‐chromosome. The mitochondrial genes clearly distinguished the Lyonia race from the Juglnas race, whereas the Tpi, Per and Ldh genealogies did not reflect the two host races. Coalescent simulations indicated gene flow at the three Z‐linked genes in both directions, whereas there was no introgression in the mitochondrial genes. The lack of introgression in mitochondrial genes suggests that female host preference is the primary force leading to the bifurcation of maternally inherited loci. Thus, the results show that a low level of gene flow coupled with the inflexible female host preference differentiates histories of divergence between maternally and biparentally inherited genes in this host race system.     

WIDESPREAD HOST‐DEPENDENT HYBRID UNFITNESS IN THE PEA APHID SPECIES COMPLEX

Linking adaptive divergence to hybrid unfitness is necessary to understand the ecological factors contributing to reproductive isolation and speciation. To date, this link has been demonstrated in few model systems, most of which encompass ecotypes that occupy relatively early stages in the speciation process. Here we extend these studies by assessing how host‐plant adaptation conditions hybrid fitness in the pea aphid, Acyrthosiphon pisum. We made crosses between and within five pea aphid biotypes adapted to different host plants and representing various stages of divergence within the complex. Performance of F1 hybrids and nonhybrids was assessed on a “universal” host that is favorable to all pea aphid biotypes in laboratory conditions. Although hybrids performed equally well as nonhybrids on the universal host, their performance was much lower than nonhybrids on the natural hosts of their parental populations. Hence, hybrids, rather than being intrinsically deficient, are maladapted to their parents’ hosts. Interestingly, the impact of this maladaptation was stronger in certain hybrids from crosses involving the most divergent biotype, suggesting that host‐dependent postzygotic isolation has continued to evolve late in divergence. Even though host‐independent deficiencies are not excluded, hybrid maladaptation to parental hosts supports the hypothesis of ecological speciation in this complex.     

Hitching a lift on the road to speciation

Understanding how speciation can take place in the presence of homogenizing gene flow remains a major challenge in evolutionary biology. In the early stages of ecological speciation, reproductive isolation between populations occupying different habitats is expected to be concentrated around genes for local adaptation. These genomic regions will show high divergence while gene exchange in other regions of the genome should continue relatively unimpaired, resulting in low levels of differentiation. The problem is to explain how speciation progresses from this point towards complete reproductive isolation, allowing genome‐wide divergence. A new study by Via and West (2008) on speciation between host races of the pea aphid, Acyrthosiphon pisum, introduces the mechanism of ‘divergence hitchhiking’ which can generate large ‘islands of differentiation’ and facilitate the build‐up of linkage disequilibrium, favouring increased reproductive isolation. This idea potentially removes a major stumbling block to speciation under continuous gene flow.     

Host‐associated divergence in sympatric host races of the leaf beetle Lochmaea capreae: implications for local adaptation and reproductive isolation

Ecological specialization is widely recognized as a major determinant of the emergence and maintenance of biodiversity. We studied two critical facets of specialization – local adaptation and habitat choice – in the host races of the leaf beetle Lochmaea capreae on willow and birch. Our results revealed that there is asymmetric disruptive selection for host use traits, and host races achieved different adaptive sets of life history traits through association with their host plant. Beetles from each host race exhibited food and oviposition preference for their own host plant. Reciprocal transplant displayed significant variation in host acceptance and performance: all families from the willow race rejected the alternative host plant before initiation of feeding and all died on this host plant. By contrast, all families from the birch race accepted willow for feeding, but they consumed less and performed less well. Intriguingly, families that performed well on birch also performed well on willow, suggesting positive genetic correlation rather than genetic trade‐offs. Our results suggest that the major proximal determinant of host specialization in the willow race is the behavioural acceptance of a plant rather than the toxicity of the food resource. However, in the birch race a combination of behavioural host acceptance and performance may play a role in specialization. Our study sheds light on the mechanisms by which divergent host adaptation might influence the evolution of reproductive isolation between herbivorous populations.     

Genomic evidence for divergence with gene flow in host races of the larch budmoth

Ecological divergence in the face of gene flow has recently become implicated as a potentially important cause of speciation and adaptive radiation. Here, we develop a genomic approach to test for divergent selection in sympatric host races of the larch budmoth Zeiraphera diniana (Lepidoptera: Tortricidae). We analysed hundreds of amplified fragment length polymorphism markers in 92 individuals in sympatric and allopatric populations, and in two backcross broods used to map the markers to individual chromosomes. The results directly confirm the existence of natural hybridization and demonstrate strong heterogeneity between chromosomes in terms of molecular divergence between host races (the average level of divergence was FST = 0.216). However, genomic heterogeneity was not found when we analysed divergence between geographically separated populations of the same host race. We conclude that the variance of the level of sympatric divergence among chromosomes is the footprint of divergent selection acting on a few linkage groups, combined with appreciable gene flow that homogenizes between-race variation at the remaining linkage groups. These results, coupled with other recent multilocus analyses of sister species pairs, demonstrate that selection-driven sympatric phase of genetic divergence in the presence of gene flow is a likely feature of speciation.     

Population differentiation and genetic variation in performance on eight hosts in the pea aphid complex

Phytophagous insects frequently use multiple host-plant species leading to the evolution of specialized host-adapted populations and sometimes eventually to speciation. Some insects are confronted with a large number of host-plant species, which may provide complex routes of gene flow between host-adapted populations. The pea aphid (Acyrthosiphon pisum) attacks a broad range of plants in the Fabaceae and it is known that populations on Trifolium pratense and Medicago sativa can be highly specialized at exploiting these species. To find out whether adaptation to a broad range of co-occurring hosts has occurred, we tested the performance of pea aphid clones collected from eight host-plant genera on all of these plants in a reciprocal transfer experiment. We provide evidence for pervasive host-plant specialization. The high performance of all aphid clones on Vicia faba suggests that this host plant could be a site of gene flow between different populations that could limit further host-associated divergence. The genetic variance in host-plant usage was partitioned into within- and among-population components, which represent different levels of host adaptation. Little evidence of within-population trade-offs in performance on different plant species was found.     

Evolutionary history of aphid-plant associations and their role in aphid diversification

Aphids are intimately linked with their host plants that constitute their only food resource and habitat, and thus impose considerable selective pressure on their evolution. It is therefore commonly assumed that host plants have greatly influenced the diversification of aphids. Here, we review what is known about the role of host plant association on aphid speciation by examining both macroevolutionary and population-level studies. Phylogenetic studies conducted at different taxonomic levels show that, as in many phytophagous insect groups, the radiation of angiosperms has probably favoured the major Tertiary diversification of aphids. These studies also highlight many aphid lineages constrained to sets of related host plants, suggesting strong evolutionary commitment in aphids’ host plant choice, but they fail to document cospeciation events between aphid and host lineages. Instead, phylogenies of several aphid genera reveal that divergence events are often accompanied by host shifts, and suggest, without constituting a formal demonstration, that aphid speciation could be a consequence of adaptation to new hosts. Experimental and field studies below the species level support reproductive isolation between host races as partly due to divergent selection by their host plants. Selected traits are mainly feeding performances and life cycle adaptations to plant phenology. Combined with behavioural preference for favourable host species, these divergent adaptations can induce pre- and post-zygotic barriers between host-specialized aphid populations. However, the hypothesis of host-driven speciation is seldom tested formally and must be weighed against overlooked explanations involving geographic isolation and non-ecological reproductive barriers in the process of speciation.     

Non-parallel morphological divergence following colonization of a new host plant

Adaptation to new ecological niches is known to spur population diversification and may lead to speciation if gene flow is ceased. While adaptation to the same ecological niche is expected to be parallel, it is more difficult to predict whether selection against maladaptive hybridization in secondary sympatry results in parallel divergence also in traits that are not directly related to the ecological niches. Such parallelisms in response to selection for reproductive isolation can be identified through estimating parallelism in reproductive character displacement across different zones of secondary contact. Here, we use a host shift in the phytophagous peacock fly Tephritis conura, with both host races represented in two geographically separate areas East and West of the Baltic Sea to investigate convergence in morphological adaptations. We asked (i) if there are consistent morphological adaptations to a host plant shift and (ii) if the response to secondary sympatry with the alternate host race is parallel across contact zones. We found surprisingly low and variable, albeit significant, divergence between host races. Only one trait, the length of the female ovipositor, which serves an important function in the interaction with the hosts, was consistently different between host races. Instead, co-existence with the other host race significantly affected the degree of morphological divergence, but the divergence was largely driven by different traits in different contact zones. Thus, local stochastic fixation or reinforcement could generate trait divergence, and additional evidence is needed to conclude whether divergence is locally adaptive.