THE EFFICACY OF DIVERGENCE HITCHHIKING IN GENERATING GENOMIC ISLANDS DURING ECOLOGICAL SPECIATION

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

The genetic mosaic suggests a new role for hitchhiking in ecological speciation

Early in ecological speciation, the genomically localized effects of divergent selection cause heterogeneity among loci in divergence between incipient species. We call this pattern of genomic variability in divergence the 'genetic mosaic of speciation'. Previous studies have used F(ST) outliers as a way to identify divergently selected genomic regions, but the nature of the relationship between outlier loci and quantitative trait loci (QTL) involved in reproductive isolation has not yet been quantified. Here, we show that F(ST) outliers between a pair of incipient species are significantly clustered around QTL for traits that cause ecologically based reproductive isolation. Around these key QTL, extensive 'divergence hitchhiking' occurs because reduced inter-race mating and negative selection decrease the opportunity for recombination between chromosomes bearing different locally adapted QTL alleles. Divergence hitchhiking is likely to greatly increase the opportunity for speciation in populations that are sympatric, regardless of whether initial divergence was sympatric or allopatric. Early in ecological speciation, analyses of population structure, gene flow or phylogeography based on different random or arbitrarily chosen neutral markers should be expected to conflict--only markers in divergently selected genomic regions will reveal the evolutionary history of adaptive divergence and ecologically based reproductive isolation. Species retain mosaic genomes for a very long time, and gene exchange in hybrid zones can vary dramatically among loci. However, in hybridizing species, the genomic regions that affect ecologically based reproductive isolation are difficult to distinguish from regions that have diverged for other reasons.     

Establishment of new mutations under divergence and genome hitchhiking

Theoretical models addressing genome-wide patterns of divergence during speciation are needed to help us understand the evolutionary processes generating empirical patterns. Here, we examine a critical issue concerning speciation-with-gene flow: to what degree does physical linkage (r < 0.5) of new mutations to already diverged genes aid the build-up of genomic islands of differentiation? We used simulation and analytical approaches to partition the probability of establishment for a new divergently selected mutation when the mutation (i) is the first to arise in an undifferentiated genome (the direct effect of selection), (ii) arises unlinked to any selected loci (r = 0.5), but within a genome that has some already diverged genes (the effect of genome-wide reductions in gene flow for facilitating divergence, which we term 'genome hitchhiking'), and (iii) arises in physical linkage to a diverged locus (divergence hitchhiking). We find that the strength of selection acting directly on a new mutation is generally the most important predictor for establishment, with divergence and genomic hitchhiking having smaller effects. We outline the specific conditions under which divergence and genome hitchhiking can aid mutation establishment. The results generate predictions about genome divergence at different points in the speciation process and avenues for further work.     

Spatially explicit models of divergence and genome hitchhiking

Strong barriers to genetic exchange can exist at divergently selected loci, whereas alleles at neutral loci flow more readily between populations, thus impeding divergence and speciation in the face of gene flow. However, ‘divergence hitchhiking’ theory posits that divergent selection can generate large regions of differentiation around selected loci. ‘Genome hitchhiking’ theory suggests that selection can also cause reductions in average genome‐wide rates of gene flow, resulting in widespread genomic divergence (rather than divergence only around specific selected loci). Spatial heterogeneity is ubiquitous in nature, yet previous models of genetic barriers to gene flow have explored limited combinations of spatial and selective scenarios. Using simulations of secondary contact of populations, we explore barriers to gene flow in various selective and spatial contexts in continuous, two‐dimensional, spatially explicit environments. In general, the effects of hitchhiking are strongest in environments with regular spatial patterning of starkly divergent habitat types. When divergent selection is very strong, the absence of intermediate habitat types increases the effects of hitchhiking. However, when selection is moderate or weak, regular (vs. random) spatial arrangement of habitat types becomes more important than the presence of intermediate habitats per se. We also document counterintuitive processes arising from the stochastic interplay between selection, gene flow and drift. Our results indicate that generalization of results from two‐deme models requires caution and increase understanding of the genomic and geographic basis of population divergence.     

Divergence hitchhiking and the spread of genomic isolation during ecological speciation-with-gene-flow

In allopatric populations, geographical separation simultaneously isolates the entire genome, allowing genetic divergence to accumulate virtually anywhere in the genome. In sympatric populations, however, the strong divergent selection required to overcome migration produces a genetic mosaic of divergent and non-divergent genomic regions. In some recent genome scans, each divergent genomic region has been interpreted as an independent incidence of migration/selection balance, such that the reduction of gene exchange is restricted to a few kilobases around each divergently selected gene. I propose an alternative mechanism, 'divergence hitchhiking' (DH), in which divergent selection can reduce gene exchange for several megabases around a gene under strong divergent selection. Not all genes/markers within a DH region are divergently selected, yet the entire region is protected to some degree from gene exchange, permitting genetic divergence from mechanisms other than divergent selection to accumulate secondarily. After contrasting DH and multilocus migration/selection balance (MM/SB), I outline a model in which genomic isolation at a given genomic location is jointly determined by DH and genome-wide effects of the progressive reduction in realized migration, then illustrate DH using data from several pairs of incipient species in the wild.     

Genomics of Divergence along a Continuum of Parapatric Population Differentiation

The patterns of genomic divergence during ecological speciation are shaped by a combination of evolutionary forces. Processes such as genetic drift, local reduction of gene flow around genes causing reproductive isolation, hitchhiking around selected variants, variation in recombination and mutation rates are all factors that can contribute to the heterogeneity of genomic divergence. On the basis of 60 fully sequenced three-spined stickleback genomes, we explore these different mechanisms explaining the heterogeneity of genomic divergence across five parapatric lake and river population pairs varying in their degree of genetic differentiation. We find that divergent regions of the genome are mostly specific for each population pair, while their size and abundance are not correlated with the extent of genome-wide population differentiation. In each pair-wise comparison, an analysis of allele frequency spectra reveals that 25–55% of the divergent regions are consistent with a local restriction of gene flow. Another large proportion of divergent regions (38–75%) appears to be mainly shaped by hitchhiking effects around positively selected variants. We provide empirical evidence that alternative mechanisms determining the evolution of genomic patterns of divergence are not mutually exclusive, but rather act in concert to shape the genome during population differentiation, a first necessary step towards ecological speciation.     

Localizing FST outliers on a QTL map reveals evidence for large genomic regions of reduced gene exchange during speciation‐with‐gene‐flow

Populations that maintain phenotypic divergence in sympatry typically show a mosaic pattern of genomic divergence, requiring a corresponding mosaic of genomic isolation (reduced gene flow). However, mechanisms that could produce the genomic isolation required for divergence‐with‐gene‐flow have barely been explored, apart from the traditional localized effects of selection and reduced recombination near centromeres or inversions. By localizing F_ST outliers from a genome scan of wild pea aphid host races on a Quantitative Trait Locus (QTL) map of key traits, we test the hypothesis that between‐population recombination and gene exchange are reduced over large ‘divergence hitchhiking’ (DH) regions. As expected under divergence hitchhiking, our map confirms that QTL and divergent markers cluster together in multiple large genomic regions. Under divergence hitchhiking, the nonoutlier markers within these regions should show signs of reduced gene exchange relative to nonoutlier markers in genomic regions where ongoing gene flow is expected. We use this predicted difference among nonoutliers to perform a critical test of divergence hitchhiking. Results show that nonoutlier markers within clusters of F_ST outliers and QTL resolve the genetic population structure of the two host races nearly as well as the outliers themselves, while nonoutliers outside DH regions reveal no population structure, as expected if they experience more gene flow. These results provide clear evidence for divergence hitchhiking, a mechanism that may dramatically facilitate the process of speciation‐with‐gene‐flow. They also show the power of integrating genome scans with genetic analyses of the phenotypic traits involved in local adaptation and population divergence.     

Divergent selection and heterogeneous genomic divergence

Levels of genetic differentiation between populations can be highly variable across the genome, with divergent selection contributing to such heterogeneous genomic divergence. For example, loci under divergent selection and those tightly physically linked to them may exhibit stronger differentiation than neutral regions with weak or no linkage to such loci. Divergent selection can also increase genome‐wide neutral differentiation by reducing gene flow (e.g. by causing ecological speciation), thus promoting divergence via the stochastic effects of genetic drift. These consequences of divergent selection are being reported in recently accumulating studies that identify: (i) ‘outlier loci’ with higher levels of divergence than expected under neutrality, and (ii) a positive association between the degree of adaptive phenotypic divergence and levels of molecular genetic differentiation across population pairs [‘isolation by adaptation’ (IBA)]. The latter pattern arises because as adaptive divergence increases, gene flow is reduced (thereby promoting drift) and genetic hitchhiking increased. Here, we review and integrate these previously disconnected concepts and literatures. We find that studies generally report 5–10% of loci to be outliers. These selected regions were often dispersed across the genome, commonly exhibited replicated divergence across different population pairs, and could sometimes be associated with specific ecological variables. IBA was not infrequently observed, even at neutral loci putatively unlinked to those under divergent selection. Overall, we conclude that divergent selection makes diverse contributions to heterogeneous genomic divergence. Nonetheless, the number, size, and distribution of genomic regions affected by selection varied substantially among studies, leading us to discuss the potential role of divergent selection in the growth of regions of differentiation (i.e. genomic islands of divergence), a topic in need of future investigation.     

Isolation by adaptation in Neochlamisus leaf beetles: host-related selection promotes neutral genomic divergence

This study tests how divergent natural selection promotes genomic differentiation during ecological speciation. Specifically, we use adaptive ecological divergence (here, population divergence in host plant use and preference) as a proxy for selection strength and evaluate the correlation between levels of adaptive and genetic differentiation across pairwise population comparisons. Positive correlations would reveal the pattern predicted by our hypothesis, that of 'isolation by adaptation' (IBA). Notably, IBA is predicted not only for selected loci but also for neutral loci. This may reflect the effects of divergent selection on neutral loci that are 'loosely linked' to divergently selected loci or on geneflow restriction that facilitates genetic drift at all loci, including neutral loci that are completely unlinked to those evolving under divergent selection. Here, we evaluate IBA in maple- and willow-associated populations of Neochlamisus bebbianae leaf beetles. To do so, we collected host preference data to construct adaptive divergence indices and used AFLPs (amplified fragment length polymorphisms) and mitochondrial sequences to quantify genetic differentiation. Partial Mantel tests showed significant IBA in 'pooled' analyses of putatively neutral and of putatively selected ('outlier') AFLP loci. This pattern was also recovered in 12% of 'locus-specific' analyses that separately evaluated genetic differentiation at individual neutral loci. These results provided evidence for widespread effects of selection on neutral genomic divergence. Our collective findings indicate that host-related selection may play important roles in the population genomic differentiation of both neutral and selected gene regions in herbivorous insects.     

Bayesian estimation of genomic clines

We developed a Bayesian genomic cline model to study the genetic architecture of adaptive divergence and reproductive isolation between hybridizing lineages. This model quantifies locus‐specific patterns of introgression with two cline parameters that describe the probability of locus‐specific ancestry as a function of genome‐wide admixture. ‘Outlier’ loci with extreme patterns of introgression relative to most of the genome can be identified. These loci are potentially associated with adaptive divergence or reproductive isolation. We simulated genetic data for admixed populations that included neutral introgression, as well as loci that were subject to directional, epistatic or underdominant selection, and analysed these data using the Bayesian genomic cline model. Under many demographic conditions, underdominance or directional selection had detectable and predictable effects on cline parameters, and ‘outlier’ loci were greatly enriched for genetic regions affected by selection. We also analysed previously published genetic data from two transects through a hybrid zone between Mus domesticus and M. musculus. We found considerable variation in rates of introgression across the genome and particularly low rates of introgression for two X‐linked markers. There were similarities and differences in patterns of introgression between the two transects, which likely reflects a combination of stochastic variability because of genetic drift and geographic variation in the genetic architecture of reproductive isolation. By providing a robust framework to quantify and compare patterns of introgression among genetic regions and populations, the Bayesian genomic cline model will advance our understanding of the genetics of reproductive isolation and the speciation process.     

The genomics of adaptation,divergence and speciation: a congealing theory

In this issue, Flaxman et al. ( 2014 ) report the results of sophisticated whole‐genome simulations of speciation with gene flow, enhancing our understanding of the process by building on previous single‐locus, multilocus and analytical works. Their findings provide us with new insights about how genomes can diverge and the importance of statistical and chromosomal linkage in facilitating reproductive isolation. The authors characterize the conditions under which, even with high gene flow and weak divergent selection, reproductive isolation between populations can occur due to the emergent stochastic process of genomewide congealing, where numerous statistically or physically linked loci of small effect allow selection to limit effective migration rates. The initial congealing event can occur within a broad range conditions, and once initiated, the self‐reinforcing process leads to rapid divergence and ultimately two reproductively isolated populations. Flaxman et al.'s ( 2014 ) work is a valuable contribution to our understanding of speciation with gene flow and in making a more predictive field of evolutionary genomics and speciation.     

The hypothesis of immune testing of partners—coordinated adaptations and changes in mating preferences

Reproductive isolation plays the key role in speciation. According to the prevailing ideas, the main speciation mechanism is gradual accumulation of genetic differences in isolated populations (allopatric phase of speciation) based on mutations, selection, and genetic drift. In this case, reproductive isolation emerges as an occasional byproduct of adaptation to different conditions (ecological speciation) or accumulation of random changes in the gene pool resulting from long-term isolation. Pure sympatric speciation assumes isolation as a direct product of selection (divergent or disruptive selection) that favors individuals selectively mating with their likes. A third possibility is substantiated below. We believe that isolation can be a regular and determined product rather than occasional byproduct of divergence. It can rely on the friend/foe discrimination mechanisms, some of which can be “immune-based” and compare the partner’s and own properties (signaling molecules, pheromones, and other antigens in a broad sense). Antigens of the major histocompatibility complex (MHC) can play a substantial role in such testing of potential mates.     

The genomics of speciation-with-gene-flow

The emerging field of speciation genomics is advancing our understanding of the evolution of reproductive isolation from the individual gene to a whole-genome perspective. In this new view it is important to understand the conditions under which 'divergence hitchhiking' associated with the physical linkage of gene regions, versus 'genome hitchhiking' associated with reductions in genome-wide rates of gene flow caused by selection, can enhance speciation-with-gene-flow. We describe here a theory predicting four phases of speciation, defined by changes in the relative effectiveness of divergence and genome hitchhiking, and review empirical data in light of the theory. We outline future directions, emphasizing the need to couple next-generation sequencing with selection, transplant, functional genomics, and mapping studies. This will permit a natural history of speciation genomics that will help to elucidate the factors responsible for population divergence and the roles that genome structure and different forms of hitchhiking play in facilitating the genesis of new biodiversity.     

THE BIOLOGY OF SPECIATION

Since Darwin published the “Origin,” great progress has been made in our understanding of speciation mechanisms. The early investigations by Mayr and Dobzhansky linked Darwin's view of speciation by adaptive divergence to the evolution of reproductive isolation, and thus provided a framework for studying the origin of species. However, major controversies and questions remain, including: When is speciation nonecological? Under what conditions does geographic isolation constitute a reproductive isolating barrier? and How do we estimate the “importance” of different isolating barriers? Here, we address these questions, providing historical background and offering some new perspectives. A topic of great recent interest is the role of ecology in speciation. “Ecological speciation” is defined as the case in which divergent selection leads to reproductive isolation, with speciation under uniform selection, polyploid speciation, and speciation by genetic drift defined as “nonecological.” We review these proposed cases of nonecological speciation and conclude that speciation by uniform selection and polyploidy normally involve ecological processes. Furthermore, because selection can impart reproductive isolation both directly through traits under selection and indirectly through pleiotropy and linkage, it is much more effective in producing isolation than genetic drift. We thus argue that natural selection is a ubiquitous part of speciation, and given the many ways in which stochastic and deterministic factors may interact during divergence, we question whether the ecological speciation concept is useful. We also suggest that geographic isolation caused by adaptation to different habitats plays a major, and largely neglected, role in speciation. We thus provide a framework for incorporating geographic isolation into the biological species concept (BSC) by separating ecological from historical processes that govern species distributions, allowing for an estimate of geographic isolation based upon genetic differences between taxa. Finally, we suggest that the individual and relative contributions of all potential barriers be estimated for species pairs that have recently achieved species status under the criteria of the BSC. Only in this way will it be possible to distinguish those barriers that have actually contributed to speciation from those that have accumulated after speciation is complete. We conclude that ecological adaptation is the major driver of reproductive isolation, and that the term “biology of speciation,” as proposed by Mayr, remains an accurate and useful characterization of the diversity of speciation mechanisms.     

Impacts of climate warming on hybrid zone movement: Geographically diffuse and biologically porous “species borders”

Abstract The ecology and evolutionary biology of insect–plant associations has realized extensive attention, especially during the past 60 years. The classifications (categorical designations) of continuous variation in biodiversity, ranging from global patterns (e.g., latitudinal gradients in species richness/diversity and degree of herbivore feeding specialization) to localized insect–plant associations that span the biospectrum from polyphenisms, polymorphisms, biotypes, demes, host races, to cryptic species, remain academically contentious. Semantic and biosystematic (taxonomical) disagreements sometimes detract from more important ecological and evolutionary processes that drive diversification, the dynamics of gene flow and local extinctions. This review addresses several aspects of insect specialization, host‐associated divergence and ecological (including “hybrid”) speciation, with special reference to the climate warming impacts on species borders of hybridizing swallowtail butterflies (Papilionidae). Interspecific hybrid introgression may result in collapse of multi‐species communities or increase species numbers via homoploid hybrid speciation. We may see diverging, merging, or emerging genotypes across hybrid zones, all part of the ongoing processes of evolution. Molecular analyses of genetic mosaics and genomic dynamics with “divergence hitchhiking”, combined with ecological, ethological and physiological studies of “species porosity”, have already begun to unveil some answers for some important ecological/evolutionary questions. (i) How rapidly can host‐associated divergence lead to new species (and why doesn't it always do so, e.g., resulting in “incomplete” speciation)? (ii) How might “speciation genes” function, and how/where would we find them? (iii) Can oscillations from specialists to generalists and back to specialists help explain global diversity in herbivorous insects? (iv) How could recombinant interspecific hybridization lead to divergence and speciation? From ancient phytochemically defined angiosperm affiliations to recent and very local geographical mosaics, the Papilionidae (swallowtail butterflies) have provided a model for enhanced understanding of ecological patterns and evolutionary processes, including host‐associated genetic divergence, genomic mosaics, genetic hitchhiking and sex‐linked speciation genes. Apparent homoploid hybrid speciation in Papilio appears to have been catalyzed by climate warming‐induced interspecific introgression of some, but not all, species diagnostic traits, reflecting strong divergent selection (discordant), especially on the Z (= X) chromosome. Reproductive isolation of these novel recombinant hybrid genotypes appears to be accomplished via a delayed post‐diapause emergence or temporal isolation, and is perhaps aided by the thermal landscape. Changing thermal landscapes appear to have created (and may destroy) novel recombinant hybrid genotypes and hybrid species.     

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.     

Recombinant hybrids retain heterozygosity at many loci: new insights into the genomics of reproductive isolation in Populus

The maintenance of species barriers in the face of gene flow is often thought to result from strong selection against intermediate genotypes, thereby preserving genetic differentiation. Most speciation genomic studies thus aim to identify exceptionally divergent loci between populations, but divergence will be affected by many processes other than reproductive isolation (RI) and speciation. Through genomic studies of recombinant hybrids sampled in the wild, genetic variation associated with RI can be observed in situ, because selection against incompatible genotypes will leave detectable patterns of variation in the hybrid genomes. To better understand the mechanisms directly involved in RI, we investigated three natural ‘replicate’ hybrid zones between two divergent Populus species via locus‐specific patterns of ancestry across recombinant hybrid genomes. As expected, genomic patterns in hybrids and their parental species were consistent with the presence of underdominant selection at several genomic regions. Surprisingly, many loci displayed greatly increased between‐species heterozygosity in recombinant hybrids despite striking genetic differentiation between the parental genomes, the opposite of what would be expected with selection against intermediate genotypes. Only a limited, reproducible set of genotypic combinations was present in hybrid genomes across localities. In the absence of clearly delimited ‘hybrid habitats’, our results suggest that complex epistatic interactions within genomes play an important role in advanced stages of RI between these ecologically divergent forest trees. This calls for more genomic studies that test for unusual patterns of genomic ancestry in hybridizing species.     

Interpreting the genomic landscape of speciation: a road map for finding barriers to gene flow

Speciation, the evolution of reproductive isolation among populations, is continuous, complex, and involves multiple, interacting barriers. Until it is complete, the effects of this process vary along the genome and can lead to a heterogeneous genomic landscape with peaks and troughs of differentiation and divergence. When gene flow occurs during speciation, barriers restricting gene flow locally in the genome lead to patterns of heterogeneity. However, genomic heterogeneity can also be produced or modified by variation in factors such as background selection and selective sweeps, recombination and mutation rate variation, and heterogeneous gene density. Extracting the effects of gene flow, divergent selection and reproductive isolation from such modifying factors presents a major challenge to speciation genomics. We argue one of the principal aims of the field is to identify the barrier loci involved in limiting gene flow. We first summarize the expected signatures of selection at barrier loci, at the genomic regions linked to them and across the entire genome. We then discuss the modifying factors that complicate the interpretation of the observed genomic landscape. Finally, we end with a road map for future speciation research: a proposal for how to account for these modifying factors and to progress towards understanding the nature of barrier loci. Despite the difficulties of interpreting empirical data, we argue that the availability of promising technical and analytical methods will shed further light on the important roles that gene flow and divergent selection have in shaping the genomic landscape of speciation.     

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.     

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.