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.

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.

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.

Bones in motion: Ontogeny of craniofacial development in sympatric arctic charr morphs.

BACKGROUND: The impressive diversity in the feeding apparatus often seen among related fish species clearly reflects differences in feeding modes and habitat utilization. Such variation can also be found within species. One example of such intraspecific diversity is the Arctic charr in Lake Thingvallavatn, where four distinct morphs coexist: two limnetic, with evenly protruding jaws, and two benthic, with subterminal lower jaws. We used these recently evolved morphs to study the role of ontogenetic variation in shaping craniofacial diversity. RESULTS: The segmental development of the pharyngeal arches and the order of events in craniofacial development is the same as has been described for teleosts, emphasizing the conserved nature of this process. However, our morphometric analyses reveal differences between morphs. Hatching is accompanied by increase in size and allometric shape changes in Arctic charr. Ontogenetic trajectories of craniofacial shape also differ significantly between morphs. CONCLUSIONS: The results point to hatching as a significant developmental event in Arctic charr and possibly other fishes. Also, the developmental origins of limnetic and benthic specializations in the craniofacial elements of Arctic charr may stem from events around hatching. This calls for investigations of the mechanisms and consequences of hatching, in the context of development and evolution. Developmental Dynamics 244:1168-1178, 2015. © 2015 Wiley Periodicals, Inc.