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

Modulation of bioelectric cues in the evolution of flying fishes

Changes to allometry, or the relative proportions of organs and tissues within organisms, is a common means for adaptive character change in evolution. However, little is understood about how relative size is specified during development and shaped during evolution. Here, through a phylogenomic analysis of genome-wide variation in 35 species of flying fishes and relatives, we identify genetic signatures in both coding and regulatory regions underlying the convergent evolution of increased paired fin size and aerial gliding behaviors. To refine our analysis, we intersected convergent phylogenomic signatures with mutants with altered fin size identified in distantly related zebrafish. Through these paired approaches, we identify a surprising role for an L-type amino acid transporter, lat4a, and the potassium channel, kcnh2a, in the regulation of fin proportion. We show that interaction between these genetic loci in zebrafish closely phenocopies the observed fin proportions of flying fishes. The congruence of experimental and phylogenomic findings point to conserved, non-canonical signaling integrating bioelectric cues and amino acid transport in the establishment of relative size in development and evolution.

Form and Function of the Vertebrate and Invertebrate Blood-Brain Barriers

The need to protect neural tissue from toxins or other substances is as old as neural tissue itself. Early recognition of this need has led to more than a century of investigation of the blood-brain barrier (BBB). Many aspects of this important neuroprotective barrier have now been well established, including its cellular architecture and barrier and transport functions. Unsurprisingly, most research has had a human orientation, using mammalian and other animal models to develop translational research findings. However, cell layers forming a barrier between vascular spaces and neural tissues are found broadly throughout the invertebrates as well as in all vertebrates. Unfortunately, previous scenarios for the evolution of the BBB typically adopt a classic, now discredited 'scala naturae' approach, which inaccurately describes a putative evolutionary progression of the mammalian BBB from simple invertebrates to mammals. In fact, BBB-like structures have evolved independently numerous times, complicating simplistic views of the evolution of the BBB as a linear process. Here, we review BBBs in their various forms in both invertebrates and vertebrates, with an emphasis on the function, evolution, and conditional relevance of popular animal models such as the fruit fly and the zebrafish to mammalian BBB research.

The Emerging Phylogenetic Perspective on the Evolution of Actinopterygian Fishes

The emergence of a new phylogeny of ray-finned fishes at the turn of the twenty-first century marked a paradigm shift in understanding the evolutionary history of half of living vertebrates. We review how the new ray-finned fish phylogeny radically departs from classical expectations based on morphology. We focus on evolutionary relationships that span the backbone of ray-finned fish phylogeny, from the earliest divergences among teleosts and nonteleosts to the resolution of major lineages of Percomorpha. Throughout, we feature advances gained by the new phylogeny toward a broader understanding of ray-finned fish evolutionary history and the implications for topics that span from the genetics of human health to reconsidering the concept of living fossils. Additionally, we discuss conceptual challenges that involve reconciling taxonomic classification with phylogenetic relationships and propose an alternate higher-level classification for Percomorpha. Our review highlights remaining areas of phylogenetic uncertainty and opportunities for comparative investigations empowered by this new phylogenetic perspective on ray-finned fishes.

Phylogeny of the damselfishes (Pomacentridae) and patterns of asymmetrical diversification in body size and feeding ecology

The damselfishes (family Pomacentridae) inhabit near-shore communities in tropical and temperature oceans as one of the major lineages in coral reef fish assemblages. Our understanding of their evolutionary ecology, morphology and function has often been advanced by increasingly detailed and accurate molecular phylogenies. Here we present the next stage of multi-locus, molecular phylogenetics for the group based on analysis of 12 nuclear and mitochondrial gene sequences from 345 of the 422 damselfishes. The resulting well-resolved phylogeny helps to address several important questions about higher-level damselfish relationships, their evolutionary history and patterns of divergence. A time-calibrated phylogenetic tree yields a root age for the family of 55.5 mya, refines the age of origin for a number of diverse genera, and shows that ecological changes during the Eocene-Oligocene transition provided opportunities for damselfish diversification. We explored the idea that body size extremes have evolved repeatedly among the Pomacentridae, and demonstrate that large and small body sizes have evolved independently at least 40 times and with asymmetric rates of transition among size classes. We tested the hypothesis that transitions among dietary ecotypes (benthic herbivory, pelagic planktivory and intermediate omnivory) are asymmetric, with higher transition rates from intermediate omnivory to either planktivory or herbivory. Using multistate hidden-state speciation and extinction models, we found that both body size and dietary ecotype are significantly associated with patterns of diversification across the damselfishes, and that the highest rates of net diversification are associated with medium body size and pelagic planktivory. We also conclude that the pattern of evolutionary diversification in feeding ecology, with frequent and asymmetrical transitions between feeding ecotypes, is largely restricted to the subfamily Pomacentrinae in the Indo-West Pacific. Trait diversification patterns for damselfishes across a fully resolved phylogeny challenge many recent general conclusions about the evolution of reef fishes.

Genome Mining-Based Discovery of Blenny Fish-Derived Peptides Targeting the Mouse κ-Opioid Receptor

Over the past years, peptides have attracted increasing interest for G protein-coupled receptor (GPCR) drug discovery and development. Peptides occupy a unique chemical space that is not easily accessible for small molecules and antibodies and provide advantages over these ligand classes such as lower toxicity and higher selectivity. The κ-opioid receptor (KOR) is a prototypic GPCR and an appealing therapeutic target for the development of safer and more effective analgesics. Recently, peptides have emerged as analgesic drug candidates with improved side effect profiles. We have previously identified plant-derived peptides, which activate KOR. Based on this precedent, here we relied on publicly available databases to discover novel KOR peptide ligands by genome mining. Using human preprodynorphin as a query, we identified blenny fish-derived peptides, referred to as blenniorphins, capable of binding to and activating KOR with nanomolar affinity and potency, respectively. Additionally, the blenniorphins altered β-arrestin-2 recruitment at the KOR. Our study demonstrates the utility of genome mining to identify peptide GPCR ligands with intriguing pharmacological properties and unveils the potential of blenny fishes as a source for novel KOR ligands.

Genome assembly and annotation of the tambaqui (Colossoma macropomum): an emblematic fish of the Amazon River Basin

Colossoma macropomum, known as "tambaqui", is the largest Characiformes fish in the Amazon River Basin and a leading species in Brazilian aquaculture and fisheries. Good quality meat and excellent adaptability to culture systems are some of its remarkable farming features. To support studies into the genetics and genomics of the tambaqui, we have produced the first high-quality genome for the species. We combined Illumina and PacBio sequencing technologies to generate a reference genome, assembled with 39× coverage of long reads and polished to a consensus quality value (QV) of 36 with 130× coverage of short reads. The genome was assembled into 1269 scaffolds (a total of 1,221,847,006 bases), with a scaffold N50 size of 40 Mb, where 93% of all assembled bases were placed in the largest 54 scaffolds corresponding to the diploid karyotype of the tambaqui. Furthermore, the NCBI Annotation Pipeline annotated genes, pseudogenes, and non-coding transcripts using the RefSeq database as evidence, guaranteeing a high-quality annotation. A Genome Data Viewer for the tambaqui was produced, which will benefit groups interested in exploring the unique genomic features of the species. The availability of a highly accurate genome assembly for tambaqui provides the foundation for the discovery of novel ecological and evolutionary insights, and is a helpful resource for aquaculture.

Analysis of Islet-1, Nkx2.1, Pax6, and Orthopedia in the forebrain of the sturgeon Acipenser ruthenus identifies conserved prosomeric characteristics

The distribution patterns of a set of conserved brain developmental regulatory transcription factors were analyzed in the forebrain of the basal actinopterygian fish Acipenser ruthenus, consistent with the prosomeric model. In the telencephalon, the pallium was characterized by ventricular expression of Pax6. In the subpallium, the combined expression of Nkx2.1/Islet-1 (Isl1) allowed to propose ventral and dorsal areas, as the septo-pallidal (Nkx2.1/Isl1+) and striatal derivatives (Isl1+), respectively, and a dorsal portion of the striatal derivatives, ventricularly rich in Pax6 and devoid of Isl1 expression. Dispersed Orthopedia (Otp) cells were found in the supracommissural and posterior nuclei of the ventral telencephalon, related to the medial portion of the amygdaloid complex. The preoptic area was identified by the Nkx2.1/Isl1 expression. In the alar hypothalamus, an Otp-expressing territory, lacking Nkx2.1/Isl1, was identified as the paraventricular domain. The adjacent subparaventricular domain (Spa) was subdivided in a rostral territory expressing Nkx2.1 and an Isl1+ caudal one. In the basal hypothalamus, the tuberal region was defined by the Nkx2.1/Isl1 expression and a rostral Otp-expressing domain was identified. Moreover, the Otp/Nkx2.1 combination showed an additional zone lacking Isl1, tentatively identified as the mamillary area. In the diencephalon, both Pax6 and Isl1 defined the prethalamic domain, and within the basal prosomere 3, scattered Pax6- and Isl1-expressing cells were observed in the posterior tubercle. Finally, a small group of Pax6 cells was observed in the pretectal area. These results improve the understanding of the forebrain evolution and demonstrate that its basic bauplan is present very early in the vertebrate lineage.

Exon-capture data and locus screening provide new insights into the phylogeny of flatfishes (Pleuronectoidei)

There is an extensive collection of literature on the taxonomy and phylogenetics of flatfishes (Pleuronectiformes) that extends over two centuries, but consensus on many of their evolutionary relationships remains elusive. Phylogenetic uncertainty stems from highly divergent results derived from morphological and genetic characters, and between various molecular datasets. Deciphering relationships is complicated by rapid diversification early in the Pleuronectiformes tree and an abundance of studies that incompletely and inconsistently sample taxa and genetic markers. We present phylogenies based on a genome-wide dataset (4,434 nuclear markers via exon-capture) and wide taxon sampling (86 species spanning 12 of 16 families) of the largest flatfish suborder (Pleuronectoidei). Nine different subsets of the data and two tree construction approaches (eighteen phylogenies in total) are remarkably consistent with other recent molecular phylogenies, and show strong support for the monophyly of all families included except Pleuronectidae. Analyses resolved a novel phylogenetic hypothesis for the family Rhombosoleidae as being within the Pleuronectoidea rather than the Soleoidea, and failed to support the subfamily Hippoglossinae as a monophyletic group. Our results were corroborated with evidence from previous phylogenetic studies to outline regions of persistent phylogenetic uncertainty and identify groups in need of further phylogenetic inference.

Holosteans contextualize the role of the teleost genome duplication in promoting the rise of evolutionary novelties in the ray-finned fish innate immune system

Over 99% of ray-finned fishes (Actinopterygii) are teleosts, a clade that comprises half of all living vertebrate species that have diversified across virtually all fresh and saltwater ecosystems. This ecological breadth raises the question of how the immunogenetic diversity required to persist under heterogeneous pathogen pressures evolved. The teleost genome duplication (TGD) has been hypothesized as the evolutionary event that provided the substrate for rapid genomic evolution and innovation. However, studies of putative teleost-specific innate immune receptors have been largely limited to comparisons either among teleosts or between teleosts and distantly related vertebrate clades such as tetrapods. Here we describe and characterize the receptor diversity of two clustered innate immune gene families in the teleost sister lineage: Holostei (bowfin and gars). Using genomic and transcriptomic data, we provide a detailed investigation of the phylogenetic history and conserved synteny of gene clusters encoding diverse immunoglobulin domain-containing proteins (DICPs) and novel immune-type receptors (NITRs). These data demonstrate an ancient linkage of DICPs to the major histocompatibility complex (MHC) and reveal an evolutionary origin of NITR variable-joining (VJ) exons that predate the TGD by at least 50 million years. Further characterizing the receptor diversity of Holostean DICPs and NITRs illuminates a sequence diversity that rivals the diversity of these innate immune receptor families in many teleosts. Taken together, our findings provide important historical context for the evolution of these gene families that challenge prevailing expectations concerning the consequences of the TGD during actinopterygiian evolution.

Excluding loci with substitution saturation improves inferences from phylogenomic data

The historical signal in nucleotide sequences becomes eroded over time by substitutions occurring repeatedly at the same sites. This phenomenon, known as substitution saturation, is recognized as one of the primary obstacles to deep-time phylogenetic inference using genome-scale data sets. We present a new test of substitution saturation and demonstrate its performance in simulated and empirical data. For some of the 36 empirical phylogenomic data sets that we examined, we detect substitution saturation in around 50% of loci. We found that saturation tends to be flagged as problematic in loci with highly discordant phylogenetic signals across sites. Within each data set, the loci with smaller numbers of informative sites are more likely to be flagged as containing problematic levels of saturation. The entropy saturation test proposed here is sensitive to high evolutionary rates relative to the evolutionary timeframe, while also being sensitive to several factors known to mislead phylogenetic inference, including short internal branches relative to external branches, short nucleotide sequences, and tree imbalance. Our study demonstrates that excluding loci with substitution saturation can be an effective means of mitigating the negative impact of multiple substitutions on phylogenetic inferences. [Phylogenetic model performance; phylogenomics; substitution model; substitution saturation; test statistics.]

The bowfin genome illuminates the developmental evolution of ray-finned fishes

The bowfin (Amia calva) is a ray-finned fish that possesses a unique suite of ancestral and derived phenotypes, which are key to understanding vertebrate evolution. The phylogenetic position of bowfin as a representative of neopterygian fishes, its archetypical body plan and its unduplicated and slowly evolving genome make bowfin a central species for the genomic exploration of ray-finned fishes. Here we present a chromosome-level genome assembly for bowfin that enables gene-order analyses, settling long-debated neopterygian phylogenetic relationships. We examine chromatin accessibility and gene expression through bowfin development to investigate the evolution of immune, scale, respiratory and fin skeletal systems and identify hundreds of gene-regulatory loci conserved across vertebrates. These resources connect developmental evolution among bony fishes, further highlighting the bowfin's importance for illuminating vertebrate biology and diversity in the genomic era.

A brief review of vertebrate sex evolution with a pledge for integrative research: towards 'sexomics'

Triggers and biological processes controlling male or female gonadal differentiation vary in vertebrates, with sex determination (SD) governed by environmental factors or simple to complex genetic mechanisms that evolved repeatedly and independently in various groups. Here, we review sex evolution across major clades of vertebrates with information on SD, sexual development and reproductive modes. We offer an up-to-date review of divergence times, species diversity, genomic resources, genome size, occurrence and nature of polyploids, SD systems, sex chromosomes, SD genes, dosage compensation and sex-biased gene expression. Advances in sequencing technologies now enable us to study the evolution of SD at broader evolutionary scales, and we now hope to pursue a sexomics integrative research initiative across vertebrates. The vertebrate sexome comprises interdisciplinary and integrated information on sexual differentiation, development and reproduction at all biological levels, from genomes, transcriptomes and proteomes, to the organs involved in sexual and sex-specific processes, including gonads, secondary sex organs and those with transcriptional sex-bias. The sexome also includes ontogenetic and behavioural aspects of sexual differentiation, including malfunction and impairment of SD, sexual differentiation and fertility. Starting from data generated by high-throughput approaches, we encourage others to contribute expertise to building understanding of the sexomes of many key vertebrate species. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part I)'.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Phylogenomic Subsampling and the Search for Phylogenetically Reliable Loci

Phylogenomic subsampling is a procedure by which small sets of loci are selected from large genome-scale data sets and used for phylogenetic inference. This step is often motivated by either computational limitations associated with the use of complex inference methods or as a means of testing the robustness of phylogenetic results by discarding loci that are deemed potentially misleading. Although many alternative methods of phylogenomic subsampling have been proposed, little effort has gone into comparing their behavior across different data sets. Here, I calculate multiple gene properties for a range of phylogenomic data sets spanning animal, fungal, and plant clades, uncovering a remarkable predictability in their patterns of covariance. I also show how these patterns provide a means for ordering loci by both their rate of evolution and their relative phylogenetic usefulness. This method of retrieving phylogenetically useful loci is found to be among the top performing when compared with alternative subsampling protocols. Relatively common approaches such as minimizing potential sources of systematic bias or increasing the clock-likeness of the data are found to fare worse than selecting loci at random. Likewise, the general utility of rate-based subsampling is found to be limited: loci evolving at both low and high rates are among the least effective, and even those evolving at optimal rates can still widely differ in usefulness. This study shows that many common subsampling approaches introduce unintended effects in off-target gene properties and proposes an alternative multivariate method that simultaneously optimizes phylogenetic signal while controlling for known sources of bias.

Integrating Sequence Capture and Restriction Site-Associated DNA Sequencing to Resolve Recent Radiations of Pelagic Seabirds

The diversification of modern birds has been shaped by a number of radiations. Rapid diversification events make reconstructing the evolutionary relationships among taxa challenging due to the convoluted effects of incomplete lineage sorting (ILS) and introgression. Phylogenomic data sets have the potential to detect patterns of phylogenetic incongruence, and to address their causes. However, the footprints of ILS and introgression on sequence data can vary between different phylogenomic markers at different phylogenetic scales depending on factors such as their evolutionary rates or their selection pressures. We show that combining phylogenomic markers that evolve at different rates, such as paired-end double-digest restriction site-associated DNA (PE-ddRAD) and ultraconserved elements (UCEs), allows a comprehensive exploration of the causes of phylogenetic discordance associated with short internodes at different timescales. We used thousands of UCE and PE-ddRAD markers to produce the first well-resolved phylogeny of shearwaters, a group of medium-sized pelagic seabirds that are among the most phylogenetically controversial and endangered bird groups. We found that phylogenomic conflict was mainly derived from high levels of ILS due to rapid speciation events. We also documented a case of introgression, despite the high philopatry of shearwaters to their breeding sites, which typically limits gene flow. We integrated state-of-the-art concatenated and coalescent-based approaches to expand on previous comparisons of UCE and RAD-Seq data sets for phylogenetics, divergence time estimation, and inference of introgression, and we propose a strategy to optimize RAD-Seq data for phylogenetic analyses. Our results highlight the usefulness of combining phylogenomic markers evolving at different rates to understand the causes of phylogenetic discordance at different timescales. [Aves; incomplete lineage sorting; introgression; PE-ddRAD-Seq; phylogenomics; radiations; shearwaters; UCEs.].

Macroevolutionary foundations of a recently evolved innate immune defense

Antagonistic interactions between hosts and parasites may drive the evolution of novel host defenses, or new parasite strategies. Host immunity is therefore one of the fastest evolving traits. But where do the novel immune traits come from? Here, we test for phylogenetic conservation in a rapidly evolving immune trait—peritoneal fibrosis. Peritoneal fibrosis is a costly defense against a specialist tapeworm, Schistocephalus solidus (Cestoda), expressed in some freshwater populations of threespine stickleback fish (Gasterosteus aculeatus, Perciformes). We asked whether stickleback fibrosis is a derived species‐specific trait or an ancestral immune response that was widely distributed across ray‐finned fish (Actinopterygii) only to be employed by threespine stickleback against the specialist parasite. We combined literature review on peritoneal fibrosis with a comparative experiment using either parasite‐specific, or nonspecific, immune challenge in deliberately selected species across fish tree of life. We show that ray‐finned fish are broadly, but not universally, able to induce peritoneal fibrosis when challenged with a generic stimulus (Alum adjuvant). The experimental species were, however, largely indifferent to the tapeworm antigen homogenate. Peritoneal fibrosis, thus, appears to be a common and deeply conserved fish immune response that was co‐opted by stickleback to adapt to a new selective challenge.

Comparing Ultraconserved Elements and Exons for Phylogenomic Analyses of Middle American Cichlids: When Data Agree to Disagree

Choosing among types of genomic markers to be used in a phylogenomic study can have a major influence on the cost, design, and results of a study. Yet few attempts have been made to compare categories of next-generation sequence markers limiting our ability to compare the suitability of these different genomic fragment types. Here, we explore properties of different genomic markers to find if they vary in the accuracy of component phylogenetic trees and to clarify the causes of conflict obtained from different data sets or inference methods. As a test case, we explore the causes of discordance between phylogenetic hypotheses obtained using a novel data set of ultraconserved elements (UCEs) and a recently published exon data set of the cichlid tribe Heroini. Resolving relationships among heroine cichlids has historically been difficult, and the processes of colonization and diversification in Middle America and the Greater Antilles are not yet well understood. Despite differences in informativeness and levels of gene tree discordance between UCEs and exons, the resulting phylogenomic hypotheses generally agree on most relationships. The independent data sets disagreed in areas with low phylogenetic signal that were overwhelmed by incomplete lineage sorting and nonphylogenetic signals. For UCEs, high levels of incomplete lineage sorting were found to be the major cause of gene tree discordance, whereas, for exons, nonphylogenetic signal is most likely caused by a reduced number of highly informative loci. This paucity of informative loci in exons might be due to heterogeneous substitution rates that are problematic to model (i.e., computationally restrictive) resulting in systematic errors that UCEs (being less informative individually but more uniform) are less prone to. These results generally demonstrate the robustness of phylogenomic methods to accommodate genomic markers with different biological and phylogenetic properties. However, we identify common and unique pitfalls of different categories of genomic fragments when inferring enigmatic phylogenetic relationships.

Seadragon genome analysis provides insights into its phenotype and sex determination locus

The iconic phenotype of seadragons includes leaf-like appendages, a toothless tubular mouth, and male pregnancy involving incubation of fertilized eggs on an open "brood patch." We de novo-sequenced male and female genomes of the common seadragon (Phyllopteryx taeniolatus) and its closely related species, the alligator pipefish (Syngnathoides biaculeatus). Transcription profiles from an evolutionary novelty, the leaf-like appendages, show that a set of genes typically involved in fin development have been co-opted as well as an enrichment of transcripts for potential tissue repair and immune defense genes. The zebrafish mutants for scpp5, which is lost in all syngnathids, were found to lack or have deformed pharyngeal teeth, supporting the hypothesis that the loss of scpp5 has contributed to the loss of teeth in syngnathids. A putative sex-determining locus encoding a male-specific amhr2y gene shared by common seadragon and alligator pipefish was identified.

A Y-linked anti-Müllerian hormone type-II receptor is the sex-determining gene in ayu, Plecoglossus altivelis

Whole-genome duplication and genome compaction are thought to have played important roles in teleost fish evolution. Ayu (or sweetfish), Plecoglossus altivelis, belongs to the superorder Stomiati, order Osmeriformes. Stomiati is phylogenetically classified as sister taxa of Neoteleostei. Thus, ayu holds an important position in the fish tree of life. Although ayu is economically important for the food industry and recreational fishing in Japan, few genomic resources are available for this species. To address this problem, we produced a draft genome sequence of ayu by whole-genome shotgun sequencing and constructed linkage maps using a genotyping-by-sequencing approach. Syntenic analyses of ayu and other teleost fish provided information about chromosomal rearrangements during the divergence of Stomiati, Protacanthopterygii and Neoteleostei. The size of the ayu genome indicates that genome compaction occurred after the divergence of the family Osmeridae. Ayu has an XX/XY sex-determination system for which we identified sex-associated loci by a genome-wide association study by genotyping-by-sequencing and whole-genome resequencing using wild populations. Genome-wide association mapping using wild ayu populations revealed three sex-linked scaffolds (total, 2.03 Mb). Comparison of whole-genome resequencing mapping coverage between males and females identified male-specific regions in sex-linked scaffolds. A duplicate copy of the anti-Müllerian hormone type-II receptor gene (amhr2bY) was found within these male-specific regions, distinct from the autosomal copy of amhr2. Expression of the Y-linked amhr2 gene was male-specific in sox9b-positive somatic cells surrounding germ cells in undifferentiated gonads, whereas autosomal amhr2 transcripts were detected in somatic cells in sexually undifferentiated gonads of both genetic males and females. Loss-of-function mutation for amhr2bY induced male to female sex reversal. Taken together with the known role of Amh and Amhr2 in sex differentiation, these results indicate that the paralog of amhr2 on the ayu Y chromosome determines genetic sex, and the male-specific amh-amhr2 pathway is critical for testicular differentiation in ayu.

Phylogenetic analyses of ray-finned fishes (Actinopterygii) using collagen type I protein sequences

Ray-finned fishes (Actinopterygii) are the largest and most diverse group of vertebrates, comprising over half of all living vertebrate species. Phylogenetic relationships between ray-finned fishes have historically pivoted on the study of morphology, which has notoriously failed to resolve higher order relationships, such as within the percomorphs. More recently, comprehensive genomic analyses have provided further resolution of actinopterygian phylogeny, including higher order relationships. Such analyses are rightfully regarded as the 'gold standard' for phylogenetics. However, DNA retrieval requires modern or well-preserved tissue and is less likely to be preserved in archaeological or fossil specimens. By contrast, some proteins, such as collagen, are phylogenetically informative and can survive into deep time. Here, we test the utility of collagen type I amino acid sequences for phylogenetic estimation of ray-finned fishes. We estimate topology using Bayesian approaches and compare the congruence of our estimated trees with published genomic phylogenies. Furthermore, we apply a Bayesian molecular clock approach and compare estimated divergence dates with previously published genomic clock analyses. Our collagen-derived trees exhibit 77% of node positions as congruent with recent genomic-derived trees, with the majority of discrepancies occurring in higher order node positions, almost exclusively within the Percomorpha. Our molecular clock trees present divergence times that are fairly comparable with genomic-based phylogenetic analyses. We estimate the mean node age of Actinopteri at ∼293 million years (Ma), the base of Teleostei at ∼211 Ma and the radiation of percomorphs beginning at ∼141 Ma (∼350 Ma, ∼250-283 Ma and ∼120-133 Ma in genomic trees, respectively). Finally, we show that the average rate of collagen (I) sequence evolution is 0.9 amino acid substitutions for every million years of divergence, with the α3 (I) sequence evolving the fastest, followed by the α2 (I) chain. This is the quickest rate known for any vertebrate group. We demonstrate that phylogenetic analyses using collagen type I amino acid sequences generate tangible signals for actinopterygians that are highly congruent with recent genomic-level studies. However, there is limited congruence within percomorphs, perhaps due to clade-specific functional constraints acting upon collagen sequences. Our results provide important insights for future phylogenetic analyses incorporating extinct actinopterygian species via collagen (I) sequencing.

Changes in Nkx2.1, Sox2, Bmp4, and Bmp16 expression underlying the lung-to-gas bladder evolutionary transition in ray-finned fishes

The key to understanding the evolutionary origin and modification of phenotypic traits is revealing the responsible underlying developmental genetic mechanisms. An important organismal trait of ray-finned fishes is the gas bladder, an air-filled organ that, in most fishes, functions for buoyancy control, and is homologous to the lungs of lobe-finned fishes. The critical morphological difference between lungs and gas bladders, which otherwise share many characteristics, is the general direction of budding during development. Lungs bud ventrally and the gas bladder buds dorsally from the anterior foregut. We investigated the genetic underpinnings of this ventral-to-dorsal shift in budding direction by studying the expression patterns of known lung genes (Nkx2.1, Sox2, and Bmp4) during the development of lungs or gas bladder in three fishes: bichir, bowfin, and zebrafish. Nkx2.1 and Sox2 show reciprocal dorsoventral expression patterns during tetrapod lung development and are important regulators of lung budding; their expression during bichir lung development is conserved. Surprisingly, we find during gas bladder development, Nkx2.1 and Sox2 expression are inconsistent with the hypothesis that they regulate the direction of gas bladder budding. Bmp4 is expressed ventrally during lung development in bichir, akin to the pattern during mouse lung development. During gas bladder development, Bmp4 is not expressed. However, Bmp16, a paralogue of Bmp4, is expressed dorsally in the developing gas bladder of bowfin. Bmp16 is present in the known genomes of Actinopteri (ray-finned fishes excluding bichir) but absent from mammalian genomes. We hypothesize that Bmp16 was recruited to regulate gas bladder development in the Actinopteri in place of Bmp4.

Evolutionary Rate Variation Among Lineages in Gene Trees has a Negative Impact on Species-Tree Inference

Phylogenetic analyses of genomic data provide a powerful means of reconstructing the evolutionary relationships among organisms, yet such analyses are often hindered by conflicting phylogenetic signals among loci. Identifying the signals that are most influential to species-tree estimation can help to inform the choice of data for phylogenomic analysis. We investigated this in an analysis of 30 phylogenomic data sets. For each data set, we examined the association between several branch-length characteristics of gene trees and the distance between these gene trees and the corresponding species trees. We found that the distance of each gene tree to the species tree inferred from the full data set was positively associated with variation in root-to-tip distances and negatively associated with mean branch support. However, no such associations were found for gene-tree length, a measure of the overall substitution rate at each locus. We further explored the usefulness of the best-performing branch-based characteristics for selecting loci for phylogenomic analyses. We found that loci that yield gene trees with high variation in root-to-tip distances have a disproportionately distant signal of tree topology compared with the complete data sets. These results suggest that rate variation across lineages should be taken into consideration when exploring and even selecting loci for phylogenomic analysis.[Branch support; data filtering; nucleotide substitution model; phylogenomics; substitution rate; summary coalescent methods.]

Orthologous Divergence and Paralogous Anticonvergence in Molecular Evolution of Triplicated Green Opsin Genes in Medaka Fish, Genus Oryzias

Gene duplication of green (RH2) opsin genes and their spectral differentiation are well documented in many teleost fish. However, their evolutionary divergence or conservation patterns among phylogenetically close but ecologically diverse species is not well explored. Medaka fish (genus Oryzias) are broadly distributed in fresh and brackish waters of Asia, with many species being laboratory-housed and feasible for genetic studies. We previously showed that a Japan strain (HNI) of medaka (Oryzias latipes) possessed three RH2 opsin genes (RH2-A, RH2-B, and RH2-C) encoding spectrally divergent photopigments. Here, we examined the three RH2 opsin genes from six Oryzias species representing three species groups: the latipes, the celebensis, and the javanicus. Photopigment reconstitution revealed that the peak absorption spectra (λ_max) of RH2-A were divergent among the species (447–469 nm), whereas those of RH2-B and RH2-C were conservative (516–519 and 486–493 nm, respectively). For the RH2-A opsins, the largest spectral shift was detected in the phylogenetic branch leading to the latipes group. A single amino acid replacement T94C explained most of the spectral shift. For RH2-B and -C opsins, we detected tracts of gene conversion between the two genes homogenizing them. Nevertheless, several amino acid differences were maintained. We showed that the spectral difference between the two opsins was attributed to largely the E/Q amino acid difference at the site 122 and to several sites with individually small spectral effects. These results depict dynamism of spectral divergence of orthologous and paralogous green opsin genes in phylogenetically close but ecologically diverse species exemplified by medaka.

Comparing diversification rates in lakes, rivers, and the sea

The diversity of species inhabiting freshwater relative to marine habitats is striking, given that freshwater habitats encompass <1% of Earth's water. The most commonly proposed explanation for this pattern is that freshwater habitats are more fragmented than marine habitats, allowing more opportunities for allopatric speciation and thus increased diversification rates in freshwater. However, speciation may be generally faster in sympatry than in allopatry, as illustrated by lacustrine radiations such as African cichlids. Such differences between rivers and lakes may be important to consider when comparing diversification broadly among freshwater and marine groups. Here I compared diversification rates of teleost fishes in marine, riverine and lacustrine habitats. I found that lakes had faster speciation and net diversification rates than other aquatic habitats. However, most freshwater diversity arose in rivers. Surprisingly, riverine and marine habitats had similar rates of net diversification on average. Biogeographic models suggest that lacustrine habitats are evolutionarily unstable, explaining the dearth of lacustrine species in spite of their rapid diversification. Collectively, these results suggest that strong diversification rate differences are unlikely to explain the freshwater paradox. Instead, this pattern may be attributable to the comparable amount of time spent in riverine and marine habitats over the 200-million-year history of teleosts.

Deep in the systematics of Camallanidae (Nematoda): using integrative taxonomy to better understand the phylogeny and consistency of diagnostic traits

Due to conflicts between classic and molecular systematics of Camallanidae, different data types were used for the first time, to better understand the evolutionary history and taxa consistency within this family. Genetic [18S and 28S rDNA; cytochrome c oxidase subunit I (COI) mtDNA], morphological and life history traits were used to infer phylogenies using Bayesian inference, reconstructed from separated and concatenated datasets. The consistency of tree and morphological traits was evaluated using the consistency index. Characters were mapped on the trees and the phylogenetic informativeness of genetic markers was estimated. Phylogenetic informativeness of 18S provided better resolution for outer nodes, COI for inners and 28S had an intermediate profile. New sequences for two camallanid species were obtained. Phylogenies of genetic and concatenated data largely agreed, showing more divergence in the COI dataset, due to its higher mutation rate vs stable morphology for diagnosing higher taxa. No genus sustained monophyly. The lack of autapomorphy and phylogenetic proximity supported the partition of Batrachocamallanus as synonym of Procamallanus and Spirocamallanus, which should not be considered as subgenera. Although traits of buccal capsule, male tail, habitat, host and biogeographic were highly consistent, intrinsic patterns varied according to different taxa assemblages. Morphological systematics of Camallanidae, based on buccal capsule, is artificial for certain taxa.

Repeated evolution of circadian clock dysregulation in cavefish populations

Circadian rhythms are nearly ubiquitous throughout nature, suggesting they are critical for survival in diverse environments. Organisms inhabiting largely arrhythmic environments, such as caves, offer a unique opportunity to study the evolution of circadian rhythms in response to changing ecological pressures. Populations of the Mexican tetra, Astyanax mexicanus, have repeatedly invaded caves from surface rivers, where individuals must contend with perpetual darkness, reduced food availability, and limited fluctuations in daily environmental cues. To investigate the molecular basis for evolved changes in circadian rhythms, we investigated rhythmic transcription across multiple independently-evolved cavefish populations. Our findings reveal that evolution in a cave environment has led to the repeated disruption of the endogenous biological clock, and its entrainment by light. The circadian transcriptome shows widespread reductions and losses of rhythmic transcription and changes to the timing of the activation/repression of core-transcriptional clock. In addition to dysregulation of the core clock, we find that rhythmic transcription of the melatonin regulator aanat2 and melatonin rhythms are disrupted in cavefish under darkness. Mutants of aanat2 and core clock gene rorca disrupt diurnal regulation of sleep in A. mexicanus, phenocopying circadian modulation of sleep and activity phenotypes of cave populations. Together, these findings reveal multiple independent mechanisms for loss of circadian rhythms in cavefish populations and provide a platform for studying how evolved changes in the biological clock can contribute to variation in sleep and circadian behavior.

Temporal dynamics of teleost populations during the Pleistocene: a report from publicly available genome data

Background

Global climate oscillation, as a selection dynamic, is an ecologically important element resulting in global biodiversity. During the glacial geological periods, most organisms suffered detrimental selection pressures (such as food shortage and habitat loss) and went through population declines. However, during the mild interglacial periods, many species re-flourished. These temporal dynamics of effective population sizes (N_e) provide essential information for understanding and predicting evolutionary outcomes during historical and ongoing global climate changes.

Results

Using high-quality genome assemblies and corresponding sequencing data, we applied the Pairwise Sequentially Markovian Coalescent (PSMC) method to quantify N_e changes of twelve representative teleost species from approximately 10 million years ago (mya) to 10 thousand years ago (kya). These results revealed multiple rounds of population contraction and expansion in most of the examined teleost species during the Neogene and the Quaternary periods. We observed that 83% (10/12) of the examined teleosts had experienced a drastic decline in N_e before the last glacial period (LGP, 110–12 kya), slightly earlier than the reported pattern of N_e changes in 38 avian species. In comparison with the peaks, almost all of the examined teleosts maintained long-term lower N_e values during the last few million years. This is consistent with increasingly dramatic glaciation during this period.

Conclusion

In summary, these findings provide a more comprehensive understanding of the historical N_e changes in teleosts. Results presented here could lead to the development of appropriate strategies to protect species in light of ongoing global climate changes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07816-7.

Evolutionary Modifications Are Moderate in the Astroglial System of Actinopterygii as Revealed by GFAP Immunohistochemistry

The present paper is the first comparative study on the astroglia of several actinopterygian species at different phylogenetical positions, teleosts (16 species), and non-teleosts (3 species), based on the immunohistochemical staining of GFAP (glial fibrillary acidic protein), the characteristic cytoskeletal intermediary filament protein, and immunohistochemical marker of astroglia. The question was, how the astroglial architecture reflexes the high diversity of this largest vertebrate group. The actinopterygian telencephalon has a so-called ‘eversive’ development in contrast to the ‘evagination’ found in sarcopterygii (including tetrapods). Several brain parts either have no equivalents in tetrapod vertebrates (e.g., torus longitudinalis, lobus inferior, lobus nervi vagi), or have rather different shapes (e.g., the cerebellum). GFAP was visualized applying DAKO polyclonal anti-GFAP serum. The study was focused mainly on the telencephalon (eversion), tectum (visual orientation), and cerebellum (motor coordination) where the evolutionary changes were most expected, but the other areas were also investigated. The predominant astroglial elements were tanycytes (long, thin, fiber-like cells). In the teleost telencephala a ‘fan-shape’ re-arrangement of radial glia reflects the eversion. In bichir, starlet, and gar, in which the eversion is less pronounced, the ‘fan-shape’ re-arrangement did not form. In the tectum the radial glial processes were immunostained, but in Ostariophysi and Euteleostei it did not extend into their deep segments. In the cerebellum Bergmann-like glia was found in each group, including non-teleosts, except for Cyprinidae. The vagal lobe was uniquely enlarged and layered in Cyprininae, and had a corresponding layered astroglial system, which left almost free of GFAP the zones of sensory and motor neurons. In conclusion, despite the diversity and evolutionary alterations of Actinopterygii brains, the diversity of the astroglial architecture is moderate. In contrast to Chondrichthyes and Amniotes; in Actinopterygii true astrocytes (stellate-shaped extraependymal cells) did not appear during evolution, and the expansion of GFAP-free areas was limited.

Accelerated Diversification Explains the Exceptional Species Richness of Tropical Characoid Fishes

The Neotropics harbor the most species-rich freshwater fish fauna on the planet, but the timing of that exceptional diversification remains unclear. Did the Neotropics accumulate species steadily throughout their long history, or attain their remarkable diversity recently? Biologists have long debated the relative support for these museum and cradle hypotheses, but few phylogenies of megadiverse tropical clades have included sufficient taxa to distinguish between them. We used 1,288 ultraconserved element loci (UCE) spanning 293 species, 211 genera and 21 families of characoid fishes to reconstruct a new, fossil-calibrated phylogeny and infer the most likely diversification scenario for a clade that includes a third of Neotropical fish diversity. This phylogeny implies paraphyly of the traditional delimitation of Characiformes because it resolves the largely Neotropical Characoidei as the sister lineage of Siluriformes (catfishes), rather than the African Citharinodei. Time-calibrated phylogenies indicate an ancient origin of major characoid lineages and reveal a much more recent emergence of most characoid species. Diversification rate analyses infer increased speciation and decreased extinction rates during the Oligocene at around 30 million years ago (Ma) during a period of mega-wetland formation in the proto-Orinoco-Amazonas. Three species-rich and ecomorphologically diverse lineages (Anostomidae, Serrasalmidae, and Characidae) that originated more than 60 Ma in the Paleocene experienced particularly notable bursts of Oligocene diversification and now account collectively for 68% of the approximately 2,150 species of Characoidei. In addition to paleogeographic changes, we discuss potential accelerants of diversification in these three lineages. While the Neotropics accumulated a museum of ecomorphologically diverse characoid lineages long ago, this geologically dynamic region also cradled a much more recent birth of remarkable species-level diversity.

Contingency and chance erase necessity in the experimental evolution of ancestral proteins

The roles of chance, contingency, and necessity in evolution are unresolved because they have never been assessed in a single system or on timescales relevant to historical evolution. We combined ancestral protein reconstruction and a new continuous evolution technology to mutate and select proteins in the B-cell lymphoma-2 (BCL-2) family to acquire protein–protein interaction specificities that occurred during animal evolution. By replicating evolutionary trajectories from multiple ancestral proteins, we found that contingency generated over long historical timescales steadily erased necessity and overwhelmed chance as the primary cause of acquired sequence variation; trajectories launched from phylogenetically distant proteins yielded virtually no common mutations, even under strong and identical selection pressures. Chance arose because many sets of mutations could alter specificity at any timepoint; contingency arose because historical substitutions changed these sets. Our results suggest that patterns of variation in BCL-2 sequences – and likely other proteins, too – are idiosyncratic products of a particular and unpredictable course of historical events.

One of the most fundamental and unresolved questions in evolutionary biology is whether the outcomes of evolution are predictable. Is the diversity of life we see today the expected result of organisms adapting to their environment throughout history (also known as natural selection) or the product of random chance? Or did chance events early in history shape the paths that evolution could take next, determining the biological forms that emerged under natural selection much later?

These questions are hard to study because evolution happened only once, long ago. To overcome this barrier, Xie, Pu, Metzger et al. developed an experimental approach that can evolve reconstructed ancestral proteins that existed deep in the past. Using this method, it is possible to replay evolution multiple times, from various historical starting points, under conditions similar to those that existed long ago. The end products of the evolutionary trajectories can then be compared to determine how predictable evolution actually is.

Xie, Pu, Metzger et al. studied proteins belonging to the BCL-2 family, which originated some 800 million years ago. These proteins have diversified greatly over time in both their genetic sequences and their ability to bind to specific partner proteins called co-regulators. Xie, Pu, Metzger et al. synthesized BCL-2 proteins that existed at various times in the past. Each ancestral protein was then allowed to evolve repeatedly under natural selection to acquire the same co-regulator binding functions that evolved during history.

At the end of each evolutionary trajectory, the genetic sequence of the resulting BCL-2 proteins was recorded. This revealed that the outcomes of evolution were almost completely unpredictable: trajectories initiated from the same ancestral protein produced proteins with very different sequences, and proteins launched from different ancestral starting points were even more dissimilar.

Further experiments identified the mutations in each trajectory that caused changes in coregulator binding. When these mutations were introduced into other ancestral proteins, they did not yield the same change in function. This suggests that early chance events influenced each protein’s evolution in an unpredictable way by opening and closing the paths available to it in the future.

This research expands our understanding of evolution on a molecular level whilst providing a new experimental approach for studying evolutionary drivers in more detail. The results suggest that BCL-2 proteins, in all their various forms, are unique products of a particular, unpredictable course of history set in motion by ancient chance events.

The deep sea is a hot spot of fish body shape evolution

Deep-sea fishes have long captured our imagination with striking adaptations to life in the mysterious abyss, raising the possibility that this cold, dark ocean region may be a key hub for physiological and functional diversification. We explore this idea through an analysis of body shape evolution across ocean depth zones in over 3000 species of marine teleost fishes. We find that the deep ocean contains twice the body shape disparity of shallow waters, driven by elevated rates of evolution in traits associated with locomotion. Deep-sea fishes display more frequent adoption of forms suited to slow and periodic swimming, whereas shallow living species are concentrated around shapes conferring strong, sustained swimming capacity and manoeuvrability. Our results support long-standing impressions of the deep sea as an evolutionary hotspot for fish body shape evolution and highlight that factors like habitat complexity and ecological interactions are potential drivers of this adaptive diversification.

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.

Whole genome sequencing of a snailfish from the Yap Trench (~7,000 m) clarifies the molecular mechanisms underlying adaptation to the deep sea

Hadal environments (depths below 6,000 m) are characterized by extremely high hydrostatic pressures, low temperatures, a scarce food supply, and little light. The evolutionary adaptations that allow vertebrates to survive in this extreme environment are poorly understood. Here, we constructed a high-quality reference genome for Yap hadal snailfish (YHS), which was captured at a depth of ~7,000 m in the Yap Trench. The final YHS genome assembly was 731.75 Mb, with a contig N50 of 0.75 Mb and a scaffold N50 of 1.26 Mb. We predicted 24,329 protein-coding genes in the YHS genome, and 24,265 of these genes were successfully functionally annotated. Phylogenetic analyses suggested that YHS diverged from a Mariana Trench snailfish approximately 0.92 million years ago. Many genes associated with DNA repair show evidence of positive selection and have expanded copy numbers in the YHS genome, possibly helping to maintain the integrity of DNA under increased hydrostatic pressure. The levels of trimethylamine N-oxide (TMAO), a potent protein stabilizer, are much higher in the muscles of YHS than in those of shallow-water fish. This difference is perhaps due to the five copies of the TMAO-generating enzyme flavin-containing monooxygenase-3 gene (fmo3) in the YHS genome and the abundance of trimethylamine (TMA)-generating bacteria in the YHS gut. Thus, the high TMAO content might help YHS adapt to high hydrostatic pressure by improving protein stability. Additionally, the evolutionary features of the YHS genes encoding sensory-related proteins are consistent with the scarce food supply and darkness in the hadal environments. These results clarify the molecular mechanisms underlying the adaptation of hadal organisms to the deep-sea environment and provide valuable genomic resources for in-depth investigations of hadal biology.

Hadal environments (depths below 6,000 m) are characterized by extremely high hydrostatic pressures, low temperatures, a scarce food supply, and little light. Fish are the only vertebrates inhabiting the hadal zone, and hadal snailfishes have been found in at least five geographically separated marine trenches. However, the genetic mechanisms that allow vertebrates to live in such extreme conditions are not well understood. Here, we constructed a high-quality reference genome for Yap hadal snailfish (YHS) captured at a depth of ~7,000 m in the Yap Trench, using long reads obtained by Pacific Biosciences Sequel sequencing. Comparative genomic analyses revealed that many genes associated with DNA repair show evidence of positive selection and have expanded copy numbers in the YHS genome, which potentially reflect the difficulty of maintaining DNA integrity under high hydrostatic pressure. Moreover, the five copies of the trimethylamine N-oxide (TMAO)-generating enzyme flavin-containing monooxygenase-3 gene (fmo3) and the abundance of trimethylamine (TMA)-generating bacteria in the YHS gut could provide enough TMAO to improve protein stability under hadal conditions. In addition, characteristics of the YHS sensory system genes were consistent with the scarce food supply and darkness in the hadal zone. Our results provide new insights into the molecular mechanisms underlying the adaptation of hadal organisms to the deep-sea environment and valuable genomic resources that will help further clarify hadal adaptations.

Chromosome-level assembly of southern catfish (silurus meridionalis) provides insights into visual adaptation to nocturnal and benthic lifestyles

The Southern catfish (Silurus meridionalis) is a nocturnal and benthic freshwater fish endemic to the Yangtze River and its tributaries. In this study, we constructed a chromosome-level draft genome of S. meridionalis using 69.7-Gb Nanopore long reads and 49.5-Gb Illumina short reads. The genome assembly was 741.2 Mb in size with a contig N50 of 13.19 Mb. An additional 116.4 Gb of Bionano and 77.4 Gb of Hi-C data were applied to assemble contigs into scaffolds and further into 29 chromosomes, resulting in a 738.9-Mb genome with a scaffold N50 of 28.04 Mb. A total of 22,965 protein-coding genes were predicted from the genome with 22,519 (98.06%) genes functionally annotated. Comparative genomic and transcriptomic analyses revealed a rod-dominated visual system which was responsible for scotopic vision. The absence of cone opsins SWS1 and SWS2 resulted in the lack of ultraviolet and blue violet sensitivity. Mutations at key amino acid sites of RH1.1, RH1.2 and RH2 resulted in spectral tuning good for dim light vision and narrow colour vision. A higher expression level of rod phototransduction genes than that of cone genes and higher rod-to-cone ratio led to higher optical sensitivity under dim light conditions. In addition, analysis of the genes involved in eye morphogenesis and development revealed the loss of some conserved noncoding elements, which might be associated with the small eyes in catfish. Together, our study provides important clues for the adaptation of the catfish visual system to the nocturnal and benthic lifestyles. The draft genome of S. meridionalis represents a valuable resource for studies of the molecular mechanisms of ecological adaptation.

Bilateral visual projections exist in non-teleost bony fish and predate the emergence of tetrapods

In most vertebrates, camera-style eyes contain retinal ganglion cell neurons that project to visual centers on both sides of the brain. However, in fish, ganglion cells were thought to innervate only the contralateral side, suggesting that bilateral visual projections appeared in tetrapods. Here we show that bilateral visual projections exist in non-teleost fishes and that the appearance of ipsilateral projections does not correlate with terrestrial transition or predatory behavior. We also report that the developmental program that specifies visual system laterality differs between fishes and mammals, as the Zic2 transcription factor, which specifies ipsilateral retinal ganglion cells in tetrapods, appears to be absent from fish ganglion cells. However, overexpression of human ZIC2 induces ipsilateral visual projections in zebrafish. Therefore, the existence of bilateral visual projections likely preceded the emergence of binocular vision in tetrapods.

Development of the caudal-fin skeleton reveals multiple convergent fusions within Atherinomorpha

Background

The caudal fin of teleosts is a highly diverse morphological structure and a valuable source of information for comparative analyses. Within the Atherinomorpha a high variation of conditions of the caudal-fin skeleton can be found. These range from complex but basal configurations to simple yet derived configurations. When comparing atherinomorph taxa, it is often difficult to decide on the homology of skeletal elements if only considering adult specimens. However, observing the development of caudal-fin skeletons allows one to evaluate complex structures, reveal homologies and developmental patterns, and even reconstruct the grundplan of the examined taxa.

Results

We studied the development of the caudal-fin skeleton in different atheriniform, beloniform and cyprinodontiform species using cleared and stained specimens. Subsequently we compared the development to find similarities and differences in terms of 1) which structures are formed and 2) which structures fuse during ontogeny. For many structures, i.e., the parhypural, the epural(s), the haemal and neural spines of the preural centra and the uroneural, there were either no or only minor differences visible between the three taxa. However, the development of the hypurals revealed a high variation of fusions within different taxa that partly occurred independently in atheriniforms, beloniforms and cyprinodontiforms. Moreover, comparing the development of the ural centra exposed two ways of formation of the compound centrum: 1) in atheriniforms and the beloniforms Oryzias and Hyporhamphus limbatus two ural centra develop and fuse during ontogeny while 2) in cyprinodontiforms and Exocoetidae (Beloniformes) only a single ural centrum is formed during ontogeny.

Conclusions

We were able to reconstruct the grundplan of the developmental pattern of the caudal-fin skeleton of the Atheriniformes, Beloniformes and Cyprinodontiformes as well as their last common ancestors. We found two developmental modes of the compound centrum within the Atherinomorpha, i.e., the fusion of two developing ural centra in atheriniforms and beloniforms and the development of only one ural centrum in cyprinodontiforms. Further differences and similarities for the examined taxa are discussed, resulting in the hypothesis that the caudal-fin development of a last common ancestor to all atherinomorphs is very much similar to that of extant atheriniforms.

Phylogenomics and Historical Biogeography of Seahorses, Dragonets, Goatfishes, and Allies (Teleostei: Syngnatharia): Assessing Factors Driving Uncertainty in Biogeographic Inferences

The charismatic trumpetfishes, goatfishes, dragonets, flying gurnards, seahorses, and pipefishes encompass a recently defined yet extraordinarily diverse clade of percomorph fishes-the series Syngnatharia. This group is widely distributed in tropical and warm-temperate regions, with a great proportion of its extant diversity occurring in the Indo-Pacific. Because most syngnatharians feature long-range dispersal capabilities, tracing their biogeographic origins is challenging. Here, we applied an integrative phylogenomic approach to elucidate the evolutionary biogeography of syngnatharians. We built upon a recently published phylogenomic study that examined ultraconserved elements by adding 62 species (total 169 species) and one family (Draconettidae), to cover ca. 25% of the species diversity and all 10 families in the group. We inferred a set of time-calibrated trees and conducted ancestral range estimations. We also examined the sensitivity of these analyses to phylogenetic uncertainty (estimated from multiple genomic subsets), area delimitation, and biogeographic models that include or exclude the jump-dispersal parameter (j). Of the three factors examined, we found that the j parameter has the strongest effect in ancestral range estimates, followed by number of areas defined, and tree topology and divergence times. After accounting for these uncertainties, our results reveal that syngnatharians originated in the ancient Tethys Sea ca. 87 Ma (84-94 Ma; Late Cretaceous) and subsequently occupied the Indo-Pacific. Throughout syngnatharian history, multiple independent lineages colonized the eastern Pacific (6-8 times) and the Atlantic (6-14 times) from their center of origin, with most events taking place following an east-to-west route prior to the closure of the Tethys Seaway ca. 12-18 Ma. Ultimately, our study highlights the importance of accounting for different factors generating uncertainty in macroevolutionary and biogeographic inferences.

A fossil-calibrated time-tree of all Australian freshwater fishes

Australian freshwater fishes are a relatively species-poor assemblage, mostly comprising groups derived from older repeated freshwater invasions by marine ancestors, plus a small number of Gondwanan lineages. These taxa are both highly endemic and highly threatened, but a comprehensive phylogeny for Australian freshwater fishes is lacking. This has hampered efforts to study their phylogenetic diversity, distribution of extinction risk, speciation rates, and rates of trait evolution. Here, we present a comprehensive dated phylogeny of 412 Australian fishes. We include all formally recognized freshwater species plus a number of genetically distinct subpopulations, species awaiting formal description, and predominantly brackish-water species that sometimes enter fresh water. The phylogeny was inferred using maximum-likelihood analysis of a multilocus data set comprising six mitochondrial and three nuclear genes from 326 taxa. We inferred the evolutionary timescale using penalized likelihood, then used a statistical approach to add 86 taxa for which no molecular data were available. The time-tree inferred in our study will provide a useful resource for macroecological studies of Australian freshwater fishes by enabling corrections for phylogenetic non-independence in evolutionary and ecological comparative analyses.

The American Paddlefish Genome Provides Novel Insights into Chromosomal Evolution and Bone Mineralization in Early Vertebrates

Sturgeons and paddlefishes (Acipenseriformes) occupy the basal position of ray-finned fishes, although they have cartilaginous skeletons as in Chondrichthyes. This evolutionary status and their morphological specializations make them a research focus, but their complex genomes (polyploidy and the presence of microchromosomes) bring obstacles and challenges to molecular studies. Here, we generated the first high-quality genome assembly of the American paddlefish (Polyodon spathula) at a chromosome level. Comparative genomic analyses revealed a recent species-specific whole-genome duplication event, and extensive chromosomal changes, including head-to-head fusions of pairs of intact, large ancestral chromosomes within the paddlefish. We also provide an overview of the paddlefish SCPP (secretory calcium-binding phosphoprotein) repertoire that is responsible for tissue mineralization, demonstrating that the earliest flourishing of SCPP members occurred at least before the split between Acipenseriformes and teleosts. In summary, this genome assembly provides a genetic resource for understanding chromosomal evolution in polyploid nonteleost fishes and bone mineralization in early vertebrates.

Chromosome-level genome assembly of burbot (Lota lota) provides insights into the evolutionary adaptations in freshwater

The burbot (Lota lota) is the only member of the order Gadiformes adapted solely to freshwater. This species has the widest longitudinal range among freshwater fish worldwide. Burbot serves as a good model for studies on adaptive genome evolution from marine to freshwater environments. However, a high-quality reference genome of burbot has not yet been released. Here, the first chromosome-level genome of burbot was constructed using PacBio long sequencing and Hi-C technology. A total of 95.24 Gb polished PacBio sequences were generated, and the preliminary genome assembly was 575.83 Mb in size with a contig N50 size of 2.15 Mb. The assembled sequences were anchored to 22 pseudochromosomes by using Hi-C data. The final assembled genome after Hi-C correction was 575.92 Mb, with a contig N50 of 2.01 Mb and a scaffold N50 of 22.10 Mb. A total of 22,067 protein-coding genes were predicted, 94.82% of which were functionally annotated. Phylogenetic analyses indicated that burbot diverged with the Atlantic cod approximately 43.8 million years ago. In addition, 377 putative genes that appear to be under positive selection in burbot were identified. These positively selected genes might be involved in the adaptation to the freshwater environment. These genome data provide an invaluable resource for the ecological and evolutionary study of the order Gadiformes.

Relative Model Fit Does Not Predict Topological Accuracy in Single-Gene Protein Phylogenetics

It is regarded as best practice in phylogenetic reconstruction to perform relative model selection to determine an appropriate evolutionary model for the data. This procedure ranks a set of candidate models according to their goodness of fit to the data, commonly using an information theoretic criterion. Users then specify the best-ranking model for inference. Although it is often assumed that better-fitting models translate to increase accuracy, recent studies have shown that the specific model employed may not substantially affect inferences. We examine whether there is a systematic relationship between relative model fit and topological inference accuracy in protein phylogenetics, using simulations and real sequences. Simulations employed site-heterogeneous mechanistic codon models that are distinct from protein-level phylogenetic inference models, allowing us to investigate how protein models performs when they are misspecified to the data, as will be the case for any real sequence analysis. We broadly find that phylogenies inferred across models with vastly different fits to the data produce highly consistent topologies. We additionally find that all models infer similar proportions of false-positive splits, raising the possibility that all available models of protein evolution are similarly misspecified. Moreover, we find that the parameter-rich GTR (general time reversible) model, whose amino acid exchangeabilities are free parameters, performs similarly to models with fixed exchangeabilities, although the inference precision associated with GTR models was not examined. We conclude that, although relative model selection may not hinder phylogenetic analysis on protein data, it may not offer specific predictable improvements and is not a reliable proxy for accuracy.

Testing the utility of alternative metrics of branch support to address the ancient evolutionary radiation of tunas, stromateoids, and allies (Teleostei: Pelagiaria)

The use of high-throughput sequencing technologies to produce genome-scale datasets was expected to settle some long-standing controversies across the Tree of Life, particularly in areas where short branches occur at deep timescales. Instead, these datasets have often yielded many well-supported but conflicting topologies, and highly variable gene-tree distributions. A variety of branch-support metrics beyond the nonparametric bootstrap are now available to assess how robust a phylogenetic hypothesis may be, as well as new methods to quantify gene-tree discordance. We applied multiple branch support metrics to an ancient group of marine fishes (Teleostei: Pelagiaria) whose interfamilial relationships have proven difficult to resolve due to a rapid accumulation of lineages very early in its history. We analyzed hundreds of loci including published UCE data and newly generated exonic data along with their flanking regions to represent all 16 extant families for more than 150 out of 284 valid species in the group. Branch support was lower for interfamilial relationships (except the SH-like aLRT and aBayes methods) regardless of the type of marker used. Several nodes that were highly supported with bootstrap had very low site and gene-tree concordance, revealing underlying conflict. Despite this conflict, we were able to identify four consistent interfamilial clades, each comprised of two or three families. Combining exons with their flanking regions also produced increased branch lengths in the deep branches of the pelagiarian tree. Our results demonstrate the limitations of employing current metrics of branch support and species-tree estimation when assessing the confidence of ancient evolutionary radiations and emphasize the necessity to embrace alternative measurements to explore phylogenetic uncertainty and discordance in phylogenomic datasets.

Fish models for investigating nutritional regulation of embryonic development

In recent decades, biologist have focused on the spatiotemporal regulation and function of genes to understand embryogenesis. It is clear that maternal diet impacts fetal development but how nutrients, like lipids and vitamins, modify developmental programs is not completely understood. Fish are useful research organisms for such investigations. Most species of fish produce eggs that develop outside the mother, dependent on a finite amount of yolk to form and grow. The developing embryo is a closed system that can be readily biochemically analyzed, easily visualized, and manipulated to understand the role of nutrients in tissue specification, organogenesis, and growth. Natural variation in yolk composition observed across fish species may be related to unique developmental strategies. In this review, we discuss the reasons that teleost fishes are powerful models to understand nutritional control of development and highlight three species that are particularly valuable for future investigations: the zebrafish, Danio rerio, the African Killifish, Nothobranchius furzeri, and the Mexican tetra, Astyanax mexicanus. This review is a part of a special issue on nutritional, hormonal, and metabolic drivers of development.

Anatomy and evolution of bioluminescent organs in the slimeheads (Teleostei: Trachichthyidae)

Bacterial bioluminescent organs in fishes have a diverse range of tissues of origin, patterns of compartmentalization, and associated light-conducting structures. The morphology of the perianal, bacterial bioluminescent organ of Aulotrachichthys prosthemius was described previously, but the light organ in other species of slimeheads, family Trachichthyidae, is poorly known. Here, we describe the anatomy of the bioluminescent organs in trachichthyids and places the evolution of this light-producing system in the context of a new phylogeny of the Trachichthyoidei to test the hypothesis that bioluminescence evolved twice in the suborder and that the light-producing component derives from the perianal ectoderm. We use gross and histological examination to provide the first description of the bioluminescent organ of Paratrachichthys and four additional species of Aulotrachichthys. Observations also strongly suggest the presence of a perianal bioluminescent organ in Sorosichthys ananasa. The updated phylogeny of the Trachichthyoidei is the first to combine morphological and DNA-sequence (11-gene fragments) evidence, and supports a monophyletic Trachichthyidae with component subfamilies Hoplostethinae and Trachichthyinae, supporting continued recognition of the family Anoplogastridae. All bioluminescent trachichthyoids share a similar bioluminescent-organ structure with elongate chambers filled with bacteria and connected to collecting ducts that, in turn, connect to superficial ducts that lead to and have lining epithelia continuous with the epidermis. In the context of the phylogeny, the bioluminescent organ of trachichthyids is inferred to have evolved as an elaboration of the proctodeum in the ancestor of Aulotrachichthys, Paratrachichthys, and Sorosichthys independently from the structurally similar cephalic bioluminescent organs in Anomalopidae and Monocentridae.