Shell morphology and color of the subtidal whelk Buccinum undatum exhibit fine-scaled spatial patterns

Publisher's version (útgefin grein)Geographical patterns in morphology can be the result of divergence among populations due to neutral or selective changes and/or phenotypic plasticity in response to different environments. Marine gastropods are ideal subjects on which to explore these patterns, by virtue of the remarkable intraspecific variation in life‐history traits and morphology often observed across relatively small spatial scales. The ubiquitous N‐Atlantic common whelk (Buccinum undatum) is well known for spatial variation in life‐history traits and morphology. Previous studies on genetic population structure have revealed that it exhibits significant differentiation across geographic distances. Within Breiðafjörður Bay, a large and shallow bay in W‐Iceland, genetic differentiation was demonstrated between whelks from sites separated by just 20 km. Here, we extended our previous studies on the common whelk in Breiðafjörður Bay by quantifying phenotypic variation in shell morphology and color throughout the Bay. We sought to test whether trait differentiation is dependent on geographic distance and/or environmental variability. Whelk in Breiðafjörður Bay displayed fine‐scale patterns of spatial variation in shape, thickness, and color diversity. Differentiation increased with increasing distance between populations, indicating that population connectivity is limited. Both shape and color varied along a gradient from the inner part of the bay in the east to the outer part in the west. Whelk shells in the innermost part of Breiðafjörður Bay were thick with an elongate shell, round aperture, and low color diversity, whereas in the outer part of the bay the shells were thinner, rounder, with a more elongate aperture and richer color diversity. Significant site‐specific difference in shell traits of the common whelk in correlation with environmental variables indicates the presence of local ecotypes and limited demographic connectivity.The project was funded by a research grant no. 141302051 from the Icelandic Centre for Research and a PhD grant from the University of Iceland Research Fund to Hildur Magnúsdóttir.Peer Reviewe

The mitochondrial genome of common whelk Buccinum undatum (Neogastropoda: Buccinidae)

Publisher's version (útgefin grein)The complete mitogenome of Buccinum undatum obtained from transcriptome analysis is presented here. The mitogenome is 15,265 bp in length and comprises 13 protein-coding genes, 2 ribosomal subunit genes, and 22 transfer RNAs. The complete genome was used to perform a phylogenetic analysis together with other Buccinoidea representatives with mitogenome data from GenBank.This work is supported by the Icelandic Research Council, under grant number 141302051.Peer Reviewe

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

Publisher's version (útgefin grein)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.We thank John Postlethwait for his valuable comments on the manuscript. This research was funded by the Icelandic Research Fund, Rannis (grant number 141360 to CAL et al., and grant number 173814–051 to SVB).Peer Reviewe

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.he project was funded by The Icelandic Center for Research (grant number: 100204011) to Sigurður S. Snorrason, Arnar Pálsson, Zophonías O. Jónsson and Bjarni K. Kristjánsson. The University of Iceland Doctoral Fund to Jóhannes Guðbrandsson and University of Iceland research fund to Arnar Pálsson, Sigurður S. Snorrason and Zophonías O. Jónsson. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Peer Reviewe

Extensive genetic differentiation between recently evolved sympatric Arctic charr morphs

Publisher's version (útgefin grein)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.This project was supported by The Icelandic Center for Research (RANNIS #100204011) to SSS and coworkers, The University of Iceland Doctoral Fund to JG and University of Iceland research fund to AP, SSS and ZOJ.Peer reviewe

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.This project was supported by The Icelandic Centre for Research (RANNIS/IRF, grant 100204) and The University of Iceland Research Fund.Peer Reviewe

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.This project was supported by The Icelandic Centre for Research (RANNIS/IRF, Grant 100204), The University of Iceland Research Fund and the Eimskip University Fund.Peer Reviewe

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 MATIN SMALL LETTER Y WITH ACUTEvatn, 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 (h(2) = 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 MATIN SMALL LETTER Y WITH ACUTEvatn 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.ISSN:0962-1083ISSN:1365-294