Comparative transcriptomics reveals domestication-associated features of Atlantic salmon lipid metabolism.

Domestication of animals imposes strong targeted selection for desired traits but can also result in unintended selection due to new domestic environments. Atlantic salmon (Salmo salmar) was domesticated in the 1970s and has subsequently been selected for faster growth in systematic breeding programmes. More recently, salmon aquaculture has replaced fish oils (FOs) with vegetable oils (VOs) in feed, radically changing the levels of essential long-chain polyunsaturated fatty acids (LC-PUFAs). Our aim here was to study the impact of domestication on metabolism and explore the hypothesis that the shift to VO diets has unintentionally selected for a domestication-specific lipid metabolism. We conducted a 96-day feeding trial of domesticated and wild salmon fed diets based on FOs, VOs or phospholipids, and compared transcriptomes and fatty acids in tissues involved in lipid absorption (pyloric caeca) and lipid turnover and synthesis (liver). Domesticated salmon had faster growth and higher gene expression in glucose and lipid metabolism compared to wild fish, possibly linked to differences in regulation of circadian rhythm pathways. Only the domesticated salmon increased expression of LC-PUFA synthesis genes when given VOs. This transcriptome response difference was mirrored at the physiological level, with domesticated salmon having higher LC-PUFA levels but lower 18:3n-3 and 18:2n-6 levels. In line with this, the VO diet decreased growth rate in wild but not domesticated salmon. Our study revealed a clear impact of domestication on transcriptomic regulation linked to metabolism and suggests that unintentional selection in the domestic environment has resulted in evolution of stronger compensatory mechanisms to a diet low in LC-PUFAs.

Genomic regions and signaling pathways associated with indicator traits for feed efficiency in juvenile Atlantic salmon (Salmo salar).

BACKGROUND: One objective of this study was to identify putative quantitative trait loci (QTL) that affect indicator phenotypes for growth, nitrogen, and carbon metabolism in muscle, liver, and adipose tissue, and for feed efficiency. Another objective was to perform an RNAseq analysis (184 fish from all families), to identify genes that are associated with carbon and nitrogen metabolism in the liver. The material consisted of a family experiment that was performed in freshwater and included 2281 individuals from 23 full-sib families. During the 12-day feed conversion test, families were randomly allocated to family tanks (50 fish per tank and 2 tanks per family) and fed a fishmeal-based diet labeled with the stable isotopes 15N and 13C at inclusion levels of 2 and 1%, respectively. RESULTS: Using a linear mixed-model algorithm, a QTL for pre-smolt growth was identified on chromosome 9 and a QTL for carbon metabolism in the liver was identified on chromosome 12 that was closely related to feed conversion ratio on a tank level. For the indicators of feed efficiency traits that were derived from the stable isotope ratios (15N and 13C) of muscle tissue and growth, no convincing QTL was detected, which suggests that these traits are polygenic. The transcriptomic analysis showed that high carbon and nitrogen metabolism was associated with individuals that convert protein from the feed more efficiently, primarily due to higher expression of the proteasome, lipid, and carbon metabolic pathways in liver. In addition, we identified seven transcription factors that were associated with carbon and nitrogen metabolism and located in the identified QTL regions. CONCLUSIONS: Analyses revealed one QTL associated with pre-smolt growth and one QTL for carbon metabolism in the liver. Both of these traits are associated with feed efficiency. However, more accurate mapping of the putative QTL will require a more diverse family material. In this experiment, fish that have a high carbon and nitrogen metabolism in the liver converted protein from the feed more efficiently, potentially because of a higher expression of the proteasome, lipid, and carbon metabolic pathways in liver. Within the QTL regions, we detected seven transcription factors that were associated with carbon and nitrogen metabolism.

SalMotifDB: a tool for analyzing putative transcription factor binding sites in salmonid genomes.

BACKGROUND: Recently developed genome resources in Salmonid fish provides tools for studying the genomics underlying a wide range of properties including life history trait variation in the wild, economically important traits in aquaculture and the evolutionary consequences of whole genome duplications. Although genome assemblies now exist for a number of salmonid species, the lack of regulatory annotations are holding back our mechanistic understanding of how genetic variation in non-coding regulatory regions affect gene expression and the downstream phenotypic effects. RESULTS: We present SalMotifDB, a database and associated web and R interface for the analysis of transcription factors (TFs) and their cis-regulatory binding sites in five salmonid genomes. SalMotifDB integrates TF-binding site information for 3072 non-redundant DNA patterns (motifs) assembled from a large number of metazoan motif databases. Through motif matching and TF prediction, we have used these multi-species databases to construct putative regulatory networks in salmonid species. The utility of SalMotifDB is demonstrated by showing that key lipid metabolism regulators are predicted to regulate a set of genes affected by different lipid and fatty acid content in the feed, and by showing that our motif database explains a significant proportion of gene expression divergence in gene duplicates originating from the salmonid specific whole genome duplication. CONCLUSIONS: SalMotifDB is an effective tool for analyzing transcription factors, their binding sites and the resulting gene regulatory networks in salmonid species, and will be an important tool for gaining a better mechanistic understanding of gene regulation and the associated phenotypes in salmonids. SalMotifDB is available at https://salmobase.org/apps/SalMotifDB .

Transcriptional development of phospholipid and lipoprotein metabolism in different intestinal regions of Atlantic salmon (Salmo salar) fry.

BACKGROUND: It has been suggested that the high phospholipid (PL) requirement in Atlantic salmon (Salmo salar) fry is due to insufficient intestinal de-novo synthesis causing low lipoprotein (LP) production and reduced transport capacity of dietary lipids. However, in-depth ontogenetic analysis of intestinal PL and LP synthesis with the development of salmon has yet to be performed. Therefore, in this paper we used RNA-Seq technology to investigate the expression of genes involved in PL synthesis and LP formation throughout early developmental stages and associate insufficient expression of synthesis pathways in salmon fry with its higher dietary PL requirement. There was a special focus on the understanding homologous genes, especially those from salmonid-specific fourth vertebrate whole-genome duplication (Ss4R), and their contribution to salmonid specific features of regulation of PL metabolic pathways. Salmon fry were sampled at 0.16 g (1 day before first-feeding), 2.5 and 10 g stages of development and transcriptomic analysis was applied separately on stomach, pyloric caeca and hindgut of the fish. RESULTS: In general, we found up-regulated pathways involved in synthesis of phosphatidylcholine (PtdCho), phosphatidylethanolamine (PtdEtn), and LP in pyloric caeca of salmon between 0.16 and 10 g. Thirteen differentially expressed genes (q < 0.05) in these pathways were highly up-regulated in 2.5 g salmon compared to 0.16 g, while only five more differentially expressed (q < 0.05) genes were found when the fish grew up to 10 g. Different homologous genes were found dominating in stomach, pyloric caeca and hindgut. However, the expression of dominating genes in pathways of PL and LP synthesis were much higher in pyloric caeca than stomach and hindgut. Salmon-specific homologous genes (Ss4R) had similar expression during development, while other homologs had more diverged expression. CONCLUSIONS: The up-regulation of the de-novo PtdCho and PtdEtn pathways confirm that salmon have decreasing requirement for dietary PL as the fish develops. The similar expressions between Ss4R homologous genes suggest that the functional divergence of these genes was incomplete compared to homologs derived from other genome duplication. The results of the present study have provided new information on the molecular mechanisms of phospholipid synthesis and lipoprotein formation in fish.

Stable Core Gut Microbiota across the Freshwater-to-Saltwater Transition for Farmed Atlantic Salmon.

Gut microbiota associations through habitat transitions are fundamentally important yet poorly understood. One such habitat transition is the migration from freshwater to saltwater for anadromous fish, such as salmon. The aim of the current work was therefore to determine the freshwater-to-saltwater transition impact on the gut microbiota in farmed Atlantic salmon, with dietary interventions resembling freshwater and saltwater diets with respect to fatty acid composition. Using deep 16S rRNA gene sequencing and quantitative PCR, we found that the freshwater-to-saltwater transition had a major association with the microbiota composition and quantity, while diet did not show significant associations with the microbiota. In saltwater there was a 100-fold increase in bacterial quantity, with a relative increase of Firmicutes and a relative decrease of both Actinobacteria and Proteobacteria Irrespective of an overall shift in microbiota composition from freshwater to saltwater, we identified three core clostridia and one Lactobacillus-affiliated phylotype with wide geographic distribution that were highly prevalent and co-occurring. Taken together, our results support the importance of the dominating bacteria in the salmon gut, with the freshwater microbiota being immature. Due to the low number of potentially host-associated bacterial species in the salmon gut, we believe that farmed salmon can represent an important model for future understanding of host-bacterium interactions in aquatic environments.IMPORTANCE Little is known about factors affecting the interindividual distribution of gut bacteria in aquatic environments. We have shown that there is a core of four highly prevalent and co-occurring bacteria irrespective of feed and freshwater-to-saltwater transition. The potential host interactions of the core bacteria, however, need to be elucidated further.

SalmoBase: an integrated molecular data resource for Salmonid species.

BACKGROUND: Salmonids are ray-finned fishes which constitute 11 genera and at least 70 species including Atlantic salmon, whitefishes, graylings, rainbow trout, and char. The common ancestor of all Salmonidae experienced a whole genome duplication (WGD) ~80 million years ago, resulting in an autotetraploid genome. Genomic rediplodization is still going on in salmonid species, providing an unique system for studying evolutionary consequences of whole genome duplication. In recent years, high quality genome sequences of Atlantic salmon and Rainbow trout has been established, due to their scientific and commercial values. In this paper we introduce SalmoBase ( http://www.salmobase.org/ ), a tool for making molecular resources for salmonids public available in a framework of visualizations and analytic tools. RESULTS: SalmoBase has been developed as a part of the ELIXIR.NO project. Currently, SalmoBase contains molecular resources for Atlantic salmon and Rainbow trout. Data can be accessed through BLAST, Genome Browser (GBrowse), Genetic Variation Browser (GVBrowse) and Gene Expression Browser (GEBrowse). CONCLUSIONS: To the best of our knowledge, SalmoBase is the first database which integrates salmonids data and allow users to study salmonids in an integrated framework. The database and its tools (e.g., comparative genomics tools, synteny browsers) will be expanded as additional public resources describing other Salmonidae genomes become available.

Extracting functional trends from whole genome duplication events using comparative genomics.

BACKGROUND: The number of species with completed genomes, including those with evidence for recent whole genome duplication events has exploded. The recently sequenced Atlantic salmon genome has been through two rounds of whole genome duplication since the divergence of teleost fish from the lineage that led to amniotes. This quadrupoling of the number of potential genes has led to complex patterns of retention and loss among gene families. RESULTS: Methods have been developed to characterize the interplay of duplicate gene retention processes across both whole genome duplication events and additional smaller scale duplication events. Further, gene expression divergence data has become available as well for Atlantic salmon and the closely related, pre-whole genome duplication pike and methods to describe expression divergence are also presented. These methods for the characterization of duplicate gene retention and gene expression divergence that have been applied to salmon are described. CONCLUSIONS: With the growth in available genomic and functional data, the opportunities to extract functional inference from large scale duplicates using comparative methods have expanded dramatically. Recently developed methods that further this inference for duplicated genes have been described.

Gene family expansion and functional diversification of chitinase and chitin synthase genes in Atlantic salmon (Salmo salar).

Chitin is one of the most abundant polysaccharides in nature, forming important structures in insects, crustaceans, and fungal cell walls. Vertebrates on the other hand are generally considered "nonchitinous" organisms, despite having highly conserved chitin metabolism-associated genes. Recent work has revealed that the largest group of vertebrates, the teleosts, have the potential to both synthesize and degrade endogenous chitin. Yet, little is known about the genes and proteins responsible for these dynamic processes. Here, we used comparative genomics, transcriptomics, and chromatin accessibility data to characterize the repertoire, evolution, and regulation of genes involved in chitin metabolism in teleosts, with a particular focus on Atlantic salmon. Reconstruction of gene family phylogenies provides evidence for an expansion of teleost and salmonid chitinase and chitin synthase genes after multiple whole-genome duplications. Analyses of multi-tissue gene expression data demonstrated a strong bias of gastrointestinal tract expression for chitin metabolism genes, but with different spatial and temporal tissue specificities. Finally, we integrated transcriptomes from a developmental time series of the gastrointestinal tract with chromatin accessibility data to identify putative transcription factors responsible for regulating chitin metabolism gene expression (CDX1 and CDX2) as well as tissue-specific divergence in the regulation of gene duplicates (FOXJ2). The findings presented here support the hypothesis that chitin metabolism genes in teleosts play a role in developing and maintaining a chitin-based barrier in the teleost gut and provide a basis for further investigations into the molecular basis of this barrier.

Functional validation of transposable element-derived cis-regulatory elements in Atlantic salmon.

Transposable elements (TEs) are hypothesized to play important roles in shaping genome evolution following whole-genome duplications (WGDs), including rewiring of gene regulation. In a recent analysis, duplicate gene copies that had evolved higher expression in liver following the salmonid WGD ∼100 million years ago were associated with higher numbers of predicted TE-derived cis-regulatory elements (TE-CREs). Yet, the ability of these TE-CREs to recruit transcription factors (TFs) in vivo and impact gene expression remains unknown. Here, we evaluated the gene-regulatory functions of 11 TEs using luciferase promoter reporter assays in Atlantic salmon (Salmo salar) primary liver cells. Canonical Tc1-Mariner elements from intronic regions showed no or small repressive effects on transcription. However, other TE-CREs upstream of transcriptional start sites increased expression significantly. Our results question the hypothesis that TEs in the Tc1-Mariner superfamily, which were extremely active following WGD in salmonids, had a major impact on regulatory rewiring of gene duplicates, but highlights the potential of other TEs in post-WGD rewiring of gene regulation in the Atlantic salmon genome.

Genetic mapping of DArT markers in the Festuca-Lolium complex and their use in freezing tolerance association analysis.

Species belonging to the Festuca-Lolium complex are important forage and turf species and as such, have been studied intensively. However, their out-crossing nature and limited availability of molecular markers make genetic studies difficult. Here, we report on saturation of F. pratensis and L. multiflorum genetic maps using Diversity Array Technology (DArT) markers and the DArTFest array.The 530 and 149 DArT markers were placed on genetic maps of L. multiflorum and F. pratensis, respectively, with overlap of 20 markers, which mapped in both species. The markers were sequenced and comparative sequence analysis was performed between L. multiflorum, rice and Brachypodium. The utility of the DArTFest array was then tested on a Festulolium population FuRs0357 in an integrated analysis using the DArT marker map positions to study associations between markers and freezing tolerance. Ninety six markers were significantly associated with freezing tolerance and five of these markers were genetically mapped to chromosomes 2, 4 and 7. Three genomic loci associated with freezing tolerance in the FuRs0357 population co-localized with chromosome segments and QTLs previously identified to be associated with freezing tolerance. The present work clearly confirms the potential of the DArTFest array in genetic studies of the Festuca-Lolium complex. The annotated DArTFest array resources could accelerate further studies and improvement of desired traits in Festuca-Lolium species.

Diet and Life Stage-Associated Lipidome Remodeling in Atlantic Salmon.

Salmon is an important source of long-chain highly unsaturated fatty acids (LC-HUFAs) such as 22:6n-3 [docosahexaenoic acid (DHA)]. In the present study, we conducted two identical experiments on salmon in freshwater (FW) and seawater (SW) stages, with a diet switch from fish oil (high in LC-HUFA) to vegetable oil (low in LC-HUFA) and vice versa. Our aim was to investigate the diet and life stage-specific features of lipid uptake (gut), processing (liver), and deposition (muscle). The lipid composition changed much faster in the gut of SW fish relative to FW fish, suggesting that the former had a higher rate of lipid absorption and transport. SW fish also had higher expression of phospholipid synthesis and lipoprotein formation genes in the gut, whereas FW fish had higher expression of lipid synthesis genes in the liver. All phospholipids except PC-44:12 and PE-44:12 were less abundant in SW, suggesting that SW fish have a higher requirement for DHA.

Photoperiod-dependent developmental reprogramming of the transcriptional response to seawater entry in Atlantic salmon (Salmo salar).

The developmental transition of juvenile salmon from a freshwater resident morph (parr) to a seawater (SW) migratory morph (smolt), known as smoltification, entails a reorganization of gill function to cope with the altered water environment. Recently, we used RNAseq to characterize the breadth of transcriptional change which takes place in the gill in the FW phase of smoltification. This highlighted the importance of extended exposure to short, winter-like photoperiods (SP) followed by a subsequent increase in photoperiod for completion of transcriptional reprogramming in FW and efficient growth following transfer to SW. Here, we extend this analysis to examine the consequences of this photoperiodic history-dependent reprogramming for subsequent gill responses upon exposure to SW. We use RNAseq to analyze gill samples taken from fish raised on the photoperiod regimes we used previously and then challenged by SW exposure for 24 hours. While fish held on constant light (LL) throughout were able to hypo-osmoregulate during a 24 hours SW challenge, the associated gill transcriptional response was highly distinctive from that in fish which had experienced a 7-week period of exposure to SP followed by a return to LL (SPLL) and had consequently acquired the characteristics of fully developed smolts. Fish transferred from LL to SP, and then held on SP for the remainder of the study was unable to hypo-osmoregulate, and the associated gill transcriptional response to SW exposure featured many transcripts apparently regulated by the glucocorticoid stress axis and by the osmo-sensing transcription factor NFAT5. The importance of these pathways for the gill transcriptional response to SW exposure appears to diminish as a consequence of photoperiod mediated induction of the smolt phenotype, presumably reflecting preparatory developmental changes taking place during this process.

Immunologic Profiling of the Atlantic Salmon Gill by Single Nuclei Transcriptomics.

Anadromous salmonids begin life adapted to the freshwater environments of their natal streams before a developmental transition, known as smoltification, transforms them into marine-adapted fish. In the wild, smoltification is a photoperiod-regulated process, involving radical remodeling of gill function to cope with the profound osmotic and immunological challenges of seawater (SW) migration. While prior work has highlighted the role of specialized "mitochondrion-rich" cells (MRCs) and accessory cells (ACs) in delivering this phenotype, recent RNA profiling experiments suggest that remodeling is far more extensive than previously appreciated. Here, we use single-nuclei RNAseq to characterize the extent of cytological changes in the gill of Atlantic salmon during smoltification and SW transfer. We identify 20 distinct cell clusters, including known, but also novel gill cell types. These data allow us to isolate cluster-specific, smoltification-associated changes in gene expression and to describe how the cellular make-up of the gill changes through smoltification. As expected, we noted an increase in the proportion of seawater mitochondrion-rich cells, however, we also identify previously unknown reduction of several immune-related cell types. Overall, our results provide fresh detail of the cellular complexity in the gill and suggest that smoltification triggers unexpected immune reprogramming.

Did gene family expansions during the Eocene-Oligocene boundary climate cooling play a role in Pooideae adaptation to cool climates?

Adaptation to cool environments is a common feature in the core group of the grass subfamily Pooideae (Triticeae and Poeae). This suggest an ancient evolutionary origin of low temperature stress tolerance dating back prior to the initiation of taxonomic divergence of core Pooideae species. Viewing the Pooideae evolution in a palaeo-climatic perspective reveals that taxonomic divergence of the core Pooideae group initiated shortly after a global super-cooling period at the Eocene-Oligocene boundary (approximately 33.5-26 Ma). This global climate cooling altered distributions of plants and animals and must have imposed selection pressure for improved low temperature stress responses. Lineage-specific gene family expansions are known to be involved in adaptation to new environmental stresses. In Pooideae, two gene families involved in low temperature stress response, the C-repeat binding factor (CBF) and fructosyl transferase (FT) gene families, has undergone lineage-specific expansions. We investigated the timing of these gene family expansions by molecular dating and found that Pooideae-specific expansion events in CBF and FT gene families took place during Eocene-Oligocene super-cooling period. We hypothesize that the E-O super-cooling exerted selection pressure for improved low temperature stress response and frost tolerance in a core Pooideae ancestor, and that those individuals with multiple copies of CBF and FT genes were favoured.

Publisher Correction: Sex-dependent dominance maintains migration supergene in rainbow trout.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Sex-dependent dominance maintains migration supergene in rainbow trout.

Males and females often differ in their fitness optima for shared traits that have a shared genetic basis, leading to sexual conflict. Morphologically differentiated sex chromosomes can resolve this conflict and protect sexually antagonistic variation, but they accumulate deleterious mutations. However, how sexual conflict is resolved in species that lack differentiated sex chromosomes is largely unknown. Here we present a chromosome-anchored genome assembly for rainbow trout (Oncorhynchus mykiss) and characterize a 55-Mb double-inversion supergene that mediates sex-specific migratory tendency through sex-dependent dominance reversal, an alternative mechanism for resolving sexual conflict. The double inversion contains key photosensory, circadian rhythm, adiposity and sex-related genes and displays a latitudinal frequency cline, indicating environmentally dependent selection. Our results show sex-dependent dominance reversal across a large autosomal supergene, a mechanism for sexual conflict resolution capable of protecting sexually antagonistic variation while avoiding the homozygous lethality and deleterious mutations associated with typical heteromorphic sex chromosomes.

The evolution and functional divergence of the beta-carotene oxygenase gene family in teleost fish--exemplified by Atlantic salmon.

In mammals, two carotenoid cleaving oxygenases are known; beta-carotene 15,15'-monooxygenase (BCMO1) and beta-carotene 9',10'-oxygenase (BCO2). BCMO1 is a key enzyme in vitamin A synthesis by symmetrically cleaving beta-carotene into 2 molecules of all-trans-retinal, while BCO2 is responsible for asymmetric cleavage of a broader range of carotenoids. Here, we show that the Atlantic salmon beta-carotene oxygenase (bco) gene family contains 5 members, three bco2 and two bcmo1 paralogs. Using public sequence databases, multiple bco genes were also found in several additional teleost species. Phylogenetic analysis indicates that bco2a and bco2b originate from the teleost fish specific genome duplication (FSGD or 3R), while the third and more distant paralog, bco2 like, might stem from a prior duplication event in the teleost lineage. The two bcmo1 paralogs (bcmo1 and bcmo1 like) appear to be the result of an ancient duplication event that took place before the divergence of ray-finned (Actinopterygii) and lobe-finned fish (Sarcopterygii), with subsequent nonfunctionalization and loss of one Sarcopterygii paralog. Gene expression analysis of the bcmo1 and bco2 paralogs in Atlantic salmon reveals regulatory divergence with tissue specific expression profiles, suggesting that the beta-carotene oxygenase subtypes have evolved functional divergences. We suggest that teleost fish have evolved and maintained an extended repertoire of beta-carotene oxygenases compared to the investigated Sarcopterygii species, and hypothesize that the main driver behind this functional divergence is the exposure to a diverse set of carotenoids in the aquatic environment.