The distribution of divergent mitochondrial DNA lineages of threespine stickleback (Gasterosteus aculeatus) in the northeastern Pacific Basin: post-glacial dispersal and lake accessibility

Aim This study furthers the documentation of the geographical distribution of two divergent (c. 3%) mitochondrial DNA clades in the threespine stickleback (Gasterosteus aculeatus) and tests the hypotheses that the northeastern Pacific distribution has been influenced by post‐glacial colonization and lake elevation and that clade identity is associated with certain morphological attributes such as reduction in body armour. Location Lakes and nearshore marine environments of the eastern Pacific Basin from southcentral Alaska to southeastern British Columbia, (BC) Canada. Methods Restriction enzyme analysis of polymerase chain reaction‐amplified mitochondrial DNA fragments (cytochrome b) from a total of 45 new populations combined with existing data for a further 45 populations. Lake elevation data were collected for 78 localities and tested for an association with mtDNA clade by contingency table analyses. Morphological data were collected on sticklebacks from eight samples representing four lake‐stream systems and tested for differentiation among populations with different mtDNA clade identities using analyses of variance. Results We extend the known distribution of the haplotypes diagnostic of the Trans‐North Pacific Clade (TNPC) southward to mid‐Vancouver Island and, for the first time, on mainland BC, in other island populations far from putative refugia, and in nearby anadromous populations. A morphological analysis indicated that the mainland population with the TNPC was not characterized by reduced spine or lateral plate (‘armour’) traits that characterize some putative relict populations on the Queen Charlotte Islands. We found a significant association between lake elevation and the presence of the TNPC; the TNPC was present more often in lakes located at or lower than 42 m than in higher elevation lakes. Main conclusions Our data support the hypothesis that post‐glacial colonization by TNPC‐bearing marine sticklebacks and aspects of lake ‘accessibility’ were important in determining the distribution of mtDNA clades in the eastern Pacific Ocean basin. More generally, our study demonstrates how processes acting both across immense geographic scales (e.g. pan‐Pacific dispersal) and local scales (lake accessibility contingent on timing and extent of isostatic rebound) may interact to explain biogeographical patterns.     

Divergent Macroparasite Infections in Parapatric Swiss Lake-Stream Pairs of Threespine Stickleback (Gasterosteus aculeatus)

Spatial heterogeneity in diversity and intensity of parasitism is a typical feature of most host-parasite interactions, but understanding of the evolutionary implications of such variation is limited. One possible outcome of infection heterogeneities is parasite-mediated divergent selection between host populations, ecotypes or species which may facilitate the process of ecological speciation. However, very few studies have described infections in population-pairs along the speciation continuum from low to moderate or high degree of genetic differentiation that would address the possibility of parasite-mediated divergent selection in the early stages of the speciation process. Here we provide an example of divergent parasitism in freshwater fish ecotypes by examining macroparasite infections in threespine stickleback (Gasterosteus aculeatus) of four Swiss lake systems each harbouring parapatric lake-stream ecotype pairs. We demonstrate significant differences in infections within and between the pairs that are driven particularly by the parasite taxa transmitted to fish from benthic invertebrates. The magnitude of the differences tended to correlate positively with the extent of neutral genetic differentiation between the parapatric lake and stream populations of stickleback, whereas no such correlation was found among allopatric populations from similar or contrasting habitats. This suggests that genetic differentiation is unrelated to the magnitude of parasite infection contrasts when gene flow is constrained by geographical barriers while in the absence of physical barriers, genetic differentiation and the magnitude of differences in infections tend to be positively correlated.     

Differences in rheotactic responses contribute to divergent habitat use between parapatric lake and stream threespine stickleback

Migration among populations is widely thought to undermine adaptive divergence, assuming gene flow arises from random movement of individuals. If individuals instead differ in dispersal behavior, phenotype‐dependent dispersal can reduce the effective rate of gene flow or even facilitate divergence. For example, parapatric populations of lake and stream stickleback tend to actively avoid dispersing into the adjoining habitat. However, the behavioral basis of this nonrandom dispersal was previously unknown. Here, we show that lake and stream stickleback exhibit divergent rheotactic responses (behavioral response to currents). During the breeding season, wild‐caught inlet stream stickleback were better than lake fish at maintaining position in currents, faced upstream more, and spent more time in low‐current areas. As a result, stream fish expended significantly less energy in currents than did lake fish. These divergent rheotactic responses likely contribute to divergent habitat use by lake and stream stickleback. Although rheotactic differences were absent in nonbreeding fish, divergent behavior of breeding‐season fish may suffice for assortative mating by breeding location. The resulting reproductive isolation between lake and stream fish may explain the fine‐scale evolutionary differentiation in parapatric stickleback populations.     

PHENOTYPE-DEPENDENT NATIVE HABITAT PREFERENCE FACILITATES DIVERGENCE BETWEEN PARAPATRIC LAKE AND STREAM STICKLEBACK

Adaptive divergence between adjoining populations reflects a balance between the diversifying effect of divergent selection and the potentially homogenizing effect of gene flow. In most models of migration-selection balance, gene flow is assumed to reflect individuals' inherent capacity to disperse, without regard to the match between individuals' phenotypes and the available habitats. However, habitat preferences can reduce dispersal between contrasting habitats, thereby alleviating migration load and facilitating adaptive divergence. We tested whether habitat preferences contribute to adaptive divergence in a classic example of migration-selection balance: parapatric lake and stream populations of three-spine stickleback ( Gasterosteus aculeatus ). Using a mark-transplant-recapture experiment on morphologically divergent parapatric populations, we showed that 90% of lake and stream stickleback returned to their native habitat, reducing migration between habitats by 76%. Furthermore, we found that dispersal into a nonnative habitat was phenotype dependent. Stream fish moving into the lake were morphologically more lake-like than those returning to the stream (and the converse for lake fish entering the stream). The strong native habitat preference documented here increases the extent of adaptive divergence between populations two- to fivefold relative to expectations with random movement. These results illustrate the potential importance of adaptive habitat choice in driving parapatric divergence.     

Quantifying the constraining influence of gene flow on adaptive divergence in the lake-stream threespine stickleback system

The constraining effect of gene flow on adaptive divergence is often inferred but rarely quantified. We illustrate ways of doing so using stream populations of threespine stickleback (Gasterosteus aculeatus) that experience different levels of gene flow from a parapatric lake population. In the Misty Lake watershed (British Columbia, Canada), the inlet stream population is morphologically divergent from the lake population, and presumably experiences little gene flow from the lake. The outlet stream population, however, shows an intermediate phenotype and may experience more gene flow from the lake. We first used microsatellite data to demonstrate that gene flow from the lake is low into the inlet but high into the outlet, and that gene flow from the lake remains relatively constant with distance along the outlet. We next combined gene flow data with morphological and habitat data to quantify the effect of gene flow on morphological divergence. In one approach, we assumed that inlet stickleback manifest well-adapted phenotypic trait values not constrained by gene flow. We then calculated the deviation between the observed and expected phenotypes for a given habitat in the outlet. In a second approach, we parameterized a quantitative genetic model of adaptive divergence. Both approaches suggest a large impact of gene flow, constraining adaptation by 80-86% in the outlet (i.e., only 14-20% of the expected morphological divergence in the absence of gene flow was observed). Such approaches may be useful in other taxa to estimate how important gene flow is in constraining adaptive divergence in nature.     

Along the speciation continuum in sticklebacks

Speciation can be viewed as a continuum, potentially divisible into several states: (1) continuous variation within panmictic populations, (2) partially discontinuous variation with minor reproductive isolation, (3) strongly discontinuous variation with strong but reversible reproductive isolation and (4) complete and irreversible reproductive isolation. Research on sticklebacks (Gasterosteidae) reveals factors that influence progress back and forth along this continuum, as well as transitions between the states. Most populations exist in state 1, even though some of these show evidence of disruptive selection and positive assortative mating. Transitions to state 2 seem to usually involve strong divergent selection coupled with at least a bit of geographic separation, such as parapatry (e.g. lake and stream pairs and mud and lava pairs) or allopatry (e.g. different lakes). Transitions to state 3 can occur when allopatric or parapatric populations that evolved under strong divergent selection come into secondary contact (most obviously the sympatric benthic and limnetic pairs), but might also occur between populations that remained in parapatry or allopatry. Transitions to state 4 might be decoupled from these selective processes, because the known situations of complete, or nearly complete, reproductive isolation (Japan Sea and Pacific Ocean pair and the recognized gasterosteid species) are always associated with chromosomal rearrangements and environment‐independent genetic incompatibilities. Research on sticklebacks has thus revealed complex and shifting interactions between selection, adaptation, mutation and geography during the course of speciation.     

Asymmetric reproductive barriers and mosaic reproductive isolation: insights from Misty lake–stream stickleback

Ecological speciation seems to occur readily but is clearly not ubiquitous – and the relative contributions of different reproductive barriers remain unclear in most systems. We here investigate the potential importance of selection against migrants in lake/stream stickleback (Gasterosteus aculeatus) from the Misty Lake system, Canada. This system is of particular interest because one population contrast (Lake vs. Outlet stream) shows very low genetic and morphological divergence, whereas another population contrast (Lake vs. Inlet stream) shows dramatic genetic and morphological divergence apparently without strong and symmetric reproductive barriers. To test whether selection against migrants might solve this “conundrum of missing reproductive isolation”, we performed a fully factorial reciprocal transplant experiment using 225 individually marked stickleback collected from the wild. Relative fitness of the different ecotypes (Lake, Inlet, and Outlet) was assessed based on survival and mass change in experimental enclosures. We found that Inlet fish performed poorly in the lake (selection against migrants in that direction), whereas Lake fish outperformed Inlet fish in all environments (no selection against migrants in the opposite direction). As predicted from their phenotypic and genetic similarity, Outlet and Lake fish performed similarly in all environments. These results suggest that selection against migrants is asymmetric and, together with previous work, indicates that multiple reproductive barriers contribute to reproductive isolation. Similar mosaic patterns of reproductive isolation are likely in other natural systems.     

Repeated Lake-Stream Divergence in Stickleback Life History within a Central European Lake Basin

Life history divergence between populations inhabiting ecologically distinct habitats might be a potent source of reproductive isolation, but has received little attention in the context of speciation. We here test for life history divergence between threespine stickleback inhabiting Lake Constance (Central Europe) and multiple tributary streams. Otolith analysis shows that lake fish generally reproduce at two years of age, while their conspecifics in all streams have shifted to a primarily annual life cycle. This divergence is paralleled by a striking and consistent reduction in body size and fecundity in stream fish relative to lake fish. Stomach content analysis suggests that life history divergence might reflect a genetic or plastic response to pelagic versus benthic foraging modes in the lake and the streams. Microsatellite and mitochondrial markers further reveal that life history shifts in the different streams have occurred independently following the colonization by Lake Constance stickleback, and indicate the presence of strong barriers to gene flow across at least some of the lake-stream habitat transitions. Given that body size is known to strongly influence stickleback mating behavior, these barriers might well be related to life history divergence.     

Morphological divergence of North‐European nine‐spined sticklebacks (Pungitius pungitius): signatures of parallel evolution

Parallel evolution is characterised by repeated, independent occurrences of similar phenotypes in a given habitat type, in different parts of the species distribution area. We studied body shape and body armour divergence between five marine, four lake, and ten pond populations of nine‐spined sticklebacks [Pungitius pungitius (Linnaeus, 1758)] in Fennoscandia. We hypothesized that marine and lake populations (large water bodies, diverse fish fauna) would be similar, whereas sticklebacks in isolated ponds (small water bodies, simple fish fauna) would be divergent. We found that pond fish had deeper bodies, shorter caudal peduncles, and less body armour (viz. shorter/absent pelvic spines, reduced/absent pelvic girdle, and reduced number of lateral plates) than marine fish. Lake fish were intermediate, but more similar to marine than to pond fish. Results of our common garden experiment concurred with these patterns, suggesting a genetic basis for the observed divergence. We also found large variation among populations within habitat types, indicating that environmental variables other than those related to gross habitat characteristics might also influence nine‐spined stickleback morphology. Apart from suggesting parallel evolution of morphological characteristics of nine‐spined sticklebacks in different habitats, the results also show a number of similarities to the evolution of three‐spined stickleback (Gasterosteus aculeatus Linnaeus, 1758) morphology. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 101 , 403–416.     

Genome divergence during evolutionary diversification as revealed in replicate lake-stream stickleback population pairs

Evolutionary diversification is often initiated by adaptive divergence between populations occupying ecologically distinct environments while still exchanging genes. The genetic foundations of this divergence process are largely unknown and are here explored through genome scans in multiple independent lake-stream population pairs of threespine stickleback. We find that across the pairs, overall genomic divergence is associated with the magnitude of divergence in phenotypes known to be under divergent selection. Along this same axis of increasing diversification, genomic divergence becomes increasingly biased towards the centre of chromosomes as opposed to the peripheries. We explain this pattern by within-chromosome variation in the physical extent of hitchhiking, as recombination is greatly reduced in chromosome centres. Correcting for this effect suggests that a great number of genes distributed widely across the genome are involved in the divergence into lake vs. stream habitats. Analyzing additional allopatric population pairs, however, reveals that strong divergence in some genomic regions has been driven by selection unrelated to lake-stream ecology. Our study highlights a major contribution of large-scale variation in recombination rate to generating heterogeneous genomic divergence and indicates that elucidating the genetic basis of adaptive divergence might be more challenging than currently recognized.     

Evolutionary history of an adaptive radiation in species pairs of threespine sticklebacks (Gasterosteus): insights from mitochondrial DNA

Species pairs of threespine stickleback, Gasterosteus aculeatus, co-exist in several lakes in the Strait of Georgia, southwestern British Columbia. One species, ‘benthics’ is robust-bodied and is morphologically and behaviourally specialized for benthivory. The other species, ‘limnetics’ is specialized for planktivory in open-water habitats of the lakes. We examined mitochondrial DNA restriction site variation in benthic and limnetic sticklebacks as well as in solitary freshwater, anadromous (sea-run), and marine populations to test: (i) if benthic and limnetic pairs have evolved only once or multiple times (parallel evolution) and (ii) if the species have evolved sympatrically, or allopatrically from ‘double invasions’ of lakes by ancestral anadromous/marine sticklebacks. Stickleback mtDNA comprised a single clade with a low (mean = 0.40%) degree of sequence divergence among the 77 haplotypes resolved. Most nucleotide diversity (97%) was found within (rather than among) populations of anadromous/marine sticklebacks whereas most diversity (77%) was found among populations in freshwater sticklebacks. Significant differences in haplotype frequencies were found between benthics and limnetics in three of the four species pair lakes examined, but in all cases the pairs within lakes were characterized by unique assemblages of closely related haplotypes. Hierarchical clustering of divergence estimates suggested that comparable species from different lakes have originated independently in all lakes because in no case did comparable species from different lakes cluster together. Divergent species within lakes tended to be more closely related to one another than to species in other lakes and there were two cases were benthics and limnetics within a particular lake were monophyletic. In two of the four two-species lakes, limnetics were less divergent from putative ancestral anadromous/marine stickleback as predicted by the double invasion hypothesis, but in the two other lakes benthics were less divergent. Our data argue strongly that the species pairs have evolved independently in each lake were they now co-exist. Further, in two lakes our data are consistent with the species having evolved by sympatric divergence, but allopatric divergence followed by introgression of mtDNA that has obscured ancestral relationships cannot be discounted completely. Finally, despite remaining uncertainty about the geography of speciation, the species appear to have evolved in the face of gene flow arguing that natural selection acting on trophic ecology has been a major component of ecological speciation in sticklebacks.     

Demography and genome divergence of lake and stream populations of an East African cichlid fish

Disentangling the processes and mechanisms underlying adaptive diversification is facilitated by the comparative study of replicate population pairs that have diverged along a similar environmental gradient. Such a setting is realized in a cichlid fish from southern Lake Tanganyika, Astatotilapia burtoni, which occurs within the lake proper as well as in various affluent rivers. Previously, we demonstrated that independent lake and stream populations show similar adaptations to the two habitat regimes. However, little is known about the evolutionary and demographic history of the A. burtoni populations in question and the patterns of genome divergence among them. Here, we apply restriction site‐associated DNA sequencing (RADseq) to examine the evolutionary history, the population structure and genomic differentiation of lake and stream populations in A. burtoni. A phylogenetic reconstruction based on genome‐wide molecular data largely resolved the evolutionary relationships among populations, allowing us to re‐evaluate the independence of replicate lake–stream population clusters. Further, we detected a strong pattern of isolation by distance, with baseline genomic divergence increasing with geographic distance and decreasing with the level of gene flow between lake and stream populations. Genome divergence patterns were heterogeneous and inconsistent among lake‐stream population clusters, which is explained by differences in divergence times, levels of gene flow and local selection regimes. In line with the latter, we only detected consistent outlier loci when the most divergent lake–stream population pair was excluded. Several of the thus identified candidate genes have inferred functions in immune and neuronal systems and show differences in gene expression between lake and stream populations.     

Adaptive divergence and the balance between selection and gene flow: lake and stream stickleback in the Misty system

We investigated the interplay between natural selection and gene flow in the adaptive divergence of threespine stickleback (Gasterosteus aculeatus) that reside parapatrically in lakes and streams. Within the Misty Lake system (Vancouver Island, British Columbia), stickleback from the inlet stream (flowing into the lake) have fewer gill rakers and deeper bodies than stickleback from the lake--differences thought to facilitate foraging (benthic macroinvertebrates in the stream vs. zooplankton in the open water of the lake). Common-garden experiments demonstrated that these differences have a genetic basis. Reciprocal transplant enclosure experiments showed that lake and inlet stickleback grow best in their home environments (although differences were subtle and often not significant). Release-recapture experiments in the inlet showed that lake fish are less well-suited than inlet fish for life in the stream (higher mortality or emigration in lake fish). Morphological divergence in the wild and under common rearing was greater between the lake and the inlet than between the lake and the outlet. Genetic divergence (mitochondrial DNA and microsatellites) was greatest between the lake and the upper inlet (1.8 km upstream from the lake), intermediate between the lake and the lower inlet (0.9 km upstream), and least between the lake and the outlet stream (1.2 km downstream). Relative levels of gene flow estimated from genetic data showed the inverse pattern. The negative association between morphological divergence and gene flow is consistent with the expectation that gene flow can constrain adaptation. Estimated absolute levels of gene flow also implied a constraint on adaptation in the outlet but not the inlet. Our results suggest that natural selection promotes the adaptive divergence of lake and stream stickleback. but that the magnitude of divergence can be constrained by gene flow.     

Trans-Arctic dispersals and the evolution of a circumpolar marine fish species complex, the capelin (Mallotus villosus)

Trans-Arctic dispersals and population and range expansions during the Pleistocene enhanced opportunities for evolutionary diversification and contributed to the process of speciation within the capelin, a northern marine-fish complex exhibiting a circumpolar distribution. Capelin is composed of four highly divergent and geographically discrete mitochondrial DNA (mtDNA) clades (609 bp; cytochrome b). Two clades occur in the North Atlantic, one associated with Canadian Atlantic waters, including Hudson Bay, and the second distributed from West Greenland to the Barents Sea. Two additional clades occur in the Arctic and northeast Pacific Oceans, representing the most recent divergence within the capelin phylogenetic tree. Judged from mtDNA diversity, capelin populations comprising all clades experienced at least one demographic and spatial reduction-expansion episode during recent Pleistocene glaciations that imprinted their molecular architecture. The large contemporary populations in the northeast Pacific and Arctic Oceans exhibited significant genetic structure whereas no such structure was detected in the equally extensive North Atlantic clades. All clades are characterized by one or two prevalent mtDNA haplotypes distributed over the entire range of the clade. Assuming a Pacific ancestor for capelin, we infer that capelin dispersed on two separate occasions to the North Atlantic. A more recent event resulted in the isolation of eastern Pacific and Arctic clades, with the Arctic clade positioned for a potential third Atlantic invasion, as revealed by the presence of this clade in the Labrador Sea. The Labrador Sea is a potential contact zone for three of the four capelin clades.     

Partitioning the effects of isolation by distance,environment, and physical barriers on genomic divergence between parapatric threespine stickleback

Genetic divergence between populations is shaped by a combination of drift, migration, and selection, yielding patterns of isolation‐by‐distance (IBD) and isolation‐by‐environment (IBE). Unfortunately, IBD and IBE may be confounded when comparing divergence across habitat boundaries. For instance, parapatric lake and stream threespine stickleback (Gasterosteus aculeatus) may have diverged due to selection against migrants (IBE), or mere spatial separation (IBD). To quantitatively partition the strength of IBE and IBD, we used recently developed population genetic software (BEDASSLE) to analyze partial genomic data from three lake‐stream clines on Vancouver Island. We find support for IBD within each of three outlet streams (unlike prior studies of lake‐stream stickleback). In addition, we find evidence for IBE (controlling for geographic distance): the genetic effect of habitat is equivalent to geographic separation of ～1.9 km of IBD. Remarkably, of our three lake‐stream pairs, IBE is strongest where migration between habitats is easiest. Such microgeographic genetic divergence would require exceptionally strong divergent selection, which multiple experiments have failed to detect. Instead, we propose that nonrandom dispersal (e.g., habitat choice) contributes to IBE. Supporting this conclusion, we show that the few migrants between habitats are a nonrandom subset of the phenotype distribution of the source population.     

Adaptive divergence between lake and stream populations of an East African cichlid fish

Divergent natural selection acting in different habitats may build up barriers to gene flow and initiate speciation. This speciation continuum can range from weak or no divergence to strong genetic differentiation between populations. Here, we focus on the early phases of adaptive divergence in the East African cichlid fish Astatotilapia burtoni, which occurs in both Lake Tanganyika (LT) and inflowing rivers. We first assessed the population structure and morphological differences in A. burtoni from southern LT. We then focused on four lake–stream systems and quantified body shape, ecologically relevant traits (gill raker and lower pharyngeal jaw) as well as stomach contents. Our study revealed the presence of several divergent lake–stream populations that rest at different stages of the speciation continuum, but show the same morphological and ecological trajectories along the lake–stream gradient. Lake fish have higher bodies, a more superior mouth position, longer gill rakers and more slender pharyngeal jaws, and they show a plant/algae and zooplankton‐biased diet, whereas stream fish feed more on snails, insects and plant seeds. A test for reproductive isolation between closely related lake and stream populations did not detect population‐assortative mating. Analyses of F1 offspring reared under common garden conditions indicate that the detected differences in body shape and gill raker length do not constitute pure plastic responses to different environmental conditions, but also have a genetic basis. Taken together, the A. burtoni lake–stream system constitutes a new model to study the factors that enhance and constrain progress towards speciation in cichlid fishes.     

VARIABLE PROGRESS TOWARD ECOLOGICAL SPECIATION IN PARAPATRY: STICKLEBACK ACROSS EIGHT LAKE-STREAM TRANSITIONS

Divergent selection between contrasting habitats can sometimes drive adaptive divergence and the evolution of reproductive isolation in the face of initially high gene flow. "Progress" along this ecological speciation pathway can range from minimal divergence to full speciation. We examine this variation for threespine stickleback fish that evolved independently across eight lake-stream habitat transitions. By quantifying stickleback diets, we show that lake-stream transitions usually coincide with limnetic-benthic ecotones. By measuring genetically based phenotypes, we show that these ecotones often generate adaptive divergence in foraging morphology. By analyzing neutral genetic markers (microsatellites), we show that adaptive divergence is often associated with the presence of two populations maintaining at least partial reproductive isolation in parapatry. Coalescent-based simulations further suggest that these populations have diverged with gene flow within a few thousand generations, although we cannot rule out the possibility of phases of allopatric divergence. Finally, we find striking variation among the eight lake-stream transitions in progress toward ecological speciation. This variation allows us to hypothesize that progress is generally promoted by strong divergent selection and limited dispersal across the habitat transitions. Our study thus makes a case for ecological speciation in a parapatric context, while also highlighting variation in the outcome.     

Parallel and nonparallel ecological,morphological and genetic divergence in lake–stream stickleback from a single catchment

Parallel phenotypic evolution in similar environments has been well studied in evolutionary biology; however, comparatively little is known about the influence of determinism and historical contingency on the nature, extent and generality of this divergence. Taking advantage of a novel system containing multiple lake–stream stickleback populations, we examined the extent of ecological, morphological and genetic divergence between three‐spined stickleback present in parapatric environments. Consistent with other lake–stream studies, we found a shift towards a deeper body and shorter gill rakers in stream fish. Morphological shifts were concurrent with changes in diet, indicated by both stable isotope and stomach contents analysis. Performing a multivariate test for shared and unique components of evolutionary response to the distance gradient from the lake, we found a strong signature of parallel adaptation. Nonparallel divergence was also present, attributable mainly to differences between river locations. We additionally found evidence of genetic substructuring across five lake–stream transitions, indicating that some level of reproductive isolation occurs between populations in these habitats. Strong correlations between pairwise measures of morphological, ecological and genetic distance between lake and stream populations supports the hypothesis that divergent natural selection between habitats drives adaptive divergence and reproductive isolation. Lake–stream stickleback divergence in Lough Neagh provides evidence for the deterministic role of selection and supports the hypothesis that parallel selection in similar environments may initiate parallel speciation.     

Parallel and nonparallel aspects of ecological, phenotypic, and genetic divergence across replicate population pairs of lake and stream stickleback

Parallel (or convergent) evolution provides strong evidence for a deterministic role of natural selection: similar phenotypes evolve when independent populations colonize similar environments. In reality, however, independent populations in similar environments always show some differences: some nonparallel evolution is present. It is therefore important to explicitly quantify the parallel and nonparallel aspects of trait variation, and to investigate the ecological and genetic explanations for each. We performed such an analysis for threespine stickleback (Gasterosteus aculeatus) populations inhabiting lake and stream habitats in six independent watersheds. Morphological traits differed in the degree to which lake-stream divergence was parallel across watersheds. Some aspects of this variation were correlated with ecological variables related to diet, presumably reflecting the strength and specifics of divergent selection. Furthermore, a genetic scan revealed some markers that diverged between lakes and streams in many of the watersheds and some that diverged in only a few watersheds. Moreover, some of the lake-stream divergence in genetic markers was associated within some of the lake-stream divergence in morphological traits. Our results suggest that parallel evolution, and deviations from it, are primarily the result of natural selection, which corresponds in only some respects to the dichotomous habitat classifications frequently used in such studies.     

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.