Polymorphism in lake trout in Great Bear Lake: intra-lake morphological diversification at two spatial scales

Widespread admixture blurs population structure and confounds Lake Trout (Salvelinus namaycush) conservation even in the genomic era

Intraspecific variation is important for species' long-term persistence in changing environments. Conservation strategies targeting intraspecific variation often rely on the identification of management or policy units below the species level based on biological differences among populations. To identify management units, this paper examines intraspecific divergence of Lake Trout (Salvelinus namaycush) in Great Slave Lake (GSL), Canada, using low-coverage whole-genome sequencing data. Specifically, we evaluate genetic differentiation and assess the relationship with morphological, mitochondrial, and putatively adaptive divergence. We show that at least three and up to five genetically distinct Lake Trout populations co-occur in GSL and exhibit differences in spatial distribution and body size, with signatures of selection. However, admixture was widespread (60% of the fish) and population structure was shallow (average FST = 0.022). These findings highlight that, even in the era of whole genome sequencing, identifying discrete management units to implement conservation efforts and policy can remain challenging in systems where gene flow among genetically and ecologically distinct populations is ubiquitous. To give more recognition to this complexity, management efforts should also focus on the area where adaptive variation is evident, and evolutionary acts are at play, to better protect species' resilience and adaptive potential in some natural systems.

Winter ecology of specialist and generalist morphs of European whitefish, Coregonus lavaretus, in subarctic northern Europe

European whitefish is a model species for adaptive radiation of fishes in temperate and subarctic lakes. In northern Europe the most commonly observed morphotypes are a generalist (LSR) morph and a pelagic specialist (DR) morph. The evolution of a pelagic specialist morph is something of an enigma, however, as this region is characterized by long, dark winters with pelagic primary production limited to a brief window in late summer. We conducted the first winter-based study of polymorphic whitefish populations to determine the winter ecology of both morphs, and we combined seasonal diet and stable isotope analysis with several proxies of condition in three polymorphic whitefish populations. The generalist LSR morph fed on benthic and pelagic prey in summer but was solely reliant on benthic prey in winter. This was associated with a noticeable but moderate reduction in condition, lipid content and stomach fullness in winter relative to summer. In contrast, the DR whitefish occupied a strict pelagic niche in both seasons. A significant reduction in pelagic prey during winter resulted in severe decrease in condition, lipid content and stomach fullness in DR whitefish in winter relative to summer, with the pelagic morph apparently approaching starvation in winter. We suggest that this divergent approach to seasonal foraging is associated with the divergent life-history traits of both morphs.

A chromosome-anchored genome assembly for Lake Trout (Salvelinus namaycush)

Here, we present an annotated, chromosome‐anchored, genome assembly for Lake Trout (Salvelinus namaycush) – a highly diverse salmonid species of notable conservation concern and an excellent model for research on adaptation and speciation. We leveraged Pacific Biosciences long‐read sequencing, paired‐end Illumina sequencing, proximity ligation (Hi‐C) sequencing, and a previously published linkage map to produce a highly contiguous assembly composed of 7378 contigs (contig N50 = 1.8 Mb) assigned to 4120 scaffolds (scaffold N50 = 44.975 Mb). Long read sequencing data were generated using DNA from a female double haploid individual. 84.7% of the genome was assigned to 42 chromosome‐sized scaffolds and 93.2% of Benchmarking Universal Single Copy Orthologues were recovered, putting this assembly on par with the best currently available salmonid genomes. Estimates of genome size based on k‐mer frequency analysis were highly similar to the total size of the finished genome, suggesting that the entirety of the genome was recovered. A mitochondrial genome assembly was also produced. Self‐versus‐self synteny analysis allowed us to identify homeologs resulting from the salmonid specific autotetraploid event (Ss4R) as well as regions exhibiting delayed rediploidization. Alignment with three other salmonid genomes and the Northern Pike (Esox lucius) genome also allowed us to identify homologous chromosomes in related taxa. We also generated multiple resources useful for future genomic research on Lake Trout, including a repeat library and a sex‐averaged recombination map. A novel RNA sequencing data set for liver tissue was also generated in order to produce a publicly available set of annotations for 49,668 genes and pseudogenes. Potential applications of these resources to population genetics and the conservation of native populations are discussed.

Morphological and genetic concordance of cutthroat trout (Oncorhynchus clarkii) diversification from western North America

The cutthroat trout (Oncorhynchus clarkii (Richardson, 1836)) is one of the most widely distributed species of freshwater fish in western North America. Occupying a diverse range of habitats, they exhibit significant phenotypic variability that is often recognized by intraspecific taxonomy. Recent molecular phylogenies have described phylogenetic diversification across cutthroat trout populations, but no study has provided a range-wide morphological comparison of taxonomic divisions. In this study, we used linear- and geometric-based morphometrics to determine if phylogenetic and subspecies divisions correspond to morphological variation in cutthroat trout, using replicate populations from throughout the geographic range of the species. Our data indicate significant morphological divergence of intraspecific categories in some, but not all, cutthroat trout subspecies. We also compare morphological distance measures with distance measures of mtDNA sequence divergence. DNA sequence divergence was positively correlated with morphological distance measures, indicating that morphologically more similar subspecies have lower sequence divergence in comparison to morphologically distant subspecies. Given these results, integrating both approaches to describing intraspecific variation may be necessary for developing a comprehensive conservation plan in wide-ranging species.

Among-individual diet variation within a lake trout ecotype: Lack of stability of niche use

In a polyphenic species, differences in resource use are expected among ecotypes, and homogeneity in resource use is expected within an ecotype. Yet, using a broad resource spectrum has been identified as a strategy for fishes living in unproductive northern environments, where food is patchily distributed and ephemeral. We investigated whether specialization of trophic resources by individuals occurred within the generalist piscivore ecotype of lake trout from Great Bear Lake, Canada, reflective of a form of diversity. Four distinct dietary patterns of resource use within this lake trout ecotype were detected from fatty acid composition, with some variation linked to spatial patterns within Great Bear Lake. Feeding habits of different groups within the ecotype were not associated with detectable morphological or genetic differentiation, suggesting that behavioral plasticity caused the trophic differences. A low level of genetic differentiation was detected between exceptionally large-sized individuals and other piscivore individuals. We demonstrated that individual trophic specialization can occur within an ecotype inhabiting a geologically young system (8,000-10,000 yr BP), a lake that sustains high levels of phenotypic diversity of lake trout overall. The characterization of niche use among individuals, as done in this study, is necessary to understand the role that individual variation can play at the beginning of differentiation processes.

"And if you gaze long into an abyss, the abyss gazes also into thee": four morphs of Arctic charr adapting to a depth gradient in Lake Tinnsjøen

The origin of species is a central topic in biology. Ecological speciation might be a driver in adaptive radiation, providing a framework for understanding mechanisms, level, and rate of diversification. The Arctic charr Salvelinus alpinus L. is a polymorphic species with huge morphological and life‐history diversity in Holarctic water systems. We studied adaptive radiation of Arctic charr in the 460‐m‐deep Lake Tinnsjøen to (a) document eco‐morphology and life‐history traits of morphs, (b) estimate reproductive isolation of morphs, and (c) illuminate Holarctic phylogeography and lineages colonizing Lake Tinnsjøen. We compared Lake Tinnsjøen with four Norwegian outgroup populations. Four field‐assigned morphs were identified in Lake Tinnsjøen: the planktivore morph in all habitats except deep profundal, the dwarf morph in shallow‐moderate profundal, the piscivore morph mainly in shallow‐moderate profundal, and a new undescribed abyssal morph in the deep profundal. Morphs displayed extensive life‐history variation in age and size. A moderate‐to‐high concordance was observed among morphs and four genetic clusters from microsatellites. mtDNA suggested two minor endemic clades in Lake Tinnsjøen originating from one widespread colonizing clade in the Holarctic. All morphs were genetically differentiated at microsatellites (F_ST: 0.12–0.20), associated with different mtDNA clade frequencies. Analyses of outgroup lakes implied colonization from a river below Lake Tinnsjøen. Our findings suggest postglacial adaptive radiation of one colonizing mtDNA lineage with niche specialization along a depth–temperature–productivity–pressure gradient. Concordance between reproductive isolation and habitats of morphs implies ecological speciation as a mechanism. Particularly novel is the extensive morph diversification with depth into the often unexplored deepwater profundal habitat, suggesting we may have systematically underestimated biodiversity in lakes. In a biological conservation framework, it is imperative to protect endemic below‐species‐level biodiversity, particularly so since within‐species variation comprises an extremely important component of the generally low total biodiversity observed in the northern freshwater systems.

Microhabitat contributes to microgeographic divergence in threespine stickleback

Since the New Synthesis, most migration-selection balance theory has predicted that there should be negligible differentiation over small spatial scales (relative to dispersal), because gene flow should erode any effect of divergent selection. Nevertheless, there are classic examples of microgeographic divergence, which theory suggests can arise under specific conditions: exceptionally strong selection, phenotypic plasticity in philopatric individuals, or nonrandom dispersal. Here, we present evidence of microgeographic morphological variation within lake and stream populations of threespine stickleback (Gasterosteus aculeatus). It seems reasonable to assume that a given lake or stream population of fish is well-mixed. However, we found this assumption to be untenable. We examined trap-to-trap variation in 34 morphological traits measured on stickleback from 16 lakes and 16 streams. Most traits varied appreciably among traps within populations. Both between-trap distance and microhabitat characteristics such as depth and substrate explained some of the within-population morphological variance. Microhabitat was also associated with genotype at particular loci but there was no genetic isolation by distance, implying that heritable habitat preferences may contribute to microgeographic variation. Our study adds to growing evidence that microgeographic divergence can occur across small spatial scales within individuals' daily dispersal neighborhood where gene flow is expected to be strong.

Parallel selection on ecologically relevant gene functions in the transcriptomes of highly diversifying salmonids

Background

Salmonid fishes are characterised by a very high level of variation in trophic, ecological, physiological, and life history adaptations. Some salmonid taxa show exceptional potential for fast, within-lake diversification into morphologically and ecologically distinct variants, often in parallel; these are the lake-resident charr and whitefish (several species in the genera Salvelinus and Coregonus). To identify selection on genes and gene categories associated with such predictable diversifications, we analysed 2702 orthogroups (4.82 Mbp total; average 4.77 genes/orthogroup; average 1783 bp/orthogroup). We did so in two charr and two whitefish species and compared to five other salmonid lineages, which do not evolve in such ecologically predictable ways, and one non-salmonid outgroup.

Results

All selection analyses are based on Coregonus and Salvelinus compared to non-diversifying taxa. We found more orthogroups were affected by relaxed selection than intensified selection. Of those, 122 were under significant relaxed selection, with trends of an overrepresentation of serine family amino acid metabolism and transcriptional regulation, and significant enrichment of behaviour-associated gene functions. Seventy-eight orthogroups were under significant intensified selection and were enriched for signalling process and transcriptional regulation gene ontology terms and actin filament and lipid metabolism gene sets. Ninety-two orthogroups were under diversifying/positive selection. These were enriched for signal transduction, transmembrane transport, and pyruvate metabolism gene ontology terms and often contained genes involved in transcriptional regulation and development. Several orthogroups showed signs of multiple types of selection. For example, orthogroups under relaxed and diversifying selection contained genes such as ap1m2, involved in immunity and development, and slc6a8, playing an important role in muscle and brain creatine uptake. Orthogroups under intensified and diversifying selection were also found, such as genes syn3, with a role in neural processes, and ctsk, involved in bone remodelling.

Conclusions

Our approach pinpointed relevant genomic targets by distinguishing among different kinds of selection. We found that relaxed, intensified, and diversifying selection affect orthogroups and gene functions of ecological relevance in salmonids. Because they were found consistently and robustly across charr and whitefish and not other salmonid lineages, we propose these genes have a potential role in the replicated ecological diversifications.

Evaluating a pattern of ecological character displacement: charr jaw morphology and diet diverge in sympatry versus allopatry across catchments in Hokkaido, Japan

Similar species that overlap in sympatry may diverge in characters related to resource use as a result of evolution or phenotypic plasticity. Dolly Varden charr (Salvelinus malma) and whitespotted charr (S. leucomaenis) overlap along streams in Hokkaido, Japan, and compete by interference for invertebrate drift-foraging positions. Previous research has shown that as drift declines during summer, Dolly Varden shift foraging modes to capture benthic prey, a behaviour facilitated by their subterminal jaw morphology. We compare body and jaw morphology of Dolly Varden in sympatry vs. allopatry in two locations to test for character displacement. Statistical analysis showed significant divergence in characters related to foraging, which was correlated with variation in individual charr diets. Dolly Varden in sympatry had shorter heads and lower jaws than in allopatry, and even within sites charr with these characteristics fed less on drifting terrestrial invertebrates but more on benthic aquatic invertebrates. Those in allopatry had longer heads and lower jaws, and fed more on terrestrial invertebrates. The close proximity of sites in one stream suggests that Dolly Varden may display phenotypic plasticity similar to other charr, allowing rapid responses in morphology to the presence of competitors. These morphological shifts probably help them maintain positive fitness when competing with whitespotted charr in Hokkaido streams.

From top to bottom: Do Lake Trout diversify along a depth gradient in Great Bear Lake, NT, Canada?

Depth is usually considered the main driver of Lake Trout intraspecific diversity across lakes in North America. Given that Great Bear Lake is one of the largest and deepest freshwater systems in North America, we predicted that Lake Trout intraspecific diversity to be organized along a depth axis within this system. Thus, we investigated whether a deep-water morph of Lake Trout co-existed with four shallow-water morphs previously described in Great Bear Lake. Morphology, neutral genetic variation, isotopic niches, and life-history traits of Lake Trout across depths (0-150 m) were compared among morphs. Due to the propensity of Lake Trout with high levels of morphological diversity to occupy multiple habitat niches, a novel multivariate grouping method using a suite of composite variables was applied in addition to two other commonly used grouping methods to classify individuals. Depth alone did not explain Lake Trout diversity in Great Bear Lake; a distinct fifth deep-water morph was not found. Rather, Lake Trout diversity followed an ecological continuum, with some evidence for adaptation to local conditions in deep-water habitat. Overall, trout caught from deep-water showed low levels of genetic and phenotypic differentiation from shallow-water trout, and displayed higher lipid content (C:N ratio) and occupied a higher trophic level that suggested an potential increase of piscivory (including cannibalism) than the previously described four morphs. Why phenotypic divergence between shallow- and deep-water Lake Trout was low is unknown, especially when the potential for phenotypic variation should be high in deep and large Great Bear Lake. Given that variation in complexity of freshwater environments has dramatic consequences for divergence, variation in the complexity in Great Bear Lake (i.e., shallow being more complex than deep), may explain the observed dichotomy in the expression of intraspecific phenotypic diversity between shallow- vs. deep-water habitats. The ambiguity surrounding mechanisms driving divergence of Lake Trout in Great Bear Lake should be seen as reflective of the highly variable nature of ecological opportunity and divergent natural selection itself.

Genetic and phenotypic variation along an ecological gradient in lake trout Salvelinus namaycush

Background

Adaptive radiation involving a colonizing phenotype that rapidly evolves into at least one other ecological variant, or ecotype, has been observed in a variety of freshwater fishes in post-glacial environments. However, few studies consider how phenotypic traits vary with regard to neutral genetic partitioning along ecological gradients. Here, we present the first detailed investigation of lake trout Salvelinus namaycush that considers variation as a cline rather than discriminatory among ecotypes. Genetic and phenotypic traits organized along common ecological gradients of water depth and geographic distance provide important insights into diversification processes in a lake with high levels of human disturbance from over-fishing.

Results

Four putative lake trout ecotypes could not be distinguished using population genetic methods, despite morphological differences. Neutral genetic partitioning in lake trout was stronger along a gradient of water depth, than by locality or ecotype. Contemporary genetic migration patterns were consistent with isolation-by-depth. Historical gene flow patterns indicated colonization from shallow to deep water. Comparison of phenotypic (Pst) and neutral genetic variation (Fst) revealed that morphological traits related to swimming performance (e.g., buoyancy, pelvic fin length) departed more strongly from neutral expectations along a depth gradient than craniofacial feeding traits. Elevated phenotypic variance with increasing water depth in pelvic fin length indicated possible ongoing character release and diversification. Finally, differences in early growth rate and asymptotic fish length across depth strata may be associated with limiting factors attributable to cold deep-water environments.

Conclusion

We provide evidence of reductions in gene flow and divergent natural selection associated with water depth in Lake Superior. Such information is relevant for documenting intraspecific biodiversity in the largest freshwater lake in the world for a species that recently lost considerable genetic diversity and is now in recovery. Unknown is whether observed patterns are a result of an early stage of incipient speciation, gene flow-selection equilibrium, or reverse speciation causing formerly divergent ecotypes to collapse into a single gene pool.

Electronic supplementary material

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

Multiple generalist morphs of Lake Trout: Avoiding constraints on the evolution of intraspecific divergence?

A generalist strategy, as an adaptation to environmental heterogeneity, is common in Arctic freshwater systems, often accompanied, however, by intraspecific divergence that promotes specialization in niche use. To better understand how resources may be partitioned in a northern system that supports intraspecific diversity of Lake Trout, trophic niches were compared among four shallow‐water morphotypes in Great Bear Lake (N65^° 56′ 39″, W120^° 50′ 59″). Bayesian mixing model analyses of stable isotopes of carbon and nitrogen were conducted on adult Lake Trout. Major niche overlap in resource use among four Lake Trout morphotypes was found within littoral and pelagic zones, which raises the question of how such polymorphism can be sustained among opportunistic generalist morphotypes. Covariances of our morphological datasets were tested against δ¹³C and δ¹⁵N values. Patterns among morphotypes were mainly observed for δ¹⁵N. This link between ecological and morphological differentiation suggested that selection pressure(s) operate at the trophic level (δ¹⁵N), independent of habitat, rather than along the habitat‐foraging opportunity axis (δ¹³C). The spatial and temporal variability of resources in Arctic lakes, such as Great Bear Lake, may have favored the presence of multiple generalists showing different degrees of omnivory along a weak benthic–pelagic gradient. Morphs 1–3 had more generalist feeding habits using both benthic and pelagic habitats than Morph 4, which was a top‐predator specialist in the pelagic habitat. Evidence for frequent cannibalism in Great Bear Lake was found across all four morphotypes and may also contribute to polymorphism. We suggest that the multiple generalist morphs described here from Great Bear Lake are a unique expression of diversity due to the presumed constraints on the evolution of generalists and contrast with the development of multiple specialists, the standard response to intraspecific divergence.

Striking Phenotypic Variation yet Low Genetic Differentiation in Sympatric Lake Trout (Salvelinus namaycush)

The study of population differentiation in the context of ecological speciation is commonly assessed using populations with obvious discreteness. Fewer studies have examined diversifying populations with occasional adaptive variation and minor reproductive isolation, so factors impeding or facilitating the progress of early stage differentiation are less understood. We detected non-random genetic structuring in lake trout (Salvelinus namaycush) inhabiting a large, pristine, postglacial lake (Mistassini Lake, Canada), with up to five discernible genetic clusters having distinctions in body shape, size, colouration and head shape. However, genetic differentiation was low (F_ST = 0.017) and genetic clustering was largely incongruent between several population- and individual-based clustering approaches. Genotype- and phenotype-environment associations with spatial habitat, depth and fish community structure (competitors and prey) were either inconsistent or weak. Striking morphological variation was often more continuous within than among defined genetic clusters. Low genetic differentiation was a consequence of relatively high contemporary gene flow despite large effective population sizes, not migration-drift disequilibrium. Our results suggest a highly plastic propensity for occupying multiple habitat niches in lake trout and a low cost of morphological plasticity, which may constrain the speed and extent of adaptive divergence. We discuss how factors relating to niche conservatism in this species may also influence how plasticity affects adaptive divergence, even where ample ecological opportunity apparently exists.

Assortative mating but no evidence of genetic divergence in a species characterized by a trophic polymorphism

Disruptive selection is a process that can result in multiple subgroups within a population, which is referred to as diversification. Foraging-related diversification has been described in many taxa, but many questions remain about the contribution of such diversification to reproductive isolation and potentially sympatric speciation. Here, we use stable isotope analysis of diet and morphological analysis of body shape to examine phenotypic divergence between littoral and pelagic foraging ecomorphs in a population of pumpkinseed sunfish (Lepomis gibbosus). We then examine reproductive isolation between ecomorphs by comparing the isotopic compositions of nesting males to eggs from their nests (a proxy for maternal diet) and use nine microsatellite loci to examine genetic divergence between ecomorphs. Our data support the presence of distinct foraging ecomorphs in this population and indicate that there is significant positive assortative mating based on diet. We did not find evidence of genetic divergence between ecomorphs, however, indicating that isolation is either relatively recent or is not strong enough to result in genetic divergence at the microsatellite loci. Based on our findings, pumpkinseed sunfish represent a system in which to further explore the mechanisms by which natural and sexual selection contribute to diversification, prior to the occurrence of sympatric speciation.