Relationships between lake ecology and morphological characters in Icelandic Arctic charr, Salvelinus alpinus

The role of neutral and adaptive evolutionary processes on patterns of genetic diversity across small cave-dwelling populations of Icelandic Arctic charr (Salvelinus alpinus)

Understanding the adaptability of small populations in the face of environmental change is a central problem in evolutionary biology. Solving this problem is challenging because neutral evolutionary processes that operate on historical and contemporary timescales can override the effects of selection in small populations. We assessed the effects of isolation by colonization (IBC), isolation by dispersal limitation (IBDL) as reflected by a pattern of isolation by distance (IBD), and isolation by adaptation (IBA) and the roles of genetic drift and gene flow on patterns of genetic differentiation among 19 cave-dwelling populations of Icelandic Arctic charr (Salvelinus alpinus). We detected evidence of IBC based on the genetic affinity of nearby cave populations and the genetic relationships between the cave populations and the presumed ancestral population in the lake. A pattern of IBD was evident regardless of whether high-level genetic structuring (IBC) was taken into account. Genetic signatures of bottlenecks and lower genetic diversity in smaller populations indicate the effect of drift. Estimates of gene flow and fish movement suggest that gene flow is limited to nearby populations. In contrast, we found little evidence of IBA as patterns of local ecological and phenotypic variation showed little association with genetic differentiation among populations. Thus, patterns of genetic variation in these small populations likely reflect localized gene flow and genetic drift superimposed onto a larger-scale structure that is largely a result of colonization history. Our simultaneous assessment of the effects of neutral and adaptive processes in a tractable and replicated system has yielded novel insights into the evolution of small populations on both historical and contemporary timescales and over a smaller spatial scale than is typically studied.

Evaluating the correlation between genome-wide diversity and the release of plastic phenotypic variation in experimental translocations to novel natural environments

Phenotypic reaction norms are often shaped and constrained by selection and are important for allowing organisms to respond to environmental change. However, selection cannot constrain reaction norms for environmental conditions that populations have not experienced. Consequently, cryptic neutral genetic variation for the reaction norm can accumulate such that a release of phenotypic variation occurs upon exposure to novel14 conditions. Most genomic diversity behaves as if functionally neutral. Therefore, genome-wide diversity metrics may correlate with levels of cryptic genetic variation and, as a result, exhibit a positive relationship with a release of phenotypic variation in novel environments. To test this hypothesis, we conducted translocations of juvenile brook trout (Salvelinus fontinalis) from 12 populations to novel uninhabited ponds that represented a gradient of environmental conditions. We assessed reaction norms for morphological traits (body size and four morphometric relative warps) across pond environmental gradients and evaluated the effect of genome-wide heterozygosity on phenotypic variability. All traits displayed plastic reaction norms. Overall, we found some evidence that a release of phenotypic variation consistent with cryptic genetic variation can occur in novel environmental conditions. However, the extent to which this release correlated with average genome-wide diversity was limited to only one of five traits examined. Our results suggest a limited link between genomic diversity26 and the accumulation of cryptic genetic variation in reaction norms. Similarly, reaction norms were constrained for many of the morphological traits examined. Past conditions may have constrained reaction norms in the putatively novel environments despite significant deviations from contemporary source population habitat. Additionally, as a generalist colonizing species brook trout may exhibit plastic phenotypes across a wide range of environmental conditions.

Temporal trends, lake-to-lake variation, and climate effects on Arctic char (Salvelinus alpinus) mercury concentrations from six High Arctic lakes in Nunavut, Canada

Climate warming and mercury (Hg) are concurrently influencing Arctic ecosystems, altering their functioning and threatening food security. Non-anadromous Arctic char (Salvelinus alpinus) in small lakes were used to biomonitor these two anthropogenic stressors, because this iconic Arctic species is a long-lived top predator in relatively simple food webs, and yet population characteristics vary greatly, reflecting differences between lake systems. Mercury concentrations in six landlocked Arctic char populations on Cornwallis Island, Nunavut have been monitored as early as 1989, providing a novel dataset to examine differences in muscle [Hg] among char populations, temporal trends, and the relationship between climate patterns and Arctic char [Hg]. We found significant lake-to-lake differences in length-adjusted Arctic char muscle [Hg], which varied by up to 9-fold. Arctic char muscle [Hg] was significantly correlated to dissolved and particulate organic carbon concentrations in water; neither watershed area or vegetation cover explained differences. Three lakes exhibited significant temporal declines in length-adjusted [Hg] in Arctic char; the other three lakes had no significant trends. Though precipitation, temperature, wind speed, and sea ice duration were tested, no single climate variable was significantly correlated to length-adjusted [Hg] across populations. However, Arctic char Hg in Resolute Lake exhibited a significant correlation with sea ice duration, which is likely closely linked to lake ice duration, and which may impact Hg processing in lakes. Additionally, Arctic char [Hg] in Amituk Lake was significantly correlated to snow fall, which may be linked to Hg deposition. The lack of consistent temporal trends in neighboring char populations indicates that currently, within lake processes are the strongest drivers of [Hg] in char in the study lakes and potentially in other Arctic lakes, and that the influence of climate change will likely vary from lake to lake.

Natural selection for body shape in resource polymorphic Icelandic Arctic charr

Resource polymorphisms exhibit remarkable intraspecific diversity and in many cases are expected to be maintained by diversifying selection. Phenotypic trade-offs can constrain morphologically intermediate individuals from effectively exploiting both alternate resources, resulting in ecological barriers to gene flow. Determining if and how phenotypic trade-offs cause fitness variation in the wild is challenging because of phenotypic and environmental correlations associated with alternative resource strategies. We investigated multiple pathways through which morphology could affect organismal performance, as measured by growth rate, and whether these effects generate diversifying selection in polymorphic Icelandic Arctic charr (Salvelinus alpinus) populations. We considered direct effects of morphology on growth and indirect effects via trophic resource use, estimated by stable isotopic signatures, and via parasitism associated with trophic resources. We sampled over 3 years in (lakes) Thingvallavatn and Vatnshlíðarvatn using the extended selection gradient path analytical approach and estimating size-dependent mortality. We found evidence for diversifying selection only in Thingvallavatn: more streamlined and terminally mouthed planktivore charr experienced greater growth, with the opposite pattern in small benthic charr. However, this effect was mediated by parasitism and nontrophic pathways, rather than trophic performance as often expected. Detection of between-morph differences in the presence (Vatnshlíðarvatn) and direction (Thingvallavatn) of size-dependent mortality, together with nontrophic effects of shape, suggests that a morphological trophic performance explanation for polymorphism is insufficient. This rare insight into selection during early diversification suggests that a complex of interacting local factors must be considered to understand how phenotype influences fitness, despite morphological variation reflecting intuitive trade-off explanations.

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.

Differential gene expression during early development in recently evolved and sympatric Arctic charr morphs

Phenotypic differences between closely related taxa or populations can arise through genetic variation or be environmentally induced, leading to altered transcription of genes during development. Comparative developmental studies of closely related species or variable populations within species can help to elucidate the molecular mechanisms related to evolutionary divergence and speciation. Studies of Arctic charr (Salvelinus alpinus) and related salmonids have revealed considerable phenotypic variation among populations and in Arctic charr many cases of extensive variation within lakes (resource polymorphism) have been recorded. One example is the four Arctic charr morphs in the ∼10,000 year old Lake Thingvallavatn, which differ in numerous morphological and life history traits. We set out to investigate the molecular and developmental roots of this polymorphism by studying gene expression in embryos of three of the morphs reared in a common garden set-up. We performed RNA-sequencing, de-novo transcriptome assembly and compared gene expression among morphs during an important timeframe in early development, i.e., preceding the formation of key trophic structures. Expectedly, developmental time was the predominant explanatory variable. As the data were affected by some form of RNA-degradation even though all samples passed quality control testing, an estimate of 3'-bias was the second most common explanatory variable. Importantly, morph, both as an independent variable and as interaction with developmental time, affected the expression of numerous transcripts. Transcripts with morph effect, separated the three morphs at the expression level, with the two benthic morphs being more similar. However, Gene Ontology analyses did not reveal clear functional enrichment of transcripts between groups. Verification via qPCR confirmed differential expression of several genes between the morphs, including regulatory genes such as AT-Rich Interaction Domain 4A (arid4a) and translin (tsn). The data are consistent with a scenario where genetic divergence has contributed to differential expression of multiple genes and systems during early development of these sympatric Arctic charr morphs.

Variation in the magnitude of morphological and dietary differences between individuals among populations of small benthic Arctic charr in relation to ecological factors

The early stages of intraspecific diversity are important for the evolution of diversification and speciation. Early stages of diversification can be seen in individual specialization, where individuals consume only a portion of the diet of the population as a whole, and how such specialization is related to phenotypic diversity within populations. Here, we study the strength of the relationship between morphological and dietary distances among individuals in eighteen populations of Icelandic small benthic charr. We furthermore studied if the strength of the relationship could be related to variation in local ecological factors these populations inhabit. In all the populations studied, there was a clear relationship between morphological and dietary distances, indicating that fish that had similar morphology were at the same time-consuming similar food items. Our findings show a systematic variation in the relationship between morphology and diet at early stages of diversification in a highly specialized small benthic charr morph. The results show the importance of fine scale comparisons within populations and furthermore the value that systematic comparisons among populations under parallel evolution can contribute toward our increased understanding of evolutionary and ecological processes.

Ecosystem size predicts eco-morphological variability in a postglacial diversification

Identifying the processes by which new phenotypes and species emerge has been a long‐standing effort in evolutionary biology. Young adaptive radiations provide a model to study patterns of morphological and ecological diversification in environmental context. Here, we use the recent radiation (ca. 12k years old) of the freshwater fish Arctic charr (Salvelinus alpinus) to identify abiotic and biotic environmental factors associated with adaptive morphological variation. Arctic charr are exceptionally diverse, and in postglacial lakes there is strong evidence of repeated parallel evolution of similar morphologies associated with foraging. We measured head depth (a trait reflecting general eco‐morphology and foraging ecology) of 1,091 individuals across 30 lake populations to test whether fish morphological variation was associated with lake bathymetry and/or ecological parameters. Across populations, we found a significant relationship between the variation in head depth of the charr and abiotic environmental characteristics: positively with ecosystem size (i.e., lake volume, surface area, depth) and negatively with the amount of littoral zone. In addition, extremely robust‐headed phenotypes tended to be associated with larger and deeper lakes. We identified no influence of co‐existing biotic community on Arctic charr trophic morphology. This study evidences the role of the extrinsic environment as a facilitator of rapid eco‐morphological diversification.

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.

The developmental transcriptome of contrasting Arctic charr (Salvelinus alpinus) morphs

Species and populations with parallel evolution of specific traits can help illuminate how predictable adaptations and divergence are at the molecular and developmental level. Following the last glacial period, dwarfism and specialized bottom feeding morphology evolved rapidly in several landlocked Arctic charrSalvelinus alpinuspopulations in Iceland.  To study the genetic divergence between small benthic morphs and limnetic morphs, we conducted RNA-sequencing charr embryos at four stages in early development. We studied two stocks with contrasting morphologies: the small benthic (SB) charr from Lake Thingvallavatn and Holar aquaculture (AC) charr.The data reveal significant differences in expression of several biological pathways during charr development. There was also an expression difference between SB- and AC-charr in genes involved in energy metabolism and blood coagulation genes. We confirmed differing expression of five genes in whole embryos with qPCR, includinglysozymeandnatterin-likewhich was previously identified as a fish-toxin of a lectin family that may be a putative immunopeptide. We also verified differential expression of 7 genes in the developing head that associated consistently with benthic v.s.limnetic morphology (studied in 4 morphs). Comparison of single nucleotide polymorphism (SNP) frequencies reveals extensive genetic differentiation between the SB and AC-charr (~1300 with more than 50% frequency difference). Curiously, three derived alleles in the otherwise conserved 12s and 16s mitochondrial ribosomal RNA genes are found in benthic charr.The data implicate multiple genes and molecular pathways in divergence of small benthic charr and/or the response of aquaculture charr to domestication. Functional, genetic and population genetic studies on more freshwater and anadromous populations are needed to confirm the specific loci and mutations relating to specific ecological traits in Arctic charr.

The developmental transcriptome of contrasting Arctic charr ( Salvelinus alpinus) morphs

Species and populations with parallel evolution of specific traits can help illuminate how predictable adaptations and divergence are at the molecular and developmental level. Following the last glacial period, dwarfism and specialized bottom feeding morphology evolved rapidly in several landlocked Arctic charr Salvelinus alpinus populations in Iceland.  

To study the genetic divergence between small benthic morphs and limnetic morphs, we conducted RNA-sequencing charr embryos at four stages in early development. We studied two stocks with contrasting morphologies: the small benthic (SB) charr from Lake Thingvallavatn and Holar aquaculture (AC) charr.

The data reveal significant differences in expression of several biological pathways during charr development. There was also an expression difference between SB- and AC-charr in genes involved in energy metabolism and blood coagulation genes. We confirmed differing expression of five genes in whole embryos with qPCR, including lysozyme and natterin-like which was previously identified as a fish-toxin of a lectin family that may be a putative immunopeptide. We also verified differential expression of 7 genes in the developing head that associated consistently with benthic v.s.limnetic morphology (studied in 4 morphs). Comparison of single nucleotide polymorphism (SNP) frequencies reveals extensive genetic differentiation between the SB and AC-charr (~1300 with more than 50% frequency difference). Curiously, three derived alleles in the otherwise conserved 12s and 16s mitochondrial ribosomal RNA genes are found in benthic charr.

The data implicate multiple genes and molecular pathways in divergence of small benthic charr and/or the response of aquaculture charr to domestication. Functional, genetic and population genetic studies on more freshwater and anadromous populations are needed to confirm the specific loci and mutations relating to specific ecological traits in Arctic charr.

Morphological variation among wild populations of Chinese rare minnow (Gobiocypris rarus): deciphering the role of evolutionary processes

Gobiocypris rarus Ye et Fu ( 1983 ) is an endemic cyprinid fish in China, and is considered to be an endangered species. From a conservation perspective, its population structure is interesting. In the present study, morphological variation, including morphometric and meristic traits, was assessed among wild samples collected in the upper Yangtze River basin. There were no significant meristic differences between sexes or among populations, except for scales in lateral line (LS). However, there were significant morphometric differences not only between sexes but also among populations. In discriminant function analysis, the first four discriminant functions explained 75.5% and 78% of the among-population morphometric variation for males and females, respectively. Almost all the truss network morphometric traits significantly contributed to population discrimination. By using all of the truss network morphometric traits, overall random assignments of male and female individuals into their original population were both 73.5%. In addition, the degree of differentiation in phenotypic traits (PST) significantly exceeds that in neutral molecular markers (FST). However, no significant correlation between PST and FST was found in males or females. Overall, these results suggest that two evolutionary processes, including phenotypic plasticity and natural selection, may contribute to the morphological patterns observed in G. rarus.

Differentiation at the MHCIIα and Cath2 loci in sympatric Salvelinus alpinus resource morphs in Lake Thingvallavatn

Northern freshwater fish may be suitable for the genetic dissection of ecological traits because they invaded new habitats after the last ice age (∼10.000 years ago). Arctic charr (Salvelinus alpinus) colonizing streams and lakes in Iceland gave rise to multiple populations of small benthic morphotypes, often in sympatry with a pelagic morphotype. Earlier studies have revealed significant, but subtle, genetic differentiation between the three most common morphs in Lake Thingvallavatn. We conducted a population genetic screen on four immunological candidate genes Cathelicidin 2 (Cath2), Hepcidin (Hamp), Liver expressed antimicrobial peptide 2a (Leap-2a), and Major Histocompatibility Complex IIα (MHCIIα) and a mitochondrial marker (D-loop) among the three most common Lake Thingvallavatn charr morphs. Significant differences in allele frequencies were found between morphs at the Cath2 and MHCIIα loci. No such signal was detected in the D-loop nor in the other two immunological genes. In Cath2 the small benthic morph deviated from the other two (FST  = 0.13), one of the substitutions detected constituting an amino acid replacement polymorphism in the antimicrobial peptide. A more striking difference was found in the MHCIIα. Two haplotypes were very common in the lake, and their frequency differed greatly between the morphotypes (from 22% to 93.5%, FST  = 0.67). We then expanded our study by surveying the variation in Cath2 and MHCIIα in 9 Arctic charr populations from around Iceland. The populations varied greatly in terms of allele frequencies at Cath2, but the variation did not correlate with morphotype. At the MHCIIα locus, the variation was nearly identical to the variation in the two benthic morphs of Lake Thingvallavatn. The results are consistent with a scenario where parts of the immune systems have diverged substantially among Arctic charr populations in Iceland, after colonizing the island ∼10.000 years ago.

Parallel and non-parallel morphological divergence among foraging specialists in European whitefish (Coregonus lavaretus)

Parallel phenotypic evolution occurs when independent populations evolve similar traits in response to similar selective regimes. However, populations inhabiting similar environments also frequently show some phenotypic differences that result from non-parallel evolution. In this study, we quantified the relative importance of parallel evolution to similar foraging regimes and non-parallel lake-specific effects on morphological variation in European whitefish (Coregonus lavaretus). We found evidence for both lake-specific morphological characteristics and parallel morphological divergence between whitefish specializing in feeding on profundal and littoral resources in three separate lakes. Foraging specialists expressed similar phenotypes in different lakes in both overall body shape and selected measured morphological traits. The morphology of the two whitefish specialists resembled that predicted from other fish species, supporting the conclusion of an adaptive significance of the observed morphological characteristics. Our results indicate that divergent natural selection resulting from foraging specialization is driving and/or maintaining the observed parallel morphological divergence. Whitefish in this study may represent an early stage of divergence towards the evolution of specialized morphs.

Fine-scale parallel patterns in diversity of small benthic Arctic charr (Salvelinus alpinus) in relation to the ecology of lava/groundwater habitats

It is critical to study factors that are important for origin and maintenance of biological diversity. A comparative approach involving a large number of populations is particularly useful. We use this approach to study the relationship between ecological factors and phenotypic diversity in Icelandic Arctic charr (Salvelinus alpinus). Numerous populations of small benthic charr have evolved in lava springs in Iceland. These charr appear morphologically similar, but differ in important morphological features related to feeding. We found a clear relationship between diversity in morphology, diet, and ecological factors among populations. In particular, there were clear differences in morphology and diet between fish coming from habitats where the lava spring flowed on as a stream compared to habitats where the lava spring flowed into a pond. Our study shows that ecological factors are important for the origin and maintenance of biological diversity. The relationship between phenotype and ecological factors are observed on a fine scale, when comparing numerous populations that are phenotypically similar. This strongly suggests that for understanding, managing, and conserving biological diversity important ecological variables have to be taken into the account.