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.     

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.     

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.     

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.     

Divergent selection and heterogeneous genomic divergence

Levels of genetic differentiation between populations can be highly variable across the genome, with divergent selection contributing to such heterogeneous genomic divergence. For example, loci under divergent selection and those tightly physically linked to them may exhibit stronger differentiation than neutral regions with weak or no linkage to such loci. Divergent selection can also increase genome‐wide neutral differentiation by reducing gene flow (e.g. by causing ecological speciation), thus promoting divergence via the stochastic effects of genetic drift. These consequences of divergent selection are being reported in recently accumulating studies that identify: (i) ‘outlier loci’ with higher levels of divergence than expected under neutrality, and (ii) a positive association between the degree of adaptive phenotypic divergence and levels of molecular genetic differentiation across population pairs [‘isolation by adaptation’ (IBA)]. The latter pattern arises because as adaptive divergence increases, gene flow is reduced (thereby promoting drift) and genetic hitchhiking increased. Here, we review and integrate these previously disconnected concepts and literatures. We find that studies generally report 5–10% of loci to be outliers. These selected regions were often dispersed across the genome, commonly exhibited replicated divergence across different population pairs, and could sometimes be associated with specific ecological variables. IBA was not infrequently observed, even at neutral loci putatively unlinked to those under divergent selection. Overall, we conclude that divergent selection makes diverse contributions to heterogeneous genomic divergence. Nonetheless, the number, size, and distribution of genomic regions affected by selection varied substantially among studies, leading us to discuss the potential role of divergent selection in the growth of regions of differentiation (i.e. genomic islands of divergence), a topic in need of future investigation.     

Selection and genomic differentiation during ecological speciation: isolating the contributions of host association via a comparative genome scan of Neochlamisus bebbianae leaf beetles

This study uses a comparative genome scan to evaluate the contributions of host plant related divergent selection to genetic differentiation and ecological speciation in maple- and willow-associated populations of Neochlamisus bebbianae leaf beetles. For each of 15 pairwise population comparisons, we identified "outlier loci" whose strong differentiation putatively reflects divergent selection. Of 447 AFLP loci, 15% were outliers across multiple population comparisons, and low linkage disequilibrium indicated that these outliers derived from multiple regions of the genome. Outliers were further classified as "host-specific" if repeatedly observed in "different-host" population comparisons but never in "same-host" comparisons. Outliers exhibiting the opposite pattern were analogously classified as "host-independent." Host-specific outliers represented 5% of all loci and were more frequent than host-independent outliers, thus revealing a large role for host-adaptation in population genomic differentiation. Evidence that host-related selection can promote divergence despite gene flow was provided by population trees. These were structured by host-association when datasets included host-specific outliers, but not when based on neutral loci, which united sympatric populations. Lastly, three host-specific outliers were highly differentiated in all nine different-host comparisons. Because host-adaptation promotes reproductive isolation in these beetles, these loci provide promising candidate gene regions for future molecular studies of ecological speciation.     

Genomic regions repeatedly involved in divergence among plant‐specialized pea aphid biotypes

Understanding the genetic bases of biological diversification is a long‐standing goal in evolutionary biology. Here, we investigate whether replicated cases of adaptive divergence involve the same genomic regions in the pea aphid, Acyrthosiphon pisum, a large complex of genetically differentiated biotypes, each specialized on different species of legumes. A previous study identified genomic regions putatively involved in host‐plant adaptation and/or reproductive isolation by performing a hierarchical genome scan in three biotypes. This led to the identification of 11 F_ST outliers among 390 polymorphic microsatellite markers. In this study, the outlier status of these 11 loci was assessed in eight biotypes specialized on other host plants. Four of the 11 previously identified outliers showed greater genetic differentiation among these additional biotypes than expected under the null hypothesis of neutral evolution (α < 0.01). Whether these hotspots of genomic divergence result from adaptive events, intrinsic barriers or reduced recombination is discussed.     

Heterogeneous genomic differentiation between walking-stick ecotypes: "isolation by adaptation" and multiple roles for divergent selection

Genetic differentiation can be highly variable across the genome. For example, loci under divergent selection and those tightly linked to them may exhibit elevated differentiation compared to neutral regions. These represent "outlier loci" whose differentiation exceeds neutral expectations. Adaptive divergence can also increase genome-wide differentiation by promoting general barriers to neutral gene flow, thereby facilitating genomic divergence via genetic drift. This latter process can yield a positive correlation between adaptive phenotypic divergence and neutral genetic differentiation (described here as "isolation-by-adaptation"). Here, we examine both these processes by combining an AFLP genome scan of two host plant ecotypes of Timema cristinae walking-sticks with existing data on adaptive phenotypic divergence and ecological speciation in these insects. We found that about 8% of loci are outliers in multiple population comparisons. Replicated comparisons between population-pairs using the same versus different host species revealed that 1-2% of loci are subject to host-related selection specifically. Locus-specific analyses revealed that up to 10% of putatively neutral (nonoutlier) AFLP loci exhibit significant isolation-by-adaptation. Our results suggest that selection may affect differentiation directly, via linkage, or by facilitating genetic drift. They thus illustrate the varied and sometimes nonintuitive contributions of selection to heterogeneous genomic differentiation.     

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.     

Isolation by adaptation in Neochlamisus leaf beetles: host-related selection promotes neutral genomic divergence

This study tests how divergent natural selection promotes genomic differentiation during ecological speciation. Specifically, we use adaptive ecological divergence (here, population divergence in host plant use and preference) as a proxy for selection strength and evaluate the correlation between levels of adaptive and genetic differentiation across pairwise population comparisons. Positive correlations would reveal the pattern predicted by our hypothesis, that of 'isolation by adaptation' (IBA). Notably, IBA is predicted not only for selected loci but also for neutral loci. This may reflect the effects of divergent selection on neutral loci that are 'loosely linked' to divergently selected loci or on geneflow restriction that facilitates genetic drift at all loci, including neutral loci that are completely unlinked to those evolving under divergent selection. Here, we evaluate IBA in maple- and willow-associated populations of Neochlamisus bebbianae leaf beetles. To do so, we collected host preference data to construct adaptive divergence indices and used AFLPs (amplified fragment length polymorphisms) and mitochondrial sequences to quantify genetic differentiation. Partial Mantel tests showed significant IBA in 'pooled' analyses of putatively neutral and of putatively selected ('outlier') AFLP loci. This pattern was also recovered in 12% of 'locus-specific' analyses that separately evaluated genetic differentiation at individual neutral loci. These results provided evidence for widespread effects of selection on neutral genomic divergence. Our collective findings indicate that host-related selection may play important roles in the population genomic differentiation of both neutral and selected gene regions in herbivorous insects.     

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.     

Patterns of divergence across the geographic and genomic landscape of a butterfly hybrid zone associated with a climatic gradient

Hybrid zones are a valuable tool for studying the process of speciation and for identifying the genomic regions undergoing divergence and the ecological (extrinsic) and nonecological (intrinsic) factors involved. Here, we explored the genomic and geographic landscape of divergence in a hybrid zone between Papilio glaucus and Papilio canadensis. Using a genome scan of 28,417 ddRAD SNPs, we identified genomic regions under possible selection and examined their distribution in the context of previously identified candidate genes for ecological adaptations. We showed that differentiation was genomewide, including multiple candidate genes for ecological adaptations, particularly those involved in seasonal adaptation and host plant detoxification. The Z chromosome and four autosomes showed a disproportionate amount of differentiation, suggesting genes on these chromosomes play a potential role in reproductive isolation. Cline analyses of significantly differentiated genomic SNPs, and of species‐diagnostic genetic markers, showed a high degree of geographic coincidence (81%) and concordance (80%) and were associated with the geographic distribution of a climate‐mediated developmental threshold (length of the growing season). A relatively large proportion (1.3%) of the outliers for divergent selection were not associated with candidate genes for ecological adaptations and may reflect the presence of previously unrecognized intrinsic barriers between these species. These results suggest that exogenous (climate‐mediated) and endogenous (unknown) clines may have become coupled and act together to reinforce reproductive isolation. This approach of assessing divergence across both the genomic and geographic landscape can provide insight about the interplay between the genetic architecture of reproductive isolation and endogenous and exogenous selection.     

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.     

Population genetic structure and secondary symbionts in host-associated populations of the pea aphid complex

Polyphagous insect herbivores experience different selection pressures on their various host plant species. How this affects population divergence and speciation may be influenced by the bacterial endosymbionts that many harbor. Here, we study the population structure and symbiont community of the pea aphid (Acyrthosiphon pisum), which feeds on a range of legume species and is known to form genetically differentiated host-adapted populations. Aphids were collected from eight legume genera in England and Germany. Extensive host plant associated differentiation was observed with this collection of pea aphids comprising nine genetic clusters, each of which could be associated with a specific food plant. Compared to host plant, geography contributed little to genetic differentiation. The genetic clusters were differentiated to varying degrees, but this did not correlate with their degree of divergence in host use. We surveyed the pea aphid clones for the presence of six facultative (secondary) bacterial endosymbionts and found they were nonrandomly distributed across the aphid genetic clusters and this distribution was similar in the two countries. Aphid clones on average carried 1.4 species of secondary symbiont with those associated with Lathyrus having significantly fewer. The results are interpreted in the light of the evolution of specialization and ecological speciation.     

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction‐site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small N_e (range = 2.1–89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome‐wide divergence. Nonetheless, outlier tests identified 3.6–6.6% of loci as high F_ST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high F_ST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small N_e in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.     

Comparative genome scan detects host‐related divergent selection in the grasshopper Hesperotettix viridis

In this study, we used a comparative genome scan to examine patterns of population differentiation with respect to host plant use in Hesperotettix viridis, a Nearctic oligophagous grasshopper locally specialized on various Asteraceae including Solidago, Gutierrezia, and Ericameria. We identified amplified fragment length polymorphism (AFLP) loci with significantly elevated F_ST (outlier loci) in multiple different‐host and same‐host comparisons of populations while controlling for geographic distance. By comparing the number and identities of outlier loci in different‐host vs. same‐host comparisons, we found evidence of host plant‐related divergent selection for some population comparisons (Solidago‐ vs. Gutierrezia‐feeders), while other comparisons (Ericameria‐ vs. Gutierrezia‐feeders) failed to demonstrate a strong role for host association in population differentiation. In comparisons of Solidago‐ vs. Gutierrezia‐feeding populations, a relatively high number of outlier loci observed repeatedly in different‐host comparisons (35% of all outliers and 2.7% of all 625 AFLP loci) indicated a significant role for host‐related selection in contributing to overall genomic differentiation in this grasshopper. Mitochondrial DNA sequence data revealed a star‐shaped phylogeny with no host‐ or geography‐related structure, low nucleotide diversity, and high haplotype diversity, suggesting a recent population expansion. mtDNA data do not suggest a long period of isolation in separate glacial refugia but are instead more compatible with a single glacial refugium and more recent divergence in host use. Our study adds to research documenting heterogeneity in differentiation across the genome as a consequence of divergent natural selection, a phenomenon that may occur as part of the process of ecological speciation.     

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.     

Integrating QTL mapping and genome scans towards the characterization of candidate loci under parallel selection in the lake whitefish (Coregonus clupeaformis)

As natural selection must act on underlying genetic variation, discovering the number and location of loci under the influence of selection is imperative towards understanding adaptive divergence in evolving populations. Studies employing genome scans have hypothesized that the action of divergent selection should reduce gene flow at the genomic locations implicated in adaptation and speciation among natural populations, yet once 'outlier' patterns of variation have been identified the function and role of such loci needs to be confirmed. We integrated adaptive QTL mapping and genomic scans among diverging sympatric pairs of the lake whitefish (Coregonus clupeaformis) species complex in order to test the hypothesis that differentiation between dwarf and normal ecotypes at growth-associated QTL was maintained by directional selection. We found evidence of significantly high levels of molecular divergence among eight growth QTL where two of the strongest candidate loci under the influence of directional selection exhibited parallel reductions of gene flow over multiple populations.