SELECTION AND GENOMIC DIFFERENTIATION DURING ECOLOGICAL SPECIATION: ISOLATING THE CONTRIBUTIONS OF HOST ASSOCIATION VIA A COMPARATIVE GENOME SCAN OF NEOCHLAMISUS BEBBIANAE LEAF BEETLES

Postmating isolation and genetically variable host use in ecologically divergent host forms of Neochlamisus bebbianae leaf beetles

Ecological speciation studies have more thoroughly addressed premating than postmating reproductive isolation. This study examines multiple postmating barriers between host forms of Neochlamisus bebbianae leaf beetles that specialize on Acer and Salix trees. We demonstrate cryptic isolation and reduced hybrid fitness via controlled matings of these host forms. These findings reveal host-associated postmating isolation, although a nonecological, 'intrinsic' basis for these patterns cannot be ruled out. Host preference and performance results among cross types further suggest sex-linked maternal effects on these traits, whereas family effects indicate their genetic basis and associated variation. Genes of major effect appear to influence these traits. Together with previous findings of premating isolation and adaptive differentiation in sympatry, our results meet many assumptions of 'speciation with gene flow' models. Here, such gene flow is likely asymmetric, with consequences for the dynamics of future ecological divergence and potential ecological speciation of these host forms.

Comparing geographical genetic differentiation between candidate and noncandidate loci for adaptation strengthens support for parallel ecological divergence in the marine snail Littorina saxatilis

The Galician sympatric ecotypes of Littorina saxatilis have been proposed as a model system for studying parallel ecological speciation. Such a model system makes a clear prediction: candidate loci (for divergent adaptation) should present a higher level of geographical differentiation than noncandidate (neutral) loci. We used 2356 amplified fragment length polymorphisms (AFLPs) and four microsatellite loci to identify candidate loci for ecological adaptation using the F(ST) outlier method. Three per cent of the studied AFLP loci were identified as candidate loci associated with adaptation, after multitest adjustments, thus contributing to ecotype differentiation (candidate loci were not detected within ecotypes). Candidate and noncandidate loci were analysed separately at four different F(ST) partitions: differences between ecotypes (overall and local), differences between localities and micro-geographical differences within ecotypes. The magnitude of F(ST) differed between candidate and noncandidate loci for all partitions except in the case of micro-geographical differentiation within ecotypes, and the microsatellites (putatively neutral) showed an identical pattern to noncandidate loci. Thus, variation in candidate loci is determined partially independent by divergent natural selection (in addition to stochastic forces) at each locality, while noncandidate loci are exclusively driven by stochastic forces. These results support the evolutionary history described for these particular populations, considered to be a clear example of incomplete sympatric ecological speciation.

Evidence for selection in response to radiation exposure: Pinus sylvestris in the Chernobyl exclusion zone

Changes of genetic structures due to viability selection are likely to occur in populations exposed to rapidly and extremely changing environmental conditions after catastrophic events. However, very little is known about the extent of selective responses and in particular the proportion of the genome involved in putatively adaptive reactions for non-model plants. We used amplified fragment length polymorphisms (AFLPs) in order to investigate genetic differences between pine (Pinus sylvestris) trees which were partially exposed to extreme environmental conditions. Genetic variation patterns of pines exposed to high radiation in the Chernobyl exclusion zone with or without phenotypic stress symptoms were compared to control trees with a similar origin. Six percent of the investigated loci (15 of 222 loci) were identified as candidates for selective responses. Moderate differentiation was observed between groups of trees showing either weak or strong phenotypic responses to high radiation levels.

Broken barriers: human-induced changes to gene flow and introgression in animals: an examination of the ways in which humans increase genetic exchange among populations and species and the consequences for biodiversity

We identify two processes by which humans increase genetic exchange among groups of individuals: by affecting the distribution of groups and dispersal patterns across a landscape, and by affecting interbreeding among sympatric or parapatric groups. Each of these processes might then have two different effects on biodiversity: changes in the number of taxa through merging or splitting of groups, and the extinction/extirpation of taxa through effects on fitness. We review the various ways in which humans are affecting genetic exchange, and highlight the difficulties in predicting the impacts on biodiversity. Gene flow and hybridization are crucially important evolutionary forces influencing biodiversity. Humans alter natural patterns of genetic exchange in myriad ways, and these anthropogenic effects are likely to influence the genetic integrity of populations and species. We argue that taking a gene-centric view towards conservation will help resolve issues pertaining to conservation and management. Editor's suggested further reading in BioEssays A systemic view of biodiversity and its conservation: Processes, interrelationships, and human culture Abstract.

Scanning the European corn borer (Ostrinia spp.) genome for adaptive divergence between host-affiliated sibling species

It has recently been shown that the European corn borer, a major pest of maize crops, is actually composed of two genetically differentiated and reproductively isolated taxa, which are found in sympatry over a wide geographical range in Eurasia. Each taxon is adapted to specific host plants: Ostrinia nubilalis feeds mainly on maize, while O. scapulalis feeds mainly on hop or mugwort. Here, we present a genome scan approach as a first step towards an integrated molecular analysis of the adaptive genomic divergence between O. nubilalis and O. scapulalis. We analysed 609 AFLP marker loci in replicate samples of sympatric populations of Ostrinia spp. collected on maize, hop and mugwort, in France. Using two genome scan methods based on the analysis of population differentiation, we found a set of 28 outlier loci that departed from the neutral expectation in one or the other method (of which a subset of 14 loci were common to both methods), which showed a significantly increased differentiation between O. nubilalis and O. scapulalis, when compared to the rest of the genome. A subset of 12 outlier loci were sequenced, of which 7 were successfully re-amplified as target candidate loci. Three of these showed homology with annotated lepidopteran sequences from public nucleotide databases.

Genomic scans detect signatures of selection along a salinity gradient in populations of the intertidal seaweed Fucus serratus on a 12 km scale

Detecting natural selection in wild populations is a central challenge in evolutionary biology and genomic scans are an important means of detecting allele frequencies that deviate from neutral expectations among marker loci. We used nine anonymous and 15 EST-linked microsatellites, 362 AFLP loci, and several neutrality tests, to identify outlier loci when comparing four populations of the seaweed Fucus serratus spaced along a 12km intertidal shore with a steep salinity gradient. Under criteria of at least two significant tests in at least two population pairs, three EST-derived and three anonymous loci revealed putative signatures of selection. Anonymous locus FsB113 was a consistent outlier when comparing least saline to fully marine sites. Locus F37 was an outlier when comparing the least saline to more saline areas, and was annotated as a polyol transporter/putative mannitol transporter - an important sugar-alcohol associated with osmoregulation by brown algae. The remaining loci could not be annotated using six different data bases. Exclusion of microsatellite outlier loci did not change either the degree or direction of differentiation among populations. In one outlier test, the number of AFLP outlier loci increased as the salinity differences between population pairs increased (up to 14); only four outliers were detected with the second test and only one was consistent with both tests. Consistency may be improved with a much more rigorous approach to replication and/or may be dependent upon the class of marker used.

Multiple approaches to detect outliers in a genome scan for selection in ocellated lizards (Lacerta lepida) along an environmental gradient

Identification of loci with adaptive importance is a key step to understand the speciation process in natural populations, because those loci are responsible for phenotypic variation that affects fitness in different environments. We conducted an AFLP genome scan in populations of ocellated lizards (Lacerta lepida) to search for candidate loci influenced by selection along an environmental gradient in the Iberian Peninsula. This gradient is strongly influenced by climatic variables, and two subspecies can be recognized at the opposite extremes: L. lepida iberica in the northwest and L. lepida nevadensis in the southeast. Both subspecies show substantial morphological differences that may be involved in their local adaptation to the climatic extremes. To investigate how the use of a particular outlier detection method can influence the results, a frequentist method, DFDIST, and a Bayesian method, BayeScan, were used to search for outliers influenced by selection. Additionally, the spatial analysis method was used to test for associations of AFLP marker band frequencies with 54 climatic variables by logistic regression. Results obtained with each method highlight differences in their sensitivity. DFDIST and BayeScan detected a similar proportion of outliers (3-4%), but only a few loci were simultaneously detected by both methods. Several loci detected as outliers were also associated with temperature, insolation or precipitation according to spatial analysis method. These results are in accordance with reported data in the literature about morphological and life-history variation of L. lepida subspecies along the environmental gradient.

The efficacy of divergence hitchhiking in generating genomic islands during ecological speciation

Genes under divergent selection flow less readily between populations than other loci. This observation has led to verbal "divergence hitchhiking" models of speciation in which decreased interpopulation gene flow surrounding loci under divergent selection can generate large regions of differentiation within the genome (genomic islands). The efficacy of this model in promoting speciation depends on the size of the region affected by divergence hitchhiking. Empirical evidence is mixed, with examples of both large and small genomic islands. To address these empirical discrepancies and to formalize the theory, we present mathematical models of divergence hitchhiking, which examine neutral differentiation around selected sites. For a single locus under selection, regions of differentiation do not extend far along a chromosome away from a selected site unless both effective population sizes and migration rates are low. When multiple loci are considered, regions of differentiation can be larger. However, with many loci under selection, genome-wide divergence occurs and genomic islands are erased. The results show that divergence hitchhiking can generate large regions of differentiation, but that the conditions under which this occurs are limited. Thus, speciation may often require multifarious selection acting on many, isolated and physically unlinked genes. How hitchhiking promotes further adaptive divergence warrants consideration.

A high incidence of selection on physiologically important genes in the three-spined stickleback, Gasterosteus aculeatus

Genome scan approaches to detect footprints of directional selection in the genomes of wild animal and plant populations have become popular tools to study local adaptation and speciation at the molecular level. Most studies thus far have used random molecular markers and found footprints of directional selection at, on average, 5% (range: 1-15%) of the examined loci. We focused on physiologically important genes that exhibit transcriptional responses to specific environmental or developmental conditions and assessed if these genes have been subject to directional selection and are responsible for local adaptation in the three-spined stickleback (Gasterosteus aculeatus). Using microsatellite markers located within or closely linked to (<6 kb) target genes, we investigated footprints of directional selection for 157 genes with known physiological functions in three marine and six freshwater populations. A high incidence (16.6%) of footprints of directional selection for these genes was revealed by four different outlier tests. In a subset of four populations screened with both physiologically important and random genes, footprints of directional selection were more frequent in physiologically important genes (13.4%) as compared with random genes (2.4%). In general, our findings indicate strong selective pressures on physiologically important genes, suggesting that these genes have significant functions in evolutionary adaptation to environmental heterogeneity.

Genetic population structure of sympatric and allopatric populations of Baltic ciscoes (Coregonus albula complex, Teleostei, Coregonidae)

Background

Teleost fishes of the Coregonidae are good model systems for studying postglacial evolution, adaptive radiation and ecological speciation. Of particular interest is whether the repeated occurrence of sympatric species pairs results from in-situ divergence from a single lineage or from multiple invasions of one or more different lineages. Here, we analysed the genetic structure of Baltic ciscoes (Coregonus albula complex), examining 271 individuals from 8 lakes in northern Germany using 1244 polymorphic AFLP loci. Six lakes had only one population of C. albula while the remaining two lakes had C. albula as well as a sympatric species (C. lucinensis or C. fontanae).

Results

AFLP demonstrated a significant population structure (Bayesian θ^B = 0.22). Lower differentiation between allopatric (θ^B = 0.028) than sympatric (0.063-0.083) populations contradicts the hypothesis of a sympatric origin of taxa, and there was little evidence for stocking or ongoing hybridization. Genome scans found only three loci that appeared to be under selection in both sympatric population pairs, suggesting a low probability of similar mechanisms of ecological segregation. However, removal of all non-neutral loci decreased the genetic distance between sympatric pairs, suggesting recent adaptive divergence at a few loci. Sympatric pairs in the two lakes were genetically distinct from the six other C. albula populations, suggesting introgression from another lineage may have influenced these two lakes. This was supported by an analysis of isolation-by-distance, where the drift-gene flow equilibrium observed among allopatric populations was disrupted when the sympatric pairs were included.

Conclusions

While the population genetic data alone can not unambiguously uncover the mode of speciation, our data indicate that multiple lineages may be responsible for the complex patterns typically observed in Coregonus. Relative differences within and among lakes raises the possibility that multiple lineages may be present in northern Germany, thus understanding the postglacial evolution and speciation in the C. albula complex requires a large-scale phylogenetic analysis of several potential founder lineages.

Sympatric ecological speciation meets pyrosequencing: sampling the transcriptome of the apple maggot Rhagoletis pomonella

Background

The full power of modern genetics has been applied to the study of speciation in only a small handful of genetic model species - all of which speciated allopatrically. Here we report the first large expressed sequence tag (EST) study of a candidate for ecological sympatric speciation, the apple maggot Rhagoletis pomonella, using massively parallel pyrosequencing on the Roche 454-FLX platform. To maximize transcript diversity we created and sequenced separate libraries from larvae, pupae, adult heads, and headless adult bodies.

Results

We obtained 239,531 sequences which assembled into 24,373 contigs. A total of 6810 unique protein coding genes were identified among the contigs and long singletons, corresponding to 48% of all known Drosophila melanogaster protein-coding genes. Their distribution across GO classes suggests that we have obtained a representative sample of the transcriptome. Among these sequences are many candidates for potential R. pomonella "speciation genes" (or "barrier genes") such as those controlling chemosensory and life-history timing processes. Furthermore, we identified important marker loci including more than 40,000 single nucleotide polymorphisms (SNPs) and over 100 microsatellites. An initial search for SNPs at which the apple and hawthorn host races differ suggested at least 75 loci warranting further work. We also determined that developmental expression differences remained even after normalization; transcripts expected to show different expression levels between larvae and pupae in D. melanogaster also did so in R. pomonella. Preliminary comparative analysis of transcript presences and absences revealed evidence of gene loss in Drosophila and gain in the higher dipteran clade Schizophora.

Conclusions

These data provide a much needed resource for exploring mechanisms of divergence in this important model for sympatric ecological speciation. Our description of ESTs from a substantial portion of the R. pomonella transcriptome will facilitate future functional studies of candidate genes for olfaction and diapause-related life history timing, and will enable large scale expression studies. Similarly, the identification of new SNP and microsatellite markers will facilitate future population and quantitative genetic studies of divergence between the apple and hawthorn-infesting host races.

Genome scan to assess the respective role of host-plant and environmental constraints on the adaptation of a widespread insect

Background

The evolutionary success of phytophagous insects could result from their adaptation to different host-plants. Alternatively, the diversification of widespread species might be driven by adaptation along environmental gradients. To disentangle the respective roles of host-plant versus abiotic environmental variables acting on the genome of an oligophagous insect, we performed a genome scan using 83 unlinked AFLP markers on larvae of the large pine weevil collected on two host-plants (pine and spruce) in four forestry regions across Europe.

Results

At this large geographic scale, the global genetic differentiation was low and there was no isolation by distance pattern, suggesting that migration is overwhelming genetic drift in this species. In this context, the widely used frequentist methods to detect outliers (e.g. Dfdist), which assume migration - drift equilibrium are not the most appropriate approach. The implementation of a recently developed Bayesian approach, conceived to detect outliers even in non-equilibrium situations, consistently detected 9 out of 83 loci as outliers. Eight of these were validated as outliers by multiple logistic regressions: six correlated with environmental variables, one with host-plant and one with the interaction between environmental variables and host-plant.

Conclusion

These results suggest a relatively greater importance of abiotic environmental variables, as opposed to factors linked with the host-plant, in shaping genetic differentiation across the genome in this species. Logistic regression allows the nature of factors involved in locus-specific selection to be precisely identified and represents another step forward in the process of identifying adaptive loci.

Asymmetric introgression between sympatric molestus and pipiens forms of Culex pipiens (Diptera: Culicidae) in the Comporta region, Portugal

Background

Culex pipiens L. is the most widespread mosquito vector in temperate regions. This species consists of two forms, denoted molestus and pipiens, that exhibit important behavioural and physiological differences. The evolutionary relationships and taxonomic status of these forms remain unclear. In northern European latitudes molestus and pipiens populations occupy different habitats (underground vs. aboveground), a separation that most likely promotes genetic isolation between forms. However, the same does not hold in southern Europe where both forms occur aboveground in sympatry. In these southern habitats, the extent of hybridisation and its impact on the extent of genetic divergence between forms under sympatric conditions has not been clarified. For this purpose, we have used phenotypic and genetic data to characterise Cx. pipiens collected aboveground in Portugal. Our aims were to determine levels of genetic differentiation and the degree of hybridisation between forms occurring in sympatry, and to relate these with both evolutionary and epidemiological tenets of this biological group.

Results

Autogeny and stenogamy was evaluated in the F1 progeny of 145 individual Cx. pipiens females. Bayesian clustering analysis based on the genotypes of 13 microsatellites revealed two distinct genetic clusters that were highly correlated with the alternative traits that define pipiens and molestus. Admixture analysis yielded hybrid rate estimates of 8-10%. Higher proportions of admixture were observed in pipiens individuals suggesting that more molestus genes are being introgressed into the pipiens form than the opposite.

Conclusion

Both physiological/behavioural and genetic data provide evidence for the sympatric occurrence of molestus and pipiens forms of Cx. pipiens in the study area. In spite of the significant genetic differentiation between forms, hybridisation occurs at considerable levels. The observed pattern of asymmetric introgression probably relates to the different mating strategies adopted by each form. Furthermore, the differential introgression of molestus genes into the pipiens form may induce a more opportunistic biting behaviour in the latter thus potentiating its capacity to act as a bridge-vector for the transmission of arboviral infections.

Ecologically dependent postmating isolation between sympatric host forms of Neochlamisus bebbianae leaf beetles

Ecological speciation is the promotion of reproductive isolation via the divergent adaptation of populations to alternative environments. A prediction peculiar to ecological speciation is that hybrids between such populations should be adapted poorly to parental environments, yielding reduced fitness and postmating isolation. However, F(1) analyses alone cannot demonstrate that ecological ("extrinsic") factors contribute to such isolation. Rather, this requires documenting a "switch" in the relative fitnesses of reciprocal backcrosses between environments. Specifically, each backcross should exhibit higher fitness in the environment of its pure parent, with which it shares the most genes, including environment-specific ones. In contrast, because genetic proportions are expected to be similar for all backcrosses ( approximately (3/4) from one parental type and approximately (1/4) from the other), the more general genetic incompatibilities responsible for "intrinsic" isolation predict no such environment-specific fitness switches. Thus, although intrinsic isolation may contribute to the fitness reduction and variation underlying such patterns, it offers an insufficient explanation for them. Here, we present a quantitative genetic "backcross" analysis of sympatric Neochlamisus bebbianae leaf beetle populations adapted to maple versus willow host plants. Results statistically supported ecological speciation predictions, notably the switch in relative fitness for backcross types, the expected rank order of cross type fitnesses, and appreciable extrinsic isolation. We additionally documented genetic variation in host-associated fitness, ruled out nongenetic maternal effects, and discuss the maintenance of ecological differentiation in sympatry. In summary, our study provides a rare and strongly supported demonstration of genetically based, ecologically dependent postmating isolation during ecological speciation.

Chromosomal inversions and species differences: when are genes affecting adaptive divergence and reproductive isolation expected to reside within inversions?

Many factors can promote speciation, and one which has received much attention is chromosomal inversions. A number of models propose that the recombination suppressing effects of inversions facilitate the maintenance of differences between interbreeding populations in genes affecting adaptive divergence and reproductive isolation. These models predict that such genes will disproportionately reside within inversions, rather than in collinear regions. This hypothesis has received some support, but exceptions exist. Additionally, the effects of known low levels of recombination within inversions on these models are uninvestigated. Here, simulations are used to compare the maintenance of genetic differences between populations following secondary contact and hybridization in different inversion models. We compare regions with no recombination within them to regions with low recombination and to collinear regions with free recombination. Our most general finding is that the low levels of recombination within an inversion often result in the loss of accentuated divergence in inverted regions compared to collinear ones. We conclude that inversions can facilitate the maintenance of species differences under some conditions, but that large or qualitative differences between inverted and collinear regions need not occur. We also find that strong selection facilitates maintenance of divergence in a manner analogous to inversions.

Detecting loci under selection in a hierarchically structured population

Patterns of genetic diversity between populations are often used to detect loci under selection in genome scans. Indeed, loci involved in local adaptations should show high F(ST) values, whereas loci under balancing selection should rather show low F(ST) values. Most tests of selection based on F(ST) use a null distribution generated under a simple island model of population differentiation. Although this model has been shown to be robust, many species have a more complex genetic structure, with some populations sharing a recent ancestry or due to the presence of barriers to gene flow between different parts of a species range. In this paper, we propose the use of a hierarchical island model, in which demes exchange more migrants within groups than between groups, to generate the joint distribution of genetic diversity within and between populations. We show that tests not accounting for a hierarchical structure, when it exists, do generate a large excess of false positive loci, whereas the hierarchical island model is robust to uncertainties about the exact number of groups and demes per group in the system. Our approach also explicitly takes into account the mutational process, and does not just rely on allele frequencies, which is important for short tandem repeat (STR) data. An application to human and stickleback STR data sets reveals a much lower number of significant loci than previously obtained under a non-hierarchical model. The elimination of false positive loci from genome scans should allow us to better determine on which specific class of genes selection is operating.

Genetic differentiation between resistance phenotypes in the phytophagous flea beetle, Phyllotreta nemorum

The flea beetle Phyllotreta nemorum L. (Coleoptera: Chrysomelidae) is genetically polymorphic for resistance against the defences of one of its host plants, Barbarea vulgaris R.Br. (Brassicales: Brassicaceae). Whereas resistant flea beetles are able to use B. vulgaris as well as other cruciferous plants as food, non-resistant beetles cannot survive on B. vulgaris. This limitation to host plant use of non-resistant beetles could potentially lead to asymmetric gene flow and some degree of genetic isolation between the different resistance-genotypes. Therefore, we studied the extent of genetic differentiation at neutral allozyme loci between samples of flea beetles that were collected at different locations and first tested for resistance phenotype. Since earlier work has shown a weak, but significant, effect of geographical distance between the samples on their genetic differentiation, in the present study variation at the neutral allozyme loci in P. nemorum was partitioned between geographical distance and resistance-phenotype. Both sources independently contributed statistically significantly to population differentiation. Thus, there appears to be a limitation to genetic exchange between the resistant and non-resistant flea beetles when corrections are made for their geographic differentiation. This is consistent with the presence of some degree of host race formation in this flea beetle.

Geographic variation in adaptation at the molecular level: a case study of plant immunity genes

Natural selection imposed by interacting species frequently varies among geographic locations and can lead to local adaptation, where alternative phenotypes are found in different populations. Little is known, however, about whether geographically variable selection acting on traits that mediate species interactions is consistent or strong enough to influence patterns of nucleotide variation at individual loci. To investigate this question, we examined patterns of nucleotide diversity and population structure at 16 plant innate immunity genes, with putative functions in defending plants against pathogens or herbivores, from six populations of teosinte (Zea mays ssp. parviglumis). Specifically, we tested whether patterns of population structure and within-population diversity at immunity genes differed from patterns found at nonimmunity (reference) loci and from neutral expectations derived from coalescent simulations of structured populations. For the majority of genes, we detected no strong evidence of geographically variable selection. However, in the wound-induced serine protease inhibitor (wip1), which inhibits the hydrolysis of dietary proteins in insect herbivores, one population showed unusually high levels of genetic differentiation, very low levels of nucleotide polymorphism, and was fixed for a novel replacement substitution in the active site of the protein. Taken together, these data suggest that wip1 experienced a recent selective sweep in one geographic region; this pattern may reflect local adaptation or an ongoing species-wide sweep. Overall, our results indicate that a signature of local adaptation at the molecular level may be uncommon-particularly for traits that are under complex genetic control.