Search

This is the peer reviewed version of the following article: Aykanat, T., Rasmussen, M., Ozerov, M., Niemelä, E., Paulin, L., Vähä, J.-P., ...Primmer, C. R. (2020). Life‐history genomic regions explain differences in Atlantic salmon marine diet specialization. Journal of Animal Ecology, ?(?), 1-15., which has been published in final form at https://doi.org/10.1111/1365-2656.13324. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Use of Self-Archived Versions. 1. Animals employ various foraging strategies along their ontogeny to acquire energy, and with varying degree of efficiencies, to support growth, maturation and subsequent reproduction events. Individuals that can efficiently acquire energy early are more likely to mature at an earlier age, as a result of faster energy gain which can fuel maturation and reproduction. 2. We aimed to test the hypothesis that heritable resource acquisition variation that covaries with efficiency along the ontogeny would influence maturation timing of individuals. 3. To test this hypothesis, we utilized Atlantic salmon as a model which exhibits a simple, hence trackable, genetic control of maturation age. We then monitored the variation in diet acquisition (quantified as stomach fullness and composition) of individuals with different ages, and linked it with genomic regions (haploblocks) that were previously identified to be associated with age‐at‐maturity. 4. Consistent with the hypothesis, we demonstrated that one of the life‐history genomic regions tested (six6) was indeed associated with age‐dependent differences in stomach fullness. Prey composition was marginally linked to six6, and suggestively (but non‐significantly) to vgll3 genomic regions. We further showed Atlantic salmon switched to the so‐called ‘feast and famine’ strategy along the ontogeny, where older age groups exhibited heavier stomach content, but that came at the expense of running on empty more often. 5. These results suggest genetic variation underlying resource utilization may explain the genetic basis of age structure in Atlantic salmon. Given that ontogenetic diet has a genetic component and the strong spatial diversity associated with these genomic regions, we predict populations with diverse maturation age will have diverse evolutionary responses to future changes in marine food web structures.

The role of female choice in sexual selection is well established, including the recognition that females choose their mates based on multiple cues. These cues may include intrinsic aspects of a male's phenotype as well as aspects of the environment associated with the male. The role of the spatial location of a potential mate has been well studied in territorial vertebrates. However, despite their role as laboratory models for studies of sexual selection, the potential for insects to choose their mates on the basis of location has scarcely been studied. We studied a natural population of individually tagged crickets (Gryllus campestris) in a meadow in Northern Spain. Adults typically move between burrows every few days, allowing us to examine how pairing success of males can be predicted by the burrow they occupy, independent of their own characteristics. We observed the entirety of ten independent breeding seasons to provide replication and to determine whether the relative importance of these factors is stable across years. We find that both male ID and the ID his burrow affect the likelihood that he is paired with a female, but the burrow has a consistently greater influence. Furthermore, the two factors interact: the relative attractiveness of an individual male depends on which burrow he occupies. Our finding demonstrates a close interaction between naturally and sexually selected traits. It also demonstrates that mate choice studies may benefit from considering not only obvious secondary sexual traits, but also more cryptic traits such as microhabitat choice. We show that female insects choosing to cohabit with a male place more importance on where he lives than on who he is, but the combination of the two is more important still. We know that female birds often choose a male based on the quality of his territory; our 10 years of observing crickets moving around a Spanish meadow to share burrows with members of the opposite sex, reveals insects can do the same.

A beta-class carbonic anhydrase (CA, EC 4.2.1.1) was cloned from the genome of the Monogenean platyhelminth Gyrodactylus salaris, a parasite of Atlantic salmon. The new enzyme, GsaCA beta has a significant catalytic activity for the physiological reaction, CO2 + H2O (sic) HCO3- + H+ with a k(cat) of 1.1 x 10(5) s(-1) and a k(cat)/K-m of 7.58 x 10(6) M-1 x s(-1). This activity was inhibited by acetazolamide (K-I of 0.46 mu M), a sulphonamide in clinical use, as well as by selected inorganic anions and small molecules. Most tested anions inhibited GsaCA beta at millimolar concentrations, but sulfamide (K-I of 81 mu M), N,N-diethyldithiocarbamate (K-I of 67 mu M) and sulphamic acid (K-I of 6.2 mu M) showed a rather efficient inhibitory action. There are currently very few non-toxic agents effective in combating this parasite. GsaCA beta is subsequently proposed as a new drug target for which effective inhibitors can be designed.

Organisms utilize varying lipid resource allocation strategies as a means to survive seasonal environmental changes and life-history stage transitions. In Atlantic salmon, a certain lipid threshold is needed to initiate sexual maturation. Because of this, an individual’s maturation schedule may be affected by changes in temperature and food availability across the seasons that create natural fluctuations of lipid reserves.Recent studies have found a genome region, including the genevgll3,that explains a large proportion of variation for size and age at maturity. Vgll3encodes a transcription co-factor that acts as an inhibitor of adipogenesis in mice and also affects condition factor and other phenotypes in juvenile salmon. However, even with many studies investigating varying temperature effects, there is a lack of temporal studies examining the effects of seasonality on such phenotypes, nor have the effects ofvgll3genotype on condition factor and maturation in different temperatures at different life stages.Here, we investigate the influence of different larval and juvenile incubation temperatures,vgll3genotype and their interactions on juvenile salmon phenotypes including body condition, and sexual maturation rate. We reared Atlantic salmon for 2 years in varying temperatures with an average 1.76 °C difference between warm and cold treatments in four different larval-juvenile phase treatment groups (Warm-Warm, Warm-Cold, Cold-Warm, and Cold-Cold) until the first occurrence of maturation in males.We found no effect of larval temperature on the measured phenotypes or maturation rate, suggesting the occurrence of growth compensation over the course of the experiment. Agreeing with previous studies, an increased maturation rate was observed in individuals of the warm juvenile temperature treatment.In addition, we observed differences in condition factor associated withvgll3genotype, wherebyvgll3*EEindividuals (the genotype associated with early maturation) had a less variable condition factor across the seasons compared to thevgll3*LL(associated with late maturation) individuals.This result suggests avgll3influence on resource acquisition and allocation strategies, possibly linked with the early maturation process, with individuals carrying the early maturationvgll3genotype having a higher early maturation rate and a higher condition factor in the spring.

A β-class carbonic anhydrase (CA, EC 4.2.1.1) present in the genome of the Monogenean platyhelminth Gyrodactylus salaris, a fish parasite, GsaCAβ, has been investigated for its inhibitory effects with a panel of sulphonamides and sulfamates, some of which in clinical use. Several effective GsaCAβ inhibitors were identified, belonging to simple heterocyclic sulphonamides, the deacetylated precursors of acetazolamide and methazolamide (KIsof 81.9-139.7 nM). Many other simple benezene sulphonamides and clinically used agents, such as acetazolamide, methazolamide, ethoxzolamide, dorzolamide, benzolamide, sulthiame and hydrochlorothiazide showed inhibition constants

Age at maturity can vary dramatically within some species. We show that, when raised in controlled conditions, Atlantic salmon age at maturity is affected greatly by temperature as well as a recently discovered age-at-maturity gene. The influence of temperature on age at maturity differed between two populations, but the effect of the age-at-maturity gene was similar in both temperatures and populations.Age at maturity is a key life history trait involving a trade-off between survival risk and reproductive investment, and is an important factor for population structures. In ectotherms, a warming environment may have a dramatic influence on development and life history, but this influence may differ between populations. While an increasing number of studies have examined population-dependent reactions with temperature, few have investigated this in the context of maturation timing. Atlantic salmon, a species of high conservation relevance, is a good study species for this topic as it displays considerable variation in age at maturity, of which a large proportion has been associated with a genomic region including the strong candidate gene vgll3. Until now, the effect of this gene in the context of different environments and populations has not been studied. Using a large-scale common-garden experiment, we find strong effects of temperature, population-of-origin, and vgll3 genotype on maturation in 2-year-old male Atlantic salmon (Salmo salar). With a temperature difference of 1.8 degrees C, maturation probability was 4.8 times higher in the warm treatment than the cold treatment. This temperature effect was population-specific and was higher in the southern (60.48 degrees N) compared to the northern (65.01 degrees N) population. The early maturation vgll3*E allele was associated with a significantly higher maturation probability, but there was no vgll3 interaction with temperature or population. Both body condition and body mass associated with maturation. The body mass association was only present in the warm treatment. Our findings demonstrate that (i) populations can vary in their response to temperature change in terms of age at maturity, (ii) high intrinsic growth could be associated with higher thermal sensitivity for life history variation and (iii) vgll3 effects on age at maturity might be similar between populations and different thermal environments.

Many post-smolt tissues were collected and screened as part of the EU funded SALSEA-Merge (Project No. 212529) part of the Seventh Framework Programme and co-sponsored by the Atlantic Salmon Trust and the Total Foundation., Gilbey, J., Utne, K.R., Wennevik, V., Beck, A.C., Kausrud, K., Hindar, K., Garcia de Leaniz, C., Cherbonnel, C., Coughlan, J., Cross, T.F., Dillane, E., Ensing, D., García-Vázquez, E., Hole, L.R., Holm, M., Holst, J.C., Jacobsen, J.A., Jensen, A.J., Karlsson, S., Ó Maoiléidigh, N., Mork, K.A., Nielsen, E.E., Nøttestad, L., Primmer, C.R., Prodöhl, P., Prusov, S., Stevens, J.R., Thomas, K., Whelan, K., McGinnity, P., Verspoor, E.

One of the most well-known life-history continuums is the fast–slow axis, where ‘fast’ individuals mature earlier than ‘slow’ individuals. ‘Fast’ individuals are predicted to be more active than ‘slow’ individuals because high activity is required to maintain a fast life-history strategy. Recent meta-analyses revealed mixed evidence for such integration. Here, we test whether known life-history genotypes differ in activity expression by using Atlantic salmon ( Salmo salar ) as a model. In salmon, variation in Vgll3, a transcription cofactor, explains approximately 40% of variation in maturation timing. We predicted that the allele related to early maturation ( vgll3 *E) would be associated with higher activity. We used an automated surveillance system to follow approximately 1900 juveniles including both migrants and non-migrants (i.e. smolt and parr fish, respectively) in semi-natural conditions over 31 days (approx. 580 000 activity measurements). In migrants, but not in non-migrants, vgll3 explained variation in activity according to our prediction in a sex-dependent manner. Specifically, in females the vgll3 *E allele was related to increasing activity, whereas in males the vgll3 *L allele (later maturation allele) was related to increasing activity. These sex-dependent effects might be a mechanism maintaining within-population genetic life-history variation.

Publisher Copyright: Copyright © 2022 the Author(s). Plasticity can allow organisms to maintain consistent performance across a wide range of environmental conditions. However, it remains largely unknown how costly plasticity is and whether a trade-off exists between plasticity and performance under optimal conditions. Biological rates generally increase with temperature, and to counter that effect, fish use physiological plasticity to adjust their biochemical and physiological functions. Zebrafish in the wild encounter large daily and seasonal temperature fluctuations, suggesting they should display high physiological plasticity. Conversely, laboratory zebrafish have been at optimal temperatures with low thermal fluctuations for over 150 generations. We treated this domestication as an evolution experiment and asked whether this has reduced the physiological plasticity of laboratory fish compared to their wild counterparts. We measured a diverse range of phenotypic traits, from gene expression through physiology to behavior, in wild and laboratory zebrafish acclimated to 15 temperatures from 10 °C to 38 °C. We show that adaptation to the laboratory environment has had major effects on all levels of biology. Laboratory fish show reduced plasticity and are thus less able to counter the direct effects of temperature on key traits like metabolic rates and thermal tolerance, and this difference is detectable down to gene expression level. Rapid selection for faster growth in stable laboratory environments appears to have carried with it a trade-off against physiological plasticity in captive zebrafish compared with their wild counterparts.

We tested the possibility that vgll3, a gene linked with maturation age in Atlantic salmon (Salmo salar), may be associated with behaviour by measuring aggressiveness and feeding activity in 380 juveniles with different vgll3 genotypes. Contrary to our prediction, individuals with the genotype associated with later maturation (vgll3*LL) were significantly more aggressive than individuals with the genotype associated with earlier maturation (vgll3*EE). Individuals with higher aggression were also significantly lighter in colour and had higher feeding activity. Although higher aggression was associated with higher feeding activity, there was no association between feeding activity and vgll3 genotype. Increased aggression of vgll3*LL individuals was independent of their sex and size, and genotypes did not differ in their condition factor. These results imply that aggressive behaviour may have an energetic cost impairing growth and condition, especially when food cannot be monopolized. This may have implications for individual fitness and aquaculture practices.

Genetic diversity within species and populations is necessary for long-term survival and thus constitutes a key component of preserving biodiversity, but until now, it has rarely been integrated into conservation policies. Ancient and historical genetic data [ancient/historical DNA (a/hDNA)], such as those from specimens stored in natural history collections, can add a temporal dimension to conservation genetic inferences by providing baseline levels of diversity that contemporary data can be compared with and help guide conservation actions. To increase the use and impact of a/hDNA research in preserving biodiversity, genetic indicators must be explicitly included in conservation policies, the benefits and limitations of using a/hDNA need to be clearly communicated to all conservation actors, and relationships between academics, museums, conservation practitioners, and policy makers must be strengthened.

Understanding the drivers of evolution is a fundamental aim in biology. However, identifying the evolutionary impacts of human activities is challenging because of a lack of temporal data and limited knowledge of the genetic basis of most traits. Here, we identify the drivers of evolution toward maturity at an earlier age in Atlantic salmon through two types of fisheries-induced evolution acting in opposing directions: an indirect effect linked with harvest of a salmon prey species (capelin) at sea (selection against late maturation) and a direct effect due to net fishing in rivers (selection against early maturation). Because capelin are harvested as an aquaculture feed protein source, we hereby determine an indirect path by which salmon aquaculture may influence wild salmon populations.

Our ability to examine genetic variation across entire genomes has enabled many studies searching for the genetic basis of local adaptation. These studies have identified numerous loci as candidates for differential local selection; however, relatively few have examined the overlap among candidate loci identified from independent studies of the same species in different geographic areas or evolutionary lineages. We used an allelotyping approach with a 220K SNP array to characterize the population genetic structure of Atlantic salmon in north-eastern Europe and ask whether the same genomic segments emerged as outliers among populations in different geographic regions. Genome-wide data recapitulated the phylogeographic structure previously inferred from mtDNA and microsatellite markers. Independent analyses of three genetically and geographically distinct groups of populations repeatedly inferred the same 17 haploblocks to contain loci under differential local selection. The most strongly supported of these replicated haploblocks had known strong associations with life-history variation or immune response in Atlantic salmon. Our results are consistent with these genomic segments harbouring large-effect loci which have a major role in Atlantic salmon diversification and are ideal targets for validation studies.

Anadromous salmonid fishes frequently exhibit strong geographic population structuring. However, population genetic differentiation of Atlantic salmon (Salmo salar) at fine geographic scales differs across equivalent spatial extents in different regions. So far, fine-scale genetic differentiation has not been assessed in rivers of the Baltic Sea, a region that contains an evolutionarily distinct Atlantic salmon lineage. Thus, Baltic salmon are currently managed on the river level, without focus on potential genetic structure and diversity within rivers. Here, we used microsatellites to characterize the genetic structure of wild juvenile salmon sampled throughout the interconnected, northern Baltic Tornio and Kalix Rivers. We found genetic differentiation within the two rivers, but not between them: salmon in the upper reaches differed from individuals in the lower reaches, regardless of river system. Further, examining smolts migrating from the river to the sea and adults returning from the sea to spawn, we found an association between the genetic structure and seasonal migration timing. Out-migrating smolts genetically assigned to upper river reaches were older and tended to reach the sea later in the season than smolts from the lower reaches. In contrast, mature adults originating from the upper reaches returned to the river early in the season. Our observation of genetic population structuring between downstream and upstream reaches of the large Tornio and Kalix rivers, and its association with migration timing, implies that careful temporal management of the northern Baltic fisheries would help to preserve the diversity and sustainability of the wild salmon stocks of these rivers.

Temperature is a key environmental parameter affecting both the phenotypes and distributions of organisms, particularly ectotherms. Rapid organismal responses to thermal environmental changes have been described for several ectotherms; however, the underlying molecular mechanisms often remain unclear. Here, we studied whole genome cytosine methylation patterns of European grayling (Thymallus thymallus) embryos from five populations with contemporary adaptations of early life history traits at either 'colder' or 'warmer' spawning grounds. We reared fish embryos in a common garden experiment using two temperatures that resembled the 'colder' and 'warmer' conditions of the natal natural environments. Genome-wide methylation patterns were similar in populations originating from colder thermal origin subpopulations, whereas single nucleotide polymorphisms uncovered from the same data identified strong population structure among isolated populations, but limited structure among interconnected populations. This was surprising because the previously studied gene expression response among populations was mostly plastic, and mainly influenced by the developmental temperature. These findings support the hypothesis of the magnified role of epigenetic mechanisms in modulating plasticity. The abundance of consistently changing methylation loci between two warmer-to-colder thermal origin population pairs suggests that local adaptation has shaped the observed methylation patterns. The dynamic nature of the methylomes was further highlighted by genome-wide and site-specific plastic responses. Our findings support both the presence of a plastic response in a subset of CpG loci, and the evolutionary role of methylation divergence between populations adapting to contrasting thermal environments.

Sex chromosomes contribute substantially to key evolutionary processes such as speciation and adaptation. Several theories suggest that evolution could occur more rapidly on sex chromosomes, but currently our understanding of whether and how this occurs is limited. Here, we present an analysis of the great tit (Parus major) genome, aiming to detect signals of faster-Z evolution. We find mixed evidence of faster divergence on the Z chromosome than autosomes, with significantly higher divergence being found in ancestral repeats, but not at 4- or 0-fold degenerate sites. Interestingly, some 4-fold sites appear to be selectively constrained, which may mislead analyses that use these sites as the neutral reference (e.g., dN/dS). Consistent with other studies in birds, the mutation rate is significantly higher in males than females, and the long-term Z-to-autosome effective population size ratio is only 0.5, significantly lower than the expected value of 0.75. These are indicative of male-driven evolution and high variance in male reproductive success, respectively. We find no evidence for an increased efficacy of positive selection on the Z chromosome. In contrast, the Z chromosome in great tits appears to be affected by increased genetic drift, which has led to detectable signals of weakened intensity of purifying selection. These results provide further evidence that the Z chromosome often has a low effective population size, and that this has important consequences for its evolution. They also highlight the importance of considering multiple factors that can affect the rate of evolution and effective population sizes of sex chromosomes.

Progress in genome sequencing now enables the large-scale generation of reference genomes. Various international initiatives aim to generate reference genomes representing global biodiversity. These genomes provide unique insights into genomic diversity and architecture, thereby enabling comprehensive analyses of population and functional genomics, and are expected to revolutionize conservation genomics.

BackgroundThe Glanville fritillary (Melitaea cinxia) butterfly is a model system for metapopulation dynamics research in fragmented landscapes. Here, we provide a chromosome-level assembly of the butterfly's genome produced from Pacific Biosciences sequencing of a pool of males, combined with a linkage map from population crosses.ResultsThe final assembly size of 484 Mb is an increase of 94 Mb on the previously published genome. Estimation of the completeness of the genome with BUSCO indicates that the genome contains 92–94% of the BUSCO genes in complete and single copies. We predicted 14,810 genes using the MAKER pipeline and manually curated 1,232 of these gene models.ConclusionsThe genome and its annotated gene models are a valuable resource for future comparative genomics, molecular biology, transcriptome, and genetics studies on this species.

Characterizing the role of different mutational effect sizes in the evolution of fitness-related traits has been a major goal in evolutionary biology for a century. Such characterization in a diversity of systems, both model and non-model, will help to understand the genetic processes underlying fitness variation. However, well-characterized genetic architectures of such traits in wild populations remain uncommon. In this study, we used haplotype-based and multi-SNP Bayesian association methods with sequencing data for 313 individuals from wild populations to test the mutational composition of known candidate regions for sea age at maturation in Atlantic salmon (Salmo salar). We detected an association at five loci out of 116 candidates previously identified in an aquaculture strain with maturation timing in wild Atlantic salmon. We found that at four of these five loci, variation explained by the locus was predominantly driven by a single SNP suggesting the genetic architecture of this trait includes multiple loci with simple, non-clustered alleles and a locus with potentially more complex alleles. This highlights the diversity of genetic architectures that can exist for fitness-related traits. Furthermore, this study provides a useful multi-SNP framework for future work using sequencing data to characterize genetic variation underlying phenotypes in wild populations.

To mitigate effects of climate change it is important to understand species’ responses to increasing temperatures. This has often been done by studying survival or activity at temperature extremes. Before such extremes are reached, however, effects on fertility may already be apparent.Sex differences in the thermal sensitivity of fertility (TSF) could impact species persistence under climate warming because female fertility is typically more limiting to population growth than male fertility. However, little is known about sex differences in TSF.Here we first demonstrate that the mating system can strongly influence TSF using the seed beetle Callosobruchus maculatus. We exposed populations carrying artificially induced mutations to two generations of short-term experimental evolution under alternative mating systems, manipulating the opportunity for natural and sexual selection on the mutations. We then measured TSF in males and females subjected to juvenile or adult heat stress.Populations kept under natural and sexual selection had higher fitness, but similar TSF, compared to control populations kept under relaxed selection. However, females had higher TSF than males, and strikingly, this sex difference had increased over only two generations in populations evolving under sexual selection.We hypothesized that an increase in male-induced harm to females during mating had played a central role in driving this evolved sex difference, and indeed, remating under conditions limiting male harassment of females reduced both male and female TSF. Moreover, we show that manipulation of mating system parameters in C. maculatus generates intraspecific variation in the sex difference in TSF equal to that found among a diverse set of studies on insects.Our study provides a causal link between the mating system and TSF. Sexual conflict, (re)mating rates, and genetic responses to sexual selection differ among ecological settings, mating systems and species. Our study therefore also provides mechanistic understanding for the variability in previously reported TSFs which can inform future experimental assays and predictions of species responses to climate warming.

Efforts to understand the genetic underpinnings of phenotypic variation are becoming more and more frequent in molecular ecology. Such efforts often lead to the identification of candidate regions showing signals of association and/or selection. These regions may contain multiple genes and therefore validation of which genes are actually responsible for the signal is required. In Atlantic salmon (Salmo salar), a large-effect locus for maturation timing, an ecologically important trait, occurs in a genomic region including two genes, vgll3 and akap11, but data for clearly determining which of the genes (or both) contribute to the association have been lacking. Here, we take advantage of natural recombination events detected between the two candidate genes in a salmon broodstock to reduce linkage disequilibrium at the locus, thus enabling delineation of the influence of variation at these two genes on early maturation. By rearing 5,895 males to maturation age, of which 81% had recombinant vgll3/akap11 allelic combinations, we found that vgll3 single nucleotide polymorphism (SNP) variation was strongly associated with early maturation, whereas there was little or no association between akap11 SNP variation and early maturation. These findings provide strong evidence supporting vgll3 as the primary candidate gene in the chromosome 25 locus for influencing early maturation. This will help guide future research for understanding the genetic processes controlling early maturation. This also exemplifies the utility of natural recombinants to more precisely map causal variation underlying ecologically important phenotypic diversity.

Behavioural responses are often the first reaction of an organism to human-induced rapid environmental change (HIREC), yet current empirical evidence provides no consensus about the main environmental features that animals respond to behaviourally or which behaviours are responsive to HIREC. To understand how changes in behaviour can be predicted by different forms of HIREC, we conducted a meta-analysis of the existing empirical literature focusing on behavioural responses to five axes of environmental change (climate change, changes in CO2, direct human impact, changes in nutrients and biotic exchanges) in five behavioural domains (aggression, exploration, activity, boldness and sociability) across a range of taxa but with a focus on fish and bird species. Our meta-analysis revealed a general absence of directional behavioural responses to HIREC. However, the absolute magnitude of the effect sizes was large. This means that animals have strong behavioural responses to HIREC, but the responses are not clearly in any particular direction. Moreover, the absolute magnitude of the effect sizes differed between different behaviours and different forms of HIREC: Exploration responded more strongly than activity, and climate change induced the strongest behavioural responses. Model heterogeneities identified that effect sizes varied primarily because of study design, and the specific sample of individuals used in a study; phylogeny also explains significant variation in our bird model. Based on these results, we make four recommendations to further our understanding: 1) a more balanced representation of laboratory and field studies, 2) consideration of context dependency, 3) standardisation of the methods and definitions used to quantify and study behaviours and 4) consideration of the role for individual differences in behaviour. © 2021 Nordic Society Oikos. Published by John Wiley & Sons Ltd

Over the past decades, Atlantic salmon (Salmo salar, Salmonidae) has emerged as a model system for sexual maturation research, owing to the high diversity of life history strategies, knowledge of trait genetic architecture, and their high economic value. The aim of this synthesis is to summarize the current state of knowledge concerning maturation in Atlantic salmon, outline knowledge gaps, and provide a roadmap for future work. We summarize the current state of knowledge: 1) maturation in Atlantic salmon takes place over the entire life cycle, starting as early as embryo development, 2) variation in the timing of maturation promotes diversity in life history strategies, 3) ecological and genetic factors influence maturation, 4) maturation processes are sex-specific and may have fitness consequences for each sex, 5) genomic studies have identified large-effect loci that influence maturation, 6) the brain-pituitary-gonadal axis regulates molecular and physiological processes of maturation, 7) maturation is a key component of fisheries, aquaculture, conservation, and management, and 8) climate change, fishing pressure, and other anthropogenic stressors likely have major effects on salmon maturation. In the future, maturation research should focus on a broader diversity of life history stages, including early embryonic development, the marine phase and return migration. We recommend studies combining ecological and genetic approaches will help disentangle the relative contributions of effects in different life history stages to maturation. Functional validation of large-effect loci should reveal how these genes influence maturation. Finally, continued research in maturation will improve our predictions concerning how salmon may adapt to fisheries, climate change, and other future challenges.

Age at maturity is a major contributor to the diversity of life history strategies in organisms. The process of maturation is influenced by both genetics and the environment, and includes changes in levels of sex hormones and behavior, but the specific factors leading to variation in maturation timing are poorly understood. gnrh1 regulates the transcription of gonadotropin genes at pubertal onset in many species, but this gene is lacking in certain teleost species including Atlantic salmon (Salmo salar), which raises the possibility of the involvement of other important regulatory factors during this process. Earlier research has reported a strong association of alternative alleles of the vgll3 gene with maturation timing in Atlantic salmon, suggesting it as a potential candidate regulating reproductive axis genes. Here, we investigated the expression of reproductive axis genes in one-year-old Atlantic salmon males with immature gonads and different vgll3 genotypes during the spawning period. We detected strong vgll3 genotype-dependent differential expression of reproductive axis genes (such as fshb, lhb, amh and igf3) tested in the pituitary, and testis. In addition, we observed differential expression of jun (ap1) and nr5a1b (sf1), potential upstream regulators of gonadotropins in the pituitary, as well as axin2, id3, insl3, itch, ptgs2a and ptger4b, the downstream targets of amh and igf3 in the testis. Hereby, we provide evidence of strong vgll3 genotype-dependent transcriptional regulation of reproductive axis genes prior to sexual maturation and suggest alternative models for distinct actions of vgll3 genotypes on the related molecular processes.

A better understanding of the genetic and phenotypic architecture underlying life-history variation is a longstanding aim in biology. Theories suggest energy metabolism determines life-history variation by modulating resource acquisition and allocation trade-offs, but the genetic underpinnings of the relationship and its dependence on ecological conditions have rarely been demonstrated. The strong genetic determination of age-at-maturity by two unlinked genomic regions ( vgll3 and six6 ) makes Atlantic salmon ( Salmo salar ) an ideal model to address these questions. Using more than 250 juveniles in common garden conditions, we quantified the covariation between metabolic phenotypes—standard and maximum metabolic rates (SMR and MMR), and aerobic scope (AS)—and the life-history genomic regions, and tested if food availability modulates the relationships. We found that the early maturation genotype in vgll3 was associated with higher MMR and consequently AS. Additionally, MMR exhibited physiological epistasis; it was decreased when late maturation genotypes co-occurred in both genomic regions. Contrary to our expectation, the life-history genotypes had no effects on SMR. Furthermore, food availability had no effect on the genetic covariation, suggesting a lack of genotype-by-environment interactions. Our results provide insights on the key organismal processes that link energy use at the juvenile stage to age-at-maturity, indicating potential mechanisms by which metabolism and life-history can coevolve.

Compared to angiosperms, gymnosperms lag behind in the availability of assembled and annotated genomes. Most genomic analyses in gymnosperms, especially conifer tree species, rely on the use of de novo assembled transcriptomes. However, the level of allelic redundancy and transcript fragmentation in these assembled transcriptomes, and their effect on downstream applications have not been fully investigated. Here, we assessed three assembly strategies, including the utility of haploid (megagametophyte) tissue during de novo assembly as single-allele guides, for six individuals and five different tissues in Pinus sylvestris. We then contrasted haploid and diploid tissue genotype calls obtained from the assembled transcriptomes to evaluate the extent of paralog mapping. The use of the haploid tissue during assembly increased its completeness without reducing the number of assembled transcripts. Our results suggest that current strategies that rely on available genomic resources as guidance to minimize allelic redundancy are less effective than the application of strategies that cluster redundant assembled transcripts. The strategy yielding the lowest levels of allelic redundancy among the assembled transcriptomes assessed here was the generation of SuperTranscripts with Lace followed by CD-HIT clustering. However, we still observed some levels of heterozygosity (multiple gene fragments per transcript reflecting allelic redundancy) in this assembled transcriptome on the haploid tissue, indicating that further filtering is required before using these assemblies for downstream applications. We discuss the influence of allelic redundancy when these reference transcriptomes are used to select regions for probe design of exome capture baits and for estimation of population genetic diversity.

Comparing observed versus theoretically expected evolutionary responses is important for our understanding of the evolutionary process, and for assessing how species may cope with anthropogenic change. Here, we document directional selection for larger female size in Atlantic salmon, using pedigree-derived estimates of lifetime reproductive success as a fitness measure. We show the trait is heritable and, thus, capable of responding to selection. The Breeder's Equation, which predicts microevolution as the product of phenotypic selection and heritability, predicted evolution of larger size. This was at odds, however, with the observed lack of either phenotypic or genetic temporal trends in body size, a so-called "paradox of stasis." To investigate this paradox, we estimated the additive genetic covariance between trait and fitness, which provides a prediction of evolutionary change according to Robertson's secondary theorem of selection (STS) that is unbiased by missing variables. The STS prediction was consistent with the observed stasis. Decomposition of phenotypic selection gradients into genetic and environmental components revealed a potential upward bias, implying unmeasured factors that covary with trait and fitness. These results showcase the power of pedigreed, wild population studies-which have largely been limited to birds and mammals-to study evolutionary processes on contemporary timescales.

Sexual maturation timing is a life-history trait central to the balance between mortality and reproduction. Maturation may be triggered when an underlying compound trait, called liability, exceeds a threshold. In many different species and especially fishes, this liability is approximated by growth and body condition. However, environmental vs. genetic contributions either directly or via growth and body condition to maturation timing remain unclear. Uncertainty exists also because the maturation process can reverse this causality and itself affect growth and body condition. In addition, disentangling the contributions of polygenic and major loci can be important. In many fishes, males mature before females, enabling the study of associations between male maturation and maturation-unbiased female liability traits. Using 40 Atlantic salmon families, longitudinal common-garden experimentation, and quantitative genetic analyses, we disentangled environmental from polygenic and major locus (vgll3) effects on male maturation, and sex-specific growth and condition. We detected polygenic heritabilities for maturation, growth, and body condition, and vgll3 effects on maturation and body condition but not on growth. Longitudinal patterns for sex-specific phenotypic liability, and for genetic variances and correlations between sexes suggested that early growth and condition indeed positively affected maturation initiation. However, towards spawning time, causality appeared reversed for males whereby maturation affected growth negatively and condition positively via both the environmental and genetic effects. Altogether, the results indicate that growth and condition are useful traits to study liability for maturation initiation, but only until maturation alters their expression, and that vgll3 contributes to maturation initiation via condition.

Signaling centers, or organizers, regulate many aspects of embryonic morphogenesis. In the mammalian molar tooth, reiterative signaling in specialized centers called enamel knots (EKs) determines tooth patterning. Preceding the primary EK, transient epithelial thickening appears, the significance of which remains debated. Using tissue confocal fluorescence imaging with laser ablation experiments, we show that this transient thickening is an earlier signaling center, the molar initiation knot (IK), that is required for the progression of tooth development. IK cell dynamics demonstrate the hallmarks of a signaling center: cell cycle exit, condensation and eventual silencing through apoptosis. IK initiation and maturation are defined by the juxtaposition of cells with high Wnt activity to Shh-expressing non-proliferating cells, the combination of which drives the growth of the tooth bud, leading to the formation of the primary EK as an independent cell cluster. Overall, the whole development of the tooth, from initiation to patterning, is driven by the iterative use of signaling centers.

Whole-genome duplications (WGD) have been considered as springboards that potentiate lineage diversification through increasing functional redundancy. Divergence in gene regulatory elements is a central mechanism for evolutionary diversification, yet the patterns and processes governing regulatory divergence following events that lead to massive functional redundancy, such as WGD, remain largely unknown. We studied the patterns of divergence and strength of natural selection on regulatory elements in the Atlantic salmon (Salmo salar) genome, which has undergone WGD 100–80 Ma. Using ChIPmentation, we first show that H3K27ac, a histone modification typical to enhancers and promoters, is associated with genic regions, tissue-specific transcription factor binding motifs, and with gene transcription levels in immature testes. Divergence in transcription between duplicated genes from WGD (ohnologs) correlated with difference in the number of proximal regulatory elements, but not with promoter elements, suggesting that functional divergence between ohnologs after WGD is mainly driven by enhancers. By comparing H3K27ac regions between duplicated genome blocks, we further show that a longer polyploid state post-WGD has constrained regulatory divergence. Patterns of genetic diversity across natural populations inferred from resequencing indicate that recent evolutionary pressures on H3K27ac regions are dominated by largely neutral evolution. In sum, our results suggest that post-WGD functional redundancy in regulatory elements continues to have an impact on the evolution of the salmon genome, promoting largely neutral evolution of regulatory elements despite their association with transcription levels. These results highlight a case where genome-wide regulatory evolution following an ancient WGD is dominated by genetic drift.

Climate change influences population demography by altering patterns of gene flow and reproductive isolation. Direct mutation rates offer the possibility for accurate dating on the within-species level but are currently only available for a handful of vertebrate species. Here, we use the first directly estimated mutation rate in birds to study the evolutionary history of pied flycatchers ( Ficedula hypoleuca ). Using a combination of demographic inference and species distribution modelling, we show that all major population splits in this forest-dependent system occurred during periods of increased climate instability and rapid global temperature change. We show that the divergent Spanish subspecies originated during the Eemian–Weichselian transition 115–104 thousand years ago (kya), and not during the last glacial maximum (26.5–19 kya), as previously suggested. The magnitude and rates of climate change during the glacial–interglacial transitions that preceded population splits in pied flycatchers were similar to, or exceeded, those predicted to occur in the course of the current, human-induced climate crisis. As such, our results provide a timely reminder of the strong impact that episodes of climate instability and rapid temperature changes can have on species' evolutionary trajectories, with important implications for the natural world in the Anthropocene.

Variation in environmental conditions during development can lead to changes in life-history traits with long-lasting effects. Here, we study how variation in temperature and host plant (i.e. the consequences of potential maternal oviposition choices) affects a suite of life-history traits in pre-diapause larvae of the Glanville fritillary butterfly. We focus on offspring survival, larval growth rates and relative fat reserves, and pay specific attention to intraspecific variation in the responses (G × E × E). Globally, thermal performance and survival curves varied between diets of two host plants, suggesting that host modifies the temperature impact, or vice versa. Additionally, we show that the relative fat content has a host-dependent, discontinuous response to developmental temperature. This implies that a potential switch in resource allocation, from more investment in growth at lower temperatures to storage at higher temperatures, is dependent on the larval diet. Interestingly, a large proportion of the variance in larval performance is explained by differences among families, or interactions with this variable. Finally, we demonstrate that these family-specific responses to the host plant remain largely consistent across thermal environments. Together, the results of our study underscore the importance of paying attention to intraspecific trait variation in the field of evolutionary ecology.

Understanding the genetic basis of repeated evolution of the same phenotype across taxa is a fundamental aim in evolutionary biology and has applications in conservation and management. However, the extent to which interspecific life-history trait polymorphisms share evolutionary pathways remains underexplored. Here, we address this gap by studying the genetic basis of a key life-history trait, age at maturity, in four species of Pacific salmonids (genus Oncorhynchus) that exhibit intra- and interspecific variation in this trait-Chinook Salmon, Coho Salmon, Sockeye Salmon, and Steelhead Trout. We tested for associations in all four species between age at maturity and two genome regions, six6 and vgll3, that are strongly associated with the same trait in Atlantic Salmon (Salmo salar). We also conducted a genome-wide association analysis in Steelhead to assess whether additional regions were associated with this trait. We found the genetic basis of age at maturity to be heterogeneous across salmonid species. Significant associations between six6 and age at maturity were observed in two of the four species, Sockeye and Steelhead, with the association in Steelhead being particularly strong in both sexes (p = 4.46 x 10(-9) after adjusting for genomic inflation). However, no significant associations were detected between age at maturity and the vgll3 genome region in any of the species, despite its strong association with the same trait in Atlantic Salmon. We discuss possible explanations for the heterogeneous nature of the genetic architecture of this key life-history trait, as well as the implications of our findings for conservation and management.

The release of captive-bred animals into the wild is commonly practised to restore or supplement wild populations but comes with a suite of ecological and genetic consequences. Vast numbers of hatchery-reared fish are released annually, ostensibly to restore/enhance wild populations or provide greater angling returns. While previous studies have shown that captive-bred fish perform poorly in the wild relative to wild-bred conspecifics, few have measured individual lifetime reproductive success (LRS) and how this affects population productivity. Here, we analyse data on Atlantic salmon from an intensely studied catchment into which varying numbers of captive-bred fish have escaped/been released and potentially bred over several decades. Using a molecular pedigree, we demonstrate that, on average, the LRS of captive-bred individuals was only 36% that of wild-bred individuals. A significant LRS difference remained after excluding individuals that left no surviving offspring, some of which might have simply failed to spawn, consistent with transgenerational effects on offspring survival. The annual productivity of the mixed population (wild-bred plus captive-bred) was lower in years where captive-bred fish comprised a greater fraction of potential spawners. These results bolster previous empirical and theoretical findings that intentional stocking, or non-intentional escapees, threaten, rather than enhance, recipient natural populations.

A total of 120 Caspian Sea trout specimens from five streams and Lar Lake of the Lar National Park plus 27 Caspian Sea trout specimens (out-group) from the Babolrud River drainage in the Caspian Sea Basin (Iran) were analyzed using 12 microsatellite loci to assess population genetic structure and to estimate the contribution of each population to the lake-run stock. In addition to the Babolrud River population that was significantly differentiated (P

The Atlantic salmon has been studied extensively, particularly as a model for understanding the genetic and environmental contributions to the evolution and development of life history traits. Expression pattern analysis in situ, however, is mostly lacking in salmon. We examine the embryonic developmental expression of six6, a candidate gene previously identified to be associated with spawning ecotypes and age at sexual maturity, in Atlantic salmon. Six6 is a member of the sine oculis homeobox family of transcription factors and is known to regulate eye and brain development in other vertebrates. We assay the expression of this gene in embryonic Atlantic salmon Salmo salar by whole-mount in situ hybridization. In line with earlier studies in other vertebrate species, we find conserved expression in the developing brain and sensory organs, including optic and olfactory primordia. However, we also find previously unreported domains of expression that suggest additional roles in axial and appendicular development, cardiovascular, intestinal, and sensory organogenesis. Each of these systems are important in the sensory ecology of Atlantic salmon, suggesting it is plausible that six6 may have pleiotropic roles in this complex phenotype.

Structural variants (SVs) are a major source of genetic and phenotypic variation, but remain challenging to accurately type and are hence poorly characterized in most species. We present an approach for reliable SV discovery in non-model species using whole genome sequencing and report 15,483 high-confidence SVs in 492 Atlantic salmon (Salmo salar L.) sampled from a broad phylogeographic distribution. These SVs recover population genetic structure with high resolution, include an active DNA transposon, widely affect functional features, and overlap more duplicated genes retained from an ancestral salmonid autotetraploidization event than expected. Changes in SV allele frequency between wild and farmed fish indicate polygenic selection on behavioural traits during domestication, targeting brain-expressed synaptic networks linked to neurological disorders in humans. This study offers novel insights into the role of SVs in genome evolution and the genetic architecture of domestication traits, along with resources supporting reliable SV discovery in non-model species., This study presents and validates a novel approach to reliably identify structural variations (SVs) in non-model genomes using whole genome sequencing, which was used to detect 15,483 SVs in 492 Atlantic salmon, shedding light on their roles in genome evolution and the genetic architecture of domestication.

Global warming is predicted to increase both acute and prolonged thermal challenges for aquatic ectotherms. Severe short- and medium-term thermal stress over hours to days may cause mortality, while longer sub-lethal thermal challenges may cause performance declines. The inter-relationship between the responses to short, medium and longer thermal challenges is unresolved. We asked if the same individuals are tolerant to both rapid and slow warming challenges, a question that has so far received little attention. Additionally, we investigated the possibility of a thermal syndrome where individuals in a population are distributed along a warm-type to cold-type axis. We tested whether different thermal traits correlate across individuals by acclimating 200 juvenile zebrafish (Danio rerio) to sub- or supra-optimal temperatures for growth (22 and 34 degrees C) for 40 days and measuring growth and thermal tolerance at two different warming rates. We found that tolerance to rapid warming correlated with tolerance to slow warming in the 22 degrees C treatment. However, individual tolerance to neither rapid nor slow warming correlated with growth at the supra-optimal temperature. We thus find some support for a syndrome-like organisation of thermal traits, but the lack of connection between tolerance and growth performance indicates a restricted generality of a thermal syndrome. The results suggest that tolerance to rapid warming may share underlying physiological mechanisms with tolerance to slower heating, and indicate that the relevance of acute critical thermal tolerance extends beyond the rapid ramping rates used to measure them.

1. Temperature has a dramatic effect on the physiology of ectothermic animals, impacting most of their biology. When temperatures increase above optimal for an animal, their growth gradually decreases. The main mechanism behind this growth rate reduction is unknown. 2. Here, we suggest the 'aerobic scope protection' hypothesis as a mechanistic explanation for the reduction in growth. 3. After a meal, metabolic rate, and hence oxygen consumption rate, transiently increase in a process called specific dynamic action (SDA). At warmer temperatures, the SDA response usually becomes temporally compressed, leading to a higher peak oxygen consumption rate. This peak in oxygen consumption rate risks taking up much of the animal's aerobic scope (the difference between resting and maximum rates of oxygen consumption), which would leave little residual aerobic scope for other aerobic functions. 4. We propose that water-breathing ectothermic animals will protect their postprandial residual aerobic scope by reducing meal sizes in order to regulate the peak SDA response during times of warming, leading to reductions in growth. 5. This hypothesis is consistent with the published literature on fishes, and we provide predictions that can be tested.

Most Atlantic salmon (Salmo salar L.) populations follow an anadromous life cycle, spending early life in freshwater, migrating to the sea for feeding, and returning to rivers to spawn. At the end of the last ice age ~10,000 years ago, several populations of Atlantic salmon became landlocked. Comparing their genomes to their anadromous counterparts can help identify genetic variation related to either freshwater residency or anadromy. The objective of this study was to identify consistently divergent loci between anadromous and landlocked Atlantic salmon strains throughout their geographical distribution, with the long‐term aim of identifying traits relevant for salmon aquaculture, including fresh and seawater growth, omega‐3 metabolism, smoltification, and disease resistance. We used a Pool‐seq approach (n = 10–40 individuals per population) to sequence the genomes of twelve anadromous and six landlocked Atlantic salmon populations covering a large part of the Northern Hemisphere and conducted a genomewide association study to identify genomic regions having been under different selection pressure in landlocked and anadromous strains. A total of 28 genomic regions were identified and included cadm1 on Chr 13 and ppargc1a on Chr 18. Seven of the regions additionally displayed consistently reduced heterozygosity in fish obtained from landlocked populations, including the genes gpr132, cdca4, and sertad2 on Chr 15. We also found 16 regions, including igf1 on Chr 17, which consistently display reduced heterozygosity in the anadromous populations compared to the freshwater populations, indicating relaxed selection on traits associated with anadromy in landlocked salmon. In conclusion, we have identified 37 regions which may harbor genetic variation relevant for improving fish welfare and quality in the salmon farming industry and for understanding life‐history traits in fish. publishedVersion

Background Understanding genetic architecture is essential for determining how traits will change in response to evolutionary processes such as selection, genetic drift and/or gene flow. In Atlantic salmon, age at maturity is an important life history trait that affects factors such as survival, reproductive success, and growth. Furthermore, age at maturity can seriously impact aquaculture production. Therefore, characterizing the genetic architecture that underlies variation in age at maturity is of key interest. Results Here, we refine our understanding of the genetic architecture for age at maturity of male Atlantic salmon using a genome-wide association study of 11,166 males from a single aquaculture strain, using imputed genotypes at 512,397 single nucleotide polymorphisms (SNPs). All individuals were genotyped with a 50K SNP array and imputed to higher density using parents genotyped with a 930K SNP array and pedigree information. We found significant association signals on 28 of 29 chromosomes (P-values: 8.7 × 10−133–9.8 × 10−8), including two very strong signals spanning the six6 and vgll3 gene regions on chromosomes 9 and 25, respectively. Furthermore, we identified 116 independent signals that tagged 120 candidate genes with varying effect sizes. Five of the candidate genes found here were previously associated with age at maturity in other vertebrates, including humans. Discussion These results reveal a mixed architecture of large-effect loci and a polygenic component that consists of multiple smaller-effect loci, suggesting a more complex genetic architecture of Atlantic salmon age at maturity than previously thought. This more complex architecture will have implications for selection on this key trait in aquaculture and for management of wild salmon populations.

Lake-dwelling fish that form species pairs/flocks characterized by body size divergence are important model systems for speciation research. While several sources of divergent selection have been identified in these systems, their importance for driving the speciation process remains elusive. A major problem is that in retrospect, we cannot distinguish selection pressures that initiated divergence from those acting later in the process. To address this issue, we reconstructed the initial stages of speciation in European whitefish (Coregonus lavaretus) using data from 357 populations of varying age (26-10 000 years). We find that whitefish speciation is driven by a large-growing predator, the northern pike (Esox lucius). Pike initiates divergence by causing a largely plastic differentiation into benthic giants and pelagic dwarfs; ecotypes that will subsequently develop partial reproductive isolation and heritable differences in gill raker number. Using an eco-evolutionary model, we demonstrate how pike’s habitat specificity and large gape size are critical for imposing a between-habitat trade-off, causing prey to mature in a safer place or at a safer size. Thereby, we propose a novel mechanism for how predators may cause dwarf/giant speciation in lake-dwelling fish species.

In species with complex life cycles, life history theory predicts that fitness is affected by conditions encountered in previous life history stages. Here, we use a four-year pedigree to investigate if time spent in two distinct life history stages has sex-specific reproductive fitness consequences in anadromous Atlantic salmon (Salmo salar). We determined the amount of years spent in fresh water as juveniles (freshwater age, FW, measured in years), and years spent in the marine environment as adults (sea age, SW, measured in sea winters) on 264 sexually mature adults collected on a river spawning ground. We then estimated reproductive fitness as the number of offspring (reproductive success) and the number of mates (mating success) using genetic parentage analysis (>5000 offspring). Sea age is significantly and positively correlated with reproductive and mating success of both sexes whereby older and larger individuals gained the highest reproductive fitness benefits (females: 62.2% increase in offspring/SW and 34.8% increase in mate number/SW; males: 201.9% offspring/SW and 60.3% mates/SW). Younger freshwater age was significantly related to older sea age and thus increased reproductive fitness, but only among females (females: -33.9% offspring/FW and -32.4% mates/FW). This result implies that females can obtain higher reproductive fitness by transitioning to the marine environment earlier. In contrast, male mating and reproductive success was unaffected by freshwater age and more males returned at a younger age than females despite the reproductive fitness advantage of later sea age maturation. Our results show that the timing of transitions between juvenile and adult phases has a sex-specific consequence on female reproductive fitness, demonstrating a life-history trade-off between maturation and reproduction in wild Atlantic salmon.

Adaptation is thought to proceed in part through spatial and temporal changes in gene expression. Fish species such as the threespine stickleback are powerful vertebrate models to study the genetic architecture of adaptive changes in gene expression since divergent adaptation to different environments is common, they are abundant and easy to study in the wild and lab, and have well-established genetic and genomic resources. Fish gills, due to their respiratory and osmoregulatory roles, show many physiological adaptations to local water chemistry, including differences in gene expression. However, obtaining high-quality RNA using popular column-based extraction methods can be challenging from small tissue samples high in cartilage and bone such as fish gills. Here, we describe a bead-based mRNA extraction and transcriptome RNA-seq protocol that does not use purification columns. The protocol can be readily scaled according to sample size for the purposes of diverse gene expression experiments using animal or plant tissue.

Gene expression changes have been recognized as important drivers of adaptation to changing environmental conditions. Little is known about the relative roles of plastic and evolutionary responses in complex gene expression networks during the early stages of divergence. Large gene expression data sets coupled with in silico methods for identifying co-expressed modules now enable systems genetics approaches also in non-model species for better understanding of gene expression responses during early divergence. Here, we combined gene co-expression analyses with population genetics to separate plastic and population (evolutionary) effects in expression networks using small salmonid populations as a model system. We show that plastic and population effects were highly variable among the six identified modules and that the plastic effects explained larger proportion of the total eigengene expression than population effects. A more detailed analysis of the population effects using a QST - FST comparison across 16622 annotated transcripts revealed that gene expression followed neutral expectations within modules and at the global level. Furthermore, two modules showed enrichment for genes coding for early developmental traits that have been previously identified as important phenotypic traits in thermal responses in the same model system indicating that co-expression analysis can capture expression patterns underlying ecologically important traits. We suggest that module-specific responses may facilitate the flexible tuning of expression levels to local thermal conditions. Overall, our study indicates that plasticity and neutral evolution are the main drivers of gene expression variance in the early stages of thermal adaptation in this system.

Understanding the genomic basis of host-parasite adaptation is important for predicting the long-term viability of species and developing successful management practices. However, in wild populations, identifying specific signatures of parasite-driven selection often presents a challenge, as it is difficult to unravel the molecular signatures of selection driven by different, but correlated, environmental factors. Furthermore, separating parasite-mediated selection from similar signatures due to genetic drift and population history can also be difficult. Populations of Atlantic salmon (Salmo salar L.) from northern Europe have pronounced differences in their reactions to the parasitic flatworm Gyrodactylus salaris Malmberg 1957 and are therefore a good model to search for specific genomic regions underlying inter-population differences in pathogen response. We used a dense Atlantic salmon SNP array, along with extensive sampling of 43 salmon populations representing the two G. salaris response extremes (extreme susceptibility vs resistant), to screen the salmon genome for signatures of directional selection while attempting to separate the parasite effect from other factors. After combining the results from two independent genome scan analyses, 57 candidate genes potentially under positive selection were identified, out of which 50 were functionally annotated. This candidate gene set was shown to be functionally enriched for lymph node development, focal adhesion genes and anti-viral response, which suggests that the regulation of both innate and acquired immunity might be an important mechanism for salmon response to G. salaris. Overall, our results offer insights into the apparently complex genetic basis of pathogen susceptibility in salmon and highlight methodological challenges for separating the effects of various environmental factors.

A major goal in biology is to understand how evolution shapes variation in individual life histories. Genome-wide association studies have been successful in uncovering genome regions linked with traits underlying life history variation in a range of species. However, lack of functional studies of the discovered genotype-phenotype associations severely restrains our understanding how alternative life history traits evolved and are mediated at the molecular level. Here, we report a cis-regulatory mechanism whereby expression of alternative isoforms of the transcription co-factor vestigial-like 3 (vgll3) associate with variation in a key life history trait, age at maturity, in Atlantic salmon (Salmo salar). Using a common-garden experiment, we first show that vgll3 genotype associates with puberty timing in one-year-old salmon males. By way of temporal sampling of vgll3 expression in ten tissues across the first year of salmon development, we identify a pubertal transition in vgll3 expression where maturation coincided with a 66% reduction in testicular vgll3 expression. The late maturation allele was not only associated with a tendency to delay puberty, but also with expression of a rare transcript isoform of vgll3 pre-puberty. By comparing absolute vgll3 mRNA copies in heterozygotes we show that the expression difference between the early and late maturity alleles is largely cis-regulatory. We propose a model whereby expression of a rare isoform from the late allele shifts the liability of its carriers towards delaying puberty. These results exemplify the potential importance of regulatory differences as a mechanism for the evolution of life history traits., Author summary Alternative life history strategies are an important source of diversity within populations and promote the maintenance of adaptive capacity and population resilience. However, in many cases the molecular basis of different life history strategies remains elusive. Age at maturity is a key adaptive life history trait in Atlantic salmon and has a relatively simple genetic basis. Using salmon age at maturity as a model, we report a mechanism whereby different transcript isoforms of the key age at maturity gene, vestigial-like 3 (vgll3), associate with variation in the timing of male puberty. Our results show how gene regulatory differences in conjunction with variation in gene transcript structure can encode for complex alternative life histories.

Recent progress has been made in identifying genomic regions implicated in trait evolution on a microevolutionary scale in many species, but whether these are relevant over macroevolutionary time remains unclear. Here, we directly address this fundamental question using bird beak shape, a key evolutionary innovation linked to patterns of resource use, divergence and speciation, as a model trait. We design a novel approach by integrating class-wide geometric-morphometric analyses with evolutionary sequence analyses of 10,322 protein coding genes as well as 229,001 genomic regions spanning 72 species. We identify 1,434 protein coding genes and 39,806 non-coding regions for which molecular rates were significantly related to rates of bill shape evolution. We show that homologs of the identified protein coding genes as well as genes in close proximity to the identified non-coding regions are involved in craniofacial embryo development in mammals. We identify two associated pathways, BMP and Wnt signalling, both of which have repeatedly been implicated in the morphological development of avian beaks. This result suggests that our approach identifying genotype-phenotype association on a genome wide scale over macroevolutionary time is feasible. While the coding and non-coding gene sets are associated with similar pathways, the actual genes are highly distinct, with significantly reduced overlap between them and bill-related phenotype associations specific to non-coding loci. Moreover we find evidence for signatures of recent diversifying selection in the identified non-coding loci in populations of Darwin9s finches that differ notably in their beak shape morphology. Our results suggest that regulatory rather than coding changes are major drivers of morphological diversification over macroevolutionary times.

Despite recent taxonomic diversification in studies linking genotype with phenotype, follow-up studies aimed at understanding the molecular processes of such genotype-phenotype associations remain rare. The age at which an individual reaches sexual maturity is an important fitness trait in many wild species. However, the molecular mechanisms regulating maturation timing processes remain obscure. A recent genome-wide association study in Atlantic salmon (Salmo salar) identified large-effect age-at-maturity-associated chromosomal regions including genesvgll3,akap11andsix6, which have roles in adipogenesis, spermatogenesis and the hypothalamic-pituitary-gonadal (HPG) axis, respectively. Here, we determine expression patterns of these genes during salmon development and their potential molecular partners and pathways. Using Nanostring transcription profiling technology, we show development- and tissue-specific mRNA expression patterns forvgll3,akap11andsix6. Correlated expression levels ofvgll3andakap11, which have adjacent chromosomal location, suggests they may have shared regulation. Further,vgll3correlating witharhgap6andyap1, andakap11withlats1andyap1suggests that Vgll3 and Akap11 take part in actin cytoskeleton regulation. Tissue-specific expression results indicate thatvgll3andakap11paralogs have sex-dependent expression patterns in gonads. Moreover,six6correlating withslc38a6andrtn1, and Hippo signaling genes suggests that Six6 could have a broader role in the HPG neuroendrocrine and cell fate commitment regulation, respectively. We conclude that Vgll3, Akap11 and Six6 may influence Atlantic salmon maturation timing via affecting on adipogenesis and gametogenesis by regulating cell fate commitment and the HPG axis. These results may help to unravel general molecular mechanisms behind maturation.

Co-inheritance in life-history traits may result in unpredictable evolutionary trajectories if not accounted for in life-history models. Iteroparity (the reproductive strategy of reproducing more than once) in Atlantic salmon (Salmo salar) is a fitness trait with substantial variation within and among populations. In the Teno River in northern Europe, iteroparous individuals constitute an important component of many populations and have experienced a sharp increase in abundance in the last 20 years, partly overlapping with a general decrease in age structure. The physiological basis of iteroparity bears similarities to that of age at first maturity, another life-history trait with substantial fitness effects in salmon. Sea age at maturity in Atlantic salmon is controlled by a major locus around the vgll3 gene, and we used this opportunity demonstrate that these two traits are co-inherited around this genome region. The odds ratio of survival until second reproduction was up to 2.4 (1.8-3.5 90% CI) times higher for fish with the early-maturing vgll3 genotype (EE) compared to fish with the late-maturing genotype (LL). The L allele was dominant in individuals remaining only one year at sea before maturation, but the dominance was reversed, with the E allele being dominant in individuals maturing after two or more years at sea. Post hoc analysis indicated that iteroparous fish with the EE genotype had accelerated growth prior to first reproduction compared to first-time spawners, across all age groups, whereas this effect was not detected in fish with the LL genotype. These results broaden the functional link around the vgll3 genome region and help us understand constraints in the evolution of life-history variation in salmon. Our results further highlight the need to account for genetic correlations between fitness traits when predicting demographic changes in changing environments.