DNA methylation reprogramming in teleosts.

Early embryonic development is crucially important but also remarkably diverse among animal taxa. Axis formation and cell lineage specification occur due to both spatial and temporal control of gene expression. This complex system involves various signaling pathways and developmental genes such as transcription factors as well as other molecular interactants that maintain cellular states, including several types of epigenetic marks. 5mC DNA methylation, the chemical modification of cytosines in eukaryotes, represents one such mark. By influencing the compaction of chromatin (a high-order DNA structure), DNA methylation can either repress or induce transcriptional activity. Mammals exhibit a reprogramming of DNA methylation from the parental genomes in the zygote following fertilization, and later in primordial germ cells (PGCs). Whether these periods of methylation reprogramming are evolutionarily conserved, or an innovation in mammals, is an emerging question. Looking into these processes in other vertebrate lineages is thus important, and teleost fish, with their extensive species richness, phenotypic diversity, and multiple rounds of whole genome duplication, provide the perfect research playground for answering such a question. This review aims to present a concise state of the art of DNA methylation reprogramming in early development in fish by summarizing findings from different research groups investigating methylation reprogramming patterns in teleosts, while keeping in mind the ramifications of the methodology used, then comparing those patterns to reprogramming patterns in mammals.

Genetic structure and relatedness of brown trout (Salmo trutta) populations in the drainage basin of the Ölfusá river, South-Western Iceland.

Background: Lake Þingvallavatn in Iceland, a part of the river Ölfusá drainage basin, was presumably populated by brown trout soon after it formed at the end of the last Ice Age. The genetic relatedness of the brown trout in Þingvallavatn to other populations in the Ölfusá drainage basin is unknown. After the building of a dam at the outlet of the lake in 1959 brown trout catches declined, though numbers have now increased. The aim of this study was to assess effects of geographic isolation and potential downstream gene flow on the genetic structure and diversity in brown trout sampled in several locations in the western side of the watershed of River Ölfusá. We hypothesized that brown trout in Lake Þingvallavatn constituted several local spawning populations connected by occasional gene flow before the damming of the lake. We also estimated the effective population size (NE) of some of these populations and tested for signs of a recent population bottleneck in Lake Þingvallavatn. Methods: We sampled brown trout inhabiting four lakes and 12 rivers within and near the watershed of River Ölfusá by means of electro- and net- fishing. After stringent data filtering, 2,597 polymorphic loci obtained from ddRADseq data from 317 individuals were ascertained as putative neutral markers. Results: Overall, the genetic relatedness of brown trout in the Ölfusá watershed reflected the connectivity and topography of the waterways. Ancestry proportion analyses and a phylogenetic tree revealed seven distinct clusters, some of which corresponded to small populations with reduced genetic diversity. There was no evidence of downstream gene flow from Lake Þingvallavatn, although gene flow was observed from much smaller mountain populations. Most locations showed low NE values (i.e., ~14.6 on average) while the putative anadromous trout from River Sog and the spawning population from River Öxará, that flows into Lake Þingvallavatn, showed notably higher NE values (i.e., 71.2 and 56.5, respectively). No signals of recent population bottlenecks were detected in the brown trout of Lake Þingvallavatn. Discussion: This is the first time that the genetic structure and diversity of brown trout in the watershed of River Ölfusá have been assessed. Our results point towards the presence of a metapopulation in the watershed of Lake Þingvallavatn, which has been influenced by restoration efforts and is now dominated by a genetic component originated in River Öxará. Many of the locations studied represent different populations. Those that are isolated in headwater streams and lakes are genetically distinct presenting low genetic diversity, yet they can be important in increasing the genetic variation in downstream populations. These populations should be considered for conservation and direct management.

Rapid and biased evolution of canalization during adaptive divergence revealed by dominance in gene expression variability during Arctic charr early development.

Adaptive evolution may be influenced by canalization, the buffering of developmental processes from environmental and genetic perturbations, but how this occurs is poorly understood. Here, we explore how gene expression variability evolves in diverging and hybridizing populations, by focusing on the Arctic charr (Salvelinus alpinus) of Thingvallavatn, a classic case of divergence between feeding habitats. We report distinct profiles of gene expression variance for both coding RNAs and microRNAs between the offspring of two contrasting morphs (benthic/limnetic) and their hybrids reared in common conditions and sampled at two key points of cranial development. Gene expression variance in the hybrids is substantially affected by maternal effects, and many genes show biased expression variance toward the limnetic morph. This suggests that canalization, as inferred by gene expression variance, can rapidly diverge in sympatry through multiple gene pathways, which are associated with dominance patterns possibly biasing evolutionary trajectories and mitigating the effects of hybridization on adaptive evolution.

Insights into Early Ontogenesis of Salmo salar: RNA Extraction, Housekeeping Gene Validation and Transcriptional Expression of Important Primordial Germ Cell and Sex-Determination Genes.

The challenge in extracting high-quality RNA impedes the investigation of the transcriptome of developing salmonid embryos. Furthermore, the mRNA expression pattern of important PGC and SD genes during the initial embryonic development of Salmo salar is yet to be studied. So, in the present study, we aimed to isolate high-quality RNA from eggs and developing embryos to check vasa, dnd1, nanos3a, sdf1, gsdf, amh, cyp19a, dmrt1 and foxl2 expression by qPCR. Additionally, four HKGs (GAPDH, UB2L3, eEf1a and ß-actin) were validated to select the best internal control for qPCR. High-quality RNA was extracted, which was confirmed by spectrophotometer, agarose gel electrophoresis and Agilent TapeStation analysis. UB2L3 was chosen as a reference gene because it exhibited lower intra- and inter-sample variation. vasa transcripts were expressed in all the developmental stages, while dnd1 was expressed only up to 40 d°C. Nanos3a was expressed in later stages and remained at its peak for a shorter period, while sdf1 showed an irregular pattern of mRNA expression. The mRNA expression levels of SD genes were observed to be upregulated during the later stages of development, prior to hatching. This study presents a straightforward methodology for isolating high-quality RNA from salmon eggs, and the resulting transcript profiles of significant PGC and SD genes in S. salar could aid in improving our comprehension of reproductive development in this commercially important species.

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.

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.

DNA methylation differences during development distinguish sympatric morphs of Arctic charr (Salvelinus alpinus).

Changes in DNA methylation in specific coding or non-coding regions can influence development and potentially divergence in traits within species and groups. While the impact of epigenetic variation on developmental pathways associated with evolutionary divergence is the focus of intense investigation, few studies have looked at recently diverged systems. Phenotypic diversity between closely related populations of Arctic charr (Salvelinus alpinus), which diverged within the last 10,000 years, offers an interesting ecological model to address such effects. Using bisulphite sequencing, we studied general DNA methylation patterns during development in the four sympatric morphs of Arctic charr from Lake Thingvallavatn. The data revealed strong differences between developmental timepoints and between morphs (mainly along the benthic-limnetic axis), both at single CpG sites and in 1000 bp-regions. Genes located close to differentially methylated CpG sites were involved in nucleosome assembly, regulation of osteoclast differentiation, and cell-matrix adhesion. Differentially methylated regions were enriched in tRNA and rRNA sequences, and half of them were located close to transcription start sites. The expression of 14 genes showing methylation differences over time or between morphs was further investigated by qPCR and nine of these were found to be differentially expressed between morphs. Four genes (ARHGEF37-like, H3-like, MPP3 and MEGF9) showed a correlation between methylation and expression. Lastly, histone gene clusters displayed interesting methylation differences between timepoints and morphs, as well as intragenic methylation variation. The results presented here provide a motivation for further studies on the contribution of epigenetic traits, such as DNA methylation, to phenotypic diversity and developmental mechanisms.

RAD sequencing of common whelk, Buccinum undatum, reveals fine-scale population structuring in Europe and cryptic speciation within the North Atlantic.

Buccinum undatum is a subtidal gastropod that exhibits clear spatial variation in several phenotypic shell traits (color, shape, and thickness) across its North Atlantic distribution. Studies of spatial phenotypic variation exist for the species; however, population genetic studies have thus far relied on a limited set of mitochondrial and microsatellite markers. Here, we greatly expand on previous work by characterizing population genetic structure in B. undatum across the North Atlantic from SNP variation obtained by RAD sequencing. There was a high degree of genetic differentiation between Canadian and European populations (Iceland, Faroe Islands, and England) consistent with the divergence of populations in allopatry (F ST > 0.57 for all pairwise comparisons). In addition, B. undatum populations within Iceland, the Faroe Islands, and England are typified by weak but significant genetic structuring following an isolation-by-distance model. Finally, we established a significant correlation between genetic structuring in Iceland and two phenotypic traits: shell shape and color frequency. The works detailed here enhance our understanding of genetic structuring in B. undatum and establish the species as an intriguing model for future genome-wide association studies.

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.

Extensive genetic differentiation between recently evolved sympatric Arctic charr morphs.

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

A way forward with eco evo devo: an extended theory of resource polymorphism with postglacial fishes as model systems.

A major goal of evolutionary science is to understand how biological diversity is generated and altered. Despite considerable advances, we still have limited insight into how phenotypic variation arises and is sorted by natural selection. Here we argue that an integrated view, which merges ecology, evolution and developmental biology (eco evo devo) on an equal footing, is needed to understand the multifaceted role of the environment in simultaneously determining the development of the phenotype and the nature of the selective environment, and how organisms in turn affect the environment through eco evo and eco devo feedbacks. To illustrate the usefulness of an integrated eco evo devo perspective, we connect it with the theory of resource polymorphism (i.e. the phenotypic and genetic diversification that occurs in response to variation in available resources). In so doing, we highlight fishes from recently glaciated freshwater systems as exceptionally well-suited model systems for testing predictions of an eco evo devo framework in studies of diversification. Studies on these fishes show that intraspecific diversity can evolve rapidly, and that this process is jointly facilitated by (i) the availability of diverse environments promoting divergent natural selection; (ii) dynamic developmental processes sensitive to environmental and genetic signals; and (iii) eco evo and eco devo feedbacks influencing the selective and developmental environments of the phenotype. We highlight empirical examples and present a conceptual model for the generation of resource polymorphism - emphasizing eco evo devo, and identify current gaps in knowledge.

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

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

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

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.

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.

Inducers of salmon innate immunity: An in vitro and in vivo approach.

Maintaining fish health is one of the most important aims in aquaculture. Prevention of fish diseases therefore is crucial and can be achieved by various different strategies, including most often a combination of different methods such as optimal feed and fish density, as well as strengthening the immune system. Understanding the fish innate immune system and developing methods to activate it, in an effort to prevent infections in the first place, has been a goal in recent years. In this study we choose different inducers of the innate immune system and examined their effects in vitro on the salmon cell line CHSE-214. We found that the butyrate derivatives 4-phenyl butyrate (PBA) and ß-hydroxy-ß-methyl butyrate (HMB) induce the expression of various innate immune genes differentially over 24-72 h. Similarly, lipids generated from fish oils were found to have an effect on the expression of the antimicrobial peptides cathelicidin and hepcidin, as well as iNOS and the viral receptor RIG-1. Interestingly we found that vitamin D3, similar as in mammals, was able to increase cathelicidin expression in fish cells. The observed induction of these different innate immune factors correlated with antibacterial activity against Aeromonas salmonicida and antiviral activity against IPNV and ISAV in vitro. To relate this data to the in vivo situation we examined cathelicidin expression in juvenile salmon and found that salmon families vary greatly in their basal cathelicidin levels. Examining cathelicidin levels in families known to be resistant to IPNV showed that these QTL-families had lower basal levels of cathelicidin in gills, than non QTL-families. Feeding fish with HMB caused a robust increase in cathelicidin expression in gills, but not skin and this was independent of the fish being resistant to IPNV. These findings support the use of fish cell lines as a tool to develop new inducers of the fish innate immune system, but also highlight the importance of the tissue studied in vivo. Understanding the response of the innate immune system in different tissues and what effect this might have on infections and downstream cellular pathways is an interesting research topic for the future.

Investigation of the effects of estrogen on skeletal gene expression during zebrafish larval head development.

The development of craniofacial skeletal structures requires well-orchestrated tissue interactions controlled by distinct molecular signals. Disruptions in normal function of these molecular signals have been associated with a wide range of craniofacial malformations. A pathway mediated by estrogens is one of those molecular signals that plays role in formation of bone and cartilage including craniofacial skeletogenesis. Studies in zebrafish have shown that while higher concentrations of 17-ß estradiol (E 2) cause severe craniofacial defects, treatment with lower concentrations result in subtle changes in head morphology characterized with shorter snouts and flatter faces. The molecular basis for these morphological changes, particularly the subtle skeletal effects mediated by lower E 2 concentrations, remains unexplored. In the present study we address these effects at a molecular level by quantitative expression analysis of sets of candidate genes in developing heads of zebrafish larvae treated with two different E 2 concentrations. To this end, we first validated three suitable reference genes, ppia2, rpl8 and tbp, to permit sensitive quantitative real-time PCR analysis. Next, we profiled the expression of 28 skeletogenesis-associated genes that potentially respond to estrogen signals and play role in craniofacial development. We found E 2 mediated differential expression of genes involved in extracellular matrix (ECM) remodelling, mmp2/9/13, sparc and timp2a, as well as components of skeletogenic pathways, bmp2a, erf, ptch1/2, rankl, rarab and sfrp1a. Furthermore, we identified a co-expressed network of genes, including cpn1, dnajc3, esr1, lman1, rrbp1a, ssr1 and tram1 with a stronger inductive response to a lower dose of E 2 during larval head development.

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.

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.

The developmental transcriptome of contrasting Arctic charr ( Salvelinus alpinus) morphs.

Species and populations with parallel evolution of specific traits can help illuminate how predictable adaptations and divergence are at the molecular and developmental level. Following the last glacial period, dwarfism and specialized bottom feeding morphology evolved rapidly in several landlocked Arctic charr Salvelinus alpinus populations in Iceland.   To study the genetic divergence between small benthic morphs and limnetic morphs, we conducted RNA-sequencing charr embryos at four stages in early development. We studied two stocks with contrasting morphologies: the small benthic (SB) charr from Lake Thingvallavatn and Holar aquaculture (AC) charr. The data reveal significant differences in expression of several biological pathways during charr development. There was also an expression difference between SB- and AC-charr in genes involved in energy metabolism and blood coagulation genes. We confirmed differing expression of five genes in whole embryos with qPCR, including lysozyme and natterin-like which was previously identified as a fish-toxin of a lectin family that may be a putative immunopeptide. We also verified differential expression of 7 genes in the developing head that associated consistently with benthic v.s.limnetic morphology (studied in 4 morphs). Comparison of single nucleotide polymorphism (SNP) frequencies reveals extensive genetic differentiation between the SB and AC-charr (~1300 with more than 50% frequency difference). Curiously, three derived alleles in the otherwise conserved 12s and 16s mitochondrial ribosomal RNA genes are found in benthic charr. The data implicate multiple genes and molecular pathways in divergence of small benthic charr and/or the response of aquaculture charr to domestication. Functional, genetic and population genetic studies on more freshwater and anadromous populations are needed to confirm the specific loci and mutations relating to specific ecological traits in Arctic charr.

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

BACKGROUND: Understanding the molecular basis of craniofacial variation can provide insights into key developmental mechanisms of adaptive changes and their role in trophic divergence and speciation. Arctic charr (Salvelinus alpinus) is a polymorphic fish species, and, in Lake Thingvallavatn in Iceland, four sympatric morphs have evolved distinct craniofacial structures. We conducted a gene expression study on candidates from a conserved gene coexpression network, focusing on the development of craniofacial elements in embryos of two contrasting Arctic charr morphotypes (benthic and limnetic). RESULTS: Four Arctic charr morphs were studied: one limnetic and two benthic morphs from Lake Thingvallavatn and a limnetic reference aquaculture morph. The presence of morphological differences at developmental stages before the onset of feeding was verified by morphometric analysis. Following up on our previous findings that Mmp2 and Sparc were differentially expressed between morphotypes, we identified a network of genes with conserved coexpression across diverse vertebrate species. A comparative expression study of candidates from this network in developing heads of the four Arctic charr morphs verified the coexpression relationship of these genes and revealed distinct transcriptional dynamics strongly correlated with contrasting craniofacial morphologies (benthic versus limnetic). A literature review and Gene Ontology analysis indicated that a significant proportion of the network genes play a role in extracellular matrix organization and skeletogenesis, and motif enrichment analysis of conserved noncoding regions of network candidates predicted a handful of transcription factors, including Ap1 and Ets2, as potential regulators of the gene network. The expression of Ets2 itself was also found to associate with network gene expression. Genes linked to glucocorticoid signalling were also studied, as both Mmp2 and Sparc are responsive to this pathway. Among those, several transcriptional targets and upstream regulators showed differential expression between the contrasting morphotypes. Interestingly, although selected network genes showed overlapping expression patterns in situ and no morph differences, Timp2 expression patterns differed between morphs. CONCLUSION: Our comparative study of transcriptional dynamics in divergent craniofacial morphologies of Arctic charr revealed a conserved network of coexpressed genes sharing functional roles in structural morphogenesis. We also implicate transcriptional regulators of the network as targets for future functional studies.