Macroevolutionary Dynamics Research Articles

Gut Microbiota as a Trigger of Accelerated Directional Adaptive Evolution: Acquisition of Herbivory in the Context of Extracellular Vesicles, MicroRNAs and Inter-Kingdom Crosstalk.

According to a traditional view, the specific diet in vertebrates is one of the key factors structuring the composition of the gut microbiota. In this interpretation, the microbiota assumes a subordinate position, where the larger host shapes, through evolution and its fitness, the taxonomical composition of the hosted microbiota. The present contribution shows how the evolution of herbivory, framed within the new concept of holobiont, the possibility of inter-kingdom crosstalk and its epigenetic effects, could pave the way to a completely reversed interpretation: instead of being passively shaped, the microbiota can mold and shape the general host body structure to increase its fitness. Central elements to consider in this context are the inter-kingdom crosstalk, the possibility of transporting RNAs through nanovesicles in feces from parents to offspring, and the activation of epigenetic processes passed on vertically from generation to generation. The new hypothesis is that the gut microbiota could play a great role in the macroevolutionary dynamics of herbivorous vertebrates, causing directly through host-microbiota dialog of epigenetic nature (i.e., methylation, histone acetylation, etc.), major changes in the organisms phenotype. The vertical exchange of the same microbial communities from parents to offspring, the interaction of these microbes with fairly uniform genotypes, and the socially restricted groups where these processes take place, could all explain the reasons why herbivory has appeared several time (and independently) during the evolution of vertebrates. The new interpretation could also represent a key factor in understanding the convergent evolution of analogous body structures in very distant lineages.

Mass turnover and recovery dynamics of a diverse Australian continental radiation.

Trends in global and local climate history have been linked to observed macroevolutionary patterns across a variety of organisms. These climatic pressures may unilaterally or asymmetrically influence the evolutionary trajectory of clades. To test and compare signatures of changing global (Eocene-Oligocene boundary cooling) and continental (Miocene aridification) environments on a continental fauna, we investigated the macroevolutionary dynamics of one of Australia's most diverse endemic radiations, pygopodoid geckos. We generated a time-calibrated phylogeny (>90% taxon coverage) to test whether (i) asymmetrical pygopodoid tree shape may be the result of mass turnover deep in the group's history, and (ii) how Miocene aridification shaped trends in biome assemblages. We find evidence of mass turnover in pygopodoids following the isolation of the Australian continental plate ∼30 million years ago, and in contrast, gradual aridification is linked to elevated speciation rates in the young arid zone. Surprisingly, our results suggest that invasion of arid habitats was not an evolutionary end point. Instead, arid Australia has acted as a source for diversity, with repeated outward dispersals having facilitated diversification of this group. This pattern contrasts trends in richness and distribution of other Australian vertebrates, illustrating the profound effects historical biome changes have on macroevolutionary patterns.

Discrete and morphometric traits reveal contrasting patterns and processes in the macroevolutionary history of a clade of scorpions.

Many palaeontological studies have investigated the evolution of entire body plans, generally relying on discrete character-taxon matrices. In contrast, macroevolutionary studies performed by neontologists have mostly focused on morphometric traits. Although these data types are very different, some studies have suggested that they capture common patterns. Nonetheless, the tests employed to support this claim have not explicitly incorporated a phylogenetic framework and may therefore be susceptible to confounding effects due to the presence of common phylogenetic structure. We address this question using the scorpion genus Brachistosternus Pocock 1893 as case study. We make use of a time-calibrated multilocus molecular phylogeny, and compile discrete and traditional morphometric data sets, both capturing the overall morphology of the organisms. We find that morphospaces derived from these matrices are significantly different, and that the degree of discordance cannot be replicated by simulations of random character evolution. Moreover, we find strong support for contrasting modes of evolution, with discrete characters being congruent with an 'early burst' scenario whereas morphometric traits suggest species-specific adaptations to have driven morphological evolution. The inferred macroevolutionary dynamics are therefore contingent on the choice of character type. Finally, we confirm that metrics of correlation fail to detect these profound differences given common phylogenetic structure in both data sets, and that methods incorporating a phylogenetic framework and accounting for expected covariance should be favoured.

Is BAMM Flawed? Theoretical and Practical Concerns in the Analysis of Multi-Rate Diversification Models.

Bayesian analysis of macroevolutionary mixtures (BAMM) is a statistical framework that uses reversible jump Markov chain Monte Carlo to infer complex macroevolutionary dynamics of diversification and phenotypic evolution on phylogenetic trees. A recent article by Moore et al. (MEA) reported a number of theoretical and practical concerns with BAMM. Major claims from MEA are that (i) BAMM’s likelihood function is incorrect, because it does not account for unobserved rate shifts; (ii) the posterior distribution on the number of rate shifts is overly sensitive to the prior; and (iii) diversification rate estimates from BAMM are unreliable. Here, we show that these and other conclusions from MEA are generally incorrect or unjustified. We first demonstrate that MEA’s numerical assessment of the BAMM likelihood is compromised by their use of an invalid likelihood function. We then show that “unobserved rate shifts” appear to be irrelevant for biologically plausible parameterizations of the diversification process. We find that the purportedly extreme prior sensitivity reported by MEA cannot be replicated with standard usage of BAMM v2.5, or with any other version when conventional Bayesian model selection is performed. Finally, we demonstrate that BAMM performs very well at estimating diversification rate variation across the n}{}{sim}20% of simulated trees in MEA’s data set for which it is theoretically possible to infer rate shifts with confidence. Due to ascertainment bias, the remaining 80% of their purportedly variable-rate phylogenies are statistically indistinguishable from those produced by a constant-rate birth–death process and were thus poorly suited for the summary statistics used in their performance assessment. We demonstrate that inferences about diversification rates have been accurate and consistent across all major previous releases of the BAMM software. We recognize an acute need to address the theoretical foundations of rate-shift models for phylogenetic trees, and we expect BAMM and other modeling frameworks to improve in response to mathematical and computational innovations. However, we remain optimistic that that the imperfect tools currently available to comparative biologists have provided and will continue to provide important insights into the diversification of life on Earth.

Macroevolutionary synthesis of flowering plant sexual systems.

Sexual system is a key determinant of genetic variation and reproductive success, affecting evolution within populations and within clades. Much research in plants has focused on evolutionary transitions away from the most common state of hermaphroditism and toward the rare state of dioecy (separate sexes). Rather than transitions predominantly toward greater sexual differentiation, however, evolution may proceed in the direction of lesser sexual differentiation. We analyzed the macroevolutionary dynamics of sexual system in angiosperm genera that contain both dioecious and nondioecious species. Our phylogenetic analyses encompass a total of 2145 species from 40 genera. Overall, we found little evidence that rates of sexual system transitions are greater in any direction. Counting the number of inferred state changes revealed a mild prevalence of transitions away from hermaphroditism and away from dioecy, toward states of intermediate sexual differentiation. We identify genera in which future studies of sexual system evolution might be especially productive, and we discuss how integrating genetic or population-level studies of sexual system could improve the power of phylogenetic comparative analyses. Our work adds to the evidence that different selective pressures and constraints act in different groups, helping maintain the variety of sexual systems observed among plants.

Bayesian estimation of fossil phylogenies and the evolution of early to middle Paleozoic crinoids (Echinodermata)

AbstractKnowledge of phylogenetic relationships among species is fundamental to understanding basic patterns in evolution and underpins nearly all research programs in biology and paleontology. However, most methods of phylogenetic inference typically used by paleontologists do not accommodate the idiosyncrasies of fossil data and therefore do not take full advantage of the information provided by the fossil record. The advent of Bayesian ‘tip-dating’ approaches to phylogeny estimation is especially promising for paleosystematists because time-stamped comparative data can be combined with probabilistic models tailored to accommodate the study of fossil taxa. Under a Bayesian framework, the recently developed fossilized birth–death (FBD) process provides a more realistic tree prior model for paleontological data that accounts for macroevolutionary dynamics, preservation, and sampling when inferring phylogenetic trees containing fossils. In addition, the FBD tree prior allows for the possibility of sampling ancestral morphotaxa. Although paleontologists are increasingly embracing probabilistic phylogenetic methods, these recent developments have not previously been applied to the deep-time invertebrate fossil record. Here, I examine phylogenetic relationships among Ordovician through Devonian crinoids using a Bayesian tip-dating approach. Results support several clades recognized in previous analyses sampling only Ordovician taxa, but also reveal instances where phylogenetic affinities are more complex and extensive revisions are necessary, particularly among the Cladida. The name Porocrinoidea is proposed for a well-supported clade of Ordovician ‘cyathocrine’ cladids and hybocrinids. The Eucladida is proposed as a clade name for the sister group of the Flexibilia herein comprised of cladids variously considered ‘cyathocrines,’ ‘dendrocrines,’ and/or ‘poteriocrines’ by other authors.

Multispecies coalescent delimits structure, not species

The multispecies coalescent model underlies many approaches used for species delimitation. In previous work assessing the performance of species delimitation under this model, speciation was treated as an instantaneous event rather than as an extended process involving distinct phases of speciation initiation (structuring) and completion. Here, we use data under simulations that explicitly model speciation as an extended process rather than an instantaneous event and carry out species delimitation inference on these data under the multispecies coalescent. We show that the multispecies coalescent diagnoses genetic structure, not species, and that it does not statistically distinguish structure associated with population isolation vs. species boundaries. Because of the misidentification of population structure as putative species, our work raises questions about the practice of genome-based species discovery, with cascading consequences in other fields. Specifically, all fields that rely on species as units of analysis, from conservation biology to studies of macroevolutionary dynamics, will be impacted by inflated estimates of the number of species, especially as genomic resources provide unprecedented power for detecting increasingly finer-scaled genetic structure under the multispecies coalescent. As such, our work also represents a general call for systematic study to reconsider a reliance on genomic data alone. Until new methods are developed that can discriminate between structure due to population-level processes and that due to species boundaries, genomic-based results should only be considered a hypothesis that requires validation of delimited species with multiple data types, such as phenotypic and ecological information.

Diversity and Coevolutionary Dynamics in High-Dimensional Phenotype Spaces

We study macroevolutionary dynamics by extending microevolutionary competition models to long timescales. It has been shown that for a general class of competition models, gradual evolutionary change in continuous phenotypes (evolutionary dynamics) can be nonstationary and even chaotic when the dimension of the phenotype space in which the evolutionary dynamics unfold is high. It has also been shown that evolutionary diversification can occur along nonequilibrium trajectories in phenotype space. We combine these lines of thinking by studying long-term coevolutionary dynamics of emerging lineages in multidimensional phenotype spaces. We use a statistical approach to investigate the evolutionary dynamics of many different systems. We find (1) that, for a given dimension of phenotype space, the coevolutionary dynamics tend to be fast and nonstationary for an intermediate number of coexisting lineages but tend to stabilize as the evolving communities reach a saturation level of diversity and (2) that the amount of diversity at the saturation level increases rapidly (exponentially) with the dimension of phenotype space. These results have implications for theoretical perspectives on major macroevolutionary patterns such as adaptive radiation, long-term temporal patterns of phenotypic changes, and the evolution of diversity.

Sea level regulated tetrapod diversity dynamics through the Jurassic/Cretaceous interval.

Reconstructing deep time trends in biodiversity remains a central goal for palaeobiologists, but our understanding of the magnitude and tempo of extinctions and radiations is confounded by uneven sampling of the fossil record. In particular, the Jurassic/Cretaceous (J/K) boundary, 145 million years ago, remains poorly understood, despite an apparent minor extinction and the radiation of numerous important clades. Here we apply a rigorous subsampling approach to a comprehensive tetrapod fossil occurrence data set to assess the group's macroevolutionary dynamics through the J/K transition. Although much of the signal is exclusively European, almost every higher tetrapod group was affected by a substantial decline across the boundary, culminating in the extinction of several important clades and the ecological release and radiation of numerous modern tetrapod groups. Variation in eustatic sea level was the primary driver of these patterns, controlling biodiversity through availability of shallow marine environments and via allopatric speciation on land.

Nematode endoparasites do not codiversify with their stick insect hosts.

Host–parasite coevolution stems from reciprocal selection on host resistance and parasite infectivity, and can generate some of the strongest selective pressures known in nature. It is widely seen as a major driver of diversification, the most extreme case being parallel speciation in hosts and their associated parasites. Here, we report on endoparasitic nematodes, most likely members of the mermithid family, infecting different Timema stick insect species throughout California. The nematodes develop in the hemolymph of their insect host and kill it upon emergence, completely impeding host reproduction. Given the direct exposure of the endoparasites to the host's immune system in the hemolymph, and the consequences of infection on host fitness, we predicted that divergence among hosts may drive parallel divergence in the endoparasites. Our phylogenetic analyses suggested the presence of two differentiated endoparasite lineages. However, independently of whether the two lineages were considered separately or jointly, we found a complete lack of codivergence between the endoparasitic nematodes and their hosts in spite of extensive genetic variation among hosts and among parasites. Instead, there was strong isolation by distance among the endoparasitic nematodes, indicating that geography plays a more important role than host‐related adaptations in driving parasite diversification in this system. The accumulating evidence for lack of codiversification between parasites and their hosts at macroevolutionary scales contrasts with the overwhelming evidence for coevolution within populations, and calls for studies linking micro‐ versus macroevolutionary dynamics in host–parasite interactions.

Environmental drivers of crocodyliform extinction across the Jurassic/Cretaceous transition

Crocodyliforms have a much richer evolutionary history than represented by their extant descendants, including several independent marine and terrestrial radiations during the Mesozoic. However, heterogeneous sampling of their fossil record has obscured their macroevolutionary dynamics, and obfuscated attempts to reconcile external drivers of these patterns. Here, we present a comprehensive analysis of crocodyliform biodiversity through the Jurassic/Cretaceous (J/K) transition using subsampling and phylogenetic approaches and apply maximum-likelihood methods to fit models of extrinsic variables to assess what mediated these patterns. A combination of fluctuations in sea-level and episodic perturbations to the carbon and sulfur cycles was primarily responsible for both a marine and non-marine crocodyliform biodiversity decline through the J/K boundary, primarily documented in Europe. This was tracked by high extinction rates at the boundary and suppressed origination rates throughout the Early Cretaceous. The diversification of Eusuchia and Notosuchia likely emanated from the easing of ecological pressure resulting from the biodiversity decline, which also culminated in the extinction of the marine thalattosuchians in the late Early Cretaceous. Through application of rigorous techniques for estimating biodiversity, our results demonstrate that it is possible to tease apart the complex array of controls on diversification patterns in major archosaur clades.

Macroevolutionary dynamics in the early diversification of Asteraceae

Spatial and temporal differences in ecological opportunity can result in disparity of net species diversification rates and consequently uneven distribution of taxon richness across the tree of life. The largest eudicotyledonous plant family Asteraceae has a global distribution and at least 460 times more species than its South American endemic sister family Calyceraceae. In this study, diversification rate dynamics across Asteraceae are examined in light of the several hypothesized causes for the family’s evolutionary success that could be responsible for rate change. The innovations of racemose capitulum and pappus, and a whole genome duplication event occurred near the origin of the family, yet we found the basal lineages of Asteraceae that evolved in South America share background diversification rates with Calyceraceae and their Australasian sister Goodeniaceae. Instead we found diversification rates increased gradually from the origin of Asteraceae approximately 69.5Ma in the late Cretaceous through the Early Eocene Climatic Optimum at least. In contrast to earlier studies, significant rate shifts were not strongly correlated with intercontinental dispersals or polyploidization. The difference is due primarily to sampling more backbone nodes, as well as calibrations placed internally in Asteraceae that resulted in earlier divergence times than those found in most previous relaxed clock studies. Two clades identified as having transformed rate processes are the Vernonioid Clade and a clade within the Heliantheae alliance characterized by phytomelanic fruit (PF Clade) that represents an American radiation. In Africa, subfamilies Carduoideae, Pertyoideae, Gymnarrhenoideae, Cichorioideae, Corymbioideae, and Asteroideae diverged in a relatively short span of only 6.5millionyears during the Middle Eocene.

A new species of Apidium (Anthropoidea, Parapithecidae) from the Sirt Basin, central Libya: First record of Oligocene primates from Libya

A new species of Apidium is the most common primate currently known from a newly discovered site near Zallah Oasis in the Sirt Basin of central Libya. Based on current knowledge of the associated fauna, this new species of Apidium is early Oligocene in age, being roughly contemporaneous with faunas from Quarries G and V in the upper part of the Jebel Qatrani Formation in Egypt that also contain species of Apidium. A phylogenetic analysis based on dental characters indicates that the new species of Apidium from Libya is the sister group of Apidium phiomense. Apidium bowni and Apidium moustafai from the Jebel Qatrani Formation in the Fayum are similar in age to the new species of Apidium from Libya, but both of these Egyptian species are more distantly related to A. phiomense from younger stratigraphic levels in the Fayum. This phylogenetic pattern underscores the benefit of enhanced geographic sampling of the fossil record, even in cases where local records are thought to be reasonably comprehensive and well documented. Oligocene parapithecids can be partitioned into two clades corresponding to the subfamilies Parapithecinae (containing Parapithecus and Simonsius) and Qatraniinae (including Qatrania and Apidium). Climatic deterioration during the early Oligocene may have impacted the macroevolutionary dynamics of early Afro-Arabian anthropoids by fostering the fragmentation of forest habitats, thereby promoting allopatric speciation among widespread populations of Apidium and other arboreal taxa.

Continental faunal exchange and the asymmetrical radiation of carnivores.

Lineages arriving on islands may undergo explosive evolutionary radiations owing to the wealth of ecological opportunities. Although studies on insular taxa have improved our understanding of macroevolutionary phenomena, we know little about the macroevolutionary dynamics of continental exchanges. Here we study the evolution of eight Carnivora families that have migrated across the Northern Hemisphere to investigate if continental invasions also result in explosive diversification dynamics. We used a Bayesian approach to estimate speciation and extinction rates from a substantial dataset of fossil occurrences while accounting for the incompleteness of the fossil record. Our analyses revealed a strongly asymmetrical pattern in which North American lineages invading Eurasia underwent explosive radiations, whereas lineages invading North America maintained uniform diversification dynamics. These invasions into Eurasia were characterized by high rates of speciation and extinction. The radiation of the arriving lineages in Eurasia coincide with the decline of established lineages or phases of climate change, suggesting differences in the ecological settings between the continents may be responsible for the disparity in diversification dynamics. These results reveal long-term outcomes of biological invasions and show that the importance of explosive radiations in shaping diversity extends beyond insular systems and have significant impact at continental scales.

A Tale of Two Hyper-diversities: Diversification dynamics of the two largest families of lichenized fungi.

Renewed interests in macroevolutionary dynamics have led to the proliferation of studies on diversification processes in large taxonomic groups, such as angiosperms, mammals, and birds. However, such a study has yet to be conducted in lichenized fungi – an extremely successful and diverse group of fungi. Analysing the most comprehensive time-calibrated phylogenies with a new analytical method, we illustrated drastically different diversification dynamics between two hyper-diverse families of lichenized fungi, Graphidaceae and Parmeliaceae, which represent more than a fourth of the total species diversity of lichenized fungi. Despite adopting a similar nutrition mode and having a similar number of species, Graphidaceae exhibited a lower speciation rate, while Parmeliaceae showed a sharp increase in speciation rate that corresponded with the aridification during the Oligocene-Miocene transition, suggesting their adaptive radiation into a novel arid habitat.

Evolutionary time-series analysis reveals the signature of frequency-dependent selection on a female mating polymorphism.

A major challenge in evolutionary biology is understanding how stochastic and deterministic factors interact and influence macroevolutionary dynamics in natural populations. One classical approach is to record frequency changes of heritable and visible genetic polymorphisms over multiple generations. Here, we combined this approach with a maximum likelihood-based population-genetic model with the aim of understanding and quantifying the evolutionary processes operating on a female mating polymorphism in the blue-tailed damselfly Ischnura elegans. Previous studies on this color-polymorphic species have suggested that males form a search image for females, which leads to excessive mating harassment of common female morphs. We analyzed a large temporally and spatially replicated data set of between-generation morph frequency changes in I. elegans. Morph frequencies were more stable than expected from genetic drift alone, suggesting the presence of selection toward a stable equilibrium that prevents local loss or fixation of morphs. This can be interpreted as the signature of negative frequency-dependent selection maintaining the phenotypic stasis and genetic diversity in these populations. Our novel analytical approach allows the estimation of the strength of frequency-dependent selection from the morph frequency fluctuations around their inferred long-term equilibria. This approach can be extended and applied to other polymorphic organisms for which time-series data across multiple generations are available.

Causal evidence between monsoon and evolution of rhizomyine rodents.

The modern Asian monsoonal systems are currently believed to have originated around the end of the Oligocene following a crucial step of uplift of the Tibetan-Himalayan highlands. Although monsoon possibly drove the evolution of many mammal lineages during the Neogene, no evidence thereof has been provided so far. We examined the evolutionary history of a clade of rodents, the Rhizomyinae, in conjunction with our current knowledge of monsoon fluctuations over time. The macroevolutionary dynamics of rhizomyines were analyzed within a well-constrained phylogenetic framework coupled with biogeographic and evolutionary rate studies. The evolutionary novelties developed by these rodents were surveyed in parallel with the fluctuations of the Indian monsoon so as to evaluate synchroneity and postulate causal relationships. We showed the existence of three drops in biodiversity during the evolution of rhizomyines, all of which reflected elevated extinction rates. Our results demonstrated linkage of monsoon variations with the evolution and biogeography of rhizomyines. Paradoxically, the evolution of rhizomyines was accelerated during the phases of weakening of the monsoons, not of strengthening, most probably because at those intervals forest habitats declined, which triggered extinction and progressive specialization toward a burrowing existence.

BAMMtools: an R package for the analysis of evolutionary dynamics on phylogenetic trees

SummaryUnderstanding the dynamics of speciation, extinction and phenotypic evolution is a central challenge in evolutionary biology. Here, we presentBAMMtools, anrpackage for the analysis and visualization of macroevolutionary dynamics on phylogenetic trees.BAMMtools is a companion package toBAMM, an open‐source program for reversible‐jumpMCMCanalyses of diversification and trait evolution.Functions inBAMMtools operate directly on output from theBAMMprogram. The package is oriented towards reconstructing and visualizing changes in evolutionary rates through time and across clades in a Bayesian statistical framework.BAMMtools enables users to extract credible sets of diversification shifts and to identify diversification histories with the maximuma posterioriprobability. Users can compare the fit of alternative diversification models using Bayes factors and by directly comparing model posterior probabilities.By providing a robust framework for quantifying uncertainty in macroevolutionary dynamics,BAMMtools will facilitate inference on the complex mixture of processes that have shaped the distribution of species and phenotypes across the tree of life.

Analysis and Visualization of Complex Macroevolutionary Dynamics: An Example from Australian Scincid Lizards

The correlation between species diversification and morphological evolution has long been of interest in evolutionary biology. We investigated the relationship between these processes during the radiation of 250+scincid lizards that constitute Australia's most species-rich clade of terrestrial vertebrates. We generated a time-calibrated phylogenetic tree for the group that was more than 85% complete at the species level and collected multivariate morphometric data for 183 species. We reconstructed the dynamics of species diversification and trait evolution using a Bayesian statistical framework (BAMM) that simultaneously accounts for variation in evolutionary rates through time and among lineages. We extended the BAMM model to accommodate time-dependent phenotypic evolution, and we describe several new methods for summarizing and visualizing macroevolutionary rate heterogeneity on phylogenetic trees. Two major clades (Lerista, Ctenotus; >90 spp. each) are associated with high rates of species diversification relative to the background rate across Australian sphenomorphine skinks. The Lerista clade is characterized by relatively high lability of body form and has undergone repeated instances of limb reduction, but Ctenotus is characterized by an extreme deceleration in the rate of body shape evolution. We estimate that rates of phenotypic evolution decreased by more than an order of magnitude in the common ancestor of the Ctenotus clade. These results provide evidence for a modal shift in phenotypic evolutionary dynamics and demonstrate that major axes of morphological variation can be decoupled from species diversification. More generally, the Bayesian framework described here can be used to identify and characterize complex mixtures of dynamic processes on phylogenetic trees. [Bayesian; diversification; evolvability; lizard; macroevolution, punctuated equilibrium, speciation.].

Automatic Detection of Key Innovations, Rate Shifts, and Diversity-Dependence on Phylogenetic Trees

A number of methods have been developed to infer differential rates of species diversification through time and among clades using time-calibrated phylogenetic trees. However, we lack a general framework that can delineate and quantify heterogeneous mixtures of dynamic processes within single phylogenies. I developed a method that can identify arbitrary numbers of time-varying diversification processes on phylogenies without specifying their locations in advance. The method uses reversible-jump Markov Chain Monte Carlo to move between model subspaces that vary in the number of distinct diversification regimes. The model assumes that changes in evolutionary regimes occur across the branches of phylogenetic trees under a compound Poisson process and explicitly accounts for rate variation through time and among lineages. Using simulated datasets, I demonstrate that the method can be used to quantify complex mixtures of time-dependent, diversity-dependent, and constant-rate diversification processes. I compared the performance of the method to the MEDUSA model of rate variation among lineages. As an empirical example, I analyzed the history of speciation and extinction during the radiation of modern whales. The method described here will greatly facilitate the exploration of macroevolutionary dynamics across large phylogenetic trees, which may have been shaped by heterogeneous mixtures of distinct evolutionary processes.