Alex Dornburg

Nine exceptional radiations plus high turnover explain species diversity in jawed vertebrates

The uneven distribution of species richness is a fundamental and unexplained pattern of vertebrate biodiversity. Although species richness in groups like mammals, birds, or teleost fishes is often attributed to accelerated cladogenesis, we lack a quantitative conceptual framework for identifying and comparing the exceptional changes of tempo in vertebrate evolutionary history. We develop MEDUSA, a stepwise approach based upon the Akaike information criterion for detecting multiple shifts in birth and death rates on an incompletely resolved phylogeny. We apply MEDUSA incompletely to a diversity tree summarizing both evolutionary relationships and species richness of 44 major clades of jawed vertebrates. We identify 9 major changes in the tempo of gnathostome diversification; the most significant of these lies at the base of a clade that includes most of the coral-reef associated fishes as well as cichlids and perches. Rate increases also underlie several well recognized tetrapod radiations, including most modern birds, lizards and snakes, ostariophysan fishes, and most eutherian mammals. In addition, we find that large sections of the vertebrate tree exhibit nearly equal rates of origination and extinction, providing some of the first evidence from molecular data for the importance of faunal turnover in shaping biodiversity. Together, these results reveal living vertebrate biodiversity to be the product of volatile turnover punctuated by 6 accelerations responsible for >85% of all species as well as 3 slowdowns that have produced “living fossils.” In addition, by revealing the timing of the exceptional pulses of vertebrate diversification as well as the clades that experience them, our diversity tree provides a framework for evaluating particular causal hypotheses of vertebrate radiations.

Resolution of ray-finned fish phylogeny and timing of diversification

Ray-finned fishes make up half of all living vertebrate species. Nearly all ray-finned fishes are teleosts, which include most commercially important fish species, several model organisms for genomics and developmental biology, and the dominant component of marine and freshwater vertebrate faunas. Despite the economic and scientific importance of ray-finned fishes, the lack of a single comprehensive phylogeny with corresponding divergence-time estimates has limited our understanding of the evolution and diversification of this radiation. Our analyses, which use multiple nuclear gene sequences in conjunction with 36 fossil age constraints, result in a well-supported phylogeny of all major ray-finned fish lineages and molecular age estimates that are generally consistent with the fossil record. This phylogeny informs three long-standing problems: specifically identifying elopomorphs (eels and tarpons) as the sister lineage of all other teleosts, providing a unique hypothesis on the radiation of early euteleosts, and offering a promising strategy for resolution of the “bush at the top of the tree” that includes percomorphs and other spiny-finned teleosts. Contrasting our divergence time estimates with studies using a single nuclear gene or whole mitochondrial genomes, we find that the former underestimates ages of the oldest ray-finned fish divergences, but the latter dramatically overestimates ages for derived teleost lineages. Our time-calibrated phylogeny reveals that much of the diversification leading to extant groups of teleosts occurred between the late Mesozoic and early Cenozoic, identifying this period as the “Second Age of Fishes.”

A comprehensive phylogeny of birds (Aves) using targeted next-generation DNA sequencing

Although reconstruction of the phylogeny of living birds has progressed tremendously in the last decade, the evolutionary history of Neoaves—a clade that encompasses nearly all living bird species—remains the greatest unresolved challenge in dinosaur systematics. Here we investigate avian phylogeny with an unprecedented scale of data: >390,000 bases of genomic sequence data from each of 198 species of living birds, representing all major avian lineages, and two crocodilian outgroups. Sequence data were collected using anchored hybrid enrichment, yielding 259 nuclear loci with an average length of 1,523 bases for a total data set of over 7.8 × 107 bases. Bayesian and maximum likelihood analyses yielded highly supported and nearly identical phylogenetic trees for all major avian lineages. Five major clades form successive sister groups to the rest of Neoaves: (1) a clade including nightjars, other caprimulgiforms, swifts, and hummingbirds; (2) a clade uniting cuckoos, bustards, and turacos with pigeons, mesites, and sandgrouse; (3) cranes and their relatives; (4) a comprehensive waterbird clade, including all diving, wading, and shorebirds; and (5) a comprehensive landbird clade with the enigmatic hoatzin (Opisthocomus hoazin) as the sister group to the rest. Neither of the two main, recently proposed Neoavian clades—Columbea and Passerea—were supported as monophyletic. The results of our divergence time analyses are congruent with the palaeontological record, supporting a major radiation of crown birds in the wake of the Cretaceous–Palaeogene (K–Pg) mass extinction.

Phylogeny and tempo of diversification in the superradiation of spiny-rayed fishes

Spiny-rayed fishes, or acanthomorphs, comprise nearly one-third of all living vertebrates. Despite their dominant role in aquatic ecosystems, the evolutionary history and tempo of acanthomorph diversification is poorly understood. We investigate the pattern of lineage diversification in acanthomorphs by using a well-resolved time-calibrated phylogeny inferred from a nuclear gene supermatrix that includes 520 acanthomorph species and 37 fossil age constraints. This phylogeny provides resolution for what has been classically referred to as the “bush at the top” of the teleost tree, and indicates acanthomorphs originated in the Early Cretaceous. Paleontological evidence suggests acanthomorphs exhibit a pulse of morphological diversification following the end Cretaceous mass extinction; however, the role of this event on the accumulation of living acanthomorph diversity remains unclear. Lineage diversification rates through time exhibit no shifts associated with the end Cretaceous mass extinction, but there is a global decrease in lineage diversification rates 50 Ma that occurs during a period when morphological disparity among fossil acanthomorphs increases sharply. Analysis of clade-specific shifts in diversification rates reveal that the hyperdiversity of living acanthomorphs is highlighted by several rapidly radiating lineages including tunas, gobies, blennies, snailfishes, and Afro-American cichlids. These lineages with high diversification rates are not associated with a single habitat type, such as coral reefs, indicating there is no single explanation for the success of acanthomorphs, as exceptional bouts of diversification have occurred across a wide array of marine and freshwater habitats.

Ancient climate change, antifreeze, and the evolutionary diversification of Antarctic fishes

The Southern Ocean around Antarctica is among the most rapidly warming regions on Earth, but has experienced episodic climate change during the past 40 million years. It remains unclear how ancient periods of climate change have shaped Antarctic biodiversity. The origin of antifreeze glycoproteins (AFGPs) in Antarctic notothenioid fishes has become a classic example of how the evolution of a key innovation in response to climate change can drive adaptive radiation. By using a time-calibrated molecular phylogeny of notothenioids and reconstructed paleoclimate, we demonstrate that the origin of AFGP occurred between 42 and 22 Ma, which includes a period of global cooling approximately 35 Ma. However, the most species-rich lineages diversified and evolved significant ecological differences at least 10 million years after the origin of AFGPs, during a second cooling event in the Late Miocene (11.6–5.3 Ma). This pattern indicates that AFGP was not the sole trigger of the notothenioid adaptive radiation. Instead, the bulk of the species richness and ecological diversity originated during the Late Miocene and into the Early Pliocene, a time coincident with the origin of polar conditions and increased ice activity in the Southern Ocean. Our results challenge the current understanding of the evolution of Antarctic notothenioids suggesting that the ecological opportunity that underlies this adaptive radiation is not linked to a single trait, but rather to a combination of freeze avoidance offered by AFGPs and subsequent exploitation of new habitats and open niches created by increased glacial and ice sheet activity.

Molecular and fossil evidence place the origin of cichlid fishes long after Gondwanan rifting

Cichlid fishes are a key model system in the study of adaptive radiation, speciation and evolutionary developmental biology. More than 1600 cichlid species inhabit freshwater and marginal marine environments across several southern landmasses. This distributional pattern, combined with parallels between cichlid phylogeny and sequences of Mesozoic continental rifting, has led to the widely accepted hypothesis that cichlids are an ancient group whose major biogeographic patterns arose from Gondwanan vicariance. Although the Early Cretaceous (ca 135 Ma) divergence of living cichlids demanded by the vicariance model now represents a key calibration for teleost molecular clocks, this putative split pre-dates the oldest cichlid fossils by nearly 90 Myr. Here, we provide independent palaeontological and relaxed-molecular-clock estimates for the time of cichlid origin that collectively reject the antiquity of the group required by the Gondwanan vicariance scenario. The distribution of cichlid fossil horizons, the age of stratigraphically consistent outgroup lineages to cichlids and relaxed-clock analysis of a DNA sequence dataset consisting of 10 nuclear genes all deliver overlapping estimates for crown cichlid origin centred on the Palaeocene (ca 65–57 Ma), substantially post-dating the tectonic fragmentation of Gondwana. Our results provide a revised macroevolutionary time scale for cichlids, imply a role for dispersal in generating the observed geographical distribution of this important model clade and add to a growing debate that questions the dominance of the vicariance paradigm of historical biogeography.

Integrating Fossil Preservation Biases in the Selection of Calibrations for Molecular Divergence Time Estimation

The selection of fossil data to use as calibration age priors in molecular divergence time estimates inherently links neontological methods with paleontological theory. However, few neontological studies have taken into account the possibility of a taphonomic bias in the fossil record when developing approaches to fossil calibration selection. The Sppil-Rongis effect may bias the first appearance of a lineage toward the recent causing most objective calibration selection approaches to erroneously exclude appropriate calibrations or to incorporate multiple calibrations that are too young to accurately represent the divergence times of target lineages. Using turtles as a case study, we develop a Bayesian extension to the fossil selection approach developed by Marshall (2008. A simple method for bracketing absolute divergence times on molecular phylogenies using multiple fossil calibrations points. Am. Nat. 171:726-742) that takes into account this taphonomic bias. Our method has the advantage of identifying calibrations that may bias age estimates to be too recent while incorporating uncertainty in phylogenetic parameter estimates such as tree topology and branch lengths. Additionally, this method is easily adapted to assess the consistency of potential calibrations to any one calibration in the candidate pool.

THE INFLUENCE OF AN INNOVATIVE LOCOMOTOR STRATEGY ON THE PHENOTYPIC DIVERSIFICATION OF TRIGGERFISH (FAMILY: BALISTIDAE)

Innovations in locomotor morphology have been invoked as important drivers of vertebrate diversification, although the influence of novel locomotion strategies on marine fish diversification remains largely unexplored. Using triggerfish as a case study, we determine whether the evolution of the distinctive synchronization of enlarged dorsal and anal fins that triggerfish use to swim may have catalyzed the ecological diversification of the group. By adopting a comparative phylogenetic approach to quantify median fin and body shape integration and to assess the tempo of functional and morphological evolution in locomotor traits, we find that: (1) functional and morphological components of the locomotive system exhibit a strong signal of correlated evolution; (2) triggerfish partitioned locomotor morphological and functional spaces early in their history; and (3) there is no strong evidence that a pulse of lineage diversification accompanied the major episode of phenotypic diversification. Together these findings suggest that the acquisition of a distinctive mode of locomotion drove an early radiation of shape and function in triggerfish, but not an early radiation of species.

Phylogenetic informativeness reconciles ray-finned fish molecular divergence times

BackgroundDiscordance among individual molecular age estimates, or between molecular age estimates and the fossil record, is observed in many clades across the Tree of Life. This discordance is attributed to a variety of variables including calibration age uncertainty, calibration placement, nucleotide substitution rate heterogeneity, or the specified molecular clock model. However, the impact of changes in phylogenetic informativeness of individual genes over time on phylogenetic inferences is rarely analyzed. Using nuclear and mitochondrial sequence data for ray-finned fishes (Actinopterygii) as an example, we extend the utility of phylogenetic informativeness profiles to predict the time intervals when nucleotide substitution saturation results in discordance among molecular ages estimated.ResultsWe demonstrate that even with identical calibration regimes and molecular clock methods, mitochondrial based molecular age estimates are systematically older than those estimated from nuclear sequences. This discordance is most severe for highly nested nodes corresponding to more recent (i.e., Jurassic-Recent) divergences. By removing data deemed saturated, we reconcile the competing age estimates and highlight that the older mtDNA based ages were driven by nucleotide saturation.ConclusionsHomoplasious site patterns in a DNA sequence alignment can systematically bias molecular divergence time estimates. Our study demonstrates that PI profiles can provide a non-arbitrary criterion for data exclusion to mitigate the influence of homoplasy on time calibrated branch length estimates. Analyses of actinopterygian molecular clocks demonstrate that scrutiny of the time scale on which sequence data is informative is a fundamental, but generally overlooked, step in molecular divergence time estimation.

The Influence of Model Averaging on Clade Posteriors: An Example Using the Triggerfishes (Family Balistidae)

Although substantial uncertainty typically surrounds the choice of the best model in most phylogenetic analyses, little is known about how accommodating this uncertainty affects phylogenetic inference. Here we explore the influence of Bayesian model averaging on the phylogenetic inference of the triggerfishes (Family: Balistidae), a charismatic group of reef fishes. We focus on clade support as this area has received little attention and is typically one of the most important outcomes of phylogenetic studies. We present a novel phylogenetic hypothesis for the family Balistidae based on an analysis of two mitochondrial (12S, 16S) and three nuclear genes (TMO-4C4, Rhodopsin, RAG1) sampled from 26 ingroup species. Despite the presence of substantial model uncertainty in almost all partitions of our data, we found model-averaged topologies and clade posteriors to be nearly identical to those conditioned on a single model. Furthermore, statistical comparison of clade posteriors using the Wilcoxon signed-rank test revealed no significant differences. Our results suggest that although current model-selection approaches are likely to lead to overparameterization of the substitution model, the consequences of conditioning on this overparameterized model are likely to be mild. Our phylogenetic results strongly support the monophyly of the triggerfishes but suggest that the genera Balistoides and Pseudobalistes are polyphyletic. Divergence time estimation supports a Miocene origin of the crown group. Despite the presence of several young species-rich subclades, statistical analysis of temporal diversification patterns reveals no significant increase in the rate of cladogenesis across geologic time intervals.