Contrasting Patterns of Genetic Differentiation among Blackcaps (Sylvia atricapilla) with Divergent Migratory Orientations in Europe

Migratory divides are thought to facilitate behavioral, ecological, and genetic divergence among populations with different migratory routes. However, it is currently contentious how much genetic divergence is needed to maintain distinct migratory behavior across migratory divides. Here we investigate patterns of neutral genetic differentiation among Blackcap (Sylvia atricapilla) populations with different migratory strategies across Europe. We compare the level of genetic divergence of populations migrating to southwestern (SW) or southeastern (SE) wintering areas with birds wintering in the British Isles following a recently established northwesterly (NW) migration route. The migratory divide between SW and SE wintering areas can be interpreted as a result of a re-colonization process after the last glaciation. Thus we predicted greater levels of genetic differentiation among the SW/SE populations. However, a lack of genetic differentiation was found between SW and SE populations, suggesting that interbreeding likely occurs among Blackcaps with different migratory orientations across a large area; therefore the SW/SE migratory divide can be seen as diffuse, broad band and is, at best, a weak isolating barrier. Conversely, weak, albeit significant genetic differentiation was evident between NW and SW migrants breeding sympatrically in southern Germany, suggesting a stronger isolating mechanism may be acting in this population. Populations located within/near the SW/SE contact zone were the least genetically divergent from NW migrants, confirming NW migrants likely originated from within the contact zone. Significant isolation-by-distance was found among eastern Blackcap populations (i.e. SE migrants), but not among western populations (i.e. NW and SW migrants), revealing different patterns of genetic divergence among Blackcap populations in Europe. We discuss possible explanations for the genetic structure of European Blackcaps and how gene flow influences the persistence of divergent migratory behaviors.     

Integrative tracking methods elucidate the evolutionary dynamics of a migratory divide

Migratory divides, the boundary between adjacent bird populations that migrate in different directions, are of considerable interest to evolutionary biologists because of their alleged role in speciation of migratory birds. However, the small size of many passerines has traditionally limited the tools available to track populations and as a result, restricted our ability to study how reproductive isolation might occur across a divide. Here, we integrate multiple approaches by using genetic, geolocator, and morphological data to investigate a migratory divide in hermit thrushes (Catharus guttatus). First, high genetic divergence between migratory groups indicates the divide is a region of secondary contact between historically isolated populations. Second, despite low sample sizes, geolocators reveal dramatic differences in overwintering locations and migratory distance of individuals from either side of the divide. Third, a diagnostic genetic marker that proved useful for tracking a key population suggests a likely intermediate nonbreeding location of birds from the hybrid zone. This finding, combined with lower return rates from this region, is consistent with comparatively lower fitness of hybrids, which is possibly due to this intermediate migration pattern. We discuss our results in the context of reproductive isolating mechanisms associated with migration patterns that have long been hypothesized to promote divergence across migratory divides.     

Local adaptation does not lead to genome‐wide differentiation in lava flow lizards

Adaptation can occur with or without genome‐wide differentiation. If adaptive loci are linked to traits involved in reproductive isolation, genome‐wide divergence is likely, and speciation is possible. However, adaptation can also lead to phenotypic differentiation without genome‐wide divergence if levels of ongoing gene flow are high. Here, we use the replicated occurrence of melanism in lava flow lizards to assess the relationship between local adaptation and genome‐wide differentiation. We compare patterns of phenotypic and genomic divergence among lava flow and nonlava populations for three lizard species and three lava flows in the Chihuahuan Desert. We find that local phenotypic adaptation (melanism) is not typically accompanied by genome‐wide differentiation. Specifically, lava populations do not generally exhibit greater divergence from nonlava populations than expected by geography alone, regardless of whether the lava formation is 5,000 or 760,000 years old. We also infer that gene flow between lava and nonlava populations is ongoing in all lava populations surveyed. Recent work in the isolation by environment and ecological speciation literature suggests that environmentally driven genome‐wide differentiation is common in nature. However, local adaptation may often simply be local adaptation rather than an early stage of ecological speciation.     

Ecological niche dimensionality and the evolutionary diversification of stick insects

The degree of phenotypic divergence and reproductive isolation between taxon pairs can vary quantitatively, and often increases as evolutionary divergence proceeds through various stages, from polymorphism to population differentiation, ecotype and race formation, speciation, and post-speciational divergence. Although divergent natural selection promotes divergence, it does not always result in strong differentiation. For example, divergent selection can fail to complete speciation, and distinct species pairs sometimes collapse ('speciation in reverse'). Widely-discussed explanations for this variability concern genetic architecture, and the geographic arrangement of populations. A less-explored possibility is that the degree of phenotypic and reproductive divergence between taxon pairs is positively related to the number of ecological niche dimensions (i.e., traits) subject to divergent selection. Some data supporting this idea stem from laboratory experimental evolution studies using Drosophila, but tests from nature are lacking. Here we report results from manipulative field experiments in natural populations of herbivorous Timema stick insects that are consistent with this 'niche dimensionality' hypothesis. In such insects, divergent selection between host plants might occur for cryptic colouration (camouflage to evade visual predation), physiology (to detoxify plant chemicals), or both of these niche dimensions. We show that divergent selection on the single niche dimension of cryptic colouration can result in ecotype formation and intermediate levels of phenotypic and reproductive divergence between populations feeding on different hosts. However, greater divergence between a species pair involved divergent selection on both niche dimensions. Although further replication of the trends reported here is required, the results suggest that dimensionality of selection may complement genetic and geographic explanations for the degree of diversification in nature.     

Phenotypic divergence and reproductive isolation between sympatric forms of Japanese threespine sticklebacks

The threespine stickleback ( Gasterosteus aculeatus ) species complex is well suited for identifying the types of phenotypic divergence and isolating barriers that contribute to reproductive isolation at early stages of speciation. In the present study, we characterize the patterns of genetic and phenotypic divergence as well as the types of isolating barriers that are present between two sympatric pairs of threespine sticklebacks in Hokkaido, Japan. One sympatric pair consists of an anadromous and a resident freshwater form and shows divergence in body size between the forms, despite the lack of genetic differentiation between them. The second sympatric pair consists of two anadromous forms, which originated before the last glacial period and are currently reproductively isolated. These two anadromous forms have diverged in many morphological traits as well as in their reproductive behaviours. Both sexual isolation and hybrid male sterility contribute to reproductive isolation between the anadromous species pair. We discuss the shared and unique aspects of phenotypic divergence and reproductive isolation in the Japanese sympatric pairs compared with postglacial stickleback species pairs. Further studies of these divergent species pairs will provide a deeper understanding of the mechanisms of speciation in sticklebacks.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 671–685.     

Migratory divides coincide with reproductive barriers across replicated avian hybrid zones above the Tibetan Plateau

Migratory divides are proposed to be catalysts for speciation across a diversity of taxa. However, it is difficult to test the relative contributions of migratory behaviour vs. other divergent traits to reproductive isolation. Comparing hybrid zones with and without migratory divides offers a rare opportunity to directly examine the contribution of divergent migratory behaviour to reproductive barriers. We show that across replicate sampling transects of two pairs of barn swallow (Hirundo rustica) subspecies, strong reproductive isolation coincided with a migratory divide spanning 20 degrees of latitude. A third subspecies pair exhibited no evidence for a migratory divide and hybridised extensively. Within migratory divides, overwintering habitats were associated with assortative mating, implicating a central contribution of divergent migratory behaviour to reproductive barriers. The remarkable geographic coincidence between migratory divides and genetic breaks supports a long‐standing hypothesis that the Tibetan Plateau is a substantial barrier contributing to the diversity of Siberian avifauna.     

Phenotypic differentiation and pre-mating isolation between allopatric populations of Girardinichthys multiradiatus

Sexual selection may lead to reproductive isolation between populations through divergence in female mate choice, and population differentiation driven by female mate choice is expected to produce pre- but not post-mating isolation. We tested these hypotheses by looking at whether allopatric populations of the Amarillo (Girardinichthys multiradiatus), a sexually dimorphic viviparous fish with effective female choice, (i) have undergone phenotypical differentiation that may be attributed to divergence in female mate choice, and (ii) are already separated by pre- and/or post-mating reproductive barriers. We found substantial divergence in morphological traits which are the target of female mate choice, and in male courtship behaviour. Strong female preferences for homogametic males indicate substantial and symmetric pre-copulatory isolation, but the few successful heterogametic crosses produced in confinement yielded litters of the same size as those produced in homogametic matings, suggesting that post-copulatory isolation between populations is non-existent or weak. It appears that the studied populations have undergone incipient speciation with a pattern that is consistent with speciation driven by sexual selection, yet further work should assess whether divergence in female preferences has promoted male phenotypic differentiation or whether variation in male attributes has driven divergence in female preferences.     

Genome‐wide differentiation in closely related populations: the roles of selection and geographic isolation

Population divergence in geographic isolation is due to a combination of factors. Natural and sexual selection may be important in shaping patterns of population differentiation, a pattern referred to as ‘isolation by adaptation’ (IBA). IBA can be complementary to the well‐known pattern of ‘isolation by distance’ (IBD), in which the divergence of closely related populations (via any evolutionary process) is associated with geographic isolation. The barn swallow Hirundo rustica complex comprises six closely related subspecies, where divergent sexual selection is associated with phenotypic differentiation among allopatric populations. To investigate the relative contributions of selection and geographic distance to genome‐wide differentiation, we compared genotypic and phenotypic variation from 350 barn swallows sampled across eight populations (28 pairwise comparisons) from four different subspecies. We report a draft whole‐genome sequence for H. rustica, to which we aligned a set of 9493 single nucleotide polymorphisms (SNPs). Using statistical approaches to control for spatial autocorrelation of phenotypic variables and geographic distance, we find that divergence in traits related to migratory behaviour and sexual signalling, as well as geographic distance, together explain over 70% of genome‐wide divergence among populations. Controlling for IBD, we find 42% of genomewide divergence is attributable to IBA through pairwise differences in traits related to migratory behaviour and sexual signalling alone. By (i) combining these results with prior studies of how selection shapes morphological differentiation and (ii) accounting for spatial autocorrelation, we infer that morphological adaptation plays a large role in shaping population‐level differentiation in this group of closely related populations.     

Low reproductive isolation and highly variable levels of gene flow reveal limited progress towards speciation between European river and brook lampreys

Ecologically based divergent selection is a factor that could drive reproductive isolation even in the presence of gene flow. Population pairs arrayed along a continuum of divergence provide a good opportunity to address this issue. Here, we used a combination of mating trials, experimental crosses and population genetic analyses to investigate the evolution of reproductive isolation between two closely related species of lampreys with distinct life histories. We used microsatellite markers to genotype over 1000 individuals of the migratory parasitic river lamprey (Lampetra fluviatilis) and freshwater‐resident nonparasitic brook lamprey (Lampetra planeri) distributed in 10 sympatric and parapatric population pairs in France. Mating trials, parentage analyses and artificial fertilizations demonstrated a low level of reproductive isolation between species even though size‐assortative mating may contribute to isolation. Most parapatric population pairs were strongly differentiated due to the joint effects of geographic distance and barriers to migration. In contrast, we found variable levels of gene flow between sympatric populations ranging from panmixia to moderate differentiation, which indicates a gradient of divergence with some population pairs that may correspond to alternative morphs or ecotypes of a single species and others that remain partially isolated. Ecologically based divergent selection may explain these variable levels of divergence among sympatric population pairs, but incomplete genome swamping following secondary contact could have also played a role. Overall, this study illustrates how highly differentiated phenotypes can be maintained despite high levels of gene flow that limit the progress towards speciation.     

Divergence and reproductive isolation in the early stages of speciation

To understand speciation we need to identify the factors causing divergence between natural populations. The traditional approach to gaining such insights has been to focus on a particular theory and ask whether observed patterns of reproductive isolation between populations or species are consistent with the hypothesis in question. However, such studies are few and they do not allow us to compare between hypotheses, so often we cannot determine the relative contribution to divergence of different potential factors. Here, I describe a study of patterns of phenotypic divergence and premating and postmating isolation between populations of the grasshopper Chorthippus parallelus. Information on the phylogeographic relationships of the populations means that a priori predictions from existing hypotheses for the evolution of reproductive isolation can be compared with observations. I assess the relative contributions to premating isolation, postmating isolation and phenotypic divergence of long periods of allopatry, adaptation to different environments and processes associated with colonisation (such as population bottlenecks). Likelihood analysis reveals that long periods of allopatry in glacial refugia are associated with postmating reproductive isolation, but not premating isolation, which is more strongly associated with colonisation. Neither premating nor postmating isolation is higher between populations differing in potential environmental selection pressures. There are only weak correlations between patterns of genetic divergence and phenotypic divergence and no correlation between premating and postmating isolation. This suggests that the potential of a taxon to exercise mate choice may affect the types of factor that promote speciation in that group. I discuss the advantages and disadvantages of the general approach of simultaneously testing competing hypotheses for the evolution of reproductive isolation.     

Parallel genetic divergence among coastal–marine ecotype pairs of European anchovy explained by differential introgression after secondary contact

Ecophenotypic differentiation among replicate ecotype pairs within a species complex is often attributed to independent outcomes of parallel divergence driven by adaptation to similar environmental contrasts. However, the extent to which parallel phenotypic and genetic divergence patterns have emerged independently is increasingly questioned by population genomic studies. Here, we document the extent of genetic differentiation within and among two geographic replicates of the coastal and marine ecotypes of the European anchovy (Engraulis encrasicolus) gathered from Atlantic and Mediterranean locations. Using a genome‐wide data set of RAD‐derived SNPs, we show that habitat type (marine vs. coastal) is the most important component of genetic differentiation among populations of anchovy. By analysing the joint allele frequency spectrum of each coastal–marine ecotype pair, we show that genomic divergence patterns between ecotypes can be explained by a postglacial secondary contact following a long period of allopatric isolation (c. 300 kyrs). We found strong support for a model including heterogeneous migration among loci, suggesting that secondary gene flow has eroded past differentiation at different rates across the genome. Markers experiencing reduced introgression exhibited strongly correlated differentiation levels among Atlantic and Mediterranean regions. These results support that partial reproductive isolation and parallel genetic differentiation among replicate pairs of anchovy ecotypes are largely due to a common divergence history prior to secondary contact. They moreover provide comprehensive insights into the origin of a surprisingly strong fine‐scale genetic structuring in a high gene flow marine fish, which should improve stock management and conservation actions.     

Parallel allochronic divergence in a winter moth due to disruption of reproductive period by winter harshness

The disruption of reproductive timing by climatic harshness may result in the temporal isolation of conspecific populations and, ultimately, in speciation. However, whether temporal isolation alone can act as the force initiating speciation and how often the same type of climatic disruption results in the divergence of allochronic populations in a lineage are largely unknown. The reproductive period of the winter geometrid moth Inurois punctigera is separated into early and late winter in habitats with severe winters, but not in habitats with mild winters, suggesting that the reproductive season is disrupted by the harshness of the mid-winter period. Here, we show that sympatric pairs of early- and late-winter populations that differ in origin exist in different regions, suggesting a parallel divergence of reproductive timing. In each region, significant genetic differentiation exists between these early- and late-winter populations, suggesting that the temporal reproductive isolation has persisted. Moreover, we demonstrate that the temporal isolation, in comparison with geographic isolation, contributes greatly to the genetic differentiation among geographic and temporal populations by an analysis of molecular variance and by a comparison of genetic differentiations (F(ST) ) between geographic populations with and without difference in reproductive season. Our results suggest that adaptive divergence of allochronically reproducing populations has occurred independently in different regions, implying the generality of the role of temporal isolation in initiating speciation in a winter moth lineage.     

The tempo of trait divergence in geographic isolation: Avian speciation across the Marañon Valley of Peru

Geographic isolation is considered essential to most speciation events, but our understanding of what controls the pace and degree of phenotypic divergence among allopatric populations remains poor. Why do some taxa exhibit phenotypic differentiation across barriers to dispersal, whereas others do not? To test factors controlling phenotypic divergence in allopatry, we employed a comparative phylogeographic approach consisting of replicates of ecologically similar Andean bird species isolated across a major biogeographic barrier, the Marañon Valley of Peru. Our study design leverages variation among codistributed taxa in their degree of plumage, morphometric, and vocal differentiation across the Marañon to examine the tempo of phenotypic evolution. We found that substantial plumage differences between populations required roughly two million years to evolve. In contrast, morphometric trait evolution showed greater idiosyncrasy and stasis. Our results demonstrate that despite a large degree of idiosyncrasy in the relationship between genetic and phenotypic divergence across taxa and environments, comparative studies within regions may reveal predictability in the pace of phenotypic divergence. Our results also suggest that social selection is important for driving differentiation of populations found in similar environments.     

Glacial cycles drive rapid divergence of cryptic field vole species

Understanding the factors that contribute to the generation of reproductively isolated forms is a fundamental goal of evolutionary biology. Cryptic species are an especially interesting challenge to study in this context since they lack obvious morphological differentiation that provides clues to adaptive divergence that may drive reproductive isolation. Geographical isolation in refugial areas during glacial cycling is known to be important for generating genetically divergent populations, but its role in the origination of new species is still not fully understood and likely to be situation dependent. We combine analysis of 35,434 single‐nucleotide polymorphisms (SNPs) with environmental niche modeling (ENM) to investigate genomic and ecological divergence in three cryptic species formerly classified as the field vole (Microtus agrestis). The SNPs demonstrate high genomic divergence (pairwise F_ST values of 0.45–0.72) and little evidence of gene flow among the three field vole cryptic species, and we argue that genetic drift may have been a particularly important mechanism for divergence in the group. The ENM reveals three areas as potential glacial refugia for the cryptic species and differing climatic niches, although with spatial overlap between species pairs. This evidence underscores the role that glacial cycling has in promoting genetic differentiation and reproductive isolation by subdivision into disjunct distributions at glacial maxima in areas relatively close to ice sheets. Future investigation of the intrinsic barriers to gene flow between the field vole cryptic species is required to fully assess the mechanisms that contribute to reproductive isolation. In addition, the Portuguese field vole (M. rozianus) shows a high inbreeding coefficient and a restricted climatic niche, and warrants investigation into its conservation status.     

Divergent host-plant use promotes reproductive isolation among cynipid gall wasp populations

Ecological speciation occurs when reproductive isolation evolves as a consequence of divergent natural selection among environments. A direct prediction of this process is that ecologically divergent pairs of populations will exhibit greater reproductive isolation than ecologically similar pairs of populations. By comparing allopatric populations of the cynipid gall wasp Belonocnema treatae infesting Quercus virginiana and Quercus geminata, we tested the role that divergent host use plays in generating ecological divergence and sexual isolation. We found differences in body size and gall structure associated with divergent host use, but no difference in neutral genetic divergence between populations on the same or different host plant. We observed significant assortative mating between populations from alternative host plants but not between allopatric populations on the same host plant. Thus, we provide evidence that divergent host use promotes speciation among gall wasp populations.     

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.     

Phenotypic divergence of the common toad (Bufo bufo) along an altitudinal gradient: evidence for local adaptation

Variation in the environment can induce different patterns of genetic and phenotypic differentiation among populations. Both neutral processes and selection can influence phenotypic differentiation. Altitudinal phenotypic variation is of particular interest in disentangling the interplay between neutral processes and selection in the dynamics of local adaptation processes but remains little explored. We conducted a common garden experiment to study the phenotypic divergence in larval life-history traits among nine populations of the common toad (Bufo bufo) along an altitudinal gradient in France. We further used correlation among population pairwise estimates of quantitative trait (Q_ST) and neutral genetic divergence (F_ST from neutral microsatellite markers), as well as altitudinal difference, to estimate the relative role of divergent selection and neutral genetic processes in phenotypic divergence. We provided evidence for a neutral genetic differentiation resulting from both isolation by distance and difference in altitude. We found evidence for phenotypic divergence along the altitudinal gradient (faster development, lower growth rate and smaller metamorphic size). The correlation between pairwise Q_STs–F_STs and altitude differences suggested that this phenotypic differentiation was most likely driven by altitude-mediated selection rather than by neutral genetic processes. Moreover, we found different divergence patterns for larval traits, suggesting that different selective agents may act on these traits and/or selection on one trait may constrain the evolution on another through genetic correlation. Our study highlighted the need to design more integrative studies on the common toad to unravel the underlying processes of phenotypic divergence and its selective agents in the context of environmental clines.     

Flyway structure in the circumpolar greater white‐fronted goose

Dispersal and migratory behavior are influential factors in determining how genetic diversity is distributed across the landscape. In migratory species, genetic structure can be promoted via several mechanisms including fidelity to distinct migratory routes. Particularly within North America, waterfowl management units have been delineated according to distinct longitudinal migratory flyways supported by banding data and other direct evidence. The greater white‐fronted goose (Anser albifrons) is a migratory waterfowl species with a largely circumpolar distribution consisting of up to six subspecies roughly corresponding to phenotypic variation. We examined the rangewide population genetic structure of greater white‐fronted geese using mtDNA control region sequence data and microsatellite loci from 23 locales across North America and Eurasia. We found significant differentiation in mtDNA between sampling locales with flyway delineation explaining a significant portion of the observed genetic variation (~12%). This is concordant with band recovery data which shows little interflyway or intercontinental movements. However, microsatellite loci revealed little genetic structure suggesting a panmictic population across most of the Arctic. As with many high‐latitude species, Beringia appears to have played a role in the diversification of this species. A common Beringian origin of North America and Asian populations and a recent divergence could at least partly explain the general lack of structure at nuclear markers. Further, our results do not provide strong support for the various taxonomic proposals for this species except for supporting the distinctness of two isolated breeding populations within Cook Inlet, Alaska (A. a. elgasi) and Greenland (A. a. flavirostris), consistent with their subspecies status.     

Cultural inheritance drives site fidelity and migratory connectivity in a long-distance migrant

Cultural transmission is thought to be a mechanism by which migratory animals settle into habitats, but little evidence exists in wild populations because of the difficulty of following individuals over successive generations and wide geographical distances. Cultural inheritance of migration routes represents a mechanism whereby geographical isolation can arise between separate groups and could constrain individuals to potentially suboptimal sites within their range. Conversely, adopting the parental migratory route in adult life, rather than dispersing randomly, may increase an individual's reproductive success because that strategy has already been proven to allow successful breeding. We combined a pedigree of related light-bellied Brent geese (Branta bernicla hrota) with 6 years of observations of marked birds to calculate the dispersal distances of adult offspring from their parents in both Ireland and Iceland. In both countries, the majority of offspring were found to recruit into or near their parental sites, indicating migratory connectivity in the flyway. Despite this kin structure, we found no evidence of genetic differentiation using genotype data from 1127 individuals across 15 microsatellite loci. We suggest that the existence of migratory connectivity of subpopulations is far more common than previous research indicates and that cultural information may play an important role in structuring reproductive isolation among them.     

Genome-wide patterns of divergence during speciation: the lake whitefish case study

The nature, size and distribution of the genomic regions underlying divergence and promoting reproductive isolation remain largely unknown. Here, we summarize ongoing efforts using young (12 000 yr BP) species pairs of lake whitefish (Coregonus clupeaformis) to expand our understanding of the initial genomic patterns of divergence observed during speciation. Our results confirmed the predictions that: (i) on average, phenotypic quantitative trait loci (pQTL) show higher F(ST) values and are more likely to be outliers (and therefore candidates for being targets of divergent selection) than non-pQTL markers; (ii) large islands of divergence rather than small independent regions under selection characterize the early stages of adaptive divergence of lake whitefish; and (iii) there is a general trend towards an increase in terms of numbers and size of genomic regions of divergence from the least (East L.) to the most differentiated species pair (Cliff L.). This is consistent with previous estimates of reproductive isolation between these species pairs being driven by the same selective forces responsible for environment specialization. Altogether, dwarf and normal whitefish species pairs represent a continuum of both morphological and genomic differentiation contributing to ecological speciation. Admittedly, much progress is still required to more finely map and circumscribe genomic islands of speciation. This will be achieved through the use of next generation sequencing data but also through a better quantification of phenotypic traits moulded by selection as organisms adapt to new environmental conditions.