Investigating the extent of parallelism in morphological and genomic divergence among lake trout ecotypes in Lake Superior

Understanding the emergence of species through the process of ecological speciation is a central question in evolutionary biology which also has implications for conservation and management. Lake trout (Salvelinus namaycush) is renowned for the occurrence of different ecotypes linked to resource and habitat use throughout North America. We aimed to unravel the fine genetic structure of the four lake trout ecotypes in Lake Superior. A total of 486 individuals from four sites were genotyped at 6822 filtered SNPs using RADseq technology. Our results revealed different extent of morphological and genetic differentiation within the different sites. Overall, genetic differentiation was weak but significant and was on average three times higher between sites (mean F_ST = 0.016) than between ecotypes within sites (mean F_ST = 0.005) indicating higher level of gene flow or a more recent shared ancestor between ecotypes within each site than between populations of the same ecotype. Evidence of divergent selection was also found between ecotypes and/or in association with morphological variation. Outlier loci found in genes related to lipid metabolism and visual acuity were of particular interest in this context of ecotypic divergence. However, we did not find clear indication of parallelism at the genomic level, despite the presence of phenotypic parallelism among some ecotypes from different sampling sites. Overall, the occurrence of different levels of both genomic and phenotypic differentiation between ecotypes within each site with several differentiated loci linked to relevant biological functions supports the presence of a continuum of divergence in lake trout.     

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.     

Rapid antagonistic coevolution between primary and secondary sexual characters in horned beetles

Different structures may compete during development for a shared and limited pool of resources to sustain growth and differentiation. The resulting resource allocation trade-offs have the potential to alter both ontogenetic outcomes and evolutionary trajectories. However, little is known about the evolutionary causes and consequences of resource allocation trade-offs in natural populations. Here, we explore the significance of resource allocation trade-offs between primary and secondary sexual traits in shaping early morphological divergences between four recently separated populations of the horned beetle Onthophagus taurus as well as macroevolutionary divergence patterns across 10 Onthophagus species. We show that resource allocation trade-offs leave a strong signature in morphological divergence patterns both within and between species. Furthermore, our results suggest that genital divergence may, under certain circumstances, occur as a byproduct of evolutionary changes in secondary sexual traits. Given the importance of copulatory organ morphology for reproductive isolation our findings begin to raise the possibility that secondary sexual trait evolution may promote speciation as a byproduct. We discuss the implications of our results on the causes and consequences of resource allocation trade-offs in insects.     

The association of feeding behaviour with the resistance and tolerance to parasites in recently diverged sticklebacks

Divergent natural selection regimes can contribute to adaptive population divergence, but can be sensitive to human‐mediated environmental change. Nutrient loading of aquatic ecosystems, for example, might modify selection pressures by altering the abundance and distribution of resources and the prevalence and infectivity of parasites. Here, we used a mesocosm experiment to test for interactive effects of nutrient loading and parasitism on host condition and feeding ecology. Specifically, we investigated whether the common fish parasite Gyrodactylus sp. differentially affected recently diverged lake and stream ecotypes of three‐spined stickleback (Gasterosteus aculeatus). We found that the stream ecotype had a higher resistance to Gyrodactylus sp. infections than the lake ecotype, and that both ecotypes experienced a cost of parasitism, indicated by negative relationships between parasite load and both stomach fullness and body condition. Overall, our results suggest that in the early stages of adaptive population divergence of hosts, parasites can affect host resistance, body condition and diet.     

Mechanisms of a locally adaptive shift in allocation among growth,reproduction, and herbivore resistance in Mimulus guttatus*

Environmental gradients can drive adaptive evolutionary shifts in plant resource allocation among growth, reproduction, and herbivore resistance. However, few studies have attempted to connect these adaptations to underlying physiological and genetic mechanisms. Here, we evaluate potential mechanisms responsible for a coordinated locally adaptive shift between growth, reproduction, and herbivore defense in the yellow monkeyflower, Mimulus guttatus. Through manipulative laboratory experiments, we found that gibberellin (GA) growth hormones may play a role in the developmental divergence between perennial and annual ecotypes of M. guttatus. Further, we detected an interaction between a locally adaptive chromosomal inversion, DIV1, and GA addition. This finding is consistent with the inversion contributing to the evolutionary divergence between inland annual and coastal perennial ecotypes by reducing GA biosynthesis/activity in perennials. Finally, we found evidence that the DIV1 inversion is partially responsible for a coordinated shift in the divergence of growth, reproduction, and herbivore resistance traits between coastal perennial and inland annual M. guttatus. The inversion has already been established to have a substantial impact on the life‐history shift between long‐term growth and rapid reproduction. Here, we demonstrate that the DIV1 inversion also has sizable impacts on both the total abundance and composition of phytochemical compounds involved in herbivore resistance.     

A cascade of evolutionary change alters consumer-resource dynamics and ecosystem function

It is becoming increasingly clear that intraspecific evolutionary divergence influences the properties of populations, communities and ecosystems. The different ecological impacts of phenotypes and genotypes may alter selection on many species and promote a cascade of ecological and evolutionary change throughout the food web. Theory predicts that evolutionary interactions across trophic levels may contribute to hypothesized feedbacks between ecology and evolution. However, the importance of 'cascading evolutionary change' in a natural setting is unknown. In lakes in Connecticut, USA, variation in migratory behaviour and feeding morphology of a fish predator, the alewife (Alosa pseudoharengus), drives life-history evolution in a species of zooplankton prey (Daphnia ambigua). Here we evaluated the reciprocal impacts of Daphnia evolution on ecological processes in laboratory mesocosms. We show that life-history evolution in Daphnia facilitates divergence in rates of population growth, which in turn significantly alters consumer-resource dynamics and ecosystem function. These experimental results parallel trends observed in lakes. Such results argue that a cascade of evolutionary change, which has occurred over contemporary timescales, alters community and ecosystem processes.     

Parallel evolution of lake whitefish dwarf ecotypes in association with limnological features of their adaptive landscape

Sympatric fish populations observed in many north temperate lakes are among the best models to study the processes of population divergence and adaptive radiation. Despite considerable research on such systems, little is known about the associations between ecological conditions and the extent of ecotypic divergence. In this study, we examined the biotic and abiotic properties of postglacial lakes in which lake whitefish, Coregonus clupeaformis, occur as a derived dwarf ecotype in sympatry with an ancestral normal ecotype. We compared 19 limnological variables between two groups of lakes known from previous studies to harbour sympatric dwarf and normal ecotypes with high and low levels of phenotypic and genetic differentiation respectively. We found clear environmental differences between the two lake groups. Namely, oxygen was the most discriminant variable, where lakes harbouring the most divergent populations were characterized by the greatest hypolimnetic oxygen depletion. These lakes also had lower zooplankton densities and a narrower distribution of zooplantonic prey length. These results suggest that the highest differentiation between sympatric ecotypes occurs in lakes with reduced habitat and prey availability that could increase competition for resources. This in turns supports the hypothesis that parallelism in the extent of phenotypic divergence among sympatric whitefish ecotypes is associated with parallelism in adaptive landscape in terms of differences in limnological characteristics, as well as availability and structure of the zooplanktonic community.     

Lessons learned from the dog genome

Extensive genetic resources and a high-quality genome sequence position the dog as an important model species for understanding genome evolution, population genetics and genes underlying complex phenotypic traits. Newly developed genomic resources have expanded our understanding of canine evolutionary history and dog origins. Domestication involved genetic contributions from multiple populations of gray wolves probably through backcrossing. More recently, the advent of controlled breeding practices has segregated genetic variability into distinct dog breeds that possess specific phenotypic traits. Consequently, genome-wide association and selective sweep scans now allow the discovery of genes underlying breed-specific characteristics. The dog is finally emerging as a novel resource for studying the genetic basis of complex traits, including behavior.     

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.     

BAC library construction, screening and clone sequencing of lake whitefish (Coregonus clupeaformis, Salmonidae) towards the elucidation of adaptive species divergence

Genomic DNA sequences and other genomic resources are essential towards the elucidation of the genomic bases of adaptive divergence and reproductive isolation. Here, we describe the construction, characterization and screening of a nonarrayed BAC library for lake whitefish (Coregonus clupeaformis). We then show how the combined use of BAC library screening and next-generation sequencing can lead to efficient full-length assembly of candidate genes. The lake whitefish BAC library consists of 181,050 clones derived from a single heterozygous fish. The mean insert size is 92 Kb, representing 5.2 haploid genome equivalents. Ten BAC clones were isolated following a quantitative real-time PCR screening approach that targeted five previously identified candidate genes. Sequencing of these clones on a 454 GS FLX system yielded 178,000 reads with a mean length of 358 bp, for a total of 63.8 Mb. De novo assembly and annotation then allowed retrieval of contigs corresponding to each candidate gene, which also contained up- and/or downstream noncoding sequences. These results suggest that the lake whitefish BAC library combined with next-generation sequencing technologies will be key resources to achieve a better understanding of both adaptive divergence and reproductive isolation in lake whitefish species pairs as well as salmonid evolution in general.     

Flexible mate choice may contribute to ecotype assortative mating in pumpkinseed sunfish (Lepomis gibbosus)

Gene flow is expected to limit adaptive divergence, but the ecological and behavioural factors that govern gene flow are still poorly understood, particularly at the earliest stages of population divergence. Reduced gene flow through mate choice (sexual isolation) can evolve even under conditions of subtle population divergence if intermediate phenotypes have reduced fitness. We indirectly tested the hypothesis that mate choice has evolved between coexisting littoral and pelagic ecotypes of polyphenic pumpkinseed sunfish (Lepomis gibbosus) that have diverged in morphology and resource use and where intermediate phenotypes have reduced performance. We assessed the ecotype of nesting males and females using stable isotope estimates of diet and a divergent male morphological trait, oral jaw width. We found positive assortative mating between ecotypes in a common spawning habitat along exposed lake shorelines, but contrary to expectations, assortative mating was variably expressed between two sampling years. Although the factors that influence variable assortative mating remain unclear, our results are consistent with mate choice being expressed by ecotypes. Despite being variably expressed, mate choice will reduce gene flow between ecotypes and could contribute to further adaptive divergence depending on its frequency and strength in the population. Our findings add to a growing body of evidence indicating mate choice behaviour can be a plastic trait, an idea that should be more explicitly considered in empirical studies of mate choice as well as conceptual frameworks of mate choice evolution and adaptive divergence.     

Insights into the role of differential gene expression on the ecological adaptation of the snail <Emphasis Type="Italic">Littorina saxatilis</Emphasis>

Background  

In the past 40 years, there has been increasing acceptance that variation in levels of gene expression represents a major source of evolutionary novelty. Gene expression divergence is therefore likely to be involved in the emergence of incipient species, namely, in a context of adaptive radiation. In this study, a genome-wide expression profiling approach (cDNA-AFLP), validated by quantitative real-time polymerase chain reaction (qPCR) were used to get insights into the role of differential gene expression on the ecological adaptation of the marine snail Littorina saxatilis. This gastropod displays two sympatric ecotypes (RB and SU) which are becoming one of the best studied systems for ecological speciation.     

Effects of soil nitrogen on diploid advantage in fireweed,Chamerion angustifolium (Onagraceae)

In many ecosystems, plant growth and reproduction are nitrogen limited. Current and predicted increases of global reactive nitrogen could alter the ecological and evolutionary trajectories of plant populations. Nitrogen is a major component of nucleic acids and cell structures, and it has been predicted that organisms with larger genomes should require more nitrogen for growth and reproduction and be more negatively affected by nitrogen scarcities than organisms with smaller genomes. In a greenhouse experiment, we tested this hypothesis by examining whether the amount of soil nitrogen supplied differentially influenced the performance (fitness, growth, and resource allocation strategies) of diploid and autotetraploid fireweed (Chamerion angustifolium). We found that soil nitrogen levels differentially impacted cytotype performance, and in general, diploids were favored under low nitrogen conditions, but this diploid advantage disappeared under nitrogen enrichment. Specifically, when nitrogen was scarce, diploids produced more seeds and allocated more biomass toward seed production relative to investment in plant biomass or total plant nitrogen than did tetraploids. As nitrogen supplied increased, such discrepancies between cytotypes disappeared. We also found that cytotype resource allocation strategies were differentially dependent on soil nitrogen, and that whereas diploids adopted resource allocation strategies that favored current season reproduction when nitrogen was limiting and future reproduction when nitrogen was more plentiful, tetraploids adopted resource allocation strategies that favored current season reproduction under nitrogen enrichment. Together these results suggest nitrogen enrichment could differentially affect cytotype performance, which could have implications for cytotypes’ ecological and evolutionary dynamics under a globally changing climate.     

Genomics of Divergence along a Continuum of Parapatric Population Differentiation

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

Mortality dynamics of Daphnia in contrasting habitats and their role in ecological divergence

1. The Daphnia pulex‐pulicaria species complex has been proposed as an example of rapid ecological speciation, associated with divergence along the gradient of waterbody size from temporary ponds to deep, stratified lakes. However, this divergence is incomplete, and thus represents an opportunity to study ecological divergence as it is occurring. 2. Dynamics of twelve populations of Daphnia in the pulex‐pulicaria species complex were monitored over 1 year. Six temporary pond populations and six permanent lake populations were compared to evaluate demographic differences that may contribute to ecological divergence in this complex. 3. Pond populations experienced greater changes in density, which were reflected in more extreme growth rates, higher birth rates and higher mortality rates than those of lake populations. 4. Mitochondrial DNA was isolated from up to three clones from each population, the D‐loop of the control region was sequenced, and a phylogenetic tree was constructed. This tree revealed two strongly supported clades. The clades were not congruent with habitat type and nominal status, indicating that interhabitat gene flow occurs easily and that the nominal taxa are incompletely diverged. 5. Published reports of genetic life history differences in the D. pulex‐pulicaria complex are consistent with the demographic differences reported here. This suggests that ecological differences between the habitats are selectively maintaining trait differences despite the possibility for genetic exchange. Thus, these taxa may be at the inception of ecological speciation.     

Speciation with gene flow and the genetics of habitat transitions

Whether speciation can advance to completion in the face of initially high levels of gene flow is a very controversial topic in evolutionary biology. Extensive gene exchange is generally considered to homogenize populations and counteract divergence. Moreover, the role of introgressive hybridization in evolution remains largely unexplored in animals, particularly in freshwater zooplankton in which allopatric speciation is considered to be the norm. Our work investigates the genetic structure of two young ecological species: the pond species, Daphnia pulex and the lake species, Daphnia pulicaria. Phylogenetic and population genetics analyses were conducted on mitochondrial NADH dehydrogenase 5 (ND5) gene, the nuclear Lactate dehydrogenase (Ldh) gene and 21 nuclear microsatellite markers in 416 individuals from habitats with various degrees of permanence. The strong and consistent phylogenetic discordance between nuclear and mitochondrial markers suggests a complex evolutionary history of multiple independent habitat transition events that involved hybridization and introgression between lake and pond Daphnia. On the other hand, the low level of contemporary gene flow between adjacent populations indicates the presence of effective habitat isolating barriers. The Daphnia system provides strong evidence for a divergence-with-gene flow speciation model that involves multiple habitat transition events.     

The evolution of parent–offspring conflict over mate choice

In human societies, parents often have a strong influence on the mate choice of their offspring. Moreover, empirical studies show that conflict over mate choice between parents and offspring is widespread across human cultures. Here we provide the first theoretical investigation into this conflict, showing that it may result from an underlying evolutionary conflict over parental resource distribution. We present a series of evolutionary simulations in which we gradually expand a standard model of sexual selection by the stepwise addition of elements of parental involvement. In our model, females obtain resources enhancing their fecundity from both their chosen mate and their parents. Potential mates differ in their ability to provide resources and may signal this ability. Both females and their parents can develop a preference for the signal, with both preferences influencing the realized mate choice of the female. Parents may differentially allocate resources among their daughters depending on the resource-provisioning abilities of their sons-in-law. When fecundity returns on investment are diminishing, we find that parents invest most in daughters whose mates provide few resources. Subsequently, the daughters evolve to exploit this allocation rule through their mate choice, which is not in the parents' best interests. This results in a conflict over mate choice between parents and their offspring, manifested as an on-going divergence of offspring and parental preferences. We predict that the conflict should be most pronounced when fathers, as opposed to mothers, control resource allocation.     

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.     

PARALLEL EVOLUTION OF LAKE-STREAM PAIRS OF THREESPINE STICKLEBACKS (GASTEROSTEUS) INFERRED FROM MITOCHONDRIAL DNA VARIATION

Three drainage systems in British Columbia, Canada, contain divergent parapatric lake-stream pairs of threespine sticklebacks (Gasterosteus aculeatus): Drizzle and Mayer Lakes on Graham Island, Queen Charlotte Islands, and Misty Lake on northeastern Vancouver Island. Ecological and morphological differences between members of all three lake-stream pairs are strikingly similar; lake fish are melanistic and slim bodied with smaller mouths and more gill rakers than the mottled-brown and robust-bodied stream sticklebacks. We estimated the level of genetic divergence between lake and stream fish in Misty Lake and tested hypotheses of single versus multiple origins of the pairs by assaying mitochondrial DNA (mtDNA) restriction site variation in samples from the three lake systems. MtDNA analysis revealed the existence of two highly divergent lineages differing by 2.7% in sequence. One lineage predominated in Misty stream fish (73%), whereas the other lineage predominated in Misty Lake samples (96%). Comparable forms (lake or stream) in the different lakes did not cluster together in terms of mtDNA nucleotide divergence, suggesting that the pairs have had independent origins. We concluded that: (1) divergent mtDNA lineages in North Pacific sticklebacks stem from historical isolation in the two major glacial refugia proposed for the North Pacific (Beringia and Cascadia); (2) the stream and lake pair in Misty Lake are distinct gene pools; (3) the divergence between parapatric lake and stream Gasterosteus represents parallel evolution having occurred at least twice in the North Pacific; and (4) different scales of evolutionary divergence exist in North Pacific Gasterosteus, that is, a relatively ancient divergence of mtDNA clades as well as recent (i.e., postglacial) divergence of ecotypes within major clades.