Salmon Genome Research Articles

Atlantic salmon type I interferon genes revisited

Type I interferons (IFN–I) play a pivotal role in vertebrate innate immunity against viruses. This study is an analysis of IFN-I genes in an updated version of the Atlantic salmon genome published in 2021 (version Ssal_v3.1), revealing 47 IFN-I genes in the Atlantic salmon genome. The GH1 locus of chromosome (Chr) 3 harbors 9 IFNa genes, 5 IFNb genes, 6 IFNc genes, 11 IFNe genes and 1 IFNf gene. The GH2 locus on Chr6 contains 1 IFNa gene, 12 IFNc genes and 1 IFNf gene while Chr19 carries a single IFNd gene. Intraperitoneal injection of Atlantic salmon presmolts with poly I:C, a mimic of virus double-stranded RNA, significantly up-regulated IFNc genes from both Chr3 and Chr6 in heart, with lower expression in head kidney. IFNe expression increased in the heart, but not in the head kidney while IFNf was strongly up-regulated in both tissues. Antiviral activity of selected IFNs was assessed by transfection of salmon cells with IFN-expressing plasmids followed by infectious pancreatic necrosis virus infection, and by injection of fish with IFN-plasmids followed by measuring expression of the antiviral Mx1 gene. The results demonstrated that IFNc from both Chr3 and Chr6 provided full protection of cells against virus infection, whereas IFNe and IFNf showed lesser protection. IFNc from Chr3 and Chr6 along with IFNe and IFNf, up-regulated the Mx1 gene in the muscle, while only the IFNcs caused induction of Mx1 in liver. Overall, this study reveals that Atlantic salmon possesses an even more potent innate immune defense against viruses than previously understood.

Environmental influence on the Atlantic salmon transcriptome and methylome during sea lice infestations

The fish's immune response is affected by different factors, including a wide range of environmental conditions that can also disrupt or promote changes in the host-pathogen interactions. How environmental conditions modulate the salmon genome during parasitism is poorly understood here. This study aimed to explore the environmental influence on the Salmo salar transcriptome and methylome infected with the sea louse Caligus rogercresseyi. Atlantic salmon were experimentally infected with lice at two temperatures (8 and 16 °C) and salinity conditions (32 and 26PSU). Fish tissues were collected from the infected Atlantic salmon for reduced representation bisulfite sequencing (RRBS) and whole transcriptome sequencing (RNA-seq) analysis. The parasitic load was highly divergent in the evaluated environmental conditions, where the lowest lice abundance was observed in fish infected at 8 °C/26PSU. Notably, transcriptome profile differences were statistically associated with the number of alternative splicing events in fish exposed to low temperature/salinity conditions. Furthermore, the temperature significantly affected the methylation level, where high values of differential methylation regions were observed at 16 °C. Also, the association between expression levels of spliced transcripts and their methylation levels was determined, revealing significant correlations with Ferroptosis and TLR KEEG pathways. This study supports the relevance of the environmental conditions during host-parasite interactions in marine ecosystems. The discovery of alternative splicing transcripts associated with DMRs is also discussed as a novel player in fish biology.

Post‐glacial recolonization and multiple scales of secondary contact contribute to contemporary Atlantic salmon (Salmo salar) genomic variation in North America

AbstractAimIn northern environments, periods of isolation during Pleistocene glaciations and subsequent recolonization and secondary contact have had a significant influence on contemporary diversity of many species. The recent advent of high‐resolution genomic analyses allows unprecedented power to resolve genomic signatures of such events in northern species. Here, we provide the highest resolution genomic characterization of Atlantic salmon in North America to date to infer glacial refugia and the geographic scales of post‐glacial secondary contact.LocationNorth America.TaxonAtlantic salmon, Salmo salar.MethodsSamples were collected for 5455 individuals from 148 populations, encompassing the majority of the Atlantic salmon's native range in North America, from Labrador to Maine. Individuals were genotyped using a 220K single nucleotide polymorphism array aligned to the Atlantic salmon (Salmo salar) genome. Spatial genetic structure (principal component analysis, k‐means clustering, admixture) was evaluated in conjunction with genomic comparisons of these identified lineages to infer the refugia during the last glacial maximum and regions of secondary contact following recolonization.ResultsSpatial genomic analyses identified three phylogeographic groups, consistent with the northward recolonization from two southern glacial refugia in North America (a western Maritime lineage and an eastern Newfoundland and Labrador lineage), with subsequent differentiation of the eastern lineage into two separate groups. Secondary contact among these North American groups was observed within the northern Gulf of St. Lawrence, and evidence of trans‐Atlantic secondary contact was detected within the eastern Newfoundland and Labrador lineage. Comparison of groups from insular Newfoundland with those from mainland Labrador suggests genomic regions displaying high differentiation were characterized by elevated European admixture, suggesting a possible role of European secondary contact in population divergence.Main ConclusionsThese findings present the first evidence suggesting that genomic diversity in extant North American Atlantic salmon populations has resulted from allopatric isolation in two glacial refugia followed by both regional and trans‐Atlantic recolonization and secondary contact and demonstrate the power of genomic tools to resolve historical drivers of diversity in wild populations.

Investigating structural variant, indel and single nucleotide polymorphism differentiation between locally adapted Atlantic salmon populations.

Genomic structural variants (SVs) are now recognized as an integral component of intraspecific polymorphism and are known to contribute to evolutionary processes in various organisms. However, they are inherently difficult to detect and genotype from readily available short-read sequencing data, and therefore remain poorly documented in wild populations. Salmonid species displaying strong interpopulation variability in both life history traits and habitat characteristics, such as Atlantic salmon (Salmo salar), offer a prime context for studying adaptive polymorphism, but the contribution of SVs to fine-scale local adaptation has yet to be explored. Here, we performed a comparative analysis of SVs, single nucleotide polymorphisms (SNPs) and small indels (<50 bp) segregating in the Romaine and Puyjalon salmon, two putatively locally adapted populations inhabiting neighboring rivers (Québec, Canada) and showing pronounced variation in life history traits, namely growth, fecundity, and age at maturity and smoltification. We first catalogued polymorphism using a hybrid SV characterization approach pairing both short- (16X) and long-read sequencing (20X) for variant discovery with graph-based genotyping of SVs across 60 salmon genomes, along with characterization of SNPs and small indels from short reads. We thus identified 115,907 SVs, 8,777,832 SNPs and 1,089,321 short indels, with SVs covering 4.8 times more base pairs than SNPs. All three variant types revealed a highly congruent population structure and similar patterns of F ST and density variation along the genome. Finally, we performed outlier detection and redundancy analysis (RDA) to identify variants of interest in the putative local adaptation of Romaine and Puyjalon salmon. Genes located near these variants were enriched for biological processes related to nervous system function, suggesting that observed variation in traits such as age at smoltification could arise from differences in neural development. This study therefore demonstrates the feasibility of large-scale SV characterization and highlights its relevance for salmonid population genomics.

Highly efficient in vivo C-to-T base editing in Atlantic salmon (Salmo salar) – A step towards aquaculture precision breeding

The salmon industry has the potential to provide nutritious food to the expanding world population. However, the industry faces several sustainability challenges related to environmental impact and fish welfare. Genome editing presents a promising approach for precision breeding, where small, targeted changes to the fish genome may facilitate the introduction of advantageous traits like disease resistance, delayed maturation, or enhanced growth. This study presents the first application of base editing in an aquaculture fish species. By co-injecting cytidine base editor AncBE4max mRNA and gRNA targeting the pigmentation gene solute carrier family 45 member 2 in Atlantic salmon embryos, we were able to induce a specific C-to-T conversion, thereby introducing a premature stop codon and achieving an albino phenotype. High-throughput sequencing revealed an average of 50–66% correct conversion of the target base depending on the base editor mRNA concentration and up to 89% in a single individual. Our results suggest that AncBE4max provides a simple and efficient approach to making precise single nucleotide edits in the Atlantic salmon genome. Additionally, we used conventional CRISPR/Cas9 combined with oligonucleotide templates to insert point mutations at various positions up- and downstream of the Cas9 cleavage site, demonstrating that the insertion efficiency was affected by the position of the mutation. The application of genome editing tools that allow for precise changes in this species is a step towards aquaculture precision breeding and the introduction of traits promoting a more sustainable industry.

Genetic architecture of individual meiotic crossover rate and distribution in Atlantic Salmon

Meiotic recombination through chromosomal crossovers ensures proper segregation of homologous chromosomes during meiosis, while also breaking down linkage disequilibrium and shuffling alleles at loci located on the same chromosome. Rates of recombination can vary between species, but also between and within individuals, sex and chromosomes within species. Indeed, the Atlantic salmon genome is known to have clear sex differences in recombination with female biased heterochiasmy and markedly different landscapes of crossovers between males and females. In male meiosis, crossovers occur strictly in the telomeric regions, whereas in female meiosis crossovers tend to occur closer to the centromeres. However, little is known about the genetic control of these patterns and how this differs at the individual level. Here, we investigate genetic variation in individual measures of recombination in > 5000 large full-sib families of a Norwegian Atlantic salmon breeding population with high-density SNP genotypes. We show that females had 1.6 × higher crossover counts (CC) than males, with autosomal linkage maps spanning a total of 2174 cM in females and 1483 cM in males. However, because of the extreme telomeric bias of male crossovers, female recombination is much more important for generation of new haplotypes with 8 × higher intra-chromosomal genetic shuffling than males. CC was heritable in females (h2 = 0.11) and males (h2 = 0.10), and shuffling was also heritable in both sex but with a lower heritability in females (h2 = 0.06) than in males (h2 = 0.11). Inter-sex genetic correlations for both traits were close to zero, suggesting that rates and distribution of crossovers are genetically distinct traits in males and females, and that there is a potential for independent genetic change in both sexes in the Atlantic Salmon. Together, these findings give novel insights into the genetic architecture of recombination in salmonids and contribute to a better understanding of how rates and distribution of recombination may evolve in eukaryotes more broadly.

Basal and immune-responsive transcript expression of two Atlantic salmon interferon regulatory factor 2 (irf2) paralogues.

Atlantic salmon (Salmo salar) is one of the most economically important aquaculture species globally. However, disease has become a prevalent threat to this industry. A thorough understanding of the genes and molecular pathways involved in the immune responses of Atlantic salmon is imperative for selective breeding of disease-resistant broodstock, as well as developing new diets and vaccines to mitigate the impact of disease. Members of the interferon regulatory factor (IRF) family of transcription factors play roles in the induction of interferons and other cytokines involved in host immune responses to intracellular and parasitic pathogens. IRF family members also play diverse roles in other biological processes, such as stress response, reproduction and development. The current study focused on one member of the IRF family: interferon regulatory factor 2 (irf2). As previously shown, due to the genome duplication that occurred ∼80 million years ago in the salmonid lineage, there are two irf2 paralogues in the Atlantic salmon genome. In silico analyses at the cDNA and deduced amino acid levels were conducted followed by phylogenetic tree construction with IRF2 amino acid sequences from various ray-finned fishes, cartilaginous fish and tetrapods. qPCR was then used to analyze paralogue-specific irf2 constitutive expression across 17 adult tissues, as well as responses to the viral mimic pIC (i.e., synthetic double-stranded RNA analog) in cultured macrophage-like cells (in vitro) and to infection with the Gram-negative bacterium Moritella viscosa in skin samples (in vivo). The qPCR studies showed sex- and paralogue-specific differences in expression across tissues. For example, expression of both paralogues was higher in ovary than in testes; expression (considering both sexes together) was highest for irf2-1 in gonad and for irf2-2 in hindgut. Both irf2 paralogues were responsive to pIC stimulation, but varied in their induction level, with irf2-1 having an overall stronger response than irf2-2. Only one paralogue, irf2-2, was significantly responsive to M. viscosa infection. Differences in irf2-1 and irf2-2 transcript expression levels constitutively across tissues, and in response to pIC and M. viscosa, may suggest neo- or subfunctionalization of the duplicated genes. This novel information expands current knowledge and provides insight into how genome duplication events may impact host regulation of important immune markers.

Leukolectin-proteins are synthesized and secreted by lectocytes, a distinct category of fish embryonic mucus cells

Fish perivitelline fluid (PVF) is a vital extra-embryonic compartment. At hatching, PVF-contents dissolve into the hatching fluid (HF). Analysis of Atlantic salmon HF reveals nearly a hundred distinct proteins, most of which were identified by advanced mass-spectrometry. However, one entity with an apparent molecular weight 26 kDa, necessitated identification from its tryptic peptides. Subsequent cloning and sequencing revealed novel leukolectin-proteins. From bioinformatic analysis, leukolectins (LL) belong in the tectonin protein-family, with recognized functions in innate immunity. This study aims to identify LL-expressing cells in diverse fish species, and to characterize the LL-gene in order to predict bio-functions of leukolectins. LL-proteins were detected in HF from several fish species and one invertebrate, using polyclonal LL-specific IgGs. Embryonic LL-immunoreactive cells were numerous in Atlantic salmon, rainbow trout, fewer in Atlantic cod, and rare in Atlantic halibut and Oikopleura dioica. LL-immunoreactive cells were termed lectocytes, which corresponded to peridermal mucuscells stained by PAS, but unstained by eosin. Hence, lectocytes and hatching-gland cells were clearly distinguished. Northern blots revealed two salmon LL-transcripts at mid-embryogenesis. Such transcripts were detected in epithelial cells of the periderm, gills and oral cavity. LL-transcripts predominated in the periderm, while choriolysin-transcripts were dominant in the gills. No co-expression of choriolysins and LL-transcripts was detected. BAC-library screening yielded salmon LL's genestructure with 4 introns, 5 exons, TATA-box, multiple upstream putative transcription-factor bindingsites and polyadenylation site. LL-gene location on chromosome ssa17 was identified in Ssal_v3.1, the 2021version of the salmon genome. In conclusion, larvae from several fish species are outfitted with mucus enriched by LL-proteins. Mucus cells are present in embryos of all fishes, but embryonic lectocyte-numbers are far higher in species with near total larval survival. When (maternal) chorionic first-line immuno-defence is lost at hatching, leukolectin-enriched mucus may provide vital protection for larvae.

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.

Analysis of microbiota-host communication mediated by butyrate in Atlantic salmon

Butyrate is a microbiota-produced metabolite, sensed by host short-chain fatty acid receptors FFAR2 (Gpr43), FFAR3 (Gpr41), HCAR2 (Gpr109A), and Histone deacetylase (HDAC) that promotes microbiota-host crosstalk. Butyrate influences energy uptake, developmental and immune response in mammals. This microbial metabolite is produced by around 79 anaerobic genera present in the mammalian gut, yet little is known about the role of butyrate in the host-microbiota interaction in salmonid fish. To further our knowledge of this interaction, we analyzed the intestinal microbiota and genome of Atlantic salmon (Salmo salar), searching for butyrate-producing genera and host butyrate receptors. We identified Firmicutes, Proteobacteria, and Actinobacteria as the main butyrate-producing bacteria in the salmon gut microbiota. In the Atlantic salmon genome, we identified an expansion of genes orthologous to FFAR2 and HCAR2 receptors, and class I and IIa HDACs that are sensitive to butyrate. In addition, we determined the expression levels of orthologous of HCAR2 in the gut, spleen, and head-kidney, and FFAR2 in RTgutGC cells. The effect of butyrate on the Atlantic salmon immune response was evaluated by analyzing the pro and anti-inflammatory cytokines response in vitro in SHK-1 cells by RT-qPCR. Butyrate decreased the expression of the pro-inflammatory cytokine IL-1β and increased anti-inflammatory IL-10 and TGF-β cytokines. Butyrate also reduced the expression of interferon-alpha, Mx, and PKR, and decreased the viral load at a higher concentration (4 mM) in cells treated with this molecule before the infection with Infectious Pancreatic Necrosis Virus (IPNV) by mechanisms independent of FFAR2, FFAR3 and HCAR2 expression that probably inhibit HDAC. Moreover, butyrate modified phosphorylation of cytoplasmic proteins in RTgutGC cells. Our data allow us to infer that Atlantic salmon have the ability to sense butyrate produced by their gut microbiota via different specific targets, through which butyrate modulates the immune response of pro and anti-inflammatory cytokines and the antiviral response.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

A highly diverse set of novel immunoglobulin-like transcript (NILT) genes in zebrafish indicates a wide range of functions with complex relationships to mammalian receptors.

Multiple novel immunoglobulin-like transcripts (NILTs) have been identified from salmon, trout, and carp. NILTs typically encode activating or inhibitory transmembrane receptors with extracellular immunoglobulin (Ig) domains. Although predicted to provide immune recognition in ray-finned fish, we currently lack a definitive framework of NILT diversity, thereby limiting our predictions for their evolutionary origin and function. In order to better understand the diversity of NILTs and their possible roles in immune function, we identified five NILT loci in the Atlantic salmon (Salmo salar) genome, defined 86 NILT Ig domains within a 3-Mbp region of zebrafish (Danio rerio) chromosome 1, and described 41 NILT Ig domains as part of an alternative haplotype for this same genomic region. We then identified transcripts encoded by 43 different NILT genes which reflect an unprecedented diversity of Ig domain sequences and combinations for a family of non-recombining receptors within a single species. Zebrafish NILTs include a sole putative activating receptor but extensive inhibitory and secreted forms as well as membrane-bound forms with no known signaling motifs. These results reveal a higher level of genetic complexity, interindividual variation, and sequence diversity for NILTs than previously described, suggesting that this gene family likely plays multiple roles in host immunity.

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.

Functional characterization of Atlantic salmon (Salmo salar L.) PepT2 transporters.

The high‐affinity/low‐capacity system Slc15a2 (PepT2) is responsible for the reuptake of di/tripeptides from the renal proximal tubule, but it also operates in many other tissues and organs. Information regarding PepT2 in teleost fish is limited and, to date, functional data are available from the zebrafish (Danio rerio) only. Here, we report the identification of two slc15a2 genes in the Atlantic salmon (Salmo salar) genome, namely slc15a2a and slc15a2b. The two encoded PepT2 proteins share 87% identity and resemble both structurally and functionally the canonical vertebrate PepT2 system. The mRNA tissue distribution analyses reveal a widespread distribution of slc15a2a transcripts, being more abundant in the brain and gills, while slc15a2b transcripts are mainly expressed in the kidney and the distal part of the gastrointestinal tract. The function of the two transporters was investigated by heterologous expression in Xenopus laevis oocytes and two‐electrode voltage‐clamp recordings of transport and presteady‐state currents. Both PepT2a and PepT2b in the presence of Gly‐Gln elicit pH‐dependent and Na+ independent inward currents. The biophysical and kinetic analysis of the recorded currents defined the transport properties, confirming that the two Atlantic salmon PepT2 proteins behave as high‐affinity/low‐capacity transporters. The recent structures and the previous kinetic schemes of rat and human PepT2 qualitatively account for the characteristics of the two Atlantic salmon proteins. This study is the first to report on the functional expression of two PepT2‐type transporters that operate in the same vertebrate organism as a result of (a) gene duplication process(es).Key points Two slc15a2‐type genes, slc15a2a and slc15a2b coding for PepT2‐type peptide transporters were found in the Atlantic salmon. slc15a2a transcripts, widely distributed in the fish tissues, are abundant in the brain and gills, while slc15a2b transcripts are mainly expressed in the kidney and distal gastrointestinal tract.Amino acids involved in vertebrate Slc15 transport function are conserved in PepT2a and PepT2b proteins.Detailed kinetic analysis indicates that both PepT2a and PepT2b operate as high‐affinity transporters.The kinetic schemes and structures proposed for the mammalian models of PepT2 are suitable to explain the function of the two Atlantic salmon transporters.

The pink salmon genome: Uncovering the genomic consequences of a two-year life cycle.

Pink salmon (Oncorhynchus gorbuscha) adults are the smallest of the five Pacific salmon native to the western Pacific Ocean. Pink salmon are also the most abundant of these species and account for a large proportion of the commercial value of the salmon fishery worldwide. A two-year life history of pink salmon generates temporally isolated populations that spawn either in even-years or odd-years. To uncover the influence of this genetic isolation, reference genome assemblies were generated for each year-class and whole genome re-sequencing data was collected from salmon of both year-classes. The salmon were sampled from six Canadian rivers and one Japanese river. At multiple centromeres we identified peaks of Fst between year-classes that were millions of base-pairs long. The largest Fst peak was also associated with a million base-pair chromosomal polymorphism found in the odd-year genome near a centromere. These Fst peaks may be the result of a centromere drive or a combination of reduced recombination and genetic drift, and they could influence speciation. Other regions of the genome influenced by odd-year and even-year temporal isolation and tentatively under selection were mostly associated with genes related to immune function, organ development/maintenance, and behaviour.

Sulfate homeostasis in Atlantic salmon is associated with differential regulation of salmonid-specific paralogs in gill and kidney.

Sulfate (SO42‐) regulation is challenging for euryhaline species as they deal with large fluctuations of SO42‐ during migratory transitions between freshwater (FW) and seawater (SW), while maintaining a stable plasma SO42‐ concentration. Here, we investigated the regulation and potential role of sulfate transporters in Atlantic salmon during the preparative switch from SO42‐ uptake to secretion. A preparatory increase in kidney and gill sodium/potassium ATPase (Nka) enzyme activity during smolt development indicate preparative osmoregulatory changes. In contrast to gill Nka activity a transient decrease in kidney Nka after direct SW exposure was observed and may be a result of reduced glomerular filtration rates and tubular flow through the kidney. In silico analyses revealed that Atlantic salmon genome comprises a single slc13a1 gene and additional salmonid‐specific duplications of slc26a1 and slc26a6a, leading to new paralogs, namely the slc26a1a and ‐b, and slc26a6a1 and ‐a2. A kidney‐specific increase in slc26a6a1 and slc26a1a during smoltification and SW transfer, suggests an important role of these sulfate transporters in the regulatory shift from absorption to secretion in the kidney. Plasma SO42‐ in FW smolts was 0.70 mM, followed by a transient increase to 1.14 ± 0.33 mM 2 days post‐SW transfer, further decreasing to 0.69 ± 0.041 mM after 1 month in SW. Our findings support the vital role of the kidney in SO42‐ excretion through the upregulated slc26a6a1, the most likely secretory transport candidate in fish, which together with the slc26a1a transporter likely removes excess SO42‐, and ultimately enable the regulation of normal plasma SO42‐ levels in SW.

Regional Expression of npy mRNA Paralogs in the Brain of Atlantic Salmon (Salmo salar, L.) and Response to Fasting.

Neuropeptide Y (NPY) is known as a potent orexigenic signal in vertebrates, but its role in Atlantic salmon has not yet been fully established. In this study, we identified three npy paralogs, named npya1, npya2, and npyb, in the Atlantic salmon genome. In silico analysis revealed that these genes are well conserved across the vertebrate’s lineage and the mature peptide sequences shared at least 77% of identity with the human homolog. We analyzed mRNA expression of npy paralogs in eight brain regions of Atlantic salmon post-smolt, and the effect of 4 days of fasting on the npy expression level. Results show that npya1 was the most abundant paralog, and was predominantly expressed in the telencephalon, followed by the midbrain and olfactory bulb. npya2 mRNA was highly abundant in hypothalamus and midbrain, while npyb was found to be highest expressed in the telencephalon, with low mRNA expression levels detected in all the other brain regions. 4 days of fasting resulted in a significant (p < 0.05) decrease of npya1 mRNA expression in the olfactory bulb, increased npya2 mRNA expression in the midbrain and decreased npyb mRNA expression in the pituitary. In the hypothalamus, the vertebrate appetite center, expression of the npy paralogs was not significantly affected by feeding status. However, we observed a trend of increased npya2 mRNA expression (p = 0.099) following 4 days of fasting. Altogether, our findings provide a solid basis for further research on appetite and energy metabolism in Atlantic salmon.

A highly diverse set of novel immunoglobulin-like transcript (NILT) genes in zebrafish indicates a wide range of functions

Abstract Multiple novel immunoglobulin-like transcripts (NILTs) have been identified from salmon, trout and carp. NILTs typically encode activating or inhibitory transmembrane receptors with extracellular immunoglobulin (Ig) domains. Although predicted to provide some level of immune recognition in ray-finned fish, we currently lack a definitive framework of NILT diversity, thereby challenging our ability to understand the evolutionary origin and specific function of these genes. In order to better understand the diversity of NILT genes and their possible roles in immune function, we identified 5 NILT loci in the Atlantic salmon (Salmo salar) genome, defined 86 NILT Ig domains within a 3 Mbp region of zebrafish (Danio rerio) chromosome 1, and described 41 NILT Ig domains as part of an alternate haplotype for this same genomic region. We then identified transcripts that define 43 different NILT genes and reflect an unprecedented range of sequence diversity and combinatorial diversity of Ig domains for a single family of non-recombining receptors. Zebrafish NILTs include a single putative activating receptor but extensive inhibitory and secreted forms as well as membrane-bound forms with no known signaling motifs. Whole genome sequences from different genetic backgrounds reveal at least three alternative haplotypes of the zebrafish NILT gene cluster that display gene content variation. This observation indicates that different individual fish can encode different combinations of NILTs which could impact their immune function. Furthermore, these results reveal a higher level of genetic complexity and sequence diversity for NILTs than previously described, suggesting that this gene family likely plays multiple roles in host immunity.

A de novo Full-Length mRNA Transcriptome Generated From Hybrid-Corrected PacBio Long-Reads Improves the Transcript Annotation and Identifies Thousands of Novel Splice Variants in Atlantic Salmon.

Atlantic salmon (Salmo salar) is a major species produced in world aquaculture and an important vertebrate model organism for studying the process of rediploidization following whole genome duplication events (Ss4R, 80 mya). The current Salmo salar transcriptome is largely generated from genome sequence based in silico predictions supported by ESTs and short-read sequencing data. However, recent progress in long-read sequencing technologies now allows for full-length transcript sequencing from single RNA-molecules. This study provides a de novo full-length mRNA transcriptome from liver, head-kidney and gill materials. A pipeline was developed based on Iso-seq sequencing of long-reads on the PacBio platform (HQ reads) followed by error-correction of the HQ reads by short-reads from the Illumina platform. The pipeline successfully processed more than 1.5 million long-reads and more than 900 million short-reads into error-corrected HQ reads. A surprisingly high percentage (32%) represented expressed interspersed repeats, while the remaining were processed into 71 461 full-length mRNAs from 23 071 loci. Each transcript was supported by several single-molecule long-read sequences and at least three short-reads, assuring a high sequence accuracy. On average, each gene was represented by three isoforms. Comparisons to the current Atlantic salmon transcripts in the RefSeq database showed that the long-read transcriptome validated 25% of all known transcripts, while the remaining full-length transcripts were novel isoforms, but few were transcripts from novel genes. A comparison to the current genome assembly indicates that the long-read transcriptome may aid in improving transcript annotation as well as provide long-read linkage information useful for improving the genome assembly. More than 80% of transcripts were assigned GO terms and thousands of transcripts were from genes or splice-variants expressed in an organ-specific manner demonstrating that hybrid error-corrected long-read transcriptomes may be applied to study genes and splice-variants expressed in certain organs or conditions (e.g., challenge materials). In conclusion, this is the single largest contribution of full-length mRNAs in Atlantic salmon. The results will be of great value to salmon genomics research, and the pipeline outlined may be applied to generate additional de novo transcriptomes in Atlantic Salmon or applied for similar projects in other species.

Genetic Drift Dominates Genome-Wide Regulatory Evolution Following an Ancient Whole-Genome Duplication in Atlantic Salmon.

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.