Extensive genetic differentiation between recently evolved sympatric Arctic charr morphs

Diversity in the internal functional feeding elements of sympatric morphs of Arctic charr (Salvelinus alpinus)

The diversity of functional feeding anatomy is particularly impressive in fishes and correlates with various interspecific ecological specializations. Intraspecific polymorphism can manifest in divergent feeding morphology and ecology, often along a benthic-pelagic axis. Arctic charr (Salvelinus alpinus) is a freshwater salmonid known for morphological variation and sympatric polymorphism and in Lake Þingvallavatn, Iceland, four morphs of charr coexist that differ in preferred prey, behaviour, habitat use, and external feeding morphology. We studied variation in six upper and lower jaw bones in adults of these four morphs using geometric morphometrics and univariate statistics. We tested for allometric differences in bone size and shape among morphs, morph effects on bone size and shape, and divergence along the benthic-pelagic axis. We also examined the degree of integration between bone pairs. We found differences in bone size between pelagic and benthic morphs for two bones (dentary and premaxilla). There was clear bone shape divergence along a benthic-pelagic axis in four bones (dentary, articular-angular, premaxilla and maxilla), as well as allometric shape differences between morphs in the dentary. Notably for the dentary, morph explained more shape variation than bone size. Comparatively, benthic morphs possess a compact and taller dentary, with shorter dentary palate, consistent with visible (but less prominent) differences in external morphology. As these morphs emerged in the last 10,000 years, these results indicate rapid functional evolution of specific feeding structures in arctic charr. This sets the stage for studies of the genetics and development of rapid and parallel craniofacial evolution.

Genomic underpinnings of head and body shape in Arctic charr ecomorph pairs

Across its Holarctic range, Arctic charr (Salvelinus alpinus) populations have diverged into distinct trophic specialists across independent replicate lakes. The major aspect of divergence between ecomorphs is in head shape and body shape, which are ecomorphological traits reflecting niche use. However, whether the genomic underpinnings of these parallel divergences are consistent across replicates was unknown but key for resolving the substrate of parallel evolution. We investigated the genomic basis of head shape and body shape morphology across four benthivore-planktivore ecomorph pairs of Arctic charr in Scotland. Through genome-wide association analyses, we found genomic regions associated with head shape (89 SNPs) or body shape (180 SNPs) separately and 50 of these SNPs were strongly associated with both body and head shape morphology. For each trait separately, only a small number of SNPs were shared across all ecomorph pairs (3 SNPs for head shape and 10 SNPs for body shape). Signs of selection on the associated genomic regions varied across pairs, consistent with evolutionary demography differing considerably across lakes. Using a comprehensive database of salmonid QTLs newly augmented and mapped to a charr genome, we found several of the head- and body-shape-associated SNPs were within or near morphology QTLs from other salmonid species, reflecting a shared genetic basis for these phenotypes across species. Overall, our results demonstrate how parallel ecotype divergences can have both population-specific and deeply shared genomic underpinnings across replicates, influenced by differences in their environments and demographic histories.

Variation in personality shaped by evolutionary history, genotype and developmental plasticity in response to feeding modalities in the Arctic charr

Animal personality has been shown to be influenced by both genetic and environmental factors and shaped by natural selection. Currently, little is known about mechanisms influencing the development of personality traits. This study examines the extent to which personality development is genetically influenced and/or environmentally responsive (plastic). We also investigated the role of evolutionary history, assessing whether personality traits could be canalized along a genetic and ecological divergence gradient. We tested the plastic potential of boldness in juveniles of five Icelandic Arctic charr morphs (Salvelinus alpinus), including two pairs of sympatric morphs, displaying various degrees of genetic and ecological divergence from the ancestral anadromous charr, split between treatments mimicking benthic versus pelagic feeding modalities. We show that differences in mean boldness are mostly affected by genetics. While the benthic treatment led to bolder individuals overall, the environmental effect was rather weak, suggesting that boldness lies under strong genetic influence with reduced plastic potential. Finally, we found hints of differences by morphs in boldness canalization through reduced variance and plasticity, and higher consistency in boldness within morphs. These findings provide new insights on how behavioural development may impact adaptive diversification.

Rapid and biased evolution of canalization during adaptive divergence revealed by dominance in gene expression variability during Arctic charr early development

Adaptive evolution may be influenced by canalization, the buffering of developmental processes from environmental and genetic perturbations, but how this occurs is poorly understood. Here, we explore how gene expression variability evolves in diverging and hybridizing populations, by focusing on the Arctic charr (Salvelinus alpinus) of Thingvallavatn, a classic case of divergence between feeding habitats. We report distinct profiles of gene expression variance for both coding RNAs and microRNAs between the offspring of two contrasting morphs (benthic/limnetic) and their hybrids reared in common conditions and sampled at two key points of cranial development. Gene expression variance in the hybrids is substantially affected by maternal effects, and many genes show biased expression variance toward the limnetic morph. This suggests that canalization, as inferred by gene expression variance, can rapidly diverge in sympatry through multiple gene pathways, which are associated with dominance patterns possibly biasing evolutionary trajectories and mitigating the effects of hybridization on adaptive evolution.

Variation in the genomic basis of parallel phenotypic and ecological divergence in benthic and pelagic morphs of Icelandic Arctic charr (Salvelinus alpinus)

Sympatric adaptive phenotypic divergence should be underlain by genomic differentiation between sub-populations. When divergence drives similar patterns of phenotypic and ecological variation within species we expect evolution to draw on common allelic variation. We investigated divergence histories and genomic signatures of adaptive divergence between benthic and pelagic morphs of Icelandic Arctic charr. Divergence histories for each of four populations were reconstructed using coalescent modelling and 14,187 single nucleotide polymorphisms. Sympatric divergence with continuous gene flow was supported in two populations while allopatric divergence with secondary contact was supported in one population; we could not differentiate between demographic models in the fourth population. We detected parallel patterns of phenotypic divergence along benthic-pelagic evolutionary trajectories among populations. Patterns of genomic differentiation between benthic and pelagic morphs were characterized by outlier loci in many narrow peaks of differentiation throughout the genome, which may reflect the eroding effects of gene flow on nearby neutral loci. We then used genome-wide association analyses to relate both phenotypic (body shape and size) and ecological (carbon and nitrogen stable isotopes) variation to patterns of genomic differentiation. Many peaks of genomic differentiation were associated with phenotypic and ecological variation in the three highly divergent populations, suggesting a genomic basis for adaptive divergence. We detected little evidence for a parallel genomic basis of differentiation as most regions and outlier loci were not shared among populations. Our results show that adaptive divergence can have varied genomic consequences in populations with relatively recent common origins, similar divergence histories, and parallel phenotypic divergence.

DNA methylation differences during development distinguish sympatric morphs of Arctic charr (Salvelinus alpinus)

Changes in DNA methylation in specific coding or non-coding regions can influence development and potentially divergence in traits within species and groups. While the impact of epigenetic variation on developmental pathways associated with evolutionary divergence is the focus of intense investigation, few studies have looked at recently diverged systems. Phenotypic diversity between closely related populations of Arctic charr (Salvelinus alpinus), which diverged within the last 10 000 years, offers an interesting ecological model to address such effects. Using bisulfite sequencing, we studied general DNA methylation patterns during development in the four sympatric morphs of Arctic charr from Lake Thingvallavatn. The data revealed strong differences between developmental timepoints and between morphs (mainly along the benthic - limnetic axis), both at single CpG sites and in 1,000bp-regions. Genes located close to differentially methylated CpG sites were involved in nucleosome assembly, regulation of osteoclast differentiation, and cell-matrix adhesion. Differentially methylated regions were enriched in tRNA and rRNA sequences, and half of them were located close to transcription start sites. The expression of 14 genes showing methylation differences over time or between morphs was further investigated by qPCR and nine of these were found to be differentially expressed between morphs. Four genes (ARHGEF37-like, H3-like, MPP3 and MEGF9) showed a correlation between methylation and expression. Lastly, histone gene clusters displayed interesting methylation differences between timepoints and morphs, as well as intragenic methylation variation. The results presented here provide a motivation for further studies on the contribution of epigenetic traits, such as DNA methylation, to phenotypic diversity and developmental mechanisms.

Genetic Causes and Consequences of Sympatric Morph Divergence in Salmonidae: A Search for Mechanisms

Repeatedly and recently evolved sympatric morphs exhibiting consistent phenotypic differences provide natural experimental replicates of speciation. Because such morphs are observed frequently in Salmonidae, this clade provides a rare opportunity to uncover the genomic mechanisms underpinning speciation. Such insight is also critical for conserving salmonid diversity, the loss of which could have significant ecological and economic consequences. Our review suggests that genetic differentiation among sympatric morphs is largely nonparallel apart from a few key genes that may be critical for consistently driving morph differentiation. We discuss alternative levels of parallelism likely underlying consistent morph differentiation and identify several factors that may temper this incipient speciation between sympatric morphs, including glacial history and contemporary selective pressures. Our synthesis demonstrates that salmonids are useful for studying speciation and poses additional research questions to be answered by future study of this family. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

The effect of growth rate on otolith-based discrimination of cod (Gadus morhua) ecotypes

Otolith shape has previously been used to identify ecotypes within the Icelandic cod (Gadus morhua) stock, using DST profiles to validate the results. Fish otolith shape variation has repeatedly been found to be largely determined by growth rate. To examine the effect of growth rate on the relationship between otolith shape and cod ecotypes (using the Pan I genotype as a proxy for ecotype), 826 archived sagittal otoliths collected over a 58 year sampling period were retrieved, the individual growth rate calculated, and otolith shape described using both Normalized Elliptic Fourier transform and Discrete Wavelet transform. Discriminant functions of otolith shape successfully classified ecotype, whether using Fourier or Wavelet descriptors, but only when excluding a heterozygous genotype from the analysis. The otolith shape variability of this genotype lowered the classification success, while otolith shape, in turn, was significantly affected by growth rate and cohort. Growth rate differences previously reported for the ecotypes were present, but were less marked than expected and indeed, growth rate variance attributable to ecotype identity was dwarfed by cohort- and location-related variance in growth. Such a strong effect of growth rate suggests that cod ecotype discrimination based on otolith shape is sensitive to both temporal and spatial variations in growth, which can mask the effect of ecotype-related growth rate differences on otolith shape.

Multivariate analysis of morphology, behaviour, growth and developmental timing in hybrids brings new insights into the divergence of sympatric Arctic charr morphs

BACKGROUND

Studying the development of fitness related traits in hybrids from populations diverging in sympatry is a fundamental approach to understand the processes of speciation. However, such traits are often affected by covariance structures that complicate the comprehension of these processes, especially because the interactive relationships between traits of different nature (e.g. morphology, behaviour, life-history) remain largely unknown in this context. In a common garden setup, we conducted an extensive examination of a large suit of traits putatively involved in the divergence of two morphs of Arctic charr (Salvelinus alpinus), and investigated the consequences of potential patterns of trait covariance on the phenotype of their hybrids. These traits were measured along ontogeny and involved growth, yolk sac resorption, developmental timing (hatching and the onset of exogeneous feeding), head morphology and feeding behaviour.

RESULTS

Growth trajectories provided the strongest signal of phenotypic divergence between the two charr. Strikingly, the first-generation hybrids did not show intermediate nor delayed growth but were similar to the smallest morph, suggesting parental biases in the inheritance of growth patterns. However, we did not observe extensive multivariate trait differences between the two morphs and their hybrids. Growth was linked to head morphology (suggesting that morphological variations in early juveniles relate to simple allometric effects) but this was the only strong signal of covariance observed between all the measured traits. Furthermore, we did not report evidence for differences in overall phenotypic variance between morphs, nor for enhanced phenotypic variability in their hybrids.

CONCLUSION

Our study shed light on the multivariate aspect of development in a context of adaptive divergence. The lack of evidence for the integration of most traits into a single covariance structure suggested that phenotypic constraints may not always favour nor impede divergence toward ecological niches differing in numerous physical and ecological variables, as observed in the respective habitats of the two charr. Likewise, the role of hybridization as a disruptive agent of trait covariance may not necessarily be significant in the evolution of populations undergoing resource polymorphism.

Genomic basis of deep-water adaptation in Arctic Charr (Salvelinus alpinus) morphs

The post-glacial colonization of Gander Lake in Newfoundland, Canada, by Arctic Charr (Salvelinus alpinus) provides the opportunity to study the genomic basis of adaptation to extreme deep-water environments. Colonization of deep-water (>50 m) habitats often requires extensive adaptation to cope with novel environmental challenges from high hydrostatic pressure, low temperature, and low light, but the genomic mechanisms underlying evolution in these environments are rarely known. Here, we compare genomic divergence between a deep-water morph adapted to depths of up to 288 m and a larger, piscivorous pelagic morph occupying shallower depths. Using both a SNP array and resequencing of whole nuclear and mitochondrial genomes, we find clear genetic divergence (F_ST  = 0.11-0.15) between deep and shallow water morphs, despite an absence of morph divergence across the mitochondrial genome. Outlier analyses identified many diverged genomic regions containing genes enriched for processes such as gene expression and DNA repair, cardiac function, and membrane transport. Detection of putative copy number variants (CNVs) uncovered 385 genes with CNVs distinct to piscivorous morphs, and 275 genes with CNVs distinct to deep-water morphs, enriched for processes associated with synapse assembly. Demographic analyses identified evidence for recent and local morph divergence, and ongoing reductions in diversity consistent with postglacial colonization. Together, these results show that Arctic Charr morph divergence has occurred through genome-wide differentiation and elevated divergence of genes underlying multiple cellular and physiological processes, providing insight into the genomic basis of adaptation in a deep-water habitat following postglacial recolonization.

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.

"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.