Rhodopsin gene copies in Japanese eel originated in a teleost-specific genome duplication

Functional diversification process of opsin genes for teleost visual and pineal photoreceptions

Most vertebrates have a rhodopsin gene with a five-exon structure for visual photoreception. By contrast, teleost fishes have an intron-less rhodopsin gene for visual photoreception and an intron-containing rhodopsin (exo-rhodopsin) gene for pineal photoreception. Here, our analysis of non-teleost and teleost fishes in various lineages of the Actinopterygii reveals that retroduplication after branching of the Polypteriformes produced the intron-less rhodopsin gene for visual photoreception, which converted the parental intron-containing rhodopsin gene into a pineal opsin in the common ancestor of the Teleostei. Additional analysis of a pineal opsin, pinopsin, shows that the pinopsin gene functions as a green-sensitive opsin together with the intron-containing rhodopsin gene for pineal photoreception in tarpon as an evolutionary intermediate state but is missing in other teleost fishes, probably because of the redundancy with the intron-containing rhodopsin gene. We propose an evolutionary scenario where unique retroduplication caused a "domino effect" on the functional diversification of teleost visual and pineal opsin genes.

Molecular Dating of the Teleost Whole Genome Duplication (3R) Is Compatible With the Expectations of Delayed Rediploidization

Vertebrate evolution has been punctuated by three whole genome duplication events that have been implicated causally in phenotypic evolution, from the origin of phenotypic novelties to explosive diversification. Arguably, the most dramatic of these is the 3R whole genome duplication event associated with the origin of teleost fishes which comprise more than half of all living vertebrate species. However, tests of a causal relationship between whole genome duplication and teleost diversification have proven difficult due to the challenge of establishing the timing of these phenomena. Here we show, based on molecular clock dating of concatenated gene alignments, that the 3R whole genome duplication event occurred in the early-middle Permian (286.18 to 267.20 million years ago; Ma), 52.02 to 12.84 million years (Myr) before the divergence of crown-teleosts in the latest Permian-earliest Late Triassic (254.36 to 234.16 Ma) and long before the major pulses of teleost diversification in Ostariophysi and Percomorpha (56.37 to 100.17 Myr and at least 139.24 to 183.29 Myr later, respectively). The extent of this temporal gap between putative cause and effect precludes 3R as a deterministic driver of teleost diversification. However, these age constraints remain compatible with the expectations of a prolonged rediploidization process following whole genome duplication which, through the effects of chromosome rearrangement and gene loss, remains a viable mechanism to explain the evolution of teleost novelties and diversification.

Characterization of potential spermatogonia biomarker genes in the European eel (Anguilla anguilla)

Identification of specific molecular markers for spermatogonial stem cells in teleost is crucial for enhancing the efficacy of reproductive biotechnologies in aquaculture, such as transplantation and surrogate production in fishes. Since it is not yet possible to distinguish spermatogonial stem cells of European eel (Anguilla anguilla) using specific molecular markers, we isolated spermatogonial cells from immature European eels to find these potential markers. We attempted this by studying three candidate genes: vasa, nanos2, and dnd1. Two vasa (vasa1 and vasa2) genes, nanos2, and dnd1 were identified, characterized, and studied in the muscle, testis, and isolated spermatogonia. Our results showed that vasa1 and vasa2 had the highest levels of expression when measured by qPCR. In situ hybridization and immunochemistry assays showed that the four genes were localized explicitly in type A spermatogonia. However, vasa1 and vasa2 exhibited stronger signals in the immature testicular tissue than the other two potential markers. According to this, vasa1 and vasa2 were found to be the most effective markers for spermatogonial cells in the European eel.

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.

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.

Recent genome duplications facilitate the phenotypic diversity of Hb repertoire in the Cyprinidae

Whole-genome duplications (WGDs) are an important contributor to phenotypic innovations in evolutionary history. The diversity of blood oxygen transport traits is the perfect reflection of physiological versatility for evolutionary success among vertebrates. In this study, the evolutionary changes of hemoglobin (Hb) repertoire driven by the recent genome duplications were detected in representative Cyprinidae fish, including eight diploid and four tetraploid species. Comparative genomic analysis revealed a substantial variation in both membership composition and intragenomic organization of Hb genes in these species. Phylogenetic reconstruction analyses were conducted to characterize the evolutionary history of these genes. Data were integrated with the expression profiles of the genes during ontogeny. Our results indicated that genome duplications facilitated the phenotypic diversity of the Hb gene family; each was associated with species-specific changes in gene content via gene loss and fusion after genome duplications. This led to repeated evolutionary transitions in the ontogenic regulation of Hb gene expression. Our results revealed that genome duplications helped to generate phenotypic changes in Cyprinidae Hb systems.

Visual pigment evolution in Characiformes: The dynamic interplay of teleost whole-genome duplication, surviving opsins and spectral tuning

Vision represents an excellent model for studying adaptation, given the genotype-to-phenotype map that has been characterized in a number of taxa. Fish possess a diverse range of visual sensitivities and adaptations to underwater light, making them an excellent group to study visual system evolution. In particular, some speciose but understudied lineages can provide a unique opportunity to better understand aspects of visual system evolution such as opsin gene duplication and neofunctionalization. In this study, we showcase the visual system evolution of neotropical Characiformes and the spectral tuning mechanisms they exhibit to modulate their visual sensitivities. Such mechanisms include gene duplications and losses, gene conversion, opsin amino acid sequence and expression variation, and A₁ /A₂ -chromophore shifts. The Characiforms we studied utilize three cone opsin classes (SWS2, RH2, LWS) and a rod opsin (RH1). However, the characiform's entire opsin gene repertoire is a product of dynamic evolution by opsin gene loss (SWS1, RH2) and duplication (LWS, RH1). The LWS- and RH1-duplicates originated from a teleost specific whole-genome duplication as well as characiform-specific duplication events. Both LWS-opsins exhibit gene conversion and, through substitutions in key tuning sites, one of the LWS-paralogues has acquired spectral sensitivity to green light. These sequence changes suggest reversion and parallel evolution of key tuning sites. Furthermore, characiforms' colour vision is based on the expression of both LWS-paralogues and SWS2. Finally, we found interspecific and intraspecific variation in A₁ /A₂ -chromophores proportions, correlating with the light environment. These multiple mechanisms may be a result of the diverse visual environments where Characiformes have evolved.

Phylogeny of teleost connexins reveals highly inconsistent intra- and interspecies use of nomenclature and misassemblies in recent teleost chromosome assemblies

BACKGROUND

Based on an initial collecting of database sequences from the gap junction protein gene family (also called connexin genes) in a few teleosts, the naming of these sequences appeared variable. The reasons could be (i) that the structure in this family is variable across teleosts, or (ii) unfortunate naming. Rather clear rules for the naming of genes in fish and mammals have been outlined by nomenclature committees, including the naming of orthologous and ohnologous genes. We therefore analyzed the connexin gene family in teleosts in more detail. We covered the range of divergence times in teleosts (eel, Atlantic herring, zebrafish, Atlantic cod, three-spined stickleback, Japanese pufferfish and spotted pufferfish; listed from early divergence to late divergence).

RESULTS

The gene family pattern of connexin genes is similar across the analyzed teleosts. However, (i) several nomenclature systems are used, (ii) specific orthologous groups contain genes that are named differently in different species, (iii) several distinct genes have the same name in a species, and (iv) some genes have incorrect names. The latter includes a human connexin pseudogene, claimed as GJA4P, but which in reality is Cx39.2P (a delta subfamily gene often called GJD2like). We point out the ohnologous pairs of genes in teleosts, and we suggest a more consistent nomenclature following the outlined rules from the nomenclature committees. We further show that connexin sequences can indicate some errors in two high-quality chromosome assemblies that became available very recently.

CONCLUSIONS

Minimal consistency exists in the present practice of naming teleost connexin genes. A consistent and unified nomenclature would be an advantage for future automatic annotations and would make various types of subsequent genetic analyses easier. Additionally, roughly 5% of the connexin sequences point out misassemblies in the new high-quality chromosome assemblies from herring and cod.

Trans-omics analyses revealed differences in hormonal and nutritional status between wild and cultured female Japanese eel (Anguilla japonica)

Long-term stock decline in the Japanese eel (Anguilla japonica) is a serious issue. To reduce natural resource utilization in Japan, artificial hormonal induction of maturation and fertilization in the Japanese eel has been intensively studied. Recent experiment on feminized (by feeding a commercial diet containing estradiol-17β for first half year) cultured female eels have shown ovulation problems, which is seldom observed in captured wild female eels. Therefore, the aim of this study is to try to investigate causes of ovulation problem frequently seen in cultured female eels by comparative trans-omics analyses. The omics data showed low growth hormone and luteinizing hormone transcription levels in the brain and low sex hormone–binding globulin transcription levels in the liver of the cultured female eels. In addition, it was found that high accumulation of glucose-6-phosphate and, maltose in the liver of the cultured female eel. It was also found that docosahexaenoic (DHA) acid, eicosapentaenoic acid (EPA) and arachidonic acid (ARA) ratios in cultured female eels were quite different from wild female eels. The data suggested that ovulation problem in cultured female eels was possibly resulted from prolonged intake of a high-carbohydrate diet and/or suboptimal DHA/EPA/ARA ratios in a diet.

Rhodopsin gene evolution in early teleost fishes

Rhodopsin mediates an essential step in image capture and is tightly associated with visual adaptations of aquatic organisms, especially species that live in dim light environments (e.g., the deep sea). The rh1 gene encoding rhodopsin was formerly considered a single-copy gene in genomes of vertebrates, but increasing exceptional cases have been found in teleost fish species. The main objective of this study was to determine to what extent the visual adaptation of teleosts might have been shaped by the duplication and loss of rh1 genes. For that purpose, homologous rh1/rh1-like sequences in genomes of ray-finned fishes from a wide taxonomic range were explored using a PCR-based method, data mining of public genetic/genomic databases, and subsequent phylogenomic analyses of the retrieved sequences. We show that a second copy of the fish-specific intron-less rh1 is present in the genomes of most anguillids (Elopomorpha), Hiodon alosoides (Osteoglossomorpha), and several clupeocephalan lineages. The phylogenetic analysis and comparisons of alternative scenarios for putative events of gene duplication and loss suggested that fish rh1 was likely duplicated twice during the early evolutionary history of teleosts, with one event coinciding with the hypothesized fish-specific genome duplication and the other in the common ancestor of the Clupeocephala. After these gene duplication events, duplicated genes were maintained in several teleost lineages, whereas some were secondarily lost in specific lineages. Alternative evolutionary schemes of rh1 and comparison with previous studies of gene evolution are also reviewed.

Whole-Genome Sequencing of 84 Japanese Eels Reveals Evidence against Panmixia and Support for Sympatric Speciation

The Japanese eel (Anguilla japonica), European eel (Anguilla anguilla), and American eel (Anguilla rostrata) are migratory, catadromous, temperate zone fish sharing several common life cycle features. The population genetics of panmixia in these eel species has already been investigated. Our extensive population genetics analysis was based on 1400 Gb of whole-genome sequence (WGS) data from 84 eels. It demonstrated that a Japanese eel group from the Kuma River differed from other populations of the same species. Even after removing the potential adapted/selected single nucleotide polymorphism (SNP) data, and with very small differences (fixation index [Fst] = 0.01), we obtained results consistently indicating that panmixia does not occur in Japanese eels. The life cycle of the Japanese eel is well-established and the Kuma River is in the center of its habitat. Nevertheless, simple reproductive isolation is not the probable cause of non-panmixia in this species. We propose that the combination of spawning area subdivision, philopatry, and habitat preference/avoidance accounts for the non-panmixia in the Japanese eel population. We named this hypothesis the “reproductive isolation like subset mapping” (RISM) model. This finding may be indicative of the initial stages of sympatric speciation in these eels.