The rainbow trout genome provides novel insights into evolution after whole-genome duplication in vertebrates

Migration-related phenotypic divergence is associated with epigenetic modifications in rainbow trout

Migration is essential for the reproduction and survival of many animals, yet little is understood about its underlying molecular mechanisms. We used the salmonid Oncorhynchus mykiss to gain mechanistic insight into smoltification, which is a morphological, physiological and behavioural transition undertaken by juveniles in preparation for seaward migration. O. mykiss is experimentally tractable and displays intra- and interpopulation variation in migration propensity. Migratory individuals can produce nonmigratory progeny and vice versa, indicating a high degree of phenotypic plasticity. One potential way that phenotypic plasticity might be linked to variation in migration-related life history tactics is through epigenetic regulation of gene expression. To explore this, we quantitatively measured genome-scale DNA methylation in fin tissue using reduced representation bisulphite sequencing of F2 siblings produced from a cross between steelhead (migratory) and rainbow trout (nonmigratory) lines. We identified 57 differentially methylated regions (DMRs) between smolt and resident O. mykiss juveniles. DMRs were high in magnitude, with up to 62% differential methylation between life history types, and over half of the gene-associated DMRs were in transcriptional regulatory regions. Many of the DMRs encode proteins with activity relevant to migration-related transitions (e.g. circadian rhythm pathway, nervous system development, protein kinase activity). This study provides the first evidence of a relationship between epigenetic variation and life history divergence associated with migration-related traits in any species.

Epithelial Cadherin Determines Resistance to Infectious Pancreatic Necrosis Virus in Atlantic Salmon

Infectious pancreatic necrosis virus (IPNV) is the cause of one of the most prevalent diseases in farmed Atlantic salmon (Salmo salar). A quantitative trait locus (QTL) has been found to be responsible for most of the genetic variation in resistance to the virus. Here we describe how a linkage disequilibrium-based test for deducing the QTL allele was developed, and how it was used to produce IPN-resistant salmon, leading to a 75% decrease in the number of IPN outbreaks in the salmon farming industry. Furthermore, we describe how whole-genome sequencing of individuals with deduced QTL genotypes was used to map the QTL down to a region containing an epithelial cadherin (cdh1) gene. In a coimmunoprecipitation assay, the Cdh1 protein was found to bind to IPNV virions, strongly indicating that the protein is part of the machinery used by the virus for internalization. Immunofluorescence revealed that the virus colocalizes with IPNV in the endosomes of homozygous susceptible individuals but not in the endosomes of homozygous resistant individuals. A putative causal single nucleotide polymorphism was found within the full-length cdh1 gene, in phase with the QTL in all observed haplotypes except one; the absence of a single, all-explaining DNA polymorphism indicates that an additional causative polymorphism may contribute to the observed QTL genotype patterns. Cdh1 has earlier been shown to be necessary for the internalization of certain bacteria and fungi, but this is the first time the protein is implicated in internalization of a virus.

Discovery and characterization of single nucleotide polymorphisms in coho salmon, Oncorhynchus kisutch

Molecular population genetic analyses have become an integral part of ecological investigation and population monitoring for conservation and management. Microsatellites have been the molecular marker of choice for such applications over the last several decades, but single nucleotide polymorphism (SNP) markers are rapidly expanding beyond model organisms. Coho salmon (Oncorhynchus kisutch) is native to the north Pacific Ocean and its tributaries, where it is the focus of intensive fishery and conservation activities. As it is an anadromous species, coho salmon typically migrate across multiple jurisdictional boundaries, complicating management and requiring shared data collection methods. Here, we describe the discovery and validation of a suite of novel SNPs and associated genotyping assays which can be used in the genetic analyses of this species. These assays include 91 that are polymorphic in the species and one that discriminates it from a sister species, Chinook salmon. We demonstrate the utility of these SNPs for population assignment and phylogeographic analyses, and map them against the draft trout genome. The markers constitute a large majority of all SNP markers described for coho salmon and will enable both population- and pedigree-based analyses across the southern part of the species native range.

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

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

Identification of single-nucleotide polymorphism markers associated with cortisol response to crowding in rainbow trout

Understanding stress responses is essential for improving animal welfare and increasing agriculture production efficiency. Previously, we reported microsatellite markers associated with quantitative trait loci (QTL) affecting plasma cortisol response to crowding in rainbow trout. In this study, our main objectives were to identify single-nucleotide polymorphism (SNP) markers associated with cortisol response to crowding in rainbow trout using both GWAS (genome-wide association studies) and QTL mapping methods and to employ rapidly expanding genomic resources for rainbow trout toward the identification of candidate genes affecting this trait. A three-generation F2 mapping family (2008052) was genotyped using RAD-seq (restriction-site-associated DNA sequencing) to identify 4874 informative SNPs. GWAS identified 26 SNPs associated with cortisol response to crowding whereas QTL mapping revealed two significant QTL on chromosomes Omy8 and Omy12, respectively. Positional candidate genes were identified using marker sequences to search the draft genome assembly of rainbow trout. One of the genes in the QTL interval on Omy12 is a putative serine/threonine protein kinase gene that was differentially expressed in the liver in response to handling and confinement stress in our previous study. A homologue of this gene was differentially expressed in zebrafish embryos exposed to diclofenac, a nonsteroidal anti-inflammatory drug (NSAID) and an environmental toxicant. NSAIDs have been shown to affect the cortisol response in rainbow trout; therefore, this gene is a good candidate based on its physical position and expression. However, the reference genome resources currently available for rainbow trout require continued improvement as demonstrated by the unmapped SNPs and the putative assembly errors detected in this study.

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

Species and populations with parallel evolution of specific traits can help illuminate how predictable adaptations and divergence are at the molecular and developmental level. Following the last glacial period, dwarfism and specialized bottom feeding morphology evolved rapidly in several landlocked Arctic charr Salvelinus alpinus populations in Iceland.  

To study the genetic divergence between small benthic morphs and limnetic morphs, we conducted RNA-sequencing charr embryos at four stages in early development. We studied two stocks with contrasting morphologies: the small benthic (SB) charr from Lake Thingvallavatn and Holar aquaculture (AC) charr.

The data reveal significant differences in expression of several biological pathways during charr development. There was also an expression difference between SB- and AC-charr in genes involved in energy metabolism and blood coagulation genes. We confirmed differing expression of five genes in whole embryos with qPCR, including lysozyme and natterin-like which was previously identified as a fish-toxin of a lectin family that may be a putative immunopeptide. We also verified differential expression of 7 genes in the developing head that associated consistently with benthic v.s.limnetic morphology (studied in 4 morphs). Comparison of single nucleotide polymorphism (SNP) frequencies reveals extensive genetic differentiation between the SB and AC-charr (~1300 with more than 50% frequency difference). Curiously, three derived alleles in the otherwise conserved 12s and 16s mitochondrial ribosomal RNA genes are found in benthic charr.

The data implicate multiple genes and molecular pathways in divergence of small benthic charr and/or the response of aquaculture charr to domestication. Functional, genetic and population genetic studies on more freshwater and anadromous populations are needed to confirm the specific loci and mutations relating to specific ecological traits in Arctic charr.

Genome-wide identification, evolution and expression analysis of nuclear receptor superfamily in Nile tilapia, Oreochromis niloticus

The nuclear receptor (NR) superfamily, which is divided into 7 subfamilies, constitutes one of the largest classes of transcription factors. In this study, through comprehensive database search, we identified all NRs (including 4 novel members) from the tilapia (75), common carp (137), zebrafish (73), fugu (73), tetraodon (72), stickleback (70), medaka (69), coelacanth (55), spotted gar (51) and elephant shark (50). For 21 NRs, two duplicates were found in teleosts, while only one in tetrapods. These duplicates, except those of DAX1, SHP and GCNF found in the elephant shark, were derived from 3R (third round of genome duplication). The linkage duplication of 5 syntenic blocks (comprising 14 duplicated NR couples) in teleosts further supported their 3R origin. Based on transcriptome data from adult tilapia, 53 NRs were found to be expressed in more than one tissue (brain, head kidney, heart, liver, kidney, muscle, ovary and testis), and 4 were tissue-specific, indicating their essential roles in the corresponding tissue. Based on the XX and XY gonadal transcriptome data from four developmental stages, 65 NRs were detected in gonads, with 21, 31, 11 and 29 expressed sexual dimorphically at 5, 30, 90 and 180days after hatching, respectively. The expression of four selected genes was examined by in situ hybridization (ISH) and quantitative PCR (qPCR) to validate the spatial and temporal expression profiles of NRs. Comparative analyses of the expression profiles of duplicated NRs revealed divergence in gene expression as well as gene function. Our results demonstrated that NRs may play important roles in sex determination and gonadal development in teleosts.

Mapping and validation of a major QTL affecting resistance to pancreas disease (salmonid alphavirus) in Atlantic salmon (Salmo salar)

Pancreas disease (PD), caused by a salmonid alphavirus (SAV), has a large negative economic and animal welfare impact on Atlantic salmon aquaculture. Evidence for genetic variation in host resistance to this disease has been reported, suggesting that selective breeding may potentially form an important component of disease control. The aim of this study was to explore the genetic architecture of resistance to PD, using survival data collected from two unrelated populations of Atlantic salmon; one challenged with SAV as fry in freshwater (POP 1) and one challenged with SAV as post-smolts in sea water (POP 2). Analyses of the binary survival data revealed a moderate-to-high heritability for host resistance to PD in both populations (fry POP 1 h²~0.5; post-smolt POP 2 h²~0.4). Subsets of both populations were genotyped for single nucleotide polymorphism markers, and six putative resistance quantitative trait loci (QTL) were identified. One of these QTL was mapped to the same location on chromosome 3 in both populations, reaching chromosome-wide significance in both the sire- and dam-based analyses in POP 1, and genome-wide significance in a combined analysis in POP 2. This independently verified QTL explains a significant proportion of host genetic variation in resistance to PD in both populations, suggesting a common underlying mechanism for genetic resistance across lifecycle stages. Markers associated with this QTL are being incorporated into selective breeding programs to improve PD resistance.

Divergent functions of fibroblast growth factor receptor-like 1 genes in grass carp (Ctenopharyngodon idella)

Fibroblast growth factor receptor-like 1 (FGFRL1) is a novel FGF receptor (FGFR) lacking an intracellular tyrosine kinase domain. FGFRs control the proliferation, differentiation and migration of cells in various tissues. However the functions of FGFRL1 in teleost fish are currently unknown. In this study, we report the identification of two fgfrl1 genes in grass carp (Ctenopharyngodon idella) that share 56% amino acid sequence identity. Both fgfrl1a and 1b were transcribed throughout embryogenesis, and mRNA levels were particularly high during somitogenesis. Using in situ hybridization, fgfrl1a transcripts were detected in notochord, somites, brain and eye at 14, 24 and 36 h post fertilization (hpf). In contrast, fgfrl1b was transcribed mainly in the endoderm at 14 hpf, in the gut and proctodeum at 24 hpf, and in the lens, pharyngeal arch and proctodeum at 36 hpf. In adult fish, fgfrl1a was abundantly expressed in heart, brain and muscle, while fgfrl1b was expressed strongly in eye, muscle and gill. Furthermore, both genes were significantly (p<0.05) up-regulated in muscle and brain during starvation and returned to normal levels rapidly after re-feeding. Exogenous treatment with different doses of human growth hormone down-regulated the expression of both genes in brain and muscle (p<0.05). These results suggest that Fgfrl1a and 1b play divergent roles in regulating growth and development in grass carp.

The genetic architecture of growth and fillet traits in farmed Atlantic salmon (Salmo salar)

BACKGROUND

Performance and quality traits such as harvest weight, fillet weight and flesh color are of economic importance to the Atlantic salmon aquaculture industry. The genetic factors underlying these traits are of scientific and commercial interest. However, such traits are typically polygenic in nature, with the number and size of QTL likely to vary between studies and populations. The aim of this study was to investigate the genetic basis of several growth and fillet traits measured at harvest in a large farmed salmon population by using SNP markers. Due to the marked heterochiasmy in salmonids, an efficient two-stage mapping approach was applied whereby QTL were detected using a sire-based linkage analysis, a sparse SNP marker map and exploiting low rates of recombination, while a subsequent dam-based analysis focused on the significant chromosomes with a denser map to confirm QTL and estimate their position.

RESULTS

The harvest traits all showed significant heritability, ranging from 0.05 for fillet yield up to 0.53 for the weight traits. In the sire-based analysis, 1695 offspring with trait records and their 20 sires were successfully genotyped for the SNPs on the sparse map. Chromosomes 13, 18, 19 and 20 were shown to harbor genome-wide significant QTL affecting several growth-related traits. The QTL on chr. 13, 18 and 20 were detected in the dam-based analysis using 512 offspring from 10 dams and explained approximately 6-7 % of the within-family variation in these traits.

CONCLUSIONS

We have detected several QTL affecting economically important complex traits in a commercial salmon population. Overall, the results suggest that the traits are relatively polygenic and that QTL tend to be pleiotropic (affecting the weight of several components of the harvested fish). Comparison of QTL regions across studies suggests that harvest trait QTL tend to be relatively population-specific. Therefore, the application of marker or genomic selection for improvement in these traits is likely to be most effective when the discovery population is closely related to the selection candidates (e.g. within-family genomic selection).

A comparison between egg trancriptomes of cod and salmon reveals species-specific traits in eggs for each species

Fish in use in aquaculture display large variation in gamete biology. To reach better understanding around this issue, this study aims at identifying if species specific "egg life history traits" can be hidden in the unfertilized egg. This was done by investigating egg transcriptome differences between Atlantic salmon and Atlantic cod. Salmon and cod eggs were selected due to their largely differencing phenotypes. An oligo microarray analysis was performed on ovulated eggs from cod (n = 8) and salmon (n = 7). The arrays were normalized to a similar spectrum for both arrays. Both arrays were re-annotated with SWISS-Prot and KEGG genes to retrieve an official gene symbol and an orthologous KEGG annotation, in salmon and cod arrays this represented 14,009 and 7,437 genes respectively. The probe linked to the highest gene expression for that particular KEGG annotation was used to compare expression between species. Differential expression was calculated for genes that had an annotation with score >300, resulting in a total of 2,457 KEGG annotations (genes) being differently expressed between the species (FD > 2). This analysis revealed that immune, signal transduction and excretory related pathways were overrepresented in salmon compared to cod. The most overrepresented pathways in cod were related to regulation of genetic information processing and metabolism. To conclude this analysis clearly point at some distinct transcriptome repertoires for cod and salmon and that these differences may explain some of the species-specific biological features for salmon and cod eggs.

New insights into the nutritional regulation of gluconeogenesis in carnivorous rainbow trout (Oncorhynchus mykiss): a gene duplication trail

The rainbow trout (Oncorhynchus mykiss) is considered to be a strictly carnivorous fish species that is metabolically adapted for high catabolism of proteins and low utilization of dietary carbohydrates. This species consequently has a "glucose-intolerant" phenotype manifested by persistent hyperglycemia when fed a high-carbohydrate diet. Gluconeogenesis in adult fish is also poorly, if ever, regulated by carbohydrates, suggesting that this metabolic pathway is involved in this specific phenotype. In this study, we hypothesized that the fate of duplicated genes after the salmonid-specific 4th whole genome duplication (Ss4R) may have led to adaptive innovation and that their study might provide new elements to enhance our understanding of gluconeogenesis and poor dietary carbohydrate use in this species. Our evolutionary analysis of gluconeogenic genes revealed that pck1, pck2, fbp1a, and g6pca were retained as singletons after Ss4r, while g6pcb1, g6pcb2, and fbp1b ohnolog pairs were maintained. For all genes, duplication may have led to sub- or neofunctionalization. Expression profiles suggest that the gluconeogenesis pathway remained active in trout fed a no-carbohydrate diet. When trout were fed a high-carbohydrate diet (30%), most of the gluconeogenic genes were non- or downregulated, except for g6pbc2 ohnologs, whose RNA levels were surprisingly increased. This study demonstrates that Ss4R in trout involved adaptive innovation via gene duplication and via the outcome of the resulting ohnologs. Indeed, maintenance of ohnologous g6pcb2 pair may contribute in a significant way to the glucose-intolerant phenotype of trout and may partially explain its poor use of dietary carbohydrates.

Temperature during early development has long-term effects on microRNA expression in Atlantic cod

Background

Environmental temperature has serious implications in life cycle of aquatic ectotherms. Understanding the molecular mechanisms of temperature acclimation and adaptation of marine organisms is of the uttermost importance for ecology, fisheries, and aquaculture, as it allows modeling the effects of global warming on population dynamics. Regulatory molecules are major modulators of acclimation and adaptation; among them, microRNAs (miRNAs) are versatile and substantial contributors to regulatory networks of development and adaptive plasticity. However, their role in thermal plasticity is poorly known. We have asked whether the temperature and its shift during the early ontogeny (embryonic and larval development) affect the miRNA repertoire of Atlantic cod (Gadus morhua), and if thermal experience has long-term consequences in the miRNA profile.

Results

We characterized miRNA during different developmental stages and in juvenile tissues using next generation sequencing. We identified 389 putative miRNA precursor loci, 120 novel precursor miRNAs, and 281 mature miRNAs. Some miRNAs showed stage- or tissue-enriched expression and miRNAs, such as the miR-17 ~ 92 cluster, myomiRs (miR-206), neuromiRs (miR-9, miR-124), miR-130b, and miR-430 showed differential expression in different temperature regimes. Long-term effect of embryonic incubation temperature was revealed on expression of some miRNAs in juvenile pituitary (miR-449), gonad (miR-27c, miR-30c, and miR-200a), and liver (let-7 h, miR-7a, miR-22, miR-34c, miR-132a, miR-192, miR-221, miR-451, miR-2188, and miR-7550), but not in brain. Some of differentially expressed miRNAs in the liver were confirmed using LNA-based rt-qPCR. The effect of temperature on methylation status of selected miRNA promoter regions was mostly inconclusive.

Conclusions

Temperature elevation by several degrees during embryonic and larval developmental stages significantly alters the miRNA profile, both short-term and long-term. Our results suggest that a further rise in seas temperature might affect life history of Atlantic cod.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1503-7) contains supplementary material, which is available to authorized users.

Functional characterization of a BCL10 isoform in the rainbow trout Oncorhynchus mykiss

•

The complexes formed by BCL10, MALT1 and CARMA proteins are key regulators of NF-κB activation.

•

We report the functional characterization of tBCL10, a BCL10 isoform from the trout Oncorhynchus mykiss.

•

tBCL10 can functionally replace the human protein.

•

The rainbow trout Oncorhynchus mykiss can serve as a model organism to study this pathway.

The complexes formed by BCL10, MALT1 and members of the family of CARMA proteins have recently been the focus of much attention because they represent a key mechanism for regulating activation of the transcription factor NF-κB. Here, we report the functional characterization of a novel isoform of BCL10 in the trout Oncorhynchus mykiss, which we named tBCL10. tBCL10 dimerizes, binds to components of the CBM complex and forms cytoplasmic filaments. Functionally, tBCL10 activates NF-κB transcription factor and is inhibited by the deubiquitinating enzyme A20. Finally, depletion experiments indicate that tBCL10 can functionally replace the human protein. This work demonstrates the evolutionary conservation of the mechanism of NF-κB activation through the CBM complex, and indicates that the rainbow trout O.mykiss can serve as a model organism to study this pathway.

Identification and expression analysis of two interleukin-23α (p19) isoforms, in rainbow trout Oncorhynchus mykiss and Atlantic salmon Salmo salar

Interleukin (IL)-23 is a heterodimeric IL-12 family cytokine composed of a p19 α-chain, linked to a p40 β-chain that is shared with IL-12. IL-23 is distinguished functionally from IL-12 by its ability to induce the production of IL-17, and differentiation of Th17 cells in mammals. Three isoforms of p40 (p40a, p40b and p40c) have been found in some 3R teleosts. Salmonids also possess three p40 isoforms (p40b1, p40b2 and p40c) although p40a is missing, and two copies (paralogues) of p40b are present that have presumably been retained following the 4R duplication in this fish lineage. Teleost p19 has been discovered recently in zebrafish, but to date there is limited information on expression and modulation of this molecule. In this report we have cloned two p19 paralogues (p19a and p19b) in salmonids, suggesting that a salmonid can possess six potential IL-23 isoforms. Whilst Atlantic salmon has two active p19 genes, the rainbow trout p19b gene may have been pseudogenized. The salmonid p19 translations share moderate identities (22.8-29.9%) to zebrafish and mammalian p19 molecules, but their identity was supported by structural features, a conserved 4 exon/3 intron gene organisation, and phylogenetic tree analysis. The active salmonid p19 genes are highly expressed in blood and gonad. Bacterial (Yersinia ruckeri) and viral infection in rainbow trout induces the expression of p19a, suggesting pathogen-specific induction of IL-23 isoforms. Trout p19a expression was also induced by PAMPs (poly IC and peptidoglycan) and the proinflammatory cytokine IL-1β in primary head kidney macrophages. These data may indicate diverse functional roles of trout IL-23 isoforms in regulating the immune response in fish.

Transcriptome assembly, gene annotation and tissue gene expression atlas of the rainbow trout

Efforts to obtain a comprehensive genome sequence for rainbow trout are ongoing and will be complemented by transcriptome information that will enhance genome assembly and annotation. Previously, transcriptome reference sequences were reported using data from different sources. Although the previous work added a great wealth of sequences, a complete and well-annotated transcriptome is still needed. In addition, gene expression in different tissues was not completely addressed in the previous studies. In this study, non-normalized cDNA libraries were sequenced from 13 different tissues of a single doubled haploid rainbow trout from the same source used for the rainbow trout genome sequence. A total of ~1.167 billion paired-end reads were de novo assembled using the Trinity RNA-Seq assembler yielding 474,524 contigs > 500 base-pairs. Of them, 287,593 had homologies to the NCBI non-redundant protein database. The longest contig of each cluster was selected as a reference, yielding 44,990 representative contigs. A total of 4,146 contigs (9.2%), including 710 full-length sequences, did not match any mRNA sequences in the current rainbow trout genome reference. Mapping reads to the reference genome identified an additional 11,843 transcripts not annotated in the genome. A digital gene expression atlas revealed 7,678 housekeeping and 4,021 tissue-specific genes. Expression of about 16,000–32,000 genes (35–71% of the identified genes) accounted for basic and specialized functions of each tissue. White muscle and stomach had the least complex transcriptomes, with high percentages of their total mRNA contributed by a small number of genes. Brain, testis and intestine, in contrast, had complex transcriptomes, with a large numbers of genes involved in their expression patterns. This study provides comprehensive de novo transcriptome information that is suitable for functional and comparative genomics studies in rainbow trout, including annotation of the genome.

MicroRNA expression profiles from eggs of different qualities associated with post-ovulatory ageing in rainbow trout (Oncorhynchus mykiss)

BACKGROUND

Egg quality is an important aspect in rainbow trout farming. Post-ovulatory aging is one of the most important factors affecting egg quality. MicroRNAs (miRNAs) are the major regulators in various biological processes and their expression profiles could serve as reliable biomarkers for various pathological and physiological conditions. The objective of this study was to identify miRNAs that are associated with egg qualities in rainbow trout using post-ovulatory aged eggs.

RESULTS

Egg samples from females on day 1, day 7, and day 14 post-ovulation (D1PO, D7PO and D14PO), which had the fertilization rates of 91.8%, 73.4% and less than 50%, respectively, were collected and small RNAs isolated from these samples were subjected to deep sequencing using the Illumina platform. The massive sequencing produced 27,342,477, 26,910,438 and 29,185,371 reads from the libraries of D1PO, D7PO and D14PO eggs, respectively. A three-way comparison of the miRNAs indicated that the egg samples shared 392 known and 236 novel miRNAs, and a total of 414, 481, and 470 known and 243, 298, and 296 novel miRNAs were identified from D1PO, D7PO and D14PO eggs, respectively. Four known miRNAs (omy-miR-193b-3p, omy-miR-203c-3p, omy-miR-499-5p and omy-miR-7550-3p) and two novel miRNAs (omy-miR-nov-95-5p and omy-miR-nov-112-5p) showed significantly higher expression in D1PO eggs relative to D14PO eggs as revealed by both deep sequencing and real time quantitative PCR analysis. GO analysis of the predicted target genes of these differentially expressed miRNAs revealed significantly enriched GO terms that are related to stress response, cell death, DNA damage, ATP generation, signal transduction and transcription regulation.

CONCLUSIONS

Results indicate that post-ovulatory ageing affects miRNA expression profiles in rainbow trout eggs, which can in turn impact egg quality. Further characterization of the differentially expressed miRNAs and their target genes may provide valuable information on the role of these miRNAs in controlling egg quality, and ultimately lead to the development of biomarkers for prediction of egg quality in rainbow trout.

Immunoglobulin isotypes in Atlantic salmon, Salmo salar

There are three major immunoglobulin (Ig) isotypes in salmonid fish: IgM, IgD and IgT, defined by the heavy chains μ, δ and τ, respectively. As a result of whole genome duplication in the ancestor of the salmonid fish family, Atlantic salmon (Salmo salar) possess two highly similar Ig heavy chain gene complexes (A and B), comprising two μ genes, two δ genes, three intact τ genes and five τ pseudogenes. The μA and μB genes correspond to two distinct sub-populations of serum IgM. The IgM-B sub-variant has a characteristic extra cysteine near the C-terminal part of the heavy chain and exhibits a higher degree of polymer disulfide cross-linking compared to IgM-A. The IgM-B:IgM-A ratio in serum is typically 60:40, but skewed ratios are also observed. The IgT isotype appears to be specialized to mucosal immune responses in salmonid fish. The concentration of IgT in serum is 100 to 1000 times lower than IgM. Secreted forms of IgD have been detected in rainbow trout, but not yet in Atlantic salmon.

The pH sensitivity of Aqp0 channels in tetraploid and diploid teleosts

Water homeostasis and the structural integrity of the vertebrate lens is partially mediated by AQP0 channels. Emerging evidence indicates that external pH may be involved in channel gating. Here we show that a tetraploid teleost, the Atlantic salmon, retains 4 aqp0 genes (aqp0a1, -0a2, -0b1, and -0b2), which are highly, but not exclusively, expressed in the lens. Functional characterization reveals that, although each paralog permeates water efficiently, the permeability is respectively shifted to the neutral, alkaline, or acidic pH in Aqp0a1, -0a2, and -0b1, whereas that of Aqp0b2 is not regulated by external pH. Mutagenesis studies demonstrate that Ser(38), His(39), and His(40) residues in the extracellular transmembrane domain of α-helix 2 facing the water pore are critical for the pH modulation of water transport. To validate these findings, we show that both zebrafish Aqp0a and -0b are functional water channels with respective pH sensitivities toward alkaline or acid pH ranges and that an N-terminal allelic variant (Ser(19)) of Aqp0b exists that abolishes water transport in Xenopus laevis oocytes. The data suggest that the alkaline pH sensitivity is a conserved trait in teleost Aqp0 a-type channels, whereas mammalian AQP0 and some teleost Aqp0 b-type channels display an acidic pH permeation preference.

Gonad specific genes in Atlantic salmon (Salmon salar L.): characterization of tdrd7-2, dazl-2, piwil1 and tdrd1 genes

Atlantic salmon is a commercially important species. Understanding key processes in their life history, such as germ cell development is essential for further improvements within salmon farming. Since salmonids have undergone an additional whole genome duplication compared to many other fish species, they possess more gene paralogues. Therefore, data on gene expression and function from other species may not apply for salmon. Our aim was to study the spatial and tissue specific expression of genes known from model species to be essential for germ cell development, to identify germ cell specific factors in salmon. Based on homology with other species, selected genes were predicted in the salmon genome assembly. Gene expression was measured by PCR in a variety of juvenile salmon tissues. For genes expressed exclusively in gonads we measured the expression in the same tissues as well as in eggs, embryos and larvae by qPCR. Finally, we revealed the cellular localization of the gonad specific mRNAs by in situ hybridization (ISH). Several of the selected genes (tdrd7, cxcr4b and dazl), were found in more than one copy (indicated by a number following the gene name) in the salmon genome. Expression of tdrd7-2, dazl-2, piwil1 and tdrd1 was detected exclusively in the testis and ovary of juvenile salmon, and transcripts of tdrd7-2, dazl-2 and piwil1 were localized within male and female germ cells. While tdrd7-2, piwil1 and tdrd1 were expressed in unfertilized eggs and all embryo and larval stages measured, dazl-2 was expressed in unfertilized eggs and embryos until the onset of gastrulation. This study shows that several of the genes known from model species to be essential for germ cell development, display paralogues in salmon with dissimilar and similar expression patterns in comparison to other species. Transcripts of tdrd7-2, dazl-2, piwil1 and tdrd1 are detected exclusively in gonads of juveniles and are found among maternal RNA of eggs and subsequent embryos. This information is valuable for further studies aiming at understanding salmon germ cell development.

Genetic evidence for ecological divergence in kokanee salmon

The evolution of locally adapted phenotypes among populations that experience divergent selective pressures is a central mechanism for generating and maintaining biodiversity. Recently, the advent of high-throughput DNA sequencing technology has provided tools for investigating the genetic basis of this process in natural populations of nonmodel organisms. Kokanee, the freshwater form of sockeye salmon (Oncorhynchus nerka), occurs as two reproductive ecotypes, which differ in spawning habitat (tributaries vs. shorelines); however, outside of the spawning season the two ecotypes co-occur in many lakes and lack diagnostic morphological characteristics. We used restriction site-associated DNA (RAD) sequencing to identify 6145 SNPs and genotype kokanee from multiple spawning sites in Okanagan Lake (British Columbia, Canada). Outlier tests revealed 18 loci putatively under divergent selection between ecotypes, all of which exhibited temporally stable allele frequencies within ecotypes. Six outliers were annotated to sequences in the NCBI database, two of which matched genes associated with early development. There was no evidence for neutral genetic differentiation; however, outlier loci demonstrated significant structure with respect to ecotype and had high assignment accuracy in mixed composition simulations. The absence of neutral structure combined with a small number of highly divergent outlier loci is consistent with theoretical predictions for the early stages of ecological divergence. These outlier loci were then applied to a realistic fisheries scenario in which additional RAD sequencing was used to genotype kokanee collected by trawl in Okanagan Lake, providing preliminary evidence that this approach may be an effective tool for conservation and management.

Applications in the search for genomic selection signatures in fish

Selection signatures are genomic regions harboring DNA sequences functionally involved in the genetic variation of traits subject to selection. Selection signatures have been intensively studied in recent years because of their relevance to evolutionary biology and their potential association with genes that control phenotypes of interest in wild and domestic populations. Selection signature research in fish has been confined to a smaller scale, due in part to the relatively recent domestication of fish species and limited genomic resources such as molecular markers, genetic mapping, DNA sequences, and reference genomes. However, recent genomic technology advances are paving the way for more studies that may contribute to the knowledge of genomic regions underlying phenotypes of biological and productive interest in fish.

Whole-body transcriptome of selectively bred, resistant-, control-, and susceptible-line rainbow trout following experimental challenge with Flavobacterium psychrophilum

Genetic improvement for enhanced disease resistance in fish is an increasingly utilized approach to mitigate endemic infectious disease in aquaculture. In domesticated salmonid populations, large phenotypic variation in disease resistance has been identified but the genetic basis for altered responsiveness remains unclear. We previously reported three generations of selection and phenotypic validation of a bacterial cold water disease (BCWD) resistant line of rainbow trout, designated ARS-Fp-R. This line has higher survival after infection by either standardized laboratory challenge or natural challenge as compared to two reference lines, designated ARS-Fp-C (control) and ARS-Fp-S (susceptible). In this study, we utilized 1.1 g fry from the three genetic lines and performed RNA-seq to measure transcript abundance from the whole body of naive and Flavobacterium psychrophilum infected fish at day 1 (early time-point) and at day 5 post-challenge (onset of mortality). Sequences from 24 libraries were mapped onto the rainbow trout genome reference transcriptome of 46,585 predicted protein coding mRNAs that included 2633 putative immune-relevant gene transcripts. A total of 1884 genes (4.0% genome) exhibited differential transcript abundance between infected and mock-challenged fish (FDR < 0.05) that included chemokines, complement components, tnf receptor superfamily members, interleukins, nod-like receptor family members, and genes involved in metabolism and wound healing. The largest number of differentially expressed genes occurred on day 5 post-infection between naive and challenged ARS-Fp-S line fish correlating with high bacterial load. After excluding the effect of infection, we identified 21 differentially expressed genes between the three genetic lines. In summary, these data indicate global transcriptome differences between genetic lines of naive animals as well as differentially regulated transcriptional responses to infection.

The Genome 10K Project: a way forward

The Genome 10K Project was established in 2009 by a consortium of biologists and genome scientists determined to facilitate the sequencing and analysis of the complete genomes of 10,000 vertebrate species. Since then the number of selected and initiated species has risen from ∼26 to 277 sequenced or ongoing with funding, an approximately tenfold increase in five years. Here we summarize the advances and commitments that have occurred by mid-2014 and outline the achievements and present challenges of reaching the 10,000-species goal. We summarize the status of known vertebrate genome projects, recommend standards for pronouncing a genome as sequenced or completed, and provide our present and future vision of the landscape of Genome 10K. The endeavor is ambitious, bold, expensive, and uncertain, but together the Genome 10K Consortium of Scientists and the worldwide genomics community are moving toward their goal of delivering to the coming generation the gift of genome empowerment for many vertebrate species.

Basal Gnathostomes provide unique insights into the evolution of vitamin B12 binders

The uptake and transport of vitamin B12 (cobalamin; Cbl) in mammals involves a refined system with three evolutionarily related transporters: transcobalamin 1 (Tcn1), transcobalamin 2 (Tcn2), and the gastric intrinsic factor (Gif). Teleosts have a single documented binder with intermediate features to the human counterparts. Consequently, it has been proposed that the expansion of Cbl binders occurred after the separation of Actinopterygians. Here, we demonstrate that the diversification of this gene family took place earlier in gnathostome ancestry. Our data indicates the presence of single copy orthologs of the Sarcopterygii/Tetrapoda duplicates Tcn1 and Gif, and Tcn2, in Chondrichthyes. In addition, a highly divergent Cbl binder was found in the Elasmobranchii. We unveil a complex scenario forged by genome, tandem duplications and lineage-specific gene loss. Our findings suggest that from an ancestral transporter, exhibiting large spectrum and high affinity binding, highly specific Cbl transporters emerged through gene duplication and mutations at the binding pocket.

High-throughput sequencing of microRNA transcriptome and expression assay in the sturgeon, Acipenser schrenckii

Sturgeons are considered as living fossils and have very high evolutionary, economical and conservation values. The multiploidy of sturgeon that has been caused by chromosome duplication may lead to the emergence of new microRNAs (miRNAs) involved in the ploidy and physiological processes. In the present study, we performed the first sturgeon miRNAs analysis by RNA-seq high-throughput sequencing combined with expression assay of microarray and real-time PCR, and aimed to discover the sturgeon-specific miRNAs, confirm the expressed pattern of miRNAs and illustrate the potential role of miRNAs-targets on sturgeon biological processes. A total of 103 miRNAs were identified, including 58 miRNAs with strongly detected signals (signal >500 and P≤0.01), which were detected by microarray. Real-time PCR assay supported the expression pattern obtained by microarray. Moreover, co-expression of 21 miRNAs in all five tissues and tissue-specific expression of 16 miRNAs implied the crucial and particular function of them in sturgeon physiological processes. Target gene prediction, especially the enriched functional gene groups (369 GO terms) and pathways (37 KEGG) regulated by 58 miRNAs (P<0.05), illustrated the interaction of miRNAs and putative mRNAs, and also the potential mechanism involved in these biological processes. Our new findings of sturgeon miRNAs expand the public database of transcriptome information for this species, contribute to our understanding of sturgeon biology, and also provide invaluable data that may be applied in sturgeon breeding.

Evolutionary history and epigenetic regulation of the three paralogous pax7 genes in rainbow trout

The extraordinary muscle growth potential of teleost fish, particular those of the Salmoninae clade, elicits questions about the regulation of the relatively highly conserved transcription factors of the myogenic program. The pseudotetraploid nature of the salmonid genome adds another layer of regulatory complexity that must be reconciled with epigenetic data to improve our understanding of the achievement of lifelong muscle growth in these fish. We identify three paralogous pax7 genes (pax7a1, pax7a2 and pax7b) in the rainbow trout genome. During in vitro myogenesis, pax7a1 transcripts remain stable, whereas pax7a2 and pax7b mRNAs increase in abundance, similarly to myogenin mRNAs but in contrast to the expression pattern of the mammalian ortholog. We also profile the distribution of repressive H3K27me3 and H3K9me3 and permissive H3K4me3 marks during in vitro myogenesis across these loci and find that pax7a2 expression is associated with decreased H3K27 trimethylation, whereas pax7b expression is correlated with decreased H3K9me3 and H3K27me3. These data link the unique differential expression of pax7 paralogs with epigenetic histone modifications in a vertebrate species displaying growth divergent from that of mammals and highlight an important divergence in the regulatory mechanisms of pax7 expression among vertebrates. The system described here provides a more comprehensive picture of the combinatorial control mechanisms orchestrating skeletal muscle growth in a salmonid, leading to a better understanding of myogenesis in this species and across Vertebrata more generally.

The lineage-specific evolution of aquaporin gene clusters facilitated tetrapod terrestrial adaptation

A major physiological barrier for aquatic organisms adapting to terrestrial life is dessication in the aerial environment. This barrier was nevertheless overcome by the Devonian ancestors of extant Tetrapoda, but the origin of specific molecular mechanisms that solved this water problem remains largely unknown. Here we show that an ancient aquaporin gene cluster evolved specifically in the sarcopterygian lineage, and subsequently diverged into paralogous forms of AQP2, -5, or -6 to mediate water conservation in extant Tetrapoda. To determine the origin of these apomorphic genomic traits, we combined aquaporin sequencing from jawless and jawed vertebrates with broad taxon assembly of >2,000 transcripts amongst 131 deuterostome genomes and developed a model based upon Bayesian inference that traces their convergent roots to stem subfamilies in basal Metazoa and Prokaryota. This approach uncovered an unexpected diversity of aquaporins in every lineage investigated, and revealed that the vertebrate superfamily consists of 17 classes of aquaporins (Aqp0 - Aqp16). The oldest orthologs associated with water conservation in modern Tetrapoda are traced to a cluster of three aqp2-like genes in Actinistia that likely arose >500 Ma through duplication of an aqp0-like gene present in a jawless ancestor. In sea lamprey, we show that aqp0 first arose in a protocluster comprised of a novel aqp14 paralog and a fused aqp01 gene. To corroborate these findings, we conducted phylogenetic analyses of five syntenic nuclear receptor subfamilies, which, together with observations of extensive genome rearrangements, support the coincident loss of ancestral aqp2-like orthologs in Actinopterygii. We thus conclude that the divergence of sarcopterygian-specific aquaporin gene clusters was permissive for the evolution of water conservation mechanisms that facilitated tetrapod terrestrial adaptation.

Genetics and genomics of disease resistance in salmonid species

Infectious and parasitic diseases generate large economic losses in salmon farming. A feasible and sustainable alternative to prevent disease outbreaks may be represented by genetic improvement for disease resistance. To include disease resistance into the breeding goal, prior knowledge of the levels of genetic variation for these traits is required. Furthermore, the information from the genetic architecture and molecular factors involved in resistance against diseases may be used to accelerate the genetic progress for these traits. In this regard, marker assisted selection and genomic selection are approaches which incorporate molecular information to increase the accuracy when predicting the genetic merit of selection candidates. In this article we review and discuss key aspects related to disease resistance in salmonid species, from both a genetic and genomic perspective, with emphasis in the applicability of disease resistance traits into breeding programs in salmonids.

Association mapping of disease resistance traits in rainbow trout using restriction site associated DNA sequencing

Recent advances in genotyping-by-sequencing have enabled genome-wide association studies in nonmodel species including those in aquaculture programs. As with other aquaculture species, rainbow trout and steelhead (Oncorhynchus mykiss) are susceptible to disease and outbreaks can lead to significant losses. Fish culturists have therefore been pursuing strategies to prevent losses to common pathogens such as Flavobacterium psychrophilum (the etiological agent for bacterial cold water disease [CWD]) and infectious hematopoietic necrosis virus (IHNV) by adjusting feed formulations, vaccine development, and selective breeding. However, discovery of genetic markers linked to disease resistance offers the potential to use marker-assisted selection to increase resistance and reduce outbreaks. For this study we sampled juvenile fish from 40 families from 2-yr classes that either survived or died after controlled exposure to either CWD or IHNV. Restriction site−associated DNA sequencing produced 4661 polymorphic single-nucleotide polymorphism loci after strict filtering. Genotypes from individual survivors and mortalities were then used to test for association between disease resistance and genotype at each locus using the program TASSEL. After we accounted for kinship and stratification of the samples, tests revealed 12 single-nucleotide polymorphism markers that were highly associated with resistance to CWD and 19 markers associated with resistance to IHNV. These markers are candidates for further investigation and are expected to be useful for marker assisted selection in future broodstock selection for various aquaculture programs.

The development and characterization of a 57K single nucleotide polymorphism array for rainbow trout

In this study, we describe the development and characterization of the first high-density single nucleotide polymorphism (SNP) genotyping array for rainbow trout. The SNP array is publically available from a commercial vendor (Affymetrix). The SNP genotyping quality was high, and validation rate was close to 90%. This is comparable to other farm animals and is much higher than previous smaller scale SNP validation studies in rainbow trout. High quality and integrity of the genotypes are evident from sample reproducibility and from nearly 100% agreement in genotyping results from other methods. The array is very useful for rainbow trout aquaculture populations with more than 40 900 polymorphic markers per population. For wild populations that were confounded by a smaller sample size, the number of polymorphic markers was between 10 577 and 24 330. Comparison between genotypes from individual populations suggests good potential for identifying candidate markers for populations' traceability. Linkage analysis and mapping of the SNPs to the reference genome assembly provide strong evidence for a wide distribution throughout the genome with good representation in all 29 chromosomes. A total of 68% of the genome scaffolds and contigs were anchored through linkage analysis using the SNP array genotypes, including ~20% of the genome assembly that has not been previously anchored to chromosomes.

Characterization of the rainbow trout spleen transcriptome and identification of immune-related genes

Resistance against diseases affects profitability of rainbow trout. Limited information is available about functions and mechanisms of teleost immune pathways. Immunogenomics provides powerful tools to determine disease resistance genes/gene pathways and develop genetic markers for genomic selection. RNA-Seq sequencing of the rainbow trout spleen yielded 93,532,200 reads (100 bp). High quality reads were assembled into 43,047 contigs. 26,333 (61.17%) of the contigs had hits to the NR protein database and 7024 (16.32%) had hits to the KEGG database. Gene ontology showed significant percentages of transcripts assigned to binding (51%), signaling (7%), response to stimuli (9%) and receptor activity (4%) suggesting existence of many immune-related genes. KEGG annotation revealed 2825 sequences belonging to "organismal systems" with the highest number of sequences, 842 (29.81%), assigned to immune system. A number of sequences were identified for the first time in rainbow trout belonging to Toll-like receptor signaling (35), B cell receptor signaling pathway (44), T cell receptor signaling pathway (56), chemokine signaling pathway (73), Fc gamma R-mediated phagocytosis (52), leukocyte transendothelial migration (60) and NK cell mediated cytotoxicity (42). In addition, 51 transcripts were identified as spleen-specific genes. The list includes 277 full-length cDNAs. The presence of a large number of immune-related genes and pathways similar to other vertebrates suggests that innate and adaptive immunity in fish are conserved. This study provides deep-sequence data of rainbow trout spleen transcriptome and identifies many new immune-related genes and full-length cDNAs. This data will help identify allelic variations suitable for genomic selection and genetic manipulation in aquaculture.

Temporally isolated lineages of pink salmon reveal unique signatures of selection on distinct pools of standing genetic variation

A species' genetic diversity bears the marks of evolutionary processes that have occurred throughout its history. However, robust detection of selection in wild populations is difficult and often impeded by lack of replicate tests. Here, we investigate selection in pink salmon (Oncorhynchus gorbuscha) using genome scans coupled with inference from a haploid-assisted linkage map. Pink salmon have a strict 2-year semelparous life history which has resulted in temporally isolated (allochronic) lineages that remain sympatric through sharing of spawning habitats in alternate years. The lineages differ in a range of adaptive traits, suggesting different genetic backgrounds. We used genotyping by sequencing of haploids to generate a high-density linkage map with 7035 loci and screened an existing panel of 8036 loci for signatures of selection. The linkage map enabled identification of novel genomic regions displaying signatures of parallel selection shared between lineages. Furthermore, 24 loci demonstrated divergent selection and differences in genetic diversity between lineages, suggesting that adaptation in the 2 lineages has arisen from different pools of standing genetic variation. Findings have implications for understanding asynchronous population abundances as well as predicting future ecosystem impacts from lineage-specific responses to climate change.

Enigmatic orthology relationships between Hox clusters of the African butterfly fish and other teleosts following ancient whole-genome duplication

Numerous ancient whole-genome duplications (WGD) have occurred during eukaryote evolution. In vertebrates, duplicated developmental genes and their functional divergence have had important consequences for morphological evolution. Although two vertebrate WGD events (1R/2R) occurred over 525 Ma, we have focused on the more recent 3R or TGD (teleost genome duplication) event which occurred approximately 350 Ma in a common ancestor of over 26,000 species of teleost fishes. Through a combination of whole genome and bacterial artificial chromosome clone sequencing we characterized all Hox gene clusters of Pantodon buchholzi, a member of the early branching teleost subdivision Osteoglossomorpha. We find 45 Hox genes organized in only five clusters indicating that Pantodon has suffered more Hox cluster loss than other known species. Despite strong evidence for homology of the five Pantodon clusters to the four canonical pre-TGD vertebrate clusters (one HoxA, two HoxB, one HoxC, and one HoxD), we were unable to confidently resolve 1:1 orthology relationships between four of the Pantodon clusters and the eight post-TGD clusters of other teleosts. Phylogenetic analysis revealed that many Pantodon genes segregate outside the conventional "a" and "b" post-TGD orthology groups, that extensive topological incongruence exists between genes physically linked on a single cluster, and that signal divergence causes ambivalence in assigning 1:1 orthology in concatenated Hox cluster analyses. Out of several possible explanations for this phenomenon we favor a model which keeps with the prevailing view of a single TGD prior to teleost radiation, but which also considers the timing of diploidization after duplication, relative to speciation events. We suggest that although the duplicated hoxa clusters diploidized prior to divergence of osteoglossomorphs, the duplicated hoxb, hoxc, and hoxd clusters concluded diploidization independently in osteoglossomorphs and other teleosts. We use the term "tetralogy" to describe the homology relationship which exists between duplicated sequences which originate through a shared WGD, but which diploidize into distinct paralogs from a common allelic pool independently in two lineages following speciation.

Targeted mutagenesis in Atlantic salmon (Salmo salar L.) using the CRISPR/Cas9 system induces complete knockout individuals in the F0 generation

Understanding the biological function behind key proteins is of great concern in Atlantic salmon, both due to a high commercial importance and an interesting life history. Until recently, functional studies in salmonids appeared to be difficult. However, the recent discovery of targeted mutagenesis using the CRISPR/Cas9 (clustered regularly interspaced palindromic repeats/CRISPR-associated) system enables performing functional studies in Atlantic salmon to a great extent. We used the CRISPR/Cas9 system to target two genes involved in pigmentation, tyrosinase (tyr) and solute carrier family 45, member 2 (slc45a2). Embryos were assayed for mutation rates at the 17 somite stage, where 40 and 22% of all injected embryos showed a high degree of mutation induction for slc45a2 and tyr, respectively. At hatching this mutation frequency was also visible for both targeted genes, displaying a graded phenotype ranging from complete lack of pigmentation to partial loss and normal pigmentation. CRISPRslc45a2/Cas9 injected embryos showing a complete lack of pigmentation or just a few spots of pigments also lacked wild type sequences when assaying more than 80 (slc45a2) sequence clones from whole embryos. This indicates that CRISPR/Cas9 can induce double-allelic knockout in the F0 generation. However, types and frequency of indels might affect the phenotype. Therefore, the variation of indels was assayed in the graded pigmentation phenotypes produced by CRISPR/Cas9-slc45a2. The results show a tendency for fewer types of indels formed in juveniles completely lacking pigmentation compared to juveniles displaying partial pigmentation. Another interesting observation was a high degree of the same indel type in different juveniles. This study shows for the first time successful use of the CRISPR/Cas9 technology in a marine cold water species. Targeted double-allelic mutations were obtained and, though the level of mosaicism has to be considered, we demonstrate that F0 fish can be used for functional studies in Atlantic salmon.

Origins and functional diversification of salinity-responsive Na(+) , K(+) ATPase α1 paralogs in salmonids

The Salmoniform whole-genome duplication is hypothesized to have facilitated the evolution of anadromy, but little is known about the contribution of paralogs from this event to the physiological performance traits required for anadromy, such as salinity tolerance. Here, we determined when two candidate, salinity-responsive paralogs of the Na(+) , K(+) ATPase α subunit (α1a and α1b) evolved and studied their evolutionary trajectories and tissue-specific expression patterns. We found that these paralogs arose during a small-scale duplication event prior to the Salmoniform, but after the teleost, whole-genome duplication. The 'freshwater paralog' (α1a) is primarily expressed in the gills of Salmoniformes and an unduplicated freshwater sister species (Esox lucius) and experienced positive selection in the freshwater ancestor of Salmoniformes and Esociformes. Contrary to our predictions, the 'saltwater paralog' (α1b), which is more widely expressed than α1a, did not experience positive selection during the evolution of anadromy in the Coregoninae and Salmonine. To determine whether parallel mutations in Na(+) , K(+) ATPase α1 may contribute to salinity tolerance in other fishes, we studied independently evolved salinity-responsive Na(+) , K(+) ATPase α1 paralogs in Anabas testudineus and Oreochromis mossambicus. We found that a quarter of the mutations occurring between salmonid α1a and α1b in functionally important sites also evolved in parallel in at least one of these species. Together, these data argue that paralogs contributing to salinity tolerance evolved prior to the Salmoniform whole-genome duplication and that strong selection and/or functional constraints have led to parallel evolution in salinity-responsive Na(+) , K(+) ATPase α1 paralogs in fishes.

A new model army: Emerging fish models to study the genomics of vertebrate Evo-Devo

Many fields of biology--including vertebrate Evo-Devo research--are facing an explosion of genomic and transcriptomic sequence information and a multitude of fish species are now swimming in this "genomic tsunami." Here, we first give an overview of recent developments in sequencing fish genomes and transcriptomes that identify properties of fish genomes requiring particular attention and propose strategies to overcome common challenges in fish genomics. We suggest that the generation of chromosome-level genome assemblies--for which we introduce the term "chromonome"--should be a key component of genomic investigations in fish because they enable large-scale conserved synteny analyses that inform orthology detection, a process critical for connectivity of genomes. Orthology calls in vertebrates, especially in teleost fish, are complicated by divergent evolution of gene repertoires and functions following two rounds of genome duplication in the ancestor of vertebrates and a third round at the base of teleost fish. Second, using examples of spotted gar, basal teleosts, zebrafish-related cyprinids, cavefish, livebearers, icefish, and lobefin fish, we illustrate how next generation sequencing technologies liberate emerging fish systems from genomic ignorance and transform them into a new model army to answer longstanding questions on the genomic and developmental basis of their biodiversity. Finally, we discuss recent progress in the genetic toolbox for the major fish models for functional analysis, zebrafish, and medaka, that can be transferred to many other fish species to study in vivo the functional effect of evolutionary genomic change as Evo-Devo research enters the postgenomic era.

The genome and linkage map of the northern pike (Esox lucius): conserved synteny revealed between the salmonid sister group and the Neoteleostei

The northern pike is the most frequently studied member of the Esociformes, the closest order to the diverse and economically important Salmoniformes. The ancestor of all salmonids purportedly experienced a whole-genome duplication (WGD) event, making salmonid species ideal for studying the early impacts of genome duplication while complicating their use in wider analyses of teleost evolution. Studies suggest that the Esociformes diverged from the salmonid lineage prior to the WGD, supporting the use of northern pike as a pre-duplication outgroup. Here we present the first genome assembly, reference transcriptome and linkage map for northern pike, and evaluate the suitability of this species to provide a representative pre-duplication genome for future studies of salmonid and teleost evolution. The northern pike genome sequence is composed of 94,267 contigs (N50 = 16,909 bp) contained in 5,688 scaffolds (N50 = 700,535 bp); the total scaffolded genome size is 878 million bases. Multiple lines of evidence suggest that over 96% of the protein-coding genome is present in the genome assembly. The reference transcriptome was constructed from 13 tissues and contains 38,696 transcripts, which are accompanied by normalized expression data in all tissues. Gene-prediction analysis produced a total of 19,601 northern pike-specific gene models. The first-generation linkage map identifies 25 linkage groups, in agreement with northern pike's diploid karyotype of 2N = 50, and facilitates the placement of 46% of assembled bases onto linkage groups. Analyses reveal a high degree of conserved synteny between northern pike and other model teleost genomes. While conservation of gene order is limited to smaller syntenic blocks, the wider conservation of genome organization implies the northern pike exhibits a suitable approximation of a non-duplicated Protacanthopterygiian genome. This dataset will facilitate future studies of esocid biology and empower ongoing examinations of the Atlantic salmon and rainbow trout genomes by facilitating their comparison with other major teleost groups.

Comparative mapping between Coho Salmon (Oncorhynchus kisutch) and three other salmonids suggests a role for chromosomal rearrangements in the retention of duplicated regions following a whole genome duplication event

Whole genome duplication has been implicated in evolutionary innovation and rapid diversification. In salmonid fishes, however, whole genome duplication significantly pre-dates major transitions across the family, and re-diploidization has been a gradual process between genomes that have remained essentially collinear. Nevertheless, pairs of duplicated chromosome arms have diverged at different rates from each other, suggesting that the retention of duplicated regions through occasional pairing between homeologous chromosomes may have played an evolutionary role across species pairs. Extensive chromosomal arm rearrangements have been a key mechanism involved in re-dipliodization of the salmonid genome; therefore, we investigated their influence on degree of differentiation between homeologs across salmon species. We derived a linkage map for coho salmon and performed comparative mapping across syntenic arms within the genus Oncorhynchus, and with the genus Salmo, to determine the phylogenetic relationship between chromosome arrangements and the retention of undifferentiated duplicated regions. A 6596.7 cM female coho salmon map, comprising 30 linkage groups with 7415 and 1266 nonduplicated and duplicated loci, respectively, revealed uneven distribution of duplicated loci along and between chromosome arms. These duplicated regions were conserved across syntenic arms across Oncorhynchus species and were identified in metacentric chromosomes likely formed ancestrally to the divergence of Oncorhynchus from Salmo. These findings support previous studies in which observed pairings involved at least one metacentric chromosome. Re-diploidization in salmon may have been prevented or retarded by the formation of metacentric chromosomes after the whole genome duplication event and may explain lineage-specific innovations in salmon species if functional genes are found in these regions.

Extensive local gene duplication and functional divergence among paralogs in Atlantic salmon

Many organisms can generate alternative phenotypes from the same genome, enabling individuals to exploit diverse and variable environments. A prevailing hypothesis is that such adaptation has been favored by gene duplication events, which generate redundant genomic material that may evolve divergent functions. Vertebrate examples of recent whole-genome duplications are sparse although one example is the salmonids, which have undergone a whole-genome duplication event within the last 100 Myr. The life-cycle of the Atlantic salmon, Salmo salar, depends on the ability to produce alternating phenotypes from the same genome, to facilitate migration and maintain its anadromous life history. Here, we investigate the hypothesis that genome-wide and local gene duplication events have contributed to the salmonid adaptation. We used high-throughput sequencing to characterize the transcriptomes of three key organs involved in regulating migration in S. salar: Brain, pituitary, and olfactory epithelium. We identified over 10,000 undescribed S. salar sequences and designed an analytic workflow to distinguish between paralogs originating from local gene duplication events or from whole-genome duplication events. These data reveal that substantial local gene duplications took place shortly after the whole-genome duplication event. Many of the identified paralog pairs have either diverged in function or become noncoding. Future functional genomics studies will reveal to what extent this rich source of divergence in genetic sequence is likely to have facilitated the evolution of extreme phenotypic plasticity required for an anadromous life-cycle.

Sex hormone-binding globulins characterization and gonadal gene expression during sex differentiation in the rainbow trout, Oncorhynchus mykiss

Sex hormone-binding globulin (SHBG) binds androgens and estrogens in the blood of many vertebrates, including teleost fish. In mammals, SHBG is synthetized in the liver and secreted into the blood. In fish, shbga also exhibits a hepatic expression. In salmonids, in which the gene has been duplicated, the recently discovered shbgb gene exhibits a predominantly ovarian expression. The present work aimed at gaining new insight into shbgb gene structure and expression during gonadal sex differentiation, a steroid-sensitive process, and Shbgb protein structure and binding characteristics; specifically, rainbow trout (Oncorhynchus mykiss) shbgb was analyzed. shbgb structure was analyzed in silico while expression was characterized during gonadal sex differentiation using all-male and all-female populations. We observed that shbgb gene and cognate-protein structures are similar to homologs previously described in zebrafish and mammals. The shbgb gene is predominantly expressed in differentiating female gonads, with increased expression around the end of ovarian differentiation. In the ovary, shbgb mRNA was detected in a subset of somatic cells surrounding the ovarian lamellae. Furthermore, Shbgb binds steroids with a higher selectivity than Shbga, exhibiting a higher affinity for estradiol compared to Shbga. In conclusion, Shbgb binding characteristics are clearly different from those of Shbga. Shbgb is expressed in the differentiating ovary during a period when the synthesis and action of testosterone and estradiol must be tightly regulated. This strongly suggests that Shbgb participates in the regulation of steroid metabolism and/or mediation, that is, needed during early gonadal development in rainbow trout.

Structural and functional divergence of gonadotropin-inhibitory hormone from jawless fish to mammals

Gonadotropin-inhibitory hormone (GnIH) was discovered as a novel hypothalamic peptide that inhibits gonadotropin release in the quail. The presence of GnIH-homologous peptides and its receptors (GnIHRs) have been demonstrated in various vertebrate species including teleosts, suggesting that the GnIH-GnIHR family is evolutionarily conserved. In avian and mammalian brain, GnIH neurons are localized in the hypothalamic nuclei and their neural projections are widely distributed. GnIH acts on the pituitary and gonadotropin-releasing hormone neurons to inhibit reproductive functions by decreasing gonadotropin release and synthesis. In addition, GnIH-GnIHR signaling is regulated by various factors, such as environmental cues and stress. However, the function of fish GnIH orthologs remains inconclusive because the physiological properties of fish GnIH peptides are debatable. This review summarizes the current research progress in GnIH-GnIHR signaling and their physiological functions in vertebrates with special emphasis on non-mammalian vertebrate species.