Effects of population, family, and diet on craniofacial morphology of Icelandic Arctic charr (Salvelinus alpinus)

Complex and Dynamic Gene-by-Age and Gene-by-Environment Interactions Underlie Functional Morphological Variation in Adaptive Divergence in Arctic Charr (Salvelinus alpinus)

The evolution of adaptive phenotypic divergence requires heritable genetic variation. However, it is underappreciated that trait heritability is molded by developmental processes interacting with the environment. We hypothesized that the genetic architecture of divergent functional traits was dependent on age and foraging environment. Thus, we induced plasticity in full-sib families of Arctic charr (Salvelinus alpinus) morphs from two Icelandic lakes by mimicking prey variation in the wild. We characterized variation in body shape and size at two ages and investigated their genetic architecture with quantitative trait locus (QTL) analysis. Age had a greater effect on body shape than diet in most families, suggesting that development strongly influences phenotypic variation available for selection. Consistent with our hypothesis, multiple QTL were detected for all traits and their location depended on age and diet. Many of the genome-wide QTL were located within a subset of duplicated chromosomal regions suggesting that ancestral whole genome duplication events have played a role in the genetic control of functional morphological variation in the species. Moreover, the detection of two body shape QTL after controlling for the effects of age provides additional evidence for genetic variation in the plastic response of morphological traits to environmental variation. Thus, functional morphological traits involved in phenotypic divergence are molded by complex genetic interactions with development and environment.

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.

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.

The genetic basis of coordinated plasticity across functional units in a Lake Malawi cichlid mapping population

Adaptive radiations are often stereotypical, as populations repeatedly specialize along conserved environmental axes. Phenotypic plasticity may be similarly stereotypical, as individuals respond to environmental cues. These parallel patterns of variation, which are often consistent across traits, have led researchers to propose that plasticity can facilitate predictable patterns of evolution along environmental gradients. This "flexible stem" model of evolution raises questions about the genetic nature of plasticity, including how complex is the genetic basis for plasticity? Is plasticity across traits mediated by many distinct loci, or few "global" regulators? To address these questions, we reared a hybrid cichlid mapping population on alternate diet regimes mimicking an important environmental axis. We show that plasticity across an array of ecologically relevant traits is generally morphologically integrated, such that traits respond in a coordinated manner, especially those with overlapping function. Our genetic data are more ambiguous. While our mapping experiment provides little evidence for global genetic regulators of plasticity, these data do contain a genetic signal for the integration of plasticity across traits. Overall, our data suggest a compromise between genetic modularity, whereby plasticity may evolve independently across traits, and low level but widespread genetic integration, establishing the potential for plasticity to experience coordinated evolution.

A way forward with eco evo devo: an extended theory of resource polymorphism with postglacial fishes as model systems

A major goal of evolutionary science is to understand how biological diversity is generated and altered. Despite considerable advances, we still have limited insight into how phenotypic variation arises and is sorted by natural selection. Here we argue that an integrated view, which merges ecology, evolution and developmental biology (eco evo devo) on an equal footing, is needed to understand the multifaceted role of the environment in simultaneously determining the development of the phenotype and the nature of the selective environment, and how organisms in turn affect the environment through eco evo and eco devo feedbacks. To illustrate the usefulness of an integrated eco evo devo perspective, we connect it with the theory of resource polymorphism (i.e. the phenotypic and genetic diversification that occurs in response to variation in available resources). In so doing, we highlight fishes from recently glaciated freshwater systems as exceptionally well-suited model systems for testing predictions of an eco evo devo framework in studies of diversification. Studies on these fishes show that intraspecific diversity can evolve rapidly, and that this process is jointly facilitated by (i) the availability of diverse environments promoting divergent natural selection; (ii) dynamic developmental processes sensitive to environmental and genetic signals; and (iii) eco evo and eco devo feedbacks influencing the selective and developmental environments of the phenotype. We highlight empirical examples and present a conceptual model for the generation of resource polymorphism - emphasizing eco evo devo, and identify current gaps in knowledge.

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.

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.

Differential expression of the aryl hydrocarbon receptor pathway associates with craniofacial polymorphism in sympatric Arctic charr

Background

The developmental basis of craniofacial morphology hinges on interactions of numerous signalling systems. Extensive craniofacial variation in the polymorphic Arctic charr, a member of the salmonid family, from Lake Thingvallavatn (Iceland), offers opportunities to find and study such signalling pathways and their key regulators, thereby shedding light on the developmental pathways, and the genetics of trophic divergence.

Results

To identify genes involved in the craniofacial differences between benthic and limnetic Arctic charr, we used transcriptome data from different morphs, spanning early development, together with data on craniofacial expression patterns and skeletogenesis in model vertebrate species. Out of 20 genes identified, 7 showed lower gene expression in benthic than in limnetic charr morphs. We had previously identified a conserved gene network involved in extracellular matrix (ECM) organization and skeletogenesis, showing higher expression in developing craniofacial elements of benthic than in limnetic Arctic charr morphs. The present study adds a second set of genes constituting an expanded gene network with strong, benthic–limnetic differential expression. To identify putative upstream regulators, we performed knowledge-based motif enrichment analyses on the regulatory sequences of the identified genes which yielded potential binding sites for a set of known transcription factors (TFs). Of the 8 TFs that we examined using qPCR, two (Ahr2b and Ap2) were found to be differentially expressed between benthic and limnetic charr. Expression analysis of several known AhR targets indicated higher activity of the AhR pathway during craniofacial development in benthic charr morphotypes.

Conclusion

These results suggest a key role of the aryl hydrocarbon receptor (AhR) pathway in the observed craniofacial differences between distinct charr morphotypes.

Electronic supplementary material

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

Hidden genetic variation evolves with ecological specialization: the genetic basis of phenotypic plasticity in Arctic charr ecomorphs

The genetic variance that determines phenotypic variation can change across environments through developmental plasticity and in turn play a strong role in evolution. Induced changes in genotype-phenotype relationships should strongly influence adaptation by exposing different sets of heritable variation to selection under some conditions, while also hiding variation. Therefore, the heritable variation exposed or hidden from selection is likely to differ among habitats. We used ecomorphs from two divergent populations of Arctic charr (Salvelinus alpinus) to test the prediction that genotype-phenotype relationships would change in relation to environment. If present over several generations this should lead to divergence in genotype-phenotype relationships under common conditions, and to changes in the amount and type of hidden genetic variance that can evolve. We performed a common garden experiment whereby two ecomorphs from each of two Icelandic lakes were reared under conditions that mimicked benthic and limnetic prey to induce responses in craniofacial traits. Using microsatellite based genetic maps, we subsequently detected QTL related to these craniofacial traits. We found substantial changes in the number and type of QTL between diet treatments and evidence that novel diet treatments can in some cases provide a higher number of QTL. These findings suggest that selection on phenotypic variation, which is both genetically and environmentally determined, has shaped the genetic architecture of adaptive divergence in Arctic charr. However, while adaptive changes are occurring in the genome there also appears to be an accumulation of hidden genetic variation for loci not expressed in the contemporary environment.