Mass extinction in tetraodontiform fishes linked to the Palaeocene-Eocene thermal maximum

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.

No phylogenetic evidence for angiosperm mass extinction at the Cretaceous-Palaeogene (K-Pg) boundary

The Cretaceous-Palaeogene mass extinction event (K-Pg) witnessed upwards of 75% of animal species going extinct, most notably among these are the non-avian dinosaurs. A major question in macroevolution is whether this extinction event influenced the rise of flowering plants (angiosperms). The fossil record suggests that the K-Pg event had a strong regional impact on angiosperms with up to 75% species extinctions, but only had a minor impact on the extinction rates of major lineages (families and orders). Phylogenetic evidence for angiosperm extinction dynamics through time remains unexplored. By analysing two angiosperm mega-phylogenies containing approximately 32 000-73 000 extant species, here we show relatively constant extinction rates throughout geological time and no evidence for a mass extinction at the K-Pg boundary. Despite high species-level extinction observed in the fossil record, our results support the macroevolutionary resilience of angiosperms to the K-Pg mass extinction event via survival of higher lineages.

Impacts of Taxon-Sampling Schemes on Bayesian Tip Dating Under the Fossilized Birth-Death Process

Evolutionary timescales can be inferred by molecular-clock analyses of genetic data and fossil evidence. Bayesian phylogenetic methods such as tip dating provide a powerful framework for inferring evolutionary timescales, but the most widely used priors for tree topologies and node times often assume that present-day taxa have been sampled randomly or exhaustively. In practice, taxon sampling is often carried out so as to include representatives of major lineages, such as orders or families. We examined the impacts of different densities of diversified sampling on Bayesian tip dating on unresolved fossilized birth-death (FBD) trees, in which fossil taxa are topologically constrained but their exact placements are averaged out. We used synthetic data generated by simulations of nucleotide sequence evolution, fossil occurrences, and diversified taxon sampling. Our analyses under the diversified-sampling FBD process show that increasing taxon-sampling density does not necessarily improve divergence-time estimates. However, when informative priors were specified for the root age or when tree topologies were fixed to those used for simulation, the performance of tip dating on unresolved FBD trees maintains its accuracy and precision or improves with taxon-sampling density. By exploring three situations in which models are mismatched, we find that including all relevant fossils, without pruning off those that are incompatible with the diversified-sampling FBD process, can lead to underestimation of divergence times. Our reanalysis of a eutherian mammal data set confirms some of the findings from our simulation study, and reveals the complexity of diversified taxon sampling in phylogenomic data sets. In highlighting the interplay of taxon-sampling density and other factors, the results of our study have practical implications for using Bayesian tip dating to infer evolutionary timescales across the Tree of Life. [Bayesian tip dating; eutherian mammals; fossilized birth-death process; phylogenomics; taxon sampling.].

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.

Convergence is Only Skin Deep: Craniofacial Evolution in Electric Fishes from South America and Africa (Apteronotidae and Mormyridae)

Apteronotidae and Mormyridae are species-rich clades of weakly electric fishes from Neotropical and Afrotropical freshwaters, respectively, known for their high morphological disparity and often regarded as a classic example of convergent evolution. Here, we use CT-imaging and 3D geometric morphometrics to quantify disparity in craniofacial morphologies, and to test the hypothesis of convergent skull-shape evolution in a phylogenetic context. For this study, we examined 391 specimens representing 78 species of Apteronotidae and Mormyridae including 30 of 37 (81%) of all valid genera with the goal to sample most of the craniofacial disparity known in these clades. We found no overlap between Apteronotidae and Mormyridae in the skull-shape morphospace using PCA and a common landmark scheme, and therefore no instances of complete phenotypic convergence. Instead, we found multiple potential instances of incomplete convergence, and at least one parallel shift among electric fish clades. The greatest components of shape variance in both families are the same as observed for most vertebrate clades: heterocephaly (i.e., opposite changes in relative sizes of the snout and braincase regions of the skull), and heterorhynchy (i.e., dorsoventral changes in relative snout flexion and mouth position). Mormyrid species examined here exhibit less craniofacial disparity than do apteronotids, potentially due to constraints associated with a larger brain size, ecological constraints related to food-type availability. Patterns of craniofacial evolution in these two clades depict a complex story of phenotypic divergence and convergence in which certain superficial similarities of external morphology obscure deeper osteological and presumably developmental differences of skull form and function. Among apteronotid and mormyrid electric fishes, craniofacial convergence is only skin deep.

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.

Evolutionary dynamics of the OR gene repertoire in teleost fishes: evidence of an association with changes in olfactory epithelium shape

Teleost fishes perceive their environment through a range of sensory modalities, among which olfaction often plays an important role. Richness of the olfactory repertoire depends on the diversity of receptors coded by homologous genes classified into four families: OR, TAAR, VR1, and VR2. Herein, we focus on the OR gene repertoire. While independent large contractions of the OR gene repertoire associated with ecological transitions have been found in mammals, little is known about the diversity of the OR gene repertoire and its evolution in teleost fishes, a group that includes more than 34,000 living species. We analyzed genomes of 163 species representing diversity in this large group. We found a large range of variation in the number of functional OR genes, from 15 in the Broad-nose Pipefish Syngnathus typhle and the Ocean Sunfish Mola mola, to 429 in the Zig-zag Eel Mastacembelus armatus. The number of OR genes was higher in species when a multilamellar olfactory rosette was present. Moreover, the number of lamellae was correlated with the richness of the OR gene repertoire. While a slow and balanced birth-and-death process generally drives the evolution of the OR gene repertoire, we inferred several episodes of high rates of gene loss, sometimes followed by large gains in the number of OR genes. These gains coincide with morphological changes of the olfactory organ and suggest a strong functional association between changes in the morphology and the evolution of the OR gene repertoire.

A Simulation-Based Evaluation of Tip-Dating Under the Fossilized Birth-Death Process

Bayesian molecular dating is widely used to study evolutionary timescales. This procedure usually involves phylogenetic analysis of nucleotide sequence data, with fossil-based calibrations applied as age constraints on internal nodes of the tree. An alternative approach is tip-dating, which explicitly includes fossil data in the analysis. This can be done, for example, through the joint analysis of molecular data from present-day taxa and morphological data from both extant and fossil taxa. In the context of tip-dating, an important development has been the fossilized birth-death process, which allows non-contemporaneous tips and sampled ancestors while providing a model of lineage diversification for the prior on the tree topology and internal node times. However, tip-dating with fossils faces a number of considerable challenges, especially, those associated with fossil sampling and evolutionary models for morphological characters. We conducted a simulation study to evaluate the performance of tip-dating using the fossilized birth-death model. We simulated fossil occurrences and the evolution of nucleotide sequences and morphological characters under a wide range of conditions. Our analyses of these data show that the number and the maximum age of fossil occurrences have a greater influence than the degree of among-lineage rate variation or the number of morphological characters on estimates of node times and the tree topology. Tip-dating with the fossilized birth-death model generally performs well in recovering the relationships among extant taxa but has difficulties in correctly placing fossil taxa in the tree and identifying the number of sampled ancestors. The method yields accurate estimates of the ages of the root and crown group, although the precision of these estimates varies with the probability of fossil occurrence. The exclusion of morphological characters results in a slight overestimation of node times, whereas the exclusion of nucleotide sequences has a negative impact on inference of the tree topology. Our results provide an overview of the performance of tip-dating using the fossilized birth-death model, which will inform further development of the method and its application to key questions in evolutionary biology.

Post-Cretaceous bursts of evolution along the benthic-pelagic axis in marine fishes

Ecological opportunity arising in the aftermath of mass extinction events is thought to be a powerful driver of evolutionary radiations. Here, we assessed how the wake of the Cretaceous-Palaeogene (K-Pg) mass extinction shaped diversification dynamics in a clade of mostly marine fishes (Carangaria), which comprises a disparate array of benthic and pelagic dwellers including some of the most astonishing fish forms (e.g. flatfishes, billfishes, remoras, archerfishes). Analyses of lineage diversification show time-heterogeneous rates of lineage diversification in carangarians, with highest rates reached during the Palaeocene. Likewise, a remarkable proportion of Carangaria's morphological variation originated early in the history of the group and in tandem with a marked incidence of habitat shifts. Taken together, these results suggest that all major lineages and body plans in Carangaria originated in an early burst shortly after the K-Pg mass extinction, which ultimately allowed the occupation of newly released niches along the benthic-pelagic habitat axis.

Two pulses of morphological diversification in Pacific pelagic fishes following the Cretaceous-Palaeogene mass extinction

Molecular phylogenies suggest some major radiations of open-ocean fish clades occurred roughly coincident with the Cretaceous-Palaeogene (K/Pg) boundary, however the timing and nature of this diversification is poorly constrained. Here, we investigate evolutionary patterns in ray-finned fishes across the K/Pg mass extinction 66 million years ago (Ma), using microfossils (isolated teeth) preserved in a South Pacific sediment core spanning 72-43 Ma. Our record does not show significant turnover of fish tooth morphotypes at the K/Pg boundary: only two of 48 Cretaceous tooth morphotypes disappear at the event in the South Pacific, a rate no different from background extinction. Capture-mark-recapture analysis finds two pulses of origination in fish tooth morphotypes following the mass extinction. The first pulse, at approximately 64 Ma, included short-lived teeth, as well as forms that contribute to an expansion into novel morphospace. A second pulse, centred at approximately 58 Ma, produced morphotype novelty in a different region of morphospace from the first pulse, and contributed significantly to Eocene tooth morphospace occupation. There was no significant increase in origination rates or expansion into novel morphospace during the early or middle Eocene, despite a near 10-fold increase in tooth abundance during that interval. Our results suggest that while the K/Pg event had a minor impact on fish diversity in terms of extinction, the removal of the few dominant Cretaceous morphotypes triggered a sequence of origination events allowing fishes to rapidly diversify morphologically, setting the stage for exceptional levels of ray-finned fish diversity in the Cenozoic.

Little lasting impact of the Paleocene-Eocene Thermal Maximum on shallow marine molluscan faunas

Global warming, acidification, and oxygen stress at the Paleocene-Eocene Thermal Maximum (PETM) are associated with severe extinction in the deep sea and major biogeographic and ecologic changes in planktonic and terrestrial ecosystems, yet impacts on shallow marine macrofaunas are obscured by the incompleteness of shelf sections. We analyze mollusk assemblages bracketing (but not including) the PETM and find few notable lasting impacts on diversity, turnover, functional ecology, body size, or life history of important clades. Infaunal and chemosymbiotic taxa become more common, and body size and abundance drop in one clade, consistent with hypoxia-driven selection, but within-clade changes are not generalizable across taxa. While an unrecorded transient response is still possible, the long-term evolutionary impact is minimal. Adaptation to already-warm conditions and slow release of CO₂ relative to the time scale of ocean mixing likely buffered the impact of PETM climate change on shelf faunas.