Charting the Course of Pinniped Evolution: insights from molecular phylogeny and fossil record integration

Pinnipeds (seals, sea lions, walruses, and their fossil relatives) are one of the most successful mammalian clades to live in the oceans. Despite a well-resolved molecular phylogeny and a global fossil record, a complete understanding of their macroevolutionary dynamics remains hampered by a lack of formal analyses that combine these two rich sources of information. We used a meta-analytic approach to infer the most densely sampled pinniped phylogeny to-date (36 recent and 93 fossil taxa) and used phylogenetic paleobiological methods to study their diversification dynamics and biogeographic history. Pinnipeds mostly diversified at constant rates. Walruses however experienced rapid turnover in which extinction rates ultimately exceeded speciation rates from 12-6 Ma, possibly due to changing sea-levels and/or competition with otariids (eared seals). Historical biogeographic analyses including fossil data allowed us to confidently identify the North Pacific and the North Atlantic (plus or minus Paratethys) as the ancestral ranges of Otarioidea (eared seals + walrus) and crown phocids (earless seals), respectively. Yet, despite the novel addition of stem pan-pinniped taxa, the region of origin for Pan-Pinnipedia remained ambiguous. These results suggest further avenues of study in pinnipeds and provide a framework for investigating other groups with substantial extinct and extant diversity.

Bayesian tip-dated phylogeny and biogeography of Cissampelideae (Menispermaceae): Mitigating the effects of homoplastic morphological characters

The integration of morphological and molecular data is essential to understand the affinities of fossil taxa and spatio-temporal evolutionary processes of organisms. However, homoplastic morphological characters can mislead the placement of fossil taxa and impact downstream analyses. Here, we provide an example of how to mitigate effectively the effect of morphological homoplasy on the placement of fossil taxa and biogeographic inferences of Cissampelideae. We assembled three data types, morphological data only, morphological data with a molecular scaffold and combined morphological and molecular data. By removing high-level homoplastic morphological data or reweighting the morphological characters, we conducted 15 parsimony, 12 undated Bayesian and four dated Bayesian analyses. Our results show that the 14 selected Cissampelideae fossil taxa are placed poorly when based only on morphological data, but the addition of molecular scaffold and combination of morphological and molecular data greatly improve the resolution of fossil nodes. We raise the monotypic Stephania subg. Botryodiscia to generic status and discover that three fossils previously assigned to Stephania should be members of Diploclisia. The Bayesian tip-dated tree recovered by removing homoplastic morphological characters with a Rescaled Consistency Index <0.25 has the highest stratigraphic fit and consequently generates more reasonable biogeographic reconstruction for Cissampelideae. Cissampelideae began to diversify in Asia in the latest Cretaceous and subsequently dispersed to South America around the Cretaceous-Palaeogene boundary. Two dispersal events from Asia to Africa occurred in the Early Eocene and the Late Eocene-Late Oligocene, respectively. These findings provide guidelines and practical methods for mitigating the effects of homoplastic morphological characters on fossil placements and Bayesian tip-dating, as well as insights into the past tropical floristic exchanges among different continents.

Reconciling fossils with phylogenies reveals the origin and macroevolutionary processes explaining the global cycad biodiversity

The determinants of biodiversity patterns can be understood using macroevolutionary analyses. The integration of fossils into phylogenies offers a deeper understanding of processes underlying biodiversity patterns in deep time. Cycadales are considered a relict of a once more diverse and globally distributed group but are restricted to low latitudes today. We still know little about their origin and geographic range evolution. Combining molecular data for extant species and leaf morphological data for extant and fossil species, we study the origin of cycad global biodiversity patterns through Bayesian total-evidence dating analyses. We assess the ancestral geographic origin and trace the historical biogeography of cycads with a time-stratified process-based model. Cycads originated in the Carboniferous on the Laurasian landmass and expanded in Gondwana in the Jurassic. Through now-vanished continental connections, Antarctica and Greenland were crucial biogeographic crossroads for cycad biogeography. Vicariance is an essential speciation mode in the deep and recent past. Their latitudinal span increased in the Jurassic and restrained toward subtropical latitudes in the Neogene in line with biogeographic inferences of high-latitude extirpations. We show the benefits of integrating fossils into phylogenies to estimate ancestral areas of origin and to study evolutionary processes explaining the global distribution of present-day relict groups.

The dynamic adaptive landscape of cetacean body size

Adaptive landscapes are central to evolutionary theory, forming a conceptual bridge between micro- and macroevolution.^1,2,3,4 Evolution by natural selection across an adaptive landscape should drive lineages toward fitness peaks, shaping the distribution of phenotypic variation within and among clades over evolutionary timescales.⁵ The location and breadth of these peaks in phenotypic space can also evolve,⁴ but whether phylogenetic comparative methods can detect such patterns has largely remained unexplored.⁶ Here, we characterize the global and local adaptive landscape for total body length in cetaceans (whales, dolphins, and relatives), a trait that spans 5 orders of magnitude, across their ∼53 million year evolutionary history. Using phylogenetic comparative methods, we analyze shifts in long-term mean body length⁷ and directional changes in average trait values⁸ for 345 living and fossil cetacean taxa. Remarkably, we find that the global macroevolutionary adaptive landscape of cetacean body length is relatively flat, with very few peak shifts occurring after cetaceans entered the oceans. Local peaks are more numerous and manifest as trends along branches linked to specific adaptations. These results contrast with previous studies using only extant taxa,⁹ highlighting the vital role of fossil data for understanding macroevolution.^10,11,12 Our results indicate that adaptive peaks are dynamic and are associated with subzones of local adaptations, creating moving targets for species adaptation. In addition, we identify limits in our ability to detect some evolutionary patterns and processes and suggest that multiple approaches are required to characterize complex hierarchical patterns of adaptation in deep time.

Cross-disciplinary information for understanding macroevolution

Many different macroevolutionary models can produce the same observations. Despite efforts in building more complex and realistic models, it may still be difficult to distinguish the processes that have generated the biodiversity we observe. In this opinion we argue that we can make new progress by reaching out across disciplines, relying on independent data and theory to constrain macroevolutionary inference. Using mainly paleontological insights and data, we illustrate how we can eliminate less plausible or implausible models, and/or parts of parameter space, while applying comparative phylogenetic approaches. We emphasize that such cross-disciplinary insights and data can be drawn between many other disciplines relevant to macroevolution. We urge cross-disciplinary training, and collaboration using common-use databases as a platform for increasing our understanding.

Deep time extinction of largest insular ant predators and the first fossil Neoponera (Formicidae: Ponerinae) from Miocene age Dominican amber

BACKGROUND

Ponerine ants are almost exclusively predatory and comprise many of the largest known ant species. Within this clade, the genus Neoponera is among the most conspicuous Neotropical predators. We describe the first fossil member of this lineage: a worker preserved in Miocene-age Dominican amber from Hispaniola.

RESULTS

Neoponera vejestoria sp. nov. demonstrates a clear case of local extinction-there are no known extant Neoponera species in the Greater Antilles. The species is attributable to an extant and well-defined species group in the genus, which suggests the group is older than previously estimated. Through CT scan reconstruction and linear morphometrics, we reconstruct the morphospace of extant and fossil ants to evaluate the history and evolution of predatory taxa in this island system.

CONCLUSIONS

The fossil attests to a shift in insular ecological community structure since the Miocene. The largest predatory taxa have undergone extinction on the island, but their extant relatives persist throughout the Neotropics. Neoponera vejestoria sp. nov. is larger than all other predatory ant workers known from Hispaniola, extant or extinct. Our results empirically demonstrate the loss of a functional niche associated with body size, which is a trait long hypothesized to be related to extinction risk.

Systematics of Miocene apes: State of the art of a neverending controversy

Hominoids diverged from cercopithecoids during the Oligocene in Afro-Arabia, initially radiating in that continent and subsequently dispersing into Eurasia. From the Late Miocene onward, the geographic range of hominoids progressively shrank, except for hominins, which dispersed out of Africa during the Pleistocene. Although the overall picture of hominoid evolution is clear based on available fossil evidence, many uncertainties persist regarding the phylogeny and paleobiogeography of Miocene apes (nonhominin hominoids), owing to their sparse record, pervasive homoplasy, and the decimated current diversity of this group. We review Miocene ape systematics and evolution by focusing on the most parsimonious cladograms published during the last decade. First, we provide a historical account of the progress made in Miocene ape phylogeny and paleobiogeography, report an updated classification of Miocene apes, and provide a list of Miocene ape species-locality occurrences together with an analysis of their paleobiodiversity dynamics. Second, we discuss various critical issues of Miocene ape phylogeny and paleobiogeography (hylobatid and crown hominid origins, plus the relationships of Oreopithecus) in the light of the highly divergent results obtained from cladistic analyses of craniodental and postcranial characters separately. We conclude that cladistic efforts to disentangle Miocene ape phylogeny are potentially biased by a long-branch attraction problem caused by the numerous postcranial similarities shared between hylobatids and hominids-despite the increasingly held view that they are likely homoplastic to a large extent, as illustrated by Sivapithecus and Pierolapithecus-and further aggravated by abundant missing data owing to incomplete preservation. Finally, we argue that-besides the recovery of additional fossils, the retrieval of paleoproteomic data, and a better integration between cladistics and geometric morphometrics-Miocene ape phylogenetics should take advantage of total-evidence (tip-dating) Bayesian methods of phylogenetic inference combining morphologic, molecular, and chronostratigraphic data. This would hopefully help ascertain whether hylobatid divergence was more basal than currently supported.

Teleonomy: Revisiting a Proposed Conceptual Replacement for Teleology

The concept of teleonomy has been attracting renewed attention recently. This is based on the idea that teleonomy provides a useful conceptual replacement for teleology, and even that it constitutes an indispensable resource for thinking biologically about purposes. However, both these claims are open to question. We review the history of teleological thinking from Greek antiquity to the modern period to illuminate the tensions and ambiguities that emerged when forms of teleological reasoning interacted with major developments in biological thought. This sets the stage for an examination of Pittendrigh's (Adaptation, natural selection, and behavior. In: Roe A, Simpson GG (eds) Behavior and evolution. Yale University Press, New Haven, pp 390-416, 1958) introduction of "teleonomy" and its early uptake in the work of prominent biologists. We then explore why teleonomy subsequently foundered and consider whether the term may yet have significance for discussions of goal-directedness in evolutionary biology and philosophy of science. This involves clarifying the relationship between teleonomy and teleological explanation, as well as asking how the concept of teleonomy impinges on research at the frontiers of evolutionary theory.

Combining palaeontological and neontological data shows a delayed diversification burst of carcharhiniform sharks likely mediated by environmental change

Estimating deep-time species-level diversification processes remains challenging. Both the fossil record and molecular phylogenies allow the estimation of speciation and extinction rates, but each type of data may still provide an incomplete picture of diversification dynamics. Here, we combine species-level palaeontological (fossil occurrences) and neontological (molecular phylogenies) data to estimate deep-time diversity dynamics through process-based birth-death models for Carcharhiniformes, the most speciose shark order today. Despite their abundant fossil record dating back to the Middle Jurassic, only a small fraction of extant carcharhiniform species is recorded as fossils, which impedes relying only on the fossil record to study their recent diversification. Combining fossil and phylogenetic data, we recover a complex evolutionary history for carcharhiniforms, exemplified by several variations in diversification rates with an early low diversity period followed by a Cenozoic radiation. We further reveal a burst of diversification in the last 30 million years, which is partially recorded with fossil data only. We also find that reef expansion and temperature change can explain variations in speciation and extinction through time. These results pinpoint the primordial importance of these environmental variables in the evolution of marine clades. Our study also highlights the benefit of combining the fossil record with phylogenetic data to address macroevolutionary questions.

Relative performance of Bayesian morphological clock and parsimony methods for phylogenetic reconstructions: Insights from the case of Myomiminae and Dryomyinae glirid rodents

Extinct organisms provide crucial information about the origin and time of origination of extant groups. The importance of morphological phylogenetics for rigorously dating the tree of life is now widely recognized and has been revitalized by methodological developments such as the application of tip-dating Bayesian approaches. Traditionally, molecular clocks have been node calibrated. However, node calibrations are often unsatisfactory because they do not allow the fossil age to inform about phylogenetic hypothesis. The introduction of tip calibrations allow fossil species to be included alongside their living relatives, and the absence of molecular sequence data for these taxa to be remedied by supplementing the sequence alignments for living taxa with phenotype character matrices for both living and fossil taxa. Therefore, only phylogenetic analyses that take into account morphological characters can incorporate both fossil and extant species. Herein we present an unprecedented morphological dataset for a vast group of glirid rodents, to which different phylogenetic methodologies have been applied. We have compared the tree topologies resulting from traditional parsimony and Bayesian phylogenetic approaches and calculate stratigraphic congruence indices for each. Bayesian tip-dated clock methods seem to outperform parsimony with our dataset. The strict consensus tree recovered by tip dating invalidates the classic classification and allows dates to be proposed for the divergence and origin of the different clades.

Integrative Phylogenetics: Tools for Palaeontologists to Explore the Tree of Life

The modern era of analytical and quantitative palaeobiology has only just begun, integrating methods such as morphological and molecular phylogenetics and divergence time estimation, as well as phenotypic and molecular rates of evolution. Calibrating the tree of life to geological time is at the nexus of many disparate disciplines, from palaeontology to molecular systematics and from geochronology to comparative genomics. Creating an evolutionary time scale of the major events that shaped biodiversity is key to all of these fields and draws from each of them. Different methodological approaches and data employed in various disciplines have traditionally made collaborative research efforts difficult among these disciplines. However, the development of new methods is bridging the historical gap between fields, providing a holistic perspective on organismal evolutionary history, integrating all of the available evidence from living and fossil species. Because phylogenies with only extant taxa do not contain enough information to either calibrate the tree of life or fully infer macroevolutionary dynamics, phylogenies should preferably include both extant and extinct taxa, which can only be achieved through the inclusion of phenotypic data. This integrative phylogenetic approach provides ample and novel opportunities for evolutionary biologists to benefit from palaeontological data to help establish an evolutionary time scale and to test core macroevolutionary hypotheses about the drivers of biological diversification across various dimensions of organisms.

The Strait of Gibraltar is an ineffective palaeogeographic barrier for some flightless darkling beetles (Coleoptera: Tenebrionidae: Pimelia)

The geographic distribution of a species is shaped by its biology and by environmental and palaeogeographic factors that interact at different spatial-temporal scales, which leads to distributions and diversification patterns observed between and within lineages. The darkling beetle genus Pimelia has been diversifying for more than 31.2 Mya showing different colonization patterns after the opening of the Gibraltar Strait 5 Mya. Three of the 14 subgenera of Pimelia have populations on both sides of the Strait. Through extensive sampling and the analysis of three molecular markers, we determine levels of intra- and interspecific genetic variation, identify evolutionary lineages in subgenera, estimate their temporal origin and distribution ranges and discuss the historical basis for the geographic and diversification patterns of Pimelia around the Strait. This single geographical feature acted both as a barrier and as a dispersal route for different Pimelia species. The Strait has represented a strong barrier for the subgenus Magrebmelia since the Middle Miocene. However, the subgenera Amblyptera and Amblypteraca share repetitive signatures of post-Messinian colonization across the Strait, possibly driven by stochastic or ‘catastrophic’ events such as tsunamis. Our demographic analyses support Wallace’s hypothesis on insect dispersal stochasticity. Some taxonomic changes, including the designation of a lectotype for Pimelia maura, are also proposed.

Molecules and fossils tell distinct yet complementary stories of mammal diversification

Reconstructing the tempo at which biodiversity arose is a fundamental goal of evolutionary biologists, yet the relative merits of evolutionary-rate estimates are debated based on whether they are derived from the fossil record or time-calibrated phylogenies (timetrees) of living species. Extinct lineages unsampled in timetrees are known to "pull" speciation rates downward, but the temporal scale at which this bias matters is unclear. To investigate this problem, we compare mammalian diversification-rate signatures in a credible set of molecular timetrees (n = 5,911 species, ∼70% from DNA) to those in fossil genus durations (n = 5,320). We use fossil extinction rates to correct or "push" the timetree-based (pulled) speciation-rate estimates, finding a surge of speciation during the Paleocene (∼66-56 million years ago, Ma) between the Cretaceous-Paleogene (K-Pg) boundary and the Paleocene-Eocene Thermal Maximum (PETM). However, about two-thirds of the K-Pg-to-PETM originating taxa did not leave modern descendants, indicating that this rate signature is likely undetectable from extant lineages alone. For groups without substantial fossil records, thankfully all is not lost. Pushed and pulled speciation rates converge starting ∼10 Ma and are equal at the present day when recent evolutionary processes can be estimated without bias using species-specific "tip" rates of speciation. Clade-wide moments of tip rates also enable enriched inference, as the skewness of tip rates is shown to approximate a clade's extent of past diversification-rate shifts. Molecular timetrees need fossil-correction to address deep-time questions, but they are sufficient for shallower time questions where extinctions are fewer.

Ecological Filtering and Exaptation in the Evolution of Marine Snakes

AbstractConvergent evolution is often attributed to adaptation of form to function, but it can also result from ecological filtering, exaptation, or nonaptation. Testing among these possibilities is critical to understanding how and why morphological similarities emerge independently in multiple lineages. To address this challenge, we combined multiple preexisting phylogenetic methods to jointly estimate the habitats and morphologies of lineages within a phylogeny. We applied this approach to the invasions of snakes into the marine realm. We utilized a data set for 1,243 extant snake species consisting of newly compiled biome occupancy information and preexisting data on reproductive strategy, body mass, and environmental temperature and elevation. We find evidence for marine clades arising from a variety of aquatic and terrestrial habitats. Furthermore, there is strong evidence of ecological filtering for nonmarine ancestors that were already viviparous, had slightly larger-than-average body sizes, and lived in environments with higher-than-average temperatures and lower-than-average elevations. In aggregate, similarities among independent lineages of marine snakes result from a combination of exaptation and strong ecological filtering. Strong barriers to entry of new habitats appear to be more important than common adaptations following invasions for producing similarities among independent lineages invading a shared, novel habitat.

Bayesian Morphological Clock versus Parsimony: An insight into the relationships and dispersal events of postvacuum Cricetidae (Rodentia, Mammalia)

Establishing an evolutionary timeline is fundamental for tackling a great variety of topics in evolutionary biology, including the reconstruction of patterns of historical biogeography, coevolution and diversification. However, the tree of life is pruned by extinction and molecular data cannot be gathered for extinct lineages. Until recently methodological challenges have prevented the application of tip-dating Bayesian approaches in morphology-based fossil-only datasets. Herein, we present a morphological dataset for a group of cricetid rodents to which we apply an array of methods fairly new in palaeontology that can be used by palaeontologists for the analysis of entirely extinct clades. We compare the tree topologies obtained by traditional parsimony, time-calibrated and non-calibrated Bayesian inference phylogenetic approaches and calculate stratigraphic congruence indices for each. Bayesian tip-dated clock methods outperform parsimony in the case of our dataset, which includes highly homoplastic morphological characters. Regardless, all three topologies support the monophyly of Megacricetodontinae, Democricetodontinae and Cricetodontinae. Dispersal and speciation events inferred through Bayesian Binary Markov chain Monte Carlo and biodiversity analyses provide evidence for a correlation between biogeographic events, climatic changes and diversification in cricetids.

Empirical and Methodological Challenges to the Model-Based Inference of Diversification Rates in Extinct Clades

Changes in speciation and extinction rates are key to the dynamics of clade diversification, but attempts to infer them from phylogenies of extant species face challenges. Methods capable of synthesizing information from extant and fossil species have yielded novel insights into diversification rate variation through time, but little is known about their behavior when analyzing entirely extinct clades. Here, we use empirical and simulated data to assess how two popular methods, PyRate and Fossil BAMM, perform in this setting. We inferred the first tip-dated trees for ornithischian dinosaurs, and combined them with fossil occurrence data to test whether the clade underwent an end-Cretaceous decline. We then simulated phylogenies and fossil records under empirical constraints to determine whether macroevolutionary and preservation rates can be teased apart under paleobiologically realistic conditions. We obtained discordant inferences about ornithischian macroevolution including a long-term speciation rate decline (BAMM), mostly flat rates with a steep diversification drop (PyRate) or without one (BAMM), and episodes of implausibly accelerated speciation and extinction (PyRate). Simulations revealed little to no conflation between speciation and preservation, but yielded spuriously correlated speciation and extinction estimates while time-smearing tree-wide shifts (BAMM) or overestimating their number (PyRate). Our results indicate that the small phylogenetic datasets available to vertebrate paleontologists and the assumptions made by current model-based methods combine to yield potentially unreliable inferences about the diversification of extinct clades. We provide guidelines for interpreting the results of the existing approaches in light of their limitations, and suggest how the latter may be mitigated.

Species diversification and ecomorphological evolution in the radiation of tanagers (Passeriformes: Thraupidae)

Ecological opportunity is hypothesized to cause an early burst of species diversification and trait evolution followed by a slowdown in diversification rates as niches are filled. Nonetheless, few studies have tested these predictions empirically with ecomorphological data at the large spatial scales relevant to most of biodiversity. Tanagers (Passeriformes: Thraupidae), the largest family of songbirds, show an early burst of species diversification and provide an excellent opportunity to test one of the hallmarks of adaptive radiation: rapid ecomorphological evolution. Here, we test for an early-burst pattern of a resource-exploiting trait (bill morphology) across the radiation of tanagers using a time-calibrated molecular phylogeny and high-resolution three-dimensional surface scans of bill structure from museum study skins. Using recently developed methods of multivariate trait evolution, we find evidence for a rapid burst of bill shape evolution early in the radiation of tanagers, followed by a subsequent decrease in rates toward the present. Likewise, we show that morphological disparity is distributed among (rather than within) subclades, indicating that most of the observed bill shape disparity evolved early in the radiation of tanagers and has slowed through time. The diversification dynamics of tanagers match patterns expected from adaptive radiation and the filling of ecomorphospace.

The comparative method is not macroevolution: across-species evidence for within-species process

It is common for studies that employ the comparative method for the study of adaptation, i.e. documentation of potentially adaptive across-species patterns of trait-environment or trait-trait correlation, to be designated as "macroevolutionary." Authors are justified in using "macroevolution" in this way by appeal to definitions such as "evolution above the species level." I argue that regarding the comparative method as "macroevolutionary" is harmful because it hides in serious ways the true causal content of hypotheses tested with the comparative method. The comparative method is a means of testing hypotheses of adaptation and their alternatives. Adaptation is a population level phenomenon, involving heritable interindividual variation that is associated with fitness differences. For example, given heritable intrapopulational variation, more streamlined individuals in populations of fast-moving aquatic animals have higher locomotory efficiency and thus better survivorship and more resources directed to reproduction than less streamlined ones. Direct evidence consistent with this population-level scenario includes the observation that many unrelated species of fast-moving aquatic animals have similar streamlined shapes, an example of the comparative method. Crucial to note in this example is that although the data are observed across species, the comparative method for studying adaptation tests hypotheses regarding standard population-level natural selection with no content that can be construed as "macro." Even less "macro," individual-level developmental dynamics can limit or bias the range of variants available for selection. Calling any of these studies "macroevolutionary" implies that some additional process is at work, shrouding the need to test adaptation hypotheses and study the range of variants that can be produced in development.

Ecological Interactions and Macroevolution: A New Field with Old Roots

Linking interspecific interactions (e.g., mutualism, competition, predation, parasitism) to macroevolution (evolutionary change on deep timescales) is a key goal in biology. The role of species interactions in shaping macroevolutionary trajectories has been studied for centuries and remains a cutting-edge topic of current research. However, despite its deep historical roots, classic and current approaches to this topic are highly diverse. Here, we combine historical and contemporary perspectives on the study of ecological interactions in macroevolution, synthesizing ideas across eras to build a zoomed-out picture of the big questions at the nexus of ecology and macroevolution. We discuss the trajectory of this important and challenging field, dividing research into work done before the 1970s, research between 1970 and 2005, and work done since 2005. We argue that in response to long-standing questions in paleobiology, evidence accumulated to date has demonstrated that biotic interactions (including mutualism) can influence lineage diversification and trait evolution over macroevolutionary timescales, and we outline major open questions for future research in the field.

Lineage Identification Affects Estimates of Evolutionary Mode in Marine Snails

In order to study evolutionary pattern and process, we need to be able to accurately identify species and the evolutionary lineages from which they are derived. Determining the concordance between genetic and morphological variation of living populations, and then directly comparing extant and fossil morphological data, provides a robust approach for improving our identification of lineages through time. We investigate genetic and shell morphological variation in extant species of Penion marine snails from New Zealand, and extend this analysis into deep time using fossils. We find that genetic and morphological variation identify similar patterns and support most currently recognized extant species. However, some taxonomic over-splitting is detected due to shell size being a poor trait for species delimitation, and we identify incorrect assignment of some fossil specimens. We infer that a single evolutionary lineage (Penion sulcatus) has existed for 22 myr, with most aspects of shell shape and shell size evolving under a random walk. However, by removing samples previously classified as the extinct species P. marwicki, we instead detect morphological stasis for one axis of shell shape variation. This result demonstrates how lineage identification can change our perception of evolutionary pattern and process. [Genotyping by sequencing; geometric morphometrics; morphological evolution; Neogastropoda; phenotype; speciation; stasis.].

Evolutionary history of zoogeographical regions surrounding the Tibetan Plateau

The Tibetan Plateau (TP) and surrounding regions have one of the most complex biotas on Earth. However, the evolutionary history of these regions in deep time is poorly understood. Here, we quantify the temporal changes in beta dissimilarities among zoogeographical regions during the Cenozoic using 4,966 extant terrestrial vertebrates and 1,278 extinct mammal genera. We identify ten present-day zoogeographical regions and find that they underwent a striking change over time. Specifically, the fauna on the TP was close to the Oriental realm in deep time but became more similar to the Palearctic realms more recently. The present-day zoogeographical regions generally emerged during the Miocene/Pliocene boundary (ca. 5 Ma). These results indicate that geological events such as the Indo-Asian Collision, the TP uplift, and the aridification of the Asian interior underpinned the evolutionary history of the zoogeographical regions surrounding the TP over different time periods.

Revisiting a Landmark Study System: No Evidence for a Punctuated Mode of Evolution in Metrarabdotos

Is speciation generally a "special time" in morphological evolution, or are lineage-splitting events just "more of the same" where the end product happens to be two separate lineages? Data on evolutionary dynamics during anagenetic and cladogenetic events among closely related lineages within a clade are rare, but the fossil record of the bryozoan genus Metrarabdotos is considered a textbook example of a clade where speciation causes rapid evolutionary change against a backdrop of morphological stasis within lineages. Here, we point to some methodological and measurement theoretical issues in the original work on Metrarabdotos. We then reanalyze a subset of the original data that can be meaningfully investigated using quantitative statistical approaches similar to those used in the original studies. We consistently fail to find variation in the evolutionary process during within-lineage evolution compared with cladogenetic events: the rates of evolution, the strength of selection, and the directions traveled in multivariate morphospace are not different when comparing evolution within lineages and at speciation events in Metrarabdotos, and genetic drift cannot be excluded as a sufficient explanation for the morphological differentiation within lineages and during speciation. Although widely considered the best example of a punctuated mode of evolution, morphological divergence and speciation are not linked in Metrarabdotos.

A genome-skimmed phylogeny of a widespread bryozoan family, Adeonidae

BACKGROUND

Understanding the phylogenetic relationships among species is one of the main goals of systematic biology. Simultaneously, credible phylogenetic hypotheses are often the first requirement for unveiling the evolutionary history of traits and for modelling macroevolutionary processes. However, many non-model taxa have not yet been sequenced to an extent such that statistically well-supported molecular phylogenies can be constructed for these purposes. Here, we use a genome-skimming approach to extract sequence information for 15 mitochondrial and 2 ribosomal operon genes from the cheilostome bryozoan family, Adeonidae, Busk, 1884, whose current systematics is based purely on morphological traits. The members of the Adeonidae are, like all cheilostome bryozoans, benthic, colonial, marine organisms. Adeonids are also geographically widely-distributed, often locally common, and are sometimes important habitat-builders.

RESULTS

We successfully genome-skimmed 35 adeonid colonies representing 6 genera (Adeona, Adeonellopsis, Bracebridgia, Adeonella, Laminopora and Cucullipora). We also contributed 16 new, circularised mitochondrial genomes to the eight previously published for cheilostome bryozoans. Using the aforementioned mitochondrial and ribosomal genes, we inferred the relationships among these 35 samples. Contrary to some previous suggestions, the Adeonidae is a robustly supported monophyletic clade. However, the genera Adeonella and Laminopora are in need of revision: Adeonella is polyphyletic and Laminopora paraphyletically forms a clade with some Adeonella species. Additionally, we assign a sequence clustering identity using cox1 barcoding region of 99% at the species and 83% at the genus level.

CONCLUSIONS

We provide sequence data, obtained via genome-skimming, that greatly increases the resolution of the phylogenetic relationships within the adeonids. We present a highly-supported topology based on 17 genes and substantially increase availability of circularised cheilostome mitochondrial genomes, and highlight how we can extend our pipeline to other bryozoans.

On trends and patterns in macroevolution: Williston's law and the branchiostegal series of extant and extinct osteichthyans

BACKGROUND

The branchiostegal series consists of an alignment of bony elements in the posterior portion of the skull of osteichthyan vertebrates. We trace the evolution of the number of elements in a comprehensive survey that includes 440 extant and 66 extinct species. Using a newly updated actinopterygian tree in combination with phylogenetic comparative analyses, we test whether osteichthyan branchiostegals follow an evolutionary trend under 'Williston's law', which postulates that osteichthyan lineages experienced a reduction of bony elements over time.

RESULTS

We detected no overall macroevolutionary trend in branchiostegal numbers, providing no support for 'Williston's law'. This result is robust to the subsampling of palaeontological data, but the estimation of the model parameters is much more ambiguous.

CONCLUSIONS

We find substantial evidence for a macroevolutionary dynamic favouring an 'early burst' of trait evolution over alternative models. Our study highlights the challenges of accurately reconstructing macroevolutionary dynamics even with large amounts of data about extant and extinct taxa.

Embracing heterogeneity: coalescing the Tree of Life and the future of phylogenomics

Building the Tree of Life (ToL) is a major challenge of modern biology, requiring advances in cyberinfrastructure, data collection, theory, and more. Here, we argue that phylogenomics stands to benefit by embracing the many heterogeneous genomic signals emerging from the first decade of large-scale phylogenetic analysis spawned by high-throughput sequencing (HTS). Such signals include those most commonly encountered in phylogenomic datasets, such as incomplete lineage sorting, but also those reticulate processes emerging with greater frequency, such as recombination and introgression. Here we focus specifically on how phylogenetic methods can accommodate the heterogeneity incurred by such population genetic processes; we do not discuss phylogenetic methods that ignore such processes, such as concatenation or supermatrix approaches or supertrees. We suggest that methods of data acquisition and the types of markers used in phylogenomics will remain restricted until a posteriori methods of marker choice are made possible with routine whole-genome sequencing of taxa of interest. We discuss limitations and potential extensions of a model supporting innovation in phylogenomics today, the multispecies coalescent model (MSC). Macroevolutionary models that use phylogenies, such as character mapping, often ignore the heterogeneity on which building phylogenies increasingly rely and suggest that assimilating such heterogeneity is an important goal moving forward. Finally, we argue that an integrative cyberinfrastructure linking all steps of the process of building the ToL, from specimen acquisition in the field to publication and tracking of phylogenomic data, as well as a culture that values contributors at each step, are essential for progress.

Inferring Diversification Rate Variation From Phylogenies With Fossils

Time-calibrated phylogenies of living species have been widely used to study the tempo and mode of species diversification. However, it is increasingly clear that inferences about species diversification-extinction rates in particular-can be unreliable in the absence of paleontological data. We introduce a general framework based on the fossilized birth-death process for studying speciation-extinction dynamics on phylogenies of extant and extinct species. The model assumes that phylogenies can be modeled as a mixture of distinct evolutionary rate regimes and that a hierarchical Poisson process governs the number of such rate regimes across a tree. We implemented the model in BAMM, a computational framework that uses reversible jump Markov chain Monte Carlo to simulate a posterior distribution of macroevolutionary rate regimes conditional on the branching times and topology of a phylogeny. The implementation, we describe can be applied to paleontological phylogenies, neontological phylogenies, and to phylogenies that include both extant and extinct taxa. We evaluate performance of the model on data sets simulated under a range of diversification scenarios. We find that speciation rates are reliably inferred in the absence of paleontological data. However, the inclusion of fossil observations substantially increases the accuracy of extinction rate estimates. We demonstrate that inferences are relatively robust to at least some violations of model assumptions, including heterogeneity in preservation rates and misspecification of the number of occurrences in paleontological data sets.

Clade-specific diversification dynamics of marine diatoms since the Jurassic

Diatoms are one of the most abundant and diverse groups of phytoplankton and play a major role in marine ecosystems and the Earth's biogeochemical cycles. Here we combine DNA metabarcoding data from the Tara Oceans expedition with palaeoenvironmental data and phylogenetic models of diversification to analyse the diversity dynamics of marine diatoms. We reveal a primary effect of variation in carbon dioxide partial pressure (pCO2) on early diatom diversification, followed by a major burst of diversification in the late Eocene epoch, after which diversification is chiefly affected by sea level, an influx of silica availability and competition with other planktonic groups. Our results demonstrate a remarkable heterogeneity of diversification dynamics across diatoms and suggest that a changing climate will favour some clades at the expense of others.

Key questions and challenges in angiosperm macroevolution

Contents Summary 1170 I. Introduction 1170 II. Six key questions 1172 III. Three key challenges 1177 IV. Conclusions 1181 Acknowledgements 1182 References 1183 SUMMARY: The origin and rapid diversification of angiosperms (flowering plants) represent one of the most intriguing topics in evolutionary biology. Despite considerable progress made in complementary fields over the last two decades (paleobotany, phylogenetics, ecology, evo-devo, genomics), many important questions remain. For instance, what has been the impact of mass extinctions on angiosperm diversification? Are the angiosperms an adaptive radiation? Has morphological evolution in angiosperms been gradual or pulsed? We propose that the recent and ongoing revolution in macroevolutionary methods provides an unprecedented opportunity to explore long-standing questions that probably hold important clues to understand present-day biodiversity. We present six key questions that explore the origin and diversification of angiosperms. We also identify three key challenges to address these questions: (1) the development of new integrative models that include diversification, multiple intrinsic and environmental traits, biogeography and the fossil record all at once, whilst accounting for sampling bias and heterogeneity of macroevolutionary processes through time and among lineages; (2) the need for large and standardized synthetic databases of morphological variation; and (3) continuous effort on sampling the fossil record, but with a revolution in current paleobotanical practice.

Rethinking Living Fossils

Biologists would be mistaken if they relegated living fossils to paleontological inquiry or assumed that the concept is dead. It is now used to describe entities ranging from viruses to higher taxa, despite recent warnings of misleading inferences. Current work on character evolution illustrates how analyzing living fossils and stasis in terms of parts (characters) and wholes (e.g., organisms and lineages) advances our understanding of prolonged stasis at many hierarchical levels. Instead of viewing the concept's task as categorizing living fossils, we show how its primary role is to mark out what is in need of explanation, accounting for the persistence of both molecular and morphological traits. Rethinking different conceptions of living fossils as specific hypotheses reveals novel avenues for research that integrate phylogenetics, ecological and evolutionary modeling, and evo-devo to produce a more unified theoretical outlook.

Automated Integration of Trees and Traits: A Case Study Using Paired Fin Loss Across Teleost Fishes

Data synthesis required for large-scale macroevolutionary studies is challenging with the current tools available for integration. Using a classic question regarding the frequency of paired fin loss in teleost fishes as a case study, we sought to create automated methods to facilitate the integration of broad-scale trait data with a sizable species-level phylogeny. Similar to the evolutionary pattern previously described for limbs, pelvic and pectoral fin reduction and loss are thought to have occurred independently multiple times in the evolution of fishes. We developed a bioinformatics pipeline to identify the presence and absence of pectoral and pelvic fins of 12,582 species. To do this, we integrated a synthetic morphological supermatrix of phenotypic data for the pectoral and pelvic fins for teleost fishes from the Phenoscape Knowledgebase (two presence/absence characters for 3047 taxa) with a species-level tree for teleost fishes from the Open Tree of Life project (38,419 species). The integration method detailed herein harnessed a new combined approach by utilizing data based on ontological inference, as well as phylogenetic propagation, to reduce overall data loss. Using inference enabled by ontology-based annotations, missing data were reduced from 98.0% to 85.9%, and further reduced to 34.8% by phylogenetic data propagation. These methods allowed us to extend the data to an additional 11,293 species for a total of 12,582 species with trait data. The pectoral fin appears to have been independently lost in a minimum of 19 lineages and the pelvic fin in 48. Though interpretation is limited by lack of phylogenetic resolution at the species level, it appears that following loss, both pectoral and pelvic fins were regained several (3) to many (14) times respectively. Focused investigation into putative regains of the pectoral fin, all within one clade (Anguilliformes), showed that the pectoral fin was regained at least twice following loss. Overall, this study points to specific teleost clades where strategic phylogenetic resolution and genetic investigation will be necessary to understand the pattern and frequency of pectoral fin reversals.

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

Integrative evolutionary analyses based upon fossil and extant species provide a powerful approach for understanding past diversification events and for assessing the tempo of evolution across the Tree of Life. Herein, we demonstrate the importance of integrating fossil and extant species for inferring patterns of lineage diversification that would otherwise be masked in analyses that examine only one source of evidence. We infer the phylogeny and macroevolutionary history of the Tetraodontiformes (triggerfishes, pufferfishes and allies), a group with one of the most extensive fossil records among fishes. Our analyses combine molecular and morphological data, based on an expanded matrix that adds newly coded fossil species and character states. Beyond confidently resolving the relationships and divergence times of tetraodontiforms, our diversification analyses detect a major mass-extinction event during the Palaeocene-Eocene Thermal Maximum (PETM), followed by a marked increase in speciation rates. This pattern is consistently obtained when fossil and extant species are integrated, whereas examination of the fossil occurrences alone failed to detect major diversification changes during the PETM. When taking into account non-homogeneous models, our analyses also detect a rapid lineage diversification increase in one of the groups (tetraodontoids) during the middle Miocene, which is considered a key period in the evolution of reef fishes associated with trophic changes and ecological opportunity. In summary, our analyses show distinct diversification dynamics estimated from phylogenies and the fossil record, suggesting that different episodes shaped the evolution of tetraodontiforms during the Cenozoic.

Development, Trait Evolution, and the Evolution of Development in Trilobites

Trilobites offer one of the best fossil records of any arthropod group. This is due to a number of factors, most notably the combination of (1) having inhabited areas where organisms are more likely to be buried and ultimately fossilized; and (2) having had a highly biomineralized exoskeleton more likely to survive the stresses of fossilization. This biomineralized exoskeleton was also morphologically complex, bearing traits that had ecological significance, and was present throughout postembryonic development, from larval to adult stages. Because the morphology of the exoskeleton changed gradually across molts during development, it is possible to reconstruct ontogenetic series for many species. Over the last decade, studies have documented both variation in modularity among closely related species and conserved developmental patterns among modules. In the latter case, trait evolution could still occur through modification of rates of morphological change along otherwise conserved ontogenetic trajectories. At the clade level, the pattern of expression and release of new exoskeletal segments during post-embryonic development was generally conserved across most species, but the relative timing of different segmentation events could vary, and developmental traits appear to have been relatively labile across the clade's evolutionary history. Most recently, comparative analyses indicate that the association between segmentation events and the timing of shifts in the rate of ontogenetic shape change varies across species. Despite these advances, we still know relatively little about how development constrained or contributed to trait evolution in trilobites, and almost nothing about the origin of novel traits in trilobites. A major (but removable) obstacle is the current lack of well-supported trilobite phylogenies that span higher taxonomic levels.

The fossilized birth-death model for the analysis of stratigraphic range data under different speciation modes

A birth-death-sampling model gives rise to phylogenetic trees with samples from the past and the present. Interpreting "birth" as branching speciation, "death" as extinction, and "sampling" as fossil preservation and recovery, this model - also referred to as the fossilized birth-death (FBD) model - gives rise to phylogenetic trees on extant and fossil samples. The model has been mathematically analyzed and successfully applied to a range of datasets on different taxonomic levels, such as penguins, plants, and insects. However, the current mathematical treatment of this model does not allow for a group of temporally distinct fossil specimens to be assigned to the same species. In this paper, we provide a general mathematical FBD modeling framework that explicitly takes "stratigraphic ranges" into account, with a stratigraphic range being defined as the lineage interval associated with a single species, ranging through time from the first to the last fossil appearance of the species. To assign a sequence of fossil samples in the phylogenetic tree to the same species, i.e., to specify a stratigraphic range, we need to define the mode of speciation. We provide expressions to account for three common speciation modes: budding (or asymmetric) speciation, bifurcating (or symmetric) speciation, and anagenetic speciation. Our equations allow for flexible joint Bayesian analysis of paleontological and neontological data. Furthermore, our framework is directly applicable to epidemiology, where a stratigraphic range is the observed duration of infection of a single patient, "birth" via budding is transmission, "death" is recovery, and "sampling" is sequencing the pathogen of a patient. Thus, we present a model that allows for incorporation of multiple observations through time from a single patient.

Combined Analysis of Extant Rhynchonellida (Brachiopoda) using Morphological and Molecular Data

Independent molecular and morphological phylogenetic analyses have often produced discordant results for certain groups which, for fossil-rich groups, raises the possibility that morphological data might mislead in those groups for which we depend upon morphology the most. Rhynchonellide brachiopods, with more than 500 extinct genera but only 19 extant genera represented today, provide an opportunity to explore the factors that produce contentious phylogenetic signal across datasets, as previous phylogenetic hypotheses generated from molecular sequence data bear little agreement with those constructed using morphological characters. Using a revised matrix of 66 morphological characters, and published ribosomal DNA sequences, we performed a series of combined phylogenetic analyses to identify conflicting phylogenetic signals. We completed a series of parsimony-based and Bayesian analyses, varying the data used, the taxa included, and the models used in the Bayesian analyses. We also performed simulation-based sensitivity analyses to assess whether the small size of the morphological data partition relative to the molecular data influenced the results of the combined analyses. In order to compare and contrast a large number of phylogenetic analyses and their resulting summary trees, we developed a measure for the incongruence between two topologies and simultaneously ignore any differences in phylogenetic resolution. Phylogenetic hypotheses generated using only morphological characters differed among each other, and with previous analyses, whereas molecular-only and combined Bayesian analyses produced extremely similar topologies. Characters historically associated with traditional classification in the Rhynchonellida have very low consistency indices on the topology preferred by the combined Bayesian analyses. Overall, this casts doubt on the use of morphological systematics to resolve relationships among the crown rhynchonellide brachiopods. However, expanding our dataset to a larger number of extinct taxa with intermediate morphologies is necessary to exclude the possibility that the morphology of extant taxa is not dominated by convergence along long branches.

Approaches to Macroevolution: 2. Sorting of Variation, Some Overarching Issues, and General Conclusions

Approaches to macroevolution require integration of its two fundamental components, within a hierarchical framework. Following a companion paper on the origin of variation, I here discuss sorting within an evolutionary hierarchy. Species sorting-sometimes termed species selection in the broad sense, meaning differential origination and extinction owing to intrinsic biological properties-can be split into strict-sense species selection, in which rate differentials are governed by emergent, species-level traits such as geographic range size, and effect macroevolution, in which rates are governed by organism-level traits such as body size; both processes can create hitchhiking effects, indirectly causing the proliferation or decline of other traits. Several methods can operationalize the concept of emergence, so that rigorous separation of these processes is increasingly feasible. A macroevolutionary tradeoff, underlain by the intrinsic traits that influence evolutionary dynamics, causes speciation and extinction rates to covary in many clades, resulting in evolutionary volatility of some clades and more subdued behavior of others; the few clades that break the tradeoff can achieve especially prolific diversification. In addition to intrinsic biological traits at multiple levels, extrinsic events can drive the waxing and waning of clades, and the interaction of traits and events are difficult but important to disentangle. Evolutionary trends can arise in many ways, and at any hierarchical level; descriptive models can be fitted to clade trajectories in phenotypic or functional spaces, but they may not be diagnostic regarding processes, and close attention must be paid to both leading and trailing edges of apparent trends. Biotic interactions can have negative or positive effects on taxonomic diversity within a clade, but cannot be readily extrapolated from the nature of such interactions at the organismic level. The relationships among macroevolutionary currencies through time (taxonomic richness, morphologic disparity, functional variety) are crucial for understanding the nature of evolutionary diversification. A novel approach to diversity-disparity analysis shows that taxonomic diversifications can lag behind, occur in concert with, or precede, increases in disparity. Some overarching issues relating to both the origin and sorting of clades and phenotypes include the macroevolutionary role of mass extinctions, the potential differences between plant and animal macroevolution, whether macroevolutionary processes have changed through geologic time, and the growing human impact on present-day macroevolution. Many challenges remain, but progress is being made on two of the key ones: (a) the integration of variation-generating mechanisms and the multilevel sorting processes that act on that variation, and (b) the integration of paleontological and neontological approaches to historical biology.

Phylogenetic classification of bony fishes

BACKGROUND

Fish classifications, as those of most other taxonomic groups, are being transformed drastically as new molecular phylogenies provide support for natural groups that were unanticipated by previous studies. A brief review of the main criteria used by ichthyologists to define their classifications during the last 50 years, however, reveals slow progress towards using an explicit phylogenetic framework. Instead, the trend has been to rely, in varying degrees, on deep-rooted anatomical concepts and authority, often mixing taxa with explicit phylogenetic support with arbitrary groupings. Two leading sources in ichthyology frequently used for fish classifications (JS Nelson's volumes of Fishes of the World and W. Eschmeyer's Catalog of Fishes) fail to adopt a global phylogenetic framework despite much recent progress made towards the resolution of the fish Tree of Life. The first explicit phylogenetic classification of bony fishes was published in 2013, based on a comprehensive molecular phylogeny ( www.deepfin.org ). We here update the first version of that classification by incorporating the most recent phylogenetic results.

RESULTS

The updated classification presented here is based on phylogenies inferred using molecular and genomic data for nearly 2000 fishes. A total of 72 orders (and 79 suborders) are recognized in this version, compared with 66 orders in version 1. The phylogeny resolves placement of 410 families, or ~80% of the total of 514 families of bony fishes currently recognized. The ordinal status of 30 percomorph families included in this study, however, remains uncertain (incertae sedis in the series Carangaria, Ovalentaria, or Eupercaria). Comments to support taxonomic decisions and comparisons with conflicting taxonomic groups proposed by others are presented. We also highlight cases were morphological support exist for the groups being classified.

CONCLUSIONS

This version of the phylogenetic classification of bony fishes is substantially improved, providing resolution for more taxa than previous versions, based on more densely sampled phylogenetic trees. The classification presented in this study represents, unlike any other, the most up-to-date hypothesis of the Tree of Life of fishes.

Approaches to Macroevolution: 1. General Concepts and Origin of Variation

Approaches to macroevolution require integration of its two fundamental components, i.e. the origin and the sorting of variation, in a hierarchical framework. Macroevolution occurs in multiple currencies that are only loosely correlated, notably taxonomic diversity, morphological disparity, and functional variety. The origin of variation within this conceptual framework is increasingly understood in developmental terms, with the semi-hierarchical structure of gene regulatory networks (GRNs, used here in a broad sense incorporating not just the genetic circuitry per se but the factors controlling the timing and location of gene expression and repression), the non-linear relation between magnitude of genetic change and the phenotypic results, the evolutionary potential of co-opting existing GRNs, and developmental responsiveness to nongenetic signals (i.e. epigenetics and plasticity), all requiring modification of standard microevolutionary models, and rendering difficult any simple definition of evolutionary novelty. The developmental factors underlying macroevolution create anisotropic probabilities-i.e., an uneven density distribution-of evolutionary change around any given phenotypic starting point, and the potential for coordinated changes among traits that can accommodate change via epigenetic mechanisms. From this standpoint, "punctuated equilibrium" and "phyletic gradualism" simply represent two cells in a matrix of evolutionary models of phenotypic change, and the origin of trends and evolutionary novelty are not simply functions of ecological opportunity. Over long timescales, contingency becomes especially important, and can be viewed in terms of macroevolutionary lags (the temporal separation between the origin of a trait or clade and subsequent diversification); such lags can arise by several mechanisms: as geological or phylogenetic artifacts, or when diversifications require synergistic interactions among traits, or between traits and external events. The temporal and spatial patterns of the origins of evolutionary novelties are a challenge to macroevolutionary theory; individual events can be described retrospectively, but a general model relating development, genetics, and ecology is needed. An accompanying paper (Jablonski in Evol Biol 2017) reviews diversity dynamics and the sorting of variation, with some general conclusions.