DIVERSIFICATION AND GENE FLOW IN NASCENT LINEAGES OF ISLAND AND MAINLAND NORTH AMERICAN TREE SQUIRRELS (TAMIASCIURUS)

Pleistocene climate cycles and glaciations had profound impacts on taxon diversification in the Boreal Forest Biome. Using population genetic analyses with multilocus data, we examined diversification, isolation, and hybridization in two sibling species of tree squirrels (Tamiasciurus douglasii and Tamiasciurus hudsonicus) with special attention to the geographically and genetically enigmatic population of T. hudsonicus on Vancouver Island, Canada. The two species differentiated only about 500,000 years ago, in the Late Pleistocene. The island population is phylogenetically nested within T. hudsonicus according to our nuclear analysis but within T. douglasii according to mitochondrial DNA. This conflict is more likely due to historical hybridization than to incomplete lineage sorting, and it appears that bidirectional gene flow occurred between the island population and both species on the mainland. This interpretation of our genetic analyses is consistent with our bioclimatic modeling, which demonstrates that both species were able to occupy this region throughout the Late Pleistocene. The divergence of the island population 40,000 years ago suggests that tree squirrels persisted in a refugium on Vancouver Island at the last glacial maximum, 20,000 years ago. Our observations demonstrate how Pleistocene climate change and habitat shifts have created incipient divergence in the presence of gene flow.     

Clinal colour variation within a panmictic population of tree squirrels,Tamiasciurus douglasii (Rodentia: Sciuridae), across an ecological gradient

Local adaptation occurs when a population in a heterogeneous environment experiences divergent ecological selection but only if selection is stronger than the homogenizing effects of gene flow. The forest environments of Oregon vary along a physical and biotic gradient from a wet, closed‐canopy forest near the coast to a drier open‐canopy forest eastward across the Cascade Mountains. The present study explores patterns of local adaptation in Douglas squirrels (Tamiasciurus douglasii) in relation to these transitions in forest structure and ecology. We test for the presence of morphological clines in relation to gene flow and, more specifically, whether any such character clines correspond with environmental clines. We sampled animals at six locations (10 specimens each) and evaluated environmental parameters across a 240‐km west‐to‐east transect. Population structure analysis of 18 microsatellite loci indicates a single, panmictic squirrel population across the entire transect. Coalescent‐based estimates show bidirectional gene flow at similar west–east intensities between squirrels in coastal and interior forests. Of the four skull traits examined, none shows a significant clinal transition. By contrast, ventral fur colour shows a strong clinal transition, from deep‐orange in coastal forest to whitish–yellow in the interior forest. This pattern of phenotypic divergence coincides with the gradient in tree‐canopy cover. Ventral fur colour of T. douglasii exemplifies a gradation of continuous phenotypic variation maintained despite ongoing gene flow in a panmictic population. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 536–546.     

Genomic variation across two barn swallow hybrid zones reveals traits associated with divergence in sympatry and allopatry

Hybrid zones are geographic regions where isolating barriers between divergent populations are challenged by admixture. Identifying factors that facilitate or inhibit hybridization in sympatry can illuminate the processes that maintain those reproductive barriers. We analysed patterns of hybridization and phenotypic variation across two newly discovered hybrid zones between three subspecies of barn swallow (Hirundo rustica). These subspecies differ in ventral coloration and wing length, traits that are targets of sexual and natural selection, respectively, and are associated with genome‐wide differentiation in allopatry. We tested the hypothesis that the degree of divergence in these traits is associated with the extent of hybridization in secondary contact. We applied measures of population structure based on >23,000 SNPs to confirm that named subspecies correspond to distinct genomic clusters, and assessed coincidence between geographic clines for ancestry and phenotype. Although gene flow was ongoing across both hybrid zones and pairwise F_ST between subspecies was extremely low, we found striking differences in the extent of hybridization. In the more phenotypically differentiated subspecies pair, clines for ancestry, wing length and ventral coloration were steep and coincident, suggestive of strong isolation and, potentially, selection associated with phenotype. In the less phenotypically differentiated pair, gene flow and phenotypic variation occurred over a wide geographic span, indicative of weaker isolation. Traits associated with genome‐wide differentiation in allopatry may thus also contribute to isolation in sympatry. We discuss potentially important additional roles for evolutionary history and ecology in shaping variation in the extent hybridization between closely related pairs of subspecies.     

RECENT ECOLOGICAL SELECTION ON REGULATORY DIVERGENCE IS SHAPING CLINAL VARIATION IN SENECIO ON MOUNT ETNA

The hybrid zone on Mount Etna (Sicily) between Senecio aethnensis and Senecio chrysanthemifolius (two morphologically and physiologically distinct species) is a classic example of an altitudinal cline. Hybridization at intermediate altitudes and gradients in phenotypic and life‐history traits occur along altitudinal transects of the volcano. The cline is considered to be a good example of ecological selection with species differences arising by divergent selection opposing gene flow. However, the possibility that the cline formed from recent secondary contact following an allopatric phase is difficult to exclude. We demonstrate a recent split between S. aethnensis and S. chrysanthemifolius (as recent as ～32,000 years ago) and sufficient gene flow (2Nm > 1) to have prevented divergence (implicating a role for diversifying selection in the maintenance of the cline). Differentially expressed genes between S. aethnensis and S. chrysanthemifolius exhibit significantly higher genetic divergence relative to “expression invariant” controls, suggesting that species differences may in part be mediated by divergent selection on differentially expressed genes involved with altitude‐related adaptation. The recent split time and the absence of fixed differences between these two ecologically distinct species suggest the rapid evolution to an altitudinal cline involving selection on both sequence and expression variation.     

Geographic variation in sexual dimorphism in the barn owl Tyto alba: a role for direct selection or genetic correlation?

Geographic variation in sexually selected traits is commonly attributed to geographic variation in the net benefit accrued from bearing such traits. Although natural and sexual selection are potentially important in shaping geographic variation, genetic constraints may also play a role. Although a genetic correlation between two traits may itself be the outcome of natural or sexual selection, it may indirectly reinforce the establishment and maintenance of cline variation with respect to one particular trait when across the cline different values of other traits are selected. Using the barn owl Tyto alba, a species in which the plumage of females is more reddish‐brown and more marked with black spots than that of males, I report results that are consistent with the hypothesis that both direct selection and genetic constraints may help establish and maintain cline variation in sexual dichromatism. In this species, inter‐individual variation in plumage coloration and spottiness has a genetic basis, and these traits are not sensitive to the environment. Data, based on the measurement of skin specimens, is consistent with the hypothesis that the stronger European cline variation in male spottiness than in female spottiness depends on the combined effects of (1) the similar cline variation in male and female plumage coloration and (2) the more intense phenotypic correlation between plumage coloration and spottiness in males (darker birds are more heavily spotted in the two sexes, but especially males) which is a general feature among the globally distributed barn owls. In northern Europe, male and female T. a. guttata are reddish‐brown and heavily spotted, and in southern Europe male and female T. a. alba are white, but only females display many spots. Here, I discuss the relative importance of direct selection, genetic correlation and the post‐ice age invasion of Europe by T. alba, in generating sex‐specific cline variation in plumage spottiness and non‐sex‐specific cline variation in plumage coloration.     

Heterogeneity and concordance in locus‐specific differentiation and introgression between species of towhees

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well‐differentiated bird species in Mexico (Pipilo maculatus and P. ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping‐by‐sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41 000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome‐wide expectation using a Bayesian genomic cline model. We also estimated locus‐specific F_ST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus‐specific introgression in the two independent hybrid zones (6–20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.     

Isolation by distance,not incipient ecological speciation,explains genetic differentiation in an Andean songbird (Aves: Furnariidae: Cranioleuca antisiensis,Line‐cheeked Spinetail) despite near threefold body size change across an environmental gradient

During the process of ecological speciation, reproductive isolation results from divergent natural selection and leads to a positive correlation between genetic divergence and adaptive phenotypic divergence, that is, isolation by adaptation (IBA). In natural populations, phenotypic differentiation is often autocorrelated with geographic distance, making IBA difficult to distinguish from the neutral expectation of isolation by distance (IBD). We examined these two alternatives in a dramatic case of clinal phenotypic variation in an Andean songbird, the Line‐cheeked Spinetail (Cranioleuca antisiensis). At its geographic extremes, this species shows a near threefold difference in body mass (11.5 to 31.0 g) with marked plumage differences. We analysed phenotypic, environmental and genetic data (5,154 SNPs) from 172 individuals and 19 populations sampled along its linear distribution in the Andes. We found that body mass was tightly correlated with environmental temperature, consistent with local adaptation as per Bergmann's rule. Using a P_ST–F_ST analysis, we found additional support for natural selection driving body mass differentiation, but these results could also be explained by environment‐mediated phenotypic plasticity. When we assessed the relative support for patterns of IBA and IBD using variance partitioning, we found that IBD was the best explanation for genetic differentiation along the cline. Adaptive phenotypic or environmental divergence can reduce gene flow, a pattern interpreted as evidence of ecological speciation's role in diversification. Our results provide a counterexample to this interpretation. Despite conditions conducive to ecological speciation, our results suggest that dramatic size and environmental differentiation within C. antisiensis are not limiting gene flow.     

Local adaptation within a hybrid species

Ecological divergence among populations may be strongly influenced by their genetic background. For instance, genetic admixture through introgressive hybridization or hybrid speciation is likely to affect the genetic variation and evolvability of phenotypic traits. We studied geographic variation in two beak dimensions and three other phenotypic traits of the Italian sparrow (Passer italiae), a young hybrid species formed through interbreeding between house sparrows (P. domesticus) and Spanish sparrows (P. hispaniolensis). We found that beak morphology was strongly influenced by precipitation regimes and that it appeared to be the target of divergent selection within Italian sparrows. Interestingly, however, the degree of parental genetic contribution in the hybrid species had no effect on phenotypic beak variation. Moreover, beak height divergence may mediate genetic differentiation between populations, consistent with isolation-by-adaptation within this hybrid species. The study illustrates how hybrid species may be relatively unconstrained by their admixed genetic background, allowing them to adapt rapidly to environmental variation.     

Selection on outlier loci and their association with adaptive phenotypes in Littorina saxatilis contact zones

A fundamental issue in speciation research is to evaluate phenotypic variation and the genomics driving the evolution of reproductive isolation between sister taxa. Above all, hybrid zones are excellent study systems for researchers to examine the association of genetic differentiation, phenotypic variation and the strength of selection. We investigated two contact zones in the marine gastropod Littorina saxatilis and utilized landmark‐based geometric morphometric analysis together with amplified fragment length polymorphism (AFLP) markers to assess phenotypic and genomic divergence between ecotypes under divergent selection. From genetic markers, we calculated the cline width, linkage disequilibrium and the average effective selection on a locus. Additionally, we conducted an association analysis linking the outlier loci and phenotypic variation between ecotypes and show that a proportion of outlier loci are associated with key adaptive phenotypic traits.     

Effects of assortative mate choice on the genomic and morphological structure of a hybrid zone between two bird subspecies

Phenotypic differentiation plays an important role in the formation and maintenance of reproductive barriers. In some cases, variation in a few key aspects of phenotype can promote and maintain divergence; hence, the identification of these traits and their associations with patterns of genomic divergence is crucial for understanding the patterns and processes of population differentiation. We studied hybridization between the alba and personata subspecies of the white wagtail (Motacilla alba), and quantified divergence and introgression of multiple morphological traits and 19,437 SNP loci on a 3,000 km transect. Our goal was to identify traits that may contribute to reproductive barriers and to assess how variation in these traits corresponds to patterns of genome‐wide divergence. Variation in only one trait—head plumage patterning—was consistent with reproductive isolation. Transitions in head plumage were steep and occurred over otherwise morphologically and genetically homogeneous populations, whereas cline centres for other traits and genomic ancestry were displaced over 100 km from the head cline. Field observational data show that social pairs mated assortatively by head plumage, suggesting that these phenotypes are maintained by divergent mating preferences. In contrast, variation in all other traits and genetic markers could be explained by neutral diffusion, although weak ecological selection cannot be ruled out. Our results emphasize that assortative mating may maintain phenotypic differences independent of other processes shaping genome‐wide variation, consistent with other recent findings that raise questions about the relative importance of mate choice, ecological selection and selectively neutral processes for divergent evolution.     

Discordant patterns of introgression across a narrow hybrid zone between two cryptic lineages of an Iberian endemic newt

The study of natural hybrid zones can illuminate aspects of lineage divergence and speciation in morphologically cryptic taxa. We studied a hybrid zone between two highly divergent but morphologically similar lineages (south‐western and south‐eastern) of the Iberian endemic Bosca's newt (Lissotriton boscai) in SW Iberia with a multilocus dataset (microsatellites, nuclear and mitochondrial genes). STRUCTURE and NEWHYBRIDS analyses retrieved few admixed individuals, which classified as backcrosses involving parental individuals of the south‐western lineage. Our results show asymmetric introgression of mtDNA beyond the contact from this lineage into the south‐eastern lineage. Analysis of nongeographic introgression patterns revealed asymmetries in the direction of introgression, but except for mtDNA, we did not find evidence for nonconcordant introgression patterns across nuclear loci. Analysis of a 150‐km transect across the hybrid zone showed broadly coincident cline widths (ca. 3.2–27.9 km), and concordant cline centres across all markers, except for mtDNA that is displaced ca. 60 km northward. Results from ecological niche modelling show that the hybrid zone is in a climatically homogenous area with suitable habitat for the species, suggesting that contact between the two lineages is unlikely to occur further south as their distributions are currently separated by an extensive area of unfavourable habitat. Taken together, our findings suggest the genetic structure of this hybrid zone results from the interplay of historical (biogeographic) and population‐level processes. The narrowness and coincidence of genetic clines can be explained by weak selection against hybrids and reflect a degree of reproductive isolation that is consistent with cryptic speciation.     

Hybridization promotes color polymorphism in the aposematic harlequin poison frog,Oophaga histrionica

Whether hybridization can be a mechanism that drives phenotypic diversity is a widely debated topic in evolutionary biology. In poison frogs (Dendrobatidae), assortative mating has been invoked to explain how new color morphs persist despite the expected homogenizing effects of natural selection. Here, we tested the complementary hypothesis that new morphs arise through hybridization between different color morphs. Specifically, we (1) reconstructed the phylogenetic relationships among the studied populations of a dart‐poison frog to provide an evolutionary framework, (2) tested whether microsatellite allele frequencies of one putative hybrid population of the polymorphic frog O. histrionica are intermediate between O. histrionica and O. lehmanni, and (3) conducted mate‐choice experiments to test whether putatively intermediate females prefer homotypic males over males from the other two populations. Our findings are compatible with a hybrid origin for the new morph and emphasize the possibility of hybridization as a mechanism generating variation in polymorphic species. Moreover, because coloration in poison frogs is aposematic and should be heavily constrained, our findings suggest that hybridization can produce phenotypic novelty even in systems where phenotypes are subject to strong stabilizing selection.     

Structure and dynamics of hybrid zones at different stages of speciation in the common vole (Microtus arvalis)

The genetic structure and dynamics of hybrid zones provide crucial information for understanding the processes and mechanisms of evolutionary divergence and speciation. In general, higher levels of evolutionary divergence between taxa are more likely to be associated with reproductive isolation and may result in suppressed or strongly restricted hybridization. In this study, we examined two secondary contact zones between three deep evolutionary lineages in the common vole (Microtus arvalis). Differences in divergence times between the lineages can shed light on different stages of reproductive isolation and thus provide information on the ongoing speciation process in M. arvalis. We examined more than 800 individuals for mitochondrial (mtDNA), Y‐chromosome and autosomal markers and used assignment and cline analysis methods to characterize the extent and direction of gene flow in the contact zones. Introgression of both autosomal and mtDNA markers in a relatively broad area of admixture indicates selectively neutral hybridization between the least‐divergent lineages (Central and Eastern) without evidence for partial reproductive isolation. In contrast, a very narrow area of hybridization, shifts in marker clines and the quasi‐absence of Y‐chromosome introgression support a moving hybrid zone and unidirectional selection against male hybrids between the lineages with older divergence (Central and Western). Data from a replicate transect further support non‐neutral processes in this hybrid zone and also suggest a role for landscape history in the movement and shaping of geneflow profiles.     

Drift versus selection as drivers of phenotypic divergence at small spatial scales: The case of Belgjarskógur threespine stickleback

Divergence in phenotypic traits is facilitated by a combination of natural selection, phenotypic plasticity, gene flow, and genetic drift, whereby the role of drift is expected to be particularly important in small and isolated populations. Separating the components of phenotypic divergence is notoriously difficult, particularly for multivariate phenotypes. Here, we assessed phenotypic divergence of threespine stickleback (Gasterosteus aculeatus) across 19 semi‐interconnected ponds within a small geographic region (~7.5 km²) using comparisons of multivariate phenotypic divergence (PST), neutral genetic (FST), and environmental (EST) variation. We found phenotypic divergence across the ponds in a suite of functionally relevant phenotypic traits, including feeding, defense, and swimming traits, and body shape (geometric morphometric). Comparisons of PSTs with FSTs suggest that phenotypic divergence is predominantly driven by neutral processes or stabilizing selection, whereas phenotypic divergence in defensive traits is in accordance with divergent selection. Comparisons of population pairwise PSTs with ESTs suggest that phenotypic divergence in swimming traits is correlated with prey availability, whereas there were no clear associations between phenotypic divergence and environmental difference in the other phenotypic groups. Overall, our results suggest that phenotypic divergence of these small populations at small geographic scales is largely driven by neutral processes (gene flow, drift), although environmental determinants (natural selection or phenotypic plasticity) may play a role.     

Phenotypic and genetic introgression across a moving woodpecker hybrid zone

In hybrid zones in which two divergent taxa come into secondary contact and interbreed, selection can maintain phenotypic diversity despite widespread genetic introgression. Red‐breasted (Sphyrapicus ruber) and red‐naped (S. nuchalis) sapsuckers meet and hybridize along a narrow contact zone that stretches from northern California to southern British Columbia. We found strong evidence for changes in the structure of this hybrid zone across time, with significant temporal shifts in allele frequencies and in the proportions of parental phenotypes across the landscape. In addition to these shifts, we found that differences in plumage predict genetic differences (R² = 0.80), suggesting that plumage is a useful proxy for assessing ancestry. We also found a significant bimodal distribution of hybrids across the contact zone, suggesting that premating barriers may be driving reproductive isolation, perhaps as a result of assortative mating based on plumage differences. However, despite evidence of selection and strong patterns of population structure between parental samples, we found only weak patterns of genetic divergence. Using museum specimens and genomic data, this study of sapsuckers provides insight into the ways in which phenotypic and genetic structure have changed over a 40‐year period, as well as insight into the mechanisms that may contribute to the maintenance of the hybrid zone over time.     

Admixture mapping of male nuptial colour and body shape in a recently formed hybrid population of threespine stickleback

Despite recent progress, we still know relatively little about the genetic architecture that underlies adaptation to divergent environments. Determining whether the genetic architecture of phenotypic adaptation follows any predictable patterns requires data from a wide variety of species. However, in many organisms, genetic studies are hindered by the inability to perform genetic crosses in the laboratory or by long generation times. Admixture mapping is an approach that circumvents these issues by taking advantage of hybridization that occurs between populations or species in the wild. Here, we demonstrate the utility of admixture mapping in a naturally occurring hybrid population of threespine sticklebacks (Gasterosteus aculeatus) from Enos Lake, British Columbia. Until recently, this lake contained two species of sticklebacks adapted to divergent habitats within the lake. This benthic–limnetic species pair diverged in a number of phenotypes, including male nuptial coloration and body shape, which were previously shown to contribute to reproductive isolation between them. However, recent ecological disturbance has contributed to extensive hybridization between the species, and there is now a single, admixed population within Enos Lake. We collected over 500 males from Enos Lake and found that most had intermediate nuptial colour and body shape. By genotyping males with nuptial colour at the two extremes of the phenotypic distribution, we identified seven genomic regions on three chromosomes associated with divergence in male nuptial colour. These genomic regions are also associated with variation in body shape, suggesting that tight linkage and/or pleiotropy facilitated adaptation to divergent environments in benthic–limnetic species pairs.     

Hybridization between Townsend's Dendroica townsendi and black‐throated green warblers D. virens in an avian suture zone

Hybrid zones between species provide natural systems for the study of processes involved in divergence, reproductive isolation and speciation. Townsend's Dendroica townsendi and black‐throated green D. virens warblers are phenotypically and genetically divergent groups that occur in western and eastern North America respectively, with potential for range contact in the Rocky Mountains of British Columbia, where other west–east avian pairs come into contact. Although one potential hybrid (a phenotypic Townsend's warbler with the black‐throated green mitochondrial DNA) has been previously reported, there have been no studies of interactions between the taxa in potential areas of sympatry. To determine whether interbreeding between these species is a regular occurrence we examined variation in individuals across the area of putative range overlap. Analysis of plumage, morphology, and mitochondrial (COI) and nuclear molecular markers (CHD1Z and numt‐Dco1) shows surprisingly extensive hybridization between these species, with at least 38% of individuals in the hybrid zone being either hybrids or backcrosses. Each of the traits displays a sigmoidal cline centred along the eastern slope of the Rocky Mountains (molecular cline centres averaging 50 km east of the crest of the Rockies, ranging from 41 to 56 km). The clines are narrow (average molecular cline width is 60 km, ranging from 40 to 87 km) relative to the dispersal distance of related warbler species, suggesting that selection is maintaining the hybrid zone; we discuss possible sources of selection. Given the narrowness of the zone we recommend the two forms should continue to be treated as separate taxonomic species. Townsend's warblers also form an extensively studied hybrid zone with their more closely related southern relative, the hermit warbler D. occidentalis. The combined system of three discrete forms separated by narrow hybrid zones provides an excellent system for the study of hybridization, reproductive isolation and speciation.     

Parallel plumage colour evolution and introgressive hybridization in wheatears

Genetic and phenotypic mosaics, in which various phenotypes and different genomic regions show discordant patterns of species or population divergence, offer unique opportunities to study the role of ancestral and introgressed genetic variation in phenotypic evolution. Here, we investigated the evolution of discordant phenotypic and genetic divergence in a monophyletic clade of four songbird taxa—pied wheatear (O. pleschanka), Cyprus wheatear (Oenanthe cypriaca), and western and eastern subspecies of black‐eared wheatear (O. h. hispanica and O. h. melanoleuca). Phenotypically, black back and neck sides distinguish pied and Cyprus wheatears from the white‐backed/necked black‐eared wheatears. Meanwhile, mitochondrial variation only distinguishes western black‐eared wheatear. In the absence of nuclear genetic data, and given frequent hybridization among eastern black‐eared and pied wheatear, it remains unclear whether introgression is responsible for discordance between mitochondrial divergence patterns and phenotypic similarities, or whether plumage coloration evolved in parallel. Multispecies coalescent analyses of about 20,000 SNPs obtained from RAD data mapped to a draft genome assembly resolve the species tree, provide evidence for the parallel evolution of colour phenotypes and establish western and eastern black‐eared wheatears as independent taxa that should be recognized as full species. The presence of the entire admixture spectrum in the Iranian hybrid zone and the detection of footprints of introgression from pied into eastern black‐eared wheatear beyond the hybrid zone despite strong geographic structure of ancestry proportions furthermore suggest a potential role for introgression in parallel plumage colour evolution. Our results support the importance of standing heterospecific and/or ancestral variation in phenotypic evolution.     

Genetic and morphological divergence between Littorina fabalis ecotypes in Northern Europe

Low dispersal marine intertidal species facing strong divergent selective pressures associated with steep environmental gradients have a great potential to inform us about local adaptation and reproductive isolation. Among these, gastropods of the genus Littorina offer a unique system to study parallel phenotypic divergence resulting from adaptation to different habitats related with wave exposure. In this study, we focused on two Littorina fabalis ecotypes from Northern European shores and compared patterns of habitat‐related phenotypic and genetic divergence across three different geographic levels (local, regional and global). Geometric morphometric analyses revealed that individuals from habitats moderately exposed to waves usually present a larger shell size with a wider aperture than those from sheltered habitats. The phenotypic clustering of L. fabalis by habitat across most locations (mainly in terms of shell size) support an important role of ecology in morphological divergence. A genome scan based on amplified fragment length polymorphisms (AFLPs) revealed a heterogeneous pattern of differentiation across the genome between populations from the two different habitats, suggesting ecotype divergence in the presence of gene flow. The contrasting patterns of genetic structure between nonoutlier and outlier loci, and the decreased sharing of outlier loci with geographic distance among locations are compatible with parallel evolution of phenotypic divergence, with an important contribution of gene flow and/or ancestral variation. In the future, model‐based inference studies based on sequence data across the entire genome will help unravelling these evolutionary hypotheses, improving our knowledge about adaptation and its influence on diversification within the marine realm.     

The effect of hybrid transgression on environmental tolerance in experimental yeast crosses

Evidence is rapidly accumulating that hybridization generates adaptive variation. Transgressive segregation in hybrids could promote the colonization of new environments. Here, we use an assay to select hybrid genotypes that can proliferate in environmental conditions beyond the conditions tolerated by their parents, and we directly compete them against parental genotypes in habitats across environmental clines. We made 45 different hybrid swarms by crossing yeast strains (both Saccharomyces cerevisiae and S. paradoxus) with different genetic and phenotypic divergence. We compared the ability of hybrids and parents to colonize seven types of increasingly extreme environmental clines, representing both natural and novel challenges (mimicking pollution events). We found that a significant majority of hybrids had greater environmental ranges compared to the average of both their parents’ ranges (mid‐parent transgression), but only a minority of hybrids had ranges exceeding their best parent (best‐parent transgression). Transgression was affected by the specific strains involved in the cross and by the test environment. Genetic and phenotypic crossing distance predicted the extent of transgression in only two of the seven environments. We isolated a set of potentially transgressive hybrids selected at the extreme ends of the clines and found that many could directly outcompete their parents across whole clines and were between 1.5‐ and 3‐fold fitter on average. Saccharomyces yeast is a good model for quantitative and replicable experimental speciation studies, which may be useful in a world where hybridization is becoming increasingly common due to the relocation of plants and animals by humans.