Estimating Age-Dependent Extinction: Contrasting Evidence from Fossils and Phylogenies

The Angiosperm Terrestrial Revolution buffered ants against extinction

With ~14,000 extant species, ants are ubiquitous and of tremendous ecological importance. They have undergone remarkable diversification throughout their evolutionary history. However, the drivers of their diversity dynamics are not well quantified or understood. Previous phylogenetic analyses have suggested patterns of diversity dynamics associated with the Angiosperm Terrestrial Revolution (ATR), but these studies have overlooked valuable information from the fossil record. To address this gap, we conducted a comprehensive analysis using a large dataset that includes both the ant fossil record (~24,000 individual occurrences) and neontological data (~14,000 occurrences), and tested four hypotheses proposed for ant diversification: co-diversification, competitive extinction, hyper-specialization, and buffered extinction. Taking into account biases in the fossil record, we found three distinct diversification periods (the latest Cretaceous, Eocene, and Oligo-Miocene) and one extinction period (Late Cretaceous). The competitive extinction hypothesis between stem and crown ants is not supported. Instead, we found support for the co-diversification, buffered extinction, and hyper-specialization hypotheses. The environmental changes of the ATR, mediated by the angiosperm radiation, likely played a critical role in buffering ants against extinction and favoring their diversification by providing new ecological niches, such as forest litter and arboreal nesting sites, and additional resources. We also hypothesize that the decline and extinction of stem ants during the Late Cretaceous was due to their hyper-specialized morphology, which limited their ability to expand their dietary niche in changing environments. This study highlights the importance of a holistic approach when studying the interplay between past environments and the evolutionary trajectories of organisms.

How Important Is Budding Speciation for Comparative Studies?

The acknowledgment of evolutionary dependence among species has fundamentally changed how we ask biological questions. Phylogenetic models became the standard approach for studies with 3 or more lineages, in particular those using extant species. Most phylogenetic comparative methods (PCMs) translate relatedness into covariance, meaning that evolutionary changes before lineages split should be interpreted together whereas after the split lineages are expected to change independently. This clever realization has shaped decades of research. Here, we discuss one element of the comparative method often ignored or assumed as unimportant: if nodes of a phylogeny represent the dissolution of the ancestral lineage into two new ones or if the ancestral lineage can survive speciation events (i.e., budding). Budding speciation is often reported in paleontological studies, due to the nature of the evidence for budding in the fossil record, but it is surprisingly absent in comparative methods. Here, we show that many PCMs assume that divergence happens as a symmetric split, even if these methods do not explicitly mention this assumption. We discuss the properties of trait evolution models for continuous and discrete traits and their adequacy under a scenario of budding speciation. We discuss the effects of budding speciation under a series of plausible evolutionary scenarios and show when and how these can influence our estimates. We also propose that long-lived lineages that have survived through a series of budding speciation events and given birth to multiple new lineages can produce evolutionary patterns that challenge our intuition about the most parsimonious history of trait changes in a clade. We hope our discussion can help bridge comparative approaches in paleontology and neontology as well as foster awareness about the assumptions we make when we use phylogenetic trees.

Drivers of the artiodactyl turnover in insular western Europe at the Eocene-Oligocene Transition

Simultaneously investigating the effects of abiotic and biotic factors on diversity dynamics is essential to understand the evolutionary history of clades. The Grande Coupure corresponds to a major faunal turnover at the Eocene-Oligocene transition (EOT) (~34.1 to 33.55 Mya) and is defined in western Europe as an extinction of insular European mammals coupled with the arrival of crown clades from Asia. Here, we focused on the species-rich group of endemic European artiodactyls to determine the drivers of the Grande Coupure during the major environmental disruptions at the EOT. Using Bayesian birth-death models, we analyzed an original high-resolution fossil dataset (90 species, >2,100 occurrences) from southwestern France (Quercy area) and estimated the regional diversification and diversity dynamics of endemic and immigrant artiodactyls. We show that the endemic artiodactyl radiation was mainly related to the Eocene tropical conditions, combined with biotic controls on speciation and clade-related diversity dependence. We further highlight that the major environmental changes at the transition (77% of species became extinct) and the concurrent increase in seasonality in Europe during the Oligocene were likely the main drivers of their decline. Surprisingly, our results do not support the widely-held hypothesis of active competition between endemic and immigrant artiodactyls but rather suggest a passive or opportunistic replacement by immigrants, which is further supported by morphological clustering of specific ecological traits across the Eocene-Oligocene transition. Our analyses provide insights into the evolutionary and ecological processes driving the diversification and decline of mammalian clades during a major biological and climatic crisis.

Global impact and selectivity of the Cretaceous-Paleogene mass extinction among sharks, skates, and rays

The Cretaceous-Paleogene event was the last mass extinction event, yet its impact and long-term effects on species-level marine vertebrate diversity remain largely uncharacterized. We quantified elasmobranch (sharks, skates, and rays) speciation, extinction, and ecological change resulting from the end-Cretaceous event using >3200 fossil occurrences and 675 species spanning the Late Cretaceous-Paleocene interval at global scale. Elasmobranchs declined by >62% at the Cretaceous-Paleogene boundary and did not fully recover in the Paleocene. The end-Cretaceous event triggered a heterogeneous pattern of extinction, with rays and durophagous species reaching the highest levels of extinction (>72%) and sharks and nondurophagous species being less affected. Taxa with large geographic ranges and/or those restricted to high-latitude settings show higher survival. The Cretaceous-Paleogene event drastically altered the evolutionary history of marine ecosystems.

Multiple drivers and lineage-specific insect extinctions during the Permo-Triassic

The Permo-Triassic interval encompasses three extinction events including the most dramatic biological crisis of the Phanerozoic, the latest Permian mass extinction. However, their drivers and outcomes are poorly quantified and understood for terrestrial invertebrates, which we assess here for insects. We find a pattern with three extinctions: the Roadian/Wordian (≈266.9 Ma; extinction of 64.5% insect genera), the Permian/Triassic (≈252 Ma; extinction of 82.6% insect genera), and the Ladinian/Carnian boundaries (≈237 Ma; extinction of 74.8% insect genera). We also unveil a heterogeneous effect of these extinction events across the major insect clades. Because extinction events have impacted Permo-Triassic ecosystems, we investigate the influence of abiotic and biotic factors on insect diversification dynamics and find that changes in floral assemblages are likely the strongest drivers of insects' responses throughout the Permo-Triassic. We also assess the effect of diversity dependence between three insect guilds; an effect ubiquitously found in current ecosystems. We find that herbivores held a central position in the Permo-Triassic interaction network. Our study reveals high levels of insect extinction that profoundly shaped the evolutionary history of the most diverse non-microbial lineage.

The Shape of Phylogenies Under Phase-Type Distributed Times to Speciation and Extinction

Phylogenetic trees describe relationships between extant species, but beyond that their shape and their relative branch lengths can provide information on broader evolutionary processes of speciation and extinction. However, currently many of the most widely used macro-evolutionary models make predictions about the shapes of phylogenetic trees that differ considerably from what is observed in empirical phylogenies. Here, we propose a flexible and biologically plausible macroevolutionary model for phylogenetic trees where times to speciation or extinction events are drawn from a Coxian phase-type (PH) distribution. First, we show that different choices of parameters in our model lead to a range of tree balances as measured by Aldous’ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β statistic. In particular, we demonstrate that it is possible to find parameters that correspond well to empirical tree balance. Next, we provide a natural extension of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β statistic to sets of trees. This extension produces less biased estimates of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β compared to using the median \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β values from individual trees. Furthermore, we derive a likelihood expression for the probability of observing an edge-weighted tree under a model with speciation but no extinction. Finally, we illustrate the application of our model by performing both absolute and relative goodness-of-fit tests for two large empirical phylogenies (squamates and angiosperms) that compare models with Coxian PH distributed times to speciation with models that assume exponential or Weibull distributed waiting times. In our numerical analysis, we found that, in most cases, models assuming a Coxian PH distribution provided the best fit.

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.

Evolutionary legacies in contemporary tetrapod imperilment

The Tree of Life will be irrevocably reshaped as anthropogenic extinctions continue to unfold. Theory suggests that lineage evolutionary dynamics, such as age since origination, historical extinction filters and speciation rates, have influenced ancient extinction patterns - but whether these factors also contribute to modern extinction risk is largely unknown. We examine evolutionary legacies in contemporary extinction risk for over 4000 genera, representing ~30,000 species, from the major tetrapod groups: amphibians, birds, turtles and crocodiles, squamate reptiles and mammals. We find consistent support for the hypothesis that extinction risk is elevated in lineages with higher recent speciation rates. We subsequently test, and find modest support for, a primary mechanism driving this pattern: that rapidly diversifying clades predominantly comprise range-restricted, and extinction-prone, species. These evolutionary patterns in current imperilment may have important consequences for how we manage the erosion of biological diversity across the Tree of Life.

Dinosaur biodiversity declined well before the asteroid impact, influenced by ecological and environmental pressures

The question why non-avian dinosaurs went extinct 66 million years ago (Ma) remains unresolved because of the coarseness of the fossil record. A sudden extinction caused by an asteroid is the most accepted hypothesis but it is debated whether dinosaurs were in decline or not before the impact. We analyse the speciation-extinction dynamics for six key dinosaur families, and find a decline across dinosaurs, where diversification shifted to a declining-diversity pattern ~76 Ma. We investigate the influence of ecological and physical factors, and find that the decline of dinosaurs was likely driven by global climate cooling and herbivorous diversity drop. The latter is likely due to hadrosaurs outcompeting other herbivores. We also estimate that extinction risk is related to species age during the decline, suggesting a lack of evolutionary novelty or adaptation to changing environments. These results support an environmentally driven decline of non-avian dinosaurs well before the asteroid impact.

Evolutionary drivers, morphological evolution and diversity dynamics of a surviving mammal clade: cainotherioids at the Eocene-Oligocene transition

The Eocene-Oligocene transition (EOT) represents a period of global environmental changes particularly marked in Europe and coincides with a dramatic biotic turnover. Here, using an exceptional fossil preservation, we document and analyse the diversity dynamics of a mammal clade, Cainotherioidea (Artiodactyla), that survived the EOT and radiated rapidly immediately after. We infer their diversification history from Quercy Konzentrat-Lagerstätte (south-west France) at the species level using Bayesian birth-death models. We show that cainotherioid diversity fluctuated through time, with extinction events at the EOT and in the late Oligocene, and a major speciation burst in the early Oligocene. The latter is in line with our finding that cainotherioids had a high morphological adaptability following environmental changes throughout the EOT, which probably played a key role in the survival and evolutionary success of this clade in the aftermath. Speciation is positively associated with temperature and continental fragmentation in a time-continuous way, while extinction seems to synchronize with environmental change in a punctuated way. Within-clade interactions negatively affected the cainotherioid diversification, while inter-clade competition might explain their final decline during the late Oligocene. Our results provide a detailed dynamic picture of the evolutionary history of a mammal clade in a context of global change.

Re-evaluation of the "law of constant extinction" for ruminants at different taxonomical scales

The "law of constant extinction," proposed by Van Valen, states that long and short-lived taxa have equal chances of going extinct. This pattern of age-independent extinction was originally inferred using the fossil record of several different taxa and relied on survivorship curves built from the literal reading of the fossil record. Van Valen's seminal work was mostly done at higher taxonomic levels, hence its prevalence at the species level could not be directly inferred. The surprisingly few subsequent studies done at the species level have challenged the prevalence of age-independent extinction, but those have, for the most part, failed to explicitly incorporate inherent biases of the fossil record. Using a recent Bayesian framework that accounts for several of those biases, including the fact that very short-living lineages might never make to the record itself, we showed that Ruminantia species present age-dependent extinction, where extinction probability decreases with species age. An analysis at the genus level suggested age-independent extinction but further examination suggested that the pattern might be more complex than previously reported by Van Valen. Our results indicate that different taxonomic levels may present different extinction regimes, which could justify the development of new macroevolutionary theory and methods.

A Total-Group Phylogenetic Metatree for Cetacea and the Importance of Fossil Data in Diversification Analyses

Phylogenetic trees provide a powerful framework for testing macroevolutionary hypotheses, but it is becoming increasingly apparent that inferences derived from extant species alone can be highly misleading. Trees incorporating living and extinct taxa are are needed to address fundamental questions about the origins of diversity and disparity but it has proved challenging to generate robust, species-rich phylogenies that include large numbers of fossil taxa. As a result, most studies of diversification dynamics continue to rely on molecular phylogenies. Here, we extend and apply a recently developed meta-analytic approach for synthesizing previously published phylogenetic studies to infer a well-resolved set of species level, time-scaled phylogenetic hypotheses for extinct and extant cetaceans (whales, dolphins and allies). Our trees extend sampling from the ∼ 90 extant species to over 500 living and extinct species, and therefore allow for more robust inference of macroevolutionary dynamics. While the diversification scenarios we recover are broadly concordant with those inferred from molecular phylogenies they differ in critical ways, notably in the relative contributions of extinction and speciation rate shifts in driving rapid radiations. The metatree approach provides the most immediate route for generating higher level phylogenies of extinct taxa, and opens the door to re-evaluation of macroevolutionary hypotheses derived only from extant taxa.

Branching patterns in phylogenies cannot distinguish diversity-dependent diversification from time-dependent diversification

One of the primary goals of macroevolutionary biology has been to explain general trends in long‐term diversity patterns, including whether such patterns correspond to an upscaling of processes occurring at lower scales. Reconstructed phylogenies often show decelerated lineage accumulation over time. This pattern has often been interpreted as the result of diversity‐dependent (DD) diversification, where the accumulation of species causes diversification to decrease through niche filling. However, other processes can also produce such a slowdown, including time dependence without diversity dependence. To test whether phylogenetic branching patterns can be used to distinguish these two mechanisms, we formulated a time‐dependent, but diversity‐independent model that matches the expected diversity through time of a DD model. We simulated phylogenies under each model and studied how well likelihood methods could recover the true diversification mode. Standard model selection criteria always recovered diversity dependence, even when it was not present. We correct for this bias by using a bootstrap method and find that neither model is decisively supported. This implies that the branching pattern of reconstructed trees contains insufficient information to detect the presence or absence of diversity dependence. We advocate that tests encompassing additional data, for example, traits or range distributions, are needed to evaluate how diversity drives macroevolutionary trends.

Using paleo-archives to safeguard biodiversity under climate change

Strategies for 21st-century environmental management and conservation under global change require a strong understanding of the biological mechanisms that mediate responses to climate- and human-driven change to successfully mitigate range contractions, extinctions, and the degradation of ecosystem services. Biodiversity responses to past rapid warming events can be followed in situ and over extended periods, using cross-disciplinary approaches that provide cost-effective and scalable information for species' conservation and the maintenance of resilient ecosystems in many bioregions. Beyond the intrinsic knowledge gain such integrative research will increasingly provide the context, tools, and relevant case studies to assist in mitigating climate-driven biodiversity losses in the 21st century and beyond.

The past and future human impact on mammalian diversity

To understand the current biodiversity crisis, it is crucial to determine how humans have affected biodiversity in the past. However, the extent of human involvement in species extinctions from the Late Pleistocene onward remains contentious. Here, we apply Bayesian models to the fossil record to estimate how mammalian extinction rates have changed over the past 126,000 years, inferring specific times of rate increases. We specifically test the hypothesis of human-caused extinctions by using posterior predictive methods. We find that human population size is able to predict past extinctions with 96% accuracy. Predictors based on past climate, in contrast, perform no better than expected by chance, suggesting that climate had a negligible impact on global mammal extinctions. Based on current trends, we predict for the near future a rate escalation of unprecedented magnitude. Our results provide a comprehensive assessment of the human impact on past and predicted future extinctions of mammals.

Fossil and phylogenetic analyses reveal recurrent periods of diversification and extinction in dictyopteran insects

Variations of speciation and extinction rates determine the fate of clades through time. Periods of high diversification and extinction (possibly mass-extinction events) can punctuate the evolutionary history of various clades, but they remain loosely defined for many biological groups, especially nonmarine invertebrates like insects. Here, we examine whether the cockroaches, mantises and termites (altogether included in Dictyoptera) have experienced episodic pulses of speciation or extinction and how these pulses may be associated with environmental fluctuations or mass extinctions. We relied on molecular phylogeny and fossil data to shed light on the times and rates at which dictyopterans diversified. The diversification of Dictyoptera has alternated between (i) periods of high diversification in the late Carboniferous, Early-Middle Triassic, Early Cretaceous and middle Palaeogene, and (ii) periods of high extinction rates particularly at the Permian-Triassic boundary, but not necessarily correlated with the major global biodiversity crises as in the mid-Cretaceous. This study advocates the importance of analyzing, when possible, both molecular phylogeny and fossil data to unveil diversification and extinction periods for a given group. The causes and consequences of extinction must be studied beyond mass-extinction events alone to gain a broader understanding of how clades wax and wane.

A 450 million years long latitudinal gradient in age-dependent extinction

Leigh Van Valen famously stated that under constant conditions extinction probability is independent of species age. To test this 'law of constant extinction', we developed a new method using deep learning to infer age‐dependent extinction and analysed 450 myr of marine life across 21 invertebrate clades. We show that extinction rate significantly decreases with age in > 90% of the cases, indicating that most species died out soon after their appearance while those which survived experienced ever decreasing extinction risk. This age‐dependent extinction pattern is stronger towards the Equator and holds true when the potential effects of mass extinctions and taxonomic inflation are accounted for. These results suggest that the effect of biological interactions on age‐dependent extinction rate is more intense towards the tropics. We propose that the latitudinal diversity gradient and selection at the species level account for this exceptional, yet little recognised, macroevolutionary and macroecological pattern.

Climate cooling and clade competition likely drove the decline of lamniform sharks

Understanding heterogeneity in species richness between closely related clades is a key research question in ecology and evolutionary biology. Multiple hypotheses have been proposed to interpret such diversity contrasts across the tree of life, with most studies focusing on speciation rates to explain clades' evolutionary radiations, while often neglecting extinction rates. Here we study a notorious biological model as exemplified by the sister relationships between mackerel sharks (Lamniformes, 15 extant species) and ground sharks (Carcharhiniformes, ∼290 extant species). Using a comprehensive fossil dataset, we found that the diversity dynamics of lamniforms waxed and waned following repeated cycles of radiation phases and declining phases. Radiation phases peaked up to 3 times the current diversity in the early Late Cretaceous. In the last 20 million years, the group declined to its present-day diversity. Along with a higher extinction risk for young species, we further show that this declining pattern is likely attributed to a combination of abiotic and biotic factors, with a cooling-driven extinction (negative correlation between temperature and extinction) and clade competition with some ground sharks. Competition from multiple clades successively drove the demise and replacement of mackerel sharks due to a failure to originate facing the rise of ground sharks, particularly since the Eocene. These effects came from ecologically similar carcharhiniform species inhibiting diversification of medium- and large-sized lamniforms. These results imply that the interplay between abiotic and biotic drivers had a substantial role in extinction and speciation, respectively, which determines the sequential rise and decline of marine apex predators.

Closing the gap between palaeontological and neontological speciation and extinction rate estimates

Measuring the pace at which speciation and extinction occur is fundamental to understanding the origin and evolution of biodiversity. Both the fossil record and molecular phylogenies of living species can provide independent estimates of speciation and extinction rates, but often produce strikingly divergent results. Despite its implications, the theoretical reasons for this discrepancy remain unknown. Here, we reveal a conceptual and methodological basis able to reconcile palaeontological and molecular evidence: discrepancies are driven by different implicit assumptions about the processes of speciation and species evolution in palaeontological and neontological analyses. We present the "birth-death chronospecies" model that clarifies the definition of speciation and extinction processes allowing for a coherent joint analysis of fossil and phylogenetic data. Using simulations and empirical analyses we demonstrate not only that this model explains much of the apparent incongruence between fossils and phylogenies, but that differences in rate estimates are actually informative about the prevalence of different speciation modes.

History is written by the victors: The effect of the push of the past on the fossil record

Survivorship biases can generate remarkable apparent rate heterogeneities through time in otherwise homogeneous birth‐death models of phylogenies. They are a potential explanation for many striking patterns seen in the fossil record and molecular phylogenies. One such bias is the “push of the past”: clades that survived a substantial length of time are likely to have experienced a high rate of early diversification. This creates the illusion of a secular rate slow‐down through time that is, rather, a reversion to the mean. An extra effect increasing early rates of lineage generation is also seen in large clades. These biases are important but relatively neglected influences on many aspects of diversification patterns in the fossil record and elsewhere, such as diversification spikes after mass extinctions and at the origins of clades; they also influence rates of fossilization, changes in rates of phenotypic evolution and even molecular clocks. These inevitable features of surviving and/or large clades should thus not be generalized to the diversification process as a whole without additional study of small and extinct clades, and raise questions about many of the traditional explanations of the patterns seen in the fossil record.

TreeSimGM: Simulating phylogenetic trees under general Bellman-Harris models with lineage-specific shifts of speciation and extinction in R

Understanding macroevolutionary processes using phylogenetic trees is a challenging and complex process that draws on mathematics, computer science and biology. Given the development of complex mathematical models and the growing computational processing power, simulation tools are becoming increasingly popular.In order to simulate phylogenetic trees, most evolutionary biologists are forced to build their own algorithms or use existing tools built on different platforms and/or as standalone programmes. The absence of a simulation tool accommodating for user-chosen model specifications limits, amongst others, model testing and pipelining with approximate Bayesian computation methods or other subsequent statistical analysis.We introduce "TreeSimGM," an r-package simulation tool for phylogenetic trees under a general Bellman and Harris model. This package allows the user to specify any desired probability distribution for the waiting times until speciation and extinction (e.g. age-dependent speciation/extinction). Upon speciation, the user can specify whether one descendant species corresponds to the ancestor species inheriting its age or whether both descendant species are new species of age 0. Moreover, it is possible to scale the waiting time to speciation/extinction for newly formed species. Thus, "TreeSimGM" not only allows the user to simulate stochastic phylogenetic trees assuming several popular existing models, such as the Yule model, the constant-rate birth-death model, and proportional to distinguishable arrangement models, but it also allows the user to formulate new models for exploration. A short explanation of the supported models and a few examples of how to use our package are presented here.As an r-package, "TreeSimGM" allows flexible and powerful stochastic phylogenetic tree simulations. Moreover, it facilitates the pipelining of outputs or inputs with other functions in r. "TreeSimGM" contributes to the tools available to the r community in the fields of ecology and evolution, is freely available under the GPL-2 licence and can be downloaded at https://cran.r-project.org/web/packages/TreeSimGM.