Integrating QTL mapping and genome scans towards the characterization of candidate loci under parallel selection in the lake whitefish (Coregonus clupeaformis)

As natural selection must act on underlying genetic variation, discovering the number and location of loci under the influence of selection is imperative towards understanding adaptive divergence in evolving populations. Studies employing genome scans have hypothesized that the action of divergent selection should reduce gene flow at the genomic locations implicated in adaptation and speciation among natural populations, yet once 'outlier' patterns of variation have been identified the function and role of such loci needs to be confirmed. We integrated adaptive QTL mapping and genomic scans among diverging sympatric pairs of the lake whitefish (Coregonus clupeaformis) species complex in order to test the hypothesis that differentiation between dwarf and normal ecotypes at growth-associated QTL was maintained by directional selection. We found evidence of significantly high levels of molecular divergence among eight growth QTL where two of the strongest candidate loci under the influence of directional selection exhibited parallel reductions of gene flow over multiple populations.     

CORRELATED TROPHIC SPECIALIZATION AND GENETIC DIVERGENCE IN SYMPATRIC LAKE WHITEFISH ECOTYPES (COREGONUS CLUPEAFORMIS): SUPPORT FOR THE ECOLOGICAL SPECIATION HYPOTHESIS

There is ample empirical evidence that phenotypic diversification in an adaptive radiation is the outcome of divergent natural selection related to differential resource use. In contrast, the role of ecological forces in favoring and maintaining reproductive isolation in nature remains poorly understood. If the same forces driving phenotypic divergence are also responsible for speciation, one would predict a correlation between the extent of trophic specialization (reflecting variable intensity of divergent natural selection) and that of reproductive isolation being reached in a given environment. We tested this hypothesis by comparing the extent of morphological and genetic differentiation between sympatric dwarf and normal whitefish ecotypes (Coregonus sp.) from six lakes of the St. John River basin (eastern Canada and northern Maine). Eight meristic variables, 19 morphometric variables, and six microsatellite loci were used to quantify morphological and genetic differentiation, respectively. Dwarf and normal ecotypes in each lake differed primarily by traits related to trophic specialization, but the extent of differentiation varied among lakes. Significant but variable genetic divergence between ecotypes within lakes was also observed. A negative correlation was observed between the extent of gene flow between ecotypes within a lake and that of their morphological differentiation in trophic-related traits. The extent of reproductive isolation reached between dwarf and normal whitefish ecotypes appears to be driven by the potential for occupying distinct trophic niches and, thus, by the same selective forces driving tropic specialization in each lake. These results therefore support the hypothesis of ecological speciation.     

Generic scan using AFLP markers as a means to assess the role of directional selection in the divergence of sympatric whitefish ecotypes

Under the ecological theory of adaptive radiation, adaptation and reproductive isolation are thought to evolve as a result of divergent natural selection. Accordingly, elucidating the genetic basis of these processes is essential toward understanding the role of selection in shaping biological diversity. In this respect, the number of genes that evolved by selection remains contentious. To address this issue, the pattern of genetic differentiation obtained using 440 AFLP loci was compared with that expected under neutrality in four sympatric pairs of lake whitefish ecotypes that evolved adaptive phenotypic differences associated with the exploitation of distinct ecological niches. On average, 14 loci showed restricted gene flow relative to neutral expectation, suggesting a role of directional selection on their divergence. Among all loci that are most likely under directional selection, six exhibited parallel patterns of divergence, which provided further support for the role of selection in driving their divergence. Overall, these results indicate that only a small proportion of scored AFLP loci (between 1.4% and 3.2%) might be linked to genes implicated in the adaptive radiation of lake whitefish.     

Role of epibenthic resource opportunities in the parallel evolution of lake whitefish species pairs (Coregonus sp.)

Parallel evolution of a dwarf and normal whitefish has been documented in six post-glacial lakes. Here, we relate the structure and seasonal variations of the epibenthic invertebrate communities to the extent of phenotypic differentiation in these species pairs. The highest phenotypic differentiation occurs in lakes characterized by less overlap in size distribution between limnetic and epibenthic prey which could represent enhanced ecological opportunities for trophic specialization and adaptive divergence. Differences in community assemblages and seasonal variation of biotic and abiotic conditions may also play a role. Accumulating evidence indicates that strong directional selection acting on dwarf whitefish may be more important than divergent selection acting on both sympatric forms in driving whitefish phenotypic divergence and ultimately, ecological speciation. Along with Landry et al. (2007), this study supports the general hypothesis that parallelism in divergence among sympatric dwarf and normal whitefish is associated with parallelism in limnological adaptive landscape.     

Combining the analyses of introgressive hybridisation and linkage mapping to investigate the genetic architecture of population divergence in the lake whitefish (Coregonus clupeaformis, Mitchill)

Adaptation and reproductive isolation, the engines of biological diversity, are still elusive when discussing the genetic bases of speciation. Namely, the number of genes and magnitude of selection acting positively or negatively on genomic traits implicated in speciation is contentious. Here, we describe the first steps of an ongoing research program aimed at understanding the genetic bases of population divergence and reproductive isolation in the lake whitefish (Coregonus clupeaformis). A preliminary linkage map originating from a hybrid cross between dwarf and normal ecotypes is presented, whereby some of the segregating AFLP markers were found to be conserved among natural populations. Maximum-likelihood was used to estimate hybrid indices from non-diagnostic markers at 998 AFLP loci. This allowed identification of the most likely candidate loci that have been under the influence of selection during the natural hybridisation of whitefish originating from different glacial races. As some of these loci could be identified on the linkage map, the possibility that selection of traits in natural populations may eventually be correlated to specific chromosomal regions was demonstrated. The future prospects and potential of these approaches to elucidate the genetic bases of adaptation and reproductive isolation among sympatric ecotypes of lake whitefish is discussed.     

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.     

Genomic insights into the vulnerability of sympatric whitefish species flocks

The erosion of habitat heterogeneity can reduce species diversity directly but can also lead to the loss of distinctiveness of sympatric species through speciation reversal. We know little about changes in genomic differentiation during the early stages of these processes, which can be mediated by anthropogenic perturbation. Here, we analyse three sympatric whitefish species (Coregonus spp) sampled across two neighbouring and connected Swiss pre‐alpine lakes, which have been differentially affected by anthropogenic eutrophication. Our data set comprises 16,173 loci genotyped across 138 whitefish using restriction‐site associated DNA sequencing (RADseq). Our analysis suggests that in each of the two lakes, the population of a different, but ecologically similar, whitefish species declined following a recent period of eutrophication. Genomic signatures consistent with hybridization are more pronounced in the more severely impacted lake. Comparisons between sympatric pairs of whitefish species with contrasting ecology, where one is shallow benthic and the other one more profundal pelagic, reveal genomic differentiation that is largely correlated along the genome, while differentiation is uncorrelated between pairs of allopatric provenance with similar ecology. We identify four genomic loci that provide evidence of parallel divergent adaptation between the shallow benthic species and the two different more profundal species. Functional annotations available for two of those loci are consistent with divergent ecological adaptation. Our genomic analysis indicates the action of divergent natural selection between sympatric whitefish species in pre‐alpine lakes and reveals the vulnerability of these species to anthropogenic alterations of the environment and associated adaptive landscape.     

THE GENETIC ARCHITECTURE OF REPRODUCTIVE ISOLATION DURING SPECIATION‐WITH‐GENE‐FLOW IN LAKE WHITEFISH SPECIES PAIRS ASSESSED BY RAD SEQUENCING

During speciation‐with‐gene‐flow, effective migration varies across the genome as a function of several factors, including proximity of selected loci, recombination rate, strength of selection, and number of selected loci. Genome scans may provide better empirical understanding of the genome‐wide patterns of genetic differentiation, especially if the variance due to the previously mentioned factors is partitioned. In North American lake whitefish (Coregonus clupeaformis), glacial lineages that diverged in allopatry about 60,000 years ago and came into contact 12,000 years ago have independently evolved in several lakes into two sympatric species pairs (a normal benthic and a dwarf limnetic). Variable degrees of reproductive isolation between species pairs across lakes offer a continuum of genetic and phenotypic divergence associated with adaptation to distinct ecological niches. To disentangle the complex array of genetically based barriers that locally reduce the effective migration rate between whitefish species pairs, we compared genome‐wide patterns of divergence across five lakes distributed along this divergence continuum. Using restriction site associated DNA (RAD) sequencing, we combined genetic mapping and population genetics approaches to identify genomic regions resistant to introgression and derive empirical measures of the barrier strength as a function of recombination distance. We found that the size of the genomic islands of differentiation was influenced by the joint effects of linkage disequilibrium maintained by selection on many loci, the strength of ecological niche divergence, as well as demographic characteristics unique to each lake. Partial parallelism in divergent genomic regions likely reflected the combined effects of polygenic adaptation from standing variation and independent changes in the genetic architecture of postzygotic isolation. This study illustrates how integrating genetic mapping and population genomics of multiple sympatric species pairs provide a window on the speciation‐with‐gene‐flow mechanism.     

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.     

Experimental evidence for reduced hybrid viability between dwarf and normal ecotypes of lake whitefish (Coregonus clupeaformis Mitchill)

Forces driving the evolution of reproductive isolation among natural populations, as well as the mechanisms involved to maintain it, are still poorly understood. Because sympatric fish ecotypes mainly differ in phenotypic traits associated with occupying distinct trophic niches, it is generally believed that reproductive isolation is mainly driven by ecological divergent selection, excluding genome incompatibility as a basis for postmating isolation. We did cross experiments between dwarf and normal ecotypes of lake whitefish (Coregonus clupeaformis Mitchill) originating from distinct glacial refugia to test the hypothesis that their geographical isolation during the Pleistocene may have led to sufficient genetic divergence for the development of reproductive isolation between them before their secondary contact. Similar fertilization success in pure and hybrid crosses indicated the absence of gametic incompatibility between the two ecotypes. In contrast, daily embryonic mortality rates were 2.4–4.7 times higher in reciprocal hybrid crosses compared to pure crosses, which supports our working hypothesis. These results, along with previous morphological and population genetic studies, indicate that both genetic and ecological mechanisms may jointly act to promote speciation among northern freshwater fish ecotypes.     

Genomic signatures of divergent selection and speciation patterns in a ‘natural experiment’, the young parallel radiations of Nicaraguan crater lake cichlid fishes

Divergent selection is the main driving force in sympatric ecological speciation and may also play a strong role in divergence between allopatric populations. Characterizing the genome‐wide impact of divergent selection often constitutes a first step in unravelling the genetic bases underlying adaptation and ecological speciation. The Midas cichlid fish (Amphilophus citrinellus) species complex in Nicaragua is a powerful system for studying evolutionary processes. Independent colonizations of isolated young crater lakes by Midas cichlid populations from the older and great lakes of Nicaragua resulted in the repeated evolution of adaptive radiations by intralacustrine sympatric speciation. In this study we performed genome scans on two repeated radiations of crater lake species and their great lake source populations (1030 polymorphic AFLPs, n ～ 30 individuals per species). We detected regions under divergent selection (0.3% in the crater lake Xiloá flock and 1.7% in the older crater lake Apoyo radiation) that might be responsible for the sympatric diversifications. We find no evidence that the same genomic regions have been involved in the repeated evolution of parallel adaptations across crater lake flocks. However, there is some genetic parallelism apparent (seven out of 51 crater lake to great lake outlier loci are shared; 13.7%) that is associated with the allopatric divergence of both crater lake flocks. Interestingly, our results suggest that the number of outlier loci involved in sympatric and allopatric divergence increases over time. A phylogeny based on the AFLP data clearly supports the monophyly of both crater lake species flocks and indicates a parallel branching order with a primary split along the limnetic‐benthic axis in both radiations.     

MULTIPLE MODES OF SPECIATION INVOLVED IN THE PARALLEL EVOLUTION OF SYMPATRIC MORPHOTYPES OF LAKE WHITEFISH (COREGONUS CLUPEAFORMIS,SALMONIDAE)

We performed a phylogenetic analysis of mtDNA variation among seven sympatric pairs of dwarf and normal morphotypes of whitefish from northern Québec and the St. John River drainage to address three questions relevant to understanding their radiation. Are all sympatric pairs reproductively isolated? Do phylogenetic analyses confirm that sympatric whitefish morphotypes found in eastern North America represent the outcome of polyphyletic evolutionary events? If so, did all sympatric pairs from the St. John River drainage originate from the same scenario of allopatric divergence and secondary contact? The hypothesis of genetic differentiation was supported for all sympatric pairs from the St. John River drainage, whereas lack of mtDNA diversity precluded any test of reproductive isolation for northern Québec populations. Patterns of mtDNA variation confirmed that dwarf and normal morphotypes evolved in parallel among independent, yet closely related, lineages, thus providing indirect evidence for the role of natural selection in promoting phenotypic radiation in whitefish. Patterns of mtDNA diversity among sympatric pairs of the St. John River indicated a complex picture of whitefish evolution that implied sympatric divergence and multiple allopatric divergence/secondary contact events on a small geographic scale. These results suggests that ecological opportunities, namely trophic niche availability, may promote population divergence in whitefish.     

Parallelism in the oxygen transport system of the lake whitefish: the role of physiological divergence in ecological speciation

In North America, populations of lake whitefish (Coregonus clupeaformis) have evolved sympatric 'dwarf' and 'normal' ecotypes that are associated with distinct trophic niches within lakes. Trophic specialization should place diverging physiological demands on individuals, and thus, genes and phenotypes associated with energy production represent ideal candidates for studies of adaptation. Here, we test for the parallel divergence of traits involved in oxygen transport in dwarf and normal lake whitefish from Québec, Canada and Maine, USA. We observed significant differences in red blood cell morphology between the ecotypes. Specifically, dwarfs exhibited larger nuclei and a higher nucleus area/total cell area than normal whitefish in all of the lakes examined. In addition, isoelectric focusing gels revealed variation in the haemoglobin protein components found in whitefish. Dwarf and normal whitefish exhibited a similar number of protein components, but the composition of these components differed, with dwarf whitefish bearing a greater proportion of cathodic components compared to the normals. Furthermore, dwarf whitefish showed significant haemoglobin gene upregulation in the brain compared with the levels shown in normals. Together, our results indicate that metabolic traits involved in oxygen transport differ between the whitefish ecotypes and the strong parallel patterns of divergence observed across lakes implicates ecologically driven selection pressures. We discuss the function of these traits in relation to the differing trophic niches occupied by the whitefish and the potential contributions of trait plasticity and genetic divergence to energetic adaptation.     

Ecological speciation in postglacial European whitefish: rapid adaptive radiations into the littoral,pelagic, and profundal lake habitats

Understanding how a monophyletic lineage of a species diverges into several adaptive forms has received increased attention in recent years, but the underlying mechanisms in this process are still under debate. Postglacial fishes are excellent model organisms for exploring this process, especially the initial stages of ecological speciation, as postglacial lakes represent replicated discrete environments with variation in available niches. Here, we combine data of niche utilization, trophic morphology, and 17 microsatellite loci to investigate the diversification process of three sympatric European whitefish morphs from three northern Fennoscandian lakes. The morphological divergence in the gill raker number among the whitefish morphs was related to the utilization of different trophic niches and was associated with reproductive isolation within and across lakes. The intralacustrine comparison of whitefish morphs showed that these systems represent two levels of adaptive divergence: (1) a consistent littoral–pelagic resource axis; and (2) a more variable littoral–profundal resource axis. The results also indicate that the profundal whitefish morph has diverged repeatedly from the ancestral littoral whitefish morph in sympatry in two different watercourses. In contrast, all the analyses performed revealed clustering of the pelagic whitefish morphs across lakes suggesting parallel postglacial immigration with the littoral whitefish morph into each lake. Finally, the analyses strongly suggested that the trophic adaptive trait, number of gill rakers, was under diversifying selection in the different whitefish morphs. Together, the results support a complex evolutionary scenario where ecological speciation acts, but where both allopatric (colonization history) and sympatric (within watercourse divergence) processes are involved.