Microevolutionary change in wild stickleback: Using integrative time-series data to infer responses to selection.

A central goal in evolutionary biology is to understand how different evolutionary processes cause trait change in wild populations. However, quantifying evolutionary change in the wild requires linking trait change to shifts in allele frequencies at causal loci. Nevertheless, datasets that allow for such tests are extremely rare and existing theoretical approaches poorly account for the evolutionary dynamics that likely occur in ecological settings. Using a decade-long integrative phenome-to-genome time-series dataset on wild threespine stickleback (Gasterosteus aculeatus), we identified how different modes of selection (directional, episodic, and balancing) drive microevolutionary change in correlated traits over time. Most strikingly, we show that feeding traits changed by as much 25% across 10 generations which was driven by changes in the genetic architecture (i.e., in both genomic breeding values and allele frequencies at genetic loci for feeding traits). Importantly, allele frequencies at genetic loci related to feeding traits changed at a rate greater than expected under drift, suggesting that the observed change was a result of directional selection. Allele frequency dynamics of loci related to swimming traits appeared to be under fluctuating selection evident in periodic population crashes in this system. Our results show that microevolutionary change in a wild population is characterized by different modes of selection acting simultaneously on different traits, which likely has important consequences for the evolution of correlated traits. Our study provides one of the most thorough descriptions to date of how microevolutionary processes result in trait change in a natural population.

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.

Genome-phenotype-environment associations identify signatures of selection in a panmictic population of threespine stickleback.

Adaptive genetic divergence occurs when selection imposed by the environment causes the genomic component of the phenotype to differentiate. However, genomic signatures of natural selection are usually identified without information on which trait is responding to selection by which selective agent(s). Here, we integrate whole-genome sequencing with phenomics and measures of putative selective agents to assess the extent of adaptive divergence in threespine stickleback occupying the highly heterogeneous lake Mývatn, NE Iceland. We find negligible genome wide divergence, yet multiple traits (body size, gill raker structure and defence traits) were divergent along known ecological gradients (temperature, predatory bird densities and water depth). SNP based heritability of all measured traits was high (h2  = 0.42-0.65), indicating adaptive potential for all traits. Environment-association analyses further identified thousands of loci putatively involved in selection, related to genes linked to, for instance, neuron development and protein phosphorylation. Finally, we found that loci linked to water depth were concurrently associated with pelvic spine length variation - supporting the conclusion that divergence in pelvic spine length occurred in the face of gene flow. Our results suggest that whilst there is substantial genetic variation in the traits measured, phenotypic divergence of Mývatn stickleback is mostly weakly associated with environmental gradients, potentially as a result of substantial gene flow. Our study illustrates the value of integrative studies that combine genomic assays of multivariate trait variation with landscape genomics.

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.

Steam caps in geothermal reservoirs can be monitored using seismic noise interferometry.

Harvesting geothermal energy often leads to a pressure drop in reservoirs, decreasing their profitability and promoting the formation of steam caps. While steam caps are valuable energy resources, they also alter the reservoir thermodynamics. Accurately measuring the steam fraction in reservoirs is essential for both operational and economic perspectives. However, steam content estimations are very limited both in space and time since current methods rely on direct measurements within production wells. Besides, these estimations normally present large uncertainties. Here, we present a pioneering method for indirectly sampling the steam content in the subsurface using the ever-present seismic background noise. We observe a consistent annual velocity drop in the Hengill geothermal field (Iceland) and establish a correlation between the velocity drop and steam buildup using in-situ borehole data. This application opens new avenues to track the evolution of any gas reservoir in the crust with a surface-based and cost-effective method.

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.

Microbial intrusion and seasonal dynamics in the groundwater microbiome of a porous basaltic rock aquifer used as municipal water reservoir.

Groundwater is a key resource for safe drinking water supply. Yet unconfined aquifers can be vulnerable to microbial contamination during extreme weather events that lead to surface runoff. The current study characterises the groundwater microbiome of a porous basaltic rock aquifer in South-West Iceland used for drinking water extraction and analyses the microbial community dynamics during surface runoff. The groundwater microbial community sampled from 12 wells across the extraction area contained over 745 prokaryotic genera and was phylogenetically similar between wells and most seasons, representing a diverse but homogenous ecosystem. The largest seasonal variation in the microbial community composition was detected during a period of concurrent snow melt and high precipitation leading to surface runoff. This period was characterised by an increased abundance of soil-associated taxa in the groundwater microbiome and specifically of taxa assigned to Aeromonas and Bacillus. A field experiment simulating high surface runoff around a groundwater well confirmed the increased abundance of surface soil microorganisms in the well water, indicating vulnerability of groundwater towards surface microbial intrusion during extreme weather events. As such events are likely to increase due to climate change, novel water management tools such as microbial community analysis could help ensure drinking water safety.

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.

Bacterial diversity in Icelandic cold spring sources and in relation to the groundwater amphipod Crangonyx islandicus.

Crangonyx islandicus is a groundwater amphipod endemic to Iceland, considered to have survived the Ice Ages in subglacial refugia. Currently the species is found in spring sources in lava fields along the tectonic plate boundary of the country. The discovery of a groundwater species in this inaccessible habitat indicates a hidden ecosystem possibly based on chemoautotrophic microorganisms as primary producers. To explore this spring ecosystem, we assessed its microbial diversity and analysed whether and how the diversity varied between the amphipods and the spring water, and if was dependent on environmental factors and geological settings. Isolated DNA from spring water and from amphipods was analysed using metabarcoding methods, targeting the 16S rRNA gene. Two genera of bacteria, Halomonas and Shewanella were dominating in the amphipod samples in terms of relative abundance, but not in the groundwater samples where Flavobacterium, Pseudomonas and Alkanindiges among others were dominating. The richness of the bacteria taxa in the microbial community of the groundwater spring sources was shaped by pH level and the beta diversity was shaped by geographic locations.