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.     

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.     

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.     

The genetic architecture of ecological speciation and the association with signatures of selection in natural lake whitefish (Coregonus sp. Salmonidae) species pairs

Adaptive evolutionary change is contingent on variation and selection; thus, understanding adaptive divergence and ultimately speciation requires information on both the genetic basis of adaptive traits as well as an understanding of the role of divergent natural selection on those traits. The lake whitefish (Coregonus clupeaformis) consists of several sympatric "dwarf" (limnetic) and normal (benthic) species pairs that co-inhabit northern postglacial lakes. These young species pairs have evolved independently and display parallelism in life history, behavioral, and morphological divergence associated with the use of distinct trophic resources. We identified phenotype-environment associations and determined the genetic architecture and the role of selection modulating population genetic divergence in sympatric dwarf and normal lake whitefish. The genetic architecture of 9 adaptive traits was analyzed in 2 hybrid backcrosses individually phenotyped throughout their life history. Significant quantitative trait loci (QTL) were associated with swimming behavior (habitat selection and predator avoidance), growth rate, morphology (condition factor and gill rakers), and life history (onset of maturity and fecundity). Genome scans among 4 natural sympatric pairs, using loci segregating in the map, revealed a signature of selection for 24 loci. Loci exhibiting a signature of selection were associated with QTL relative to other regions of the genome more often than expected by chance alone. Two parallel QTL outliers for growth and condition factor exhibited segregation distortion in both mapping families, supporting the hypothesis that adaptive divergence contributing to parallel reductions of gene flow among natural populations may cause genetic incompatibilities. Overall, these findings offer evidence that the genetic architecture of ecological speciation is associated with signatures of selection in nature, providing strong support for the hypothesis that divergent natural selection is currently maintaining adaptive differentiation and promoting ecological speciation in lake whitefish species pairs.     

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.     

Contrasting patterns of mitochondrial DNA and microsatellite introgressive hybridization between lineages of lake whitefish (Coregonus clupeaformis); relevance for speciation

We performed a combined analysis of mitochondrial DNA (mtDNA) and microsatellite loci among lake whitefish (Coregonus clupeaformis) populations in order to assess the levels of congruence between both types of markers in defining patterns of genetic structuring, introgressive hybridization and inferring population origins in the hybrid zone of the St. John River basin. A second objective was to test the hypothesis that secondary contact between glacial lineages always resulted in the occurrence of sympatric dwarf and normal whitefish ecotypes. Fish were sampled from 35 populations and polymorphism was screened at mtDNA and six microsatellite loci for a total of 688 and 763 whitefish, respectively. Four lakes harbouring a single whitefish population of normal ecotype admixed with mtDNA haplotypes of different lineages were found. This confirmed that secondary contact between whitefish evolutionary lineages did not always result in the persistence of reproductively isolated ecotypes. Microsatellites further supported the definition of distinct glacial lineages by identifying lineage-specific allelic size groups. They also further supported the hypothesis that ecotypes originated from either a single founding lineage (sympatric divergence) or following secondary contacts between lineages (allopatric divergence), depending on the lake. In general, however, the pattern of population differentiation and introgressive hybridization observed at microsatellites was in sharp contrast with that depicted by mtDNA variation. Both factorial correspondence analysis and analysis of admixture proportion revealed a much more pronounced pattern of introgressive hybridization than depicted by mtDNA analyses. Variable levels of introgression indicated that environmental differences may be as important as the historical contingency of secondary contact in explaining the persistence of sympatric ecotypes and the differential pattern of introgressive hybridization among lakes. Whitefish populations from the St. John River basin hybrid zone represent a rare illustration of a continuum of both morphological and genetic differentiation within a given taxon, spanning from complete introgression to possibly complete reproductive isolation, depending on lakes. Thus, each lake may be viewed as a different temporal snapshot taken throughout the gradual process 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.     

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.     

Molecular population divergence and sexual selection on morphology in the banded demoiselle (Calopteryx splendens)

The importance of sexual selection in population divergence is of much interest, mainly because it is thought to cause reproductive isolation and hence could lead to speciation. Sexually selected traits have been hypothesized to diverge faster between populations than other traits, presumably because of differences in the strength, mechanism or dynamics of selection. We investigated this by quantifying population divergence in eight morphological characters in 12 south Swedish populations of a sexually dimorphic damselfly, the banded demoiselle (Calopteryx splendens). The morphological characters included a secondary sexual character, the male melanized wing spot, which has an important function in both inter- and intrasexual selection. In addition, we investigated molecular population divergence, revealed by amplified fragment length polymorphism (AFLP) analysis. Molecular population divergence was highly significant among these Northern European populations (overall F(st)=0.054; pairwise population F(st)'s ranged from approximately 0 to 0.13). We found evidence for isolation-by-distance (r=0.70) for the molecular markers and a significant correlation between molecular and phenotypic population divergence (r=0.39). One interpretation is that population divergence for the AFLP loci are affected by genetic drift, but is also indirectly influenced by selection, due to linkage with loci for the phenotypic traits. Field estimates of sexual and natural selection from two of the populations revealed fairly strong sexual selection on wing spot length, indicating that this trait has the potential to rapidly diverge, provided that variation is heritable and the observed selection is chronic.     

The role of local ecology during hybridization at the initial stages of ecological speciation in a marine snail

Hybrid zones of ecologically divergent populations are ideal systems to study the interaction between natural selection and gene flow during the initial stages of speciation. Here, we perform an amplified fragment length polymorphism (AFLP) genome scan in parallel hybrid zones between divergent ecotypes of the marine snail Littorina saxatilis, which is considered a model case for the study of ecological speciation. Ridged‐Banded (RB) and Smooth‐Unbanded (SU) ecotypes are adapted to different shore levels and microhabitats, although they present a sympatric distribution at the mid‐shore where they meet and mate (partially assortatively). We used shell morphology, outlier and nonoutlier AFLP loci from RB, SU and hybrid specimens captured in sympatry to determine the level of phenotypic and genetic introgression. We found different levels of introgression at parallel hybrid zones and nonoutlier loci showed more gene flow with greater phenotypic introgression. These results were independent from the phylogeography of the studied populations, but not from the local ecological conditions. Genetic variation at outlier loci was highly correlated with phenotypic variation. In addition, we used the relationship between genetic and phenotypic variation to estimate the heritability of morphological traits and to identify potential Quantitative Trait Loci to be confirmed in future crosses. These results suggest that ecology (exogenous selection) plays an important role in this hybrid zone. Thus, ecologically based divergent natural selection is responsible, simultaneously, for both ecotype divergence and hybridization. On the other hand, genetic introgression occurs only at neutral loci (nonoutliers). In the future, genome‐wide studies and controlled crosses would give more valuable information about this process of speciation in the face of gene flow.     

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.     

The genetic basis of intrinsic and extrinsic post-zygotic reproductive isolation jointly promoting speciation in the lake whitefish species complex (Coregonus clupeaformis)

Understanding the genetic architecture of reproductive barriers and the evolutionary forces that drove their divergence represents a considerable challenge towards understanding speciation. The objective of this study was to determine the genetic basis of intrinsic and extrinsic post-zygotic isolation in diverging populations of dwarf and normal lake whitefish with allopatric glacial origins. We found that the rate of embryonic mortality was 5.3-6.5 times higher in dwarf-normal hybrid backcrosses during development than in F1 dwarf and normal crosses. When comparing embryos that died during development against larvae that successfully hatched, patterns of Mendelian segregation at 101 loci whose linkage is known identified 13 loci distributed over seven linkage groups that exhibited significant shifts in segregation ratios leading to significant segregation distortion at these loci in the surviving progeny. Controlled crosses and quantitative trait loci analysis revealed a significant genetic basis for developmental time until emergence, a trait critical to fish larval survival in nature. Hatching backcross progeny exhibited asynchronous emergence and transgressive segregation, suggesting that extrinsic post-zygotic isolation may select against hybridization in specific environmental contexts. Evidence of a genetic basis for increased embryonic mortality followed by asynchronous emergence indicated that intrinsic and extrinsic mechanisms are not mutually exclusive in the formation and maintenance of reproductive isolation, but may be jointly promoting population divergence and ultimately speciation.     

Parasite diversity,patterns of MHC II variation and olfactory based mate choice in diverging three-spined stickleback ecotypes

Ecological speciation has been the subject of intense research in evolutionary biology but the genetic basis of the actual mechanism driving reproductive isolation has rarely been identified. The extreme polymorphism of the major histocompatibility complex (MHC), probably maintained by parasite-mediated selection, has been proposed as a potential driver of population divergence. We performed an integrative field and experimental study using three-spined stickleback river and lake ecotypes. We characterized their parasite load and variation at MHC class II loci. Fish from lakes and rivers harbor contrasting parasite communities and populations possess different MHC allele pools that could be the result of a combined action of genetic drift and parasite-mediated selection. We show that individual MHC class II diversity varies among populations and is lower in river ecotypes. Our results suggest the action of homogenizing selection within habitat type and diverging selection between habitat types. Finally, reproductive isolation was suggested by experimental evidence: in a flow channel design females preferred assortatively the odor of their sympatric male. This demonstrates the role of olfactory cues in maintaining reproductive isolation between diverging fish ecotypes.     

Linkage maps of the dwarf and Normal lake whitefish (Coregonus clupeaformis) species complex and their hybrids reveal the genetic architecture of population divergence

Elucidating the genetic architecture of population divergence may reveal the evolution of reproductive barriers and the genomic regions implicated in the process. We assembled genetic linkage maps for the dwarf and Normal lake whitefish species complex and their hybrids. A total of 877 AFLP loci and 30 microsatellites were positioned. The homology of mapped loci between families supported the existence of 34 linkage groups (of 40n expected) exhibiting 83% colinearity among linked loci between these two families. Classes of AFLP markers were not randomly distributed among linkage groups. Both AFLP and microsatellites exhibited deviations from Mendelian expectations, with 30.4% exhibiting significant segregation distortion across 28 linkage groups of the four linkage maps in both families (P < 0.00001). Eight loci distributed over seven homologous linkage groups were significantly distorted in both families and the level of distortion, when comparing homologous loci of the same phase between families, was correlated (Spearman R = 0.378, P = 0.0021). These results suggest that substantial divergence incurred during allopatric glacial separation and subsequent sympatric ecological specialization has resulted in several genomic regions that are no longer complementary between dwarf and Normal populations issued from different evolutionary glacial lineages.     

Selection on hybrids of ecologically divergent ecotypes of a marine snail: the relative importance of exogenous and endogenous barriers

Unravelling the form of selection acting on hybrids of ecotypes undergoing ecological speciation is essential to understand the mechanisms behind the evolution of reproductive isolation in the face of gene flow. Shell phenotype is known to be affected by natural selection and is involved in the fitness of the marine snail Littorina saxatilis. Here, we studied the association between shell traits and fitness in hybrids in order to determine the relative role of exogenous and endogenous selection in this hybrid zone of L. saxatilis. We show that directional selection is the predominant mode of selection among hybrids. We also show its heterogeneity, affecting different shell traits, within populations at the level of the microhabitat. Therefore, endogenous selection mechanisms are most probably lacking in this hybrid zone and exogenous barriers (pre‐ and post‐zygotic) are possibly one of the main forces behind the evolution of barriers to gene flow between these ecologically divergent ecotypes. This study shows how this barrier might represent an important type of reproductive isolation within ecological speciation, and this should be taken into account in future studies of speciation in hybrid zones. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 391–400.     

Genetic differentiation among European whitefish ecotypes based on microsatellite data

The amount of genetic differentiation at DNA microsatellite loci in European whitefish (Coregonus lavaretus) was assessed among ecotypes, populations and run-timing types. The magnitude of genetic changes potentially caused by hatchery broodstock rearing were also compared with those observed in corresponding natural populations. A total of 35 populations were studied, including 33 Coregonus lavaretus populations and two samples of Coregonus peled. Five of the six whitefish ecotypes in Finland were represented within C. lavaretus populations. Genetic diversity among C. lavaretus populations proved to be high compared to two C. peled populations. The genetic D(A) distance between these two species was as high as 0.86. The genetic differentiation among ecotypes was generally low and thus gives support for the hypothesis of one native European whitefish species in Scandinavia. Among the ecotypes the northern, large sparsely-rakered, bottom-dwelling whitefish was most unique. Thus, observed genetic differences in quantitative traits have either developed independently of phylogenetic lineages, or have mixed and later changed according to environments and selection pressures. Overall genetic distances between the anadromous whitefish populations along the Finnish coast, especially in the Bothnian Bay area, were small. Populations of this area have been heavily influenced by human activities, and they also have the highest probability of mixing by natural means. In two cases, the Rivers Iijoki and Tornionjoki, statistically significant genetic differences could be observed between summer- and autumn-run spawning-time types. Wild populations had slightly higher allelic diversity than hatchery-reared populations of corresponding rivers. Although some reduction in genetic diversity during hatchery rearing is possible, it is an important aid in maintaining endangered populations.     

昆虫同域物种形成机制及检验

从生物学、生态和遗传的角度阐述昆虫同域物种形成过程中涉及到的可能性机制。昆虫同域种的分化与作用于同域初始种群的歧化选择密切相关,歧化选择间接导致种群生态特征和遗传特征的分化,促进同域近缘种群间的生殖隔离。同域物种形成的过程中涉及到性状替换、性选择、同型交配等机制。寄主专化型多见于昆虫同域种的分化过程中,一般以植食性昆虫为主。有关昆虫同域物种形成的检验机制有多种,归纳起来主要包括同型交配的检验、遗传漂流的量化、遗传分化程度和连锁不平衡(LD)的检测、杂交后代适合度的估算等。目前发现在许多昆虫种类中存在同域物种形成的可能性,但是有关其隔离机制并没有得到充分的解释。     

Quantitative genetic inheritance of morphological divergence in a lake-stream stickleback ecotype pair: implications for reproductive isolation

Ecological selection against hybrids between populations occupying different habitats might be an important component of reproductive isolation during the initial stages of speciation. The strength and directionality of this barrier to gene flow depends on the genetic architecture underlying divergence in ecologically relevant phenotypes. We here present line cross analyses of inheritance for two key foraging-related morphological traits involved in adaptive divergence between stickleback ecotypes residing parapatrically in lake and stream habitats within the Misty Lake watershed (Vancouver Island, Canada). One main finding is the striking genetic dominance of the lake phenotype for body depth. Selection associated with this phenotype against first- and later-generation hybrids should therefore be asymmetric, hindering introgression from the lake to the stream population but not vice versa. Another main finding is that divergence in gill raker number is inherited additively and should therefore contribute symmetrically to reproductive isolation. Our study suggests that traits involved in adaptation might contribute to reproductive isolation qualitatively differently, depending on their mode of inheritance.