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

Genome sequencing of Coryphaenoides yaquinae reveals convergent and lineage-specific molecular evolution in deep-sea adaptation

Abyssal (3501-6500 m) and hadal (>6500 m) fauna evolve under harsh abiotic stresses, characterized by high hydrostatic pressure, darkness and food shortage, providing unique opportunities to investigate mechanisms underlying environmental adaptation. Genomes of several hadal species have recently been reported. However, the genetic adaptation of deep sea species across a broad spectrum of ocean depths has yet to be thoroughly investigated, due to the challenges imposed by collecting the deep sea species. To elucidate the correlation between genetic innovation and vertical distribution, we generated a chromosome-level genome assembly of the macrourids Coryphaenoides yaquinae, which is widely distributed in the abyssal/hadal zone ranging from 3655 to 7259 m in depth. Genomic comparisons among shallow, abyssal and hadal-living species identified idiosyncratic and convergent genetic alterations underlying the extraordinary adaptations of deep-sea species including light perception, circadian regulation, hydrostatic pressure and hunger tolerance. The deep-sea fishes (Coryphaenoides Sp. and Pseudoliparis swirei) venturing into various ocean depths independently have undergone convergent amino acid substitutions in multiple proteins such as rhodopsin 1, pancreatic and duodenal homeobox 1 and melanocortin 4 receptor which are known or verified in zebrafish to be related with vision adaptation and energy expenditure. Convergent evolution events were also identified in heat shock protein 90 beta family member 1 and valosin-containing protein genes known to be related to hydrostatic pressure adaptation specifically in fishes found around the hadal range. The uncovering of the molecular convergence among the deep-sea species shed new light on the common genetic innovations required for deep-sea adaptation by the fishes.

A complete biosynthetic pathway of the long-chain polyunsaturated fatty acids in an amphidromous fish, ayu sweetfish Plecoglossus altivelis (Stomiati; Osmeriformes)

The biosynthetic capability of the long-chain polyunsaturated fatty acids (LC-PUFA) in teleosts are highly diversified due to evolutionary events such as gene loss and subsequent neo- and/or sub-functionalisation of enzymes encoded by existing genes. In the present study, we have comprehensively characterised genes potentially involved in LC-PUFA biosynthesis, namely one front-end desaturase (fads2) and eight fatty acid elongases (elovl1a, elovl1b, elovl4a, elovl4b, elovl5, elovl7, elovl8a and elovl8b) from an amphidromous teleost, Ayu sweetfish, Plecoglossus altivelis. Functional analysis confirmed Fads2 with Δ6, Δ5 and Δ8 desaturase activities towards multiple PUFA substrates and several Elovl enzymes exhibited elongation capacities towards C18-20 or C18-22 PUFA substrates. Consequently, P. altivelis possesses a complete enzymatic capability to synthesise physiologically important LC-PUFA including arachidonic acid (ARA, 20:4n-6), eicosapentaenoic acid (EPA, 20:5n-3) and docosahexaenoic acid (DHA, 22:6n-3) from their C18 precursors. Interestingly, the loss of elovl2 gene in P. altivelis was corroborated by genomic and phylogenetic analyses. However, this constraint would possibly be overcome by the function of alternative Elovl enzymes, such as Elovl1b, which has not hitherto been functionally characterised in teleosts. The present study contributes novel insights into LC-PUFA biosynthesis in the relatively understudied teleost group, Osmeriformes (Stomiati), thereby enhancing our understanding of the complement of LC-PUFA biosynthetic genes within teleosts.

The Hox code responsible for the patterning of the anterior vertebrae in zebrafish

The vertebral column is a characteristic structure of vertebrates. Genetic studies in mice have shown that Hox-mediated patterning plays a key role in specifying discrete anatomical regions of the vertebral column. Expression pattern analyses in several vertebrate embryos have provided correlative evidence that the anterior boundaries of Hox expression coincide with distinct anatomical vertebrae. However, because functional analyses have been limited to mice, it remains unclear which Hox genes actually function in vertebral patterning in other vertebrates. In this study, various zebrafish Hox mutants were generated for loss-of-function phenotypic analysis to functionally decipher the Hox code responsible for the zebrafish anterior vertebrae between the occipital and thoracic vertebrae. We found that Hox genes in HoxB- and HoxC-related clusters participate in regulating the morphology of the zebrafish anterior vertebrae. In addition, medaka hoxc6a was found to be responsible for anterior vertebral identity, as in zebrafish. Based on phenotypic similarities with Hoxc6 knockout mice, our results suggest that the Hox patterning system, including at least Hoxc6, may have been functionally established in the vertebral patterning of the common ancestor of ray-finned and lobe-finned fishes.

Scoutknife: A naïve, whole genome informed phylogenetic robusticity metric

Background: The phylogenetic bootstrap, first proposed by Felsenstein in 1985, is a critically important statistical method in assessing the robusticity of phylogenetic datasets. Core to its concept was the use of pseudo sampling - assessing the data by generating new replicates derived from the initial dataset that was used to generate the phylogeny. In this way, phylogenetic support metrics could overcome the lack of perfect, infinite data. With infinite data, however, it is possible to sample smaller replicates directly from the data to obtain both the phylogeny and its statistical robusticity in the same analysis. Due to the growth of whole genome sequencing, the depth and breadth of our datasets have greatly expanded and are set to only expand further. With genome-scale datasets comprising thousands of genes, we can now obtain a proxy for infinite data. Accordingly, we can potentially abandon the notion of pseudo sampling and instead randomly sample small subsets of genes from the thousands of genes in our analyses. Methods: We introduce Scoutknife, a jackknife-style subsampling implementation that generates 100 datasets by randomly sampling a small number of genes from an initial large-gene dataset to jointly establish both a phylogenetic hypothesis and assess its robusticity. We assess its effectiveness by using 18 previously published datasets and 100 simulation studies. Results: We show that Scoutknife is conservative and informative as to conflicts and incongruence across the whole genome, without the need for subsampling based on traditional model selection criteria. Conclusions: Scoutknife reliably achieves comparable results to selecting the best genes on both real and simulation datasets, while being resistant to the potential biases caused by selecting for model fit. As the amount of genome data grows, it becomes an even more exciting option to assess the robusticity of phylogenetic hypotheses.

Molecular Dating of the Teleost Whole Genome Duplication (3R) Is Compatible With the Expectations of Delayed Rediploidization

Vertebrate evolution has been punctuated by three whole genome duplication events that have been implicated causally in phenotypic evolution, from the origin of phenotypic novelties to explosive diversification. Arguably, the most dramatic of these is the 3R whole genome duplication event associated with the origin of teleost fishes which comprise more than half of all living vertebrate species. However, tests of a causal relationship between whole genome duplication and teleost diversification have proven difficult due to the challenge of establishing the timing of these phenomena. Here we show, based on molecular clock dating of concatenated gene alignments, that the 3R whole genome duplication event occurred in the early-middle Permian (286.18 to 267.20 million years ago; Ma), 52.02 to 12.84 million years (Myr) before the divergence of crown-teleosts in the latest Permian-earliest Late Triassic (254.36 to 234.16 Ma) and long before the major pulses of teleost diversification in Ostariophysi and Percomorpha (56.37 to 100.17 Myr and at least 139.24 to 183.29 Myr later, respectively). The extent of this temporal gap between putative cause and effect precludes 3R as a deterministic driver of teleost diversification. However, these age constraints remain compatible with the expectations of a prolonged rediploidization process following whole genome duplication which, through the effects of chromosome rearrangement and gene loss, remains a viable mechanism to explain the evolution of teleost novelties and diversification.

Characterization and expression profiles of toll-like receptor genes (TLR2 and TLR5) in immune tissues of hybrid yellow catfish under bacterial infection

The yellow catfish (Pelteobagrus fulvidraco) is one of the most economically important freshwater species in Asia. However, pathogenic bacterial infections often cause high rates of mortality and economic losses in practical aquaculture. Previous studies in mammals have shown that Toll-like receptor 2 (TLR2) and Toll-like receptor 5 (TLR5) are involved in the recognition of cell wall components such as lipopolysaccharides and flagella of various bacteria, thereby acting as key regulators in the innate immunity response. However, TLR2 and TLR5 in yellow catfish have not been characterized. In the present study, TLR2 and TLR5 were examined through comparative genomic approaches. The gene structure, collinearity, protein spatial structure, and phylogenetic relationships were compared with those in multiple representative vertebrates. Meanwhile, quantitative real-time PCR was conducted to explore transcriptional changes in TLR2 and TLR5 in immune tissues after infection with exogenous A. hydrophila and E. tarda. The results demonstrated the presence of TLR2 and TLR5 in yellow catfish. However, a systematic analysis showed that TLR2 was not associated with the arrangement of diverse neighboring genes. The expression of hybrid yellow catfish TLR2 transcripts in multiple tissues (including liver, spleen, kidney, and intestine) was significantly up-regulated after infection with A. hydrophila and E. tarda, suggesting that hybrid yellow catfish TLR2 and TLR5 may participate in the immune process. Taken together, the results indicate that TLR2 and TLR5 are conserved in terms of evolution and possess significant antibacterial activity as well as regulatory properties in immune-related tissues and thus play key roles in host defense against pathogen invasion.

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.

Phylogenomic analysis and molecular identification of true fruit flies

The family Tephritidae in the order Diptera, known as true fruit flies, are agriculturally important insect pests. However, the phylogenetic relationships of true fruit flies, remain controversial. Moreover, rapid identification of important invasive true fruit flies is essential for plant quarantine but is still challenging. To this end, we sequenced the genome of 16 true fruit fly species at coverage of 47-228×. Together with the previously reported genomes of nine species, we reconstructed phylogenetic trees of the Tephritidae using benchmarking universal single-copy ortholog (BUSCO), ultraconserved element (UCE) and anchored hybrid enrichment (AHE) gene sets, respectively. The resulting trees of 50% taxon-occupancy dataset for each marker type were generally congruent at 88% nodes for both concatenation and coalescent analyses. At the subfamily level, both Dacinae and Trypetinae are monophyletic. At the species level, Bactrocera dorsalis is more closely related to Bactrocera latifrons than Bactrocera tryoni. This is inconsistent with previous conclusions based on mitochondrial genes but consistent with recent studies based on nuclear data. By analyzing these genome data, we screened ten pairs of species-specific primers for molecular identification of ten invasive fruit flies, which PCR validated. In summary, our work provides draft genome data of 16 true fruit fly species, addressing the long-standing taxonomic controversies and providing species-specific primers for molecular identification of invasive fruit flies.

Teleost Hox code defines regional identities competent for the formation of dorsal and anal fins

The dorsal and anal fins can vary widely in position and length along the anterior-posterior axis in teleost fishes. However, the molecular mechanisms underlying the diversification of these fins remain unknown. Here, we used genetic approaches in zebrafish and medaka, in which the relative positions of the dorsal and anal fins are opposite, to demonstrate the crucial role of hox genes in the patterning of the teleost posterior body, including the dorsal and anal fins. By the CRISPR-Cas9-induced frameshift mutations and positional cloning of spontaneous dorsalfinless medaka, we show that various hox mutants exhibit the absence of dorsal or anal fins, or a stepwise posterior extension of these fins, with vertebral abnormalities. Our results indicate that multiple hox genes, primarily from hoxc-related clusters, encompass the regions responsible for the dorsal and anal fin formation along the anterior-posterior axis. These results further suggest that shifts in the anterior boundaries of hox expression which vary among fish species, lead to diversification in the position and size of the dorsal and anal fins, similar to how modulations in Hox expression can alter the number of anatomically distinct vertebrae in tetrapods. Furthermore, we show that hox genes responsible for dorsal fin formation are different between zebrafish and medaka. Our results suggest that a novel mechanism has occurred during teleost evolution, in which the gene network responsible for fin formation might have switched to the regulation downstream of other hox genes, leading to the remarkable diversity in the dorsal fin position.

Genetic Diversity and Population Structure of Rhodeus uyekii in the Republic of Korea Revealed by Microsatellite Markers from Whole Genome Assembly

This study is the first report to characterize the Rhodus uyekii genome and study the development of microsatellite markers and their markers applied to the genetic structure of the wild population. Genome assembly was based on PacBio HiFi and Illumina HiSeq paired-end sequencing, resulting in a draft genome assembly of R. uyekii. The draft genome was assembled into 2652 contigs. The integrity assessment of the assemblies indicates that the quality of the draft assemblies is high, with 3259 complete BUSCOs (97.2%) in the database of Verbrata. A total of 31,166 predicted protein-coding genes were annotated in the protein database. The phylogenetic tree showed that R. uyekii is a close but distinct relative of Onychostoma macrolepis. Among the 10 fish genomes, there were significant gene family expansions (8-2387) and contractions (16-2886). The average number of alleles amplified by the 21 polymorphic markers ranged from 6 to 23, and the average PIC value was 0.753, which will be useful for evolutionary and genetic analysis. Using population genetic analysis, we analyzed genetic diversity and the genetic structures of 120 individuals from 6 populations. The average number of alleles per population ranged from 7.6 to 9.9, observed heterozygosity ranged from 0.496 to 0.642, and expected heterozygosity ranged from 0.587 to 0.783. Discriminant analysis of principal components According to the analysis method, the population was divided into three populations (BS vs. DC vs. GG, GC, MS, DC). In conclusion, our study provides a useful resource for comparative genomics, phylogeny, and future population studies of R. uyekii.

Colonization of the ocean floor by jawless vertebrates across three mass extinctions

BACKGROUND

The deep (> 200 m) ocean floor is often considered to be a refugium of biodiversity; many benthic marine animals appear to share ancient common ancestry with nearshore and terrestrial relatives. Whether this pattern holds for vertebrates is obscured by a poor understanding of the evolutionary history of the oldest marine vertebrate clades. Hagfishes are jawless vertebrates that are either the living sister to all vertebrates or form a clade with lampreys, the only other surviving jawless fishes.

RESULTS

We use the hagfish fossil record and molecular data for all recognized genera to construct a novel hypothesis for hagfish relationships and diversification. We find that crown hagfishes persisted through three mass extinctions after appearing in the Permian ~ 275 Ma, making them one of the oldest living vertebrate lineages. In contrast to most other deep marine vertebrates, we consistently infer a deep origin of continental slope occupation by hagfishes that dates to the Paleozoic. Yet, we show that hagfishes have experienced marked body size diversification over the last hundred million years, contrasting with a view of this clade as morphologically stagnant.

CONCLUSION

Our results establish hagfishes as ancient members of demersal continental slope faunas and suggest a prolonged accumulation of deep sea jawless vertebrate biodiversity.

FishTEDB 2.0: an update fish transposable element (TE) database with new functions to facilitate TE research

We launched the initial version of FishTEDB in 2018, which aimed to establish an open-source, user-friendly, data-rich transposable element (TE) database. Over the past 5 years, FishTEDB 1.0 has gained approximately 10 000 users, accumulating more than 450 000 interactions. With the unveiling of extensive fish genome data and the increasing emphasis on TE research, FishTEDB needs to extend the richness of data and functions. To achieve the above goals, we introduced 33 new fish species to FishTEDB 2.0, encompassing a wide array of fish belonging to 48 orders. To make the updated database more functional, we added a genome browser to visualize the positional relationship between TEs and genes and the estimated TE insertion time in different species. In conclusion, we released a new version of the fish TE database, FishTEDB 2.0, designed to assist researchers in the future study of TE functions and promote the progress of biological theories related to TEs. Database URL: https://www.fishtedb.com/.

Genetic adaptations for the oceanic success of fish eggs

Genetic adaptations of organisms living in extreme environments are fundamental to our understanding of where life can evolve. Water is the single limiting parameter in this regard, yet when released in the oceans, the single-celled eggs of marine bony fishes (teleosts) have no means of acquiring it. They are strongly hyposmotic to seawater and lack osmoregulatory systems. Paradoxically, modern teleosts successfully release vast quantities of eggs in the extreme saline environment and recorded the most explosive radiation in vertebrate history. Here, we highlight key genetic adaptations that evolved to solve this paradox by filling the pre-ovulated eggs with water. The degree of water acquisition is uniquely prevalent to marine teleosts, permitting the survival and oceanic dispersal of their eggs.

Deploy the proboscis!: Functional morphology and kinematics of a novel form of extreme jaw protrusion in the hingemouth, Phractolaemus ansorgii (Gonorynchiformes)

Premaxillary protrusion and the performance advantages it confers are implicated in the success of diverse lineages of teleost fishes, such as Cypriniformes and Acanthomorpha. Although premaxillary protrusion has evolved independently at least five times within bony fishes, much of the functional work investigating this kinesis relates to mechanisms found only in these two clades. Few studies have characterized feeding mechanisms in less-diverse premaxilla-protruding lineages and fewer yet have investigated the distinctive anatomy underlying jaw kinesis in these lineages. Here, we integrated dissection, clearing and staining, histology, micro-CT, and high-speed videography to investigate an isolated and independent origin of jaw protrusion in the hingemouth, Phractolaemus ansorgii, which employs a complex arrangement of bones, musculature, and connective tissues to feed on benthic detritus via a deployable proboscis. Our goals were to provide an integrative account of the underlying architecture of P. ansorgii's feeding apparatus and to assess the functional consequences of this drastic deviation from the more typical teleost condition. Phractolaemus ansorgii's cranial anatomy is distinct from all other fishes in that its adducted lower jaw is caudally oriented, and it possesses a mouth at the terminal end of an elongated, tube-like proboscis that is unique in its lack of skeletal support from the oral jaws. Instead, its mouth is supported primarily by hyaline-cell cartilage and other rigid connective tissues, and features highly flexible lips that are covered in rows of keratinous unculi. Concomitant changes to the adductor musculature likely allow for the flexibility to protrude the mouth dorsally and ventrally as observed during different feeding behaviors, while the intrinsic compliance of the lips allows for more effective scraping of irregular surfaces. From our feeding videos, we find that P. ansorgii is capable of modulating the distance of protrusion, with maximum anterior protrusion exceeding 30% of head length. This represents a previously undescribed example of extreme jaw protrusion on par with many acanthomorph species. Protrusion is much slower in P. ansorgii-reaching an average speed of 2.74 cm/s-compared to acanthomorphs feeding on elusive prey or even benthivorous cypriniforms. However, this reorganization of cranial anatomy may reflect a greater need for dexterity to forage more precisely in multiple directions and on a wide variety of surface textures. Although this highly modified mechanism may have limited versatility over evolutionary timescales, it has persisted in solitude within Gonorynchiformes, representing a novel functional solution for benthic feeding in tropical West African rivers.

High Diversity of Long Terminal Repeat Retrotransposons in Compact Vertebrate Genomes: Insights from Genomes of Tetraodontiformes

This study aimed to investigate the evolutionary profile (including diversity, activity, and abundance) of retrotransposons (RTNs) with long terminal repeats (LTRs) in ten species of Tetraodontiformes. These species, Arothron firmamentum, Lagocephalus sceleratus, Pao palembangensis, Takifugu bimaculatus, Takifugu flavidus, Takifugu ocellatus, Takifugu rubripes, Tetraodon nigroviridis, Mola mola, and Thamnaconus septentrionalis, are known for having the smallest genomes among vertebrates. Data mining revealed a high diversity and wide distribution of LTR retrotransposons (LTR-RTNs) in these compact vertebrate genomes, with varying abundances among species. A total of 819 full-length LTR-RTN sequences were identified across these genomes, categorized into nine families belonging to four different superfamilies: ERV (Orthoretrovirinae and Epsilon retrovirus), Copia, BEL-PAO, and Gypsy (Gmr, Mag, V-clade, CsRN1, and Barthez). The Gypsy superfamily exhibited the highest diversity. LTR family distribution varied among species, with Takifugu bimaculatus, Takifugu flavidus, Takifugu ocellatus, and Takifugu rubripes having the highest richness of LTR families and sequences. Additionally, evidence of recent invasions was observed in specific tetraodontiform genomes, suggesting potential transposition activity. This study provides insights into the evolution of LTR retrotransposons in Tetraodontiformes, enhancing our understanding of their impact on the structure and evolution of host genomes.

pipesnake: generalized software for the assembly and analysis of phylogenomic datasets from conserved genomic loci

MOTIVATION

Phylogenetics has moved into the era of genomics, incorporating enormous volumes of data to study questions at both shallow and deep scales. With this increase in information, phylogeneticists need new tools and skills to manipulate and analyze these data. To facilitate these tasks and encourage reproducibility, the community is increasingly moving toward automated workflows.

RESULTS

Here we present pipesnake, a phylogenomics pipeline written in Nextflow for the processing, assembly, and phylogenetic estimation of genomic data from short-read sequences. pipesnake is an easy to use and efficient software package designed for this next era in phylogenetics.

AVAILABILITY AND IMPLEMENTATION

pipesnake is publicly available on GitHub at https://github.com/AusARG/pipesnake and accompanied by documentation and a wiki/tutorial.

Chromosome-level genome assembly and population genomic analysis reveal evolution and local adaptation in common hairfin anchovy (Setipinna tenuifilis)

Understanding the genetic structure and the factors associated with adaptive diversity has significant implications for the effective management of wild populations under threat from overfishing and climate change. The common hairfin anchovy (Setipinna tenuifilis) is an economically and ecologically important pelagic fish species, spanning a broad latitudinal gradient along marginal seas of the Northwest Pacific. In this study, we constructed the first reference genome of S. tenuifilis using PacBio long reads and high-resolution chromosome conformation capture (Hi-C) technology. The assembled genome was 798.38 Mb with a contig N50 of 1.43 Mb and a scaffold N50 of 32.42 Mb, which were anchored onto 24 pseudochromosomes. A total of 22,019 genes were functionally annotated, which accounted for 95.27% of the predicted protein-coding genes. Chromosomal collinearity analysis revealed chromosome fusion or fission events in Clupeiformes species. Three genetic groups of S. tenuifilis were revealed along the Chinese coast using restriction site-associated DNA sequencing (RADseq). We investigated the influence of four bioclimatic variables as potential drivers of adaptive divergence in S. tenuifilis, suggesting that these environmental variables, especially sea surface temperature, may play important roles as drivers of spatially varying selection for S. tenuifilis. We also identified candidate functional genes underlying adaptive mechanisms and ecological tradeoffs using redundancy analysis (RDA) and BayeScan analysis. In summary, this study sheds light on the evolution and spatial patterns of genetic variation of S. tenuifilis, providing a valuable genomic resource for further biological and genetic studies on this species and other closely related Clupeiformes.

The promise and pitfalls of synteny in phylogenomics

Reconstructing the tree of life remains a central goal in biology. Early methods, which relied on small numbers of morphological or genetic characters, often yielded conflicting evolutionary histories, undermining confidence in the results. Investigations based on phylogenomics, which use hundreds to thousands of loci for phylogenetic inquiry, have provided a clearer picture of life's history, but certain branches remain problematic. To resolve difficult nodes on the tree of life, 2 recent studies tested the utility of synteny, the conserved collinearity of orthologous genetic loci in 2 or more organisms, for phylogenetics. Synteny exhibits compelling phylogenomic potential while also raising new challenges. This Essay identifies and discusses specific opportunities and challenges that bear on the value of synteny data and other rare genomic changes for phylogenomic studies. Synteny-based analyses of highly contiguous genome assemblies mark a new chapter in the phylogenomic era and the quest to reconstruct the tree of life.

Morphological Diversity and Evolution of Jaw Morphologies in Zeiform Fishes (Teleostei, Paracanthopterygii)

Zeiformes (dories, tinselfishes, and oreos) are primarily benthopelagic acanthomorph fishes, distributed between 50 and 1000 m depth on continental slopes and on flanks of oceanic islands and seamounts. Among the interesting morphological adaptations of zeiform fishes are their unique and highly protrusible jaws involving premaxillae with long ascending processes and a four-bar linkage, including mobile palatines that pivot on their posterior articulation. This adaptation for increased jaw protrusion has enabled zeiform fishes to capture elusive prey more efficiently and is arguably a major factor in their morphological diversity and evolutionary success. This study examines the evolution of zeiform jaw morphologies using 3D landmark-based multivariate morphometrics as well as phylomorphospace analysis. Results show that the descendants of the zeiform ancestor branched rapidly early in their history, retaining conservative jaw morphologies during this early branching, but subsequently strongly diverged in many of the resulting lineages. Results from this study are compared with earlier research based on overall body form, demonstrating that morphological variation within Zeiformes arose along at least two distinct trajectories: body form and jaw morphology. Variation among genera in body form is not associated with variation among the same genera in jaw morphology, and vice versa. Hypotheses to explain the apparent decoupling of body shape and jaw morphology are addressed along with avenues for further study to better understand the morphological evolution of these iconic fishes.

Evolution of the enhancer-rich regulatory region of the gene for the cell-type specific transcription factor POU1F1

Precise spatio-temporal expression of genes in organogenesis is regulated by the coordinated interplay of DNA elements such as promoter and enhancers present in the regulatory region of a given locus. POU1F1 transcription factor plays a crucial role in the development of somatotrophs, lactotrophs and thyrotrophs in the anterior pituitary gland, and in maintaining high expression of growth hormone, prolactin and TSH. In mouse, expression of POU1F1 is controlled by a region fenced by two CTCF sites, containing 5 upstream enhancer elements, designated E-A (5' to 3'). Elements C, B and A correspond to elements shown previously to play a role in pituitary development and hormonal expression; functional roles for elements E and D have not been reported. We performed comparative sequence analysis of this regulatory region and discovered that three elements, B, C and E, are present in all vertebrate groups except Agnatha. One very long (>2 kb) element (A) is unique to mammals suggesting a specific change in regulation of the gene in this group. Using DNA accessibility assay (ATAC-seq) we showed that conserved elements in anterior pituitary of four non-mammals are open, suggesting functionality as regulatory elements. We showed that, in many non-mammalian vertebrates, an additional upstream exon closely follows element E, leading to alternatively spliced transcripts. Here, element E functions as an alternative promoter, but in mammals this feature is lost, suggesting conversion of alternative promoter to enhancer. Our work shows that regulation of POU1F1 changed markedly during the course of vertebrate evolution, use of a low number of enhancer elements combined with alternative promoters in non-mammalian vertebrates being replaced by use of a unique combination of regulatory units in mammals. Most importantly, our work suggests that evolutionary conversion of alternate promoter to enhancer could be one of the evolutionary mechanisms of enhancer birth.

Single Cell RNA Sequencing Provides Clues for the Developmental Genetic Basis of Syngnathidae's Evolutionary Adaptations

Seahorses, pipefishes, and seadragons are fishes from the family Syngnathidae that have evolved extraordinary traits including male pregnancy, elongated snouts, loss of teeth, and dermal bony armor. The developmental genetic and cellular changes that led to the evolution of these traits are largely unknown. Recent syngnathid genomes revealed suggestive gene content differences and provide the opportunity for detailed genetic analyses. We created a single cell RNA sequencing atlas of Gulf pipefish embryos to understand the developmental basis of four traits: derived head shape, toothlessness, dermal armor, and male pregnancy. We completed marker gene analyses, built genetic networks, and examined spatial expression of select genes. We identified osteochondrogenic mesenchymal cells in the elongating face that express regulatory genes bmp4, sfrp1a, and prdm16. We found no evidence for tooth primordia cells, and we observed re-deployment of osteoblast genetic networks in developing dermal armor. Finally, we found that epidermal cells expressed nutrient processing and environmental sensing genes, potentially relevant for the brooding environment. The examined pipefish evolutionary innovations are composed of recognizable cell types, suggesting derived features originate from changes within existing gene networks. Future work addressing syngnathid gene networks across multiple stages and species is essential for understanding how their novelties evolved.

An investigation of gill and blood carbonic anhydrase characteristics in three basal actinopterygian species: alligator gar (Atractosteus spatula), white sturgeon (Acipenser transmontanus) and Senegal bichir (Polypterus senegalus)

Many teleosts possess a unique set of respiratory characteristics allowing enhanced oxygen unloading to the tissues during stress. This system comprises three major components: highly pH sensitive haemoglobins (large Bohr and Root effects), rapid red blood cell (RBC) intracellular pH (pHi) protection, and a heterogeneous distribution of membrane-bound plasma-accessible carbonic anhydrase (paCA; absence in the gills). The first two components have received considerable research effort; however, the evolutionary loss of branchial paCA has received little attention. In the current study, we investigated the availability of branchial membrane-bound CA, along with several other CA-related characteristics in species belonging to three basal actinopterygian groups: the Lepisosteiformes, Acipenseriformes and Polypteriformes to assess the earlier hypothesis that Root effect haemoglobins constrain branchial paCA availability. We present the first evidence suggesting branchial membrane-bound CA presence in a basal actinopterygian species: the Senegal bichir (Polypterus senegalus) and show that like the teleosts, white sturgeon (Acipenser transmontanus) and alligator gar (Atractosteus spatula) do not possess branchial membrane-bound CA. We discuss the varying respiratory strategies for these species and propose that branchial paCA may have been lost much earlier than previously thought, likely in relation to the changes in haemoglobin buffer capacity associated with the increasing magnitude of the Bohr effect. The findings described here represent an important advancement in our understanding of the evolution of the unique system of enhanced oxygen unloading thought to be present in most teleosts, a group that encompasses half of all vertebrates.

Toward transparent taxonomy: an interactive web-tool for evaluating competing taxonomic arrangements

Informative and consistent taxonomy above the species level is essential to communication about evolution, biodiversity and conservation, and yet the practice of taxonomy is considered opaque and subjective by non-taxonomist scientists and the public alike. While various proposals have tried to make the basis for the ranking and inclusiveness of taxa more transparent and objective, widespread adoption of these ideas has lagged. Here, we present TaxonomR, an interactive online decision-support tool to evaluate alternative taxonomic classifications. This tool implements an approach that quantifies the criteria commonly used in taxonomic treatments and allows the user to interactively manipulate weightings for different criteria to compare scores for taxonomic groupings under those weights. We use the butterfly taxon Argynnis to demonstrate how different weightings applied to common taxonomic criteria result in fundamentally different genus-level classifications that are predominantly used in different continents and geographic regions. These differences are objectively compared and quantified using TaxonomR to evaluate the kinds of criteria that have been emphasized in earlier classifications, and the nature of the support for current alternative taxonomic arrangements. The main role of TaxonomR is to make taxonomic decisions transparent via an explicit prioritization scheme. TaxonomR is not a prescriptive application. Rather, it aims to be a tool for facilitating our understanding of alternative taxonomic classifications that can, in turn, potentially support global harmony in biodiversity assessments through evidence-based discussion and community-wide resolution of historically entrenched taxonomic tensions.

Lessons from assembling UCEs: A comparison of common methods and the case of Clavinomia (Halictidae)

Sequence data assembly is a foundational step in high-throughput sequencing, with untold consequences for downstream analyses. Despite this, few studies have interrogated the many methods for assembling phylogenomic UCE data for their comparative efficacy, or for how outputs may be impacted. We study this by comparing the most commonly used assembly methods for UCEs in the under-studied bee lineage Nomiinae and a representative sampling of relatives. Data for 63 UCE-only and 75 mixed taxa were assembled with five methods, including ABySS, HybPiper, SPAdes, Trinity and Velvet, and then benchmarked for their relative performance in terms of locus capture parameters and phylogenetic reconstruction. Unexpectedly, Trinity and Velvet trailed the other methods in terms of locus capture and DNA matrix density, whereas SPAdes performed favourably in most assessed metrics. In comparison with SPAdes, the guided-assembly approach HybPiper generally recovered the highest quality loci but in lower numbers. Based on our results, we formally move Clavinomia to Dieunomiini and render Epinomia once more a subgenus of Dieunomia. We strongly advise that future studies more closely examine the influence of assembly approach on their results, or, minimally, use better-performing assembly methods such as SPAdes or HybPiper. In this way, we can move forward with phylogenomic studies in a more standardized, comparable manner.

Macroevolutionary dynamics of gene family gain and loss along multicellular eukaryotic lineages

The gain and loss of genes fluctuate over evolutionary time in major eukaryotic clades. However, the full profile of these macroevolutionary trajectories is still missing. To give a more inclusive view on the changes in genome complexity across the tree of life, here we recovered the evolutionary dynamics of gene family gain and loss ranging from the ancestor of cellular organisms to 352 eukaryotic species. We show that in all considered lineages the gene family content follows a common evolutionary pattern, where the number of gene families reaches the highest value at a major evolutionary and ecological transition, and then gradually decreases towards extant organisms. This supports theoretical predictions and suggests that the genome complexity is often decoupled from commonly perceived organismal complexity. We conclude that simplification by gene family loss is a dominant force in Phanerozoic genomes of various lineages, probably underpinned by intense ecological specializations and functional outsourcing.

Structure-activity relationships for alkyl-phenanthrenes support two independent but interacting synergistic models for PAC mixture potency

Multiple lines of evidence at whole animal, cellular and molecular levels implicate polycyclic aromatic compounds (PACs) with three rings as drivers of crude oil toxicity to developing fish. Phenanthrene (P0) and its alkylated homologs (C1- through C4-phenanthrenes) comprise the most prominent subfraction of tricyclic PACs in crude oils. Among this family, P0 has been studied intensively, with more limited detail available for the C4-phenanthrene 1-methyl-7-isopropyl-phenanthrene (1-M,7-IP, or retene). While both compounds are cardiotoxic, P0 impacts embryonic cardiac function and development through direct blockade of K+ and Ca2+ currents that regulate cardiomyocyte contractions. In contrast, 1-M,7-IP dysregulates aryl hydrocarbon receptor (AHR) activation in developing ventricular cardiomyocytes. Although no other compounds have been assessed in detail across the larger family of alkylated phenanthrenes, increasing alkylation might be expected to shift phenanthrene family member activity from K+/Ca2+ ion current blockade to AHR activation. Using embryos of two distantly related fish species, zebrafish and Atlantic haddock, we tested 14 alkyl-phenanthrenes in both acute and latent developmental cardiotoxicity assays. All compounds were cardiotoxic, and effects were resolved into impacts on multiple, highly specific aspects of heart development or function. Craniofacial defects were clearly linked to developmental cardiotoxicity. Based on these findings, we suggest a novel framework to delineate the developmental toxicity of petrogenic PAC mixtures in fish, which incorporates multi-mechanistic pathways that produce interactive synergism at the organ level. In addition, relationships among measured embryo tissue concentrations, cytochrome P4501A mRNA induction, and cardiotoxic responses suggest a two-compartment toxicokinetic model that independently predicts high potency of PAC mixtures through classical metabolic synergism. These two modes of synergism, specific to the sub-fraction of phenanthrenes, are sufficient to explain the high embryotoxic potency of crude oils, independent of as-yet unmeasured compounds in these complex environmental mixtures.

Solving an enigma in the tree of life, at the origins of teleost fishes

Tracing the phylogenetic relationships between species is one of the fundamental objectives of evolutionary biology. Since Charles Darwin's seminal work in the 19th century, considerable progress has been made towards establishing a tree of life that summarises the evolutionary history of species. Nevertheless, substantial uncertainties still remain. Specifically, the relationships at the origins of teleost fishes have been the subject of extensive debate over the last 50 years. This question has major implications for various research fields: there are almost 30,000 species in the teleost group, which includes invaluable model organisms for biomedical, evolutionary and ecological studies. Here, we present the work in which we solved this enigma. We demonstrated that eels are more closely related to bony-tongued fishes than to the rest of teleost fishes. We achieved this by taking advantage of new genomic data and leveraging innovative phylogenetic markers. Notably, in addition to traditional molecular phylogeny methods based on the evolution of gene sequences, we also considered the evolution of gene order along the DNA molecule. We discuss the challenges and opportunities that these new markers represent for the field of molecular phylogeny, and in particular the possibilities they offer for re-examining other controversial branches in the tree of life.

Full-length transcriptome reveals the circularly polarized light response-related molecular genetic characteristics of Oratosquilla oratoria

The mantis shrimp is the only animal that can recognize circularly polarized light (CPL), but its molecular genetic characteristics are unclear. Multi-tissue level full-length (FL) transcriptome sequencing of Oratosquilla oratoria, a representative widely distributed mantis shrimp, was performed in the present study. We used comparative transcriptomics to explore the critical genes of O. oratoria selected by CPL and the GNβ gene associated with CPL signal transduction was hypothesized to be positively selected. Furthermore, the FL transcriptomes of O. oratoria compound eyes under five light conditions were sequenced and used to detect alternative splicing (AS). The ASs associated with CPL recognition mainly occurred in the LWS, ARR and TRPC regions. The number of FL transcripts with AS events and annotation information also provided evidence that O. oratoria could recognize LCPL. Additionally, 51 sequences belonging to the LWS, UV and Peropsin gene families were identified based on conserved 7tm domains. The LWS, UV and Peropsin opsins have similar 3D structures with seven domains across the cell membrane and conserved KSLRTPSN, DRY, and QAKK motifs. In conclusion, these results are undoubtedly valuable for perfecting the vision theory of O. oratoria and other mantis shrimp.

Identification, characterization and the inflammatory regulating effect of NOD1/2 in sturgeon

As an ancient species with both conservation and commercial value, Sturgeon's inflammatory regulation mechanism is a research point. Nucleotide-binding and oligomerization domain-containing proteins 1 and 2 (NOD1/2) are classical intracellular pattern recognition receptors (PRRs) in immunity of anti-bacterial infection. However, the characterization and function of NOD1/2 in Sturgeon are still unclear. In this study, we analyzed the synteny relationship of NOD1/2 genes between Acipenser ruthenus and representative fishes at the genome-level. Results showed that the ArNOD2 collinear genes pair was present in all representative fishes. The duplicated ArNOD1/2 genes were under purifying selection during evolution as indicated by their Ka/Ks values. To explore the function of NOD1/2, we further investigated their expression patterns and the effects of pathogenic infection, PAMPs treatment, and siRNA interference in Acipenser baerii, the sibling species of A. ruthenus. Results showed that both AbNOD1/2 were expressed at early developmental stages and in different tissues. Pathogenic infection in vivo and PAMPs treatment in vitro demonstrated that AbNOD1/2 could respond to pathogen stimulation. siRNA interference with AbNOD1/2 inhibited expression levels of RIPK2 and inflammatory cytokines compared to the control group after iE-DAP or MDP treatment. This study hinted that the AbNOD1/2 could stimulate the inflammatory cytokines response during evolutionary processes.

In silico insights into intra- and inter-species interactions of piscine gonadotropin hormones and receptor crosstalk

In mammals, the gonadotropins follicle-stimulating hormone (FSH) and luteinizing hormone (LH) are macromolecules secreted during specific reproductive phases and display strict specificity towards their cognate receptors. However, fish gonadotropins (GTH) and their receptors (GTHR) display diverse species-specific expression patterns, secretion patterns, and intra- and interspecies cross-activation. To uncover the molecular basis of this diversity, we generated and analyzed 29 in-silico models of intra- and inter-species combinations of sturgeon, carp, tilapia, and human gonadotropins with piscine receptors and analyzed the resulting receptor activation and signal transduction of these GTHR-GTH complexes in-vitro. Our results suggest that unlike humans, the surface charge on piscine FSH/LH β-seatbelt and N107huLHCGR/K104hFSHR homologs does not necessarily determine binding specificity. Instead, sequence and structural variations allow piscine GTHs significant conformational flexibility when binding to the receptor extracellular domain, thereby enabling cross-activation. The resulting diversity may support various reproductive strategies in different environmental niches.

Mitochondrial genome structure and composition in 70 fishes: a key resource for fisheries management in the South Atlantic

BACKGROUND

Phylogenetic gaps of public databases of reference sequences are a major obstacle for comparative genomics and management of marine resources, particularly in the Global South, where economically important fisheries and conservation flagship species often lack closely-related references. We applied target-enrichment to obtain complete mitochondrial genomes of marine ichthyofauna from the Brazilian coast selected based on economic significance, conservation status and lack of phylogenetically-close references. These included sardines (Dorosomatidae, Alosidae), mackerels (Scombridae) croakers (Sciaenidae), groupers (Epinephelidae) and snappers (Lutjanidae).

RESULTS

Custom baits were designed to enrich mitochondrial DNA across a broad phylogenetic range of fishes. Sequencing generated approximately 100k reads per sample, which were assembled in a total of 70 complete mitochondrial genomes and include fifty-two new additions to GenBank, including five species with no previous mitochondrial data. Departures from the typical gene content and order occurred in only three taxa and mostly involved tRNA gene duplications. Start-codons for all genes, except Cytochrome C Oxidase subunit I (COI), were consistently ATG, whilst a wide range of stop-codons deviated from the prevailing TAA. Phylogenetic analysis confirmed assembly accuracy and revealed signs of cryptic diversification within the Mullus genus. Lineage delimitation methods using Sardinella aurita and S. brasiliensis mitochondrial genomes support a single Operational Taxonomic Unit.

CONCLUSIONS

Target enrichment was highly efficient, providing complete novel mitochondrial genomes with little sequencing effort. These sequences are deposited in public databases to enable subsequent studies in population genetics and adaptation of Latin American fish species and serve as a vital resource for conservation and management programs that rely on molecular data for species and genus-level identification.

Chromosome-level genome assembly of the silver pomfret Pampus argenteus

Pampus argenteus (Euphrasen, 1788) is one of the major fishery species in coastal China. Pampus argenteus has a highly specialized morphology, and its declining fishery resources have encouraged massive research efforts on its aquacultural biology. In this study, we reported the first high-quality chromosome-level genome of P. argenteus obtained by integrating Illumina, PacBio HiFi, and Hi-C sequencing techniques. The final size of the genome was 518.06 Mb, with contig and scaffold N50 values of 20.47 and 22.86 Mb, respectively. The sequences were anchored and oriented onto 24 pseudochromosomes based on Hi-C data corresponding to the 24-chromatid karyotype of P. argenteus. A colinear relationship was observed between the P. argenteus genome and that of a closely related species (Scomber japonicus). A total of 24,696 protein-coding genes were identified from the genome, 98.9% of which were complete BUSCOs. This report represents the first case of high-quality chromosome-level genome assembly for P. argenteus and can provide valuable information for future evolutionary, conservation, and aquacultural research.

Conservation of Affinity Rather Than Sequence Underlies a Dynamic Evolution of the Motif-Mediated p53/MDM2 Interaction in Ray-Finned Fishes

The transcription factor and cell cycle regulator p53 is marked for degradation by the ubiquitin ligase MDM2. The interaction between these 2 proteins is mediated by a conserved binding motif in the disordered p53 transactivation domain (p53TAD) and the folded SWIB domain in MDM2. The conserved motif in p53TAD from zebrafish displays a 20-fold weaker interaction with MDM2, compared to the interaction in human and chicken. To investigate this apparent difference, we tracked the molecular evolution of the p53TAD/MDM2 interaction among ray-finned fishes (Actinopterygii), the largest vertebrate clade. Intriguingly, phylogenetic analyses, ancestral sequence reconstructions, and binding experiments showed that different loss-of-affinity changes in the canonical binding motif within p53TAD have occurred repeatedly and convergently in different fish lineages, resulting in relatively low extant affinities (KD = 0.5 to 5 μM). However, for 11 different fish p53TAD/MDM2 interactions, nonconserved regions flanking the canonical motif increased the affinity 4- to 73-fold to be on par with the human interaction. Our findings suggest that compensating changes at conserved and nonconserved positions within the motif, as well as in flanking regions of low conservation, underlie a stabilizing selection of "functional affinity" in the p53TAD/MDM2 interaction. Such interplay complicates bioinformatic prediction of binding and calls for experimental validation. Motif-mediated protein-protein interactions involving short binding motifs and folded interaction domains are very common across multicellular life. It is likely that the evolution of affinity in motif-mediated interactions often involves an interplay between specific interactions made by conserved motif residues and nonspecific interactions by nonconserved disordered regions.

A phylogenetics-based nomenclature system for steroid receptors in teleost fishes

Teleost fishes have emerged as tractable models for studying the neuroendocrine regulation of social behavior via molecular genetic techniques, such as CRISPR/Cas9 gene editing. Moreover, teleosts provide an opportunity to investigate the evolution of steroid receptors and their functions, as species within this lineage possess novel steroid receptor paralogs that resulted from a teleost-specific whole genome duplication. Although teleost fishes have grown in popularity as models for behavioral neuroendocrinology, there is not a consistent nomenclature system for steroid receptors and their genes, which may impede a clear understanding of steroid receptor paralogs and their functions. Here, we used a phylogenetic approach to assess the relatedness of protein sequences encoding steroid receptor paralogs in 18 species from 12 different orders of the Infraclass Teleostei. While most similarly named sequences grouped based on the established phylogeny of the teleost lineage, our analysis revealed several inconsistencies in the nomenclature of steroid receptor paralogs, particularly for sequences encoding estrogen receptor beta (ERβ). Based on our results, we propose a nomenclature system for teleosts in which Greek symbols refer to proteins and numbers refer to genes encoding different subtypes of steroid receptors within the five major groups of this nuclear receptor subfamily. Collectively, our results bridge a critical gap by providing a cohesive naming system for steroid receptors in teleost fishes, which will serve to improve communication, promote collaboration, and enhance our understanding of the evolution and function of steroid receptors across vertebrates.

The serotonin blocker Ketanserin reduces coral reef fish Ctenochaetus striatus aggressive behaviour during between-species social interactions

A multitude of species engages in social interactions not only with their conspecifics but also with other species. Such interspecific interactions can be either positive, like helping, or negative, like aggressive behaviour. However, the physiological mechanisms of these behaviours remain unclear. Here, we manipulated the serotonin system, a well-known neurohormone for regulating intraspecific aggressive behaviour, to investigate its role in interspecific aggression. We tested whether serotonin blockade affects the aggressive behaviour of a coral reef fish species (Ctenochaetus striatus) that engages in mutualistic interactions with another species, the cleaner fish (Labroides dimidiatus). Although this mutualistic cleaning relationship may appear positive, cleaner fish do not always cooperate and remove ectoparasites from the other coral reef fish ("clients") but tend to cheat and bite the client's protective layer of mucus. Client fish thus often apply control mechanisms, like chasing, to deter their cleaner fish partners from cheating. Our findings show that blocking serotonin receptors 5-HT2A and 5-HT2C with ketanserin reduced the client fish's aggressive behaviour towards cleaner fish, but in the context where the latter did not cheat. These results are evidence of the involvement of serotonin in regulating aggressive behaviour at the between-species social interactions level. Yet, the direction of effect we found here is the opposite of previous findings using a similar experimental set-up and ecological context but with a different client fish species (Scolopsis bilineatus). Together, it suggests that serotonin's role in aggressive behaviour is complex, and at least in this mutualistic ecological context, its function is species-dependent. This warrants, to some extent, careful interpretations from single-species studies looking into the physiological mechanisms of social behaviour.

Looking for a needle in a haystack: de novo phenotypic target identification reveals Hippo pathway-mediated miR-202 regulation of egg production

Understanding microRNA (miRNA) functions has been hampered by major difficulties in identifying their biological target(s). Currently, the main limitation is the lack of a suitable strategy to identify biologically relevant targets among a high number of putative targets. Here we provide a proof of concept of successful de novo (i.e. without prior knowledge of its identity) miRNA phenotypic target (i.e. target whose de-repression contributes to the phenotypic outcomes) identification from RNA-seq data. Using the medaka mir-202 knock-out (KO) model in which inactivation leads to a major organism-level reproductive phenotype, including reduced egg production, we introduced novel criteria including limited fold-change in KO and low interindividual variability in gene expression to reduce the list of 2853 putative targets to a short list of 5. We selected tead3b, a member of the evolutionarily-conserved Hippo pathway, known to regulate ovarian functions, due to its remarkably strong and evolutionarily conserved binding affinity for miR-202-5p. Deleting the miR-202-5p binding site in the 3' UTR of tead3b, but not of other Hippo pathway members sav1 and vgll4b, triggered a reduced egg production phenotype. This is one of the few successful examples of de novo functional assignment of a miRNA phenotypic target in vivo in vertebrates.

Disentangling historical relationships within Poeciliidae (Teleostei: Cyprinodontiformes) using ultraconserved elements

Poeciliids (Cyprinodontiformes: Poeciliidae), commonly known as livebearers, are popular fishes in the aquarium trade (e.g., guppies, mollies, swordtails) that are widely distributed in the Americas, with 274 valid species in 27 genera. This group has undergone various taxonomic changes recently, spurred by investigations using traditional genetic markers. Here we used over 1,000 ultraconserved loci to infer the relationships within Poeciliidae in the first attempt at understanding their diversification based on genome-scale data. We explore gene tree discordance and investigate potential incongruence between concatenation and coalescent inference methods. Our aim is to examine the influence of incomplete lineage sorting and reticulate evolution on the poeciliids' evolutionary history and how these factors contribute to the observed gene tree discordace. Our concatenated and coalescent phylogenomic inferences recovered four major clades within Poeciliidae. Most supra-generic level relationships we inferred were congruent with previous molecular studies, but we found some disagreements; the Middle American taxa Phallichthys and Poecilia (Mollienesia) were recovered as non-monophyletic, and unlike other recent molecular studies, we recovered Brachyrhaphis as monophyletic. Our study is the first to provide signatures of reticulate evolution in Poeciliidae at the family level; however, continued finer-scale investigations are needed to understand the complex evolutionary history of the family along with a much-needed taxonomic re-evaluation.

Modeling familial and sporadic Parkinson's disease in small fishes

The establishment of animal models for Parkinson's disease (PD) has been challenging. Nevertheless, once established, they will serve as valuable tools for elucidating the causes and pathogenesis of PD, as well as for developing new strategies for its treatment. Following the recent discovery of a series of PD causative genes in familial cases, teleost fishes, including zebrafish and medaka, have often been used to establish genetic PD models because of their ease of breeding and gene manipulation, as well as the high conservation of gene orthologs. Some of the fish lines can recapitulate PD phenotypes, which are often more pronounced than those in rodent genetic models. In addition, a new experimental teleost fish, turquoise killifish, can be used as a sporadic PD model, because it spontaneously manifests age-dependent PD phenotypes. Several PD fish models have already made significant contributions to the discovery of novel PD pathological features, such as cytosolic leakage of mitochondrial DNA and pathogenic phosphorylation in α-synuclein. Therefore, utilizing various PD fish models with distinct degenerative phenotypes will be an effective strategy for identifying emerging facets of PD pathogenesis and therapeutic modalities.

Gamete-exporting organs of vertebrates: dazed and confused

Mature gametes are transported externally for fertilization. In vertebrates, the gonads are located within the coelom. Consequently, each species has specific organs for export, which often vary according to sex. In most vertebrates, sperm ducts and oviducts develop from the Wolffian and Müllerian ducts, respectively. However, exceptions exist. Both sexes of cyclostomes, as well as females of basal teleosts, lack genital ducts but possess genital pores. In teleosts of both sexes, genital ducts are formed through the posterior extensions of gonads. These structures appear to be independent of both Wolffian and Müllerian ducts. Furthermore, the development of Wolffian and Müllerian ducts differs significantly among various vertebrates. Are these gamete-exporting organs homologous or not? A question extensively debated around the turn of the 20th century but now largely overlooked. Recent research has revealed the indispensable role of Wnt4a in genital duct development in both sexes of teleosts: zebrafish and medaka. wnt4a is an ortholog of mammalian Wnt4, which has functions in Müllerian duct formation. These results suggest a potential homology between the mammalian Müllerian ducts and genital ducts in teleosts. To investigate the homology of gamete-exporting organs in vertebrates, more detailed descriptions of their development across vertebrates, using modern cellular and genetic tools, are needed. Therefore, this review summarizes existing knowledge and unresolved questions on the structure and development of gamete-exporting organs in diverse vertebrate groups. This also underscores the need for comprehensive studies, particularly on cyclostomes, cartilaginous fishes, basal ray-finned fishes, and teleosts.

Anatomy of the endocrine pancreas in actinopterygian fishes and its phylogenetic implications

The anatomy and organisation of the endocrine pancreas in ray-finned fishes vary widely. The two main morphoanatomical character states are diffuse versus compact pancreatic tissue. The latter are called Brockmann Bodies (BBs), or principal islets. The present study is the first comprehensive survey on the anatomy of the endocrine pancreas (diffuse versus compact) across 322 actinopterygian species in 39 orders and 135 families based on literature, specimen dissections, and Magnetic Resonance Imaging (MRI). The data obtained show that large endocrine pancreatic islets (BB) have appeared several times in teleost evolution: in some ostariophysian clades and within the Salmoniformes and Neoteleostei. Acanthomorpha (spiny-rayed fishes) is the largest clade of the Neoteleostei. Within this clade, an absence of BBs is only observed in flying fishes (Exocoetidae), parrotfishes (Scarinae), and some of the scarine relatives, the Labridae. The presence of BBs in examined jellynose fish species from the Ateleopodiformes indicates support for its sister-group relationship to the Ctenosquamata (Myctophiformes + Acanthomorpha). More investigations are still needed to corroborate the presence or absence of BBs as a putative synapomorphy for a clade comprising Ateleopodiformes and Ctenosquamata.

Hypoxia-inducible factor 1 alpha protein increases without changes in mRNA during acute hypoxic exposure of the Gulf killifish, Fundulus grandis

The hypoxia inducible factor 1 (HIF1) is a central regulator of the molecular responses of animals to low oxygen. While the hypoxia-responsiveness of HIF1 is generally attributed to the stabilization of the alpha protein subunit (HIF1α) at low oxygen, several studies on fish report increased tissue levels of HIF1A mRNA during hypoxia, suggesting transcriptional regulation. In the current study, HIF1α protein and HIF1A mRNA were determined in parallel in tissues of Gulf killifish, Fundulus grandis, exposed to short-term hypoxia (24 h at 1 mg O2 l-1). HIF1α protein was higher in brain, ovary, and skeletal muscle from fish exposed to hypoxia compared with normoxic controls by 6 h, and it remained elevated in brain and ovary at 24 h. In contrast, HIF1A mRNA levels were unaffected by hypoxia in any tissue. Moreover, HIF1α protein and HIF1A mRNA levels in the same tissues were not correlated with one another during either normoxia or hypoxia. Hence, an increase in HIF1α protein does not depend upon an increase in HIF1A mRNA during acute exposure to low oxygen in this species. The results support the widely accepted mechanism of post-translational protein stabilization, rather than new transcription, during the initial response of fish to hypoxia.

The role of the 5' HoxA genes in the development of the hindgut, vent, and a novel sphincter in a derived teleost (bluebanded goby, Lythrypnus dalli)

Unique expression patterns of the 5' HoxA genes are associated with the evolution and development of novel features including claspers in cartilaginous fishes, modified pectoral fins in batoids, and the yolk sac extension in Cypriniformes. Here, we demonstrate a role for HoxA11a and HoxA13a in demarcating the hindgut in fishes of the family Gobiidae, including a novel sphincter called the intestinal rectal sphincter (IRS). Disruption of 5' HoxA expression, via manipulation of retinoic acid signaling, results in failure of the IRS and/or vent to develop. Furthermore, exposure to HoxA disruptors alters 5' HoxA expression, in association with developmental phenotypes, demonstrating a functional link between 5' HoxA expression and development of a novel feature in the bluebanded goby, Lythrypnus dalli.

Ontogenetic variation in the ecological stoichiometry of 10 fish species during early development

The chemical composition and stoichiometry of vertebrate bodies changes greatly during ontogeny as phosphorus-rich bones form, but we know little about the variation among species during early development. Such variation is important because element ratios in animal bodies influence which element limits growth and how animals contribute to nutrient cycling. We quantified ontogenetic variation from embryos through 2-3 months of age in 10 species of fish in six different families, ranging in adult size from 73 to 720 mm in length. We measured whole-body concentrations (percentage of dry mass) and ratios of carbon (C), nitrogen (N), and phosphorus (P) as fish developed. We also quantified whole-body concentrations of calcium (Ca), because Ca should reflect bone development, and RNA, which can be a major pool of body P. To account for interspecific differences in adult size, we also examined how trends changed with relative size, defined as body length divided by adult length. Ontogenetic changes in body composition and ratios were relatively similar among species and were more similar when expressed as a function of relative size compared to age. Body P increased rapidly in all species (likely because of bone development) from embryos until individuals were ~5%-8% of adult size. Body N also increased, while body C, C:N, C:P, and N:P all decreased over this period. Body Ca increased with development but was more variable among species. Body RNA was low in embryos, increased rapidly in young larvae, then decreased as fish reached 5%-8% of adult size. After fish were about 5%-8% of adult size, changes in body composition were relatively slight for all elements and ratios. These results reveal a consistency in the dynamics of body stoichiometry during early ontogeny, presumably because of similar constraints on the allocation of elements to bones and other body pools. Because most changes occur when individuals are <1 month old (<10% of adult size for that species), early ontogenetic variation in body stoichiometry may be especially important for growth limitation of individuals and ecosystem-level nutrient cycling.

Phylogenetic structure of body shape in a diverse inland ichthyofauna

Body shape is a fundamental metric of animal diversity affecting critical behavioral and ecological dynamics and conservation status, yet previously available methods capture only a fraction of total body-shape variance. Here we use structure-from-motion (SFM) 3D photogrammetry to generate digital 3D models of adult fishes from the Lower Mississippi Basin, one of the most diverse temperate-zone freshwater faunas on Earth, and 3D geometric morphometrics to capture morphologically distinct shape variables, interpreting principal components as growth fields. The mean body shape in this fauna resembles plesiomorphic teleost fishes, and the major dimensions of body-shape disparity are similar to those of other fish faunas worldwide. Major patterns of body-shape disparity are structured by phylogeny, with nested clades occupying distinct portions of the morphospace, most of the morphospace occupied by multiple distinct clades, and one clade (Acanthomorpha) accounting for over half of the total body shape variance. In contrast to previous studies, variance in body depth (59.4%) structures overall body-shape disparity more than does length (31.1%), while width accounts for a non-trivial (9.5%) amount of the total body-shape disparity.

Widespread sympatry in a species-rich clade of marine fishes (Carangoidei)

A universal paradigm describing patterns of speciation across the tree of life has been debated for decades. In marine organisms, inferring patterns of speciation using contemporary and historical patterns of biogeography is challenging due to the deficiency of species-level phylogenies and information on species' distributions, as well as conflicting relationships between species' dispersal, range size and co-occurrence. Most research on global patterns of marine fish speciation and biogeography has focused on coral reef or pelagic species. Carangoidei is an ecologically important clade of marine fishes that use coral reef and pelagic environments. We used sequence capture of 1314 ultraconserved elements (UCEs) from 154 taxa to generate a time-calibrated phylogeny of Carangoidei and its parent clade, Carangiformes. Age-range correlation analyses of the geographical distributions and divergence times of sister species pairs reveal widespread sympatry, with 73% of sister species pairs exhibiting sympatric geographical distributions, regardless of node age. Most species pairs coexist across large portions of their ranges. We also observe greater disparity in body length and maximum depth between sympatric relative to allopatric sister species. These and other ecological or behavioural attributes probably facilitate sympatry among the most closely related carangoids.

Ganoin and acrodin formation on scales and teeth in spotted gar: A vital role of enamelin in the unique process of enamel mineralization

Gars and bichirs develop scales and teeth with ancient actinopterygian characteristics. Their scale surface and tooth collar are covered with enamel, also known as ganoin, whereas the tooth cap is equipped with an enamel-like tissue, acrodin. Here, we investigated the formation and mineralization of the ganoin and acrodin matrices in spotted gar, and the evolution of the scpp5, ameloblastin (ambn), and enamelin (enam) genes, which encode matrix proteins of ganoin. Results suggest that, in bichirs and gars, all these genes retain structural characteristics of their orthologs in stem actinopterygians, presumably reflecting the presence of ganoin on scales and teeth. During scale formation, Scpp5 and Enam were initially found in the incipient ganoin matrix and the underlying collagen matrix, whereas Ambn was detected mostly in a surface region of the well-developed ganoin matrix. Although collagen is the principal acrodin matrix protein, Scpp5 was detected within the matrix. Similarities in timings of mineralization and the secretion of Scpp5 suggest that acrodin evolved by the loss of the matrix secretory stage of ganoin formation: dentin formation is immediately followed by the maturation stage. The late onset of Ambn secretion during ganoin formation implies that Ambn is not essential for mineral ribbon formation, the hallmark of the enamel matrix. Furthermore, Scpp5 resembles amelogenin that is not important for the initial formation of mineral ribbons in mammals. It is thus likely that the evolution of ENAM was vital to the origin of the unique mineralization process of the enamel matrix.

Resolving Marine-Freshwater Transitions by Diatoms Through a Fog of Gene Tree Discordance

Despite the obstacles facing marine colonists, most lineages of aquatic organisms have colonized and diversified in freshwaters repeatedly. These transitions can trigger rapid morphological or physiological change and, on longer timescales, lead to increased rates of speciation and extinction. Diatoms are a lineage of ancestrally marine microalgae that have diversified throughout freshwater habitats worldwide. We generated a phylogenomic data set of genomes and transcriptomes for 59 diatom taxa to resolve freshwater transitions in one lineage, the Thalassiosirales. Although most parts of the species tree were consistently resolved with strong support, we had difficulties resolving a Paleocene radiation, which affected the placement of one freshwater lineage. This and other parts of the tree were characterized by high levels of gene tree discordance caused by incomplete lineage sorting and low phylogenetic signal. Despite differences in species trees inferred from concatenation versus summary methods and codons versus amino acids, traditional methods of ancestral state reconstruction supported six transitions into freshwaters, two of which led to subsequent species diversification. Evidence from gene trees, protein alignments, and diatom life history together suggest that habitat transitions were largely the product of homoplasy rather than hemiplasy, a condition where transitions occur on branches in gene trees not shared with the species tree. Nevertheless, we identified a set of putatively hemiplasious genes, many of which have been associated with shifts to low salinity, indicating that hemiplasy played a small but potentially important role in freshwater adaptation. Accounting for differences in evolutionary outcomes, in which some taxa became locked into freshwaters while others were able to return to the ocean or become salinity generalists, might help further distinguish different sources of adaptive mutation in freshwater diatoms.

Evolution and phylogenetic distribution of endo-α-mannosidase

While glycans underlie many biological processes, such as protein folding, cell adhesion, and cell-cell recognition, deep evolution of glycosylation machinery remains an understudied topic. N-linked glycosylation is a conserved process in which mannosidases are key trimming enzymes. One of them is the glycoprotein endo-α-1,2-mannosidase which participates in the initial trimming of mannose moieties from an N-linked glycan inside the cis-Golgi. It is unique as the only endo-acting mannosidase found in this organelle. Relatively little is known about its origins and evolutionary history; so far it was reported to occur only in vertebrates. In this work, a taxon-rich bioinformatic survey to unravel the evolutionary history of this enzyme, including all major eukaryotic clades and a wide representation of animals, is presented. The endomannosidase was found to be more widely distributed in animals and other eukaryotes. The protein motif changes in context of the canonical animal enzyme were tracked. Additionally, the data show the two canonical vertebrate endomannosidase genes, MANEA and MANEAL, arose at the second round of the two vertebrate genome duplications and one more vertebrate paralog, CMANEAL, is uncovered. Finally, a framework where N-glycosylation co-evolved with complex multicellularity is described. A better understanding of the evolution of core glycosylation pathways is pivotal to understanding biology of eukaryotes in general, and the Golgi apparatus in particular. This systematic analysis of the endomannosidase evolution is one step toward this goal.

Marine and freshwater fishes of Alabama: a revised checklist and discussion of taxonomic issues

Checklists are fundamental and important tools for organizing information about biodiversity that provide a basis for conservation and additional scientific research. While Alabama is recognized as an aquatic biodiversity hotspot with the highest native freshwater fish diversity in the contiguous United States, we currently lack an up-to-date list of the states fishes. In particular, much has changed over the past ~20 years regarding our knowledge of fishes from Alabama and the Mobile River Basin, rendering past comprehensive treatments by Mettee et al. (1996) and Boschung and Mayden (2004) out of date. Here, we provide a revised checklist of marine and freshwater fishes known from the coastal and inland waters of Alabama that includes 463 species (335 primarily freshwater fishes, and 128 marine or diadromous fishes) in 35 orders, 78 families, and 176 genera. Extant, extirpated, and extinct species are included, as are putative candidate species. The checklist is based on prior work, searches of the literature and online sources, as well as parsing a large compilation of >140,000 fish records for Alabama and the Mobile River Basin from 37 data providers in the global Fishnet2 database (www.fishnet2.net) and >4000 marine survey records from the SEAMAP database (https://www.gsmfc.org/seamap.php). After editing and quality control checks, the final combined database contained 144,215 collection records, ~95% of which were georeferenced. We discuss the species descriptions, nomenclatural changes, and updates to marine species that account for changes to the state list, and we close with a discussion of ~13 candidate species forms that remain undescribed, which represent outstanding taxonomic issues in need of further research attention.

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.