Fish-T1K (Transcriptomes of 1,000 Fishes) Project: large-scale transcriptome data for fish evolution studies

Genetic improvement and genomic resources of important cyprinid species: status and future perspectives for sustainable production

Cyprinid species are the most cultured aquatic species around the world in terms of quantity and total value. They account for 25% of global aquaculture production and significantly contribute to fulfilling the demand for fish food. The aquaculture of these species is facing severe concerns in terms of seed quality, rising feed costs, disease outbreaks, introgression of exotic species, environmental impacts, and anthropogenic activities. Numerous researchers have explored biological issues and potential methods to enhance cyprinid aquaculture. Selective breeding is extensively employed in cyprinid species to enhance specific traits like growth and disease resistance. In this context, we have discussed the efforts made to improve important cyprinid aquaculture practices through genetic and genomic approaches. The recent advances in DNA sequencing technologies and genomic tools have revolutionized the understanding of biological research. The generation of a complete genome and other genomic resources in cyprinid species has significantly strengthened molecular-level investigations into disease resistance, growth, reproduction, and adaptation to changing environments. We conducted a comprehensive review of genomic research in important cyprinid species, encompassing genome, transcriptome, proteome, metagenome, epigenome, etc. This review reveals that considerable data has been generated for cyprinid species. However, the seamless integration of this valuable data into genetic selection programs has yet to be achieved. In the upcoming years, genomic techniques, gene transfer, genome editing tools are expected to bring a paradigm shift in sustainable cyprinid aquaculture production. The comprehensive information presented here will offer insights for the cyprinid aquaculture research community.

CoSFISH: a comprehensive reference database of COI and 18S rRNA barcodes for fish

Fish, being a crucial component of aquatic ecosystems, holds significant importance from both economic and ecological perspectives. However, the identification of fish at the species level remains challenging, and there is a lack of a taxonomically complete and comprehensive reference sequence database for fish. Therefore, we developed CoSFISH, an online fish database. Currently, the database contains 21 535 cytochrome oxidase I sequences and 1074 18S rRNA sequences of 21 589 species, belonging to 8 classes and 90 orders. We additionally incorporate online analysis tools to aid users in comparing, aligning and analyzing sequences, as well as designing primers. Users can upload their own data for analysis, in addition to using the data stored in the database directly. CoSFISH offers an extensive fish database and incorporates online analysis tools, making it a valuable resource for the study of fish diversity, phylogenetics and biological evolution. Database URL:  http://210.22.121.250:8888/CoSFISH/home/indexPage.

Transcriptome mining extends the host range of the Flaviviridae to non-bilaterians

The flavivirids (family Flaviviridae) are a group of positive-sense RNA viruses that include well-documented agents of human disease. Despite their importance and ubiquity, the timescale of flavivirid evolution is uncertain. An ancient origin, spanning millions of years, is supported by their presence in both vertebrates and invertebrates and by the identification of a flavivirus-derived endogenous viral element in the peach blossom jellyfish genome (Craspedacusta sowerbii, phylum Cnidaria), implying that the flaviviruses arose early in the evolution of the Metazoa. To date, however, no exogenous flavivirid sequences have been identified in these hosts. To help resolve the antiquity of the Flaviviridae, we mined publicly available transcriptome data across the Metazoa. From this, we expanded the diversity within the family through the identification of 32 novel viral sequences and extended the host range of the pestiviruses to include amphibians, reptiles, and ray-finned fish. Through co-phylogenetic analysis we found cross-species transmission to be the predominate macroevolutionary event across the non-vectored flavivirid genera (median, 68 per cent), including a cross-species transmission event between bats and rodents, although long-term virus-host co-divergence was still a regular occurrence (median, 23 per cent). Notably, we discovered flavivirus-like sequences in basal metazoan species, including the first associated with Cnidaria. This sequence formed a basal lineage to the genus Flavivirus and was closer to arthropod and crustacean flaviviruses than those in the tamanavirus group, which includes a variety of invertebrate and vertebrate viruses. Combined, these data attest to an ancient origin of the flaviviruses, likely close to the emergence of the metazoans 750-800 million years ago.

Aquaculture omics: An update on the current status of research and data analysis

Aquaculture is the fast-growing agricultural sector and has the ability to meet the growing demand for protein nutritional security for future population. In future aquaculture is going to be the major source of fish proteins as capture fisheries reached at its maximum. However, several challenges need to overcome such as lack of genetically improved strains/varieties, lack of species-specific feed/functional feed, round the year availability of quality fish seed, pollution of ecosystems and increased frequencies of disease occurrence etc. In recent years, the continuous development of high throughput sequencing technology has revolutionized the biological sciences and provided necessary tools. Application of 'omics' in aquaculture research have been successfully used to resolve several productive and reproductive issues and thus ensure its sustainability and profitability. To date, high quality draft genomes of over fifty fish species have been generated and successfully used to develop large number of single nucleotide polymorphism markers (SNPs), marker panels and other genomic resources etc in several aquaculture species. Similarly, transcriptome profiling and miRNAs analysis have been used in aquaculture research to identify key transcripts and expression analysis of candidate genes/miRNAs involved in reproduction, immunity, growth, development, stress toxicology and disease. Metagenome analysis emerged as a promising scientific tool to analyze the complex genomes contained within microbial communities. Metagenomics has been successfully used in the aquaculture sector to identify novel and potential pathogens, antibiotic resistance genes, microbial roles in microcosms, microbial communities forming biofloc, probiotics etc. In the current review, we discussed application of high-throughput technologies (NGS) in the aquaculture sector.

Evolutionary Genomics Reveals Multiple Functions of Arylalkylamine N-Acetyltransferase in Fish

As an important hormone, melatonin participates in endocrine regulation of diverse functions in vertebrates. Its biosynthesis is catalyzed by four cascaded enzymes, among them, arylalkylamine N-acetyltransferase (AANAT) is the most critical one. Although only single aanat gene has been identified in most groups of vertebrates, researchers including us have determined that fish have the most diverse of aanat genes (aanat1a, aanat1b, and aanat2), playing various potential roles such as seasonal migration, amphibious aerial vision, and cave or deep-sea adaptation. With the rapid development of genome and transcriptome sequencing, more and more putative sequences of fish aanat genes are going to be available. Related phylogeny and functional investigations will enrich our understanding of AANAT functions in various fish species.

RNA-Seq of Single Fish Cells - Seeking Out the Leukocytes Mediating Immunity in Teleost Fishes

The immune system is a complex and sophisticated biological system, spanning multiple levels of complexity, from the molecular level to that of tissue. Our current understanding of its function and complexity, of the heterogeneity of leukocytes, is a result of decades of concentrated efforts to delineate cellular markers using conventional methods of antibody screening and antigen identification. In mammalian models, this led to in-depth understanding of individual leukocyte subsets, their phenotypes, and their roles in health and disease. The field was further propelled forward by the development of single-cell (sc) RNA-seq technologies, offering an even broader and more integrated view of how cells work together to generate a particular response. Consequently, the adoption of scRNA-seq revealed the unexpected plasticity and heterogeneity of leukocyte populations and shifted several long-standing paradigms of immunology. This review article highlights the unprecedented opportunities offered by scRNA-seq technology to unveil the individual contributions of leukocyte subsets and their crosstalk in generating the overall immune responses in bony fishes. Single-cell transcriptomics allow identifying unseen relationships, and formulating novel hypotheses tailored for teleost species, without the need to rely on the limited number of fish-specific antibodies and pre-selected markers. Several recent studies on single-cell transcriptomes of fish have already identified previously unnoticed expression signatures and provided astonishing insights into the diversity of teleost leukocytes and the evolution of vertebrate immunity. Without a doubt, scRNA-seq in tandem with bioinformatics tools and state-of-the-art methods, will facilitate studying the teleost immune system by not only defining key markers, but also teaching us about lymphoid tissue organization, development/differentiation, cell-cell interactions, antigen receptor repertoires, states of health and disease, all across time and space in fishes. These advances will invite more researchers to develop the tools necessary to explore the immunology of fishes, which remain non-conventional animal models from which we have much to learn.

Genome-wide identification and characterization of 14-3-3 genes in fishes

The 14-3-3 family genes are highly conserved regulatory factors in eukaryotes with involvement in multiple important cellular processes. However, detailed investigations of this family in fishes are very limited. Here, a comparative genomic and transcriptomic survey were performed to investigate the 14-3-3 family in fishes. We confirmed that the numbers of 14-3-3 genes ranged from 5 to 7 in non-teleost fishes, as well as additional 14-3-3 genes (9 to 11) in teleost fishes. In addition, some special teleost fishes possess 17 to 25 14-3-3s, which undergone the fourth whole-genome duplication (WGD). We also found that six pairs of fish 14-3-3 genes were clustered with mammalian ε, γ, ς, η, τand β isotypes, respectively, while σ was absent with a potential specificity within mammals, on the basis of their phylogenetic and synteny analyses. According to our results, we inferred that the diversity of 14-3-3 genes in fishes seems to be generated from a combination of WGD and gene loss. Comparative transcriptomic analysis revealed that there are differences in tissue distribution, and we speculated that 14-3-3 genes may contribute to terrestrial adaptations in mudskippers. In addition, protein sequence alignments of 14-3-3s supported their differential roles in fishes. In summary, our present comparative genomic and transcriptomic survey will benefit for further functional investigations of these fish genes.

Comprehensive phylogeny of Konosirus punctatus (Clupeiformes: Clupeidae) based on transcriptomic data

Konosirus punctatus is an important species for the structure of marine ecosystems. Meanwhile, it is a native species in the northwest Pacific Ocean and supports important commercial fishery. In the present study, we generated the whole transcriptome of K. punctatus from combined tissues (muscle, liver, gill, heart, kidney, swim bladder and sexual gland) using Illumina RNA-seq technology and a total of 46087110 clean reads were obtained, corresponding to 6531521430 nucleotides. Meanwhile, 10000 clean reads were randomly selected and compared with NT database to examine the possible contamination. Results showed that 6754 clean reads were distributed among some species closely related with K. punctatus, indicating no-pollution. De novo assembly was performed and all clean reads were assembled to produce 71610 longest unigenes with an N50 of 906 bp. Among all the unigenes, 43974 unigenes were annotated in at least one database and 3172 unigenes were annotated in all databases. All unigenes were further analyzed to predict the gene structure and we have obtained a total of 54864 coding sequences (CDSs) and 17326 simple sequence repeats (SSRs). Saturation analyses were applied to evaluate the accuracy of gene expression and we hypothesized that the detection of gene expression might be effective. Finally, single-copy orthologous genes were applied to construct the phylogenetic relationship of K. punctatus. Results showed that K. punctatus diverged from the common ancestor with Alosa alosa, Alosa pseudoharengus and Sardina pilchardus at approx. 61.16–92.52 MYA. The present study will provide a foundational molecular information for the biological research of K. punctatus.

Identification and Characterization of Zebrafish Tlr4 Coreceptor Md-2

The zebrafish (Danio rerio) is a powerful model organism for studies of the innate immune system. One apparent difference between human and zebrafish innate immunity is the cellular machinery for LPS sensing. In amniotes, the protein complex formed by TLR4 and myeloid differentiation factor 2 (Tlr4/Md-2) recognizes the bacterial molecule LPS and triggers an inflammatory response. It is believed that zebrafish have neither Md-2 nor Tlr4; Md-2 has not been identified outside of amniotes, whereas the zebrafish tlr4 genes appear to be paralogs, not orthologs, of amniote TLR4s We revisited these conclusions. We identified a zebrafish gene encoding Md-2, ly96 Using single-cell RNA sequencing, we found that ly96 is transcribed in cells that also transcribe genes diagnostic for innate immune cells, including the zebrafish tlr4-like genes. In larval zebrafish, ly96 is expressed in a small number of macrophage-like cells. In a functional assay, zebrafish Md-2 and Tlr4ba form a complex that activates NF-κB signaling in response to LPS. In larval zebrafish ly96 loss-of-function mutations perturbed LPS-induced cytokine production but gave little protection against LPS toxicity. Finally, by analyzing the genomic context of tlr4 genes in 11 jawed vertebrates, we found that tlr4 arose prior to the divergence of teleosts and tetrapods. Thus, an LPS-sensitive Tlr4/Md-2 complex is likely an ancestral feature shared by mammals and zebrafish, rather than a de novo invention on the tetrapod lineage. We hypothesize that zebrafish retain an ancestral, low-sensitivity Tlr4/Md-2 complex that confers LPS responsiveness to a specific subset of innate immune cells.

Gouania willdenowi is a teleost fish without immunoglobulin genes

Immunoglobulin (Ig) genes encode antibodies in jawed vertebrates. They are essential elements of the adaptive immune response. Ig exists in soluble form or as part of the B cell membrane antigen receptor (BCR). Studies of Ig genes in fish genomes reveal the absence of Ig genes in Gouania willdenowi by deletion of the entire Ig locus from the canonical chromosomal region. The genes coding for integral BCR proteins, CD79a and CD79b, are also absent. Genes exist for T α/β lymphocyte receptors but not for the T γ/δ receptors. The results of the genomic analysis are independently corroborated with RNA-Seq transcriptomes from other Gobiesocidae species. From the transcriptome studies, Ig is also absent from these other Gobiesocidae species, Acyrtus sp. and Tomicodon sp. Present evidence suggests that Ig is missing from all species of the Gobiesocidae family.

Immunoglobulin T genes in Actinopterygii

In teleost fishes, there are three immunoglobulin isotypes named immunoglobulin M (IgM), D (IgD), and T (IgT). IgT was the last to be described in teleost fishes, and it is specific to them. From recent fish genomes, we identified and studied the immunoglobulin heavy chain genes in Actinopterygii. For this analysis, a custom bioinformatics and machine learning pipeline, we call CHfinder, was developed that identifies the exons coding for the CH domains of fish immunoglobulins. Some IgT in teleost and holostean fish found from this systematic study have not been previously described. Phylogenetic analysis of the deduced amino acid sequences of the IgT CH1 exons reveals they are similar to the CH1 of IgM. This analysis also shows that the other three domains (CH2, CH3, and CH4) were not the result of recent IgM duplication processes in Actinopterygii, demonstrating that it is an immunoglobulin of earlier origin. The bioinformatics program, CHfinder, is publicly available at https://github.com/compimmuno/CHfinder.

Exon probe sets and bioinformatics pipelines for all levels of fish phylogenomics

Exon markers have a long history of use in phylogenetics of ray-finned fishes, the most diverse clade of vertebrates with more than 35,000 species. As the number of published genomes increases, it has become easier to test exons and other genetic markers for signals of ancient duplication events and filter out paralogues that can mislead phylogenetic analysis. We present seven new probe sets for current target-capture phylogenomic protocols that capture 1,104 exons explicitly filtered for paralogues using gene trees. These seven probe sets span the diversity of teleost fishes, including four sets that target five hyperdiverse percomorph clades which together comprise ca. 17,000 species (Carangaria, Ovalentaria, Eupercaria, and Syngnatharia + Pelagiaria combined). We additionally included probes to capture legacy nuclear exons and mitochondrial markers that have been commonly used in fish phylogenetics (despite some exons being flagged for paralogues) to facilitate integration of old and new molecular phylogenetic matrices. We tested these probes experimentally for 56 fish species (eight species per probe set) and merged new exon-capture sequence data into an existing data matrix of 1,104 exons and 300 ray-finned fish species. We provide an optimized bioinformatics pipeline to assemble exon capture data from raw reads to alignments for downstream analysis. We show that legacy loci with known paralogues are at risk of assembling duplicated sequences with target-capture, but we also assembled many useful orthologous sequences that can be integrated with many PCR-generated matrices. These probe sets are a valuable resource for advancing fish phylogenomics because targeted exons can easily be extracted from increasingly available whole genome and transcriptome data sets, and also may be integrated with existing PCR-based exon and mitochondrial data.

Raw transcriptomics data to gene specific SSRs: a validated free bioinformatics workflow for biologists

Recent advances in next-generation sequencing technologies have paved the path for a considerable amount of sequencing data at a relatively low cost. This has revolutionized the genomics and transcriptomics studies. However, different challenges are now created in handling such data with available bioinformatics platforms both in assembly and downstream analysis performed in order to infer correct biological meaning. Though there are a handful of commercial software and tools for some of the procedures, cost of such tools has made them prohibitive for most research laboratories. While individual open-source or free software tools are available for most of the bioinformatics applications, those components usually operate standalone and are not combined for a user-friendly workflow. Therefore, beginners in bioinformatics might find analysis procedures starting from raw sequence data too complicated and time-consuming with the associated learning-curve. Here, we outline a procedure for de novo transcriptome assembly and Simple Sequence Repeats (SSR) primer design solely based on tools that are available online for free use. For validation of the developed workflow, we used Illumina HiSeq reads of different tissue samples of Santalum album (sandalwood), generated from a previous transcriptomics project. A portion of the designed primers were tested in the lab with relevant samples and all of them successfully amplified the targeted regions. The presented bioinformatics workflow can accurately assemble quality transcriptomes and develop gene specific SSRs. Beginner biologists and researchers in bioinformatics can easily utilize this workflow for research purposes.

Anemonefish, a model for Eco-Evo-Devo

Anemonefish, are a group of about 30 species of damselfish (Pomacentridae) that have long aroused the interest of coral reef fish ecologists. Combining a series of original biological traits and practical features in their breeding that are described in this paper, anemonefish are now emerging as an experimental system of interest for developmental biology, ecology and evolutionary sciences. They are small sized and relatively easy to breed in specific husbandries, unlike the large-sized marine fish used for aquaculture. Because they live in highly structured social groups in sea anemones, anemonefish allow addressing a series of relevant scientific questions such as the social control of growth and sex change, the mechanisms controlling symbiosis, the establishment and variation of complex color patterns, and the regulation of aging. Combined with the use of behavioral experiments, that can be performed in the lab or directly in the wild, as well as functional genetics and genomics, anemonefish provide an attractive experimental system for Eco-Evo-Devo.

Divergent Influenza-Like Viruses of Amphibians and Fish Support an Ancient Evolutionary Association

Influenza viruses (family Orthomyxoviridae) infect a variety of vertebrates, including birds, humans, and other mammals. Recent metatranscriptomic studies have uncovered divergent influenza viruses in amphibians, fish and jawless vertebrates, suggesting that these viruses may be widely distributed. We sought to identify additional vertebrate influenza-like viruses through the analysis of publicly available RNA sequencing data. Accordingly, by data mining, we identified the complete coding segments of five divergent vertebrate influenza-like viruses. Three fell as sister lineages to influenza B virus: salamander influenza-like virus in Mexican walking fish (Ambystoma mexicanum) and plateau tiger salamander (Ambystoma velasci), Siamese algae-eater influenza-like virus in Siamese algae-eater fish (Gyrinocheilus aymonieri) and chum salmon influenza-like virus in chum salmon (Oncorhynchus keta). Similarly, we identified two influenza-like viruses of amphibians that fell as sister lineages to influenza D virus: cane toad influenza-like virus and the ornate chorus frog influenza-like virus, in the cane toad (Rhinella marina) and ornate chorus frog (Microhyla fissipes), respectively. Despite their divergent phylogenetic positions, these viruses retained segment conservation and splicing consistent with transcriptional regulation in influenza B and influenza D viruses, and were detected in respiratory tissues. These data suggest that influenza viruses have been associated with vertebrates for their entire evolutionary history.

Initial data release and announcement of the 10,000 Fish Genomes Project (Fish10K)

Background

With more than 30,000 species, fish—including bony, jawless, and cartilaginous fish—are the largest vertebrate group, and include some of the earliest vertebrates. Despite their critical roles in many ecosystems and human society, fish genomics lags behind work on birds and mammals. This severely limits our understanding of evolution and hinders progress on the conservation and sustainable utilization of fish.

Results

Here, we announce the Fish10K project, a portion of the Earth BioGenome Project aiming to sequence 10,000 representative fish genomes in a systematic fashion within 10 years, and we officially welcome collaborators to join this effort. As a step towards this goal, we herein describe a feasible workflow for the procurement and storage of biospecimens, as well as sequencing and assembly strategies.

Conclusions

To illustrate, we present the genomes of 10 fish species from a cohort of 93 species chosen for technology development.

Liver transcriptome resources of four commercially exploited teleost species

The generation of omic resources is central to develop adequate management strategies for species with economic value. Here, we provide high-coverage RNA-seq datasets of liver tissue (containing between 80,2 and 88,4 million of paired-end reads) from four wildtype teleost species with high commercial value: Trachurus trachurus (TTR; Atlantic horse mackerel), Scomber scombrus (SSC; Atlantic mackerel), Trisopterus luscus (TLU; pout), and Micromesistius poutassou (MPO; blue whiting). A comprehensive assembly pipeline, using de novo single and multi-kmer assembly approaches, produced 64 single high-quality liver transcriptomes - 16 per species. The final assemblies, with N50 values ranging from 2,543-3,700 bp and BUSCO (Benchmarking Universal Single-Copy Orthologs) completeness values between 81.8-86.5% of the Actinopterygii gene set, were subjected to open reading frame (ORF) prediction and functional annotation. Our study provides the first transcriptomic resources for these species and offers valuable tools to evaluate both neutral and selected genetic variation among populations, and to identify candidate genes for environmental adaptation assisting in the investigation of the effects of global changes in fisheries.

Characterization of the prohormone complement in Amphiprion and related fish species integrating genome and transcriptome assemblies

The Amphiprion (anemonefish or clownfish) family of teleost fish, which is not a common model species, exhibits multiple unique characteristics, including social control of body size and protandrous sex change. The social changes in sex and body size are modulated by neuropeptide signaling pathways. These neuropeptides are formed from complex processing from larger prohormone proteins; understanding the neuropeptide complement requires information on complete prohormones sequences. Genome and transcriptome information within and across 22 teleost fish species, including 11 Amphiprion species, were assembled and integrated to achieve the first comprehensive survey of their prohormone genes. This information enabled the identification of 175 prohormone isoforms from 159 prohormone proteins across all species. This included identification of 9 CART prepropeptide genes and the loss of insulin-like 5B and tachykinin precursor 1B genes in Pomacentridae species. Transcriptome assemblies generally detected most prohormone genes but provided fewer prohormone genes than genome assemblies due to the lack of expression of prohormone genes or specific isoforms and tissue sampled. Comparisons between duplicate genes indicated that subfunctionalization, degradation, and neofunctionalization may be occurring between all copies. Characterization of the prohormone complement lays the foundation for future peptidomic investigation of the molecular basis of social physiology and behavior in the teleost fish.

Phylogenomics investigation of sparids (Teleostei: Spariformes) using high-quality proteomes highlights the importance of taxon sampling

Sparidae (Teleostei: Spariformes) are a family of fish constituted by approximately 150 species with high popularity and commercial value, such as porgies and seabreams. Although the phylogeny of this family has been investigated multiple times, its position among other teleost groups remains ambiguous. Most studies have used a single or few genes to decipher the phylogenetic relationships of sparids. Here, we conducted a thorough phylogenomic analysis using five recently available Sparidae gene-sets and 26 high-quality, genome-predicted teleost proteomes. Our analysis suggested that Tetraodontiformes (puffer fish, sunfish) are the closest relatives to sparids than all other groups used. By analytically comparing this result to our own previous contradicting finding, we show that this discordance is not due to different orthology assignment algorithms; on the contrary, we prove that it is caused by the increased taxon sampling of the present study, outlining the great importance of this aspect in phylogenomic analyses in general.

A Phylogenomic Supertree of Birds

It has long been appreciated that analyses of genomic data (e.g., whole genome sequencing or sequence capture) have the potential to reveal the tree of life, but it remains challenging to move from sequence data to a clear understanding of evolutionary history, in part due to the computational challenges of phylogenetic estimation using genome-scale data. Supertree methods solve that challenge because they facilitate a divide-and-conquer approach for large-scale phylogeny inference by integrating smaller subtrees in a computationally efficient manner. Here, we combined information from sequence capture and whole-genome phylogenies using supertree methods. However, the available phylogenomic trees had limited overlap so we used taxon-rich (but not phylogenomic) megaphylogenies to weave them together. This allowed us to construct a phylogenomic supertree, with support values, that included 707 bird species (~7% of avian species diversity). We estimated branch lengths using mitochondrial sequence data and we used these branch lengths to estimate divergence times. Our time-calibrated supertree supports radiation of all three major avian clades (Palaeognathae, Galloanseres, and Neoaves) near the Cretaceous-Paleogene (K-Pg) boundary. The approach we used will permit the continued addition of taxa to this supertree as new phylogenomic data are published, and it could be applied to other taxa as well.

Divergence, evolution and adaptation in ray-finned fish genomes

With the rapid development of next-generation sequencing technologies and bioinformatics, over 50 ray-finned fish genomes by far have been sequenced with high quality. The genomic work provides abundant genetic resources for deep understanding of divergence, evolution and adaptation in the fish genomes. They are also instructive for identification of candidate genes for functional verification, molecular breeding, and development of novel marine drugs. As an example of other omics data, the Fish-T1K project generated a big database of fish transcriptomes to integrate with these published fish genomes for potential applications. In this review, we highlight the above-mentioned recent investigations and core topics on the ray-finned fish genome research, with a main goal to obtain a deeper understanding of fish biology for theoretical and practical applications.

Transcriptome sequencing of the gill and barbel of Southern catfish (Silurus meridionalis) revealed immune responses and novel rhamnose-binding lectins (RBLs)

Southern catfish (Silurus meridionalis) is an economically important species widely cultured in China. It is well known for its fast growth, strong resistance to diseases and euryphage. However, little is known about the mechanisms for its powerful immune systems. Our Fish-T1K project has finished its first phase of 200 fish transcriptomes, with sequencing of gills in most examined fishes. In this study, we performed transcriptome sequencing of the gill and the maxillary barbel of Southern catfish, with the latter as a control. High expression of immune-related transcripts were observed in these two tissues. We observed that genes in the T cell receptor signaling pathway had higher transcription values in the gill than in the barbel. In addition, eight new rhamnose-binding lectins (RBLs) were identified and their carbohydrate recognition domains (CRDs) were classified according to the eight conserved cysteine residues and two conserved motifs (-YGR- and -DPC-). This is the first transcriptome report by high-throughput sequencing of the Southern catfish. Our genomic data and discovery of novel RBLs in this project should be able to promote better understandings of the roles of gills in immune responses and disease prevention for further aquaculture.

Whole Genome Sequencing of the Pirarucu (Arapaima gigas) Supports Independent Emergence of Major Teleost Clades

The Pirarucu (Arapaima gigas) is one of the world's largest freshwater fishes and member of the superorder Osteoglossomorpha (bonytongues), one of the oldest lineages of ray-finned fishes. This species is an obligate air-breather found in the basin of the Amazon River with an attractive potential for aquaculture. Its phylogenetic position among bony fishes makes the Pirarucu a relevant subject for evolutionary studies of early teleost diversification. Here, we present, for the first time, a draft genome version of the A. gigas genome, providing useful information for further functional and evolutionary studies. The A. gigas genome was assembled with 103-Gb raw reads sequenced in an Illumina platform. The final draft genome assembly was ∼661 Mb, with a contig N50 equal to 51.23 kb and scaffold N50 of 668 kb. Repeat sequences accounted for 21.69% of the whole genome, and a total of 24,655 protein-coding genes were predicted from the genome assembly, with an average of nine exons per gene. Phylogenomic analysis based on 24 fish species supported the postulation that Osteoglossomorpha and Elopomorpha (eels, tarpons, and bonefishes) are sister groups, both forming a sister lineage with respect to Clupeocephala (remaining teleosts). Divergence time estimations suggested that Osteoglossomorpha and Elopomorpha lineages emerged independently in a period of ∼30 Myr in the Jurassic. The draft genome of A. gigas provides a valuable genetic resource for further investigations of evolutionary studies and may also offer a valuable data for economic applications.

Comprehensive phylogeny of ray-finned fishes (Actinopterygii) based on transcriptomic and genomic data

Our understanding of phylogenetic relationships among bony fishes has been transformed by analysis of a small number of genes, but uncertainty remains around critical nodes. Genome-scale inferences so far have sampled a limited number of taxa and genes. Here we leveraged 144 genomes and 159 transcriptomes to investigate fish evolution with an unparalleled scale of data: >0.5 Mb from 1,105 orthologous exon sequences from 303 species, representing 66 out of 72 ray-finned fish orders. We apply phylogenetic tests designed to trace the effect of whole-genome duplication events on gene trees and find paralogy-free loci using a bioinformatics approach. Genome-wide data support the structure of the fish phylogeny, and hypothesis-testing procedures appropriate for phylogenomic datasets using explicit gene genealogy interrogation settle some long-standing uncertainties, such as the branching order at the base of the teleosts and among early euteleosts, and the sister lineage to the acanthomorph and percomorph radiations. Comprehensive fossil calibrations date the origin of all major fish lineages before the end of the Cretaceous.

A resource for sustainable management: De novo assembly and annotation of the liver transcriptome of the Atlantic chub mackerel, Scomber colias

Mackerels represent a valuable fishery worldwide. Their ample geographic distribution and capture levels make them an insightful model to address stock management strategies in the context of global changes. Yet, and despite recent impressive genome and transcriptome sequencing efforts from teleost species, available resources from the Scombridae family are comparatively scarce. Here, we generated the first high-quality de novo assembly of the liver transcriptome of the Atlantic chub mackerel (Scomber colias). Through the use of RNA-Seq Illumina technology, 111,124,228 clean reads were obtained for the liver transcriptome. De novo assembly resulted in 93,731 transcripts with an N50 of 1462 bp. This dataset provides an important insight into the context of fisheries management.

Characterization of viral RNA splicing using whole-transcriptome datasets from host species

RNA alternative splicing (AS) is an important post-transcriptional mechanism enabling single genes to produce multiple proteins. It has been well demonstrated that viruses deploy host AS machinery for viral protein productions. However, knowledge on viral AS is limited to a few disease-causing viruses in model species. Here we report a novel approach to characterizing viral AS using whole transcriptome dataset from host species. Two insect transcriptomes (Acheta domesticus and Planococcus citri) generated in the 1,000 Insect Transcriptome Evolution (1KITE) project were used as a proof of concept using the new pipeline. Two closely related densoviruses (Acheta domesticus densovirus, AdDNV, and Planococcus citri densovirus, PcDNV, Ambidensovirus, Densovirinae, Parvoviridae) were detected and analyzed for AS patterns. The results suggested that although the two viruses shared major AS features, dramatic AS divergences were observed. Detailed analysis of the splicing junctions showed clusters of AS events occurred in two regions of the virus genome, demonstrating that transcriptome analysis could gain valuable insights into viral splicing. When applied to large-scale transcriptomics projects with diverse taxonomic sampling, our new method is expected to rapidly expand our knowledge on RNA splicing mechanisms for a wide range of viruses.

Identifying Sequenced Eukaryotic Genomes and Transcriptomes with diArk

The diArk Eukaryotic Genome Database is a manually curated and updated repository of available eukaryotic genome and transcriptome assemblies. diArk is a key resource for researchers interested in comparative eukaryotic genomics, and the entry point to browsing sequenced eukaryotes in general and to find the most closely related species to the own organism of interest in particular. The exponentially increasing number of sequenced species demands sophisticated search and data presentation tools. In this chapter we describe how to navigate the diArk database keeping a first-time user in mind.

A Comparative Genomic and Transcriptomic Survey Provides Novel Insights into N-Acetylserotonin Methyltransferase (ASMT) in Fish

Melatonin is a multifunctional bioactive molecule that plays comprehensive physiological roles in all living organisms. N-acetylserotonin methyltransferase (ASMT, also known as hydroxyindole O-methyltransferase or HIOMT) is the final enzyme for biosynthesis of melatonin. Here, we performed a comparative genomic and transcriptomic survey to explore the ASMT family in fish. Two ASMT isotypes (ASMT1 and ASMT2) and a new ASMT-like (ASMTL) are all extracted from teleost genomes on the basis of phylogenetic and synteny analyses. We confirmed that C-terminal of the ASMTL proteins (ASMTL-ASMT) is homology to the full length of ASMT1 and ASMT2. Our results also demonstrate that the two ASMT isotypes and their distribution in teleosts seem to be the result of combinations of whole-genome duplication (WGD) and gene loss. Differences were also observed in tissue distribution and relative transcript abundances of ASMT1, ASMT2 and ASMTL through transcriptomic analysis. Protein sequence alignment and 3D structure prediction of ASMTs and ASMTL suggest differential roles for these ASMT genes. In summary, our current work provides novel insights into the ASMT genes in fish by combination of genomic and transcriptomic data.

Phylogenetic classification of bony fishes

BACKGROUND

Fish classifications, as those of most other taxonomic groups, are being transformed drastically as new molecular phylogenies provide support for natural groups that were unanticipated by previous studies. A brief review of the main criteria used by ichthyologists to define their classifications during the last 50 years, however, reveals slow progress towards using an explicit phylogenetic framework. Instead, the trend has been to rely, in varying degrees, on deep-rooted anatomical concepts and authority, often mixing taxa with explicit phylogenetic support with arbitrary groupings. Two leading sources in ichthyology frequently used for fish classifications (JS Nelson's volumes of Fishes of the World and W. Eschmeyer's Catalog of Fishes) fail to adopt a global phylogenetic framework despite much recent progress made towards the resolution of the fish Tree of Life. The first explicit phylogenetic classification of bony fishes was published in 2013, based on a comprehensive molecular phylogeny ( www.deepfin.org ). We here update the first version of that classification by incorporating the most recent phylogenetic results.

RESULTS

The updated classification presented here is based on phylogenies inferred using molecular and genomic data for nearly 2000 fishes. A total of 72 orders (and 79 suborders) are recognized in this version, compared with 66 orders in version 1. The phylogeny resolves placement of 410 families, or ~80% of the total of 514 families of bony fishes currently recognized. The ordinal status of 30 percomorph families included in this study, however, remains uncertain (incertae sedis in the series Carangaria, Ovalentaria, or Eupercaria). Comments to support taxonomic decisions and comparisons with conflicting taxonomic groups proposed by others are presented. We also highlight cases were morphological support exist for the groups being classified.

CONCLUSIONS

This version of the phylogenetic classification of bony fishes is substantially improved, providing resolution for more taxa than previous versions, based on more densely sampled phylogenetic trees. The classification presented in this study represents, unlike any other, the most up-to-date hypothesis of the Tree of Life of fishes.

From Marine Venoms to Drugs: Efficiently Supported by a Combination of Transcriptomics and Proteomics

The potential of marine natural products to become new drugs is vast; however, research is still in its infancy. The chemical and biological diversity of marine toxins is immeasurable and as such an extraordinary resource for the discovery of new drugs. With the rapid development of next-generation sequencing (NGS) and liquid chromatography-tandem mass spectrometry (LC-MS/MS), it has been much easier and faster to identify more toxins and predict their functions with bioinformatics pipelines, which pave the way for novel drug developments. Here we provide an overview of related bioinformatics pipelines that have been supported by a combination of transcriptomics and proteomics for identification and function prediction of novel marine toxins.

Targeted sequencing for high-resolution evolutionary analyses following genome duplication in salmonid fish: Proof of concept for key components of the insulin-like growth factor axis

High-throughput sequencing has revolutionised comparative and evolutionary genome biology. It has now become relatively commonplace to generate multiple genomes and/or transcriptomes to characterize the evolution of large taxonomic groups of interest. Nevertheless, such efforts may be unsuited to some research questions or remain beyond the scope of some research groups. Here we show that targeted high-throughput sequencing offers a viable alternative to study genome evolution across a vertebrate family of great scientific interest. Specifically, we exploited sequence capture and Illumina sequencing to characterize the evolution of key components from the insulin-like growth (IGF) signalling axis of salmonid fish at unprecedented phylogenetic resolution. The IGF axis represents a central governor of vertebrate growth and its core components were expanded by whole genome duplication in the salmonid ancestor ~95Ma. Using RNA baits synthesised to genes encoding the complete family of IGF binding proteins (IGFBP) and an IGF hormone (IGF2), we captured, sequenced and assembled orthologous and paralogous exons from species representing all ten salmonid genera. This approach generated 299 novel sequences, most as complete or near-complete protein-coding sequences. Phylogenetic analyses confirmed congruent evolutionary histories for all nineteen recognized salmonid IGFBP family members and identified novel salmonid-specific IGF2 paralogues. Moreover, we reconstructed the evolution of duplicated IGF axis paralogues across a replete salmonid phylogeny, revealing complex historic selection regimes - both ancestral to salmonids and lineage-restricted - that frequently involved asymmetric paralogue divergence under positive and/or relaxed purifying selection. Our findings add to an emerging literature highlighting diverse applications for targeted sequencing in comparative-evolutionary genomics. We also set out a viable approach to obtain large sets of nuclear genes for any member of the salmonid family, which should enable insights into the evolutionary role of whole genome duplication before additional nuclear genome sequences become available.