Morphological variation over ontogeny and environments in resource polymorphic 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.

The evolution of genomic organization through chromosomal rearrangements in Arctic charr (Salvelinus alpinus)

Chromosomal rearrangements (CRs) can play an important role in evolutionary diversification by preserving linkage among favourable alleles through reduced recombination and/or by reducing hybrid fitness due to genomic incompatibilities. Our goal was to determine to what extent CRs contribute to known patterns of genetic variation in Arctic charr (Salvelinus alpinus). To address this goal, we compared genetic linkage maps to identify whole arm CRs and smaller scale structural variants (SVs) such as translocations/transpositions and inversions found in groups of populations that reflect the temporal sequence of geographic isolation events. If CRs contribute to genetic differentiation, we expected that CRs would be specific to glacial lineages, geographic clusters of populations within lineages, and sympatric morphs. We detected fusions and fissions of whole chromosome arms and SV involving translocations/transpositions of the sex-determining gene (sdY) and inversions. Several CRs were shared across populations from the Arctic and Atlantic glacial lineages, Canadian and Icelandic populations within the Atlantic lineage, between two Icelandic populations and sympatric morphs within Icelandic populations, suggesting that their origin predates geographic isolation in glacial refugia. Other CRs were specific to single populations, which suggests a more recent origin of these variants in refugia, during post-glacial recolonization and/or in contemporary populations. Thus, CRs contribute relatively little to known patterns of genetic differentiation at different geographic scales but represent a pool of standing genetic variation for evolution.

Foraging-induced craniofacial plasticity is associated with an early, robust and dynamic transcriptional response

Phenotypic plasticity is the ability of a single genotype to vary its phenotype in response to the environment. Plasticity of the skeletal system in response to mechanical input is widely studied, but the timing of its transcriptional regulation is not well understood. Here, we used the cichlid feeding apparatus to examine the transcriptional dynamics of skeletal plasticity over time. Using three closely related species that vary in their ability to remodel bone and a panel of 11 genes, including well-studied skeletal differentiation markers and newly characterized environmentally sensitive genes, we examined plasticity at one, two, four and eight weeks following the onset of alternate foraging challenges. We found that the plastic species exhibited environment-specific bursts in gene expression beginning at one week, followed by a sharp decline in levels, while the species with more limited plasticity exhibited consistently low levels of gene expression. This trend held across nearly all genes, suggesting that it is a hallmark of the larger plasticity regulatory network. We conclude that plasticity of the cichlid feeding apparatus is not the result of slowly accumulating gene expression difference over time, but rather is stimulated by early bursts of environment-specific gene expression followed by a return to homeostatic levels.

Heat shocks during egg incubation led to developmental, morphological, and behavioral differences in Arctic charr (Salvelinus alpinus)

Temperature variation is affecting fish biodiversity worldwide, causing changes in geographic distribution, phenotypic structure, and even species extinction. Incubation is a critical stage for stenothermic species, which are vulnerable to large temperature fluctuations, and its effects on the phenotype at later developmental stages are understudied, despite the fact that the phenotype being essential for organism ecology and evolution. In this study, we tested the effects of heat shocks during the embryonic period on the phenotype of Arctic charr (Salvelinus alpinus). We repeatedly quantified multiple phenotypic traits, including morphology, development, and behavior, over a period of 4 months, from hatching to juvenile stage in individuals that had experienced heat shocks (+ 5°C on 24 h, seven times) during their embryonic stage and those that had not. We found that heat shocks led to smaller body size at hatching and a lower sociability. Interestingly, these effects weakened throughout the development of individuals and even reversed in the case of body size. We also found an accelerated growth rate and a higher body condition in the presence of heat shocks. Our study provides evidence that heat shocks experienced during incubation can have long-lasting effects on an individual's phenotype. This highlights the importance of the incubation phase for the development of ectothermic organisms and suggests that temperature fluctuations may have significant ecological and evolutionary implications for Arctic charr. Given the predicted increase in extreme events and the unpredictability of temperature fluctuations, it is critical to further investigate their effects on development by examining fluctuations that vary in frequency and intensity.

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 Dynamic Ontogenetic Shape Patterns of Adaptive Divergence and Sexual Dimorphism

The interplay between ecological diversification and sexual dimorphism has been largely overlooked in the literature. Sexually dimorphic species which are also undergoing adaptive radiations are ideal for filling this knowledge gap. The Arctic charr in lake Thingvallavatn is one such system: it is a sexually dimorphic species which has recently diverged along the benthic-limnetic ecological axis. In a long-running common-garden experiment we studied the shape variation throughout ontogeny of intra- and inter- morph crosses of benthic and limnetic charr from the lake. We found that shape differences between ecomorphs and sexes had a genetic component. Prior to the onset of sexual maturation, shape differences were attributable to cross type and were related to adaptations to benthic and limnetic niches, i.e., shorter lower jaws and rounder snouts in the benthic and evenly protruding snouts and pointier snouts in the limnetic. Reciprocal hybrids showed intermediate, transgressive and/or maternal morphologies. However, after the onset of sexual maturation larger morphological differences occurred between sexes than among cross types. Taken together, our results demonstrate that the interplay between ecological diversification and sexual dimorphism is complex and dynamic throughout ontogeny, and that long-term common garden experiments are immensely valuable for studying shape dynamics in different evolutionary scenarios.

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.

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.

Evolvability under climate change: Bone development and shape plasticity are heritable and correspond with performance in Arctic charr (Salvelinus alpinus)

Environmental conditions can impact the development of phenotypes and in turn the performance of individuals. Climate change, therefore, provides a pressing need to extend our understanding of how temperature will influence phenotypic variation. To address this, we assessed the impact of increased temperatures on ecologically significant phenotypic traits in Arctic charr (Salvelinus alpinus). We raised Arctic charr at 5°C and 9°C to simulate a predicted climate change scenario and examined temperature-induced variation in ossification, bone metabolism, skeletal morphology, and escape response. Fish reared at 9°C exhibited less cartilage and bone development at the same developmental stage, but also higher bone metabolism in localized regions. The higher temperature treatment also resulted in significant differences in craniofacial morphology, changes in the degree of variation, and fewer vertebrae. Both temperature regime and vertebral number affected escape response performance, with higher temperature leading to decreased latency. These findings demonstrate that climate change has the potential to impact development through multiple routes with the potential for plasticity and the release of cryptic genetic variation to have strong impacts on function through ecological performance and survival.

Developmental plasticity reveals hidden fish phenotypes and enables morphospace diversification

The establishment of a given phenotype is only one expression from a range of hidden developmental possibilities. Developmental plasticity at hidden reaction norms might elicit phenotypic diversification under new developmental environments. Current discussion benefits from empirical analyses that integrate multiple environmental stimuli to evaluate how plastic responses may shape phenotypic variation. We raised Megaleporinus macrocephalus fish in different environmental settings to address contributions of developmental plasticity for emergence of new phenotypes and subsequent morphospace diversification. Plastic morphotypes were evaluated at two complementary scales, the M. macrocephalus morphospace and the higher taxonomic level of Anostomidae family. Morphospace analyses demonstrated that developmental plasticity quickly releases distinct head morphotypes that were hidden in the parental monomorphic population. Plastic morphotypes occupied discrete and previously unfilled morphospace regions, a result obtained from comparisons with a control population and in analyses including several Anostomidae species. Plastic responses involved adjustments in shape and relative position of head bonesets, and fish raised under specific environmental combinations rescued phenotypic patterns described for different genera. Therefore, developmental plasticity possibly contributes to adaptive radiation in Anostomidae. Results illustrate how plastic responses enable morphospace diversification and contribute to 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.

Does phenotypic plasticity initiate developmental bias?

The generation of variation is paramount for the action of natural selection. Although biologists are now moving beyond the idea that random mutation provides the sole source of variation for adaptive evolution, we still assume that variation occurs randomly. In this review, we discuss an alternative view for how phenotypic plasticity, which has become well accepted as a source of phenotypic variation within evolutionary biology, can generate nonrandom variation. Although phenotypic plasticity is often defined as a property of a genotype, we argue that it needs to be considered more explicitly as a property of developmental systems involving more than the genotype. We provide examples of where plasticity could be initiating developmental bias, either through direct active responses to similar stimuli across populations or as the result of programmed variation within developmental systems. Such biased variation can echo past adaptations that reflect the evolutionary history of a lineage but can also serve to initiate evolution when environments change. Such adaptive programs can remain latent for millions of years and allow development to harbor an array of complex adaptations that can initiate new bouts of evolution. Specifically, we address how ideas such as the flexible stem hypothesis and cryptic genetic variation overlap, how modularity among traits can direct the outcomes of plasticity, and how the structure of developmental signaling pathways is limited to a few outcomes. We highlight key questions throughout and conclude by providing suggestions for future research that can address how plasticity initiates and harbors developmental bias.

Within-stream phenotypic divergence in head shape of brown trout associated with invasive brook trout

Competition with a non-native species can lead to morphological changes in native organisms induced by phenotypic plasticity, and by selection against individuals that do not adjust their morphology to the novel selection pressure. The morphological changes in native organisms are often associated with rapid behavioural responses to competition with the invader. However, knowledge of the interaction between the behaviour and morphology of native organisms competing with a non-native species remains scarce. Here, we investigated the effect of competition with non-native brook trout Salvelinus fontinalis on head shape of native brown trout Salmo trutta in a stream system where changes in diet and territorial behaviour of sympatric brown trout have previously been demonstrated. We found that sympatric brown trout had smaller eyes, shorter lower jaws and more terminal mouth than allopatric conspecifics. These differences in head shape were highly repeatable over a period of 12 months. Apparent survival indicated that the selection on head shape of brown trout was weaker in the sympatric than in the allopatric stretch of the stream. The results suggest that these changes reinforce divergences of foraging strategies between the allopatric and sympatric brown trout, which can negatively affect their population dynamics and trophic function in the food-web.

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.

Extensive genetic differentiation between recently evolved sympatric Arctic charr morphs

The availability of diverse ecological niches can promote adaptation of trophic specializations and related traits, as has been repeatedly observed in evolutionary radiations of freshwater fish. The role of genetics, environment, and history in ecologically driven divergence and adaptation, can be studied on adaptive radiations or populations showing ecological polymorphism. Salmonids, especially the Salvelinus genus, are renowned for both phenotypic diversity and polymorphism. Arctic charr (Salvelinus alpinus) invaded Icelandic streams during the glacial retreat (about 10,000 years ago) and exhibits many instances of sympatric polymorphism. Particularly, well studied are the four morphs in Lake Þingvallavatn in Iceland. The small benthic (SB), large benthic (LB), planktivorous (PL), and piscivorous (PI) charr differ in many regards, including size, form, and life history traits. To investigate relatedness and genomic differentiation between morphs, we identified variable sites from RNA-sequencing data from three of those morphs and verified 22 variants in population samples. The data reveal genetic differences between the morphs, with the two benthic morphs being more similar and the PL-charr more genetically different. The markers with high differentiation map to all linkage groups, suggesting ancient and pervasive genetic separation of these three morphs. Furthermore, GO analyses suggest differences in collagen metabolism, odontogenesis, and sensory systems between PL-charr and the benthic morphs. Genotyping in population samples from all four morphs confirms the genetic separation and indicates that the PI-charr are less genetically distinct than the other three morphs. The genetic separation of the other three morphs indicates certain degree of reproductive isolation. The extent of gene flow between the morphs and the nature of reproductive barriers between them remain to be elucidated.

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.

Gene expression in the phenotypically plastic Arctic charr (Salvelinus alpinus): A focus on growth and ossification at early stages of development

Gene expression during development shapes the phenotypes of individuals. Although embryonic gene expression can have lasting effects on developmental trajectories, few studies consider the role of maternal effects, such as egg size, on gene expression. Using qPCR, we characterize relative expression of 14 growth and/or skeletal promoting genes across embryonic development in Arctic charr (Salvelinus alpinus). We test to what extent their relative expression is correlated with egg size and size at early life‐stages within the study population. We predict smaller individuals to have higher expression of growth and skeletal promoting genes, due to less maternal resources (i.e., yolk) and prioritization of energy toward ossification. We found expression levels to vary across developmental stages and only three genes (Mmp9, Star, and Sgk1) correlated with individual size at a given developmental stage. Contrary to our hypothesis, expression of Mmp9 and Star showed a non‐linear relationship with size (at post fertilization and hatching, respectively), whilst Sgk1 was higher in larger embryos at hatching. Interestingly, these genes are also associated with craniofacial divergence of Arctic charr morphs. Our results indicate that early life‐stage variation in gene expression, concomitant to maternal effects, can influence developmental plasticity and potentially the evolution of resource polymorphism in fishes.

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.

Allometric trajectories of body and head morphology in three sympatric Arctic charr (Salvelinus alpinus (L.)) morphs

A study of body and head development in three sympatric reproductively isolated Arctic charr (Salvelinus alpinus (L.)) morphs from a subarctic lake (Skogsfjordvatn, northern Norway) revealed allometric trajectories that resulted in morphological differences. The three morphs were ecologically assigned to a littoral omnivore, a profundal benthivore and a profundal piscivore, and this was confirmed by genetic analyses (microsatellites). Principal component analysis was used to identify the variables responsible for most of the morphological variation of the body and head shape. The littoral omnivore and the profundal piscivore morph had convergent allometric trajectories for the most important head shape variables, developing bigger mouths and relatively smaller eyes with increasing head size. The two profundal morphs shared common trajectories for the variables explaining most of the body and head shape variation, namely head size relative to body size, placement of the dorsal and pelvic fins, eye size and mouth size. In contrast, the littoral omnivore and the profundal benthivore morphs were not on common allometric trajectories for any of the examined variables. The findings suggest that different selective pressures could have been working on traits related to their trophic niche such as habitat and diet utilization of the three morphs, with the two profundal morphs experiencing almost identical environmental conditions.

Niche partitioning and the role of intraspecific niche variation in structuring a guild of generalist anurans

Intra-population niche differences in generalist foragers have captured the interest of ecologists, because such individuality can have important ecological and evolutionary implications. Few researchers have investigated how these differences affect the relationships among ecologically similar, sympatric species. Using stable isotopes, stomach contents, morphology and habitat preference, we examined niche partitioning within a group of five anurans and determined whether variation within species could facilitate resource partitioning. Species partitioned their niches by trophic level and by foraging habitat. However, there was considerable intraspecific variation in trophic level, with larger individuals generally feeding at higher trophic levels. For species at intermediate trophic levels, smaller individuals overlapped in trophic level with individuals of smaller species and larger individuals overlapped with the smallest individuals from larger species. Species varied in carbon isotopes; species with enriched carbon isotope ratios foraged farther from ponds, whereas species with depleted carbon isotope values foraged closer to ponds. Our study shows that these species partition their niches by feeding at different trophic levels and foraging at different distances from ponds. The intraspecific variation in trophic level decreased the number of individuals from each species that overlapped in trophic level with individuals from other species, which can facilitate species coexistence.

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.

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.

Transcriptional dynamics of a conserved gene expression network associated with craniofacial divergence in Arctic charr

Background

Understanding the molecular basis of craniofacial variation can provide insights into key developmental mechanisms of adaptive changes and their role in trophic divergence and speciation. Arctic charr (Salvelinus alpinus) is a polymorphic fish species, and, in Lake Thingvallavatn in Iceland, four sympatric morphs have evolved distinct craniofacial structures. We conducted a gene expression study on candidates from a conserved gene coexpression network, focusing on the development of craniofacial elements in embryos of two contrasting Arctic charr morphotypes (benthic and limnetic).

Results

Four Arctic charr morphs were studied: one limnetic and two benthic morphs from Lake Thingvallavatn and a limnetic reference aquaculture morph. The presence of morphological differences at developmental stages before the onset of feeding was verified by morphometric analysis. Following up on our previous findings that Mmp2 and Sparc were differentially expressed between morphotypes, we identified a network of genes with conserved coexpression across diverse vertebrate species. A comparative expression study of candidates from this network in developing heads of the four Arctic charr morphs verified the coexpression relationship of these genes and revealed distinct transcriptional dynamics strongly correlated with contrasting craniofacial morphologies (benthic versus limnetic). A literature review and Gene Ontology analysis indicated that a significant proportion of the network genes play a role in extracellular matrix organization and skeletogenesis, and motif enrichment analysis of conserved noncoding regions of network candidates predicted a handful of transcription factors, including Ap1 and Ets2, as potential regulators of the gene network. The expression of Ets2 itself was also found to associate with network gene expression. Genes linked to glucocorticoid signalling were also studied, as both Mmp2 and Sparc are responsive to this pathway. Among those, several transcriptional targets and upstream regulators showed differential expression between the contrasting morphotypes. Interestingly, although selected network genes showed overlapping expression patterns in situ and no morph differences, Timp2 expression patterns differed between morphs.

Conclusion

Our comparative study of transcriptional dynamics in divergent craniofacial morphologies of Arctic charr revealed a conserved network of coexpressed genes sharing functional roles in structural morphogenesis. We also implicate transcriptional regulators of the network as targets for future functional studies.

Electronic supplementary material

The online version of this article (doi:10.1186/2041-9139-5-40) 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.

Wnt signalling underlies the evolution of new phenotypes and craniofacial variability in Lake Malawi cichlids

Progress towards understanding adaptive radiations at the mechanistic level is still limited with regard to the proximate molecular factors that both promote and constrain evolution. Here we focus on the craniofacial skeleton and show that expanded Wnt/β-catenin signalling early in ontogeny is associated with the evolution of phenotypic novelty and ecological opportunity in Lake Malawi cichlids. We demonstrate that the mode of action of this molecular change is to effectively lock into place an early larval phenotype, likely through accelerated rates of bone deposition. However, we demonstrate further that this change toward phenotypic novelty may in turn constrain evolutionary potential through the corresponding reduction in craniofacial plasticity at later stages of ontogeny. In all, our data implicate the Wnt pathway as an important mediator of craniofacial form and offer new insights into how developmental systems can evolve to both promote and constrain evolutionary change.

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.

Validation of reference genes for expression studies during craniofacial development in arctic charr

Arctic charr (Salvelinus alpinus) is a highly polymorphic species and in Lake Thingvallavatn, Iceland, four phenotypic morphs have evolved. These differences in morphology, especially in craniofacial structures are already apparent during embryonic development, indicating that genes important in the formation of the craniofacial features are expressed differentially between the morphs. In order to generate tools to examine these expression differences in Arctic charr, the aim of the present study was to identify reference genes for quantitative real-time PCR (qPCR). The specific aim was to select reference genes which are able to detect very small expression differences among different morphs. We selected twelve candidate reference genes from the literature, identified corresponding charr sequences using data derived from transcriptome sequencing (RNA-seq) and examined their expression using qPCR. Many of the candidate reference genes were found to be stably expressed, yet their quality-rank as reference genes varied considerably depending on the type of analysis used. In addition to commonly used software for reference gene validation, we used classical statistics to evaluate expression profiles avoiding a bias for reference genes with similar expression patterns (co-regulation). Based on these analyses we chose three reference genes, ACTB, UB2L3 and IF5A1 for further evaluation. Their consistency was assessed in an expression study of three known craniofacially expressed genes, sparc (or osteonectin), matrix metalloprotease 2 (mmp2) and sox9 (sex-determining region Y box 9 protein) using qPCR in embryo heads derived from four charr groups at three developmental time points. The three reference genes were found to be very suitable for studying expression differences between the morphotypes, enabling robust detection of small relative expression changes during charr development. Further, the results showed that sparc and mmp2 are differentially expressed in embryos of different Arctic charr morphotypes.

Behavioral, ecological and genetic differentiation in an open environment--a study of a mysid population in the Baltic Sea

Diel vertical migration (DVM) is often assumed to encompass an entire population. However, bimodal nighttime vertical distributions have been observed in various taxa. Mysid shrimp populations also display this pattern with one group concentrated in the pelagia and the other near the bottom. This may indicate alternative migratory strategies, resembling the seasonal partial migrations seen in birds, fishes and amphibians, where only a subset of the population migrates. To assess the persistence of these alternative strategies, we analyzed the nitrogen and carbon stable isotope signatures (as proxies for diet), biochemical indices (as proxies for growth condition), and genetic population divergence in the Baltic mysid Mysis salemaai collected at night in the pelagia and close to the bottom. Stable isotope signatures were significantly different between migrants (pelagic samples) and residents (benthic samples), indicating persistent diet differences, with pelagic mysids having a more uniform and carnivorous diet. Sequencing of the mitochondrial cytochrome subunit I (COI) gene showed genetic differentiation attributable to geographic location but not between benthic and pelagic groups. Divergent migration strategies were however supported by significantly lower gene flow between benthic populations indicating that these groups have a lower predisposition for horizontal migrations compared to pelagic ones. Different migration strategies did not convey measurable growth benefits as pelagic and benthic mysids had similar growth condition indices. Thus, the combination of ecological, biochemical and genetic markers indicate that this partial migration may be a plastic behavioral trait that yields equal growth benefits.

Phenotypic plasticity, heterochrony and ontogenetic repatterning during juvenile development of divergent Arctic charr (Salvelinus alpinus)

Phenotypic plasticity is a developmental process that plays a role as a source of variation for evolution. Models of adaptive divergence make the prediction that increasing ecological specialization should be associated with lower levels of plasticity. We tested for differences in the magnitude, rate and trajectory of morphological plasticity in two lake populations of Arctic charr (Salvelinus alpinus) that exhibited variation in the degree of resource polymorphism. We reared offspring on diet treatments that mimicked benthic and pelagic prey. Offspring from the more divergent population had lower levels of morphological plasticity. Allometry influenced the rate of shape change over ontogeny, with differences in rate among ecomorphs being minimal when allometric variation was removed. However, plasticity in the spatial trajectory of development was extensive across ecomorphs, both with and without the inclusion of allometric variation, suggesting that different aspects of shape development can evolve independently.

Evolution of adaptive diversity and genetic connectivity in Arctic charr (Salvelinus alpinus) in Iceland

The ecological theory of adaptive radiation predicts that the evolution of phenotypic diversity within species is generated by divergent natural selection arising from different environments and competition between species. Genetic connectivity among populations is likely also to have an important role in both the origin and maintenance of adaptive genetic diversity. Our goal was to evaluate the potential roles of genetic connectivity and natural selection in the maintenance of adaptive phenotypic differences among morphs of Arctic charr, Salvelinus alpinus, in Iceland. At a large spatial scale, we tested the predictive power of geographic structure and phenotypic variation for patterns of neutral genetic variation among populations throughout Iceland. At a smaller scale, we evaluated the genetic differentiation between two morphs in Lake Thingvallavatn relative to historically explicit, coalescent-based null models of the evolutionary history of these lineages. At the large spatial scale, populations are highly differentiated, but weakly structured, both geographically and with respect to patterns of phenotypic variation. At the intralacustrine scale, we observe modest genetic differentiation between two morphs, but this level of differentiation is nonetheless consistent with strong reproductive isolation throughout the Holocene. Rather than a result of the homogenizing effect of gene flow in a system at migration-drift equilibrium, the modest level of genetic differentiation could equally be a result of slow neutral divergence by drift in large populations. We conclude that contemporary and recent patterns of restricted gene flow have been highly conducive to the evolution and maintenance of adaptive genetic variation in Icelandic Arctic charr.