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

Environmental correlates of adaptive diversification in postglacial freshwater fishes

Determining how environmental conditions contribute to divergence among populations and drive speciation is fundamental to resolving mechanisms and understanding outcomes in evolutionary biology. Postglacial freshwater fish species in the Northern Hemisphere are ideal biological systems to explore the effects of environment on diversification in morphology, ecology, and genetics (ecomorph divergences) within lakes. To date, various environmental factors have been implicated in the presence of multiple ecomorphs within particular lakes or regions. However, concerted evidence for generalizable patterns in environmental variables associated with speciation across geographical regions and across species and genera has been lacking. Here, we aimed to identify key biotic and abiotic factors associated with ecological divergence of postglacial freshwater fish species into multiple sympatric ecomorphs, focusing on species in the well-studied, widespread, and co-distributed genera Gasterosteus, Salvelinus, and Coregonus (stickleback, charr, and whitefish, respectively). We found that the presence of multiple sympatric ecomorphs tended to be associated with increasing lake surface area, maximum depth, and nutrient availability. In addition, predation, competition, and prey availability were suggested to play a role in divergence into multiple ecomorphs, but the effects of biotic factors require further study. Although we identified several environmental factors correlated with the presence of multiple ecomorphs, there were substantial data gaps across species and regions. An improved understanding of these systems may provide insight into both generalizable environmental factors involved in speciation in other systems, and potential ecological and evolutionary responses of species complexes when these variables are altered by environmental change.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

The interaction of resource use and gene flow on the phenotypic divergence of benthic and pelagic morphs of Icelandic Arctic charr (Salvelinus alpinus)

Conceptual models of adaptive divergence and ecological speciation in sympatry predict differential resource use, phenotype-environment correlations, and reduced gene flow among diverging phenotypes. While these predictions have been assessed in past studies, connections among them have rarely been assessed collectively. We examined relationships among phenotypic, ecological, and genetic variation in Arctic charr (Salvelinus alpinus) from six Icelandic localities that have undergone varying degrees of divergence into sympatric benthic and pelagic morphs. We characterized morphological variation with geometric morphometrics, tested for differential resource use between morphs using stable isotopes, and inferred the amount of gene flow from single nucleotide polymorphisms. Analysis of stable isotopic signatures indicated that sympatric morphs showed similar difference in resource use across populations, likely arising from the common utilization of niche space within each population. Carbon isotopic signature was also a significant predictor of individual variation in body shape and size, suggesting that variation in benthic and pelagic resource use is associated with phenotypic variation. The estimated percentage of hybrids between sympatric morphs varied across populations (from 0% to 15.6%) but the majority of fish had genotypes (ancestry coefficients) characteristic of pure morphs. Despite evidence of reduced gene flow between sympatric morphs, we did not detect the expected negative relationship between divergence in resource use and gene flow. Three lakes showed the expected pattern, but morphs in the fourth showed no detectable hybridization and had relatively low differences in resource use between them. This coupled with the finding that resource use and genetic differentiation had differential effects on body shape variation across populations suggests that reproductive isolation maintains phenotypic divergence between benthic and pelagic morphs when the effects of resource use are relatively low. Our ability to assess relationships between phenotype, ecology, and genetics deepens our understanding of the processes underlying adaptive divergence in sympatry.

Isotopic niche size of Coregonus artedi (sensu lato) Increases in the presence of Mysis diluviana, expanded habitat use and phenotypic diversity

Post-glacial colonization of lakes in Algonquin Park, Ontario, Canada resulted in food webs with cisco (Coregonus artedi sensu lato) and either Mysis diluviana or Chaoborus spp. as the dominant diel migrator. Mysis as prey, its diel movements and benthic occupancy, are hypothesized to be key elements of ecological opportunity for cisco diversity in the Laurentian Great Lakes. If correct, the hypothesis strongly implies that lakes with Mysis would have greater trophic niche size and drive greater adaptive radiation of cisco forms relative to lakes without Mysis. The dichotomy in diel migrator in Algonquin Park lakes was an opportunity to assess the isotopic niche size of cisco (δ15N and δ13C) and determine if niche size expands with Mysis presence. We found the presence of Mysis is necessary to expand isotopic niche size in our study lakes. The use of habitats not typically associated with the ancestral form of cisco (e.g., benthic habitats) and phenotypic diversity (blackfin and cisco) also continue to expand niche size in Mysis-based food webs. Partial ecological speciation based on a large niche space appears to be present in one lake (Cauchon Lake) where use of alternative habitats is the only real difference in cisco. The presence of blackfin expands niche space in Cedar and Radiant Lakes. This was not matched in Hogan Lake where niche space was relatively smaller with similar forms. Possible reasons for this discrepancy may be related to the asymmetric basin of Hogan Lake and whether the two forms overlap during cool and cold-water periods of the annual temperature cycle. By comparing trophic niche size among lakes with and without Mysis, we conclude that Mysis provides a key ecological opportunity for cisco diversity in our study lakes and likely more widely.

Ecological speciation in European whitefish is driven by a large-gaped predator

Lake-dwelling fish that form species pairs/flocks characterized by body size divergence are important model systems for speciation research. Although several sources of divergent selection have been identified in these systems, their importance for driving the speciation process remains elusive. A major problem is that in retrospect, we cannot distinguish selection pressures that initiated divergence from those acting later in the process. To address this issue, we studied the initial stages of speciation in European whitefish (Coregonus lavaretus) using data from 358 populations of varying age (26-10,000 years). We find that whitefish speciation is driven by a large-growing predator, the northern pike (Esox lucius). Pike initiates divergence by causing a largely plastic differentiation into benthic giants and pelagic dwarfs: ecotypes that will subsequently develop partial reproductive isolation and heritable differences in gill raker number. Using an eco-evolutionary model, we demonstrate how pike's habitat specificity and large gape size are critical for imposing a between-habitat trade-off, causing prey to mature in a safer place or at a safer size. Thereby, we propose a novel mechanism for how predators may cause dwarf/giant speciation in lake-dwelling fish species.

Evidence for host effect on the intestinal microbiota of whitefish (Coregonus sp.) species pairs and their hybrids

Investigating relationships between microbiota and their host is essential toward a full understanding of how animal adapt to their environment. Lake Whitefish offers a powerful system to investigate processes of adaptive divergence where the dwarf, limnetic species evolved repeatedly from the normal, benthic species. We compared the transient intestinal microbiota between both species from the wild and in controlled conditions, including their reciprocal hybrids. We sequenced the 16s rRNA gene V3-V4 regions to (a) test for parallelism in the transient intestinal microbiota among sympatric pairs, (b) test for transient intestinal microbiota differences among dwarf, normal, and hybrids reared under identical conditions, and (c) compare intestinal microbiota between wild and captive whitefish. A significant host effect on microbiota taxonomic composition was observed when all lakes were analyzed together and in three of the five species pairs. In captive whitefish, host effect was also significant. Microbiota of both reciprocal hybrids fell outside of that observed in the parental forms. Six genera formed a bacterial core which was present in captive and wild whitefish, suggesting a horizontal microbiota transmission. Altogether, our results complex interactions among the host, the microbiota, and the environment, and we propose that these interactions define three distinct evolutionary paths of the intestinal microbiota.

Evolution and diversity of two cisco forms in an outlet of glacial Lake Algonquin

The diversity of Laurentian Great Lakes ciscoes (Coregonus artedi, sensu lato) arose via repeated local adaptive divergence including deepwater ciscoes that are now extirpated or threatened. The nigripinnis form, or Blackfin Cisco, is extirpated from the Great Lakes and remains only in Lake Nipigon. Putative nigripinnis populations were recently discovered in sympatry with artedi in a historical drainage system of glacial Lake Algonquin, the precursor of lakes Michigan and Huron. Given the apparent convergence on Great Lakes form, we labeled this form blackfin. Here, we test the hypothesis that nigripinnis may have colonized this area from the Great Lakes as a distinct lineage. It would then represent a relict occurrence of the historical diversity of Great Lakes ciscoes. Alternatively, blackfin could have evolved in situ in several lakes. We captured more than 600 individuals in the benthic or pelagic habitat in 14 lakes in or near Algonquin Provincial Park (Ontario, Canada). Fish were compared based on habitat, morphology, and genetic variation at 6,676 SNPs. Contrary to our expectations, both cisco and blackfin belonged to an Atlantic lineage that colonized the area from the east, not from the Great Lakes. Sympatric cisco and blackfin were closely related while fish from different lakes were genetically differentiated, strongly suggesting the repeated in situ origin of each form. Across lakes, there was a continuum of ecological, morphological, and genetic differentiation that could be associated with alternative resources and lake characteristics. This study uncovers a new component of cisco diversity in inland lakes of Canada that evolved independently from ciscoes of the Laurentian Great lakes. The diversity of cisco revealed in this study and across their Canadian range presents a challenge for designating conservation units at the intraspecific level within the framework of the Committee on the Status of Endangered Wildlife in Canada (COSEWIC).

(Non)Parallel Evolution

Parallel evolution across replicate populations has provided evolutionary biologists with iconic examples of adaptation. When multiple populations colonize seemingly similar habitats, they may evolve similar genes, traits, or functions. Yet, replicated evolution in nature or in the laboratory often yields inconsistent outcomes: Some replicate populations evolve along highly similar trajectories, whereas other replicate populations evolve to different extents or in distinct directions. To understand these heterogeneous outcomes, biologists are increasingly treating parallel evolution not as a binary phenomenon but rather as a quantitative continuum ranging from parallel to nonparallel. By measuring replicate populations’ positions along this (non)parallel continuum, we can test hypotheses about evolutionary and ecological factors that influence the extent of repeatable evolution. We review evidence regarding the manifestation of (non)parallel evolution in the laboratory, in natural populations, and in applied contexts such as cancer. We enumerate the many genetic, ecological, and evolutionary processes that contribute to variation in the extent of parallel evolution.

Modeling the Multiple Facets of Speciation-with-Gene-Flow toward Inferring the Divergence History of Lake Whitefish Species Pairs (Coregonus clupeaformis)

Parallel divergence across replicated species pairs occurring in similar environmental contrasts may arise through distinct evolutionary scenarios. Deciphering whether such parallelism actually reflects repeated parallel divergence driven by divergent selection or a single divergence event with subsequent gene flow needs to be ascertained. Reconstructing historical gene flow is therefore of fundamental interest to understand how demography and selection jointly shaped genomic divergence during speciation. Here, we use an extended modeling framework to explore the multiple facets of speciation-with-gene-flow with demo-genetic divergence models that capture both temporal and genomic variation in effective population size and migration rate. We investigate the divergence history of replicate sympatric species pairs of Lake Whitefish (normal benthic and dwarf limnetic) characterized by variable degrees of ecological divergence and reproductive isolation. Genome-wide SNPs were used to document the extent of genetic differentiation in each species pair, and 26 divergence models were fitted and compared with the unfolded joint allele frequency spectrum of each pair. We found evidence that a recent (circa 3,000–4,000 generations) asymmetrical secondary contact between expanding postglacial populations has accompanied Whitefish diversification. Our results suggest that heterogeneous genomic differentiation has emerged through the combined effects of linked selection generating variable rates of lineage sorting across the genome during geographical isolation, and heterogeneous introgression eroding divergence at different rates across the genome upon secondary contact. This study thus provides a new retrospective insight into the historical demographic and selective processes that shaped a continuum of divergence associated with ecological speciation.

Molecular evolution and expression of oxygen transport genes in livebearing fishes (Poeciliidae) from hydrogen sulfide rich springs

Hydrogen sulfide (H2S) is a natural toxicant in some aquatic environments that has diverse molecular targets. It binds to oxygen transport proteins, rendering them non-functional by reducing oxygen-binding affinity. Hence, organisms permanently inhabiting H2S-rich environments are predicted to exhibit adaptive modifications to compensate for the reduced capacity to transport oxygen. We investigated 10 lineages of fish of the family Poeciliidae that have colonized freshwater springs rich in H2S—along with related lineages from non-sulfidic environments—to test hypotheses about the expression and evolution of oxygen transport genes in a phylogenetic context. We predicted shifts in the expression of and signatures of positive selection on oxygen transport genes upon colonization of H2S-rich habitats. Our analyses indicated significant shifts in gene expression for multiple hemoglobin genes in lineages that have colonized H2S-rich environments, and three hemoglobin genes exhibited relaxed selection in sulfidic compared to non-sulfidic lineages. However, neither changes in gene expression nor signatures of selection were consistent among all lineages in H2S-rich environments. Oxygen transport genes may consequently be predictable targets of selection during adaptation to sulfidic environments, but changes in gene expression and molecular evolution of oxygen transport genes in H2S-rich environments are not necessarily repeatable across replicated lineages.

Holobionts and ecological speciation: the intestinal microbiota of lake whitefish species pairs

BACKGROUND

It is well established that symbionts have considerable impact on their host, yet the investigation of the possible role of the holobiont in the host's speciation process is still in its infancy. In this study, we compared the intestinal microbiota among five sympatric pairs of dwarf (limnetic) and normal (benthic) lake whitefish Coregonus clupeaformis representing a continuum in the early stage of ecological speciation. We sequenced the 16s rRNA gene V3-V4 regions of the intestinal microbiota present in a total of 108 wild sympatric dwarf and normal whitefish as well as the water bacterial community from five lakes to (i) test for differences between the whitefish intestinal microbiota and the water bacterial community and (ii) test for parallelism in the intestinal microbiota of dwarf and normal whitefish.

RESULTS

The water bacterial community was distinct from the intestinal microbiota, indicating that intestinal microbiota did not reflect the environment, but rather the intrinsic properties of the host microbiota. Our results revealed a strong influence of the host (dwarf or normal) on the intestinal microbiota with pronounced conservation of the core intestinal microbiota (mean ~ 44% of shared genera). However, no clear evidence for parallelism was observed, whereby non-parallel differences between dwarf and normal whitefish were observed in three of the lakes while similar taxonomic composition was observed for the two other species pairs.

CONCLUSIONS

This absence of parallelism across dwarf vs. normal whitefish microbiota highlighted the complexity of the holobiont and suggests that the direction of selection could be different between the host and its microbiota.

Convergent evolution of reduced energy demands in extremophile fish

Convergent evolution in organismal function can arise from nonconvergent changes in traits that contribute to that function. Theory predicts that low resource availability and high maintenance costs in extreme environments select for reductions in organismal energy demands, which could be attained through modifications of body size or metabolic rate. We tested for convergence in energy demands and underlying traits by investigating livebearing fish (genus Poecilia) that have repeatedly colonized toxic, hydrogen sulphide-rich springs. We quantified variation in body size and routine metabolism across replicated sulphidic and non-sulphidic populations in nature, modelled total organismal energy demands, and conducted a common-garden experiment to test whether population differences had a genetic basis. Sulphidic populations generally exhibited smaller body sizes and lower routine metabolic rates compared to non-sulphidic populations, which together caused significant reductions in total organismal energy demands in extremophile populations. Although both mechanisms contributed to variation in organismal energy demands, variance partitioning indicated reductions of body size overall had a greater effect than reductions of routine metabolism. Finally, population differences in routine metabolism documented in natural populations were maintained in common-garden reared individuals, indicating evolved differences. In combination with other studies, these results suggest that reductions in energy demands may represent a common theme in adaptation to physiochemical stressors. Selection for reduced energy demand may particularly affect body size, which has implications for life history evolution in extreme environments.

Standing chromosomal variation in Lake Whitefish species pairs: the role of historical contingency and relevance for speciation

The role of chromosome changes in speciation remains a debated topic, although demographic conditions associated with divergence should promote their appearance. We tested a potential relationship between chromosome changes and speciation by studying two Lake Whitefish (Coregonus clupeaformis) lineages that recently colonized postglacial lakes following allopatry. A dwarf limnetic species evolved repeatedly from the normal benthic species, becoming reproductively isolated. Lake Whitefish hybrids experience mitotic and meiotic instability, which may result from structurally divergent chromosomes. Motivated by this observation, we test the hypothesis that chromosome organization differs between Lake Whitefish species pairs using cytogenetics. While chromosome and fundamental numbers are conserved between the species (2n = 80, NF = 98), we observe extensive polymorphism of subtle karyotype traits. We describe intrachromosomal differences associated with heterochromatin and repetitive DNA, and test for parallelism among three sympatric species pairs. Multivariate analyses support the hypothesis that differentiation at the level of subchromosomal markers mostly appeared during allopatry. Yet we find no evidence for parallelism between species pairs among lakes, consistent with colonization effect or postcolonization differentiation. The reported intrachromosomal polymorphisms do not appear to play a central role in driving adaptive divergence between normal and dwarf Lake Whitefish. We discuss how chromosomal differentiation in the Lake Whitefish system may contribute to the destabilization of mitotic and meiotic chromosome segregation in hybrids, as documented previously. The chromosome structures detected here are still difficult to sequence and assemble, demonstrating the value of cytogenetics as a complementary approach to understand the genomic bases of speciation.

Adaptation and acclimation of traits associated with swimming capacity in Lake Whitefish (coregonus clupeaformis) ecotypes

Background

Improved performance in a given ecological niche can occur through local adaptation, phenotypic plasticity, or a combination of these mechanisms. Evaluating the relative importance of these two mechanisms is needed to better understand the cause of intra specific polymorphism. In this study, we reared populations of Lake Whitefish (Coregonus clupeaformis) representing the’normal’ (benthic form) and the ‘dwarf’ (derived limnetic form) ecotypes in two different conditions (control and swim-training) to test the relative importance of adaptation and acclimation in the differentiation of traits related to swimming capacity. The dwarf whitefish is a more active swimmer than the normal ecotype, and also has a higher capacity for aerobic energy production in the swimming musculature. We hypothesized that dwarf fish would show changes in morphological and physiological traits consistent with reductions in the energetic costs of swimming and maintenance metabolism.

Results

We found differences in traits predicted to decrease the costs of prolonged swimming and standard metabolic rate and allow for a more active lifestyle in dwarf whitefish. Dwarf whitefish evolved a more streamlined body shape, predicted to lead to a decreased drag, and a smaller brain, which may decrease their standard metabolic rate. Contrary to predictions, we also found evidence of acclimation in liver size and metabolic enzyme activities.

Conclusion

Results support the view that local adaptation has contributed to the genetically-based divergence of traits associated with swimming activity. Presence of post-zygotic barriers limiting gene flow between these ecotype pairs may have favoured repeated local adaptation to the limnetic niches.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0732-y) contains supplementary material, which is available to authorized users.

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.

Identifying designatable units for intraspecific conservation prioritization: a hierarchical approach applied to the lake whitefish species complex (Coregonus spp.)

The concept of the designatable unit (DU) affords a practical approach to identifying diversity below the species level for conservation prioritization. However, its suitability for defining conservation units in ecologically diverse, geographically widespread and taxonomically challenging species complexes has not been broadly evaluated. The lake whitefish species complex (Coregonus spp.) is geographically widespread in the Northern Hemisphere, and it contains a great deal of variability in ecology and evolutionary legacy within and among populations, as well as a great deal of taxonomic ambiguity. Here, we employ a set of hierarchical criteria to identify DUs within the Canadian distribution of the lake whitefish species complex. We identified 36 DUs based on (i) reproductive isolation, (ii) phylogeographic groupings, (iii) local adaptation and (iv) biogeographic regions. The identification of DUs is required for clear discussion regarding the conservation prioritization of lake whitefish populations. We suggest conservation priorities among lake whitefish DUs based on biological consequences of extinction, risk of extinction and distinctiveness. Our results exemplify the need for extensive genetic and biogeographic analyses for any species with broad geographic distributions and the need for detailed evaluation of evolutionary history and adaptive ecological divergence when defining intraspecific conservation units.

RAD-QTL Mapping Reveals Both Genome-Level Parallelism and Different Genetic Architecture Underlying the Evolution of Body Shape in Lake Whitefish (Coregonus clupeaformis) Species Pairs

Parallel changes in body shape may evolve in response to similar environmental conditions, but whether such parallel phenotypic changes share a common genetic basis is still debated. The goal of this study was to assess whether parallel phenotypic changes could be explained by genetic parallelism, multiple genetic routes, or both. We first provide evidence for parallelism in fish shape by using geometric morphometrics among 300 fish representing five species pairs of Lake Whitefish. Using a genetic map comprising 3438 restriction site−associated DNA sequencing single-nucleotide polymorphisms, we then identified quantitative trait loci underlying body shape traits in a backcross family reared in the laboratory. A total of 138 body shape quantitative trait loci were identified in this cross, thus revealing a highly polygenic architecture of body shape in Lake Whitefish. Third, we tested for evidence of genetic parallelism among independent wild populations using both a single-locus method (outlier analysis) and a polygenic approach (analysis of covariation among markers). The single-locus approach provided limited evidence for genetic parallelism. However, the polygenic analysis revealed genetic parallelism for three of the five lakes, which differed from the two other lakes. These results provide evidence for both genetic parallelism and multiple genetic routes underlying parallel phenotypic evolution in fish shape among populations occupying similar ecological niches.

Parallel evolution in Ugandan crater lakes: repeated evolution of limnetic body shapes in haplochromine cichlid fish

Background

The enormous diversity found in East African cichlid fishes in terms of morphology, coloration, and behavior have made them a model for the study of speciation and adaptive evolution. In particular, haplochromine cichlids, by far the most species-rich lineage of cichlids, are a well-known textbook example for parallel evolution. Southwestern Uganda is an area of high tectonic activity, and is home to numerous crater lakes. Many Ugandan crater lakes were colonized, apparently independently, by a single lineage of haplochromine cichlids. Thereby, this system could be considered a natural experiment in which one can study the interaction between geographical isolation and natural selection promoting phenotypic diversification.

Results

We sampled 13 crater lakes and six potentially-ancestral older lakes and, using both mitochondrial and microsatellite markers, discovered strong genetic and morphological differentiation whereby (a) geographically close lakes tend to be genetically more similar and (b) three different geographic areas seem to have been colonized by three independent waves of colonization from the same source population. Using a geometric morphometric approach, we found that body shape elongation (i.e. a limnetic morphology) evolved repeatedly from the ancestral deeper-bodied benthic morphology in the clear and deep crater lake habitats.

Conclusions

A pattern of strong genetic and morphological differentiation was observed in the Ugandan crater lakes. Our data suggest that body shape changes have repeatedly evolved into a more limnetic-like form in several Ugandan crater lakes after independent waves of colonization from the same source population. The observed morphological changes in crater lake cichlids are likely to result from a common selective regime.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0287-3) contains supplementary material, which is available to authorized users.

Combining next-generation sequencing and online databases for microsatellite development in non-model organisms

Next-generation sequencing (NGS) is revolutionising marker development and the rapidly increasing amount of transcriptomes published across a wide variety of taxa is providing valuable sequence databases for the identification of genetic markers without the need to generate new sequences. Microsatellites are still the most important source of polymorphic markers in ecology and evolution. Motivated by our long-term interest in the adaptive radiation of a non-model species complex of whitefishes (Coregonus spp.), in this study, we focus on microsatellite characterisation and multiplex optimisation using transcriptome sequences generated by Illumina® and Roche-454, as well as online databases of Expressed Sequence Tags (EST) for the study of whitefish evolution and demographic history. We identified and optimised 40 polymorphic loci in multiplex PCR reactions and validated the robustness of our analyses by testing several population genetics and phylogeographic predictions using 494 fish from five lakes and 2 distinct ecotypes.

Speciation in ninespine stickleback: reproductive isolation and phenotypic divergence among cryptic species of Japanese ninespine stickleback

Although similar patterns of phenotypic diversification are often observed in phylogenetically independent lineages, differences in the magnitude and direction of phenotypic divergence have been also observed among independent lineages, even when exposed to the same ecological gradients. The stickleback family is a good model with which to explore the ecological and genetic basis of parallel and nonparallel patterns of phenotypic evolution, because there are a variety of populations and species that are locally adapted to divergent environments. Although the patterns of phenotypic divergence as well as the genetic and ecological mechanisms have been well characterized in threespine sticklebacks, Gasterosteus aculeatus, we know little about the patterns of phenotypic diversification in other stickleback lineages. In eastern Hokkaido, Japan, there are three species of ninespine sticklebacks, Pungitius tymensis and the freshwater type and the brackish-water type of the P. pungitius-P. sinensis species complex. They utilize divergent habitats along coast-stream gradients of rivers. Here, we investigated genetic, ecological and phenotypic divergence among three species of Japanese ninespine sticklebacks. Divergence in trophic morphology and salinity tolerance occurred in the direction predicted by the patterns observed in threespine sticklebacks. However, the patterns of divergence in armour plate were different from those previously found in threespine sticklebacks. Furthermore, the genetic basis of plate variation may differ from that in threespine sticklebacks. Because threespine sticklebacks are well-established model for evolutionary research, the sympatric trio of ninespine sticklebacks will be an invaluable resource for ecological and genetic studies on both common and lineage-specific patterns of phenotypic diversification.

Indirect trophic interactions with an invasive species affect phenotypic divergence in a top consumer

While phenotypic responses to direct species interactions are well studied, we know little about the consequences of indirect interactions for phenotypic divergence. In this study we used lakes with and without the zebra mussel to investigate effects of indirect trophic interactions on phenotypic divergence between littoral and pelagic perch. We found a greater phenotypic divergence between littoral and pelagic individuals in lakes with zebra mussels and propose a mussel-mediated increase in pelagic and benthic resource availability as a major factor underlying this divergence. Lakes with zebra mussels contained higher densities of large plankton taxa and large invertebrates. We suggest that this augmented resource availability improved perch foraging opportunities in both the littoral and pelagic zones. Perch in both habitats could hence express a more specialized foraging morphology, leading to an increased divergence of perch forms in lakes with zebra mussels. As perch do not prey on mussels directly, we conclude that the increased divergence results from indirect interactions with the mussels. Our results hence suggest that species at lower food web levels can indirectly affect phenotypic divergence in species at the top of the food chain.

Variable extent of parallelism in respiratory, circulatory, and neurological traits across lake whitefish species pairs

Parallel adaptive radiation events provide a powerful framework for investigations of ecology's contribution to phenotypic diversification. Ecologically driven divergence has been invoked to explain the repeated evolution of sympatric dwarf and normal lake whitefish (Coregonus clupeaformis) species in multiple lakes in eastern North America. Nevertheless, links between most putatively adaptive traits and ecological variation remain poorly defined within and among whitefish species pairs. Here, we examine four species pairs for variation in gill, heart, and brain size; three traits predicted to show strong phenotypic responses to ecological divergence. In each of the species pairs, normals exhibited larger body size standardized gills compared to dwarfs – a pattern that is suggestive of a common ecological driver of gill size divergence. Within lakes, the seasonal hypoxia experienced in the benthic environment is a likely factor leading to the requirement for larger gills in normals. Interestingly, the morphological pathways used to achieve larger gills varied between species pairs from Québec and Maine, which may imply subtle non-parallelism in gill size divergence related to differences in genetic background. There was also a non-significant trend toward larger hearts in dwarfs, the more active species of the two, whereas brain size varied exclusively among the lake populations. Taken together, our results suggest that the diversification of whitefish has been driven by parallel and non-parallel ecological conditions across lakes. Furthermore, the phenotypic response to ecological variation may depend on genetic background of each population.

The evolutionary ecology of dwarfism in three-spined sticklebacks

1. Body size is a defining phenotypic trait, but the ecological causes of its evolution are poorly understood. Most studies have considered only a single putative causal agent and have failed to recognise that different environmental agents are often correlated. 2. Darwin suggested that although trait variation across populations is often associated with abiotic variation, evolution is more likely to be driven by biotic factors correlated with the abiotic variation. This hypothesis has received little explicit attention. 3. We use structural equation modelling to quantify the relative importance of abiotic (pH, metal concentrations) and biotic (competition, predation) factors in the evolution of body size in three-spined sticklebacks Gasterosteus aculeatus on the island of North Uist, Scotland. We combine phenotypic data from multiple isolated populations, detailed characterisation of their environment and a common garden experiment that establishes the genetic basis of size differences. 4. Three-spined sticklebacks on North Uist show almost unprecedented intraspecific evolution of body size that has taken place rapidly (<16,000 years). The smallest fish mature at only 7% of the mass of ancestral, anadromous fish. Dwarfism is associated with reduced abundance of a smaller competitor species, the nine-spined stickleback Pungitius pungitius, and with low pH indicative of poor resource conditions. Dwarfism also tends to occur where an important predator, the brown trout Salmo trutta, is also small. The abundance of P. pungitius and the size of S. trutta are themselves related to underlying abiotic environmental variation. 5. Despite the close association between abiotic and biotic factors across populations, our results support Darwin's hypothesis that biotic factors, associated with variation in the abiotic environment, are more important in explaining evolution than is abiotic variation per se. This study demonstrates the importance of considering the relationships between environmental variables before conclusions can be drawn about the causes of (body size) evolution on islands.

Oxidative phosphorylation gene transcription in whitefish species pairs reveals patterns of parallel and nonparallel physiological divergence

Across multiple lakes in North America, lake whitefish (Coregonus clupeaformis) have independently evolved 'dwarf' and 'normal' sympatric species pairs that exhibit pronounced phenotypic and genetic divergence. In particular, traits associated with metabolism have been shown to be highly differentiated between whitefish species. Here, we examine the transcription of genes associated with the five mitochondrial and nuclear genome-encoded oxidative phosphorylation (OXPHOS) complexes, the primary physiological mechanism responsible for the production of ATP, in whitefish species pairs from Cliff Lake and Webster Lake in Maine, USA. We observed OXPHOS gene transcription divergence between dwarf and normal whitefish in each of the two lakes, with the former exhibiting transcription upregulation for genes associated with each of the OXPHOS complexes. We also observed a significant influence of lake on transcription levels for some of the genes, indicating that inter-lake ecological or genetic differences are contributing to variation in OXPHOS gene transcription levels. Together, our results support the hypothesis that metabolic divergence is a critical adaptation involved in whitefish speciation and implicate OXPHOS gene upregulation as a factor involved in meeting the enhanced energetic demands of dwarf whitefish. Further studies are now needed to evaluate the contribution of genetically vs. plasticity driven variation in transcription associated with this critical physiological pathway.

Genome-wide patterns of divergence during speciation: the lake whitefish case study

The nature, size and distribution of the genomic regions underlying divergence and promoting reproductive isolation remain largely unknown. Here, we summarize ongoing efforts using young (12 000 yr BP) species pairs of lake whitefish (Coregonus clupeaformis) to expand our understanding of the initial genomic patterns of divergence observed during speciation. Our results confirmed the predictions that: (i) on average, phenotypic quantitative trait loci (pQTL) show higher F(ST) values and are more likely to be outliers (and therefore candidates for being targets of divergent selection) than non-pQTL markers; (ii) large islands of divergence rather than small independent regions under selection characterize the early stages of adaptive divergence of lake whitefish; and (iii) there is a general trend towards an increase in terms of numbers and size of genomic regions of divergence from the least (East L.) to the most differentiated species pair (Cliff L.). This is consistent with previous estimates of reproductive isolation between these species pairs being driven by the same selective forces responsible for environment specialization. Altogether, dwarf and normal whitefish species pairs represent a continuum of both morphological and genomic differentiation contributing to ecological speciation. Admittedly, much progress is still required to more finely map and circumscribe genomic islands of speciation. This will be achieved through the use of next generation sequencing data but also through a better quantification of phenotypic traits moulded by selection as organisms adapt to new environmental conditions.

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.