Explosive diversification of marine fishes at the Cretaceous-Palaeogene boundary

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.

Synergistic innovations enabled the radiation of anglerfishes in the deep open ocean

Major ecological transitions are thought to fuel diversification, but whether they are contingent on the evolution of certain traits called key innovations1 is unclear. Key innovations are routinely invoked to explain how lineages rapidly exploit new ecological opportunities.1,2,3 However, investigations of key innovations often focus on single traits rather than considering trait combinations that collectively produce effects of interest.4 Here, we investigate the evolution of synergistic trait interactions in anglerfishes, which include one of the most species-rich vertebrate clades in the bathypelagic, or "midnight," zone of the deep sea: Ceratioidea.5 Ceratioids are the only vertebrates that possess sexual parasitism, wherein males temporarily attach or permanently fuse to females to mate.6,7 We show that the rapid transition of ancestrally benthic anglerfishes into pelagic habitats occurred during a period of major global warming 50-35 million years ago.8,9 This transition coincided with the origins of sexual parasitism, which is thought to increase the probability of successful reproduction once a mate is found in the midnight zone, Earth's largest habitat.5,6,7 Our reconstruction of the evolutionary history of anglerfishes and the loss of immune genes support that permanently fusing clades have convergently degenerated their adaptive immunity. We find that degenerate adaptive immune genes and sexual body size dimorphism, both variably present in anglerfishes outside the ceratioid radiation, likely promoted their transition into the bathypelagic zone. These results show how traits from separate physiological, morphological, and reproductive systems can interact synergistically to drive major transitions and subsequent diversification in novel environments.

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.

Ecological niche conservatism spurs diversification in response to climate change

Lengthy debate has surrounded the theoretical and empirical science of whether climatic niche evolution is related to increased or decreased rates of biological diversification. Because species can persist for thousands to millions of years, these questions cross broad scales of time and space. Thus, short-term experiments may not provide comprehensive understanding of the system, leading to the emergence of contrasting opinions: niche evolution may increase diversity by allowing species to explore and colonize new geographic areas across which they could speciate; or, niche conservatism might augment biodiversity by supporting isolation of populations that may then undergo allopatric speciation. Here, we use a simulation approach to test how biological diversification responds to different rates and modes of niche evolution. We find that niche conservatism promotes biological diversification, whereas labile niches-whether adapting to the conditions available or changing randomly-generally led to slower diversification rates. These novel results provide a framework for understanding how Earth-life interactions produced such a diverse biota.

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.

Fifth mass extinction event triggered the diversification of the largest family of freshwater gastropods (Caenogastropoda: Truncatelloidea: Hydrobiidae)

The fifth mass extinction event (MEE) at the Cretaceous-Palaeogene (K-Pg) boundary 66 million years ago (Ma) led to massive species loss but also triggered the diversification of higher taxa. Five models have been proposed depending on whether this diversification occurred before, during or after the K-Pg boundary and the rate of species accumulation. While the effects of the K-Pg MEE on vertebrate evolution are relatively well understood, the impact on invertebrates, particularly in freshwater ecosystems, remains controversial. One example is the hyperdiverse Hydrobiidae-the most species-rich family of freshwater gastropods. Whereas some studies place its origin in the Jurassic or even Carboniferous, most fossil records postdate the K-Pg event. We therefore used robustly time-calibrated multi-locus phylogenies of >400 species representing >100 hydrobiid genera to unravel its evolutionary history and patterns of diversification. We found that the family started diversifying shortly after the K-Pg boundary (∼60 Ma; 95% highest posterior density 52-69 Ma). Lineage richness gradually increased to the present and phylogenetic diversity until ∼25 Ma. These findings suggest that diversification was not initially driven by ecological opportunity. Combining the two criteria of timing and rate of diversification, a soft-explosive diversification model of aquatic vertebrates best fits the patterns observed. We also show that most higher hydrobiid taxa (i.e. subfamilies) diversified from the Middle Oligocene to Middle Miocene (i.e. 12-28 Ma). Two of the 15 major clades delimited are described here as new subfamilies (i.e. Bullaregiinae n. subfam. and Pontobelgrandiellinae n. subfam.), whose members are restricted to subterranean waters. Our results are an important contribution to understanding how the fifth MEE has shaped evolution and patterns of biodiversity in continental aquatic systems. Given the high extinction risks faced by many hydrobiids today, they also emphasise the need to study the biodiversity of vulnerable ecosystems.

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.

Whole-genome Comparisons Identify Repeated Regulatory Changes Underlying Convergent Appendage Evolution in Diverse Fish Lineages

Fins are major functional appendages of fish that have been repeatedly modified in different lineages. To search for genomic changes underlying natural fin diversity, we compared the genomes of 36 percomorph fish species that span over 100 million years of evolution and either have complete or reduced pelvic and caudal fins. We identify 1,614 genomic regions that are well-conserved in fin-complete species but missing from multiple fin-reduced lineages. Recurrent deletions of conserved sequences in wild fin-reduced species are enriched for functions related to appendage development, suggesting that convergent fin reduction at the organismal level is associated with repeated genomic deletions near fin-appendage development genes. We used sequencing and functional enhancer assays to confirm that PelA, a Pitx1 enhancer previously linked to recurrent pelvic loss in sticklebacks, has also been independently deleted and may have contributed to the fin morphology in distantly related pelvic-reduced species. We also identify a novel enhancer that is conserved in the majority of percomorphs, drives caudal fin expression in transgenic stickleback, is missing in tetraodontiform, syngnathid, and synbranchid species with caudal fin reduction, and alters caudal fin development when targeted by genome editing. Our study illustrates a broadly applicable strategy for mapping phenotypes to genotypes across a tree of vertebrate species and highlights notable new examples of regulatory genomic hotspots that have been used to evolve recurrent phenotypes across 100 million years of fish evolution.

Diversification of the sleepers (Gobiiformes: Gobioidei: Eleotridae) and evolution of the root gobioid families

Eleotridae (sleepers) and five smaller families are the earliest diverging lineages within Gobioidei. Most inhabit freshwaters in and around the Indo-Pacific, but Eleotridae also includes species that have invaded the Neotropics as well as several inland radiations in the freshwaters of Australia, New Zealand, and New Guinea. Previous efforts to infer phylogeny of these families have been based on sets of mitochondrial or nuclear loci and have yielded uncertain resolution of clades within Eleotridae. We expand the taxon sampling of previous studies and use genomic data from nuclear ultraconserved elements (UCEs) to infer phylogeny, then calibrate the hypothesis with recently discovered fossils. Our hypothesis clarifies ambiguously resolved relationships, provides a timescale for divergences, and indicates the core crown Eleotridae diverged over a short period 24.3-26.3 Ma in the late Oligocene. Within Eleotridae, we evaluate diversification dynamics with BAMM and find evidence for an overall slowdown in diversification over the past 35 Ma, but with a sharp increase 3.5 Ma in the genus Mogurnda, a clade of brightly colored species found in the freshwaters of Australia and New Guinea.

Genomics of cold adaptations in the Antarctic notothenioid fish radiation

Numerous novel adaptations characterise the radiation of notothenioids, the dominant fish group in the freezing seas of the Southern Ocean. To improve understanding of the evolution of this iconic fish group, here we generate and analyse new genome assemblies for 24 species covering all major subgroups of the radiation, including five long-read assemblies. We present a new estimate for the onset of the radiation at 10.7 million years ago, based on a time-calibrated phylogeny derived from genome-wide sequence data. We identify a two-fold variation in genome size, driven by expansion of multiple transposable element families, and use the long-read data to reconstruct two evolutionarily important, highly repetitive gene family loci. First, we present the most complete reconstruction to date of the antifreeze glycoprotein gene family, whose emergence enabled survival in sub-zero temperatures, showing the expansion of the antifreeze gene locus from the ancestral to the derived state. Second, we trace the loss of haemoglobin genes in icefishes, the only vertebrates lacking functional haemoglobins, through complete reconstruction of the two haemoglobin gene clusters across notothenioid families. Both the haemoglobin and antifreeze genomic loci are characterised by multiple transposon expansions that may have driven the evolutionary history of these genes.

An enigmatic pelagic fish with internalized red muscle: A future regional endotherm or forever an ectotherm?

Ectothermy and endothermy in extant fishes are defined by distinct integrated suites of characters. Although only ⁓0.1% of fishes are known to have endothermic capacity, recent discoveries suggest that there may still be uncommon pelagic fish species with yet to be discovered endothermic traits. Among the most rarely encountered marine fishes, the louvar Luvarus imperialis is a remarkable example of adaptive evolution as the only extant pelagic species in the order Acanthuriformes (including surgeonfishes, tangs, unicornfishes and Moorish idol). Magnetic resonance imaging and gross necropsy did not yield evidence of cranial or visceral endothermy but revealed a central-posterior distribution of myotomal red muscle that is a mixture of the character states typifying ectotherms (lateral-posterior) and red muscle endotherms (central-anterior). Dissection of a specimen confirmed, and an osteological proxy supported, that L. imperialis has not evolved the vascular rete that is vital to retaining heat in the red muscle. The combination of presumably relying on caudal propulsion while exhibiting internal red muscle without associated retia is unique to L. imperialis among all extant fishes, raising the macroevolutionary question of whether this species - in geologic timescales - will remain an ectotherm or evolve red muscle endothermy.

Range-wide population genomics of common seadragons shows secondary contact over a former barrier and insights on illegal capture

BACKGROUND

Common seadragons (Phyllopteryx taeniolatus, Syngnathidae) are an emblem of the diverse endemic fauna of Australia's southern rocky reefs, the newly recognized "Great Southern Reef." A lack of assessments spanning this global biodiversity hotspot in its entirety is currently hampering an understanding of the factors that have contributed to its diversity. The common seadragon has a wide range across Australia's entire temperate south and includes a geogenetic break over a former land bridge, which has called its status as a single species into question. As a popular aquarium display that sells for high prices, common seadragons are also vulnerable to illegal capture.

RESULTS

Here, we provide range-wide nuclear sequences (986 variable Ultraconserved Elements) for 198 individuals and mitochondrial genomes for 140 individuals to assess species status, identify genetic units and their diversity, and trace the source of two poached individuals. Using published data of the other two seadragon species, we found that lineages of common seadragons have diverged relatively recently (< 0.63 Ma). Within common seadragons, we found pronounced genetic structure, falling into three major groups in the western, central, and eastern parts of the range. While populations across the Bassian Isthmus were divergent, there is also evidence for secondary contact since the passage opened. We found a strong cline of genetic diversity from the range center tapering symmetrically towards the range peripheries. Based on their genetic similarities, the poached individuals were inferred to have originated from around Albany in southwestern Australia.

CONCLUSIONS

We conclude that common seadragons constitute a single species with strong geographic structure but coherence through gene flow. The low genetic diversity on the east and west coasts is concerning given that these areas are projected to face fast climate change. Our results suggest that in addition to their life history, geological events and demographic expansions have all played a role in shaping populations in the temperate south. These insights are an important step towards understanding the historical determinants of the diversity of species endemic to the Great Southern Reef.

Phylogenetic and Evolutionary Comparison of Mitogenomes Reveal Adaptive Radiation of Lampriform Fishes

Lampriform fishes (Lampriformes), which primarily inhabit deep-sea environments, are large marine fishes varying from the whole-body endothermic opah to the world's longest bony fish-giant oarfish, with species morphologies varying from long and thin to deep and compressed, making them an ideal model for studying the adaptive radiation of teleost fishes. Moreover, this group is important from a phylogenetic perspective owing to their ancient origins among teleosts. However, knowledge about the group is limited, which is, at least partially, due to the dearth of recorded molecular data. This study is the first to analyze the mitochondrial genomes of three lampriform species (Lampris incognitus, Trachipterus ishikawae, and Regalecus russelii) and infer a time-calibrated phylogeny, including 68 species among 29 orders. Our phylomitogenomic analyses support the classification of Lampriformes as monophyletic and sister to Acanthopterygii; hence, addressing the longstanding controversy regarding the phylogenetic status of Lampriformes among teleosts. Comparative mitogenomic analyses indicate that tRNA losses existed in at least five Lampriformes species, which may reveal the mitogenomic structure variation associated with adaptive radiation. However, codon usage in Lampriformes did not change significantly, and it is hypothesized that the nucleus transported the corresponding tRNA, which led to function substitutions. The positive selection analysis revealed that atp8 and cox3 were positively selected in opah, which might have co-evolved with the endothermic trait. This study provides important insights into the systematic taxonomy and adaptive evolution studies of Lampriformes species.

Syngnathoid Evolutionary History and the Conundrum of Fossil Misplacement

Seahorses, pipefishes, trumpetfishes, shrimpfishes, and allies are a speciose, globally distributed clade of fishes that have evolved a large number of unusual body plans. The clade that includes all these forms, Syngnathoidei, has become a model for the study of life history evolution, population biology, and biogeography. Yet, the timeline of syngnathoid evolution has remained highly contentious. This debate is largely attributable to the nature of the syngnathoid fossil record, which is both poorly described and patchy for several major lineages. Although fossil syngnathoids have been used to calibrate molecular phylogenies, the interrelationships of extinct species and their affinities to major living syngnathoid clades have scarcely been quantitatively tested. Here, I use an expanded morphological dataset to reconstruct the evolutionary relationships and clade ages of fossil and extant syngnathoids. Phylogenies generated using different analytical methodologies are largely congruent with molecular phylogenetic trees of Syngnathoidei but consistently find novel placements for several key taxa used as fossil calibrators in phylogenomic studies. Tip-dating of the syngnathoid phylogeny finds a timeline for their evolution that differs slightly from the one inferred using molecular trees but is generally congruent with a post-Cretaceous diversification event. These results emphasize the importance of quantitatively testing the relationships of fossil species, particularly when they are critical to assessing divergence times.

Evolution of the germline mutation rate across vertebrates

The germline mutation rate determines the pace of genome evolution and is an evolving parameter itself¹. However, little is known about what determines its evolution, as most studies of mutation rates have focused on single species with different methodologies². Here we quantify germline mutation rates across vertebrates by sequencing and comparing the high-coverage genomes of 151 parent-offspring trios from 68 species of mammals, fishes, birds and reptiles. We show that the per-generation mutation rate varies among species by a factor of 40, with mutation rates being higher for males than for females in mammals and birds, but not in reptiles and fishes. The generation time, age at maturity and species-level fecundity are the key life-history traits affecting this variation among species. Furthermore, species with higher long-term effective population sizes tend to have lower mutation rates per generation, providing support for the drift barrier hypothesis³. The exceptionally high yearly mutation rates of domesticated animals, which have been continually selected on fecundity traits including shorter generation times, further support the importance of generation time in the evolution of mutation rates. Overall, our comparative analysis of pedigree-based mutation rates provides ecological insights on the mutation rate evolution in vertebrates.

A latitudinal gradient of deep-sea invasions for marine fishes

Although the tropics harbor the greatest species richness globally, recent work has demonstrated that, for many taxa, speciation rates are faster at higher latitudes. Here, we explore lability in oceanic depth as a potential mechanism for this pattern in the most biodiverse vertebrates - fishes. We demonstrate that clades with the highest speciation rates also diversify more rapidly along the depth gradient, drawing a fundamental link between evolutionary and ecological processes on a global scale. Crucially, these same clades also inhabit higher latitudes, creating a prevailing latitudinal gradient of deep-sea invasions concentrated in poleward regions. We interpret these findings in the light of classic ecological theory, unifying the latitudinal variation of oceanic features and the physiological tolerances of the species living there. This work advances the understanding of how niche lability sculpts global patterns of species distributions and underscores the vulnerability of polar ecosystems to changing environmental conditions.

Whole-genome comparisons identify repeated regulatory changes underlying convergent appendage evolution in diverse fish lineages

Fins are major functional appendages of fish that have been repeatedly modified in different lineages. To search for genomic changes underlying natural fin diversity, we compared the genomes of 36 wild fish species that either have complete or reduced pelvic and caudal fins. We identify 1,614 genomic regions that are well-conserved in fin-complete species but missing from multiple fin-reduced lineages. Recurrent deletions of conserved sequences (CONDELs) in wild fin-reduced species are enriched for functions related to appendage development, suggesting that convergent fin reduction at the organismal level is associated with repeated genomic deletions near fin-appendage development genes. We used sequencing and functional enhancer assays to confirm that PelA , a Pitx1 enhancer previously linked to recurrent pelvic loss in sticklebacks, has also been independently deleted and may have contributed to the fin morphology in distantly related pelvic-reduced species. We also identify a novel enhancer that is conserved in the majority of percomorphs, drives caudal fin expression in transgenic stickleback, is missing in tetraodontiform, s yngnathid, and synbranchid species with caudal fin reduction, and which alters caudal fin development when targeted by genome editing. Our study illustrates a general strategy for mapping phenotypes to genotypes across a tree of vertebrate species, and highlights notable new examples of regulatory genomic hotspots that have been used to evolve recurrent phenotypes during 100 million years of fish evolution.

Phylogenetics and the Cenozoic radiation of lampreys

The development of a movable jaw is one of the most important transitions in the evolutionary history of animals.¹ Jawed vertebrates rapidly diversified after appearing approximately 470 million years ago. Today, only lampreys and hagfishes represent the once dominant jawless grade^2,3,4 and comprise less than 1% of living vertebrate species. Their relationship to other vertebrates ranks among the more contentious problems in animal phylogenetics.^{5,6,7,8,9,10,11,12} Further, the phylogenetic relationships within lampreys and hagfishes remain unclear,^13,14,15 and the ages of their living lineages are largely unexplored.^16,17 Because of their importance for the genomic and developmental changes that prefigured jawed vertebrate diversity,^18,19,20,21 the evolutionary history of lampreys and hagfishes is a major frontier of organismal biology. Of these two clades, lampreys²² are more ecologically diverse, exhibiting freshwater, anadromous, and fully marine forms, as well as parasitic and nonparasitic species.^23,24 Here, we present a new phylogeny and historical biogeographic reconstruction of all living lampreys. Whereas the early diversification of this clade tracks Pangaean fragmentation, lampreys also rapidly radiated in the northern hemisphere during the mid-Cretaceous and directly after the Cretaceous-Paleogene extinction. These radiations mirrored concurrent ones in other animals and plants and coincided with changes to lamprey ecology and feeding behavior. Our results suggest that 80% of living lamprey clades appeared in the last 20 million years of Earth history. Rather than gradually accumulating since the oldest stem-group forms appeared in the early Paleozoic, living lamprey biodiversity results from diversifications extending from the Cretaceous to present.

Evolutionary time and species diversity in aquatic ecosystems worldwide

The latitudinal diversity gradient (LDG) is frequently described as the most dramatic biodiversity pattern on Earth, yet ecologists and biogeographers have failed to reach consensus on its primary cause. A key problem in explaining the LDG involves collinearity between multiple factors that are predicted to affect species richness in the same direction. In terrestrial systems, energy input, geographic area, and evolutionary time for species accumulation tend to covary positively with species richness at the largest spatial scales, such that their individual contributions to the LDG are confounded in global analyses. I review three diversity patterns from marine and freshwater systems that break this collinearity and which may thus provide stronger tests of the influence of time on global richness gradients. Specifically, I contrast biodiversity patterns along oceanic depth gradients, in geologically young versus ancient lakes, and in the north versus south polar marine biomes. I focus primarily on fishes due to greater data availability but synthesize patterns for invertebrates where possible. I find that regional-to-global species richness generally declines with depth in the oceans, despite the great age and stability of the deep-sea biome. Geologically ancient lakes generally do not contain more species than young lakes, and the Antarctic marine biome is not appreciably more species rich than the much younger Arctic marine biome. However, endemism is consistently higher in older systems. Patterns for invertebrate groups are less clear than for fishes and reflect a critical need for primary biodiversity data. In summary, the available data suggest that species richness is either decoupled from or only weakly related to the amount of time for diversification. These results suggest that energy, productivity, or geographic area are the primary drivers of large-scale diversity gradients. To the extent that marine and terrestrial diversity gradients result from similar processes, these examples provide evidence against a primary role for evolutionary time as the cause of the LDG.

Alternating regimes of shallow and deep-sea diversification explain a species-richness paradox in marine fishes

The deep sea contains a surprising diversity of life, including iconic fish groups such as anglerfishes and lanternfishes. Still, >65% of marine teleost fish species are restricted to the photic zone <200 m, which comprises less than 10% of the ocean's total volume. From a macroevolutionary perspective, this paradox may be explained by three hypotheses: 1) shallow water lineages have had more time to diversify than deep-sea lineages, 2) shallow water lineages have faster rates of speciation than deep-sea lineages, or 3) shallow-to-deep sea transition rates limit deep-sea richness. Here we use phylogenetic comparative methods to test among these three non-mutually exclusive hypotheses. While we found support for all hypotheses, the disparity in species richness is better described as the uneven outcome of alternating phases that favored shallow or deep diversification over the past 200 million y. Shallow marine teleosts became incredibly diverse 100 million y ago during a period of warm temperatures and high sea level, suggesting the importance of reefs and epicontinental settings. Conversely, deep-sea colonization and speciation was favored during brief episodes when cooling temperatures increased the efficiency of the ocean's carbon pump. Finally, time-variable ecological filters limited shallow-to-deep colonization for much of teleost history, which helped maintain higher shallow richness. A pelagic lifestyle and large jaws were associated with early deep-sea colonists, while a demersal lifestyle and a tapered body plan were typical of later colonists. Therefore, we also suggest that some hallmark characteristics of deep-sea fishes evolved prior to colonizing the deep sea.

Phylogenomics and body shape morphometrics reveal recent diversification in the goatfishes (Syngnatharia: Mullidae)

Clades of marine fishes exhibit many patterns of diversification, ranging from relatively constant throughout time to rapid changes in the rates of speciation and extinction. The goatfishes (Syngnatharia: Mullidae) are a family of marine, reef associated fishes with a relatively recent origin, distributed globally in tropical and temperate waters. Despite their abundance and economic importance, the goatfishes remain one of the few coral reef families for which the species level relationships have not been examined using genomic techniques. Here we use phylogenomic analysis of ultra-conserved elements (UCE) and exon data to resolve a well-supported, time-calibrated phylogeny for 72 species of goatfishes, supporting a recent crown age of the goatfishes at 21.9 million years ago. We used this framework to test hypotheses about the associations among body shape morphometrics, taxonomy, and phylogeny, as well as to explore relative diversification rates across the phylogeny. Body shape was strongly associated with generic-level taxonomy of goatfishes, with morphometric analyses showing evidence for high phylogenetic signal across all morphotypes. Rates of diversification in this clade reveal a recent sharp increase in lineage accumulation, with 92% of the goatfish species sampled across all clades and major body plans having originated in just the past 5 million years. We suggest that habitat diversity in the early Pliocene oceans and the generalist ecology of goatfishes are key factors in the unusual evolutionary tempo of the family Mullidae.

Paleogene emergence and evolutionary history of the Amazonian fossorial fish genus Tarumania (Teleostei: Tarumaniidae)

Tarumania walkerae is a rare fossorial freshwater fish species from the lower Rio Negro, Central Amazonia, composing the monotypic and recently described family Tarumaniidae. The family has been proposed as the sister group of Erythrinidae by both morphological and molecular studies despite distinct arrangements of the superfamily Erythrinoidea within Characiformes. Recent phylogenomic studies and time-calibrated analyses of characoid fishes have not included specimens of Tarumania in their analyses. We obtained genomic data for T. walkerae and constructed a phylogeny based on 1795 nuclear loci with 488,434 characters of ultraconserved elements (UCEs) for 108 terminals including specimens of all 22 characiform families. The phylogeny confirms the placement of Tarumaniidae as sister to Erythrinidae but differs from the morphological hypothesis in the placement of the two latter families as sister to the clade with Hemiodontidae, Cynodontidae, Serrasalmidae, Parodontidae, Anostomidae, Prochilodontidae, Chilodontidae, and Curimatidae. The phylogeny calibrated with five characoid fossils indicates that Erythrinoidea diverged from their relatives during the Late Cretaceous circa 90 Ma (108–72 Ma), and that Tarumania diverged from the most recent common ancestor of Erythrinidae during the Paleogene circa 48 Ma (66–32 Ma). The occurrence of the erythrinoid-like †Tiupampichthys in the Late Cretaceous–Paleogene formations of the El Molino Basin of Bolivia supports our hypothesis for the emergence of the modern Erythrinidae and Tarumaniidae during the Paleogene.

The rise of biting during the Cenozoic fueled reef fish body shape diversification

We demonstrate that the stunning trophic diversity of modern reef fishes is a relatively recent state driven by a dramatic transformation in representation of major feeding modes. Since the Early Cenozoic, when over 95% of teleost lineages were suction feeders, there has been a steady increase in direct biting feeding modes. A variety of novelties and jaw modifications permitted reef fishes to feed on substrate-bound prey using direct biting and grazing behaviors and opened this rich adaptive zone, which we show elevated rates of body shape evolution. Taken together, our results indicate that recent diversification of the feeding mechanism played a major role in ecologically and phenotypically shaping the modern fauna of reef fishes.

Diversity of feeding mechanisms is a hallmark of reef fishes, but the history of this variation is not fully understood. Here, we explore the emergence and proliferation of a biting mode of feeding, which enables fishes to feed on attached benthic prey. We find that feeding modes other than suction, including biting, ram biting, and an intermediate group that uses both biting and suction, were nearly absent among the lineages of teleost fishes inhabiting reefs prior to the end-Cretaceous mass extinction, but benthic biting has rapidly increased in frequency since then, accounting for about 40% of reef species today. Further, we measured the impact of feeding mode on body shape diversification in reef fishes. We fit a model of multivariate character evolution to a dataset comprising three-dimensional body shape of 1,530 species of teleost reef fishes across 111 families. Dedicated biters have accumulated over half of the body shape variation that suction feeders have in just 18% of the evolutionary time by evolving body shape ∼1.7 times faster than suction feeders. As a possible response to the ecological and functional diversity of attached prey, biters have dynamically evolved both into shapes that resemble suction feeders as well as novel body forms characterized by lateral compression and small jaws. The ascendance of species that use biting mechanisms to feed on attached prey reshaped modern reef fish assemblages and has been a major contributor to their ecological and phenotypic diversification.

FishPIE: a universal phylogenetically informative exon markers set for ray-finned fishes

Understanding the evolutionary history of the highly diverse ray-finned fishes has been challenging, and the development of more universal primers for phylogenetic analyses may help overcoming these challenges. We developed FishPIE, a nested PCR primer set of 82 phylogenetically informative exon markers, and tested it on 203 species from 31 orders of Actinopterygii. We combined orthologous sequences of the FishPIE markers obtained from published genomes and transcriptomes and constructed the phylogeny of 710 species belonging to 190 families and 60 orders. The resulting phylogenies had topologies comparable to previous phylogenomic studies. We demonstrated that the FishPIE markers could address phylogenetic questions across broad taxonomic levels. By incorporating the newly sequenced taxa, we were able to shed new light on the phylogeny of the highly diverse Cypriniformes. Thus, FishPIE holds great promise for generating genetic data for broad taxonomic groups and accelerating our understanding of the fish tree of life.

•

FishPIE is a nested PCR primer set of 82 markers for fish phylogenetic analysis

•

The markers can be broadly applied to all orders of ray-finned fishes

•

Their phylogenetic performance is comparable to that of genomic analyses

Rising from the Ashes: The Biogeographic Origins of Modern Coral Reef Fishes

During the excavation of Mayan tombs, little did the archaeologists know that the fossils they discovered in the tomb stones would fundamentally alter our understanding of the earliest origins of coral reef fishes. Located just 500 kilometers from the point where an asteroid impact reconfigured the world's biological systems 66 million years ago, we find the earliest origins of three typical reef fish groups. Their presence in Mexico just 3 million years after this impact finally reconciles the conflict between the fossil and phylogenetic evidence for the earliest origins of reef fishes. The incorporation of these fossils into a global reconstruction of fish evolutionary history reveals a new picture of the early biogeography of reef fishes, with strong Atlantic links. From locations associated with biological destruction and societal collapse, we see evidence of the origins of one of the world's most diverse and spectacular marine ecosystems: coral reefs.

The impact of paleoclimatic changes on body size evolution in marine fishes

General rules are useful tools for understanding how organisms evolve. Cope’s rule (tendency to increase in size over evolutionary time) and Bergmann’s rule (tendency to grow to larger sizes in cooler climates) both relate to body size, an important factor that affects the biology, ecology, and physiology of organisms. These rules are well studied in endotherms but remain poorly understood among ectotherms. Here, we show that paleoclimatic changes strongly shaped the trajectory of body size evolution in tetraodontiform fishes. Their body size evolution is explained by both Cope’s and Bergmann’s rules, highlighting the impact of paleoclimatic changes on aquatic organisms, which rely on their environment for temperature regulation and are likely more susceptible than terrestrial vertebrates to climatic changes.

Body size is an important species trait, correlating with life span, fecundity, and other ecological factors. Over Earth’s geological history, climate shifts have occurred, potentially shaping body size evolution in many clades. General rules attempting to summarize body size evolution include Bergmann’s rule, which states that species reach larger sizes in cooler environments and smaller sizes in warmer environments, and Cope’s rule, which poses that lineages tend to increase in size over evolutionary time. Tetraodontiform fishes (including pufferfishes, boxfishes, and ocean sunfishes) provide an extraordinary clade to test these rules in ectotherms owing to their exemplary fossil record and the great disparity in body size observed among extant and fossil species. We examined Bergmann’s and Cope’s rules in this group by combining phylogenomic data (1,103 exon loci from 185 extant species) with 210 anatomical characters coded from both fossil and extant species. We aggregated data layers on paleoclimate and body size from the species examined, and inferred a set of time-calibrated phylogenies using tip-dating approaches for downstream comparative analyses of body size evolution by implementing models that incorporate paleoclimatic information. We found strong support for a temperature-driven model in which increasing body size over time is correlated with decreasing oceanic temperatures. On average, extant tetraodontiforms are two to three times larger than their fossil counterparts, which otherwise evolved during periods of warmer ocean temperatures. These results provide strong support for both Bergmann’s and Cope’s rules, trends that are less studied in marine fishes compared to terrestrial vertebrates and marine invertebrates.

Molecular phylogenetics of the Clupeiformes based on exon-capture data and a new classification of the order

The Clupeiformes, including among others herrings, anchovies, shads and menhadens are ecologically and commercially important, yet their phylogenetic relationships are still controversial. Previous classification of Clupeiformes were based on morphological characters or lack of synapomorphic characters. More recent studies based on molecular data as well as new morphological evidence are keeping challenging their phylogenetic relations and there is still no consensus on many interrelationships within the Clupeiformes. In this study, we collected nuclear sequence data from 4,434 single-copy protein coding loci using a gene-capture method. We obtained a robust phylogeny based on 1,165 filtered loci with less than 30 % missing data. Our major findings include: 1) reconfirmation of monophyly of the Clupeiformes, that is, Denticipitidae is sister to all other clupeiforms; 2) the polyphyletic nature of dussumieriids and early branching of Spratelloididae from all other clupeoids were confirmed using datasets curated for less missing data and more balanced base composition in the respective taxa. The next branching clade is the monophyletic Engraulidae. Pristigasteridae also is monophyletic, but it was nested in the previously defined "Clupeidae". Within Pristigasteridae there is no support for monophyletic Pelloninae. Chirocentrus is close to Dussumieria and not to engraulids. The miniaturized Sundasalanx is placed close to the ehiravine Clupeonella, however, with a relatively deep split. The genus Clupea, is not part of the diverse "Clupeidae", but part of a clade containing additionally Sprattus and Etrumeus. Within the crown group clades, Alosidae and Dorosomatidae are retrieved as sister clades. Based on new fossil calibration points, we found that major lineages of the clupeiforms diverged in the late Cretaceous and early Paleogene. The extinction event at the end of the Cretaceous may have created ecological niches, which could have fueled the diversification of clupeiform fishes. Based on the strong evidence of the present study, we propose an updated classification of Clupeiformes consisting of ten families: Denticipitidae; Spratelloididae; Engraulidae (Engraulinae + Coiliinae); Clupeidae; Chirocentridae; Dussumieriidae; Pristigasteridae; Ehiravidae; Alosidae, Dorosomatidae.

Molecular phylogenetics reveals the evolutionary history of marine fishes (Actinopterygii) endemic to the subtropical islands of the Southwest Pacific

Remote oceanic islands of the Pacific host elevated levels of actinopterygian (ray-finned fishes) endemism. Characterizing the evolutionary histories of these endemics has provided insight into the generation and maintenance of marine biodiversity in many regions. The subtropical islands of Lord Howe, Norfolk, and Rangitāhua (Kermadec) in the Southwest Pacific are yet to be comprehensively studied. Here, we characterize the spatio-temporal diversification of marine fishes endemic to these Southwest Pacific islands by combining molecular phylogenies and the geographic distribution of species. We built Bayesian ultrametric trees based on open-access and newly generated sequences for five mitochondrial and ten nuclear loci, and using fossil data for time calibration. We present the most comprehensive phylogenies to date for marine ray-finned fish genera, comprising 34 species endemic to the islands, including the first phylogenetic placements for 11 endemics. Overall, our topologies confirm the species status of all endemics, including three undescribed taxa. Our phylogenies highlight the predominant affinity of these endemics with the Australian fish fauna (53%), followed by the East Pacific (15%), and individual cases where the closest sister taxon of our endemic is found in the Northwest Pacific and wider Indo-Pacific. Nonetheless, for a quarter of our focal endemics, their geographic affinity remains unresolved due to sampling gaps within their genera. Our divergence time estimates reveal that the majority of endemic lineages (67.6%) diverged after the emergence of Lord Howe (6.92 Ma), the oldest subtropical island in the Southwest Pacific, suggesting that these islands have promoted diversification. However, divergence ages of some endemics pre-date the emergence of the islands, suggesting they may have originated outside of these islands, or, in some cases, ages may be overestimated due to unsampled taxa. To fully understand the role of the Southwest Pacific subtropical islands as a 'cradle' for diversification, our study advocates for further regional surveys focused on tissue collection for DNA analysis.

Prolonged morphological expansion of spiny-rayed fishes following the end-Cretaceous

Spiny-rayed fishes (Acanthomorpha) dominate modern marine habitats and account for more than a quarter of all living vertebrate species. Previous time-calibrated phylogenies and patterns from the fossil record explain this dominance by correlating the origin of major acanthomorph lineages with the Cretaceous-Palaeogene mass extinction. Here we infer a time-calibrated phylogeny using ultraconserved elements that samples 91.4% of all acanthomorph families and investigate patterns of body shape disparity. Our results show that acanthomorph lineages steadily accumulated throughout the Cenozoic and underwent a significant expansion of among-clade morphological disparity several million years after the end-Cretaceous. These acanthomorph lineages radiated into and diversified within distinct regions of morphospace that characterize iconic lineages, including fast-swimming open-ocean predators, laterally compressed reef fishes, bottom-dwelling flatfishes, seahorses and pufferfishes. The evolutionary success of spiny-rayed fishes is the culmination of multiple species-rich and phenotypically disparate lineages independently diversifying across the globe under a wide range of ecological conditions.

Phylogenomic analysis of Syngnathidae reveals novel relationships, origins of endemic diversity and variable diversification rates

Background

Seahorses, seadragons, pygmy pipehorses, and pipefishes (Syngnathidae, Syngnathiformes) are among the most recognizable groups of fishes because of their derived morphology, unusual life history, and worldwide distribution. Despite previous phylogenetic studies and recent new species descriptions of syngnathids, the evolutionary relationships among several major groups within this family remain unresolved.

Results

Here, we provide a reconstruction of syngnathid phylogeny based on genome-wide sampling of 1314 ultraconserved elements (UCEs) and expanded taxon sampling to assess the current taxonomy and as a basis for macroevolutionary insights. We sequenced a total of 244 new specimens across 117 species and combined with published UCE data for a total of 183 species of Syngnathidae, about 62% of the described species diversity, to compile the most data-rich phylogeny to date. We estimated divergence times using 14 syngnathiform fossils, including nine fossils with newly proposed phylogenetic affinities, to better characterize current and historical biogeographical patterns, and to reconstruct diversification through time. We present a phylogenetic hypothesis that is well-supported and provides several notable insights into syngnathid evolution. We found nine non-monophyletic genera, evidence for seven cryptic species, five potentially invalid synonyms, and identified a novel sister group to the seahorses, the Indo-Pacific pipefishes Halicampus macrorhynchus and H. punctatus. In addition, the morphologically distinct southwest Pacific seahorse Hippocampus jugumus was recovered as the sister to all other non-pygmy seahorses. As found in many other groups, a high proportion of syngnathid lineages appear to have originated in the Central Indo-Pacific and subsequently dispersed to adjoining regions. Conversely, we also found an unusually high subsequent return of lineages from southern Australasia to the Central Indo-Pacific. Diversification rates rose abruptly during the Middle Miocene Climate Transition and peaked after the closure of the Tethys Sea.

Conclusions

Our results reveal a previously underappreciated diversity of syngnathid lineages. The observed biogeographic patterns suggest a significant role of the southern Australasian region as a source and sink of lineages. Shifts in diversification rates imply possible links to declining global temperatures, the separation of the Atlantic and Pacific faunas, and the environmental changes associated with these events.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01271-w.

Phylogenomic analyses highlight innovation and introgression in the continental radiations of Fagaceae across the Northern Hemisphere

Northern Hemisphere forests changed drastically in the early Eocene with the diversification of the oak family (Fagaceae). Cooling climates over the next 20 million years fostered the spread of temperate biomes that became increasingly dominated by oaks and their chestnut relatives. Here we use phylogenomic analyses of nuclear and plastid genomes to investigate the timing and pattern of major macroevolutionary events and ancient genome-wide signatures of hybridization across Fagaceae. Innovation related to seed dispersal is implicated in triggering waves of continental radiations beginning with the rapid diversification of major lineages and resulting in unparalleled transformation of forest dynamics within 15 million years following the K-Pg extinction. We detect introgression at multiple time scales, including ancient events predating the origination of genus-level diversity. As oak lineages moved into newly available temperate habitats in the early Miocene, secondary contact between previously isolated species occurred. This resulted in adaptive introgression, which may have further amplified the diversification of white oaks across Eurasia.

Fagaceae are diverse family including trees of ecological and economic importance. This phylogenomic analysis of nuclear and plastid genomes reconstructs evolutionary history and finds evidence of multiple adaptive introgression events in this important plant family.

Phylogeny, diversification, and biogeography of a hemiclonal hybrid system of native Australian freshwater fishes (Gobiiformes: Gobioidei: Eleotridae: Hypseleotris)

Background

Carp gudgeons (genus Hypseleotris) are a prominent part of the Australian freshwater fish fauna, with species distributed around the western, northern, and eastern reaches of the continent. We infer a calibrated phylogeny of the genus based on nuclear ultraconserved element (UCE) sequences and using Bayesian estimation of divergence times, and use this phylogeny to investigate geographic patterns of diversification with GeoSSE. The southeastern species have hybridized to form hemiclonal lineages, and we also resolve relationships of hemiclones and compare their phylogenetic placement in the UCE phylogeny with a hypothesis based on complete mitochondrial genomes. We then use phased SNPs extracted from the UCE sequences for population structure analysis among the southeastern species and hemiclones.

Results

Hypseleotris cyprinoides, a widespread euryhaline species known from throughout the Indo-Pacific, is resolved outside the remainder of the species. Two Australian radiations comprise the bulk of Hypseleotris, one primarily in the northwestern coastal rivers and a second inhabiting the southeastern region including the Murray–Darling, Bulloo-Bancannia and Lake Eyre basins, plus coastal rivers east of the Great Dividing Range. Our phylogenetic results reveal cytonuclear discordance between the UCE and mitochondrial hypotheses, place hemiclone hybrids among their parental taxa, and indicate that the genus Kimberleyeleotris is nested within the northwestern Hypseleotris radiation along with three undescribed species. We infer a crown age for Hypseleotris of 17.3 Ma, date the radiation of Australian species at roughly 10.1 Ma, and recover the crown ages of the northwestern (excluding H. compressa) and southeastern radiations at 5.9 and 7.2 Ma, respectively. Range-dependent diversification analyses using GeoSSE indicate that speciation and extinction rates have been steady between the northwestern and southeastern Australian radiations and between smaller radiations of species in the Kimberley region and the Arnhem Plateau. Analysis of phased SNPs confirms inheritance patterns and reveals high levels of heterozygosity among the hemiclones.

Conclusions

The northwestern species have restricted ranges and likely speciated in allopatry, while the southeastern species are known from much larger areas, consistent with peripatric speciation or allopatric speciation followed by secondary contact. Species in the northwestern Kimberley region differ in shape from those in the southeast, with the Kimberley species notably more elongate and slender than the stocky southeastern species, likely due to the different topographies and flow regimes of the rivers they inhabit.

Evolutionary transition to the ectomycorrhizal habit in the genomes of a hyperdiverse lineage of mushroom-forming fungi

The ectomycorrhizal (ECM) symbiosis has independently evolved from diverse types of saprotrophic ancestors. In this study, we seek to identify genomic signatures of the transition to the ECM habit within the hyperdiverse Russulaceae. We present comparative analyses of the genomic architecture and the total and secreted gene repertoires of 18 species across the order Russulales, of which 13 are newly sequenced, including a representative of a saprotrophic member of Russulaceae, Gloeopeniophorella convolvens. The genomes of ECM Russulaceae are characterized by a loss of genes for plant cell wall-degrading enzymes (PCWDEs), an expansion of genome size through increased transposable element (TE) content, a reduction in secondary metabolism clusters, and an association of small secreted proteins (SSPs) with TE 'nests', or dense aggregations of TEs. Some PCWDEs have been retained or even expanded, mostly in a species-specific manner. The genome of G. convolvens possesses some characteristics of ECM genomes (e.g. loss of some PCWDEs, TE expansion, reduction in secondary metabolism clusters). Functional specialization in ECM decomposition may drive diversification. Accelerated gene evolution predates the evolution of the ECM habit, indicating that changes in genome architecture and gene content may be necessary to prime the evolutionary switch.

Molecular phylogeny of the threadfin fishes (Polynemidae) using ultraconserved elements

Threadfins (Teleostei: Polynemidae) are a group of fishes named for their elongated and threadlike pectoral-fin rays. These fishes are commonly found in the world's tropical and subtropical waters, and are an economically important group for people living in these regions, with more than 100,000 t harvested in recent years. However, we do not have a detailed understanding of polynemid evolutionary history such that these fishes can be monitored, managed and conserved as an important tropical food source. Recent studies hypothesize at least one genus of threadfins is polyphyletic, and no studies have focused on generating a hypothesis of relationship for the Polynemidae using DNA sequences. In this study, we analyse a genomic dataset of ultraconserved-element and mitochondrial loci to construct a phylogeny of the Polynemidae. We recover the threadfins as a clade sister to flatfishes, with the most taxonomically rich genus, Polydactylus, being resolved as polyphyletic. When comparing our dataset to data from previous studies, we find that a few recent broad-scale phylogenies of fishes have incorporated mislabelled, misidentified or chimeric terminals into their analyses, impacting the relationships of threadfins they recover. We highlight these problematic sequences, providing revised identifications based on the data sequenced in this study. We then discuss the intrarelationships of threadfins, highlighting morphological or ecological characters that support the clades we recover.

Limitations of phylogenomic data can drive inferred speciation rate shifts

Biodiversity analyses of phylogenomic timetrees have produced many high-profile examples of shifts in the rate of speciation across the tree of life. Temporally correlated events in ecology, climate, and biogeography are frequently invoked to explain these rate shifts. In a re-examination of 15 genomic timetrees and 25 major published studies of the pattern of speciation through time, we observed an unexpected correlation between the timing of reported rate shifts and the information content of sequence alignments. Here, we show that the paucity of sequence variation and insufficient species sampling in phylogenomic datasets are the likely drivers of many inferred speciation rate shifts, rather than the proposed biological explanations. Therefore, data limitations can produce predictable but spurious signals of rate shifts even when speciation rates may be similar across taxa and time. Our results suggest that the reliable detection of speciation rate shifts requires the acquisition and assembly of long phylogenomic alignments with near-complete species sampling and accurate estimates of species richness for the clades of study.

Investigations into the ancestry of the Grape-eye Seabass (Hemilutjanus macrophthalmos) reveal novel limits and relationships for the Acropomatiformes (Teleostei: Percomorpha)

Abstract For 175 years, an unremarkable bass, the Grape-eye Seabass (Hemilutjanus macrophthalmos), has been known from coastal waters in the Eastern Pacific. To date, its phylogenetic placement and classification have been ignored. A preliminary osteological examination of Hemilutjanus hinted that it may have affinities with the Acropomatiformes. To test this hypothesis, we conducted a phylogenetic analysis using UCE and Sanger sequence data to study the placement of Hemilutjanus and the limits and relationships of the Acropomatiformes. We show that Hemilutjanus is a malakichthyid, and our results corroborate earlier studies that have resolved a polyphyletic Polyprionidae; accordingly, we describe Stereolepididae, new family, for Stereolepis. With these revisions, the Acropomatiformes is now composed of the: Acropomatidae; Banjosidae; Bathyclupeidae; Champsodontidae; Creediidae; Dinolestidae; Epigonidae; Glaucosomatidae; Hemerocoetidae; Howellidae; Lateolabracidae; Malakichthyidae; Ostracoberycidae; Pempheridae; Pentacerotidae; Polyprionidae; Scombropidae; Stereolepididae, new family; Symphysanodontidae; Synagropidae; and Schuettea. Finally, using our new hypothesis, we demonstrate that acropomatiforms repeatedly evolved bioluminescence and transitioned between shallow waters and the deep sea.

A highly contiguous nuclear genome assembly of the mandarinfish Synchiropus splendidus (Syngnathiformes: Callionymidae)

The fish order Syngnathiformes has been referred to as a collection of misfit fishes, comprising commercially important fish such as red mullets as well as the highly diverse seahorses, pipefishes, and seadragons-the well-known family Syngnathidae, with their unique adaptations including male pregnancy. Another ornate member of this order is the species mandarinfish. No less than two types of chromatophores have been discovered in the spectacularly colored mandarinfish: the cyanophore (producing blue color) and the dichromatic cyano-erythrophore (producing blue and red). The phylogenetic position of mandarinfish in Syngnathiformes, and their promise of additional genetic discoveries beyond the chromatophores, made mandarinfish an appealing target for whole-genome sequencing. We used linked sequences to create synthetic long reads, producing a highly contiguous genome assembly for the mandarinfish. The genome assembly comprises 483 Mbp (longest scaffold 29 Mbp), has an N50 of 12 Mbp, and an L50 of 14 scaffolds. The assembly completeness is also high, with 92.6% complete, 4.4% fragmented, and 2.9% missing out of 4584 BUSCO genes found in ray-finned fishes. Outside the family Syngnathidae, the mandarinfish represents one of the most contiguous syngnathiform genome assemblies to date. The mandarinfish genomic resource will likely serve as a high-quality outgroup to syngnathid fish, and furthermore for research on the genomic underpinnings of the evolution of novel pigmentation.

Divergence time estimation of Galliformes based on the best gene shopping scheme of ultraconserved elements

BACKGROUND

Divergence time estimation is fundamental to understanding many aspects of the evolution of organisms, such as character evolution, diversification, and biogeography. With the development of sequence technology, improved analytical methods, and knowledge of fossils for calibration, it is possible to obtain robust molecular dating results. However, while phylogenomic datasets show great promise in phylogenetic estimation, the best ways to leverage the large amounts of data for divergence time estimation has not been well explored. A potential solution is to focus on a subset of data for divergence time estimation, which can significantly reduce the computational burdens and avoid problems with data heterogeneity that may bias results.

RESULTS

In this study, we obtained thousands of ultraconserved elements (UCEs) from 130 extant galliform taxa, including representatives of all genera, to determine the divergence times throughout galliform history. We tested the effects of different "gene shopping" schemes on divergence time estimation using a carefully, and previously validated, set of fossils. Our results found commonly used clock-like schemes may not be suitable for UCE dating (or other data types) where some loci have little information. We suggest use of partitioning (e.g., PartitionFinder) and selection of tree-like partitions may be good strategies to select a subset of data for divergence time estimation from UCEs. Our galliform time tree is largely consistent with other molecular clock studies of mitochondrial and nuclear loci. With our increased taxon sampling, a well-resolved topology, carefully vetted fossil calibrations, and suitable molecular dating methods, we obtained a high quality galliform time tree.

CONCLUSIONS

We provide a robust galliform backbone time tree that can be combined with more fossil records to further facilitate our understanding of the evolution of Galliformes and can be used as a resource for comparative and biogeographic studies in this group.

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.

Evolution and dispersal of snakes across the Cretaceous-Paleogene mass extinction

Mass extinctions have repeatedly shaped global biodiversity. The Cretaceous-Paleogene (K-Pg) mass extinction caused the demise of numerous vertebrate groups, and its aftermath saw the rapid diversification of surviving mammals, birds, frogs, and teleost fishes. However, the effects of the K-Pg extinction on the evolution of snakes-a major clade of predators comprising over 3,700 living species-remains poorly understood. Here, we combine an extensive molecular dataset with phylogenetically and stratigraphically constrained fossil calibrations to infer an evolutionary timescale for Serpentes. We reveal a potential diversification among crown snakes associated with the K-Pg mass extinction, led by the successful colonisation of Asia by the major extant clade Afrophidia. Vertebral morphometrics suggest increasing morphological specialisation among marine snakes through the Paleogene. The dispersal patterns of snakes following the K-Pg underscore the importance of this mass extinction event in shaping Earth's extant vertebrate faunas.

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

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

Comparing diversification rates in lakes, rivers, and the sea

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

Spiny and soft-rayed fin domains in acanthomorph fish are established through a BMP-gremlin-shh signaling network

With over 18,000 species, the Acanthomorpha, or spiny-rayed fishes, form the largest and arguably most diverse radiation of vertebrates. One of the key novelties that contributed to their evolutionary success are the spiny rays in their fins that serve as a defense mechanism. We investigated the patterning mechanisms underlying the differentiation of median fin Anlagen into discrete spiny and soft-rayed domains during the ontogeny of the direct-developing cichlid fish Astatotilapia burtoni Distinct transcription factor signatures characterize these two fin domains, whereby mutually exclusive expression of hoxa13a/b with alx4a/b and tbx2b marks the spine to soft-ray boundary. The soft-ray domain is established by BMP inhibition via gremlin1b, which synergizes in the posterior fin with shh secreted from a zone of polarizing activity. Modulation of BMP signaling by chemical inhibition or gremlin1b CRISPR/Cas9 knockout induces homeotic transformations of spines into soft rays and vice versa. The expression of spine and soft-ray genes in nonacanthomorph fins indicates that a combination of exaptation and posterior expansion of an ancestral developmental program for the anterior fin margin allowed the evolution of robustly individuated spiny and soft-rayed domains. We propose that a repeated exaptation of such pattern might underly the convergent evolution of anterior spiny-fin elements across fishes.

Phylogenomics illuminates the evolution of bobtail and bottletail squid (order Sepiolida)

Bobtail and bottletail squid are small cephalopods with striking anti-predatory defensive mechanisms, bioluminescence, and complex morphology; that inhabit nektobenthic and pelagic environments around the world's oceans. Yet, the evolution and diversification of these animals remain unclear. Here, we used shallow genome sequencing of thirty-two bobtail and bottletail squids to estimate their evolutionary relationships and divergence time. Our phylogenetic analyses show that each of Sepiadariidae, Sepiolidae, and the three subfamilies of the Sepiolidae are monophyletic. We found that the ancestor of the Sepiolinae very likely possessed a bilobed light organ with bacteriogenic luminescence. Sepiolinae forms a sister group to Rossinae and Heteroteuthinae, and split into Indo-Pacific and Atlantic-Mediterranean lineages. The origin of these lineages coincides with the end of the Tethys Sea and the separation of these regions during the Eocene and the beginning of the Oligocene. We demonstrated that sepiolids radiated after the Late Cretaceous and that major biogeographic events might have shaped their distribution and speciation.

Chromosome-level genome assembly of the endangered humphead wrasse Cheilinus undulatus: Insight into the expansion of opsin genes in fishes

Wrasses are dominant components of major coral reef systems. Among wrasses, Cheilinus undulatus is an endangered species with high economic and ecological value that exhibits sex reversal of females to males, while sexual selection occurs in breeding aggregations. However, the molecular-associated mechanism underlying this remains unclear. Opsin gene diversification is regarded as a potent force in sexual selection. Here we present a genome assembly of C. undulatus, using Illumina, Nanopore and Hi-C sequencing. The 1.17 Gb genome was generated from 328 contigs with an N50 length of 16.5 Mb and anchored to 24 chromosomes. In total, 22,218 genes were functionally annotated, and 96.36% of BUSCO genes were fully represented. Transcriptomic analyses showed that 96.79% of the predicted genes were expressed. Transposons were most abundant, accounting for 39.88% of the genome, with low divergence, owing to their evolution with close species ~60.53 million years ago. In total, 567/1,826 gene families were expanded and contracted in the reconstructed phylogeny, respectively. Forty-six genes were under positive selection. Comparative genomic analyses with other fish revealed expansion of opsin SWS2B, LWS1 and Rh2. The elevated duplicates of SWS2B were generated by gene conversions via transposition of transposons followed by nonallelic homologous recombination. Amino acid substitutions of opsin paralogues occurred at key tuning sites, causing a spectral shift in maximal absorbance of visual pigment to capture functional changes. Among these opsin genes, SWS2B-3 and 4 and Rh1 are expressed in the retina. The genome sequence of C. undulatus provides valuable resources for future investigation of the conservation, evolution and behaviour of fishes.