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

Genomic signatures associated with recurrent scale loss in cyprinid fish

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

Chromosome-level genome assembly reveals adaptive evolution of the invasive Amazon sailfin catfish (Pterygoplichthys pardalis)

Catfish represents a diverse lineage with variable number of chromosomes and complex relationships with humans. Although certain species pose significant invasive threats to native fish populations, comprehensive genomic investigations into the evolutionary adaptations that contribute to their invasion success are lacking. To address this gap, our study presents a high-quality genome assembly of the Amazon sailfin catfish (Pterygoplichthys pardalis), a member of the armored catfish family, along with a comprehensive comparative genomic analysis. By utilizing conserved genomic regions across different catfish species, we reconstructed the 29 ancestral chromosomes of catfish, including two microchromosomes (28 and 29) that show different fusion and breakage patterns across species. Our analysis shows that the Amazon sailfin catfish genome is notably larger (1.58 Gb) with more than 40,000 coding genes. The genome expansion was linked to early repetitive sequence expansions and recent gene duplications. Several expanded genes are associated with immune functions, including antigen recognition domains like the Ig-v-set domain and the tandem expansion of the CD300 gene family. We also identified specific insertions in CNEs (conserved non-coding elements) near genes involved in cellular processes and neural development. Additionally, rapidly evolving and positively selected genes in the Amazon sailfin catfish genome were found to be associated with collagen formation. Moreover, we identified multiple positively selected codons in hoxb9, which may lead to functional alterations. These findings provide insights into molecular adaptations in an invasive catfish that may underlie its widespread invasion success.

Cranial kinematics and prey-type effects in Amia ocellicauda feeding strikes

Variability in the biomechanics and kinematics of prey capture in vertebrates has been studied extensively, with evidence of multiple strategies for successful feeding in many taxa. Early research into suction-feeding strikes in fishes hypothesized that fish utilize a set of pre-programmed strike kinematics that cannot be altered once initiated. However, more recent evidence has demonstrated that teleost fishes not only deploy unique strike kinematics for different prey types, but they also alter their kinematics in response to a prey item attempting to escape. It has not yet been explicitly investigated whether non-teleost actinopterygians can also modulate the strike in response to different prey types. Here, we examined the kinematics of suction strikes in bowfin, Amia ocellicauda, a holostean fish most closely related to gars. We recorded Amia feeding on both feeder fish and worms, two types of live prey differing in evasiveness, using X-ray reconstruction of moving morphology. We found significant prey-type effects on the magnitude, timing and velocity of jaw opening, hyoid arch depression and pectoral girdle motions. These prey-type effects demonstrate that the ability to modulate feeding strikes evolved early in actinopterygian fishes and is possibly the ancestral state for jawed vertebrates.

Genomics comparisons provide new insights into the evolution of karyotype and body patterns in Anguilliformes species

Anguilliformes species not only possess distinctive appearance (such as body elongation and absence of pectoral and/or pelvic fins), but also display diversity in chromosome number, supporting them as a suitable model for studying karyotype evolution and related molecular mechanisms of evolutionary body patterns. However, the ancestral chromosomes and evolutionary chromosomal reorganization in various eels have not been reported yet. The most regulatory or related genes of their distinctive appearance are still unknown. Here, we predicted an eel-based ancestral chromosome karyotype for the first time, and revealed multiple chromosomal fusion and fission events that reduced the ancestral chromosome number from haploid n = 21 to the commonly extant n = 19 within the Anguilla lineage. Moreover, we carried out a genome-wide comparison of two significant gene families including homeobox (Hox) and T-box (tbx), revealing genomic loss of some Hox genes (such as HoxB9β and HoxD13α) and variation of certain tbx gene (i.e., tbx5) may be responsible for the evolutionary development of pectoral fins. Interestingly, loss of certain secretory calcium-binding phosphoprotein (SCPP) genes was identified in various eel genomes, which possibly contribute to the common reduction of scales. Overall, our current findings provide new insights into evolutionary karyotype and body pattern evolution across diverse Anguilliformes species.

Understanding vertebrate immunity through comparative immunology

Evolutionary immunology has entered a new era. Classical studies, using just a handful of model animal species, combined with clinical observations, provided an outline of how innate and adaptive immunity work together to ensure tissue homeostasis and to coordinate the fight against infections. However, revolutionary advances in cellular and molecular biology, genomics and methods of genetic modification now offer unprecedented opportunities. They provide immunologists with the possibility to consider, at unprecedented scale, the impact of the astounding phenotypic diversity of vertebrates on immune system function. This Perspective is intended to highlight some of the many interesting, but largely unexplored, biological phenomena that are related to immune function among the roughly 60,000 existing vertebrate species. Importantly, hypotheses arising from such wide-ranging comparative studies can be tested in representative and genetically tractable species. The emerging general principles and the discovery of their evolutionarily selected variations may inspire the future development of novel therapeutic strategies for human immune disorders.

Investigating the Impact of Whole-Genome Duplication on Transposable Element Evolution in Teleost Fishes

Transposable elements (TEs) can make up more than 50% of any given vertebrate's genome, with substantial variability in TE composition among lineages. TE variation is often linked to changes in gene regulation, genome size, and speciation. However, the role that genome duplication events have played in generating abrupt shifts in the composition of the mobilome over macroevolutionary timescales remains unclear. We investigated the degree to which the teleost genome duplication (TGD) shaped the diversification trajectory of the teleost mobilome. We integrate a new high coverage genome of Polypterus bichir with data from over 100 publicly available actinopterygian genomes to assess the macroevolutionary implications of genome duplication events on TE evolution in teleosts. Our results provide no evidence for a substantial shift in mobilome composition following the TGD event. Instead, the diversity of the teleost mobilome appears to have been shaped by a history of lineage-specific shifts in composition that are not correlated with commonly evoked drivers of diversification such as body size, water column usage, or latitude. Collectively, these results provide additional evidence for an emerging perspective that TGD did not catalyze bursts of diversification and innovation in the actinopterygian mobilome.

Transcriptomic data support phylogenetic congruence and reveal genomic changes associated with the repeated evolution of annualism in aplocheiloid killifishes (Cyprinodontiformes)

Repeated evolution of novel life histories that are correlated with ecological variables offers opportunities to study convergence in genetic, developmental, and metabolic features. Nearly half of the 800 species of Aplocheiloid killifishes, a clade of teleost fishes with a circumtropical distribution, are "annual" or seasonal species that survive in ephemeral bodies of water that desiccate and are unfeasible for growth, reproduction, or survival for weeks to months every year. But the repeated evolution of adaptations that are key features of the annual life history among these fishes remains poorly known without a robust phylogenetic framework. We present a large-scale phylogenomic reconstruction of aplocheiloid killifishes evolution using newly sequenced transcriptomes obtained from a diversity of killifish lineages representing putative independent origins of annualism. Ancestral state estimation shows that developmental dormancy (diapause), a key trait of the killifish annual life cycle, may have originated up to seven times independently among African and South American lineages. To further explore the genetic basis of this unique trait, we measure changes in evolutionary rates among orthologous genes across the killifish tree of life by quantifying codon evolution using dN/dS ratios. We show that some genes have higher dN/dS ratios in lineages leading to species with annual life history. Many of them constitute key developmental genes or nuclear-encoded metabolic genes that control oxidative phosphorylation. Lastly, we compare these genes with higher ω to genes previously associated to developmental dormancy and metabolic shifts in killifishes and other vertebrates, and thereby identify molecular evolutionary signatures of repeated transitions to extreme environments.

Analysis of a shark reveals ancient, Wnt-dependent, habenular asymmetries in vertebrates

The mode of evolution of left-right asymmetries in the vertebrate habenulae remains largely unknown. Using a transcriptomic approach, we show that in a cartilaginous fish, the catshark Scyliorhinus canicula, habenulae exhibit marked asymmetries, in both their medial and lateral components. Comparisons across vertebrates suggest that those identified in lateral habenulae reflect an ancestral gnathostome trait, partially conserved in lampreys, and independently lost in tetrapods and neopterygians. Asymmetry formation involves distinct mechanisms in the catshark lateral and medial habenulae. Medial habenulae are submitted to a marked, asymmetric temporal regulation of neurogenesis, undetectable in their lateral counterparts. Conversely, asymmetry formation in lateral habenulae results from asymmetric choices of neuronal identity in post-mitotic progenitors, a regulation dependent on the repression of Wnt signaling by Nodal on the left. Based on comparisons with the mouse and the zebrafish, we propose that habenular asymmetry formation involves a recurrent developmental logic across vertebrates, which relies on conserved, temporally regulated genetic programs sequentially shaping choices of neuronal identity on both sides and asymmetrically modified by Wnt activity.

Anterior-posterior constraint on Hedgehog signaling by hhip in teleost fin elaboration

Pectoral fins, the anterior paired fins in fish, have enhanced maneuvering abilities due to morphological changes. Teleosts have fewer radial bones in their pectoral fins than basal species, resulting in more-elaborate fins. The mechanism behind this radial constraint change in teleosts is unclear. Here, we found that mutations in hhip, which encodes an antagonist of Hedgehog signaling, led to an increase in radial bones in a localized region. Expression of the Shh genes, encoding ligands of Hedgehog signaling, coincided with notable hhip expression specifically during early development. We suggest that a negative feedback effect of Hedgehog signaling by hhip regulates the constraint of the pectoral fin in zebrafish. Additionally, re-analysis of hhip-related gene expression data in zebrafish and basal species revealed that the notable hhip expression during early development is a characteristic of zebrafish that is not observed in basal species. Region-specific expression of Hox13 genes in the zebrafish pectoral fin indicated that the median region, analogous to the region with abundant radials in basal species, is expanded in hhip-/- zebrafish. These data underscore potential morphological evolution through constrained diversity.

The Chromosome-Scale Genome of Chitala ornata Illuminates the Evolution of Early Teleosts

Teleosts are the most prolific vertebrates, occupying the vast majority of aquatic environments, and their pectoral fins have undergone remarkable physiological transformations throughout their evolution. Studying early teleost fishes, such as those belonging to the Osteoglossiformes order, could offer crucial insights into the adaptive evolution of pectoral fins within this group. In this study, we have assembled a chromosomal-level genome for the Clown featherback (Chitala ornata), achieving the highest quality genome assembly for Osteoglossiformes to date, with a contig N50 of 32.78 Mb and a scaffold N50 of 40.73 Mb. By combining phylogenetic analysis, we determined that the Clown featherback diverged approximately 202 to 203 million years ago (Ma), aligning with continental separation events. Our analysis revealed the intriguing discovery that a unique deletion of regulatory elements is adjacent to the Gli3 gene, specifically in teleosts. This deletion might be tied to the specialized adaptation of their pectoral fins. Furthermore, our findings indicate that specific contractions and expansions of transposable elements (TEs) in teleosts, including the Clown featherback, could be connected to their adaptive evolution. In essence, this study not only provides a high-quality genomic resource for Osteoglossiformes but also sheds light on the evolutionary trajectory of early teleosts.

Paired fins in vertebrate evolution and ontogeny

The origin of paired appendages became one of the most important adaptations of vertebrates, allowing them to lead active lifestyles and explore a wide range of ecological niches. The basic form of paired appendages in evolution is the fins of fishes. The problem of paired appendages has attracted the attention of researchers for more than 150 years. During this time, a number of theories have been proposed, mainly based on morphological data, two of which, the Balfour-Thacher-Mivart lateral fold theory and Gegenbaur's gill arch theory, have not lost their relevance. So far, however, none of the proposed ideas has been supported by decisive evidence. The study of the evolutionary history of the appearance and development of paired appendages lies at the intersection of several disciplines and involves the synthesis of paleontological, morphological, embryological, and genetic data. In this review, we attempt to summarize and discuss the results accumulated in these fields and to analyze the theories put forward regarding the prerequisites and mechanisms that gave rise to paired fins and limbs in vertebrates.

Genomics comparisons of three chromosome-level mudskipper genome assemblies reveal molecular clues for water-to-land evolution and adaptation

INTRODUCTION

Mudskippers are a large group of amphibious fishes that have developed many morphological and physiological capacities to live on land. Genomics comparisons of chromosome-level genome assemblies of three representative mudskippers, Boleophthalmus pectinirostris (BP), Periophthalmus magnuspinnatus (PM) and P. modestus (PMO), may be able to provide novel insights into the water-to-land evolution and adaptation.

METHODS

Two chromosome-level genome assemblies for BP and PM were respectively sequenced by an integration of PacBio, Nanopore and Hi-C sequencing. A series of standard assembly and annotation pipelines were subsequently performed for both mudskippers. We also re-annotated the PMO genome, downloaded from NCBI, to obtain a redundancy-reduced annotation. Three-way comparative analyses of the three mudskipper genomes in a large scale were carried out to discover detailed genomic differences, such as different gene sizes, and potential chromosomal fission and fusion events. Comparisons of several representative gene families among the three amphibious mudskippers and some other teleosts were also performed to find some molecular clues for terrestrial adaptation.

RESULTS

We obtained two high-quality haplotype genome assemblies with 23 and 25 chromosomes for BP and PM respectively. We also found two specific chromosome fission events in PM. Ancestor chromosome analysis has discovered a common fusion event in mudskipper ancestor. This fusion was then retained in all the three mudskipper species. A loss of some SCPP (secretory calcium-binding phosphoprotein) genes were identified in the three mudskipper genomes, which could lead to reduction of scales for a part-time terrestrial residence. The loss of aanat1a gene, encoding an important enzyme (arylalkylamine N-acetyltransferase 1a, AANAT1a) for dopamine metabolism and melatonin biosynthesis, was confirmed in PM but not in PMO (as previously reported existence in BP), suggesting a better air vision of PM than both PMO and BP. Such a tiny variation within the genus Periophthalmus exemplifies to prove a step-by-step evolution for the mudskippers' water-to-land adaptation.

CONCLUSION

These high-quality mudskipper genome assemblies will become valuable genetic resources for in-depth discovery of genomic evolution for the terrestrial adaptation of amphibious fishes.

Convergent genomic signatures associated with vertebrate viviparity

BACKGROUND

Viviparity-live birth-is a complex and innovative mode of reproduction that has evolved repeatedly across the vertebrate Tree of Life. Viviparous species exhibit remarkable levels of reproductive diversity, both in the amount of care provided by the parent during gestation, and the ways in which that care is delivered. The genetic basis of viviparity has garnered increasing interest over recent years; however, such studies are often undertaken on small evolutionary timelines, and thus are not able to address changes occurring on a broader scale. Using whole genome data, we investigated the molecular basis of this innovation across the diversity of vertebrates to answer a long held question in evolutionary biology: is the evolution of convergent traits driven by convergent genomic changes?

RESULTS

We reveal convergent changes in protein family sizes, protein-coding regions, introns, and untranslated regions (UTRs) in a number of distantly related viviparous lineages. Specifically, we identify 15 protein families showing evidence of contraction or expansion associated with viviparity. We additionally identify elevated substitution rates in both coding and noncoding sequences in several viviparous lineages. However, we did not find any convergent changes-be it at the nucleotide or protein level-common to all viviparous lineages.

CONCLUSIONS

Our results highlight the value of macroevolutionary comparative genomics in determining the genomic basis of complex evolutionary transitions. While we identify a number of convergent genomic changes that may be associated with the evolution of viviparity in vertebrates, there does not appear to be a convergent molecular signature shared by all viviparous vertebrates. Ultimately, our findings indicate that a complex trait such as viviparity likely evolves with changes occurring in multiple different pathways.

An Ancestral Major Histocompatibility Complex Organization in Cartilaginous Fish: Reconstructing MHC Origin and Evolution

Cartilaginous fish (sharks, rays, and chimeras) comprise the oldest living jawed vertebrates with a mammalian-like adaptive immune system based on immunoglobulins (Ig), T-cell receptors (TCRs), and the major histocompatibility complex (MHC). Here, we show that the cartilaginous fish "adaptive MHC" is highly regimented and compact, containing (i) a classical MHC class Ia (MHC-Ia) region containing antigen processing (antigen peptide transporters and immunoproteasome) and presenting (MHC-Ia) genes, (ii) an MHC class II (MHC-II) region (with alpha and beta genes) with linkage to beta-2-microglobulin (β2m) and bromodomain-containing 2, (iii) nonclassical MHC class Ib (MHC-Ib) regions with 450 million-year-old lineages, and (iv) a complement C4 associated with the MHC-Ia region. No MHC-Ib genes were found outside of the elasmobranch MHC. Our data suggest that both MHC-I and MHC-II genes arose after the second round of whole-genome duplication (2R) on a human chromosome (huchr) 6 precursor. Further analysis of MHC paralogous regions across early branching taxa from all jawed vertebrate lineages revealed that Ig/TCR genes likely arose on a precursor of the huchr9/12/14 MHC paralog. The β2m gene is linked to the Ig/TCR genes in some vertebrates suggesting that it was present at 1R, perhaps as the donor of C1 domain to the primordial MHC gene. In sum, extant cartilaginous fish exhibit a conserved and prototypical MHC genomic organization with features found in various vertebrates, reflecting the ancestral arrangement for the jawed vertebrates.

Fish hearing revealed: Do we understand hearing in critical fishes and marine tetrapods

Hearing evolved in lampreys with a frequency range of 50-200 Hz. This hearing range is comparable to that of elasmobranchs, most non-teleosts, and lungfish. Elasmobranchs most likely use the saccule and the papilla neglecta (PN) for hearing. In non-teleosts and teleosts, lungfish, and certain tetrapods the saccule is the likely sensor for sound reception while the lagena and the PN are important for gravistatic sensing. Coelacanth and most tetrapods have a basilar papilla (BP) for hearing. In coelacanth and tetrapods, the hair cells of the BP are in contact with a basilar and a tectorial membrane. These membranes transmit mechanical vibrations. A cochlear aqueduct (CA) provides a connection between the cerebrospinal fluid that has a sodium rich space in coelacanth and tetrapods while the potassium rich endolymph is known in vertebrates. A unique feature is known in basic sarcopterygians, the intracranial joint, that never developed in actinopterygians and has been lost in lungfish and tetrapods. The BP in coelacanths is thought to generate pressure with the intracranial joint that will be transmitted to the CA. Lungs or a swim bladder are not forming in Chondrichthyes, structures that have a major impact on hearing in teleosts and tetrapods.

Conserved enhancers control notochord expression of vertebrate Brachyury

The cell type-specific expression of key transcription factors is central to development and disease. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three conserved Brachyury-controlling notochord enhancers, T3, C, and I, in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, in cis deletion of all three enhancers in mouse abolishes Brachyury/T/Tbxt expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. The three Brachyury-driving notochord enhancers are conserved beyond mammals in the brachyury/tbxtb loci of fishes, dating their origin to the last common ancestor of jawed vertebrates. Our data define the vertebrate enhancers for Brachyury/T/TBXTB notochord expression through an auto-regulatory mechanism that conveys robustness and adaptability as ancient basis for axis development.

Ancient fish lineages illuminate toll-like receptor diversification in early vertebrate evolution

Since its initial discovery over 50 years ago, understanding the evolution of the vertebrate RAG- mediated adaptive immune response has been a major area of research focus for comparative geneticists. However, how the evolutionary novelty of an adaptive immune response impacted the diversity of receptors associated with the innate immune response has received considerably less attention until recently. Here, we investigate the diversification of vertebrate toll-like receptors (TLRs), one of the most ancient and well conserved innate immune receptor families found across the Tree of Life, integrating genomic data that represent all major vertebrate lineages with new transcriptomic data from Polypteriformes, the earliest diverging ray-finned fish lineage. Our analyses reveal TLR sequences that reflect the 6 major TLR subfamilies, TLR1, TLR3, TLR4, TLR5, TLR7, and TLR11, and also currently unnamed, yet phylogenetically distinct TLR clades. We additionally recover evidence for a pulse of gene gain coincident with the rise of the RAG-mediated adaptive immune response in jawed vertebrates, followed by a period of rapid gene loss during the Cretaceous. These gene losses are primarily concentrated in marine teleost fish and synchronous with the mid Cretaceous anoxic event, a period of rapid extinction for marine species. Finally, we reveal a mismatch between phylogenetic placement and gene nomenclature for up to 50% of TLRs found in clades such as ray-finned fishes, cyclostomes, amphibians, and elasmobranchs. Collectively, these results provide an unparalleled perspective of TLR diversity and offer a ready framework for testing gene annotations in non-model species.

A single-cell atlas of West African lungfish respiratory system reveals evolutionary adaptations to terrestrialization

The six species of lungfish possess both lungs and gills and are the closest extant relatives of tetrapods. Here, we report a single-cell transcriptome atlas of the West African lungfish (Protopterus annectens). This species manifests the most extreme form of terrestrialization, a life history strategy to survive dry periods that can last for years, characterized by dormancy and reversible adaptive changes of the gills and lungs. Our atlas highlights the cell type diversity of the West African lungfish, including gene expression consistent with phenotype changes of terrestrialization. Comparison with terrestrial tetrapods and ray-finned fishes reveals broad homology between the swim bladder and lung cell types as well as shared and idiosyncratic changes of the external gills of the West African lungfish and the internal gills of Atlantic salmon. The single-cell atlas presented here provides a valuable resource for further exploration of the respiratory system evolution in vertebrates and the diversity of lungfish terrestrialization.

Abandoning the Isochore Theory Can Help Explain Genome Compositional Organization in Fish

The organization of the genome nucleotide (AT/GC) composition in vertebrates remains poorly understood despite the numerous genome assemblies available. Particularly, the origin of the AT/GC heterogeneity in amniotes, in comparison to the homogeneity in anamniotes, is controversial. Recently, several exceptions to this dichotomy were confirmed in an ancient fish lineage with mammalian AT/GC heterogeneity. Hence, our current knowledge necessitates a reevaluation considering this fact and utilizing newly available data and tools. We analyzed fish genomes in silico with as low user input as possible to compare previous approaches to assessing genome composition. Our results revealed a disparity between previously used plots of GC% and histograms representing the authentic distribution of GC% values in genomes. Previous plots heavily reduced the range of GC% values in fish to comply with the alleged AT/GC homogeneity and AT-richness of their genomes. We illustrate how the selected sequence size influences the clustering of GC% values. Previous approaches that disregarded chromosome and genome sizes, which are about three times smaller in fish than in mammals, distorted their results and contributed to the persisting confusion about fish genome composition. Chromosome size and their transposons may drive the AT/GC heterogeneity apparent on mammalian chromosomes, whereas far less in fishes.

A chromosome-level genome assembly of longnose gar, Lepisosteus osseus

Holosteans (gars and bowfins) represent the sister lineage to teleost fishes, the latter being a clade that comprises over half of all living vertebrates and includes important models for comparative genomics and human health. A major distinction between the evolutionary history of teleosts and holosteans is that all teleosts experienced a genome duplication event in their early evolutionary history. As the teleost genome duplication occurred after teleosts diverged from holosteans, holosteans have been heralded as a means to bridge teleost models to other vertebrate genomes. However, only three species of holosteans have been genome-sequenced to date, and sequencing of more species is needed to fill sequence sampling gaps and provide a broader comparative basis for understanding holostean genome evolution. Here we report the first high quality reference genome assembly and annotation of the longnose gar (Lepisosteus osseus). Our final assembly consists of 22,709 scaffolds with a total length of 945 bp with contig N50 of 116.61 kb. Using BRAKER2, we annotated a total of 30,068 genes. Analysis of the repetitive regions of the genome reveals the genome to contain 29.12% transposable elements, and the longnose gar to be the only other known vertebrate outside of the spotted gar and bowfin to contain CR1, L2, Rex1, and Babar. These results highlight the potential utility of holostean genomes for understanding the evolution of vertebrate repetitive elements, and provide a critical reference for comparative genomic studies utilizing ray-finned fish models.

Independent rediploidization masks shared whole genome duplication in the sturgeon-paddlefish ancestor

Whole genome duplication (WGD) is a dramatic evolutionary event generating many new genes and which may play a role in survival through mass extinctions. Paddlefish and sturgeon are sister lineages that both show genomic evidence for ancient WGD. Until now this has been interpreted as two independent WGD events due to a preponderance of duplicate genes with independent histories. Here we show that although there is indeed a plurality of apparently independent gene duplications, these derive from a shared genome duplication event occurring well over 200 million years ago, likely close to the Permian-Triassic mass extinction period. This was followed by a prolonged process of reversion to stable diploid inheritance (rediploidization), that may have promoted survival during the Triassic-Jurassic mass extinction. We show that the sharing of this WGD is masked by the fact that paddlefish and sturgeon lineage divergence occurred before rediploidization had proceeded even half-way. Thus, for most genes the resolution to diploidy was lineage-specific. Because genes are only truly duplicated once diploid inheritance is established, the paddlefish and sturgeon genomes are thus a mosaic of shared and non-shared gene duplications resulting from a shared genome duplication event.

The little skate genome and the evolutionary emergence of wing-like fins

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments1,2. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins-including gene expression, chromatin occupancy and three-dimensional conformation-we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

A highly diverse set of novel immunoglobulin-like transcript (NILT) genes in zebrafish indicates a wide range of functions with complex relationships to mammalian receptors

Multiple novel immunoglobulin-like transcripts (NILTs) have been identified from salmon, trout, and carp. NILTs typically encode activating or inhibitory transmembrane receptors with extracellular immunoglobulin (Ig) domains. Although predicted to provide immune recognition in ray-finned fish, we currently lack a definitive framework of NILT diversity, thereby limiting our predictions for their evolutionary origin and function. In order to better understand the diversity of NILTs and their possible roles in immune function, we identified five NILT loci in the Atlantic salmon (Salmo salar) genome, defined 86 NILT Ig domains within a 3-Mbp region of zebrafish (Danio rerio) chromosome 1, and described 41 NILT Ig domains as part of an alternative haplotype for this same genomic region. We then identified transcripts encoded by 43 different NILT genes which reflect an unprecedented diversity of Ig domain sequences and combinations for a family of non-recombining receptors within a single species. Zebrafish NILTs include a sole putative activating receptor but extensive inhibitory and secreted forms as well as membrane-bound forms with no known signaling motifs. These results reveal a higher level of genetic complexity, interindividual variation, and sequence diversity for NILTs than previously described, suggesting that this gene family likely plays multiple roles in host immunity.

Structure/Function Studies on the Activation Motif of Two Non-Mammalian Mrap1 Orthologs, and Observations on the Phylogeny of Mrap1, Including a Novel Characterization of an Mrap1 from the Chondrostean Fish, Polyodon spathula

In derived bony vertebrates, activation of the melanocortin-2 receptor (Mc2r) by its ACTH ligand requires chaperoning by the Mc2r accessory protein (Mrap1). The N-terminal domain of the non-mammalian tetrapod MRAP1 from chicken (c; Gallus gallus) has the putative activation motif, W18D19Y20I21, and the N-terminal domain in the neopterygian ray-finned fish Mrap1 from bowfin (bf; Amia calva) has the putative activation motif, Y18D19Y20I21. The current study used an alanine-substitution paradigm to test the hypothesis that only the Y20 position in the Mrap1 ortholog of these non-mammalian vertebrates is required for activation of the respective Mc2r ortholog. Instead, we found that for cMRAP1, single alanine-substitution resulted in a gradient of inhibition of activation (Y20 >> D19 = W18 > I21). For bfMrap1, single alanine-substitution also resulted in a gradient of inhibition of activation (Y20 >> D19 > I21 > Y18). This study also included an analysis of Mc2r activation in an older lineage of ray-finned fish, the paddlefish (p), Polyodon spathula (subclass Chondronstei). Currently no mrap1 gene has been found in the paddlefish genome. When pmc2r was expressed alone in our CHO cell/cAMP reporter gene assay, no activation was observed following stimulation with ACTH. However, when pmc2r was co-expressed with either cmrap1 or bfmrap1 robust dose response curves were generated. These results indicate that the formation of an Mc2r/Mrap1 heterodimer emerged early in the radiation of the bony vertebrates.

A basal actinopterygian melanocortin receptor: Molecular and functional characterization of an Mc2r ortholog from the Senegal bichir (Polypterus senegalus)

In bony vertebrates, melanocortin 2 receptor (Mc2r) specifically binds adrenocorticotropic hormone (ACTH) and is responsible for mediating anterior pituitary signaling that stimulates corticosteroid production in the adrenal gland/interrenal cells. In bony fishes Mc2r requires the chaperoning of an accessory protein (Mrap1) to traffic to the membrane surface and bind ACTH. Here, we evaluated the structure and pharmacological properties of Mc2r from the Senegal bichir (Polypterus senegalus), which represents the most basal bony fish from which an Mc2r has been pharmacologically studied to date. In our experiments, cDNA constructs of the Mc2r from the Senegal bichir (sbMc2r) and various vertebrate Mrap1s were heterologously co-expressed in Chinese hamster ovary (CHO) cells, stimulated by ACTH or melanocyte-stimulating hormone (α-MSH) ligands, and assessed using a luciferase reporter gene assay. When expressed without an Mrap1, sbMc2r was not activated by ACTH. When co-expressed with Mrap1 from either chicken (Gallus gallus) or bowfin (Amia calva), sbMc2r could be activated in a dose-dependent manner by ACTH, but not α-MSH. Co-expression of sbMrap2 with sbMc2r resulted in no detectable activation of the receptor. Collectively, these results demonstrate that sbMc2r has pharmacological properties similar to those of Mc2rs of later-evolved bony fishes, such as Mrap1 dependence and ACTH selectivity, indicating that these qualities of Mc2r function are ancestral to all bony fish Mc2rs.

Otoliths suggest lifespans more than 30 years for free-living bowfin Amia calva: Implications for fisheries management in the bowfishing era

The bowfin Amia calva is an amiid (Amiiformes) relict native to North America. It is the last surviving member of the Halecomorphi, a group of fishes that evolved more than 250 million years ago. Despite the phylogenetic significance of the amiids in vertebrate evolution, little has been published about their age and growth. Recreational bowfin harvest is currently unregulated throughout most of the USA, yet new recreational fisheries are emerging. As such, bowfin are increasingly harvested by sport bowfishing without limit, in addition to their growing commercial harvest for caviar. From 2017 to 2021 we studied a total of 81 bowfin from 11 populations across the east-west gradient of Minnesota within a narrow latitudinal margin (<50 km) of the 46th parallel north. We compared the allometry and translucence of bowfin asteriscus, lapillus and sagittal otoliths and found the lapillus otoliths provide consistent readability for age estimation despite being the smallest of the set. Size-at-age data derived from otoliths indicated that bowfin are sexually dimorphic in asymptotic length and may live up to 33 years, which is 15 years longer than previously estimated in wild populations, but comparable to what has been reported in captivity. Overall, 28% of the otolith-aged fish were estimated as older than the previously reported maximum age for wild bowfin populations. Our findings suggest that the bowfin life history may exhibit slower growth, greater longevity, and more variable recruitment than previously recognized, which sets the stage for more otolith-derived population demographics across their range and age validation study. Our results have direct implications for conservation of bowfin, especially amidst the increasing rates of exploitation during the bowfishing era.

Hidden species diversity in an iconic living fossil vertebrate

Ancient, species-poor lineages persistently occur across the Tree of life. These lineages are likely to contain unrecognized species diversity masked by the low rates of morphological evolution that characterize living fossils. Halecomorphi is a lineage of ray-finned fishes that diverged from its closest relatives before 200 Ma and is represented by only one living species in eastern North America, the bowfin, Amia calva Linnaeus. Here, we use double digest restriction-site-associated DNA sequencing and morphology to illuminate recent speciation in bowfins. Our results support the delimitation of a second living species of Amia, with the timing of diversification dating to the Plio-Pleistocene. This delimitation expands the species diversity of an ancient lineage that is integral to studies of vertebrate genomics and development, yet is facing growing conservation threats driven by the caviar fishery.

Phylogenomic analysis of the bowfin (Amia calva) reveals unrecognized species diversity in a living fossil lineage

The Bowfin (Amia calva), as currently recognized, represents the sole living member of the family Amiidae, which dates back to approximately 150 Ma. Prior to 1896, 13 species of extant Bowfins had been described, but these were all placed into a single species with no rationale or analysis given. This situation has persisted until the present day, with little attention given to re-evaluation of those previously described nominal forms. Here, we present a phylogenomic analysis based on over 21,000 single nucleotide polymorphisms (SNPs) from 94 individuals that unambiguously demonstrates the presence of at least two independent evolutionary lineages within extant Amia populations that merit species-level standing, as well as the possibility of two more. These findings not only expand the recognizable species diversity in an iconic, ancient lineage, but also demonstrate the utility of such methods in addressing previously intractable questions of molecular systematics and phylogeography in slowly evolving groups of ancient fishes.

Whole genome assembly of the armored loricariid catfish Ancistrus triradiatus highlights herbivory signatures

The catfish Ancistrus triradiatus belongs to the species-rich family Loricariidae. Loricariids display remarkable traits such as herbivory, a benthic lifestyle, the absence of scales but the presence of dermal bony plates. They are exported as ornamental fish worldwide, with escaped fishes becoming a threat locally. Although genetic and phylogenetic studies are continuously increasing and developmental genetic investigations are underway, no genome assembly has been formally proposed for Loricariidae yet. We report a high-quality genome assembly of Ancistrus triradiatus using long and short reads, and a newly assembled transcriptome. The genome assembly is composed of 9530 scaffolds, including 85.6% of ray-finned fish BUSCOs, and 26,885 predicted protein-coding genes. The genomic GC content is higher than in other catfishes, reflecting the higher metabolism associated with herbivory. The examination of the SCPP gene family indicates that the genes presumably triggering scale loss when absent, are present in the scaleless A. triradiatus, questioning their explanatory role. The analysis of the opsin gene repertoire revealed that gene losses associated to the nocturnal lifestyle of catfishes were not entirely found in A. triradiatus, as the UV-sensitive opsin 5 is present. Finally, most gene family expansions were related to immunity except the gamma crystallin gene family which controls pupil shape and sub-aquatic vision. Thus, the genome of A. triradiatus reveals that fish herbivory may be related to the photic zone habitat, conditions metabolism, photoreception and visual functions. This genome is the first for the catfish suborder Loricarioidei and will serve as backbone for future genetic, developmental and conservation studies.

Evolution of the nitric oxide synthase family in vertebrates and novel insights in gill development

Nitric oxide (NO) is an ancestral key signalling molecule essential for life and has enormous versatility in biological systems, including cardiovascular homeostasis, neurotransmission and immunity. Although our knowledge of NO synthases (Nos), the enzymes that synthesize NO in vivo, is substantial, the origin of a large and diversified repertoire of nos gene orthologues in fishes with respect to tetrapods remains a puzzle. The recent identification of nos3 in the ray-finned fish spotted gar, which was considered lost in this lineage, changed this perspective. This finding prompted us to explore nos gene evolution, surveying vertebrate species representing key evolutionary nodes. This study provides noteworthy findings: first, nos2 experienced several lineage-specific gene duplications and losses. Second, nos3 was found to be lost independently in two different teleost lineages, Elopomorpha and Clupeocephala. Third, the expression of at least one nos paralogue in the gills of developing shark, bichir, sturgeon, and gar, but not in lamprey, suggests that nos expression in this organ may have arisen in the last common ancestor of gnathostomes. These results provide a framework for continuing research on nos genes' roles, highlighting subfunctionalization and reciprocal loss of function that occurred in different lineages during vertebrate genome duplications.

Whole-genome survey and phylogenetic analysis of Gadus macrocephalus

Gadus macrocephalus (Pacific cod) is an economically important species on the northern coast of the Pacific. Although numerous studies on G. macrocephalus exist, there are few reports on its genomic data. Here, we used whole-genome sequencing data to elucidate the genomic characteristics and phylogenetic relationship of G. macrocephalus. From the 19-mer frequency distribution, the genome size was estimated to be 658.22 Mb. The heterozygosity, repetitive sequence content and GC content were approximately 0.62%, 27.50% and 44.73%, respectively. The draft genome sequences were initially assembled, yielding a total of 500,760 scaffolds (N50 = 3565 bp). A total of 789,860 microsatellite motifs were identified from the genomic data, and dinucleotide repeat was the most dominant simple sequence repeat motif. As a byproduct of whole-genome sequencing, the mitochondrial genome was assembled to investigate the evolutionary relationships between G. macrocephalus and its relatives. On the basis of 13 protein-coding gene sequences of the mitochondrial genome of Gadidae species, the maximum likelihood phylogenetic tree showed that complicated relationships and divergence times among Gadidae species. Demographic history analysis revealed changes in the G. macrocephalus population during the Pleistocene by using the pairwise sequentially Markovian coalescent model. These findings supplement the genomic data of G. macrocephalus, and make a valuable contribution to the whole-genome studies on G. macrocephalus.

Insights into chromosomal evolution and sex determination of Pseudobagrus ussuriensis (Bagridae, Siluriformes) based on a chromosome-level genome

Pseudobagrus ussuriensis is an aquaculture catfish with significant sexual dimorphism. In this study, a chromosome-level genome with a size of 741.97 Mb was assembled for female P. ussuriensis. A total of 26 chromosome-level contigs covering 97.34% of the whole-genome assembly were obtained with an N50 of 28.53 Mb and an L50 of 11. A total of 24,075 protein-coding genes were identified, with 91.54% (22,039) genes being functionally annotated. Based on the genome assembly, four chromosome evolution clusters of catfishes were identified and the formation process of P. ussuriensis chromosomes was predicted. A total of 55 sex-related quantitative trait loci (QTLs) with a phenotypic variance explained value of 100% were located on chromosome 8 (chr08). The QTLs and other previously identified sex-specific markers were located in a sex-determining region of 16.83 Mb (from 6.90 to 23.73 Mb) on chr08, which was predicted as the X chromosome. The sex-determining region comprised 554 genes, with 135 of which being differently expressed between males and females/pseudofemales, and 16 candidate sex-determining genes were screened out. The results of this study provided a useful chromosome-level genome for genetic, genomic and evolutionary studies of P. ussuriensis, and also be useful for further studies on sex-determination mechanism analysis and sex-control breeding of this fish.

Genome of the Rio Pearlfish (Nematolebias whitei), a bi-annual killifish model for Eco-Evo-Devo in extreme environments

The Rio Pearlfish, Nematolebias whitei, is a bi-annual killifish species inhabiting seasonal pools in the Rio de Janeiro region of Brazil that dry twice per year. Embryos enter dormant diapause stages in the soil, waiting for the inundation of the habitat which triggers hatching and commencement of a new life cycle. Rio Pearlfish represents a convergent, independent origin of annualism from other emerging killifish model species. While some transcriptomic datasets are available for Rio Pearlfish, thus far, a sequenced genome has been unavailable. Here, we present a high quality, 1.2 Gb chromosome-level genome assembly, genome annotations, and a comparative genomic investigation of the Rio Pearlfish as representative of a vertebrate clade that evolved environmentally cued hatching. We show conservation of 3D genome structure across teleost fish evolution, developmental stages, tissues, and cell types. Our analysis of mobile DNA shows that Rio Pearlfish, like other annual killifishes, possesses an expanded transposable element profile with implications for rapid aging and adaptation to harsh conditions. We use the Rio Pearlfish genome to identify its hatching enzyme gene repertoire and the location of the hatching gland, a key first step in understanding the developmental genetic control of hatching. The Rio Pearlfish genome expands the comparative genomic toolkit available to study convergent origins of seasonal life histories, diapause, and rapid aging phenotypes. We present the first set of genomic resources for this emerging model organism, critical for future functional genetic, and multiomic explorations of “Eco-Evo-Devo” phenotypes of resilience and adaptation to extreme environments.

Sonic hedgehog is Essential for Proximal-Distal Outgrowth of the Limb Bud in Salamanders

The developing forelimb has been a foundational model to understand how specified progenitor cells integrate genetic information to produce the tetrapod limb bauplan. Although the reigning hypothesis is that all tetrapods develop limbs in a similar manner, recent work suggests that urodeles have evolved a derived mode of limb dvelopment. Here, we demonstrate through pharmacological and genetic inactivation of Sonic hedgehog (Shh) signaling in axolotls that Shh directs expansion and survival of limb progenitor cells in addition to patterning the limb across the proximodistal and antero-posterior axis. In contrast to inactivation of Shh in mouse or chick embryos where a humerus, radius, and single digit develop, Shh crispant axolotls completely lack forelimbs. In rescuing limb development by implanting SHH-N protein beads into the nascent limb field of Shh crispants, we show that the limb field is specified in the absence of Shh and that hedgehog pathway activation is required to initiate proximodistal outgrowth. When our results are examined alongside other derived aspects of salamander limb development and placed in a phylogenetic context, a new hypothesis emerges whereby the ability for cells at an amputation plane to activate morphogenesis and regenerate a limb may have evolved uniquely in urodeles.

Diversification Rate is Associated with Rate of Molecular Evolution in Ray-Finned Fish (Actinopterygii)

Understanding the factors that drive diversification of taxa across the tree of life is a key focus of macroevolutionary research. While the effects of life history, ecology, climate and geography on diversity have been studied for many taxa, the relationship between molecular evolution and diversification has received less attention. However, correlations between rates of molecular evolution and diversification rate have been detected in a range of taxa, including reptiles, plants and birds. A correlation between rates of molecular evolution and diversification rate is a prediction of several evolutionary theories, including the evolutionary speed hypothesis which links variation in mutation rates to differences in speciation rates. If it is widespread, such correlations could also have significant practical impacts, if they are not adequately accounted for in phylogenetic inference of evolutionary rates and timescales. Ray-finned fish (Actinopterygii) offer a prime target to test for this relationship due to their extreme variation in clade size suggesting a wide range of diversification rates. We employ both a sister-pairs approach and a whole-tree approach to test for correlations between substitution rate and net diversification. We also collect life history and ecological trait data and account for potential confounding factors including body size, latitude, max depth and reef association. We find evidence to support a relationship between diversification and synonymous rates of nuclear evolution across two published backbone phylogenies, as well as weak evidence for a relationship between mitochondrial nonsynonymous rates and diversification at the genus level.

Convergent losses of SCPP genes and ganoid scales among non-teleost actinopterygians

Various secretory calcium-binding phosphoprotein (SCPP) genes are expressed in the skin and jaw during the formation of bone, teeth, and scales in osteichthyans (bony vertebrates). Among these mineralized skeletal units is the ganoid scale, found in many fossil actinopterygians (ray-finned fish) but confirmed only in Polypteriformes (bichirs, reedfish) and Lepisosteiformes (gars) among extant clades. Here, we examined SCPP genes in the genome of seven non-teleost actinopterygian species that possess or do not possess ganoid scales. As a result, 39-43 SCPP genes were identified in Polypteriformes and Lepisosteiformes, whereas 22-24 SCPP genes were found in Acipenseriformes (sturgeons, paddlefish) and Amiiformes (bowfin). Most of these genes form two clusters in the genome of Polypteriformes, Lepisosteiformes, and Amiiformes, and these two clusters are duplicated in Acipenseriformes. Despite their distant phylogenetic relationship, Polypteriformes and Lepisosteiformes retain many orthologous SCPP genes. These results imply that common ancestors of extant actinopterygians possessed a large repertoire of SCPP genes, and that many SCPP genes were lost independently in Acipenseriformes and Amiiformes. Notably, most SCPP genes originally located in one of the two SCPP gene clusters are retained in Polypteriformes and Lepisosteiformes but were secondarily lost in Acipenseriformes and Amiiformes. In Lepisosteiformes, orthologs of these lost genes show high or detectable expression levels in the skin but not in the jaw. We thus hypothesize that many SCPP genes located in this cluster are involved in the formation of ganoid scales in Polypteriformes and Lepisosteiformes, and that their orthologs and ganoid scales were convergently lost in Acipenseriformes and Amiiformes.

Efficient CRISPR Mutagenesis in Sturgeon Demonstrates Its Utility in Large, Slow-Maturing Vertebrates

In the last decade, the CRISPR/Cas9 bacterial virus defense system has been adapted as a user-friendly, efficient, and precise method for targeted mutagenesis in eukaryotes. Though CRISPR/Cas9 has proven effective in a diverse range of organisms, it is still most often used to create mutant lines in lab-reared genetic model systems. However, one major advantage of CRISPR/Cas9 mutagenesis over previous gene targeting approaches is that its high efficiency allows the immediate generation of near-null mosaic mutants. This feature could potentially allow genotype to be linked to phenotype in organisms with life histories that preclude the establishment of purebred genetic lines; a group that includes the vast majority of vertebrate species. Of particular interest to scholars of early vertebrate evolution are several long-lived and slow-maturing fishes that diverged from two dominant modern lineages, teleosts and tetrapods, in the Ordovician, or before. These early-diverging or "basal" vertebrates include the jawless cyclostomes, cartilaginous fishes, and various non-teleost ray-finned fishes. In addition to occupying critical phylogenetic positions, these groups possess combinations of derived and ancestral features not seen in conventional model vertebrates, and thus provide an opportunity for understanding the genetic bases of such traits. Here we report successful use of CRISPR/Cas9 mutagenesis in one such non-teleost fish, sterlet Acipenser ruthenus, a small species of sturgeon. We introduced mutations into the genes Tyrosinase, which is needed for melanin production, and Sonic hedgehog, a pleiotropic developmental regulator with diverse roles in early embryonic patterning and organogenesis. We observed disruption of both loci and the production of consistent phenotypes, including both near-null mutants' various hypomorphs. Based on these results, and previous work in lamprey and amphibians, we discuss how CRISPR/Cas9 F0 mutagenesis may be successfully adapted to other long-lived, slow-maturing aquatic vertebrates and identify the ease of obtaining and injecting eggs and/or zygotes as the main challenges.

FishmiRNA: An evolutionarily supported microRNA annotation and expression database for ray-finned fishes

microRNAs are important post-transcriptional regulators of gene expression involved in countless biological processes and are widely studied across metazoans. While miRNA research continues to grow, the large community of fish miRNA researchers lacks exhaustive resources consistent among species. To fill this gap, we developed FishmiRNA, an evolutionarily supported microRNA annotation and expression database for ray-finned fishes: www.fishmirna.org. The self-explanatory database contains detailed, manually-curated miRNA annotations with orthology relationships rigorously established by sequence similarity and conserved syntenies, and expression data provided for each detected mature miRNA. In just few clicks, users can download the annotation and expression database in several convenient formats either in its entirety or a subset. Simple filters and BLAST search options also permit the simultaneous exploration and visual comparison of expression data for up to any ten mature miRNAs across species and organs. FishmiRNA was specifically designed for ease of use to reach a wide audience.

On the relationship between extant innate immune receptors and the evolutionary origins of jawed vertebrate adaptive immunity

For over half a century, deciphering the origins of the genomic loci that form the jawed vertebrate adaptive immune response has been a major topic in comparative immunogenetics. Vertebrate adaptive immunity relies on an extensive and highly diverse repertoire of tandem arrays of variable (V), diversity (D), and joining (J) gene segments that recombine to produce different immunoglobulin (Ig) and T cell receptor (TCR) genes. The current consensus is that a recombination-activating gene (RAG)-like transposon invaded an exon of an ancient innate immune VJ-bearing receptor, giving rise to the extant diversity of Ig and TCR loci across jawed vertebrates. However, a model for the evolutionary relationships between extant non-recombining innate immune receptors and the V(D)J receptors of the jawed vertebrate adaptive immune system has only recently begun to come into focus. In this review, we provide an overview of non-recombining VJ genes, including CD8β, CD79b, natural cytotoxicity receptor 3 (NCR3/NKp30), putative remnants of an antigen receptor precursor (PRARPs), and the multigene family of signal-regulatory proteins (SIRPs), that play a wide range of roles in immune function. We then focus in detail on the VJ-containing novel immune-type receptors (NITRs) from ray-finned fishes, as recent work has indicated that these genes are at least 50 million years older than originally thought. We conclude by providing a conceptual model of the evolutionary origins and phylogenetic distribution of known VJ-containing innate immune receptors, highlighting opportunities for future comparative research that are empowered by this emerging evolutionary perspective.

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.

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.