The genome regulatory landscape of Atlantic salmon liver through smoltification.

The anadromous Atlantic salmon undergo a preparatory physiological transformation before seawater entry, referred to as smoltification. Key molecular developmental processes involved in this life stage transition, such as remodeling of gill functions, are known to be synchronized and modulated by environmental cues like photoperiod. However, little is known about the photoperiod influence and genome regulatory processes driving other canonical aspects of smoltification such as the large-scale changes in lipid metabolism and energy homeostasis in the developing smolt liver. Here we generate transcriptome, DNA methylation, and chromatin accessibility data from salmon livers across smoltification under different photoperiod regimes. We find a systematic reduction of expression levels of genes with a metabolic function, such as lipid metabolism, and increased expression of energy related genes such as oxidative phosphorylation, during smolt development in freshwater. However, in contrast to similar studies of the gill, smolt liver gene expression prior to seawater transfer was not impacted by photoperiodic history. Integrated analyses of gene expression, chromatin accessibility, and transcription factor (TF) binding signatures highlight chromatin remodeling and TF dynamics underlying smolt gene regulatory changes. Differential peak accessibility patterns largely matched differential gene expression patterns during smoltification and we infer that ZNF682, KLFs, and NFY TFs are important in driving a liver metabolic shift from synthesis to break down of organic compounds in freshwater. Overall, chromatin accessibility and TFBS occupancy were highly correlated to changes in gene expression. On the other hand, we identified numerous differential methylation patterns across the genome, but associated genes were not functionally enriched or correlated to observed gene expression changes across smolt development. Taken together, this work highlights the relative importance of chromatin remodeling during smoltification and demonstrates that metabolic remodeling occurs as a preadaptation to life at sea that is not to a large extent driven by photoperiod history.

Hypothalamic tanycytes as mediators of maternally programmed seasonal plasticity.

In mammals, maternal photoperiodic programming (MPP) provides a means whereby juvenile development can be matched to forthcoming seasonal environmental conditions.1,2,3,4 This phenomenon is driven by in utero effects of maternal melatonin5,6,7 on the production of thyrotropin (TSH) in the fetal pars tuberalis (PT) and consequent TSH receptor-mediated effects on tanycytes lining the 3rd ventricle of the mediobasal hypothalamus (MBH).8,9,10 Here we use LASER capture microdissection and transcriptomic profiling to show that TSH-dependent MPP controls the attributes of the ependymal region of the MBH in juvenile animals. In Siberian hamster pups gestated and raised on a long photoperiod (LP) and thereby committed to a fast trajectory for growth and reproductive maturation, the ependymal region is enriched for tanycytes bearing sensory cilia and receptors implicated in metabolic sensing. Contrastingly, in pups gestated and raised on short photoperiod (SP) and therefore following an over-wintering developmental trajectory with delayed sexual maturation, the ependymal region has fewer sensory tanycytes. Post-weaning transfer of SP-gestated pups to an intermediate photoperiod (IP), which accelerates reproductive maturation, results in a pronounced shift toward a ciliated tanycytic profile and formation of tanycytic processes. We suggest that tanycytic plasticity constitutes a mechanism to tailor metabolic development for extended survival in variable overwintering environments.

Efficient transfection of Atlantic salmon primary hepatocyte cells for functional assays and gene editing.

The expansion of genomic resources for Atlantic salmon over the past half decade has enabled efficient interrogation of genetic traits by large-scale correlation of genotype to phenotype. Moving from correlation to causation will require genotype-phenotype relationships to be tested experimentally in a cost-efficient and cell context-relevant manner. To enable such future experiments, we have developed a method for the isolation and genetic manipulation of primary hepatocytes from Atlantic salmon for use in heterologous expression, reporter assay, and gene editing experiments. We chose the liver as the tissue of interest because it is the metabolic hub and many current Atlantic salmon research projects focus on understanding metabolic processes to improve traits such as the growth rate, total fat content, and omega-3 content. We find that isolated primary hepatocytes are optimally transfected with both plasmid and ribonucleoprotein using a Neon electroporator at 1,400 V, 10 ms, and 2 pulses. Transfection efficiency with plasmid and cutting efficiency with ribonucleoprotein were optimally 46% and 60%, respectively. We also demonstrate a 26 times increase in luciferase expression under the promoter of the key liver metabolic gene, elovl5b, compared to an empty vector, in line with expected liver-specific expression. Taken together, this work provides a valuable resource enabling transfection and gene editing experiments in a context-relevant and cost-effective system.

Analyses of Genome Regulatory Evolution Following Whole-Genome Duplication Using the Phylogenetic EVE Model.

Whole-genome duplications (WGDs) are important in shaping the evolution of complex genomes, including rewiring of genome regulation. To address key questions about how WGDs impact the evolution of genome regulation, we need to understand the relative importance of selection versus drift and temporal evolutionary dynamics. One promising class of statistical models that can help address such questions are phylogenetic Ornstein-Uhlenbeck (OU) models.Here we present a computational pipeline for the comparative phylogenetic analyses of genome regulation using an OU model. We have implemented this model in R and provide a step-by-step protocol for the use of this model, including example scripts and simulated test data. We provide the nonspecialist a brief overview of how this model works and how to perform tests for signatures of selection on genome regulation as well as power simulations to aid in experimental design and interpretation of results. We believe that these resources could help polyploidy research move forward in an era of rapidly increasing functional genomics data across the tree of life.

Molecular Regulation of Biosynthesis of Long Chain Polyunsaturated Fatty Acids in Atlantic Salmon.

Salmon is a rich source of health-promoting omega-3 long chain polyunsaturated fatty acids (n-3 LC-PUFA), such as eicosapentaenoic acid (EPA, 20:5n-3) and docosahexaenoic acid (DHA, 22:6n-3). The LC-PUFA biosynthetic pathway in Atlantic salmon is one of the most studied compared to other teleosts. This has largely been due to the massive replacement of LC-PUFA-rich ingredients in aquafeeds with terrestrial plant oils devoid of these essential fatty acids (EFA) which ultimately pushed dietary content towards the minimal requirement of EFA. The practice would also reduce tissue content of n-3 LC-PUFA compromising the nutritional value of salmon to the human consumer. These necessitated detailed studies of endogenous biosynthetic capability as a contributor to these EFA. This review seeks to provide a comprehensive and concise overview of the current knowledge about the molecular genetics of PUFA biosynthesis in Atlantic salmon, highlighting the enzymology and nutritional regulation as well as transcriptional control networks. Furthermore, we discuss the impact of genome duplication on the complexity of salmon LC-PUFA pathway and highlight probable implications on endogenous biosynthetic capabilities. Finally, we have also compiled and made available a large RNAseq dataset from 316 salmon liver samples together with an R-script visualization resource to aid in explorative and hypothesis-driven research into salmon lipid metabolism.

SALARECON connects the Atlantic salmon genome to growth and feed efficiency.

Atlantic salmon (Salmo salar) is the most valuable farmed fish globally and there is much interest in optimizing its genetics and rearing conditions for growth and feed efficiency. Marine feed ingredients must be replaced to meet global demand, with challenges for fish health and sustainability. Metabolic models can address this by connecting genomes to metabolism, which converts nutrients in the feed to energy and biomass, but such models are currently not available for major aquaculture species such as salmon. We present SALARECON, a model focusing on energy, amino acid, and nucleotide metabolism that links the Atlantic salmon genome to metabolic fluxes and growth. It performs well in standardized tests and captures expected metabolic (in)capabilities. We show that it can explain observed hypoxic growth in terms of metabolic fluxes and apply it to aquaculture by simulating growth with commercial feed ingredients. Predicted limiting amino acids and feed efficiencies agree with data, and the model suggests that marine feed efficiency can be achieved by supplementing a few amino acids to plant- and insect-based feeds. SALARECON is a high-quality model that makes it possible to simulate Atlantic salmon metabolism and growth. It can be used to explain Atlantic salmon physiology and address key challenges in aquaculture such as development of sustainable feeds.

Genome-Wide Reconstruction of Rediploidization Following Autopolyploidization across One Hundred Million Years of Salmonid Evolution.

The long-term evolutionary impacts of whole-genome duplication (WGD) are strongly influenced by the ensuing rediploidization process. Following autopolyploidization, rediploidization involves a transition from tetraploid to diploid meiotic pairing, allowing duplicated genes (ohnologs) to diverge genetically and functionally. Our understanding of autopolyploid rediploidization has been informed by a WGD event ancestral to salmonid fishes, where large genomic regions are characterized by temporally delayed rediploidization, allowing lineage-specific ohnolog sequence divergence in the major salmonid clades. Here, we investigate the long-term outcomes of autopolyploid rediploidization at genome-wide resolution, exploiting a recent "explosion" of salmonid genome assemblies, including a new genome sequence for the huchen (Hucho hucho). We developed a genome alignment approach to capture duplicated regions across multiple species, allowing us to create 121,864 phylogenetic trees describing genome-wide ohnolog divergence across salmonid evolution. Using molecular clock analysis, we show that 61% of the ancestral salmonid genome experienced an initial "wave" of rediploidization in the late Cretaceous (85-106 Ma). This was followed by a period of relative genomic stasis lasting 17-39 My, where much of the genome remained tetraploid. A second rediploidization wave began in the early Eocene and proceeded alongside species diversification, generating predictable patterns of lineage-specific ohnolog divergence, scaling in complexity with the number of speciation events. Using gene set enrichment, gene expression, and codon-based selection analyses, we provide insights into potential functional outcomes of delayed rediploidization. This study enhances our understanding of delayed autopolyploid rediploidization and has broad implications for future studies of WGD events.

Chromosome-Level Genome Assembly of Chinese Sucker (Myxocyprinus asiaticus) Reveals Strongly Conserved Synteny Following a Catostomid-Specific Whole-Genome Duplication.

Fishes of the family Catostomidae ("suckers"; Teleostei: Cypriniformes) are hypothesized to have undergone an allopolyploidy event approximately 60 Ma. However, genomic evidence has previously been unavailable to assess this hypothesis. We sequenced and assembled the first chromosome-level catostomid genome, Chinese sucker (Myxocyprinus asiaticus), and present clear evidence of a catostomid-specific whole-genome duplication (WGD) event ("Cat-4R"). Our results reveal remarkably strong, conserved synteny since this duplication event, as well as between Myxocyprinus and an unduplicated outgroup, zebrafish (Danio rerio). Gene content and repetitive elements are also approximately evenly distributed across homeologous chromosomes, suggesting that both subgenomes retain some function, with no obvious bias in gene fractionation or subgenome dominance. The Cat-4R duplication provides another independent example of genome evolution following WGD in animals, in this case at the extreme end of conserved genome architecture over at least 25.2 Myr since the duplication. The M. asiaticus genome is a useful resource for researchers interested in understanding genome evolution following WGD in animals.

Comparative regulomics supports pervasive selection on gene dosage following whole genome duplication.

BACKGROUND: Whole genome duplication (WGD) events have played a major role in eukaryotic genome evolution, but the consequence of these extreme events in adaptive genome evolution is still not well understood. To address this knowledge gap, we used a comparative phylogenetic model and transcriptomic data from seven species to infer selection on gene expression in duplicated genes (ohnologs) following the salmonid WGD 80-100 million years ago. RESULTS: We find rare cases of tissue-specific expression evolution but pervasive expression evolution affecting many tissues, reflecting strong selection on maintenance of genome stability following genome doubling. Ohnolog expression levels have evolved mostly asymmetrically, by diverting one ohnolog copy down a path towards lower expression and possible pseudogenization. Loss of expression in one ohnolog is significantly associated with transposable element insertions in promoters and likely driven by selection on gene dosage including selection on stoichiometric balance. We also find symmetric expression shifts, and these are associated with genes under strong evolutionary constraints such as ribosome subunit genes. This possibly reflects selection operating to achieve a gene dose reduction while avoiding accumulation of "toxic mutations". Mechanistically, ohnolog regulatory divergence is dictated by the number of bound transcription factors in promoters, with transposable elements being one likely source of novel binding sites driving tissue-specific gains in expression. CONCLUSIONS: Our results imply pervasive adaptive expression evolution following WGD to overcome the immediate challenges posed by genome doubling and to exploit the long-term genetic opportunities for novel phenotype evolution.

Transkingdom network analysis provides insight into host-microbiome interactions in Atlantic salmon.

BACKGROUND: The Atlantic salmon gut constitutes an intriguing system for studying host-microbiota interactions due to the dramatic environmental change salmon experiences during its life cycle. Yet, little is known about the role of interactions in this system and there is a general deficit in computational methods for integrative analysis of omics data from host-microbiota systems. METHODS: We developed a pipeline to integrate host RNAseq data and microbial 16S rRNA amplicon sequencing data using weighted correlation network analysis. Networks are first inferred from each dataset separately, followed by module detections and finally robust identification of interactions via comparisons of representative module profiles. Through the use of module profiles, this network-based dimensionality reduction approach provides a holistic view into the discovery of potential host-microbiota symbionts. RESULTS: We analyzed host gene expression from the gut epithelial tissue and microbial abundances from the salmon gut in a long-term feeding trial spanning the fresh-/salt-water transition and including two feeds resembling the fatty acid compositions available in salt- and fresh-water environments, respectively. We identified several host modules with significant correlations to both microbiota modules and variables such as feed, growth and sex. Although the strongest associations largely coincided with the fresh-/salt-water transition, there was a second layer of correlations associating smaller host modules to both variables and microbiota modules. Hence, we identify extensive reprogramming of the gut epithelial transcriptome and large scale coordinated changes in gut microbiota composition associated with water type as well as evidence of host-microbiota interactions linked to feed.

Diversified regulation of circadian clock gene expression following whole genome duplication.

Across taxa, circadian control of physiology and behavior arises from cell-autonomous oscillations in gene expression, governed by a networks of so-called 'clock genes', collectively forming transcription-translation feedback loops. In modern vertebrates, these networks contain multiple copies of clock gene family members, which arose through whole genome duplication (WGD) events during evolutionary history. It remains unclear to what extent multiple copies of clock gene family members are functionally redundant or have allowed for functional diversification. We addressed this problem through an analysis of clock gene expression in the Atlantic salmon, a representative of the salmonids, a group which has undergone at least 4 rounds of WGD since the base of the vertebrate lineage, giving an unusually large complement of clock genes. By comparing expression patterns across multiple tissues, and during development, we present evidence for gene- and tissue-specific divergence in expression patterns, consistent with functional diversification of clock gene duplicates. In contrast to mammals, we found no evidence for coupling between cortisol and circadian gene expression, but cortisol mediated non-circadian regulated expression of a subset of clock genes in the salmon gill was evident. This regulation is linked to changes in gill function necessary for the transition from fresh- to sea-water in anadromous fish. Overall, this analysis emphasises the potential for a richly diversified clock gene network to serve a mixture of circadian and non-circadian functions in vertebrate groups with complex genomes.

The structural variation landscape in 492 Atlantic salmon genomes.

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.

Mid and hindgut transcriptome profiling analysis of Atlantic salmon (Salmon salar) under unpredictable chronic stress.

The intestinal epithelium is a selectively permeable barrier for nutrients, electrolytes and water, while maintaining effective protection against pathogens. Combinations of stressors throughout an animal's life, especially in agriculture and aquaculture settings, may affect the regular operativity of this organ with negative consequences for animal welfare. In the current study, we report the effects of a three-week unpredictable chronic stress (UCS) period on the intestinal morphology and transcriptome response of Atlantic salmon (Salmon salar) parr midgut and hindgut. Midgut and hindgut from both control and UCS fish were collected for histology and RNA-sequencing analysis to identify respective changes in the membrane structures and putative genes and pathways responding to UCS. Histological analysis did not show any significant effect on morphometric parameters. In the midgut, 1030 genes were differentially expressed following UCS, resulting in 279 genes which were involved in 13 metabolic pathways, including tissue repair pathways. In the hindgut, following UCS, 591 differentially expressed genes were detected with 426 downregulated and 165 upregulated. A total of 53 genes were related to three pathways. Downregulated genes include cellular senescence pathways, p53 signalling and cytokine-cytokine receptor pathways. The overall results corroborate that salmon parr were at least partly habituating to the UCS treatment. In midgut, the main upregulation was related to cell growth and repair, while in the hindgut there were indications of the activated apoptotic pathway, reduced cell repair and inhibited immune/anti-inflammatory capacity. This may be the trade-off between habituating to UCS and health resilience. This study suggests possible integrated genetic regulatory mechanisms that are tuned when farmed Atlantic salmon parr attempt to cope with UCS.

Liver slice culture as a model for lipid metabolism in fish.

Hepatic lipid metabolism is traditionally investigated in vitro using hepatocyte monocultures lacking the complex three-dimensional structure and interacting cell types essential liver function. Precision cut liver slice (PCLS) culture represents an alternative in vitro system, which benefits from retention of tissue architecture. Here, we present the first comprehensive evaluation of the PCLS method in fish (Atlantic salmon, Salmo salar L.) and validate it in the context of lipid metabolism using feeding trials, extensive transcriptomic data, and fatty acid measurements. We observe an initial period of post-slicing global transcriptome adjustment, which plateaued after 3 days in major metabolic pathways and stabilized through 9 days. PCLS fed alpha-linolenic acid (ALA) and insulin responded in a liver-like manner, increasing lipid biosynthesis gene expression. We identify interactions between insulin and ALA, where two PUFA biosynthesis genes that were induced by insulin or ALA alone, were highly down-regulated when insulin and ALA were combined. We also find that transcriptomic profiles of liver slices are exceedingly more similar to whole liver than hepatocyte monocultures, both for lipid metabolism and liver marker genes. PCLS culture opens new avenues for high throughput experimentation on the effect of "novel feed composition" and represent a promising new strategy for studying genotype-specific molecular features of metabolism.

The Chromosome-Level Genome Assembly of European Grayling Reveals Aspects of a Unique Genome Evolution Process Within Salmonids.

Salmonids represent an intriguing taxonomical group for investigating genome evolution in vertebrates due to their relatively recent last common whole genome duplication event, which occurred between 80 and 100 million years ago. Here, we report on the chromosome-level genome assembly of European grayling (Thymallus thymallus), which represents one of the earliest diverged salmonid subfamilies. To achieve this, we first generated relatively long genomic scaffolds by using a previously published draft genome assembly along with long-read sequencing data and a linkage map. We then merged those scaffolds by applying synteny evidence from the Atlantic salmon (Salmo salar) genome. Comparisons of the European grayling genome assembly to the genomes of Atlantic salmon and Northern pike (Esox lucius), the latter used as a nonduplicated outgroup, detailed aspects of the characteristic chromosome evolution process that has taken place in European grayling. While Atlantic salmon and other salmonid genomes are portrayed by the typical occurrence of numerous chromosomal fusions, European grayling chromosomes were confirmed to be fusion-free and were characterized by a relatively large proportion of paracentric and pericentric inversions. We further reported on transposable elements specific to either the European grayling or Atlantic salmon genome, on the male-specific sdY gene in the European grayling chromosome 11A, and on regions under residual tetrasomy in the homeologous European grayling chromosome pairs 9A-9B and 25A-25B. The same chromosome pairs have been observed under residual tetrasomy in Atlantic salmon and in other salmonids, suggesting that this feature has been conserved since the subfamily split.

Evolution of Cold Acclimation and Its Role in Niche Transition in the Temperate Grass Subfamily Pooideae.

The grass subfamily Pooideae dominates the grass floras in cold temperate regions and has evolved complex physiological adaptations to cope with extreme environmental conditions like frost, winter, and seasonality. One such adaptation is cold acclimation, wherein plants increase their frost tolerance in response to gradually falling temperatures and shorter days in the autumn. However, understanding how complex traits like cold acclimation evolve remains a major challenge in evolutionary biology. Here, we investigated the evolution of cold acclimation in Pooideae and found that a phylogenetically diverse set of Pooideae species displayed cold acclimation capacity. However, comparing differential gene expression after cold treatment in transcriptomes of five phylogenetically diverse species revealed widespread species-specific responses of genes with conserved sequences. Furthermore, we studied the correlation between gene family size and number of cold-responsive genes as well as between selection pressure on coding sequences of genes and their cold responsiveness. We saw evidence of protein-coding and regulatory sequence evolution as well as the origin of novel genes and functions contributing toward evolution of a cold response in Pooideae. Our results reflect that selection pressure resulting from global cooling must have acted on already diverged lineages. Nevertheless, conservation of cold-induced gene expression of certain genes indicates that the Pooideae ancestor may have possessed some molecular machinery to mitigate cold stress. Evolution of adaptations to seasonally cold climates is regarded as particularly difficult. How Pooideae evolved to transition from tropical to temperate biomes sheds light on how complex traits evolve in the light of climate changes.

The Grayling Genome Reveals Selection on Gene Expression Regulation after Whole-Genome Duplication.

Whole-genome duplication (WGD) has been a major evolutionary driver of increased genomic complexity in vertebrates. One such event occurred in the salmonid family ∼80 Ma (Ss4R) giving rise to a plethora of structural and regulatory duplicate-driven divergence, making salmonids an exemplary system to investigate the evolutionary consequences of WGD. Here, we present a draft genome assembly of European grayling (Thymallus thymallus) and use this in a comparative framework to study evolution of gene regulation following WGD. Among the Ss4R duplicates identified in European grayling and Atlantic salmon (Salmo salar), one-third reflect nonneutral tissue expression evolution, with strong purifying selection, maintained over ∼50 Myr. Of these, the majority reflect conserved tissue regulation under strong selective constraints related to brain and neural-related functions, as well as higher-order protein-protein interactions. A small subset of the duplicates have evolved tissue regulatory expression divergence in a common ancestor, which have been subsequently conserved in both lineages, suggestive of adaptive divergence following WGD. These candidates for adaptive tissue expression divergence have elevated rates of protein coding- and promoter-sequence evolution and are enriched for immune- and lipid metabolism ontology terms. Lastly, lineage-specific duplicate divergence points toward underlying differences in adaptive pressures on expression regulation in the nonanadromous grayling versus the anadromous Atlantic salmon. Our findings enhance our understanding of the role of WGD in genome evolution and highlight cases of regulatory divergence of Ss4R duplicates, possibly related to a niche shift in early salmonid evolution.

A systemic study of lipid metabolism regulation in salmon fingerlings and early juveniles fed plant oil.

In salmon farming, the scarcity of fish oil has driven a shift towards the use of plant-based oil from vegetable or seed, leading to fish feed low in long-chain PUFA (LC-PUFA) and cholesterol. Atlantic salmon has the capacity to synthesise both LC-PUFA and cholesterol, but little is known about the regulation of synthesis and how it varies throughout salmon life span. Here, we present a systemic view of lipid metabolism pathways based on lipid analyses and transcriptomic data from salmon fed contrasting diets of plant or fish oil from first feeding. We analysed four tissues (stomach, pyloric caeca, hindgut and liver) at three life stages (initial feeding 0·16 g, 2·5 g fingerlings and 10 g juveniles). The strongest response to diets higher in plant oil was seen in pyloric caeca of fingerlings, with up-regulation of thirty genes in pathways for cholesterol uptake, transport and biosynthesis. In juveniles, only eleven genes showed differential expression in pyloric caeca. This indicates a higher requirement of dietary cholesterol in fingerlings, which could result in a more sensitive response to plant oil. The LC-PUFA elongation and desaturation pathway was down-regulated in pyloric caeca, probably regulated by srebp1 genes. In liver, cholesterol metabolism and elongation and desaturation genes were both higher on plant oil. Stomach and hindgut were not notably affected by dietary treatment. Plant oil also had a higher impact on fatty acid composition of fingerlings compared with juveniles, suggesting that fingerlings have less metabolic regulatory control when primed with plant oil diet compared with juveniles.