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

The osteology of the golden grey mullet Liza aurata (Teleostei: Mugiliformes: Mugilidae) including interactive three-dimensional reconstructions

Grey mullets are remarkably characterized by their overall uniform external morphology. Identifying species as well as positioning the Mugiliformes in a phylogenetic context is rather difficult. Most recently they were placed in the newly erected Ovalentaria, but more detailed relationships to potential sister taxa were not resolved. Studying the internal morphology, especially the osteology, might provide new insights into the evolution of the Mugiliformes as well as help clarify their phylogenetic position within the Ovalentaria. A detailed osteology of the golden grey mullet Liza aurata is presented. The use of cleared and stained specimens allowed for a complete examination of bony and cartilaginous structures, and a 3D reconstruction from a μCT data set provided additional information on the positional relationships of the bones. Following this, the data obtained were compared with different mugilid species, particularly with the flathead grey mullet Mugil cephalus. Several differences between these species could be identified, such as the position of the basisphenoid, the shape of the hyomandibular and the composition of the branchial arches. These characters might help in understanding the evolutionary changes happening within the mugiliforms and will provide the basis to study this taxon in detail, finally allowing the reconstruction of the body plan of grey mullets.

Extreme genomic volatility characterizes the evolution of the immunoglobulin heavy chain locus in cyprinodontiform fishes

The evolution of the adaptive immune system has provided vertebrates with a uniquely sophisticated immune toolkit, enabling them to mount precise immune responses against a staggeringly diverse range of antigens. Like other vertebrates, teleost fishes possess a complex and functional adaptive immune system; however, our knowledge of the complex antigen-receptor genes underlying its functionality has been restricted to a small number of experimental and agricultural species, preventing systematic investigation into how these crucial gene loci evolve. Here, we analyse the genomic structure of the immunoglobulin heavy chain (IGH) gene loci in the cyprinodontiforms, a diverse and important group of teleosts present in many different habitats across the world. We reconstruct the complete IGH loci of the turquoise killifish (Nothobranchius furzeri) and the southern platyfish (Xiphophorus maculatus) and analyse their in vivo gene expression, revealing the presence of species-specific splice isoforms of transmembrane IGHM. We further characterize the IGH constant regions of 10 additional cyprinodontiform species, including guppy, Amazon molly, mummichog and mangrove killifish. Phylogenetic analysis of these constant regions suggests multiple independent rounds of duplication and deletion of the teleost-specific antibody class IGHZ in the cyprinodontiform lineage, demonstrating the extreme volatility of IGH evolution. Focusing on the cyprinodontiforms as a model taxon for comparative evolutionary immunology, this work provides novel genomic resources for studying adaptive immunity and sheds light on the evolutionary history of the adaptive immune system.

The genome of Mekong tiger perch (Datnioides undecimradiatus) provides insights into the phylogenetic position of Lobotiformes and biological conservation

Mekong tiger perch (Datnioides undecimradiatus) is an ornamental and vulnerable freshwater fish native to the Mekong basin in Indochina, belonging to the order Lobotiformes. Here, we generated 121X stLFR co-barcode clean reads and 18X Oxford Nanopore MinION reads and obtained a 595 Mb Mekong tiger perch genome, which is the first whole genome sequence in the order Lobotiformes. Based on this genome, the phylogenetic tree analysis suggested that Lobotiformes is more closely related to Sciaenidae than to Tetraodontiformes, resolving a long-time dispute. We depicted the genes involved in pigment development in Mekong tiger perch and results confirmed that the four rate-limiting genes of pigment synthesis had been retained after fish-specific genome duplication. We also estimated the demographic history of Mekong tiger perch, which showed that the effective population size suffered a continuous reduction possibly related to the contraction of immune-related genes. Our study provided a reference genome resource for the Lobotiformes, as well as insights into the phylogenetic position of Lobotiformes and biological conservation.

Basal teleosts provide new insights into the evolutionary history of teleost-duplicated aromatase

Duplicated cyp19a1 genes (cyp19a1a encoding aromatase a and cyp19a1b encoding aromatase b) have been identified in an increasing number of teleost species. Cyp19a1a is mainly expressed in the gonads, while cyp19a1b is mainly expressed in the brain, specifically in radial glial cells, as largely investigated by Kah and collaborators. The third round of whole-genome duplication that specifically occurred in the teleost lineage (TWGD or 3R) is likely at the origin of the duplicated cyp19a1 paralogs. In contrast to the situation in other teleosts, our previous studies identified a single cyp19a1 in eels (Anguilla), which are representative species of a basal group of teleosts, Elopomorpha. In the present study, using genome data mining and phylogenetic and synteny analyses, we confirmed that the whole aromatase genomic region was duplicated in eels, with most aromatase-neighboring genes being conserved in duplicate in eels, as in other teleosts. These findings suggest that specific gene loss of one of the 3R-duplicated cyp19a1 paralogs occurred in Elopomorpha after TWGD. Similarly, a single cyp19a1 gene was found in the arowana, which is a representative species of another basal group of teleosts, Osteoglossomorpha. In eels, the single cyp19a1 is expressed in both the brain and the gonads, as observed for the single CYP19A1 gene present in other vertebrates. The results of phylogenetic, synteny, closest neighboring gene, and promoter structure analyses showed that the single cyp19a1 of the basal teleosts shared conserved properties with both teleost cyp19a1a and cyp19a1b paralogs, which did not allow us to conclude which of the 3R-duplicated paralogs (cyp19a1a or cyp19a1b) was lost in Elopomorpha. Elopomorpha and Osteoglossomorpha cyp19a1 genes exhibited preserved ancestral functions, including expression in both the gonad and brain. We propose that the subfunctionalization of the 3R-duplicated cyp19a1 paralogs expressed specifically in the gonad or brain occurred in Clupeocephala, after the split of Clupeocephala from Elopomorpha and Osteoglossomorpha, which represented a driving force for the conservation of both 3R-duplicated paralogs in all extant Clupeocephala. In contrast, the functional redundancy of the undifferentiated 3R-duplicated cyp19a1 paralogs in elopomorphs and osteoglossomorphs would have favored the loss of one 3R paralog in basal teleosts.

Linking Virus Discovery to Immune Responses Visualized during Zebrafish Infections

The discovery of new viruses currently outpaces our capacity for experimental examination of infection biology. To better couple virus discovery with immunology, we genetically modified zebrafish to visually report on virus infections. After generating a strain that expresses green fluorescent protein (GFP) under an interferon-stimulated gene promoter, we repeatedly observed transgenic larvae spontaneously expressing GFP days after hatching. RNA sequencing comparisons of co-housed GFP-positive and GFP-negative zebrafish revealed a naturally occurring picornavirus that induced a canonical interferon-mediated response and hundreds of antiviral defense genes not observed following immunostimulatory treatments or experimental infections with other viruses. Among the many genes induced by picornavirus infection was a large set encoding guanosine triphosphatase (GTPase) of immunity-associated proteins (GIMAPs). The GIMAP gene family is massively expanded in fish genomes and may also play a crucial role in antiviral responses in mammals, including humans. We subsequently detected zebrafish picornavirus in publicly available sequencing data from seemingly asymptomatic zebrafish in many research institutes and found that it altered gene expression in a previous study of zebrafish development. Experiments revealed a horizontal mode of virus transmission, highlighting a system for studying the spread of picornavirus infections within and between individuals. Our study describes a naturally occurring picornavirus that elicits strong antiviral responses in zebrafish and provides new strategies for simultaneously discovering viruses and their impact on vertebrate hosts.

Facial and gill musculature of polynemid fishes, with notes on their possible relationships with sciaenids (Percomorphacea: Perciformes)

The Polynemidae is a family of primarily marine fishes with eight genera and 42 extant species. Many aspects of their morphology are largely unknown, with few reports about their osteology and barely any information on their myology. This paper describes and illustrates in detail all facial and branchial muscles of representative species of polynemids. Our analysis demonstrates the existence of several remarkable and previously unknown specializations in the polynemid musculature. The aponeurotic and completely independent origin of the pars promalaris of the adductor mandibulae is apparently unique among percomorphs. The differentiation of this section into lateral and medial subsections; the total separation of the promalaris from the retromalaris; the differentiation of the pars primordialis of the levator arcus palatini into external and internal subsections are also uncommon features of polynemids that are shared by sciaenids, thus supporting the hypothesis of a closer relationship between these families.

Large-bodied sabre-toothed anchovies reveal unanticipated ecological diversity in early Palaeogene teleosts

Many modern groups of marine fishes first appear in the fossil record during the early Palaeogene (66-40 Ma), including iconic predatory lineages of spiny-rayed fishes that appear to have originated in response to ecological roles left empty after the Cretaceous/Palaeogene extinction. The hypothesis of extinction-mediated ecological release likewise predicts that other fish groups have adopted novel predatory ecologies. Here, we report remarkable trophic innovation in early Palaeogene clupeiforms (herrings and allies), a group whose modern representatives are generally small-bodied planktivores. Two forms, the early Eocene (Ypresian) †Clupeopsis from Belgium and a new genus from the middle Eocene (Lutetian) of Pakistan, bear conspicuous features indicative of predatory ecology, including large size, long gapes and caniniform dentition. Most remarkable is the presence of a single, massive vomerine fang offset from the midline in both. Numerous features of the neurocranium, suspensorium and branchial skeleton place these taxa on the engraulid (anchovy) stem as the earliest known representatives of the clade. The identification of large-bodied, piscivorous anchovies contributes to an emerging picture of a phylogenetically diverse guild of predatory ray-finned fishes in early Palaeogene marine settings, which include completely extinct lineages alongside members of modern marine groups and taxa that are today restricted to freshwater or deep-sea environments.

Between an ocean and a high place: coastal drainage isolation generates endemic cryptic species in the Cape kurper Sandelia capensis (Anabantiformes: Anabantidae), Cape Region, South Africa

This study investigated the range-wide phylogenetics and biogeography of the Cape kurper Sandelia capensis, a primary freshwater fish endemic to and widespread within the Cape Floristic Region (CFR) of South Africa. Maximum likelihood, Bayesian phylogenetic and haplotype network analyses, based on two mitochondrial and two nuclear genes, revealed the existence of three reciprocally monophyletic, deeply divergent and allopatric clades that probably represent cryptic species. The West Coast Clade is largely confined to the Langvlei, Verlorenvlei, Berg and Diep Rivers, the Klein River Clade is endemic to the Klein River and the South Coast Clade is found everywhere else in the range of S. capensis sensu lato. It was hypothesised that divergences within S. capensis sensu lato probably occurred because of isolation of coastal drainages by persistent drainage divides or vicariance of current tributaries by the drowning of their confluences by high sea levels. The current distribution of lineages could be due to historical range expansion and gene flow via river capture or some other mode of transdivide dispersal or dispersal during periods of low sea level via palaeoriver confluences of currently isolated coastal rivers. Comparison of BEAST2 estimated divergence times with the timing of climatic, geological and geomorphological events supported long-term coastal drainage isolation, punctuated by rare transdivide dispersal events and limited palaeoriver dispersal, as the best explanation of current phylogeographic and divergence patterns in S. capensis. Hydrological barriers that block upstream passage in palaeotributaries could hypothetically explain why S. capensis failed to disperse through certain palaeoriver confluences. There were several sites where biogeographic patterns have likely been confounded by human translocation of S. capensis. Alien fish predators and water extraction may threaten the three cryptic species more severely than previously realised, due to their smaller population sizes and inhabitation of only a portion of the range previously ascribed to S. capensis sensu lato. The preponderance of cryptic diversity and endemism in the CFR suggests that additional undescribed cryptic species of obligate freshwater fishes may be found in short coastal river systems around the world, especially in regions with a history of geological stability and a narrow continental shelf.

Seeing the rainbow: mechanisms underlying spectral sensitivity in teleost fishes

Among vertebrates, teleost eye diversity exceeds that found in all other groups. Their spectral sensitivities range from ultraviolet to red, and the number of visual pigments varies from 1 to over 40. This variation is correlated with the different ecologies and life histories of fish species, including their variable aquatic habitats: murky lakes, clear oceans, deep seas and turbulent rivers. These ecotopes often change with the season, but fish may also migrate between ecotopes diurnally, seasonally or ontogenetically. To survive in these variable light habitats, fish visual systems have evolved a suite of mechanisms that modulate spectral sensitivities on a range of timescales. These mechanisms include: (1) optical media that filter light, (2) variations in photoreceptor type and size to vary absorbance and sensitivity, and (3) changes in photoreceptor visual pigments to optimize peak sensitivity. The visual pigment changes can result from changes in chromophore or changes to the opsin. Opsin variation results from changes in opsin sequence, opsin expression or co-expression, and opsin gene duplications and losses. Here, we review visual diversity in a number of teleost groups where the structural and molecular mechanisms underlying their spectral sensitivities have been relatively well determined. Although we document considerable variability, this alone does not imply functional difference per se. We therefore highlight the need for more studies that examine species with known sensitivity differences, emphasizing behavioral experiments to test whether such differences actually matter in the execution of visual tasks that are relevant to the fish.

Harnessing genomics to fast-track genetic improvement in aquaculture

Aquaculture is the fastest-growing farmed food sector and will soon become the primary source of fish and shellfish for human diets. In contrast to crop and livestock production, aquaculture production is derived from numerous, exceptionally diverse species that are typically in the early stages of domestication. Genetic improvement of production traits via well-designed, managed breeding programmes has great potential to help meet the rising seafood demand driven by human population growth. Supported by continuous advances in sequencing and bioinformatics, genomics is increasingly being applied across the broad range of aquaculture species and at all stages of the domestication process to optimize selective breeding. In the future, combining genomic selection with biotechnological innovations, such as genome editing and surrogate broodstock technologies, may further expedite genetic improvement in aquaculture.

Comprehensive Transcriptome Analysis Reveals Insights into Phylogeny and Positively Selected Genes of Sillago Species

Sillago species lives in the demersal environments and face multiple stressors, such as localized oxygen depletion, sulfide accumulation, and high turbidity. In this study, we performed transcriptome analyses of seven Sillago species to provide insights into the phylogeny and positively selected genes of this species. After de novo assembly, 82,024, 58,102, 63,807, 85,990, 102,185, 69,748, and 102,903 unigenes were generated from S. japonica, S. aeolus, S. sp.1, S. sihama, S. sp.2, S. parvisquamis, and S. sinica, respectively. Furthermore, 140 shared orthologous exon markers were identified and then applied to reconstruct the phylogenetic relationships of the seven Sillago species. The reconstructed phylogenetic structure was significantly congruent with the prevailing morphological and molecular biological view of Sillago species relationships. In addition, a total of 44 genes were identified to be positively selected, and these genes were potential participants in the stress response, material (carbohydrate, amino acid and lipid) and energy metabolism, growth and differentiation, embryogenesis, visual sense, and other biological processes. We suspected that these genes possibly allowed Sillago species to increase their ecological adaptation to multiple environmental stressors.

Transcriptome Ortholog Alignment Sequence Tools (TOAST) for phylogenomic dataset assembly

Background

Advances in next-generation sequencing technologies have reduced the cost of whole transcriptome analyses, allowing characterization of non-model species at unprecedented levels. The rapid pace of transcriptomic sequencing has driven the public accumulation of a wealth of data for phylogenomic analyses, however lack of tools aimed towards phylogeneticists to efficiently identify orthologous sequences currently hinders effective harnessing of this resource.

Results

We introduce TOAST, an open source R software package that can utilize the ortholog searches based on the software Benchmarking Universal Single-Copy Orthologs (BUSCO) to assemble multiple sequence alignments of orthologous loci from transcriptomes for any group of organisms. By streamlining search, query, and alignment, TOAST automates the generation of locus and concatenated alignments, and also presents a series of outputs from which users can not only explore missing data patterns across their alignments, but also reassemble alignments based on user-defined acceptable missing data levels for a given research question.

Conclusions

TOAST provides a comprehensive set of tools for assembly of sequence alignments of orthologs for comparative transcriptomic and phylogenomic studies. This software empowers easy assembly of public and novel sequences for any target database of candidate orthologs, and fills a critically needed niche for tools that enable quantification and testing of the impact of missing data. As open-source software, TOAST is fully customizable for integration into existing or novel custom informatic pipelines for phylogenomic inference.

Software, a detailed manual, and example data files are available through github carolinafishes.github.io

The sterlet sturgeon genome sequence and the mechanisms of segmental rediploidization

Sturgeons seem to be frozen in time. The archaic characteristics of this ancient fish lineage place it in a key phylogenetic position at the base of the ~30,000 modern teleost fish species. Moreover, sturgeons are notoriously polyploid, providing unique opportunities to investigate the evolution of polyploid genomes. We assembled a high-quality chromosome-level reference genome for the sterlet, Acipenser ruthenus. Our analysis revealed a very low protein evolution rate that is at least as slow as in other deep branches of the vertebrate tree, such as that of the coelacanth. We uncovered a whole-genome duplication that occurred in the Jurassic, early in the evolution of the entire sturgeon lineage. Following this polyploidization, the rediploidization of the genome included the loss of whole chromosomes in a segmental deduplication process. While known adaptive processes helped conserve a high degree of structural and functional tetraploidy over more than 180 million years, the reduction of redundancy of the polyploid genome seems to have been remarkably random.

A genome assembly of the sterlet, Acipenser ruthenus, reveals a whole-genome duplication early in the evolution of the entire sturgeon lineage and provides details about the rediploidization of the genome.

The Polycomb Orthologues in Teleost Fishes and Their Expression in the Zebrafish Model

The Polycomb Repressive Complex 1 (PRC1) is a chromatin-associated protein complex involved in transcriptional repression of hundreds of genes controlling development and differentiation processes, but also involved in cancer and stem cell biology. Within the canonical PRC1, members of Pc/CBX protein family are responsible for the targeting of the complex to specific gene loci. In mammals, the Pc/CBX protein family is composed of five members generating, through mutual exclusion, different PRC1 complexes with potentially distinct cellular functions. Here, we performed a global analysis of the cbx gene family in 68 teleost species and traced the distribution of the cbx genes through teleost evolution in six fish super-orders. We showed that after the teleost-specific whole genome duplication, cbx4, cbx7 and cbx8 are retained as pairs of ohnologues. In contrast, cbx2 and cbx6 are present as pairs of ohnologues in the genome of several teleost clades but as singletons in others. Furthermore, since zebrafish is a widely used vertebrate model for studying development, we report on the expression of the cbx family members during zebrafish development and in adult tissues. We showed that all cbx genes are ubiquitously expressed with some variations during early development.

A new species of deropristid trematode from the sterlet Acipenser ruthenus (Actinopterygii: Acipenseridae) and revision of superfamily affiliation of the family Deropristidae

A new species, Skrjabinopsolus nudidorsalis sp. nov. is described from the sterlet Acipenser ruthenus, caught in the River Volga basin (Russia). This species differs from previously described congeners by the absence of vitelline follicles on the dorsal side of the body. The complete 18S rRNA and partial 28S rRNA gene sequences obtained for S. nudidorsalis are the first molecular data for the family Deropristidae. The results of phylogenetic analysis indicate that Deropristidae is sister to the Monorchiidae + Lissorchiidae group. The results of the phylogenetic study contradict the current taxonomic hypothesis that Deropristidae belongs to the superfamily Lepocreadioidea and allow inclusion of this family in Monorchioidea. The morphological similarity of deropristids to other monorchioids is recognizable from the presence of a bipartite internal seminal vesicle, spinous cirrus and a voluminous, armed metraterm.

Neuroanatomical Distribution of the Serotonergic System in the Brain and Retina of Holostean Fishes, The Sister Group to Teleosts

Among actinopterygian fishes, holosteans are the phylogenetically closest group to teleosts but they have been much less studied, particularly regarding the neurochemical features of their central nervous system. The serotonergic system is one of the most important and conserved systems of neurotransmission in all vertebrates. By means of immunohistochemistry against serotonin (5-hydroxytryptamine), we have conducted a comprehensive and complete description of this system in the brain and retina of representative species of the 3 genera of holostean fishes, belonging to the only 2 extant orders, Amiiformes and Lepisosteiformes. Serotonin-immunoreactive cell groups were detected in the preoptic area, the hypothalamic paraventricular organ, the epiphysis, the pretectal region, the long and continuous column of the raphe, the spinal cord, and the inner nuclear layer of the retina. Specifically, the serotonergic cell groups in the preoptic area, the epiphysis, the pretectum, and the retina had never been identified in previous studies in this group of fishes. Widespread serotonergic innervation was observed in all main brain regions, but more abundantly in the subpallium, the hypothalamus, the habenula, the optic tectum, the so-called cerebellar nucleus, and the area postrema. The comparative analysis of these results with those in other groups of vertebrates reveals some extremely conserved features, such as the presence of serotonergic cells in the retina, the pineal organ, and the raphe column, while other characteristics, like the serotonergic populations in the preoptic area, the paraventricular organ, the pretectum, and the spinal cord are generally present in all fish groups, but have been lost in most amniotes.

Developmental and Cellular Basis of Vertical Bar Color Patterns in the East African Cichlid Fish Haplochromis latifasciatus

The East African adaptive radiations of cichlid fishes are renowned for their diversity in coloration. Yet, the developmental basis of pigment pattern formation remains largely unknown. One of the most common melanic patterns in cichlid fishes are vertical bar patterns. Here we describe the ontogeny of this conspicuous pattern in the Lake Kyoga species Haplochromis latifasciatus. Beginning with the larval stages we tracked the formation of this stereotypic color pattern and discovered that its macroscopic appearance is largely explained by an increase in melanophore density and accumulation of melanin during the first 3 weeks post-fertilization. The embryonal analysis is complemented with cytological quantifications of pigment cells in adult scales and the dermis beneath the scales. In adults, melanic bars are characterized by a two to threefold higher density of melanophores than in the intervening yellow interbars. We found no strong support for differences in other pigment cell types such as xanthophores. Quantitative PCRs for twelve known pigmentation genes showed that expression of melanin synthesis genes tyr and tyrp1a is increased five to sixfold in melanic bars, while xanthophore and iridophore marker genes are not differentially expressed. In summary, we provide novel insights on how vertical bars, one of the most widespread vertebrate color patterns, are formed through dynamic control of melanophore density, melanin synthesis and melanosome dispersal.

Global determinants of freshwater and marine fish genetic diversity

Genetic diversity is estimated to be declining faster than species diversity under escalating threats, but its spatial distribution remains poorly documented at the global scale. Theory predicts that similar processes should foster congruent spatial patterns of genetic and species diversity, but empirical studies are scarce. Using a mined database of 50,588 georeferenced mitochondrial DNA barcode sequences (COI) for 3,815 marine and 1,611 freshwater fish species respectively, we examined the correlation between genetic diversity and species diversity and their global distributions in relation to climate and geography. Genetic diversity showed a clear spatial organisation, but a weak association with species diversity for both marine and freshwater species. We found a predominantly positive relationship between genetic diversity and sea surface temperature for marine species. Genetic diversity of freshwater species varied primarily across the regional basins and was negatively correlated with average river slope. The detection of genetic diversity patterns suggests that conservation measures should consider mismatching spatial signals across multiple facets of biodiversity.

Biogeographic patterns of genetic diversity are poorly documented, especially for fish species. Here the authors show that (mitochondrial) genetic diversity has global spatial organization patterns with different environmental drivers for marine and freshwater fishes, where genetic diversity is only partly congruent with species richness.

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.

Evolution of the androgen receptor: Perspectives from human health to dancing birds

Androgenic hormones orchestrate the development and activation of diverse reproductive phenotypes across vertebrates. Although extensive work investigates how selection for these traits modifies individual elements of this signaling system (e.g., hormone or androgen receptor [AR] levels), we know less about natural variation in the AR sequence across vertebrates. Our knowledge of AR sequence mutations is largely limited to work in human patients or cell-lines, providing a framework to contextualize single mutations at the expense of evolutionary timescale. Here we unite both perspectives in a review that explores the functional significance of AR on a domain-by-domain basis, using existing knowledge to highlight how and why each region might evolve. We then examine AR sequence variation on different timescales by examining sequence variation in clades originating in the Cambrian (vertebrates; >500 mya) and Cretaceous (birds; >65 mya). In each case, we characterize how the receptor has changed over time and discuss which regions are most likely to evolve in response to selection. Overall, domains that are required for androgenic signaling to function (e.g., DNA- and ligand-binding) tend to be conserved. Meanwhile, areas that interface with co-regulatory molecules can exhibit notable variation even between closely related species. We propose that accumulating mutations in regulatory regions is one way that AR structure might act as a substrate for selection to guide the evolution of reproductive traits. By synthesizing literature across disciplines and highlighting the evolutionary potential of specific AR regions, we hope to inspire new avenues of integrative research into endocrine system evolution.

Mechanisms and Drivers for the Establishment of Life Cycle Complexity in Myxozoan Parasites

It is assumed that complex life cycles in cnidarian parasites belonging to the Myxozoa result from incorporation of vertebrates into simple life cycles exploiting aquatic invertebrates. However, nothing is known about the driving forces and implementation of this event, though it fostered massive diversification. We performed a comprehensive search for myxozoans in evolutionary ancient fishes (Chondrichthyes), and more than doubled existing 18S rDNA sequence data, discovering seven independent phylogenetic lineages. We performed cophylogenetic and character mapping methods in the largest monophyletic dataset and demonstrate that host and parasite phylogenies are strongly correlated, and that tectonic changes may explain phylogeographic clustering in recent skates and softnose skates, in the Atlantic. The most basal lineages of myxozoans inhabit the bile of chondrichthyans, an immunologically privileged site and protective niche, easily accessible from the gut via the bile duct. We hypothesize that feed-integration is a likely mechanism of host acquisition, an idea supported by feeding habits of chimaeras and ancient sharks and by multiple entries of different parasite lineages from invertebrates into the new host group. We provide exciting first insights into the early evolutionary history of ancient metazoan parasites in a host group that embodies more evolutionary distinctiveness than most other vertebrates.

Characterization of two kcnk3 genes in Nile tilapia (Oreochromis niloticus): Molecular cloning, tissue distribution, and transcriptional changes in various salinity of seawater

As one important member of the two-pore-domain potassium channel (K_2P) family, potassium channel subfamily K member 3 (KCNK3) has been reported for thermogenesis regulation, energy homeostasis, membrane potential conduction, and pulmonary hypertension in mammals. However, its roles in fishes are far less examined and published. In the present study, we identified two kcnk3 genes (kcnk3a and kcnk3b) in an euryhaline fish, Nile tilapia (Oreochromis niloticus), by molecular cloning, genomic survey and laboratory experiments to investigate their potential roles for osmoregulation. We obtained full-length coding sequences of the kcnk3a and kcnk3b genes (1209 and 1173 bp), which encode 402 and 390 amino acids, respectively. Subsequent multiple sequence alignments, putative 3D-structure model prediction, genomic survey and phylogenetic analysis confirmed that two kcnk3 paralogs are widely presented in fish genomes. Interestingly, a DNA fragment inversion of a kcnk3a cluster was found in Cypriniforme in comparison with other fishes. Quantitative real-time PCRs demonstrated that both the tilapia kcnk3 genes were detected in all the examined tissues with a similar distribution pattern, and the highest transcriptions were observed in the heart. Meanwhile, both kcnk3 genes in the gill were proved to have a similar transcriptional change pattern in response to various salinity of seawater, implying that they might be involved in osmoregulation. Furthermore, three predicted transcription factors (arid3a, arid3b, and arid5a) of both kcnk3 genes also showed a similar pattern as their target genes in response to the various salinity, suggesting their potential positive regulatory roles. In summary, we for the first time characterized the two kcnk3 genes in Nile tilapia, and demonstrated their potential involvement in osmoregulation for this economically important fish.

The vertebrate Aqp14 water channel is a neuropeptide-regulated polytransporter

Water channels (aquaporins) were originally discovered in mammals with fourteen subfamilies now identified (AQP0-13). Here we show that a functional Aqp14 subfamily phylogenetically related to AQP4-type channels exists in all vertebrate lineages except hagfishes and eutherian mammals. In contrast to the water-selective classical aquaporins, which have four aromatic-arginine constriction residues, Aqp14 proteins present five non-aromatic constriction residues and facilitate the permeation of water, urea, ammonia, H2O2 and glycerol. Immunocytochemical assays suggest that Aqp14 channels play important osmoregulatory roles in piscine seawater adaptation. Our data indicate that Aqp14 intracellular trafficking is tightly regulated by the vasotocinergic/isotocinergic neuropeptide and receptor systems, whereby protein kinase C and A transduction pathways phosphorylate highly conserved C-terminal residues to control channel plasma membrane insertion. The neuropeptide regulation of Aqp14 channels thus predates the vasotocin/vasopressin regulation of AQP2-5-6 orthologs observed in tetrapods. These findings demonstrate that vertebrate Aqp14 channels represent an ancient subfamily of neuropeptide-regulated polytransporters.

Comparative phylogeography of trans-Andean freshwater fishes based on genome-wide nuclear and mitochondrial markers

The Neotropical region represents one of the greatest biodiversity hot spots on earth. Despite its unparalleled biodiversity, regional comparative phylogeographic studies are still scarce, with most focusing on model clades (e.g. birds) and typically examining a handful of loci. Here, we apply a genome-wide comparative phylogeographic approach to test hypotheses of codiversification of freshwater fishes in the trans-Andean region. Using target capture methods, we examined exon data for over 1,000 loci combined with complete mitochondrial genomes to study the phylogeographic history of five primary fish species (>150 individuals) collected from eight major river basins in Northwestern South America and Lower Central America. To assess their patterns of genetic structure, we inferred genealogical concordance taking into account all major aspects of phylogeography (within loci, across multiple genes, across species and among biogeographic provinces). Based on phylogeographic concordance factors, we tested four a priori biogeographic hypotheses, finding support for three of them and uncovering a novel, unexpected pattern of codiversification. The four emerging inter-riverine patterns are as follows: (a) Tuira + Atrato, (b) Ranchería + Catatumbo, (c) Magdalena system and (d) Sinú + Atrato. These patterns are interpreted as shared responses to the complex uplifting and orogenic processes that modified or sundered watersheds, allowing codiversification and speciation over geological time. We also find evidence of cryptic speciation in one of the species examined and instances of mitochondrial introgression in others. These results help further our knowledge of the historical geographic factors shaping the outstanding biodiversity of the Neotropics.

Comparative Analysis of the Organization of the Catecholaminergic Systems in the Brain of Holostean Fishes (Actinopterygii/Neopterygii)

Living holosteans, comprising 8 species of bowfins and gars, form a small monophyletic group of actinopterygian fishes, which are currently considered as the sister group to the enormously numerous teleosts and have largely been neglected in neuroanatomical studies. We have studied the catecholaminergic (CAergic) systems by means of antibodies against tyrosine hydroxylase (TH) and dopamine (DA) in the brain of representative species of the 3 genera included in the 2 orders of holostean fishes: Amia calva (Amiiformes) and Lepisosteus platyrhincus, Lepisosteus oculatus, and Atractosteus spatula (Lepisosteiformes). Different groups of TH/DA-immunoreactive (ir) cells were observed in the olfactory bulb, subpallium, and preoptic area of the telencephalon. Hypothalamic groups were labeled in the suprachiasmatic nucleus, tuberal (only in A. calva), retrotuberal, and retromamillary areas; specifically, the paraventricular organ showed only DA immunoreactivity. In the diencephalon, TH/DA-ir groups were detected in the prethalamus, posterior tubercle, and pretectum. In the caudal hindbrain, the solitary tract nucleus and area postrema presented TH/DA-ir cell groups, and also the spinal cord and the retina. Only in A. calva, particular CAergic cell groups were observed in the habenula, the mesencephalic tegmentum, and in the locus coeruleus. Following a neuromeric analysis, the comparison of these results with those obtained in other classes of fishes and tetrapods shows many common traits of CAergic systems shared by most vertebrates and in addition highlights unique features of actinopterygian fishes.

Serotonin systems in three socially communicating teleost species, the grunting toadfish (Allenbatrachus grunniens), a South American marine catfish (Ariopsis seemanni), and the upside-down catfish (Synodontis nigriventris)

We investigated immunohistochemically the distribution of serotonergic cell populations in three teleost species (one toadfish, Allenbatrachus grunniens, and two catfishes, Synodontis nigriventris and Ariopsis seemanni). All three species exhibited large populations of 5-HT positive neurons in the paraventricular organ (PVO) and the dorsal (Hd) and caudal (Hc) periventricular hypothalamic zones, plus a smaller one in the periventricular pretectum, a few cells in the pineal stalk, and - only in catfishes - in the preoptic region. Furthermore, the rhombencephalic superior and inferior raphe always contained ample serotonergic cells. In each species, a neuronal mass extended into the hypothalamic lateral recess. Only in the toadfish, did this intraventricular structure contain serotonergic cells and arise from Hd, whereas in the catfishes it emerged from medially and represents the dorsal tuberal nucleus seen in other catfishes as well. Serotonergic cells in PVO, Hd and Hc were liquor-contacting. Those of the PVO extended into the midline area of the periventricular posterior tubercular nucleus in both catfishes. Dopaminergic, liquor-contacting neurons were additionally investigated using an antibody against tyrosine hydroxylase (TH) in S. nigriventris showing that TH was never co-localized with serotonin. Because TH antibodies are known to reveal mostly or only the TH1 enzyme, we hypothesize that th1-expressing dopamine cells (unlike th2-expressing ones) do not co-localize with serotonin. Since the three investigated species engage in social communication using swim bladder associated musculature, we investigated the serotonergic innervation of the hindbrain vocal or electromotor nuclei initiating the social signal. We found in all three species serotonergic fibers seemingly originating from close-by serotonergic neurons of inferior raphe or anterior spinal cord. Minor differences appear to be rather species-specific than dependent on the type of social communication.

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.

Evolution and diversity of transposable elements in fish genomes

Transposable elements (TEs) are genomic sequences that can move, multiply, and often form sizable fractions of vertebrate genomes. Fish belong to a unique group of vertebrates, since their karyotypes and genome sizes are more diverse and complex, with probably higher diversity and evolution specificity of TE. To investigate the characteristics of fish TEs, we compared the mobilomes of 39 species, and observed significant variation of TE content in fish (from 5% in pufferfish to 56% in zebrafish), along with a positive correlation between fish genome size and TE content. In different classification hierarchies, retrotransposons (class), long terminal repeat (order), as well as Helitron, Maverick, Kolobok, CMC, DIRS, P, I, L1, L2, and 5S (superfamily) were all positively correlated with fish genome size. Consistent with previous studies, our data suggested fish genomes to not always be dominated by DNA transposons; long interspersed nuclear elements are also prominent in many species. This study suggests CR1 distribution in fish genomes to be obviously regular, and provides new clues concerning important events in vertebrate evolution. Altogether, our results highlight the importance of TEs in the structure and evolution of fish genomes and suggest fish species diversity to parallel transposon content diversification.

The cephalic lateral-line system of Characiformes (Teleostei: Ostariophysi): anatomy and phylogenetic implications

The lateral-line system has been traditionally recognized as an important source of phylogenetic information for different groups of fishes. Although extensively studied in Siluriformes and Cypriniformes, the lateral-line system of Characiformes remained underexplored. In the present study, the anatomy of the cephalic lateral-line canals of characiforms is described in detail and a unifying terminology that considers the ontogeny and homologies of the components of this system is offered. Aspects of the arrangement of lateral-line canals, as well as the number, location and size of canal tubules and pores, resulted in the identification of novel putative synapomorphies for Characiformes and several of its subgroups. The study also revised synapomorphies previously proposed for different characiform families and provided comments on their observed distribution across the order based on extensive taxon sampling. Information from the ontogenetic studies of the cephalic lateral-line canal system and a proposal for the proper use of these data to detect truncations in the development of the lateral-line canals across the order is also offered.

From Africa to Antarctica: Exploring the Metabolism of Fish Heart Mitochondria Across a Wide Thermal Range

The thermal sensitivity of ectotherms is largely dictated by the impact of temperature on cellular bioenergetics, particularly on mitochondrial functions. As the thermal sensitivity of bioenergetic pathways depends on the structural and kinetic properties of its component enzymes, optimization of their collective function to different thermal niches is expected to have occurred through selection. In the present study, we sought to characterize mitochondrial phenotypic adjustments to thermal niches in eight ray-finned fish species occupying a wide range of thermal habitats by comparing the activities of key mitochondrial enzymes in their hearts. We measured the activity of four enzymes that control substrate entrance into the tricarboxylic acid (TCA) cycle: pyruvate kinase (PK), pyruvate dehydrogenase complex (PDHc), carnitine palmitoyltransferase (CPT), and hydroxyacyl-CoA dehydrogenase (HOAD). We also assayed enzymes of the electron transport system (ETS): complexes I, II, I + III, and IV. Enzymes were assayed at five temperatures (5, 10, 15, 20, and 25°C). Our results showed that the activity of CPT, a gatekeeper of the fatty acid pathway, was higher in the cold-water fish than in the warmer-adapted fish relative to the ETS (complexes I and III) when measured close to the species optimal temperatures. The activity of HOAD showed a similar pattern relative to CI + III and thermal environment. By contrast, PDHc and PK did not show the similar patterns with respect to CI + III and temperature. Cold-adapted species had high CIV activities compared to those of upstream complexes (I, II, I + III) whereas the converse was true for warm-adapted species. Our findings reveal a significant variability of heart mitochondrial organization among species that can be linked to temperature adaptation. Cold-adapted fish do not appear to compensate for PDHc activity but likely adjust fatty acids oxidation through higher activities of CPT and HOAD relative to complexes I + III.

Subfunctionalization and evolution of liver-expressed antimicrobial peptide 2 (LEAP2) isoform genes in Siberian sturgeon (Acipenser baerii), a primitive chondrostean fish species

Two liver-expressed antimicrobial peptide 2 (LEAP2) isoforms were characterized in a primitive chondrostean sturgeon species, Acipenser baerii (Acipenseriformes). A. baerii LEAP2 isoforms represented essentially common structures shared by their vertebrate orthologs at both genomic (i.e., tripartite organization) and peptide (two conserved disulfide bonds) levels. A. baerii LEAP2 isoforms (designed LEAP2AB and LEAP2C, respectively) phylogenetically occupy the most basal position in the actinopterygian lineage and represent an intermediate character between teleostean and tetrapodian LEAP2s in the sequence alignment. Molecular phylogenetic analysis including LEAP2s from extant primitive fish species indicated that the evolutionary origin of ancestral LEAP2 in vertebrate groups should date back to earlier than the actinopterygian-sarcopterygian split. Gene expression assays under both basal and stimulated conditions suggested that A. baerii LEAP2 isoforms have undergone substantial subfunctionalization in tissue distribution pattern, developmental/ontogenetic expression, and immune responses. LEAP2AB showed a predominant liver expression, while LEAP2C exhibited the highest level of expression in the intestine. LEAP2C was a more dominantly expressed isoform during embryonic development and prelarval ontogeny. The LEAP2AB isoform is more closely associated with innate immune response to microbial invasion, compared with LEAP2C, as evidenced by results from LPS, poly(I:C) and Aeromonas hydrophila challenges. Synthetic mature peptides of LEAP2AB displayed a more potent antimicrobial activity than did LEAP2C. Data from this study could be useful not only to provide deeper insights into the evolutionary mechanism of LEAP2 in the actinopterygian lineage but also to better understand the innate immunity of this commercially important chondrostean species.

A phylogenomic framework for pelagiarian fishes (Acanthomorpha: Percomorpha) highlights mosaic radiation in the open ocean

The fish clade Pelagiaria, which includes tunas as its most famous members, evolved remarkable morphological and ecological variety in a setting not generally considered conducive to diversification: the open ocean. Relationships within Pelagiaria have proven elusive due to short internodes subtending major lineages suggestive of rapid early divergences. Using a novel sequence dataset of over 1000 ultraconserved DNA elements (UCEs) for 94 of the 286 species of Pelagiaria (more than 70% of genera), we provide a time-calibrated phylogeny for this widely distributed clade. Some inferred relationships have clear precedents (e.g. the monophyly of 'core' Stromateoidei, and a clade comprising 'Gempylidae' and Trichiuridae), but others are unexpected despite strong support (e.g. Chiasmodontidae + Tetragonurus). Relaxed molecular clock analysis using node-based fossil calibrations estimates a latest Cretaceous origin for Pelagiaria, with crown-group families restricted to the Cenozoic. Estimated mean speciation rates decline from the origin of the group in the latest Cretaceous, although credible intervals for root and tip rates are broad and overlap in most cases, and there is higher-than-expected partitioning of body shape diversity (measured as fineness ratio) between clades concentrated during the Palaeocene-Eocene. By contrast, more direct measures of ecology show either no substantial deviation from a null model of diversification (diet) or patterns consistent with evolutionary constraint or high rates of recent change (depth habitat). Collectively, these results indicate a mosaic model of diversification. Pelagiarians show high morphological disparity and modest species richness compared to better-studied fish radiations in contrasting environments. However, this pattern is also apparent in other clades in open-ocean or deep-sea habitats, and suggests that comparative study of such groups might provide a more inclusive model of the evolution of diversity in fishes.

Assexon: Assembling Exon Using Gene Capture Data

Exon capture across species has been one of the most broadly applied approaches to acquire multi-locus data in phylogenomic studies of non-model organisms. Methods for assembling loci from short-read sequences (eg, Illumina platforms) that rely on mapping reads to a reference genome may not be suitable for studies comprising species across a wide phylogenetic spectrum; thus, de novo assembling methods are more generally applied. Current approaches for assembling targeted exons from short reads are not particularly optimized as they cannot (1) assemble loci with low read depth, (2) handle large files efficiently, and (3) reliably address issues with paralogs. Thus, we present Assexon: a streamlined pipeline that de novo assembles targeted exons and their flanking sequences from raw reads. We tested our method using reads from Lepisosteus osseus (4.37 Gb) and Boleophthalmus pectinirostris (2.43 Gb), which are captured using baits that were designed based on genome sequence of Lepisosteus oculatus and Oreochromis niloticus, respectively. We compared performance of Assexon to PHYLUCE and HybPiper, which are commonly used pipelines to assemble ultra-conserved element (UCE) and Hyb-seq data. A custom exon capture analysis pipeline (CP) developed by Yuan et al was compared as well. Assexon accurately assembled more than 3400 to 3800 (20%-28%) loci than PHYLUCE and more than 1900 to 2300 (8%-14%) loci than HybPiper across different levels of phylogenetic divergence. Assexon ran at least twice as fast as PHYLUCE and HybPiper. Number of loci assembled using CP was comparable with Assexon in both tests, while Assexon ran at least 7 times faster than CP. In addition, some steps of CP require the user’s interaction and are not fully automated, and this user time was not counted in our calculation. Both Assexon and CP retrieved no paralogs in the testing runs, but PHYLUCE and Hybpiper did. In conclusion, Assexon is a tool for accurate and efficient assembling of large read sets from exon capture experiments. Furthermore, Assexon includes scripts to filter poorly aligned coding regions and flanking regions, calculate summary statistics of loci, and select loci with reliable phylogenetic signal. Assexon is available at https://github.com/yhadevol/Assexon.

Draft Genome and Complete Hox-Cluster Characterization of the Sterlet (Acipenser ruthenus)

Background: Sturgeons (Chondrostei: Acipenseridae) are a group of “living fossil” fishes at a basal position among Actinopteri. They have raised great public interest due to their special evolutionary position, species conservation challenges, as well as their highly-prized eggs (caviar). The sterlet, Acipenser ruthenus, is a relatively small-sized member of sturgeons and has been widely distributing in both Europe and Asia. In this study, we performed whole genome sequencing, de novo assembly and gene annotation of the tarlet to construct its draft genome.

Findings: We finally obtained a 1.83-Gb genome assembly (BUSCO completeness of 81.6%) from a total of 316.8-Gb raw reads generated by an Illumina Hiseq 2500 platform. The scaffold N50 and contig N50 values reached 191.06 and 18.88 kb, respectively. The sterlet genome was predicted to be comprised of 42.84% repeated sequences and to contain 22,184 protein-coding genes, of which 21,112 (95.17%) have been functionally annotated with at least one hit in public databases. A genetic phylogeny demonstrated that the sterlet is situated in the basal position among ray-finned fishes and 4dTv analysis estimated that a recent whole genome duplication occurred 21.3 million years ago. Moreover, seven Hox clusters carrying 68 Hox genes were characterized in the sterlet. Phylogeny of HoxA clusters in the sterlet and American paddlefish divided these sturgeons into two groups, confirming the independence of each lineage-specific genome duplication in Acipenseridae and Polyodontidae.

Conclusions: This draft genome makes up for the lack of genomic and molecular data of the sterlet and its Hox clusters. It also provides a genetic basis for further investigation of lineage-specific genome duplication and the early evolution of ray-finned fishes.

Transcriptome analyses reveal alterations in muscle metabolism, immune responses and reproductive behavior of Japanese mantis shrimp (Oratosquilla oratoria) at different cold temperature

Low temperature reduction is thought to cause widespread effects on the physical and behavioral traits of marine organisms, which include metabolic processes, immune responses, and reproductive behavior. Crustaceans are generally considered sensitive to temperature reduction due to the lack of efficient regulators. To better understand the molecular regulatory mechanisms of crustacean exposure to cold stress, Japanese mantis shrimp (Oratosquilla oratoria) was chosen as a representative crustacean. Transcriptomic responses in O. oratoria from five temperatures (25 °C, 22 °C, 19 °C, 16 °C, and 13 °C) were studied using RNA-seq. A total of 64.91 Gb of clean transcriptomic data were generated in 10 libraries and then spliced into 52,107 unigenes with an average length of 1089 bp and an N50 length of 1872 bp. A total of 14,841 unigenes was annotated in at least one database using Blastx alignment. Compared with the control temperature (25 °C), 7, 21, 58, and 236 unigenes were significantly differentially expressed at 22 °C, 19 °C, 16 °C, and 13 °C, respectively. GO analysis showed that 6, 20, 27, and 35 terms were significantly enriched at 22 °C, 19 °C, 16 °C, and 13 °C, respectively. In addition, 2, 5, 2, and 10 significant pathways were presented at 22 °C, 19 °C, 16 °C, and 13 °C, respectively. Combining NR, GO, and KEGG annotation information, many genes significantly differentially expressed at low temperatures may be associated with metabolic processes, immune response, and reproductive behavior. Additionally, we reconstructed the phylogenetic relationship based on 366 orthologous genes and the predicted differentiation time of O. oratoria and P. vannamei range from 212.82 to 365.30 Mya. Furthermore, 16 orthologous genes were identified as PSGs and 30 orthologous genes were identified as FEGs and these adaptive genes were associated with energy metabolism, stress response and immunity, and multiple cellular processing. These results provide fundamental information about molecular mechanisms regulating cold stress response of O. oratoria.

Whole-Genome Sequencing of the Giant Devil Catfish, Bagarius yarrelli

As one economically important fish in the southeastern Himalayas, the giant devil catfish (Bagarius yarrelli) has been known for its extraordinarily large body size. It can grow up to 2 m, whereas the non-Bagarius sisorids only reach 10-30 cm. Another outstanding characteristic of Bagarius species is the salmonids-like reddish flesh color. Both body size and flesh color are interesting questions in science and also valuable features in aquaculture that worth of deep investigations. Bagarius species therefore are ideal materials for studying body size evolution and color depositions in fish muscles, and also potential organisms for extensive utilization in Asian freshwater aquaculture. In a combination of Illumina and PacBio sequencing technologies, we de novo assembled a 571-Mb genome for the giant devil catfish from a total of 153.4-Gb clean reads. The scaffold and contig N50 values are 3.1 and 1.6 Mb, respectively. This genome assembly was evaluated with 93.4% of Benchmarking Universal Single-Copy Orthologs completeness, 98% of transcripts coverage, and highly homologous with a chromosome-level-based genome of channel catfish (Ictalurus punctatus). We detected that 35.26% of the genome assembly is composed of repetitive elements. Employing homology, de novo, and transcriptome-based annotations, we annotated a total of 19,027 protein-coding genes for further use. In summary, we generated the first high-quality genome assembly of the giant devil catfish, which provides an important genomic resource for its future studies such as the body size and flesh color issues, and also for facilitating the conservation and utilization of this valuable catfish.

Deep evolutionary origin of limb and fin regeneration

Salamanders and lungfishes are the only sarcopterygians (lobe-finned vertebrates) capable of paired appendage regeneration, regardless of the amputation level. Among actinopterygians (ray-finned fishes), regeneration after amputation at the fin endoskeleton has only been demonstrated in polypterid fishes (Cladistia). Whether this ability evolved independently in sarcopterygians and actinopterygians or has a common origin remains unknown. Here we combine fin regeneration assays and comparative RNA-sequencing (RNA-seq) analysis of Polypterus and axolotl blastemas to provide support for a common origin of paired appendage regeneration in Osteichthyes (bony vertebrates). We show that, in addition to polypterids, regeneration after fin endoskeleton amputation occurs in extant representatives of 2 other nonteleost actinopterygians: the American paddlefish (Chondrostei) and the spotted gar (Holostei). Furthermore, we assessed regeneration in 4 teleost species and show that, with the exception of the blue gourami (Anabantidae), 3 species were capable of regenerating fins after endoskeleton amputation: the white convict and the oscar (Cichlidae), and the goldfish (Cyprinidae). Our comparative RNA-seq analysis of regenerating blastemas of axolotl and Polypterus reveals the activation of common genetic pathways and expression profiles, consistent with a shared genetic program of appendage regeneration. Comparison of RNA-seq data from early Polypterus blastema to single-cell RNA-seq data from axolotl limb bud and limb regeneration stages shows that Polypterus and axolotl share a regeneration-specific genetic program. Collectively, our findings support a deep evolutionary origin of paired appendage regeneration in Osteichthyes and provide an evolutionary framework for studies on the genetic basis of appendage regeneration.

The untapped potential of medaka and its wild relatives

The medaka is a fish that has served as a model organism for over a century, yet there is still much to learn about its life in the wild. Here we summarize the current knowledge, highlight recent progress and outline remaining gaps in our understanding of the natural history of medaka. It has also become clear over time that rather than being a single species, medaka comprises an entire species complex, so disentangling the species boundaries is an important goal for future research. Moreover, medaka and other ricefishes exhibit striking functional diversity, little of which has been investigated to date. As such, there are opportunities to use the resources developed for medaka to study other ricefishes, and to learn more about medaka itself in an evolutionary context.

Pharyngeal Jaws Converge by Similar Means, Not to Similar Ends, When Minnows (Cypriniformes: Leuciscidae) Adapt to New Dietary Niches

Convergent evolution is at the forefront of many form-function studies. There are many examples of multiple independent lineages evolving a similar morphology in response to similar functional demands, providing a framework for testing hypotheses of form-function evolution. However, there are numerous clades with underappreciated convergence, in which there is a perceived homogeneity in morphology. In these groups, it can be difficult to investigate causal relationships of form and function (e.g., diet influencing the evolution of jaw morphology) without the ability to disentangle phylogenetic signal from convergence. Leuciscids (Cypriniformes: Leuciscidae; formerly nested within Cyprinidae) are a species-rich clade of fishes that have diversified to occupy nearly every freshwater trophic niche, yet are considered to have relatively low morphological diversity relative to other large freshwater clades. Within the North American leuciscids, many genera contain at least one herbivore, insectivore, and larvaphage. We created 3D models from micro-computed tomography scans of 165 leuciscid species to measure functionally relevant traits within the pharyngeal jaws of these fishes. Using a published phylogeny, we tested these metrics for evolutionary integration, phylogenetic signal, and correlation with diet. Measurements of the pharyngeal jaws, muscle attachment areas, and teeth showed strong positive evolutionary correlation with each other and negative evolutionary correlation with measurements of the inter-ceratobranchial ligament (ICB ligament). Using diet data from published literature, we found extensive dietary convergence within Leuciscidae. The most common transitions we found were between herbivorous and invertivorous taxa and between insectivore types (aquatic vs. terrestrial). We document a trade-off in which herbivorous leuciscids have large teeth, short ICB ligaments, and large muscle attachment areas, whereas insectivorous leuciscids showed the opposite pattern. Inverse patterns of morphological integration between the ICB ligament the rest of the pharyngeal jaw correspond this dietary trade-off, which indicates that coordinated evolution of morphological traits contributes to functional diversity in this clade. However, these patterns only emerge in the context of phylogeny, meaning that the pharyngeal jaws of North American leuciscids converge by similar means (structural changes in response to dietary demands), but not necessarily to similar ends (absolute phenotype).

Chromosomal and Genomic Dynamics of Satellite DNAs in Characidae (Characiformes, Teleostei) Species

Satellite DNAs (satDNAs) are tandemly repeated DNA sequences with great abundance in eukaryotic genomes. A single species may carry up to hundreds of satDNA families, which is collectively called as "satellitome," each showing its own dynamics and evolution rates. In this context, all live species contain a satDNA library that may be partially or totally shared with other related species/populations. In the late few years, next-generation sequencing (NGS) and novel bioinformatic tools facilitated the massive characterization of these sequences at low costs, and consequently, comparing satDNAs between species. In this study, we characterized two novel satDNAs (MsaSat03-80 and MsaSat04-142) in three characid fish (Astyanax paranae and Astyanax fasciatus and two populations of Moenkhausia sanctaefilomenae) and mapped their chromosomal location to unveil the evolutionary dynamics of satDNA repeats in those species. Our results evidenced that MsaSat03 is present in the genomes of all analyzed species, but is clustered only in the chromosomes of M. sanctaefilomenae, exhibiting a conserved number and location of sites. Conversely, MsaSat04 sequences is restricted to M. sanctaefilomenae and shows a differential distribution between the two analyzed populations. Altogether, our analyses point to a complex history of satDNA families in characid fish and the utility of NGS data for comparative satDNA analysis.

Staging and normal table of postembryonic development of the clownfish (Amphiprion ocellaris)

Background

The clownfish Amphiprion ocellaris is one of the rare coral reef fish species that can be reared in aquaria. With relatively short embryonic and larval development, it could be used as a model species to study the impact of global changes such as temperature rise or anthropogenic threats (eg, pollution) on the postembryonic development at molecular and endocrinological levels. Establishing a developmental table allows us to standardize sampling for the scientific community willing to conduct experiments on this species on different areas: ecology, evolution, and developmental biology.

Results

Here, we describe the postembryonic developmental stages for the clownfish A. ocellaris from hatching to juvenile stages (30 days posthatching). We quantitatively followed the postembryonic growth and described qualitative traits: head, paired and unpaired fins, notochord flexion, and pigmentation changes. The occurrence of these changes over time allowed us to define seven stages, for which we provide precise descriptions.

Conclusions

Our work gives an easy system to determine A. ocellaris postembryonic stages allowing, thus, to develop this species as a model species for coral reef fishes. In light of global warming, the access to the full postembryonic development stages of coral reef fish is important to determine stressors that can affect such processes.

Seven developmental stages have been identified to describe the larval development of the clownfish Amphiprion ocellaris.

Clownfish larvae undergo two distinct developmental growth phases that correspond to growth and metamorphosis.

A dichotomous key determination has been created to assist users in identifying the various developmental stages.

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.

Stable transgenesis in Astyanax mexicanus using the Tol2 transposase system

BACKGROUND

Astyanax mexicanus is a well-established fish model system for evolutionary and developmental biology research. These fish exist as surface forms that inhabit rivers and 30 different populations of cavefish. Despite important progress in the deployment of new technologies, deep mechanistic insights into the genetic basis of evolution, development, and behavior have been limited by a lack of transgenic lines commonly used in genetic model systems.

RESULTS

Here, we expand the toolkit of transgenesis by characterizing two novel stable transgenic lines that were generated using the highly efficient Tol2 system, commonly used to generate transgenic zebrafish. A stable transgenic line consisting of the zebrafish ubiquitin promoter expresses enhanced green fluorescent protein ubiquitously throughout development in a surface population of Astyanax. To define specific cell-types, a Cntnap2-mCherry construct labels lateral line mechanosensory neurons in zebrafish. Strikingly, both constructs appear to label the predicted cell types, suggesting many genetic tools and defined promoter regions in zebrafish are directly transferrable to cavefish.

CONCLUSION

The lines provide proof-of-principle for the application of Tol2 transgenic technology in A. mexicanus. Expansion on these initial transgenic lines will provide a platform to address broadly important problems in the quest to bridge the genotype-phenotype gap.

Into Africa via docked India: a fossil climbing perch from the Oligocene of Tibet helps solve the anabantid biogeographical puzzle

The northward drift of the Indian Plate and its collision with Eurasia have profoundly impacted the evolutionary history of the terrestrial organisms, especially the ones along the Indian Ocean rim. Climbing perches (Anabantidae) are primary freshwater fishes showing a disjunct south Asian-African distribution, but with an elusive paleobiogeographic history due to the lack of fossil evidence. Here, based on an updated time-calibrated anabantiform phylogeny integrating a number of relevant fossils, the divergence between Asian and African climbing perches is estimated to have occurred in the middle Eocene (ca. 40 Ma, Ma: million years ago), a time when India had already joined with Eurasia. The key fossil lineage is †Eoanabas, the oldest anabantid known so far, from the upper Oligocene of the Tibetan Plateau. Ancestral range reconstructions suggest a Southeast Asian origin in the early Eocene (ca. 48 Ma) and subsequent dispersals to Tibet and then India for this group. Thereby we propose their westbound dispersal to Africa via the biotic bridge between India and Africa. If so, climbing perch precursors had probably followed the paleobiogeographical route of snakehead fishes, which have a slightly older divergence between African and Asian taxa. As such, our study echoes some recent molecular analyses in rejecting the previously held "Gondwana continental drift vicariance" or late Mesozoic dispersal scenarios for the climbing perches, but provides a unique biogeographical model to highlight the role of the pre-uplift Tibet and the docked India in shaping the disjunct distribution of some air-breathing freshwater fishes around the Indian Ocean.