Salmonid-specific Whole Genome Duplication Research Articles

Comprehensive Analysis of YTH Domain-Containing Genes, Encoding m6A Reader and Their Response to Temperature Stresses and Yersinia ruckeri Infection in Rainbow Trout (Oncorhynchus mykiss)

YTH domain-containing genes are important readers of N6-methyladenosine (m6A) modifications with ability to directly affect the fates of distinct RNAs in organisms. Despite their importance, little is known about YTH domain-containing genes in teleosts until now. In the present study, a total of 10 YTH domain-containing genes have been systematically identified and functionally characterized in rainbow trout (Oncorhynchus mykiss). According to the phylogenetic tree, gene structure and syntenic analysis, these YTH domain-containing genes could be classified into three evolutionary subclades, including YTHDF, YTHDC1 and YTHDC2. Of them, the copy number of OmDF1, OmDF2, OmDF3, and OmDC1 were duplicated or even triplicated in rainbow trout due to the salmonid-specific whole-genome duplication event. The three-dimensional protein structure analysis revealed that there were similar structures and the same amino acid residues that were associated with cage formation between humans and rainbow trout, implying their similar manners in binding to m6A modification. Additionally, the results of qPCR experiment indicated that the expression patterns of a few YTH domain-containing genes, especially OmDF1b, OmDF3a and OmDF3b, were significantly different in liver tissue of rainbow trout under four different temperatures (7 °C, 11 °C, 15 °C, and 19 °C). The expression levels of OmDF1a, OmDF1b and OmDC1a were obviously repressed in spleen tissue of rainbow trout at 24 h after Yersinia ruckeri infection, while increased expression was detected in OmDF3b. This study provides a systemic overview of YTH domain-containing genes in rainbow trout and reveals their biological roles in responses to temperature stress and bacterial infection.

A homology guide for Pacific salmon genus Oncorhynchus resolves patterns of ohnolog retention, resolution and local adaptation following the salmonid-specific whole-genome duplication event.

Salmonid fishes have emerged as a tractable model to study whole-genome duplications (WGDs) as this group has undergone four rounds of WGDs. While most of the salmonid genome has returned to a diploid state, a significant proportion of genes are maintained as duplicates and are referred to as ohnologs. The fact that much of the modern salmonid gene repertoire is comprised of ohnologs, while other genes have returned to their singleton state creates complications for genetic studies by obscuring homology relationships. The difficulty this creates is particularly prominent in Pacific salmonids belonging to genus Oncorhynchus who are the focus of intense genetics-based conservation and management efforts owing to the important ecological and cultural roles these fish play. To address this gap, we generated a homology guide for six species of Oncorhynchus with available genomes and used this guide to describe patterns of ohnolog retention and resolution. Overall, we find that ohnologs comprise approximately half of each species modern gene repertoires, which are functionally enriched for genes involved in DNA binding, while the less numerous singleton genes are heavily enriched in dosage-sensitive processes such as mitochondrial metabolism. Additionally, by reanalyzing published expression data from locally adapted strains of O. mykiss, we show that numerous ohnologs exhibit adaptive expression profiles; however, ohnologs are not more likely to display adaptive signatures than either paralogs or singletons. Finally, we demonstrate the utility of our homology guide by investigating the evolutionary relationship among genes highlighted as playing a role in salmonid life-history traits or gene editing targets.

Differential retention of duplicated retinoid-binding protein (crabp &rbp) genes in the rainbow trout genome after two whole genome duplications and their responses to dietary canola oil

The complexity of retinoid intracellular signaling by lipid-binding proteins is further confounded in salmonids, owing to a whole genome duplication (WGD), followed by a second salmonid-specific WGD. As such, many more cellular retinoic acid-binding protein (crabp)/retinol-binding protein (rbp) genes may exist in salmonids, compared to other vertebrate species. In this study, genomic database searches and phylogenetic analyses identified seven crabp and eight rbp genes in the rainbow trout. A parsimonious explanation for the number of crabp/rbp genes in the rainbow trout genome is duplicated genes were lost from the genome by nonfunctionalization following the salmonid-WGD, and the extant duplicate crabp/rbp genes of the much earlier teleost-WGD have been retained by either sub- or neofunctionalization of tissue-specific regulatory mechanisms. Nutritional regulation of retinoid-binding protein genes in the liver was also determined in the fish fed diets containing various amounts of canola oil (CO) for 12 weeks along with hepatic fatty acid profile, Vitamin A (VitA) level, and growth parameters. The crabp/rbp transcripts were modulated markedly by experimental diets except for rbp2a.1/rbp2a.2 at the end of the trial. This result suggests that rbp2a.1/rbp2a.2 don't exhibit obvious divergence of transcriptional regulatory mechanisms that control their nutritional regulation by dietary CO. Growth parameters didn't differ markedly among dietary treatments. An understanding of VitA homeostasis, especially the regulatory effects of retinoids on transcriptional control of target genes, may have considerable importance to basic fish development, physiology, and aquaculture and this study may lead to future research on molecular VitA metabolism especially in salmonids experienced salmonid-WGD.

Identification of another primordial CD80/86 molecule in rainbow trout: Insights into the origin and evolution of CD80 and CD86 in vertebrates

In mammals, the binding of distinct costimulatory ligands CD80 and CD86 to their receptors is essential for optimal T cell activation. Previous studies have shown that only a single CD80/86 gene exists in rainbow trout (termed rtCD80/86A) and other teleost fish, suggesting that CD80 and CD86 arose by gene duplication in the tetrapod branch, after the separation of fish and tetrapods. However, in this study, another CD80/86 gene has been cloned from rainbow trout, termed rtCD80/86B. The sequence identity between trout CD80/86 is significantly higher than that between CD80 and CD86 in tetrapods, indicating that CD80 and CD86 underwent divergent evolution in vertebrates, especially in tetrapods. Gene synteny analyses showed that the CD80 and CD86 genes are closely located in the same chromosome in tetrapods. However, CD80/86 genes are located in two distinct chromosomes in rainbow trout and Atlantic salmon, suggesting that salmonid CD80/86 genes arose by the salmonid-specific whole-genome duplication (WGD) event. Expression analysis showed that rtCD80/86A was more abundant and inducible than rtCD80/86B in various tissues, indicating the important role of rtCD80/86A in trout immunity. Interestingly, we found that head kidney B cells showed higher expression level of rtCD80/86A and rtCD80/86B when compared with the other leukocytes, suggesting a potential role for trout B cells as antigen-presenting cells (APCs).

Life-stage-associated remodelling of lipid metabolism regulation in Atlantic salmon.

Atlantic salmon migrates from rivers to sea to feed, grow and develop gonads before returning to spawn in freshwater. The transition to marine habitats is associated with dramatic changes in the environment, including water salinity, exposure to pathogens and shift in dietary lipid availability. Many changes in physiology and metabolism occur across this life-stage transition, but little is known about the molecular nature of these changes. Here, we use a long-term feeding experiment to study transcriptional regulation of lipid metabolism in Atlantic salmon gut and liver in both fresh- and saltwater. We find that lipid metabolism becomes significantly less plastic to differences in dietary lipid composition when salmon transitions to saltwater and experiences increased dietary lipid availability. Expression of genes in liver relating to lipogenesis and lipid transport decreases overall and becomes less responsive to diet, while genes for lipid uptake in gut become more highly expressed. Finally, analyses of evolutionary consequences of the salmonid-specific whole-genome duplication on lipid metabolism reveal several pathways with significantly different (p<.05) duplicate retention or duplicate regulatory conservation. We also find a limited number of cases where the whole-genome duplication has resulted in an increased gene dosage. In conclusion, we find variable and pathway-specific effects of the salmonid genome duplication on lipid metabolism genes. A clear life-stage-associated shift in lipid metabolism regulation is evident, and we hypothesize this to be, at least partly, driven by nondietary factors such as the preparatory remodelling of gene regulation and physiology prior to sea migration.

Evolution and Expression of Tissue Globins in Ray-Finned Fishes.

The globin gene family encodes oxygen-binding hemeproteins conserved across the major branches of multicellular life. The origins and evolutionary histories of complete globin repertoires have been established for many vertebrates, but there remain major knowledge gaps for ray-finned fish. Therefore, we used phylogenetic, comparative genomic and gene expression analyses to discover and characterize canonical “non-blood” globin family members (i.e., myoglobin, cytoglobin, neuroglobin, globin-X, and globin-Y) across multiple ray-finned fish lineages, revealing novel gene duplicates (paralogs) conserved from whole genome duplication (WGD) and small-scale duplication events. Our key findings were that: (1) globin-X paralogs in teleosts have been retained from the teleost-specific WGD, (2) functional paralogs of cytoglobin, neuroglobin, and globin-X, but not myoglobin, have been conserved from the salmonid-specific WGD, (3) triplicate lineage-specific myoglobin paralogs are conserved in arowanas (Osteoglossiformes), which arose by tandem duplication and diverged under positive selection, (4) globin-Y is retained in multiple early branching fish lineages that diverged before teleosts, and (5) marked variation in tissue-specific expression of globin gene repertoires exists across ray-finned fish evolution, including several previously uncharacterized sites of expression. In this respect, our data provide an interesting link between myoglobin expression and the evolution of air breathing in teleosts. Together, our findings demonstrate great-unrecognized diversity in the repertoire and expression of nonblood globins that has arisen during ray-finned fish evolution.

Regulatory divergence of homeologous Atlantic salmon elovl5 genes following the salmonid-specific whole-genome duplication

Fatty acyl elongase 5 (elovl5) is a critical enzyme in the vertebrate biosynthetic pathway which produces the physiologically essential long-chain polyunsaturated fatty acids (LC-PUFA), docosahexenoic acid (DHA), and eicosapentenoic acid (EPA) from 18 carbon fatty acids precursors. In contrast to most other vertebrates, Atlantic salmon possess two copies of elovl5 (elovl5a and elovl5b) as a result of a whole-genome duplication (WGD) which occurred at the base of the salmonid lineage. WGDs have had a major influence on vertebrate evolution, providing extra genetic material, enabling neofunctionalization to accelerate adaptation and speciation. However, little is known about the mechanisms by which such duplicated homeologous genes diverge. Here we show that homeologous Atlantic salmon elovl5a and elovl5b genes have been asymmetrically colonised by transposon-like elements. Identical locations and identities of insertions are also present in the rainbow trout duplicate elovl5 genes, but not in the nearest extant representative preduplicated teleost, the northern pike. Both elovl5 salmon duplicates possessed conserved regulatory elements that promoted Srebp1- and Srebp2-dependent transcription, and differences in the magnitude of Srebp response between promoters could be attributed to a tandem duplication of SRE and NF-Y cofactor binding sites in elovl5b. Furthermore, an insertion in the promoter region of elovl5a confers responsiveness to Lxr/Rxr transcriptional activation. Our results indicate that most, but not all, transposon mobilisation into elovl5 genes occurred after the split from the common ancestor of pike and salmon, but before more recent salmonid speciations, and that divergence of elovl5 regulatory regions have enabled neofuntionalization by promoting differential expression of these homeologous genes.

Divergent spatial regulation of duplicated fatty acid-binding protein (fabp) genes in rainbow trout (Oncorhynchus mykiss)

The increased use of plant oil as a dietary supplement with the resultant high dietary lipid loads challenges the lipid transport, metabolism and storage mechanisms in economically important aquaculture species, such as rainbow trout. Fatty acid-binding proteins (Fabp), ubiquitous in tissues highly active in fatty acid metabolism, participate in lipid uptake and transport, and overall lipid homeostasis. In the present study, searches of nucleotide sequence databases identified mRNA transcripts coded by 14 different fatty acid-binding protein (fabp) genes in rainbow trout (Oncorhynchus mykiss), which include the complete minimal suite of seven distinct fabp genes (fabp1, 2, 3, 6, 7, 10 and 11) discovered thus far in teleost fishes. Phylogenetic analyses suggest that many of these extant fabp genes in rainbow trout exist as duplicates, which putatively arose owing to the teleost-specific whole genome duplication (WGD); three pairs of duplicated fabp genes (fabp2a.1/fabp2a.2, fabp7b.1/fabp7b.2 and fabp10a.1/fabp10a.2) most likely were generated by the salmonid-specific WGD subsequent to the teleost-specific WGD; and fabp3 and fabp6 exist as single copy genes in the rainbow trout genome. Assay of the steady-state levels of fabp gene transcripts by RT-qPCR revealed: (1) steady-state transcript levels differ substantially between fabp genes and, in some instances, by as much as 30×104-fold; (2) some fabp transcripts are widely distributed in many tissues, whereas others are restricted to one or a few tissues; and (3) divergence of regulatory mechanisms that control spatial transcription of duplicated fabp genes in rainbow trout appears related to length of time since their duplication. The suite of fabp genes described here provides the foundation to investigate the role(s) of fatty acid-binding proteins in the uptake, mobilization and storage of fatty acids in cultured fish fed diets differing in lipid content, especially the use of plant oil as a dietary supplement. These nutritional dietary supplements may well lead to high lipid loads with the resultant challenges to lipid homeostasis and, thus, health of cultivated fish which may be mediated by appropriate transcriptional control of fabp genes.

A synthetic rainbow trout linkage map provides new insights into the salmonid whole genome duplication and the conservation of synteny among teleosts

BackgroundRainbow trout is an economically important fish and a suitable experimental organism in many fields of biology including genome evolution, owing to the occurrence of a salmonid specific whole-genome duplication (4th WGD). Rainbow trout is among some of the most studied teleosts and has benefited from substantial efforts to develop genomic resources (e.g., linkage maps. Here, we first generated a synthetic map by merging segregation data files derived from three independent linkage maps. Then, we used it to evaluate genome conservation between rainbow trout and three teleost models, medaka, stickleback and zebrafish and to further investigate the extent of the 4th WGD in trout genome.ResultsThe INRA linkage map was updated by adding 211 new markers. After standardization of marker names, consistency of marker assignment to linkage groups and marker orders was checked across the three different data sets and only loci showing consistent location over all or almost all of the data sets were kept. This resulted in a synthetic map consisting of 2226 markers and 29 linkage groups spanning over 3600 cM. Blastn searches against medaka, stickleback, and zebrafish genomic databases resulted in 778, 824 and 730 significant hits respectively while blastx searches yielded 505, 513 and 510 significant hits. Homology search results revealed that, for most rainbow trout chromosomes, large syntenic regions encompassing nearly whole chromosome arms have been conserved between rainbow trout and its closest models, medaka and stickleback. Large conserved syntenies were also found between the genomes of rainbow trout and the reconstructed teleost ancestor. These syntenies consolidated the known homeologous affinities between rainbow trout chromosomes due to the 4th WGD and suggested new ones.ConclusionsThe synthetic map constructed herein further highlights the stability of the teleost genome over long evolutionary time scales. This map can be easily extended by incorporating new data sets and should help future rainbow trout whole genome sequence assembly. Finally, the persistence of large conserved syntenies across teleosts should facilitate the identification of candidate genes through comparative mapping, even if the occurrence of intra-chromosomal micro-rearrangement may hinder the accurate prediction their genomic location.

Ovary‐predominant wnt4 expression during gonadal differentiation is not conserved in the rainbow trout (Oncorhynchus mykiss)

The Wnt/β-catenin pathway is crucial for ovarian differentiation in mammals, and WNT4 is an important protein that regulates this process. While the role of Wnt4 in gonadal differentiation is relatively well characterized in mammals, little is known regarding its role in teleost fish. Therefore, we investigated the potential activity of wnt4 in gonadal differentiation in rainbow trout (Oncorhynchus mykiss), focusing on the teleost and salmonid gene duplications. Phylogenetic and synteny analyses demonstrated that teleost fish possess two wnt4 genes, wnt4a and wnt4b, as a consequence of the teleost-specific whole-genome duplication (3R). In rainbow trout, we also identified an additional wnt4 gene, which is a wnt4a paralog that likely resulted from the salmonid-specific whole-genome duplication (4R). These two Wnt4a proteins (Wnt4a1 and Wnt4a2) share a high identity (>80%) with other vertebrate Wnt4 proteins, whereas Wnt4b is clearly more divergent (60% identity). During embryogenesis and adulthood, the wnt4a1/2 transcripts were expressed in various tissues, including the ovaries and testes. In contrast, wnt4b expression was restricted to the nervous system, suggesting a sub- or a neo-functionalization of this divergent paralog. During early gonadal differentiation in both males and females, the wnt4a1/2 transcripts were detected in the somatic cells surrounding the germ cells, with a slight sexual dimorphism in favor of males. These results demonstrate that, unlike mammals, rainbow trout do not display an ovary-predominant wnt4 expression profile during early gonadal differentiation.

The genetic basis of salinity tolerance traits in Arctic charr (Salvelinus alpinus)

BackgroundThe capacity to maintain internal ion homeostasis amidst changing conditions is particularly important for teleost fishes whose reproductive cycle is dependent upon movement from freshwater to seawater. Although the physiology of seawater osmoregulation in mitochondria-rich cells of fish gill epithelium is well understood, less is known about the underlying causes of inter- and intraspecific variation in salinity tolerance. We used a genome-scan approach in Arctic charr (Salvelinus alpinus) to map quantitative trait loci (QTL) correlated with variation in four salinity tolerance performance traits and six body size traits. Comparative genomics approaches allowed us to infer whether allelic variation at candidate gene loci (e.g., ATP1α1b, NKCC1, CFTR, and cldn10e) could have underlain observed variation.ResultsCombined parental analyses yielded genome-wide significant QTL on linkage groups 8, 14 and 20 for salinity tolerance performance traits, and on 1, 19, 20 and 28 for body size traits. Several QTL exhibited chromosome-wide significance. Among the salinity tolerance performance QTL, trait co-localizations occurred on chromosomes 1, 4, 7, 18 and 20, while the greatest experimental variation was explained by QTL on chromosomes 20 (19.9%), 19 (14.2%), 4 (14.1%) and 12 (13.1%). Several QTL localized to linkage groups exhibiting homeologous affinities, and multiple QTL mapped to regions homologous with the positions of candidate gene loci in other teleosts. There was no gene × environment interaction among body size QTL and ambient salinity.ConclusionsVariation in salinity tolerance capacity can be mapped to a subset of Arctic charr genomic regions that significantly influence performance in a seawater environment. The detection of QTL on linkage group 12 was consistent with the hypothesis that variation in salinity tolerance may be affected by allelic variation at the ATP1α1b locus. IGF2 may also affect salinity tolerance capacity as suggested by a genome-wide QTL on linkage group 19. The detection of salinity tolerance QTL in homeologous regions suggests that candidate loci duplicated from the salmonid-specific whole-genome duplication may have retained their function on both sets of homeologous chromosomes. Homologous affinities suggest that loci affecting salinity tolerance in Arctic charr may coincide with QTL for smoltification and salinity tolerance traits in rainbow trout. The effects of body size QTL appear to be independent of changes in ambient salinity.

Genomic organization and evolution of the Atlantic salmon hemoglobin repertoire.

BackgroundThe genomes of salmonids are considered pseudo-tetraploid undergoing reversion to a stable diploid state. Given the genome duplication and extensive biological data available for salmonids, they are excellent model organisms for studying comparative genomics, evolutionary processes, fates of duplicated genes and the genetic and physiological processes associated with complex behavioral phenotypes. The evolution of the tetrapod hemoglobin genes is well studied; however, little is known about the genomic organization and evolution of teleost hemoglobin genes, particularly those of salmonids. The Atlantic salmon serves as a representative salmonid species for genomics studies. Given the well documented role of hemoglobin in adaptation to varied environmental conditions as well as its use as a model protein for evolutionary analyses, an understanding of the genomic structure and organization of the Atlantic salmon α and β hemoglobin genes is of great interest.ResultsWe identified four bacterial artificial chromosomes (BACs) comprising two hemoglobin gene clusters spanning the entire α and β hemoglobin gene repertoire of the Atlantic salmon genome. Their chromosomal locations were established using fluorescence in situ hybridization (FISH) analysis and linkage mapping, demonstrating that the two clusters are located on separate chromosomes. The BACs were sequenced and assembled into scaffolds, which were annotated for putatively functional and pseudogenized hemoglobin-like genes. This revealed that the tail-to-tail organization and alternating pattern of the α and β hemoglobin genes are well conserved in both clusters, as well as that the Atlantic salmon genome houses substantially more hemoglobin genes, including non-Bohr β globin genes, than the genomes of other teleosts that have been sequenced.ConclusionsWe suggest that the most parsimonious evolutionary path leading to the present organization of the Atlantic salmon hemoglobin genes involves the loss of a single hemoglobin gene cluster after the whole genome duplication (WGD) at the base of the teleost radiation but prior to the salmonid-specific WGD, which then produced the duplicated copies seen today. We also propose that the relatively high number of hemoglobin genes as well as the presence of non-Bohr β hemoglobin genes may be due to the dynamic life history of salmon and the diverse environmental conditions that the species encounters.Data deposition: BACs S0155C07 and S0079J05 (fps135): GenBank GQ898924; BACs S0055H05 and S0014B03 (fps1046): GenBank GQ898925