From Local Adaptation to Ecological Speciation in Copepod Populations from Neighboring Lakes

Continental copepods have been derived from several independent invasive events from the sea, but the subsequent evolutionary processes that account for the current diversity in lacustrine environments are virtually unknown. Salinity is highly variable among lakes and constitutes a source of divergent selection driving potential reproductive isolation. We studied four populations of the calanoid copepod Leptodiaptomus cf. sicilis inhabiting four neighboring lakes with a common history (since the Late Pleistocene) located in the Oriental Basin, Mexico; one lake is shallow and varies in salinity periodically (1.4–10 g L^-1), while three are deep and permanent, with constant salinity (0.5, 1.1 and 6.5 g L^-1, respectively). We hypothesized that (1) these populations belong to a different species than L. sicilis sensu stricto and (2) are experiencing ecologically based divergence due to salinity differences. We assessed morphological and molecular (mtDNA) COI variation, as well as fitness differences and tests of reproductive isolation. Although relationships of the Mexican populations with L. sicilis s.s. could not be elucidated, we identified a clear pattern of divergent selection driven by salinity conditions. The four populations can still be considered a single biological species (sexual recognition and hybridization are still possible in laboratory conditions), but they have diverged into at least three different phenotypes: two locally adapted, specialized in the lakes of constant salinity (saline vs. freshwater), and an intermediate generalist phenotype inhabiting the temporary lake with fluctuating salinity. The specialized phenotypes are poorly suited as migrants, so prezygotic isolation due to immigrant inviability is highly probable. This implication was supported by molecular evidence that showed restricted gene flow, persistence of founder events, and a pattern of allopatric fragmentation. This study showed how ecologically based divergent selection may explain diversification patterns in lacustrine copepods.     

A gene with major phenotypic effects as a target for selection vs. homogenizing gene flow

Genes with major phenotypic effects facilitate quantifying the contribution of genetic vs. plastic effects to adaptive divergence. A classical example is Ectodysplasin (Eda), the major gene controlling lateral plate phenotype in three‐spined stickleback. Completely plated marine stickleback populations evolved repeatedly towards low‐plated freshwater populations, representing a prime example of parallel evolution by natural selection. However, many populations remain polymorphic for lateral plate number. Possible explanations for this polymorphism include relaxation of selection, disruptive selection or a balance between divergent selection and gene flow. We investigated 15 polymorphic stickleback populations from brackish and freshwater habitats in coastal North‐western Europe. At each site, we tracked changes in allele frequency at the Eda gene between subadults in fall, adults in spring and juveniles in summer. Eda genotypes were also compared for body size and reproductive investment. We observed a fitness advantage for the Eda allele for the low morph in freshwater and for the allele for the complete morph in brackish water. Despite these results, the differentiation at the Eda gene was poorly correlated with habitat characteristics. Neutral population structure was the best predictor of spatial variation in lateral plate number, suggestive of a substantial effect of gene flow. A meta‐analysis revealed that the signature of selection at Eda was weak compared to similar studies in stickleback. We conclude that a balance between divergent selection and gene flow can maintain stickleback populations polymorphic for lateral plate number and that ecologically relevant genes may not always contribute much to local adaptation, even when targeted by selection.     

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.     

Where the Lake Meets the Sea: Strong Reproductive Isolation Is Associated with Adaptive Divergence between Lake Resident and Anadromous Three-Spined Sticklebacks

Contact zones between divergent forms of the same species are often characterised by high levels of phenotypic diversity over small geographic distances. What processes are involved in generating such high phenotypic diversity? One possibility is that introgression and recombination between divergent forms in contact zones results in greater phenotypic and genetic polymorphism. Alternatively, strong reproductive isolation between forms may maintain distinct phenotypes, preventing homogenisation by gene flow. Contact zones between divergent freshwater-resident and anadromous stickleback (Gasterosteus aculeatus L.) forms are numerous and common throughout the species distribution, offering an opportunity to examine these contrasting hypotheses in greater detail. This study reports on an interesting new contact zone located in a tidally influenced lake catchment in western Ireland, characterised by high polymorphism for lateral plate phenotypes. Using neutral and QTL-linked microsatellite markers, we tested whether the high diversity observed in this contact zone arose as a result of introgression or reproductive isolation between divergent forms: we found strong support for the latter hypothesis. Three phenotypic and genetic clusters were identified, consistent with two divergent resident forms and a distinct anadromous completely plated population that migrates in and out of the system. Given the strong neutral differentiation detected between all three morphotypes (mean F_ST = 0.12), we hypothesised that divergent selection between forms maintains reproductive isolation. We found a correlation between neutral genetic and adaptive genetic differentiation that support this. While strong associations between QTL linked markers and phenotypes were also observed in this wild population, our results support the suggestion that such associations may be more complex in some Atlantic populations compared to those in the Pacific. These findings provide an important foundation for future work investigating the dynamics of gene flow and adaptive divergence in this newly discovered stickleback contact zone.     

Resource Competition and Coevolution in Sticklebacks

Threespine stickleback in young postglacial lakes provide a compelling example of coevolution between species that compete for resources. Coexisting pairs of stickleback species are highly divergent in habitat, diet, and body size and shape, whereas stickleback occurring alone in lakes are intermediate. We used experiments in ponds to test mechanisms of divergence between coexisting species. The results support the hypothesis of coevolution by resource competition between stickleback, but we found evidence that interactions with natural enemies also contribute to divergence. Natural selection arising from these interactions selects against intermediate phenotypes, included hybrids, and thus has contributed to the origin and persistence of stickleback species.     

Parallel and nonparallel aspects of ecological, phenotypic, and genetic divergence across replicate population pairs of lake and stream stickleback

Parallel (or convergent) evolution provides strong evidence for a deterministic role of natural selection: similar phenotypes evolve when independent populations colonize similar environments. In reality, however, independent populations in similar environments always show some differences: some nonparallel evolution is present. It is therefore important to explicitly quantify the parallel and nonparallel aspects of trait variation, and to investigate the ecological and genetic explanations for each. We performed such an analysis for threespine stickleback (Gasterosteus aculeatus) populations inhabiting lake and stream habitats in six independent watersheds. Morphological traits differed in the degree to which lake-stream divergence was parallel across watersheds. Some aspects of this variation were correlated with ecological variables related to diet, presumably reflecting the strength and specifics of divergent selection. Furthermore, a genetic scan revealed some markers that diverged between lakes and streams in many of the watersheds and some that diverged in only a few watersheds. Moreover, some of the lake-stream divergence in genetic markers was associated within some of the lake-stream divergence in morphological traits. Our results suggest that parallel evolution, and deviations from it, are primarily the result of natural selection, which corresponds in only some respects to the dichotomous habitat classifications frequently used in such studies.     

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.     

Phenotypic Plasticity and the Ecotypic Differentiation of Aggressive Behavior in Threespine Stickleback

The threespine stickleback fish, Gasterosteus aculeatus, has undergone a remarkable postglacial adaptive radiation in which an ancient oceanic ancestor has given rise to uncountable freshwater populations. The radiation is characterized by repeated, independent evolution of similar derived phenotypes under similar environmental conditions. A common pattern of divergence is caused by differences in habitat that favor morphological and behavioral features that enhance efficiency of feeding on plankton (limnetic ecotypes) vs. those that enhance efficiency of feeding on benthic invertebrates (benthic ecotypes). These two ecotypes exhibit consistently different patterns of courtship and of foraging and cannibalistic behavior (divergent behavioral syndromes). Here, we demonstrate that there also exist differences in aggression toward conspecifics that are likely to be characteristic of the ecotypes. We report differences in patterns of aggression toward rivals between the ecotypes and offer evidence of differences in the patterns of phenotypic plasticity (norms of reaction) for these traits across population types, and of differences in the incorporation of aggressive elements of behavior in courtship. These data support an earlier suggestion that differences in aggressive tendencies could have facilitated assortative mating between the four benthic–limnetic species pairs found in British Columbia lakes, and they demonstrate the need to evaluate divergent behavioral phenotypes in this radiation as phenotypic norms of reaction rather than as fixed traits.     

A comparison of the association between large haplotype blocks under selection and the presence/absence of inversions

Inversions may contribute to ecological adaptation and phenotypic diversity, and with the advent of “second” and “third” generation sequencing technologies, the ability to detect inversion polymorphisms has been enhanced dramatically. A key molecular consequence of an inversion is the suppression of recombination allowing independent accumulation of genetic changes between alleles over time. This may lead to the development of divergent haplotype blocks maintained by balancing selection. Thus, divergent haplotype blocks are often considered as indicating the presence of an inversion. In this paper, we first review the features of a 7.7 Mb inversion causing the Rose‐comb phenotype in chicken, as a model for how inversions evolve and directly affect phenotypes. Second, we compare the genetic basis for alternative mating strategies in ruff and timing of reproduction in Atlantic herring, both associated with divergent haplotype blocks. Alternative male mating strategies in ruff are associated with a 4.5 Mb inversion that occurred about 4 million years ago. In fact, the ruff inversion shares some striking features with the Rose‐comb inversion such as disruption of a gene at one of the inversion breakpoints and generation of a new allele by recombination between the inverted and noninverted alleles. In contrast, inversions do not appear to be a major reason for the fairly large haplotype blocks (range 10–200 kb) associated with ecological adaptation in the herring. Thus, it is important to note that divergent haplotypes may also be maintained by natural selection in the absence of structural variation.     

Sympatric and Allopatric Divergence of MHC Genes in Threespine Stickleback

Parasites can strongly affect the evolution of their hosts, but their effects on host diversification are less clear. In theory, contrasting parasite communities in different foraging habitats could generate divergent selection on hosts and promote ecological speciation. Immune systems are costly to maintain, adaptable, and an important component of individual fitness. As a result, immune system genes, such as those of the Major Histocompatability Complex (MHC), can change rapidly in response to parasite-mediated selection. In threespine stickleback (Gasterosteus aculeatus), as well as in other vertebrates, MHC genes have been linked with female mating preference, suggesting that divergent selection acting on MHC genes might influence speciation. Here, we examined genetic variation at MHC Class II loci of sticklebacks from two lakes with a limnetic and benthic species pair, and two lakes with a single species. In both lakes with species pairs, limnetics and benthics differed in their composition of MHC alleles, and limnetics had fewer MHC alleles per individual than benthics. Similar to the limnetics, the allopatric population with a pelagic phenotype had few MHC alleles per individual, suggesting a correlation between MHC genotype and foraging habitat. Using a simulation model we show that the diversity and composition of MHC alleles in a sympatric species pair depends on the amount of assortative mating and on the strength of parasite-mediated selection in adjacent foraging habitats. Our results indicate parallel divergence in the number of MHC alleles between sympatric stickleback species, possibly resulting from the contrasting parasite communities in littoral and pelagic habitats of lakes.     

Parallel evolution and ecological selection: replicated character displacement in spadefoot toads

Ecological character displacement—trait evolution stemming from selection to lessen resource competition between species—is most often inferred from a pattern in which species differ in resource-use traits in sympatry but not in allopatry, and in which sympatric populations within each species differ from conspecific allopatric populations. Yet, without information on population history, the presence of a divergent phenotype in multiple sympatric populations does not necessarily imply that there has been repeated evolution of character displacement. Instead, such a pattern may arise if there has been character displacement in a single ancestral population, followed by gene flow carrying the divergent phenotype into multiple, derived, sympatric populations. Here, we evaluate the likelihood of such historical events versus ongoing ecological selection in generating divergence in trophic morphology between multiple populations of spadefoot toad (Spea multiplicata) tadpoles that are in sympatry with a heterospecific and those that are in allopatry. We present both phylogenetic and population genetic evidence indicating that the same divergent trait, which minimizes resource competition with the heterospecific, has arisen independently in multiple sympatric populations. These data, therefore, provide strong indirect support for competition''s role in divergent trait evolution.     

Tameness correlates with domestication related traits in a Red Junglefowl intercross

Early animal domestication may have been driven by selection on tameness. Selection on only tameness can bring about correlated selection responses in other traits, not intentionally selected upon, which may be one cause of the domesticated phenotype. We predicted that genetically reduced fear towards humans in Red Junglefowl, ancestors of domesticated chickens, would be correlated to other traits included in the domesticated phenotype. Fear level was determined by a standardised behaviour test, where the reaction towards an approaching human was recorded. We first selected birds for eight generations for either high or low fear levels in this test, to create two divergent selection lines. An F3 intercross, with birds from the eighth generation as parentals, was generated to study correlations between fear-of-human scores and other unselected phenotypes, possibly caused by pleiotropy or linkage. Low fear-of-human scores were associated with higher body weight and growth rates, and with increased activity in an open field test, indicating less general fearfulness. In females, low fear-of-human scores were also associated with more efficient fear habituation and in males with an increased tendency to emit food calls in a mirror test, indicating increased social dominance. Low fear-of-human scores were also associated with smaller brain relative to body weight, and with larger cerebrum relative to total brain weight in females. All these effects are in line with the changes observed in domesticated chickens compared to their ancestors, and we conclude that tameness may have been a driving factor underlying some aspects of the domesticated phenotype.     

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.     

Correlated evolution of life-history with size at maturity in Daphnia pulicaria: patterns within and between populations

Explaining the repeated evolution of similar sets of traits under similar environmental conditions is an important issue in evolutionary biology. The extreme alternative classes of explanations for correlated suites of traits are optimal adaptation and genetic constraint resulting from pleiotropy. Adaptive explanations presume that individual traits are free to evolve to their local optima and that convergent evolution represents particularly adaptive combinations of traits. Alternatively, if pleiotropy is strong and difficult to break, strong selection on one or a few particularly important characters would be expected to result in consistent correlated evolution of associated traits. If pleiotropy is common, we predict that the pattern of divergence among populations will consistently reflect the within-population genetic architecture. To test the idea that the multivariate life-history phenotype is largely a byproduct of strong selection on body size, we imposed divergent artificial selection on size at maturity upon two populations of the cladoceran Daphnia pulicaria, chosen on the basis of their extreme divergence in body size. Overall, the trajectory of divergence between the two natural populations did not differ from that predicted by the genetic architecture within each population. However, the pattern of correlated responses suggested the presence of strong pleiotropic constraints only for adult body size and not for other life-history traits. One trait, offspring size, appears to have evolved in a way different from that expected from the within-population genetic architecture and may be under stabilizing selection.     

Divergent selection as revealed by P(ST) and QTL-based F(ST) in three-spined stickleback (Gasterosteus aculeatus) populations along a coastal-inland gradient

Three measures of divergence, estimated at nine putatively neutral microsatellite markers, 14 quantitative traits, and seven quantitative trait loci (QTL) were compared in eight populations of the three-spined stickleback (Gasterosteus aculeatus L.) living in the Scheldt river basin (Belgium). Lowland estuarine and polder populations were polymorphic for the number of lateral plates, whereas upland freshwater populations were low-plated. The number of short gill rakers and the length of dorsal and pelvic spines gradually declined along a coastal-inland gradient. Plate number, short gill rakers and spine length showed moderate to strong signals of divergent selection between lowland and upland populations in comparison between P(ST) (a phenotypic alternative for Q(ST)) and neutral F(ST). However, such comparisons rely on the unrealistic assumption that phenotypic variance equals additive genetic variance, and that nonadditive genetic effects and environmental effects can be minimized. In order to verify this assumption and to confirm the phenotypic signals of divergence, we tested for divergent selection at the underlying QTL. For plate number, strong genetic evidence for divergent selection between lowland and upland populations was obtained based on an intron marker of the Eda gene, of which the genotype was highly congruent with plate morph. Genetic evidence for divergent selection on short gill rakers was limited to some population pairs where F(ST) at only one of two QTL was detected as an outlier, although F(ST) at both loci correlated significantly with P(ST). No genetic confirmation was obtained for divergent selection on dorsal spine length, as no outlier F(ST)s were detected at dorsal spine QTL, and no significant correlations with P(ST) were observed.     

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.     

Genetic evidence for ecological divergence in kokanee salmon

The evolution of locally adapted phenotypes among populations that experience divergent selective pressures is a central mechanism for generating and maintaining biodiversity. Recently, the advent of high‐throughput DNA sequencing technology has provided tools for investigating the genetic basis of this process in natural populations of nonmodel organisms. Kokanee, the freshwater form of sockeye salmon (Oncorhynchus nerka), occurs as two reproductive ecotypes, which differ in spawning habitat (tributaries vs. shorelines); however, outside of the spawning season the two ecotypes co‐occur in many lakes and lack diagnostic morphological characteristics. We used restriction site‐associated DNA (RAD) sequencing to identify 6145 SNPs and genotype kokanee from multiple spawning sites in Okanagan Lake (British Columbia, Canada). Outlier tests revealed 18 loci putatively under divergent selection between ecotypes, all of which exhibited temporally stable allele frequencies within ecotypes. Six outliers were annotated to sequences in the NCBI database, two of which matched genes associated with early development. There was no evidence for neutral genetic differentiation; however, outlier loci demonstrated significant structure with respect to ecotype and had high assignment accuracy in mixed composition simulations. The absence of neutral structure combined with a small number of highly divergent outlier loci is consistent with theoretical predictions for the early stages of ecological divergence. These outlier loci were then applied to a realistic fisheries scenario in which additional RAD sequencing was used to genotype kokanee collected by trawl in Okanagan Lake, providing preliminary evidence that this approach may be an effective tool for conservation and management.     

Incipient speciation driven by hypertrophied lips in Midas cichlid fishes?

Sympatric speciation has been debated in evolutionary biology for decades. Although it has gained in acceptance recently, still only a handful of empirical examples are seen as valid (e.g. crater lake cichlids). In this study, we disentangle the role of hypertrophied lips in the repeated adaptive radiations of Nicaraguan crater lake cichlid fish. We assessed the role of disruptive selection and assortative mating during the early stages of divergence and found a functional trade‐off in feeding behaviour between thick‐ and thin‐lipped ecotypes, suggesting that this trait is a target of disruptive selection. Thick‐lipped fish perform better on nonevasive prey at the cost of a poorer performance on evasive prey. Using enclosures in the wild, we found that thick‐lipped fish perform significantly better in rocky than in sandy habitats. We found almost no mixed pairs during two breeding seasons and hence significant assortative mating. Genetic differentiation between ecotypes seems to be related to the time since colonization, being subtle in L. Masaya (1600 generations ago) and absent in the younger L. Apoyeque (<600 generations ago). Genome‐wide differentiation between ecotypes was higher in the old source lakes than in the young crater lakes. Our results suggest that hypertrophied lips might be promoting incipient sympatric speciation through divergent selection (ecological divergence in feeding performance) and nonrandom mating (assortative mating) in the young Nicaraguan crater lakes. Nonetheless, further manipulative experiments are needed in order to confirm the role of hypertrophied lips as the main cue for assortative mating.