Invertebrate diversity in groundwater-filled lava caves is influenced by both neutral- and niche-based processes.

Understanding which factors shape and maintain biodiversity is essential to understand how ecosystems respond to crises. Biodiversity observed in ecological communities is a result of the interaction of various factors which can be classified as either neutral- or niche-based. The importance of these processes has been debated, but many scientists believe that both processes are important. Here, we use unique ecosystems in groundwater-filled lava caves near Lake Mývatn, to examine the importance of neutral- versus niche-based factors for shaping invertebrate communities. We studied diversity in benthic and epibenthic invertebrate communities and related them to ecological variables. We hypothesized that if neutral processes are the main drivers of community structure we would not see any clear relationship between the structure of community within caves and ecological factors. If niche-based processes are important we should see clear relationships between community structure and variation in ecological variables across caves. Both communities were species poor, with low densities of invertebrates, showing the resource limited and oligotrophic nature of these systems. Unusually for Icelandic freshwater ecosystems, the benthic communities were not dominated by Chironomidae (Diptera) larvae, but rather by crustaceans, mainly Cladocera. The epibenthic communities were not shaped by environmental variables, suggesting that they may have been structured primarily by neutral processes. The benthic communities were shaped by the availability of energy, and to some extent pH, suggesting that niche-based processes were important drivers of community structure, although neutral processes may still be relevant. The results suggest that both processes are important for invertebrate communities in freshwater, and research should focus on understanding both of these processes. The ponds we studied are representative of a number of freshwater ecosystems that are extremely vulnerable for human disturbance, making it even more important to understand how their biodiversity is shaped and maintained.

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

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

Fine scale diversity in the lava: genetic and phenotypic diversity in small populations of Arctic charr Salvelinus alpinus.

BACKGROUND: A major goal in evolutionary biology is to understand the processes underlying phenotypic variation in nature. Commonly, studies have focused on large interconnected populations or populations found along strong environmental gradients. However, studies on small fragmented populations can give strong insight into evolutionary processes in relation to discrete ecological factors. Evolution in small populations is believed to be dominated by stochastic processes, but recent work shows that small populations can also display adaptive phenotypic variation, through for example plasticity and rapid adaptive evolution. Such evolution takes place even though there are strong signs of historical bottlenecks and genetic drift. Here we studied 24 small populations of the freshwater fish Arctic charr (Salvelinus alpinus) found in groundwater filled lava caves. Those populations were found within a few km2-area with no apparent water connections between them. We studied the relative contribution of neutral versus non-neutral evolutionary processes in shaping phenotypic divergence, by contrasting patterns of phenotypic and neutral genetic divergence across populations in relation to environmental measurements. This allowed us to model the proportion of phenotypic variance explained by the environment, taking in to account the observed neutral genetic structure. RESULTS: These populations originated from the nearby Lake Mývatn, and showed small population sizes with low genetic diversity. Phenotypic variation was mostly correlated with neutral genetic diversity with only a small environmental effect. CONCLUSIONS: Phenotypic diversity in these cave populations appears to be largely the product of neutral processes, fitting the classical evolutionary expectations. However, the fact that neutral processes did not explain fully the phenotypic patterns suggests that further studies can increase our understanding on how neutral evolutionary processes can interact with other forces of selection at early stages of divergence. The accessibility of these populations has provided the opportunity for long-term monitoring of individual fish, allowing tracking how the environment can influence phenotypic and genetic divergence for shaping and maintaining diversity in small populations. Such studies are important, especially in freshwater, as habitat alteration is commonly breaking populations into smaller units, which may or may not be viable.

Experimental evidence for adaptive divergence in response to a warmed habitat reveals roles for morphology, allometry and parasite resistance.

Ectotherms are expected to be particularly vulnerable to climate change-driven increases in temperature. Understanding how populations adapt to novel thermal environments will be key for informing mitigation plans. We took advantage of threespine stickleback (Gasterosteus aculeatus) populations inhabiting adjacent geothermal (warm) and ambient (cold) habitats to test for adaptive evolutionary divergence using a field reciprocal transplant experiment. We found evidence for adaptive morphological divergence, as growth (length change) in non-native habitats related to head, posterior and total body shape. Higher growth in fish transplanted to a non-native habitat was associated with morphological shape closer to native fish. The consequences of transplantation were asymmetric with cold sourced fish transplanted to the warm habitat suffering from lower survival rates and greater parasite prevalence than warm sourced fish transplanted to the cold habitat. We also found divergent shape allometries that related to growth. Our findings suggest that wild populations can adapt quickly to thermal conditions, but immediate transitions to warmer conditions may be particularly difficult.

Demographic Basis of Spatially Structured Fluctuations in a Threespine Stickleback Metapopulation.

AbstractUncovering the demographic basis of population fluctuations is a central goal of population biology. This is particularly challenging for spatially structured populations, which require disentangling synchrony in demographic rates from coupling via movement between locations. In this study, we fit a stage-structured metapopulation model to a 29-year time series of threespine stickleback abundance in the heterogeneous and productive Lake Mývatn, Iceland. The lake comprises two basins (North and South) connected by a channel through which the stickleback disperse. The model includes time-varying demographic rates, allowing us to assess the potential contributions of recruitment and survival, spatial coupling via movement, and demographic transience to the population's large fluctuations in abundance. Our analyses indicate that recruitment was only modestly synchronized between the two basins, whereas survival probabilities of adults were more strongly synchronized, contributing to cyclic fluctuations in the lake-wide population size with a period of approximately 6 years. The analyses further show that the two basins were coupled through movement, with the North Basin subsidizing the South Basin and playing a dominant role in driving the lake-wide dynamics. Our results show that cyclic fluctuations of a metapopulation can be explained in terms of the combined effects of synchronized demographic rates and spatial coupling.

Geothermal stickleback populations prefer cool water despite multigenerational exposure to a warm environment.

Given the threat of climate change to biodiversity, a growing number of studies are investigating the potential for organisms to adapt to rising temperatures. Earlier work has predicted that physiological adaptation to climate change will be accompanied by a shift in temperature preferences, but empirical evidence for this is lacking. Here, we test whether exposure to different thermal environments has led to changes in preferred temperatures in the wild. Our study takes advantage of a "natural experiment" in Iceland, where freshwater populations of threespine sticklebacks (Gasterosteus aculeatus) are found in waters warmed by geothermal activity year-round (warm habitats), adjacent to populations in ambient-temperature lakes (cold habitats). We used a shuttle-box approach to measure temperature preferences of wild-caught sticklebacks from three warm-cold population pairs. Our prediction was that fish from warm habitats would prefer higher water temperatures than those from cold habitats. We found no support for this, as fish from both warm and cold habitats had an average preferred temperature of 13°C. Thus, our results challenge the assumption that there will be a shift in ectotherm temperature preferences in response to climate change. In addition, since warm-habitat fish can persist at relatively high temperatures despite a lower-temperature preference, we suggest that preferred temperature alone may be a poor indicator of a population's adaptive potential to a novel thermal environment.

Testing the predictability of morphological evolution in contrasting thermal environments.

Gaining the ability to predict population responses to climate change is a pressing concern. Using a "natural experiment," we show that testing for divergent evolution in wild populations from contrasting thermal environments provides a powerful approach, and likely an enhanced predictive power for responses to climate change. Specifically, we used a unique study system in Iceland, where freshwater populations of threespine sticklebacks (Gasterosteus aculeatus) are found in waters warmed by geothermal activity, adjacent to populations in ambient-temperature water. We focused on morphological traits across six pairs from warm and cold habitats. We found that fish from warm habitats tended to have a deeper mid-body, a subterminally orientated jaw, steeper craniofacial profile, and deeper caudal region relative to fish from cold habitats. Our common garden experiment showed that most of these differences were heritable. Population age did not appear to influence the magnitude or type of thermal divergence, but similar types of divergence between thermal habitats were more prevalent across allopatric than sympatric population pairs. These findings suggest that morphological divergence in response to thermal habitat, despite being relatively complex and multivariate, are predictable to a degree. Our data also suggest that the potential for migration of individuals between different thermal habitats may enhance nonparallel evolution and reduce our ability to predict responses to climate change.

A warmer environment can reduce sociability in an ectotherm.

The costs and benefits of being social vary with environmental conditions, so individuals must weigh the balance between these trade-offs in response to changes in the environment. Temperature is a salient environmental factor that may play a key role in altering the costs and benefits of sociality through its effects on food availability, predator abundance, and other ecological parameters. In ectotherms, changes in temperature also have direct effects on physiological traits linked to social behaviour, such as metabolic rate and locomotor performance. In light of climate change, it is therefore important to understand the potential effects of temperature on sociality. Here, we took the advantage of a 'natural experiment' of threespine sticklebacks from contrasting thermal environments in Iceland: geothermally warmed water bodies (warm habitats) and adjacent ambient-temperature water bodies (cold habitats) that were either linked (sympatric) or physically distinct (allopatric). We first measured the sociability of wild-caught adult fish from warm and cold habitats after acclimation to a low and a high temperature. At both acclimation temperatures, fish from the allopatric warm habitat were less social than those from the allopatric cold habitat, whereas fish from sympatric warm and cold habitats showed no differences in sociability. To determine whether differences in sociability between thermal habitats in the allopatric population were heritable, we used a common garden breeding design where individuals from the warm and the cold habitat were reared at a low or high temperature for two generations. We found that sociability was indeed heritable but also influenced by rearing temperature, suggesting that thermal conditions during early life can play an important role in influencing social behaviour in adulthood. By providing the first evidence for a causal effect of rearing temperature on social behaviour, our study provides novel insights into how a warming world may influence sociality in animal populations.

Variation in egg size and offspring phenotype among and within seven Arctic charr morphs.

Maternal effects have the potential to alter early developmental processes of offspring and contribute to adaptive diversification. Egg size is a major contributor to offspring phenotype, which can influence developmental trajectories and potential resource use. However, to what extent intraspecific variation in egg size facilitates evolution of resource polymorphism is poorly understood. We studied multiple resource morphs of Icelandic Arctic charr, ranging from an anadromous morph-with a phenotype similar to the proposed ancestral phenotype-to sympatric morphs that vary in their degree of phenotypic divergence from the ancestral anadromous morph. We characterized variation in egg size and tested whether egg size influenced offspring phenotype at early life stages (i.e., timing of- and size at- hatching and first feeding [FF]). We predicted that egg size would differ among morphs and be less variable as morphs diverge away from the ancestral anadromous phenotype. We also predicted that egg size would correlate with offspring size and developmental timing. We found morphs had different egg size, developmental timing, and size at hatching and FF. Egg size increased as phenotypic proximity to the ancestral anadromous phenotype decreased, with larger eggs generally giving rise to larger offspring, especially at FF, but egg size had no effect on developmental rate. The interaction between egg size and the environment may have a profound impact on offspring fitness, where the resulting differences in early life-history traits may act to initiate and/or maintain resource morphs diversification.

Invertebrate communities in springs across a gradient in thermal regimes.

In many respects, freshwater springs can be considered as unique ecosystems on the fringe of aquatic habitats. This integrates their uniqueness in terms of stability of environmental metrics. The main objective of our study was to evaluate how environmental variables may shape invertebrate diversity and community composition in different freshwater spring types and habitats within. In order to do so, we sampled invertebrates from 49 springs in Iceland, where we included both limnocrene and rheocrene springs. At each site, samples were taken from the benthic substrate of the spring ("surface") and the upwelling groundwater at the spring source ("source"). To collect invertebrates from the spring sources we used a modified method of "electrobugging" and Surber sampler for collecting invertebrates from the surface. In total, 54 invertebrate taxa were identified, mostly Chironomidae (Diptera). Chironomid larvae also dominated in terms of abundance (67%), followed by Ostracoda (12%) and Copepoda (9%). The species composition in the surface samples differed considerably between rheocrene and limnocrene springs and was characterised by several indicator species. Alpha diversity was greater at the surface of springs than at the source, but the beta diversity was higher at the source. Diversity, as summarized by taxa richness and Shannon diversity, was negatively correlated with temperature at the surface. At the source, on the other hand, Shannon diversity increased with temperature. The community assembly in springs appears to be greatly affected by water temperature, with the source community of hot springs being more niche-assembled (i.e., affected by mechanisms of tolerance and adaptation) than the source community of cold springs, which is more dispersal-assembled (i.e., by mechanisms of drift and colonization).

Variation in parasite resistance of Arctic charr, Salvelinus alpinus, between and within sympatric morphs.

Genetic variation in resistance against parasite infections is a predominant feature in host-parasite systems. However, mechanisms maintaining genetic polymorphism in resistance in natural host populations are generally poorly known. We explored whether differences in natural infection pressure between resource-based morphs of Arctic charr (Salvelinus alpinus) have resulted in differentiation in resistance profiles. We experimentally exposed offspring of two morphs from Lake Þingvallavatn (Iceland), the pelagic planktivorous charr ("murta") and the large benthivorous charr ("kuðungableikja"), to their common parasite, eye fluke Diplostomum baeri, infecting the eye humor. We found that there were no differences in resistance between the morphs, but clear differences among families within each morph. Moreover, we found suggestive evidence of resistance of offspring within families being positively correlated with the parasite load of the father, but not with that of the mother. Our results suggest that the inherited basis of parasite resistance in this system is likely to be related to variation among host individuals within each morph rather than ecological factors driving divergent resistance profiles at morph level. Overall, this may have implications for evolution of resistance through processes such as sexual selection.

Re-identification of individuals from images using spot constellations: a case study in Arctic charr (Salvelinus alpinus).

The ability to re-identify individuals is fundamental to the individual-based studies that are required to estimate many important ecological and evolutionary parameters in wild populations. Traditional methods of marking individuals and tracking them through time can be invasive and imperfect, which can affect these estimates and create uncertainties for population management. Here we present a photographic re-identification method that uses spot constellations in images to match specimens through time. Photographs of Arctic charr (Salvelinus alpinus) were used as a case study. Classical computer vision techniques were compared with new deep-learning techniques for masks and spot extraction. We found that a U-Net approach trained on a small set of human-annotated photographs performed substantially better than a baseline feature engineering approach. For matching the spot constellations, two algorithms were adapted, and, depending on whether a fully or semi-automated set-up is preferred, we show how either one or a combination of these algorithms can be implemented. Within our case study, our pipeline both successfully identified unmarked individuals from photographs alone and re-identified individuals that had lost tags, resulting in an approximately 4% increase in our estimate of survival rate. Overall, our multi-step pipeline involves little human supervision and could be applied to many organisms.

Differences among families in craniofacial shape at early life-stages of Arctic charr (Salvelinus alpinus).

BACKGROUND: Organismal fitness can be determined at early life-stages, but phenotypic variation at early life-stages is rarely considered in studies on evolutionary diversification. The trophic apparatus has been shown to contribute to sympatric resource-mediated divergence in several taxa. However, processes underlying diversification in trophic traits are poorly understood. Using phenotypically variable Icelandic Arctic charr (Salvelinus alpinus), we reared offspring from multiple families under standardized laboratory conditions and tested to what extent family (i.e. direct genetic and maternal effects) contributes to offspring morphology at hatching (H) and first feeding (FF). To understand the underlying mechanisms behind early life-stage variation in morphology, we examined how craniofacial shape varied according to family, offspring size, egg size and candidate gene expression. RESULTS: Craniofacial shape (i.e. the Meckel's cartilage and hyoid arch) was more variable between families than within families both across and within developmental stages. Differences in craniofacial morphology between developmental stages correlated with offspring size, whilst within developmental stages only shape at FF correlated with offspring size, as well as female mean egg size. Larger offspring and offspring from females with larger eggs consistently had a wider hyoid arch and contracted Meckel's cartilage in comparison to smaller offspring. CONCLUSIONS: This study provides evidence for family-level variation in early life-stage trophic morphology, indicating the potential for parental effects to facilitate resource polymorphism.

Multigenerational exposure to elevated temperatures leads to a reduction in standard metabolic rate in the wild.

In light of global climate change, there is a pressing need to understand and predict the capacity of populations to respond to rising temperatures. Metabolic rate is a key trait that is likely to influence the ability to cope with climate change. Yet, empirical and theoretical work on metabolic rate responses to temperature changes has so far produced mixed results and conflicting predictions.Our study addresses this issue using a novel approach of comparing fish populations in geothermally warmed lakes and adjacent ambient-temperature lakes in Iceland. This unique 'natural experiment' provides repeated and independent examples of populations experiencing contrasting thermal environments for many generations over a small geographic scale, thereby avoiding the confounding factors associated with latitudinal or elevational comparisons. Using Icelandic sticklebacks from three warm and three cold habitats, we measured individual metabolic rates across a range of acclimation temperatures to obtain reaction norms for each population.We found a general pattern for a lower standard metabolic rate (SMR) in sticklebacks from warm habitats when measured at a common temperature, as predicted by Krogh's rule. Metabolic rate differences between warm- and cold-habitat sticklebacks were more pronounced at more extreme acclimation temperatures, suggesting the release of cryptic genetic variation upon exposure to novel conditions, which can reveal hidden evolutionary potential. We also found a stronger divergence in metabolic rate between thermal habitats in allopatry than sympatry, indicating that gene flow may constrain physiological adaptation when dispersal between warm and cold habitats is possible.In sum, our study suggests that fish may diverge toward a lower SMR in a warming world, but this might depend on connectivity and gene flow between different thermal habitats. A free Plain Language Summary can be found within the Supporting Information of this article.

Repeatability of Adaptive Radiation Depends on Spatial Scale: Regional Versus Global Replicates of Stickleback in Lake Versus Stream Habitats.

The repeatability of adaptive radiation is expected to be scale-dependent, with determinism decreasing as greater spatial separation among "replicates" leads to their increased genetic and ecological independence. Threespine stickleback (Gasterosteus aculeatus) provide an opportunity to test whether this expectation holds for the early stages of adaptive radiation-their diversification in freshwater ecosystems has been replicated many times. To better understand the repeatability of that adaptive radiation, we examined the influence of geographic scale on levels of parallel evolution by quantifying phenotypic and genetic divergence between lake and stream stickleback pairs sampled at regional (Vancouver Island) and global (North America and Europe) scales. We measured phenotypes known to show lake-stream divergence and used reduced representation genome-wide sequencing to estimate genetic divergence. We assessed the scale dependence of parallel evolution by comparing effect sizes from multivariate models and also the direction and magnitude of lake-stream divergence vectors. At the phenotypic level, parallelism was greater at the regional than the global scale. At the genetic level, putative selected loci showed greater lake-stream parallelism at the regional than the global scale. Generally, the level of parallel evolution was low at both scales, except for some key univariate traits. Divergence vectors were often orthogonal, highlighting possible ecological and genetic constraints on parallel evolution at both scales. Overall, our results confirm that the repeatability of adaptive radiation decreases at increasing spatial scales. We suggest that greater environmental heterogeneity at larger scales imposes different selection regimes, thus generating lower repeatability of adaptive radiation at larger spatial scales.

Drift versus selection as drivers of phenotypic divergence at small spatial scales: The case of Belgjarskógur threespine stickleback.

Divergence in phenotypic traits is facilitated by a combination of natural selection, phenotypic plasticity, gene flow, and genetic drift, whereby the role of drift is expected to be particularly important in small and isolated populations. Separating the components of phenotypic divergence is notoriously difficult, particularly for multivariate phenotypes. Here, we assessed phenotypic divergence of threespine stickleback (Gasterosteus aculeatus) across 19 semi-interconnected ponds within a small geographic region (~7.5 km2) using comparisons of multivariate phenotypic divergence (PST), neutral genetic (FST), and environmental (EST) variation. We found phenotypic divergence across the ponds in a suite of functionally relevant phenotypic traits, including feeding, defense, and swimming traits, and body shape (geometric morphometric). Comparisons of PSTs with FSTs suggest that phenotypic divergence is predominantly driven by neutral processes or stabilizing selection, whereas phenotypic divergence in defensive traits is in accordance with divergent selection. Comparisons of population pairwise PSTs with ESTs suggest that phenotypic divergence in swimming traits is correlated with prey availability, whereas there were no clear associations between phenotypic divergence and environmental difference in the other phenotypic groups. Overall, our results suggest that phenotypic divergence of these small populations at small geographic scales is largely driven by neutral processes (gene flow, drift), although environmental determinants (natural selection or phenotypic plasticity) may play a role.

Invertebrate biodiversity in cold groundwater fissures in Iceland.

Iceland has an abundance of fissures that are parallel to the Mid-Atlantic Ridge where bedrock cracks as a result of continental rifting. Some fissures penetrate the aquifer and expose the groundwater within the bedrock, becoming springs. As such, groundwater fissures have uniform and constant physical and chemical environment but they can differ greatly in morphology. In addition, there is often great variation in depth within fissures and substrate types contrast between vertical rock wall and more heterogenous horizontal bottom. The variation in morphological environment may create dissimilar habitats with unique characteristics and/or influence distribution of resources. Our objective was to study macrozoobenthos communities in cold groundwater fissures in Iceland in relation to physical habitat by comparing invertebrate diversity and density both between fissures with different morphological characteristics as well as between substrate types and depths within fissures. Samples were collected in two fissures in SW Iceland, Silfra and Flosagjá. Assemblages were similar between fissures except for higher densities of cladocerans in Flosagjá fissure. Within fissures, there was significant difference in Shannon diversity between substrate types in Flosagjá, and ostracods were found in significantly higher densities on the bottom. The distribution of all other taxa groups was homogenous in both fissures regardless of depth gradient and substrate. Invertebrates were found to be living within and around a biofilm that covered the entire substrate. These biofilm mats are made from Cyanobacteria and benthic diatoms, which are successful under low light conditions and may minimize any effect of the heterogeneous habitat creating a uniform and suitable microhabitat for invertebrates regardless of depth and substrate type.

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.

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

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