De novo European eel transcriptome provides insights into the evolutionary history of duplicated genes in teleost lineages

Superoxidase dismutases (SODs) in the European eel: Gene characterization, expression response to temperature combined with hormonal maturation and possible migratory implications

Superoxide dismutases (SODs) are antioxidant enzymes that protect cells from oxidation. Three SODs have been identified in mammals, but there is limited information in teleosts. This study investigates SODs in the European eel and their expression patterns during testis maturation. Phylogenetic and synteny analyses revealed SODs paralogs and their evolution in vertebrates. The eel possesses one SOD1 and two SOD2/3 (a and b), indicating SOD2 and SOD3 duplication in elopomorphs. SODs expression were then evaluated in various male and female tissues. SOD1 is more expressed in females, while SOD2a and SOD2b dominate brain-pituitary-gonad tissues in both sexes. SOD3a showed predominant expression in the ovary and the male livers, whereas SOD3b was found in the pituitary and brain of both sexes. The effects of different maturation protocols (standard hormonal treatment vs. same protocol preceded with cold seawater pre-treatment) on SODs expression during testis maturation were evaluated. Salinity increase at the onset of standard treatment at 20 °C, simulating early migration, upregulated SOD1, SOD2a, and SOD2b, coinciding with spermatogonia type A differentiated cells dominance. Thereafter, SOD2a and SOD3a decreased, while SOD2b increased during hormonal treatment-induced spermatogenesis. Pre-treatment with seawater at 10 °C, mimicking the conditions at the beginning of the seawater migration, downregulated SOD1 but increased SOD3a expression. Finally, the standard hormonal treatment, replicating spawning at higher temperatures, downregulated SOD1 in eels without any pre-treatment while SOD2a expression increased in pre-treated eels. This study revealed tissue-specific, sex-dependent, and maturation-related SOD expression patterns, predicting SODs dynamic expression profiles during their reproductive migration.

Evolutionary history and functional characterization of duplicated G protein-coupled estrogen receptors in European sea bass

Across vertebrates, the numerous estrogenic functions are mainly mediated by nuclear and membrane receptors, including the G protein-coupled estrogen receptor (GPER) that has been mostly associated with rapid non-genomic responses. Although Gper-mediated signalling has been characterized in only few fish species, Gpers in fish appear to present more mechanistic functionalities as those of mammals due to additional gene duplicates. In this study, we ran a thorough investigation of the fish Gper evolutionary history in light of available genomes, we carried out the functional characterization of the two gper gene duplicates of European sea bass (Dicentrarchus labrax) using luciferase reporter gene transactivation assays, validated it with natural and synthetic estrogen agonists/antagonists and applied it to other chemicals of aquaculture and ecotoxicological interest. Phylogenetic and synteny analyses of fish gper1 and gper1-like genes suggest their duplication may have not resulted from the teleost-specific whole genome duplication. We confirmed that both sbsGper isoforms activate the cAMP signalling pathway and respond differentially to distinct estrogenic compounds. Therefore, as observed for nuclear estrogen receptors, both sbsGpers duplicates retain estrogenic activity although they differ in their specificity and potency (Gper1 being more potent and more specific than Gper1-like), suggesting a more conserved role for Gper1 than for Gper1-like. In addition, Gpers were able to respond to estrogenic environmental pollutants known to interfere with estrogen signalling, such as the phytoestrogen genistein and the anti-depressant fluoxetine, a point that can be taken into account in aquatic environment pollution screenings and chemical risk assessment, complementing previous assays for sea bass nuclear estrogen receptors.

Two CcCCL19bs orchestrate an antibacterial immune response in Yellow River carp (Cyprinus carpio haematopterus)

Chemokines are a group of chemotactic cytokines with an essential role in homeostasis as well as immunity via specific G protein-coupled receptors and atypical receptors. In our study, two Yellow River carp (Cyprinus carpio haematopterus) CCL19b genes (CcCCL19bs), tentatively named CcCCL19b_a and CcCCL19b_b, were cloned. The open reading frames (ORFs) of CcCCL19b_a and CcCCL19b_b were both 333 bp that encoded a 12 kDa protein with 110 amino acid residues. CcCCL19bs contained a signal peptide and a SCY domain with four typical conserved cysteine residues. The two CcCCL19b proteins shared high similarities with each other in both secondary and three-dimensional structure. Phylogenetic analysis showed that CcCCL19bs and other CCL19bs from tetraploid cyprinid fish were clustered into one clade. CcCCL19bs were highly expressed in gill and intestine in healthy fish, and a significant up-regulation of gene expression after Aeromonas hydrophila infection and poly(I:C) stimulation was observed in gill, liver, and head kidney. Furthermore, chemotaxis and antibacterial activity of CcCCL19bs were studied. The results indicated that recombinant CcCCL19b_a and CcCCL19b_b protein (rCcCCL19b_a and rCcCCL19b_b) exhibited significant attraction to primary head kidney leukocytes (HKLs). Meanwhile, both of rCcCCL19bs could promote the proliferation of HKLs, and significantly up-regulate the expressions of IL-1β, CCR7, and IL-6, and down-regulate the expression of IL-10 in primary HKLs. In vitro, rCcCCL19bs could bind and aggregate A. hydrophila and Staphylococcus aureus. The rCcCCL19bs exhibited significant antibacterial activity against A. hydrophila, but not S. aureus. Moreover, they inhibited the growth of A. hydrophila and S. aureus. In vivo, overexpression of CcCCL19bs contributed to the bacterial clearance. These studies suggested that CcCCL19bs orchestrate an antibacterial immune response.

Independent rediploidization masks shared whole genome duplication in the sturgeon-paddlefish ancestor

Whole genome duplication (WGD) is a dramatic evolutionary event generating many new genes and which may play a role in survival through mass extinctions. Paddlefish and sturgeon are sister lineages that both show genomic evidence for ancient WGD. Until now this has been interpreted as two independent WGD events due to a preponderance of duplicate genes with independent histories. Here we show that although there is indeed a plurality of apparently independent gene duplications, these derive from a shared genome duplication event occurring well over 200 million years ago, likely close to the Permian-Triassic mass extinction period. This was followed by a prolonged process of reversion to stable diploid inheritance (rediploidization), that may have promoted survival during the Triassic-Jurassic mass extinction. We show that the sharing of this WGD is masked by the fact that paddlefish and sturgeon lineage divergence occurred before rediploidization had proceeded even half-way. Thus, for most genes the resolution to diploidy was lineage-specific. Because genes are only truly duplicated once diploid inheritance is established, the paddlefish and sturgeon genomes are thus a mosaic of shared and non-shared gene duplications resulting from a shared genome duplication event.

Identification of type II interferons and receptors in an osteoglossiform fish, the arapaima Arapaima gigas

In mammals, type II interferon (IFN; i.e. IFN-γ) signalling transduces through its specific receptors IFN-γR1 and IFN-γR2. In an osteoglossiform fish, the arapaima Arapaima gigas, three type II IFNs, IFN-γ-like, IFN-γ and IFN-γrel, and their four possible receptor subunits IFN-γR1-1, IFN-γR1-2, IFN-γR2-1 and IFN-γR2-2 were identified in this study. The three type II IFN genes are composed of four exons and three introns, and they all contain IFN-γ signature motif and signal peptide, with the presence of potential nuclear localization signal (NLS) in IFN-γ-like and IFN-γ. The IFN-γR1-1, IFN-γR1-2, IFN-γR2-1 and IFN-γR2-2 are composed of seven exons and six introns, with predicted IFN-γR1-1 and IFN-γR1-2 proteins containing JAK1 and STAT1 binding sites, and IFN-γR2-1 and IFN-γR2-2 containing JAK2 binding sites. Gene synteny analysis showed that the type II IFN and their receptor loci are duplicated in arapaima. All these genes were expressed constitutively in all organs/tissues examined, and responded to the stimulation of polyI:C. The prokaryotic recombinant IFN-γ-like, IFN-γ and IFN-γrel proteins can significantly induce the upregulation of immune-related genes in trunk kidney leucocytes. The ligand-receptor relationship analyses revealed that recombinant IFN-γ-like, IFN-γ, and IFN-γrel transduce downstream signalling through IFN-γR1-1/IFN-γR2-1, IFN-γR1-2/IFN-γR2-2, and IFN-γR1-1, respectively, in xenogeneic cells with the overexpression of original or chimeric receptors. In addition, tyrosine (Y) 366 and Y377 in the intracellular region may be essential for the function of IFN-γR1-2 and IFN-γR1-1, respectively. The finding of type II IFN system in A. gigas thus provides different knowledge in understanding the diversity and evolution of type II IFN ligand-receptor relationships in vertebrates.

A Chromosome-level assembly of the Japanese eel genome, insights into gene duplication and chromosomal reorganization

Japanese eels (Anguilla japonica) are commercially important species, harvested extensively for food. Currently, this and related species (American and European eels) are challenging to breed on a commercial basis. As a result, the wild stock is used for aquaculture. Moreover, climate change, habitat loss, water pollution, and altered ocean currents affect eel populations negatively. Accordingly, the International Union for Conservation of Nature lists Japanese eels as endangered and on its red list. Here we presented a high-quality genome assembly for Japanese eels and demonstrated that large chromosome reorganizations occurred in the events of third-round whole-genome duplications (3R-WRDs). Several chromosomal fusions and fissions have reduced the ancestral protochromosomal number of 25 to 19 in the Anguilla lineage. A phylogenetic analysis of the expanded gene families showed that the olfactory receptors (group δ and ζ genes) and voltage-gated Ca2+ channels expanded significantly. Both gene families are crucial for olfaction and neurophysiology. Additional tandem and proximal duplications occurred following 3R-WGD to acquire immune-related genes for an adaptive advantage against various pathogens. The Japanese eel assembly presented here can be used to study other Anguilla species relating to evolution and conservation.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

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

Genome-wide identification, evolutionary analysis, and antimicrobial activity prediction of CC chemokines in allotetraploid common carp, Cyprinus carpio

Chemokines are a group of secreted small molecules which are essential for cell migration in physiological and pathological conditions by binding to specific chemokine receptors. They are structurally classified into five groups, namely CXC, CC, CX3C, XC and CX. CC chemokine group is the largest one among them. In this study, we identified and characterized 61 CC chemokines from allotetraploid common carp (Cyprinus carpio). The sequence analyses showed that the majority of CC chemokines had an N-terminal signal peptide, and an SCY domain, and all CC chemokines were located in the extracellular region. Phylogenetic, evolutionary and syntenic analyses confirmed that CC chemokines were annotated as 11 different types (CCL19, CCL20, CCL25, CCL27, CCL32, CCL33, CCL34, CCL35, CCL36, CCL39, and CCL44), which exhibited unique gene arrangement pattern and chromosomal location respectively. Furthermore, genome synteny analyses between common carp and four representative teleost species indicated expansion of common carp CC chemokines resulted from the whole genome duplication (WGD) event. Additionally, the continuous evolution of gene CCL25s in teleost afforded a novel viewpoint to explain the WGD event in teleost. Then, we predicted the three-dimensional structures and probable function regions of common carp CC chemokines. All the CC chemokines core structures were constituted of an N-loop, a three-stranded β-sheet, and a C-terminal helix. Finally, 43 CC chemokines were predicted to have probable general antimicrobial activity. Their tertiary structures, cationic and amphiphilic physicochemical property supported the viewpoint. To verify the prediction, six recombinant CCL19s proteins were prepared and the antibacterial activity against Escherichia coli and Aeromonas hydrophila were verified. The results supported our prediction that rCCL19a.1s (rCCL19a.1_a, rCCL19a.1_b) and rCCL19bs (rCCL19b_a, rCCL19b_b), especially rCCL19bs, exhibited extremely significant inhibition to the growth of both E. coli and A. hydrophila. On the contrary, two rCCL19a.2s had no significant inhibitory effect. These studies suggested that CC chemokines were essential in immune system evolution and not monofunctional during pathogen infection.

An atlas of fish genome evolution reveals delayed rediploidization following the teleost whole-genome duplication

Teleost fishes are ancient tetraploids descended from an ancestral whole-genome duplication that may have contributed to the impressive diversification of this clade. Whole-genome duplications can occur via self-doubling (autopolyploidy) or via hybridization between different species (allopolyploidy). The mode of tetraploidization conditions evolutionary processes by which duplicated genomes return to diploid meiotic pairing, and subsequent genetic divergence of duplicated genes (cytological and genetic rediploidization). How teleosts became tetraploid remains unresolved, leaving a fundamental gap in the interpretation of their functional evolution. As a result of the whole-genome duplication, identifying orthologous and paralogous genomic regions across teleosts is challenging, hindering genome-wide investigations into their polyploid history. Here, we combine tailored gene phylogeny methodology together with a state-of-the-art ancestral karyotype reconstruction to establish the first high-resolution comparative atlas of paleopolyploid regions across 74 teleost genomes. We then leverage this atlas to investigate how rediploidization occurred in teleosts at the genome-wide level. We uncover that some duplicated regions maintained tetraploidy for more than 60 million years, with three chromosome pairs diverging genetically only after the separation of major teleost families. This evidence suggests that the teleost ancestor was an autopolyploid. Further, we find evidence for biased gene retention along several duplicated chromosomes, contradicting current paradigms that asymmetrical evolution is specific to allopolyploids. Altogether, our results offer novel insights into genome evolutionary dynamics following ancient polyploidizations in vertebrates.

Stimulatory Effects of Androgens on Eel Primary Ovarian Development - from Phenotypes to Genotypes

Androgens stimulate primary ovarian development in Vertebrate. Japanese eels underwent operation to sample the pre- and post-treated ovarian tissues from the same individual. Ovarian phenotypic or genotypic data were mined in a pair. A correlation between the initial ovarian status (determined by kernel density estimation (KDE), presented as a probability density of oocyte size) and the consequence of androgen (17MT) treatment (change in ovary) has been showed. The initial ovarian status appeared to be important to influence ovarian androgenic sensitivity. The initial ovary was important to the outcomes of androgen treatments, and ePAV (expression presence-absence variation) is existing in Japanese eel by analyze DEGs; core, unique, or accessory genes were identified, the sensitivities of initial ovaries were correlated with their gene expression profiles. We speculated the importance of genetic differential expression on the variations of phenotypes by 17MT, and transcriptomic approach seems to allow extracting multiple layers of genomic data.

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

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

Investigating the Transcriptomic and Expression Presence-Absence Variation Exist in Japanese Eel (Anguilla japonica), a Primitive Teleost

The pan-genome was defined as the complete gene set across strains, and it is built upon genes displaying presence-absence variations (PAVs); the pan-transcriptome is defined by recalling the pan-genome. Indeed, a PAV is reflected from the expression presence-absence variation (ePAV). In this study, treated with androgen, eels, which are a primitive fish from the basal lineage of Teleost, with different ovarian developments were chosen and submitted to RAN-sequencing. Transcriptomes were the assembly against eel genome scaffolds; a pair was the unit (the same eel before and after treatment) to analyze DEGs (differentially expressed genes); the core, unique, or accessory genes were identified, and the list of DEGs was analyzed to investigate ePAV. The results suggest that there was ePAV in Japanese eel, and the ePAV of eel was analyzed by pathway enrichment. These results signify the importance of genetic differential expression on the variations of phenotypes by androgen, and a transcriptomic approach appears to enable extracting multiple layers of genomic data.

DNA Transposon Expansion is Associated with Genome Size Increase in Mudminnows

Genome sizes of eukaryotic organisms vary substantially, with whole-genome duplications (WGD) and transposable element expansion acting as main drivers for rapid genome size increase. The two North American mudminnows, Umbra limi and Umbra pygmaea, feature genomes about twice the size of their sister lineage Esocidae (e.g., pikes and pickerels). However, it is unknown whether all Umbra species share this genome expansion and which causal mechanisms drive this expansion. Using flow cytometry, we find that the genome of the European mudminnow is expanded similarly to both North American species, ranging between 4.5 and 5.4 pg per diploid nucleus. Observed blocks of interstitially located telomeric repeats in U. limi suggest frequent Robertsonian rearrangements in its history. Comparative analyses of transcriptome and genome assemblies show that the genome expansion in Umbra is driven by the expansion of DNA transposon and unclassified repeat sequences without WGD. Furthermore, we find a substantial ongoing expansion of repeat sequences in the Alaska blackfish Dallia pectoralis, the closest relative to the family Umbridae, which might mark the beginning of a similar genome expansion. Our study suggests that the genome expansion in mudminnows, driven mainly by transposon expansion, but not WGD, occurred before the separation into the American and European lineage.

Mobile Elements in Ray-Finned Fish Genomes

Ray-finned fishes (Actinopterygii) are a very diverse group of vertebrates, encompassing species adapted to live in freshwater and marine environments, from the deep sea to high mountain streams. Genome sequencing offers a genetic resource for investigating the molecular bases of this phenotypic diversity and these adaptations to various habitats. The wide range of genome sizes observed in fishes is due to the role of transposable elements (TEs), which are powerful drivers of species diversity. Analyses performed to date provide evidence that class II DNA transposons are the most abundant component in most fish genomes and that compared to other vertebrate genomes, many TE superfamilies are present in actinopterygians. Moreover, specific TEs have been reported in ray-finned fishes as a possible result of an intricate relationship between TE evolution and the environment. The data summarized here underline the biological interest in Actinopterygii as a model group to investigate the mechanisms responsible for the high biodiversity observed in this taxon.

Vertebrate Alpha2,8-Sialyltransferases (ST8Sia): A Teleost Perspective

We identified and analyzed α2,8-sialyltransferases sequences among 71 ray-finned fish species to provide the first comprehensive view of the Teleost ST8Sia repertoire. This repertoire expanded over the course of Vertebrate evolution and was primarily shaped by the whole genome events R1 and R2, but not by the Teleost-specific R3. We showed that duplicated st8sia genes like st8sia7, st8sia8, and st8sia9 have disappeared from Tetrapods, whereas their orthologues were maintained in Teleosts. Furthermore, several fish species specific genome duplications account for the presence of multiple poly-α2,8-sialyltransferases in the Salmonidae (ST8Sia II-r1 and ST8Sia II-r2) and in Cyprinuscarpio (ST8Sia IV-r1 and ST8Sia IV-r2). Paralogy and synteny analyses provided more relevant and solid information that enabled us to reconstruct the evolutionary history of st8sia genes in fish genomes. Our data also indicated that, while the mammalian ST8Sia family is comprised of six subfamilies forming di-, oligo-, or polymers of α2,8-linked sialic acids, the fish ST8Sia family, amounting to a total of 10 genes in fish, appears to be much more diverse and shows a patchy distribution among fish species. A focus on Salmonidae showed that (i) the two copies of st8sia2 genes have overall contrasted tissue-specific expressions, with noticeable changes when compared with human co-orthologue, and that (ii) st8sia4 is weakly expressed. Multiple sequence alignments enabled us to detect changes in the conserved polysialyltransferase domain (PSTD) of the fish sequences that could account for variable enzymatic activities. These data provide the bases for further functional studies using recombinant enzymes.