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

Phylogenomic incongruence, hypothesis testing, and taxonomic sampling: The monophyly of characiform fishes

Phylogenomic studies using genome-wide datasets are quickly becoming the state of the art for systematics and comparative studies, but in many cases, they result in strongly supported incongruent results. The extent to which this conflict is real depends on different sources of error potentially affecting big datasets (assembly, stochastic, and systematic error). Here, we apply a recently developed methodology (GGI or gene genealogy interrogation) and data curation to new and published datasets with more than 1000 exons, 500 ultraconserved element (UCE) loci, and transcriptomic sequences that support incongruent hypotheses. The contentious non-monophyly of the order Characiformes proposed by two studies is shown to be a spurious outcome induced by sample contamination in the transcriptomic dataset and an ambiguous result due to poor taxonomic sampling in the UCE dataset. By exploring the effects of number of taxa and loci used for analysis, we show that the power of GGI to discriminate among competing hypotheses is diminished by limited taxonomic sampling, but not equally sensitive to gene sampling. Taken together, our results reinforce the notion that merely increasing the number of genetic loci for a few representative taxa is not a robust strategy to advance phylogenetic knowledge of recalcitrant groups. We leverage the expanded exon capture dataset generated here for Characiformes (206 species in 23 out of 24 families) to produce a comprehensive phylogeny and a revised classification of the order.

Genome Assembly for a Yunnan-Guizhou Plateau "3E" Fish, Anabarilius grahami (Regan), and Its Evolutionary and Genetic Applications

A Yunnan-Guizhou Plateau fish, the Kanglang white minnow (Anabarilius grahami), is a typical “3E” (Endangered, Endemic, and Economic) species in China. Its distribution is limited to Fuxian Lake, the nation’s second deepest lake, with a significant local economic value but a drastically declining wild population. This species has been evaluated as VU (Vulnerable) in the China Species Red List. As one of the “Four Famous Fish” in Yunnan province, the artificial breeding has been achieved since 2003. It has not only re-established its wild natural populations by reintroduction of the artificial breeding stocks, but also brought a wide and popular utilization of this species to the local fish farms. A. grahami has become one of the main native aquaculture species in Yunnan province, and the artificial production has been emerging in steady growth each year. To promote the conservation and sustainable utilization of this fish, we initiated its whole genome sequencing project using an Illumina Hiseq2500 platform. The assembled genome size of A. grahami is 1.006 Gb, accounting for 98.63% of the estimated genome size (1.020 Gb), with contig N50 and scaffold N50 values of 26.4 kb and 4.41 Mb, respectively. Approximately about 50.38% of the genome was repetitive. A total of 25,520 protein-coding genes were subsequently predicted. A phylogenetic tree based on 4,580 single-copy genes from A. grahami and 18 other cyprinids revealed three well-supported subclades within the Cyprinidae. This is the first inter-subfamily relationship of cyprinids at genome level, providing a simple yet useful framework for understanding the traditional but popular subfamily classification systems. Interestingly, a further population demography of A. grahami uncovered a historical relationship between this fish and Fuxian Lake, suggesting that range expansion or shrinkage of the habitat has had a remarkable impact on the population size of endemic plateau fishes. Additionally, a total of 33,836 simple sequence repeats (SSR) markers were identified, and 11 loci were evaluated for a preliminary genetic diversity analysis in this study, thus providing another useful genetic resource for studying this “3E” species.

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.

A phylogenomic approach to reconstruct interrelationships of main clupeocephalan lineages with a critical discussion of morphological apomorphies

BACKGROUND

Previous molecular studies on the phylogeny and classification of clupeocephalan fishes revealed numerous new taxonomic entities. For re-analysing these taxa, we perform target gene capturing and subsequent next generation sequencing of putative ortholog exons of major clupeocephalan lineages. Sequence information for the RNA bait design was derived from publicly available genomes of bony fishes. Newly acquired sequence data comprising > 800 exon sequences was subsequently used for phylogenetic reconstructions.

RESULTS

Our results support monophyletic Otomorpha comprising Alepocephaliformes. Within Ostariophysi, Gonorynchiformes are sister to a clade comprising Cypriniformes, Characiformes, Siluriformes and Gymnotiformes, where the interrelationships of Characiformes, Siluriformes and Gymnotiformes remain enigmatic. Euteleosts comprise four major clades: Lepidogalaxiiformes, Protacanthopterygii, Stomiatii, and Galaxiiformes plus Neoteleostei. The monotypic Lepidogalaxiiformes form the sister-group to all remaining euteleosts. Protacanthopterygii, comprising Argentini-, Esoci- and Salmoniformes, is sister to Stomiatii (Osmeriformes and Stomiatiformes) and Galaxiiformes plus Neoteleostei.

CONCLUSIONS

Several proposed monophyla defined by morphological apomorphies within the Clupeocephalan phylogeny are confirmed by the phylogenetic estimates presented herein. However, other morphologically described groups cannot be reconciled with molecular phylogenies. Thus, numerous morphological apomoprhies of supposed monophyla are called into question. The interpretation of suggested morphological synapomorphies of otomorph fishes is strongly affected by the inclusion of deep-sea inhabiting, and to that effect morphologically adapted Alepocephaliformes. Our revision of these potential synapomorphies, in the context that Alepocephaliformes are otomorph fishes, reveals that only a single character of the total nine characters proposed as synapomorphic for the group is clearly valid for all otomorphs. Three further characters remain possible apomorphies since their status remains unclear in the deep-sea adapted Alepocephaliformes showing developmental lag and lacking a swim bladder. Further, our analysis places Galaxiiformes as sister group to neoteleosts, which contradicts some previous molecular phylogenetic studies. This needs further investigation from a morphological perspective, as suggested synapomophies for several euteleostean lineages are challenged or still lacking. For the verification of results presented herein, a denser phylogenomic-level taxon sampling should be applied.

Retention of fatty acyl desaturase 1 (fads1) in Elopomorpha and Cyclostomata provides novel insights into the evolution of long-chain polyunsaturated fatty acid biosynthesis in vertebrates

Background

Provision of long-chain polyunsaturated fatty acids (LC-PUFA) in vertebrates occurs through the diet or via endogenous production from C₁₈ precursors through consecutive elongations and desaturations. It has been postulated that the abundance of LC-PUFA in the marine environment has remarkably modulated the gene complement and function of Fads in marine teleosts. In vertebrates two fatty acyl desaturases, namely Fads1 and Fads2, encode ∆5 and ∆6 desaturases, respectively. To fully clarify the evolutionary history of LC-PUFA biosynthesis in vertebrates, we investigated the gene repertoire and function of Fads from species placed at key evolutionary nodes.

Results

We demonstrate that functional Fads1Δ5 and Fads2∆6 arose from a tandem gene duplication in the ancestor of vertebrates, since they are present in the Arctic lamprey. Additionally, we show that a similar condition was retained in ray-finned fish such as the Senegal bichir and spotted gar, with the identification of fads1 genes in these lineages. Functional characterisation of the isolated desaturases reveals the first case of a Fads1 enzyme with ∆5 desaturase activity in the Teleostei lineage, the Elopomorpha. In contrast, in Osteoglossomorpha genomes, while no fads1 was identified, two separate fads2 duplicates with ∆6 and ∆5 desaturase activities respectively were uncovered.

Conclusions

We conclude that, while the essential genetic components involved LC-PUFA biosynthesis evolved in the vertebrate ancestor, the full completion of the LC-PUFA biosynthesis pathway arose uniquely in gnathostomes.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1271-5) contains supplementary material, which is available to authorized users.

A review of the systematic biology of fossil and living bony-tongue fishes, Osteoglossomorpha (Actinopterygii: Teleostei)

ABSTRACT The bony-tongue fishes, Osteoglossomorpha, have been the focus of a great deal of morphological, systematic, and evolutionary study, due in part to their basal position among extant teleostean fishes. This group includes the mooneyes (Hiodontidae), knifefishes (Notopteridae), the abu (Gymnarchidae), elephantfishes (Mormyridae), arawanas and pirarucu (Osteoglossidae), and the African butterfly fish (Pantodontidae). This morphologically heterogeneous group also has a long and diverse fossil record, including taxa from all continents and both freshwater and marine deposits. The phylogenetic relationships among most extant osteoglossomorph families are widely agreed upon. However, there is still much to discover about the systematic biology of these fishes, particularly with regard to the phylogenetic affinities of several fossil taxa, within Mormyridae, and the position of Pantodon. In this paper we review the state of knowledge for osteoglossomorph fishes. We first provide an overview of the diversity of Osteoglossomorpha, and then discuss studies of the phylogeny of Osteoglossomorpha from both morphological and molecular perspectives, as well as biogeographic analyses of the group. Finally, we offer our perspectives on future needs for research on the systematic biology of Osteoglossomorpha.

"Out of the Can": A Draft Genome Assembly, Liver Transcriptome, and Nutrigenomics of the European Sardine, Sardina pilchardus

Clupeiformes, such as sardines and herrings, represent an important share of worldwide fisheries. Among those, the European sardine (Sardina pilchardus, Walbaum 1792) exhibits significant commercial relevance. While the last decade showed a steady and sharp decline in capture levels, recent advances in culture husbandry represent promising research avenues. Yet, the complete absence of genomic resources from sardine imposes a severe bottleneck to understand its physiological and ecological requirements. We generated 69 Gbp of paired-end reads using Illumina HiSeq X Ten and assembled a draft genome assembly with an N50 scaffold length of 25,579 bp and BUSCO completeness of 82.1% (Actinopterygii). The estimated size of the genome ranges between 655 and 850 Mb. Additionally, we generated a relatively high-level liver transcriptome. To deliver a proof of principle of the value of this dataset, we established the presence and function of enzymes (Elovl2, Elovl5, and Fads2) that have pivotal roles in the biosynthesis of long chain polyunsaturated fatty acids, essential nutrients particularly abundant in oily fish such as sardines. Our study provides the first omics dataset from a valuable economic marine teleost species, the European sardine, representing an essential resource for their effective conservation, management, and sustainable exploitation.

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.

Phylogeny of the hillstream loach genus Mesonoemacheilus reveals widespread diversification through ancient drainage connections in the Western Ghats Biodiversity Hotspot

Rivers draining the Western Ghats (WG) mountain ranges in peninsular India harbor an exceptionally diverse, unique and evolutionarily distinct assemblage of lower vertebrates with high levels of endemism, attributed to their evolution and potentially long history of isolation during the Late Cretaceous or Early Tertiary. A molecular phylogeny of hillstream loaches of the genus Mesonoemacheilus endemic to the WG revealed the presence of four clades which we designate as 'species groups'. A consensus of various species delimitation methods indicates the likelihood of 'at least' seven more undescribed species within Mesonoemacheilus. Molecular clock analysis dates the basal clade around 38 mya in the Paleogene, and subsequent diversification in the Neogene and Quaternary periods resulting in the current genetic diversity. Biogeographic analysis suggests that vicariance events which separated the rivers on either side of the two geological barriers/gaps, the Palghat and Shencottah, in the Neogene, as well as range contractions and cladogenetic events contributed to the current patterns of diversity and distribution of this genus. Our results also provide preliminary indications on the interconnections and faunal exchange between historical river drainages in the WG region.