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.     

Model‐based demographic inference of introgression history in European whitefish species pairs'

Parallel phenotypic differentiation is generally attributed to parallel adaptive divergence as an evolutionary response to similar environmental contrasts. Such parallelism may actually originate from several evolutionary scenarios ranging from repeated parallel divergence caused by divergent selection to a unique divergence event followed by gene flow. Reconstructing the evolutionary history underlying parallel phenotypic differentiation is thus fundamental to understand the relative contribution of demography and selection on genomic divergence during speciation. In this study, we investigate the divergence history of replicate European whitefish (Coregonus lavaretus), limnetic and benthic species pairs from two lakes in Norway and two lakes in Switzerland. Demographic models accounting for semi‐permeability and linked selection were fitted to the unfolded joint allele frequency spectrum built from genome‐wide SNPs and compared to each other in each species pair. We found strong support for a model of asymmetrical post‐glacial secondary contact between glacial lineages in all four lakes. Moreover, our results suggest that heterogeneous genomic differentiation has been shaped by the joint action of linked selection accelerating lineage sorting during allopatry, and heterogeneous migration eroding divergence at different rates along the genome following secondary contact. Our analyses reveal how the interplay between demography, selection and historical contingency has influenced the levels of diversity observed in previous whitefish phylogeographic studies. This study thus provides new insights into the historical demographic and selective processes that shaped the divergence associated with ecological speciation in European whitefish.     

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.     

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.     

SEXUAL SELECTION AGAINST HYBRIDS BETWEEN SYMPATRIC STICKLEBACK SPECIES: EVIDENCE FROM A FIELD EXPERIMENT

Sexual selection against viable, fertile hybrids may contribute to reproductive isolation between recently diverged species. If so, then sexual selection may be implicated in the speciation process. Laboratory measures of the mating success of hybrids may underestimate the amount of sexual selection against them if selection pressures are habitat specific. Male F₁ hybrids between sympatric benthic and limnetic sticklebacks (Gasterosteus aculeatus complex) do not suffer a mating disadvantage when tested in the laboratory. However, in the wild males choose different microhabitats and parental females tend to be found in the same habitats as conspecific males. This sets up the opportunity for sexual selection against male hybrids because they must compete with parental males for access to parental females. To test for sexual selection against adult F₁ hybrid males, we examined their mating success in enclosures in their preferred habitat (open, unvegetated substrate) where limnetic males and females also predominate. We found significantly reduced mating success in F₁ hybrid males compared with limnetic males. Thus, sexual selection, like other mechanisms of postzygotic isolation between young sister species, may be stronger in a wild setting than in the laboratory because of habitat-specific selection pressures. Our results are consistent with, but do not confirm, a role for sexual selection in stickleback speciation.     

EXPERIMENTAL CONFIRMATION THAT BODY SIZE DETERMINES MATE PREFERENCE VIA PHENOTYPE MATCHING IN A STICKLEBACK SPECIES PAIR

Mate choice by phenotype matching, whereby individuals prefer a mate whose phenotype is similar to their own, should facilitate speciation with gene flow. This is because the genes that control mate signal (the phenotype being matched) also determine the preferred mate signal (“mate preference”). Speciation is made even easier if phenotype matching is based on a trait under divergent natural selection. In this case, assortative mating should readily evolve as a byproduct of divergent selection on the trait. Previous observational studies of assortative mating between sympatric, hybridizing threespine stickleback species (Gasterosteus aculeatus complex) suggested that phenotype matching might occur by body size, a trait under divergent natural selection. To test this, we used experimental manipulation of body size to rule out the effects of confounding variables. We found that size‐manipulated benthic and limnetic stickleback females prefer mates whose body size more closely matches their own. It is thus likely that assortative mating by phenotype matching has facilitated the origin and persistence of benthic and limnetic threespine sticklebacks in the face of gene flow.     

Incipient sympatric speciation in Midas cichlid fish from the youngest and one of the smallest crater lakes in Nicaragua due to differential use of the benthic and limnetic habitats?

Understanding how speciation can occur without geographic isolation remains a central objective in evolutionary biology. Generally, some form of disruptive selection and assortative mating are necessary for sympatric speciation to occur. Disruptive selection can arise from intraspecific competition for resources. If this competition leads to the differential use of habitats and variation in relevant traits is genetically determined, then assortative mating can be an automatic consequence (i.e., habitat isolation). In this study, we caught Midas cichlid fish from the limnetic (middle of the lake) and benthic (shore) habitats of Crater Lake Asososca Managua to test whether some of the necessary conditions for sympatric speciation due to intraspecific competition and habitat isolation are given. Lake As. Managua is very small (<900 m in diameter), extremely young (maximally 1245 years of age), and completely isolated. It is inhabited by, probably, only a single endemic species of Midas cichlids, Amphilophus tolteca. We found that fish from the limnetic habitat were more elongated than fish collected from the benthic habitat, as would be predicted from ecomorphological considerations. Stable isotope analyses confirmed that the former also exhibit a more limnetic lifestyle than the latter. Furthermore, split‐brood design experiments in the laboratory suggest that phenotypic plasticity is unlikely to explain much of the observed differences in body elongation that we observed in the field. Yet, neutral markers (microsatellites) did not reveal any genetic clustering in the population. Interestingly, demographic inferences based on RAD‐seq data suggest that the apparent lack of genetic differentiation at neutral markers could simply be due to a lack of time, as intraspecific competition may only have begun a few hundred generations ago.     

A test of hybrid growth disadvantage in wild, free-ranging species pairs of threespine stickleback (Gasterosteus aculeatus) and its implications for ecological speciation

Ecological speciation is the evolution of reproductive isolation as a direct or indirect consequence of divergent natural selection. Reduced performance of hybrids in nature is thought to be an important process by which natural selection can favor the evolution of assortative mating and drive speciation. Benthic and limnetic sympatric species of threespine stickleback (Gasterosteus aculeatus) are adapted to alternative trophic niches (bottom browsing vs. open water planktivory, respectively) and reduced feeding performance of hybrids is thought to have contributed to the evolution of reproductive isolation. We tested this “hybrid‐disadvantage hypothesis” by inferring growth rates from otoliths sampled from wild, free‐ranging benthic, limnetic, and hybrid sticklebacks in two lakes. There were significant differences in growth rate between lakes, life‐history stages, and among years (maximum P = 0.02), as well as interactions between most factors, but not between hybrid and parental species sticklebacks in most comparisons. Our results provide little evidence of a growth disadvantage in hybrid sticklebacks when free‐ranging in nature. Although trophic ecology per se may contribute less to ecological speciation than envisioned, it may act in concert with other aspects of stickleback biology, such as interactions with parasites, predators, competitors, and/or sexual selection, to present strong multifarious selection against hybrids.     

Genomic architecture of ecologically divergent body shape in a pair of sympatric crater lake cichlid fishes

Determining the genetic bases of adaptations and their roles in speciation is a prominent issue in evolutionary biology. Cichlid fish species flocks are a prime example of recent rapid radiations, often associated with adaptive phenotypic divergence from a common ancestor within a short period of time. In several radiations of freshwater fishes, divergence in ecomorphological traits — including body shape, colour, lips and jaws — is thought to underlie their ecological differentiation, specialization and, ultimately, speciation. The Midas cichlid species complex (Amphilophus spp.) of Nicaragua provides one of the few known examples of sympatric speciation where species have rapidly evolved different but parallel morphologies in young crater lakes. This study identified significant QTL for body shape using SNPs generated via ddRAD sequencing and geometric morphometric analyses of a cross between two ecologically and morphologically divergent, sympatric cichlid species endemic to crater Lake Apoyo: an elongated limnetic species (Amphilophus zaliosus) and a high‐bodied benthic species (Amphilophus astorquii). A total of 453 genome‐wide informative SNPs were identified in 240 F₂ hybrids. These markers were used to construct a genetic map in which 25 linkage groups were resolved. Seventy‐two segregating SNPs were linked to 11 QTL. By annotating the two most highly supported QTL‐linked genomic regions, genes that might contribute to divergence in body shape along the benthic–limnetic axis in Midas cichlid sympatric adaptive radiations were identified. These results suggest that few genomic regions of large effect contribute to early stage divergence in Midas cichlids.     

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.     

REPRODUCTIVE CHARACTER DISPLACEMENT OF MALE STICKLEBACK MATE PREFERENCE: REINFORCEMENT OR DIRECT SELECTION?

Reproductive isolation can evolve between species as a byproduct of adaptation to different niches, through reinforcement, and by direct selection on mating preferences. We investigated the potential role of direct selection in the reproductive isolation between sympatric species of threespine sticklebacks. Each sympatric pair consists of a small-bodied limnetic species and large-bodied benthic species. We compared the mate preferences and courtship behavior of males from one sympatric limnetic population and two allopatric populations. We used limnetic-like allopatric populations to control for the effects of ecological character displacement and adaptation to different niches on mate preferences. The sympatric limnetic males preferred the small limnetic females, whereas the allopatric limnetic-like males preferred the large benthic females, suggesting that adaptation to the limnetic niche does not automatically confer a preference for small limnetic females. This reproductive character displacement of male preference is consistent with the predictions of both reinforcement and direct selection on mate preferences. To test for direct selection, we assessed a prediction of one proposed mechanism: predation by benthic females on eggs guarded by limnetic males. The allopatric males come from populations in which there is no egg predation. Sympatric limnetic males were more aggressive toward benthic females than toward limnetic females, whereas the allopatric limnetic-like males did not treat the two types of females differently. The contrast in male behavior suggests that egg predation has shaped male preferences. Direct selection is potentially more effective than indirect selection via reinforcement, and it is likely that it has been important in building up reproductive isolation between limnetic and benthic sticklebacks.     

An experimental test of the effect of predation upon behaviour and trait correlations in the threespine stickleback

Benthic and limnetic threespine stickleback (Gasterosteus aculeatus) are a classic example of ecological speciation. Behavioural and armour divergence between these species has been predicted to be the result of divergent selection driven in part by differential predation from cutthroat trout (Oncorhynchus clarki). To experimentally test this prediction, we reared split families of benthic–limnetic hybrids in the presence or absence of trout predation. Our results show that the presence of trout had little effect upon stickleback behaviour. We then compared performance in behavioural assays among stickleback that varied in armour to test if armour morphology correlates with behaviour. Our measurements revealed trait correlations between several behaviours and components of armour morphology. Trout predation did not result in an increased correlation between traits, therefore differential trout predation between benthics and limnetics is unlikely to be the cause of these correlations. The presence of trait correlations in advanced generation hybrids suggests that pleiotropy or linkage between genes underlying behaviour and armour morphology may be greater than previously appreciated.     

Adaptive plasticity and genetic divergence in feeding efficiency during parallel adaptive radiation of whitefish (Coregonus spp.)

Parallel phenotypic divergence in replicated adaptive radiations could either result from parallel genetic divergence in response to similar divergent selection regimes or from equivalent phenotypically plastic response to the repeated occurrence of contrasting environments. In post‐glacial fish, replicated divergence in phenotypes along the benthic‐limnetic habitat axis is commonly observed. Here, we use two benthic‐limnetic species pairs of whitefish from two Swiss lakes, raised in a common garden design, with reciprocal food treatments in one species pair, to experimentally measure whether feeding efficiency on benthic prey has a genetic basis or whether it underlies phenotypic plasticity (or both). To do so, we offered experimental fish mosquito larvae, partially burried in sand, and measured multiple feeding efficiency variables. Our results reveal both, genetic divergence as well as phenotypically plastic divergence in feeding efficiency, with the phenotypically benthic species raised on benthic food being the most efficient forager on benthic prey. This indicates that both, divergent natural selection on genetically heritable traits and adaptive phenotypic plasticity, are likely important mechanisms driving phenotypic divergence in adaptive radiation.     

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.     

Adaptation and acclimation of aerobic exercise physiology in Lake Whitefish ecotypes (Coregonus clupeaformis)

The physiological mechanisms underlying local adaptation in natural populations of animals, and whether the same mechanisms contribute to adaptation and acclimation, are largely unknown. Therefore, we tested for evolutionary divergence in aerobic exercise physiology in laboratory bred, size‐matched crosses of ancestral, benthic, normal Lake Whitefish (Coregonus clupeaformis) and derived, limnetic, more actively swimming “dwarf” ecotypes. We acclimated fish to constant swimming (emulating limnetic foraging) and control conditions (emulating normal activity levels) to simultaneously study phenotypic plasticity. We found extensive divergence between ecotypes: dwarf fish generally had constitutively higher values of traits related to oxygen transport (ventricle size) and use by skeletal muscle (percent oxidative muscle, mitochondrial content), and also evolved differential plasticity of mitochondrial function (Complex I activity and flux through Complexes I–IV and IV). The effects of swim training were less pronounced than differences among ecotypes and the traits which had a significant training effect (ventricle protein content, ventricle malate dehydrogenase activity, and muscle Complex V activity) did not differ among ecotypes. Only one trait, ventricle mass, varied in a similar manner with acclimation and adaptation and followed a pattern consistent with genetic accommodation. Overall, the physiological and biochemical mechanisms underlying acclimation and adaptation to swimming activity in Lake Whitefish differ.     

Evidence for nonallopatric speciation among closely related sympatric Heliotropium species in the Atacama Desert

The genetic structure of populations of closely related, sympatric species may hold the signature of the geographical mode of the speciation process. In fully allopatric speciation, it is expected that genetic differentiation between species is homogeneously distributed across the genome. In nonallopatric speciation, the genomes may remain undifferentiated to a large extent. In this article, we analyzed the genetic structure of five sympatric species from the plant genus Heliotropium in the Atacama Desert. We used amplified fragment length polymorphisms (AFLPs) to characterize the genetic structure of these species and evaluate their genetic differentiation as well as the number of loci subject to positive selection using divergence outlier analysis (DOA). The five species form distinguishable groups in the genetic space, with zones of overlap, indicating that they are possibly not completely isolated. Among‐species differentiation accounts for 35% of the total genetic differentiation (F_ST = 0.35), and F_ST between species pairs is positively correlated with phylogenetic distance. DOA suggests that few loci are subject to positive selection, which is in line with a scenario of nonallopatric speciation. These results support the idea that sympatric species of Heliotropium sect. Cochranea are under an ongoing speciation process, characterized by a fluctuation of population ranges in response to pulses of arid and humid periods during Quaternary times.     

Patterns and processes in the evolutionary history of parrotfishes (Family Labridae)

Phylogenetic reconstruction of the evolutionary relationships among 61 of the 70 species of the parrotfish genera Chlorurus and Scarus (Family Labridae) based on mitochondrial and nuclear gene sequences retrieved 15 well‐supported clades with mid Pliocene/Pleistocene diversification. Twenty‐two reciprocally monophyletic sister‐species pairs were identified: 64% were allopatric, and the remainder were sympatric. Age of divergence was similar for allopatric and sympatric species pairs. Sympatric sister pairs displayed greater divergence in morphology, ecology, and sexually dimorphic colour patterns than did allopatric pairs, suggesting that both genetic drift in allopatric species pairs and ecologically adaptive divergence between members of sympatric pairs have played a role in diversification. Basal species typically have small geographical ranges and are restricted to geographically and ecologically peripheral reef habitats. We found little evidence that a single dominant process has driven diversification, nor did we detect a pattern of discrete, sequential stages of diversification in relation to habitat, ecology, and reproductive biology. The evolution of Chlorurus and Scarus has been complex, involving a number of speciation processes. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 107 , 529–557.     

Ecological speciation in anemone‐associated snapping shrimps (Alpheus armatus species complex)

Divergent natural selection driven by competition for limited resources can promote speciation, even in the presence of gene flow. Reproductive isolation is more likely to result from divergent selection when the partitioned resource is closely linked to mating. Obligate symbiosis and host fidelity (mating on or near the host) can provide this link, creating ideal conditions for speciation in the absence of physical barriers to dispersal. Symbiotic organisms often experience competition for hosts, and host fidelity ensures that divergent selection for a specific host or host habitat can lead to speciation and strengthen pre‐existing reproductive barriers. Here, we present evidence that diversification of a sympatric species complex occurred despite the potential for gene flow and that partitioning of host resources (both by species and by host habitat) has contributed to this diversification. Four species of snapping shrimps (Alpheus armatus, A. immaculatus, A. polystictus and A. roquensis) are distributed mainly sympatrically in the Caribbean, while the fifth species (A. rudolphi) is restricted to Brazil. All five species are obligate commensals of sea anemones with a high degree of fidelity and ecological specificity for host species and habitat. We analysed sequence data from 10 nuclear genes and the mitochondrial COI gene in 11–16 individuals from each of the Caribbean taxa and from the only available specimen of the Brazilian taxon. Phylogenetic analyses support morphology‐based species assignments and a well‐supported Caribbean clade. The Brazilian A. rudolphi is recovered as an outgroup to the Caribbean taxa. Isolation–migration coalescent analysis provides evidence for historical gene flow among sympatric sister species. Our data suggest that both selection for a novel host and selection for host microhabitat may have promoted diversification of this complex despite gene flow.