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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

The role of structural genomic variants in population differentiation and ecotype formation in Timema cristinae walking sticks

Theory predicts that structural genomic variants such as inversions can promote adaptive diversification and speciation. Despite increasing empirical evidence that adaptive divergence can be triggered by one or a few large inversions, the degree to which widespread genomic regions under divergent selection are associated with structural variants remains unclear. Here we test for an association between structural variants and genomic regions that underlie parallel host‐plant‐associated ecotype formation in Timema cristinae stick insects. Using mate‐pair resequencing of 20 new whole genomes we find that moderately sized structural variants such as inversions, deletions and duplications are widespread across the genome, being retained as standing variation within and among populations. Using 160 previously published, standard‐orientation whole genome sequences we find little to no evidence that the DNA sequences within inversions exhibit accentuated differentiation between ecotypes. In contrast, a formerly described large region of reduced recombination that harbours genes controlling colour‐pattern exhibits evidence for accentuated differentiation between ecotypes, which is consistent with differences in the frequency of colour‐pattern morphs between host‐associated ecotypes. Our results suggest that some types of structural variants (e.g., large inversions) are more likely to underlie adaptive divergence than others, and that structural variants are not required for subtle yet genome‐wide genetic differentiation with gene flow.     

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.     

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.     

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.     

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.     

Insight into the roles of selection in speciation from genomic patterns of divergence and introgression in secondary contact in venomous rattlesnakes

Investigating secondary contact of historically isolated lineages can provide insight into how selection and drift influence genomic divergence and admixture. Here, we studied the genomic landscape of divergence and introgression following secondary contact between lineages of the Western Diamondback Rattlesnake (Crotalus atrox) to determine whether genomic regions under selection in allopatry also contribute to reproductive isolation during introgression. We used thousands of nuclear loci to study genomic differentiation between two lineages that have experienced recent secondary contact following isolation, and incorporated sampling from a zone of secondary contact to identify loci that are resistant to gene flow in hybrids. Comparisons of patterns of divergence and introgression revealed a positive relationship between allelic differentiation and resistance to introgression across the genome, and greater‐than‐expected overlap between genes linked to lineage‐specific divergence and loci that resist introgression. Genes linked to putatively selected markers were related to prominent aspects of rattlesnake biology that differ between populations of Western Diamondback rattlesnakes (i.e., venom and reproductive phenotypes). We also found evidence for selection against introgression of genes that may contribute to cytonuclear incompatibility, consistent with previously observed biased patterns of nuclear and mitochondrial alleles suggestive of partial reproductive isolation due to cytonuclear incompatibilities. Our results provide a genome‐scale perspective on the relationships between divergence and introgression in secondary contact that is relevant for understanding the roles of selection in maintaining partial isolation of lineages, causing admixing lineages to not completely homogenize.     

A test of genomic modularity among life‐history adaptations promoting speciation with gene flow

Speciation with gene flow may require adaptive divergence of multiple traits to generate strong ecologically based reproductive isolation. Extensive negative pleiotropy or physical linkage of genes in the wrong phase affecting these diverging traits may therefore hinder speciation, while genetic independence or “modularity” among phenotypic traits may reduce constraints and facilitate divergence. Here, we test whether the genetics underlying two components of diapause life history, initial diapause intensity and diapause termination timing, constrain differentiation between sympatric hawthorn and apple‐infesting host races of the fly Rhagoletis pomonella through analysis of 10,256 SNPs measured via genotyping‐by‐sequencing (GBS). Loci genetically associated with diapause termination timing were mainly observed for SNPs mapping to chromosomes 1–3 in the genome, most notably for SNPs displaying higher levels of linkage disequilibrium (LD), likely due to inversions. In contrast, selection on initial diapause intensity affected loci on all five major chromosomes of the genome, specifically those showing low levels of LD. This lack of overlap in genetically associated loci suggests that the two diapause phenotypes are largely modular. On chromosome 2, however, intermediate level LD loci and a subgroup of high LD loci displayed significant negative relationships between initial diapause intensity and diapause termination time. These gene regions on chromosome 2 therefore affected both traits, while most regions were largely independent. Moreover, loci associated with both measured traits also tended to exhibit highly divergent allele frequencies between the host races. Thus, the presence of nonoverlapping genetic modules likely facilitates simultaneous, adaptive divergence for the measured life‐history components.     

Experimental evidence of genome‐wide impact of ecological selection during early stages of speciation‐with‐gene‐flow

Theory predicts that speciation‐with‐gene‐flow is more likely when the consequences of selection for population divergence transitions from mainly direct effects of selection acting on individual genes to a collective property of all selected genes in the genome. Thus, understanding the direct impacts of ecologically based selection, as well as the indirect effects due to correlations among loci, is critical to understanding speciation. Here, we measure the genome‐wide impacts of host‐associated selection between hawthorn and apple host races of Rhagoletis pomonella (Diptera: Tephritidae), a model for contemporary speciation‐with‐gene‐flow. Allele frequency shifts of 32 455 SNPs induced in a selection experiment based on host phenology were genome wide and highly concordant with genetic divergence between co‐occurring apple and hawthorn flies in nature. This striking genome‐wide similarity between experimental and natural populations of R. pomonella underscores the importance of ecological selection at early stages of divergence and calls for further integration of studies of eco‐evolutionary dynamics and genome divergence.     

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.     

Reproductive isolation between divergent races of pea aphids on two hosts. II. Selection against migrants and hybrids in the parental environments

Sympatric races of pea aphids on alfalfa and red clover are highly ecologically specialized and significantly reproductively isolated. Much of the restriction of gene flow between the specialized populations is due to habitat choice behavior of the winged colonizers (Via 1999). Here, we document additional pre- and postmating reproductive isolation through selection against migrants and hybrids in the parental environments. First, a group of randomly chosen genotypes from each race that were experimentally migrated between hosts had very low survival and reproduction on the alternate host relative to genotypes originating from that host (natives). Such selection against cross-host migrants forms a premating barrier to gene flow because it is likely to reduce migrant frequencies before the sexual forms are induced in the fall. Our reciprocal transplant experiment also shows that natural selection acts directly on individual migrants between the crops to favor host choice behavior: genotypes from each host suffered large losses of fitness when forced to migrate to the alternate host plant relative to the fitness they would have enjoyed had they been able to choose their native host. In a companion field study, sequential sampling throughout the summer in newly colonized fields of both alfalfa and clover revealed a decrease in the frequency of host-specific marker alleles characteristic of the alternate crop. These field data further support the hypothesis that selection disfavors migrants that cross between crops. Second, when two sets of F1 hybrids between the races were reciprocally tested on alfalfa and clover, both sets had significantly lower average fitness than the specialized parent in each of the two environments. This demographic selection against hybrids in the parental environments is a source of postmating reproductive isolation between the specialized races. Finally, significant genetic variation in fitness traits was seen among F1 hybrid genotypes from both crosses between alfalfa and clover specialists. Although this variation suggests that a generalized pea aphid could evolve, such generalists are not seen in field collections of these populations.     

Complex trait differentiation between host-populations of the pea aphid Acyrthosiphon pisum (Harris): implications for the evolution of ecological specialisation

Variation in traits affecting preference for, and performance on, new habitats is a key factor in the initiation of ecological specialisation and adaptive speciation. However, habitat and resource use also involves other traits whose influence on ecological and genetic divergence remains poorly understood. In the present study, we investigated the extent of variation of life-history traits among sympatric populations of the pea aphid Acyrthosiphon pisum , which shows several host races that are specialised on various plants of the family Fabaceae plants and is an established model for ecological speciation. First, we assessed the community structure of microbial partners within host populations of the pea aphid. The effect of these microbes on host fitness is uncertain, although there is growing evidence that they may modulate various important adaptive traits of their host such as plant utilisation and resistance against natural enemies. Second, we performed a multivariate analysis on several ecologically relevant features of host populations recorded in the present and previous studies (including microbial composition, colour morph, reproductive mode, and male dispersal phenotype), enabling the identification of correlations between phenotypic traits. We discuss the ecological significance of these associations of traits in relation to the habitat characteristics of pea aphid populations, and their consequences for the evolution of ecological specialisation and sympatric speciation.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 718–727.