Five major shifts of diversification through the long evolutionary history of Magnoliidae (angiosperms)
© Massoni et al.; licensee BioMed Central. 2015
Received: 10 October 2014
Accepted: 24 February 2015
Published: 18 March 2015
With 10,000 species, Magnoliidae are the largest clade of flowering plants outside monocots and eudicots. Despite an ancient and rich fossil history, the tempo and mode of diversification of Magnoliidae remain poorly known. Using a molecular data set of 12 markers and 220 species (representing >75% of genera in Magnoliidae) and six robust, internal fossil age constraints, we estimate divergence times and significant shifts of diversification across the clade. In addition, we test the sensitivity of magnoliid divergence times to the choice of relaxed clock model and various maximum age constraints for the angiosperms.
Compared with previous work, our study tends to push back in time the age of the crown node of Magnoliidae (178.78-126.82 million years, Myr), and of the four orders, Canellales (143.18-125.90 Myr), Piperales (158.11-88.15 Myr), Laurales (165.62-112.05 Myr), and Magnoliales (164.09-114.75 Myr). Although families vary in crown ages, Magnoliidae appear to have diversified into most extant families by the end of the Cretaceous. The strongly imbalanced distribution of extant diversity within Magnoliidae appears to be best explained by models of diversification with 6 to 13 shifts in net diversification rates. Significant increases are inferred within Piperaceae and Annonaceae, while the low species richness of Calycanthaceae, Degeneriaceae, and Himantandraceae appears to be the result of decreases in both speciation and extinction rates.
This study provides a new time scale for the evolutionary history of an important, but underexplored, part of the tree of angiosperms. The ages of the main clades of Magnoliidae (above the family level) are older than previously thought, and in several lineages, there were significant increases and decreases in net diversification rates. This study is a new robust framework for future investigations of trait evolution and of factors influencing diversification in this group as well as angiosperms as a whole.
Understanding the diversification history of species-rich clades is a major goal in evolutionary biology  as they provide important insights into the evolution of life on earth. With over 10,000 plant species, Magnoliidae, sensu Cantino et al. , are the largest clade of flowering plants outside monocots and eudicots and comprise four orders: Canellales (123 spp.), Laurales (3874 spp., incl. avocado), Magnoliales (2978 spp., incl. magnolias), and Piperales (3190 spp., incl. black pepper). Investigating the diversification of Magnoliidae is of primary importance to provide guidance for future hypotheses on the drivers of diversification of flowering plants as a whole. In addition, many organisms (incl. various butterfly and beetle groups) are highly dependent on this group for feeding or reproduction [3,4], and magnoliid species are an important part of tropical ecosystems [5,6]. As a consequence, investigating their tempo and mode of diversification will provide key knowledge to understand not only the evolutionary history of Magnoliidae as a whole, but also that of these related groups  and of the environments they live in .
It has often been suggested that the biology of angiosperms has a great influence on their own diversification [9,10]. The biology of Magnoliidae differs from that of eudicots and monocots in several respects. For instance, most members of Magnoliidae are mainly pollinated by beetles, flies, or thrips, while, in contrast to eudicots and monocots, bee or wind pollination are rare [3,11]. In addition, floral and vegetative morphology of Magnoliidae is highly variable and several of these variable traits have been shown to be associated with variation of diversification rates in other flowering plant lineages .
Previous molecular dating studies at large (e.g., angiosperms) or narrow (e.g., a particular family) taxonomic scale have suggested that the crown node age of Magnoliidae is older than the crown nodes of both eudicots and monocots, ranging from 240.16 to 85.07 Myr [12,13], and that several of the 19 families of Magnoliidae appear to have originated before the end of the Cretaceous [13-22]. In addition, the ages of numerous nodes within Magnoliidae have yet to be estimated . In parallel, few studies have focused on the diversification of the group, being either at the scale of all angiosperms  or within families [25,26].
In the present study, we provide for the first time a complete time scale for Magnoliidae at the familial level. We then use this new time scale to test for significant diversification rate shifts during the evolutionary history of the clade. To do so, we take advantage of a recent improvement in the phylogeny of the group [16,27-35] and a new revision of the rich fossil record of Magnoliidae providing reliable minimum-age calibration points . In order to conduct these analyses for nodes above and below the family level, we used a dataset including more than 75% of the existing genera and 12 molecular markers from the three plant genomes (6 plastid, 4 mitochondrial, 2 nuclear).
All runs in each of the five BEAST analyses with different calibration scheme (angio-130, −140, −150, −170 and −200) converged effectively. The concatenated post-burn-in samples of the four runs in each analyses showed final effective sampling sizes (ESS) of the likelihood far above 200 [see Additional file 1]. Except for the position of Siparunaceae (Laurales) and Magnoliaceae (Magnoliales), which remained unresolved in all analyses, the great majority of the relationships among the families received high support values [see Additional file 2] and is similar in all analyses [see Additional file 3].
Conducting BEAST analyses without molecular data highlighted the non-uniformity of the temporal constraints effectively used in the analyses (results not shown). However, no effective temporal priors violated the paleontological information used to calibrate the present analyses. In addition, the temporal posterior distributions of calibrated nodes obtained with the molecular data differ from the effective temporal prior distributions obtained without these data. These mismatches for calibrated nodes suggest that molecular data are well informative.
Age estimates for selected nodes of Magnoliidae
95% CI (mean or PL)
95% IC (mean or PL)
95% CI (mean or PL)
95% CI (mean or PL)
95% CI (mean or PL)
(Canellales + Piperales)
(Laurales + Magnoliales)
Tempo of Magnoliidae evolution
We provide here for the first time a near complete temporal framework for the evolution of Magnoliidae above the generic level. Applying different maximum age constraints to the crown node of angiosperms highlighted the sensitivity of inferred ages estimated for the deepest nodes of Magnoliidae (Figure 1). Although previous studies have recognized a long gap in the fossil record of the angiosperm stem lineage [13,21,37], there is no indisputable argument in favor of one particular maximum age constraint for the crown node of angiosperms. Here we will systematically consider the whole range of age estimates obtained across our five maximum age calibration schemes and both molecular dating approaches used.
We restricted ourselves to conservative fossil age constraints based on phylogenetic analyses of fossil and extant taxa . This led us to consider only 10 out of more than 100 described fossils putatively belonging to Magnoliidae . However, our new age estimates for families and orders of Magnoliidae appear, in general, to be compatible with the putative fossil record attributed to each of these taxa. Future investigations of the phylogenetic placement of other magnoliid fossils will likely provide additional minimum age constraints, which could push back in time several young estimates supported in the present study, and decrease the size of the credibility intervals. This study further illustrates that, just as paleontological dating, molecular dating estimates are often associated with large uncertainties. The latter approach has to be seen as an attempt to reduce the range of the most likely ages for nodes constrained by the age of the fossil record securely placed and dated, and to evaluate the probability of the ages of nodes for which there is no direct fossil record. Because all the ambiguity of our current knowledge has to be taken into account, the molecular dating approach cannot provide exact secure ages, except for exceptionally fossil-rich clades.
Based on this most reliable knowledge in the field of paleobotany of Magnoliidae  and chrono-stratigraphy , combined with a comprehensive sample of taxa and molecular markers, our data support an origin of extant lineages of Magnoliidae (crown node of the group) during the Jurassic or Lower Cretaceous between 178.78 and 126.82 Myr. This result tend to be older than previously suggested [12,15,21,22,37,40-42]. For example, the angiosperm-wide study of Bell et al. , in which 33 magnoliid genera were sampled, estimated the crown node of Magnoliidae to be between 138 and 108 Myr. In another, more recent angiosperm-wide analysis sampling 34 species of Magnoliidae and using a PL approach, Zanne et al.  estimated the age of the crown node of Magnoliidae to be 147 Myr without uncertainty associated. Our new estimates for the crown-group age of orders also tend to be older than previously inferred [13,15,16,22,37] (Figure 1). However, some studies have supported much older ages for these nodes [13,14,21,43]. At the familial level, in addition to providing the first estimates for the ages of crown-group Canellaceae (52.34-11.33 Myr), Saururaceae (80.79-46.69 Myr), and Winteraceae (56.73-12.11 Myr; Table 1), the crown nodes of families are also generally dated to an older range of dates than previously estimated [16-20,26,44], with the exception of Atherospermataceae (Laurales) for which we found a younger age than previously suggested . The large amount of combinations of parameters influencing age estimates limits straightforward explanations of these differences with previous studies (e.g., topology used ; fossil species and priors used to define and model age constraints ; molecular dating method ). Despite variation among the ages obtained for the crown ages of families, our results confirm that Magnoliidae diversified into morphologically distinct clades (now identified as families) by the end of the Cretaceous (66 Myr, Figure 1), but probably earlier than thought before [16-20,26,44]. This provides arguments in favor of an earlier diversification of angiosperms altogether, with the presence of more flowering plant lineages in the Cretaceous than previously thought.
Finally, taking into account the whole range of age estimates obtained across our five different calibration schemes, leads to long ranges of age estimates that are compatible with biogeographic scenarios previously suggested in Magnoliidae [17,26,44,45,48-52]. However, viewed in isolation, the different assumptions about the maximum age of angiosperms lead to alternative time scales for the evolutionary history of the families of Magnoliidae. Therefore, biogeographic scenarios will need to be reevaluated in the context of explicit assumptions on the age of the angiosperms.
Mode of diversification of Magnoliidae
The new timetree obtained for Magnoliidae in this study allowed us to detect an average of nine significant diversification rate shifts across lineages, implying that diversification has not been a homogeneous process throughout the history of the clade (Figure 2). Explaining the present day diversity of species-rich clades has generally required several diversification rate shifts, but previous studies using the same statistical approach in other large angiosperm clades found, in general, fewer shifts than in the present analysis. Within eudicots, Beaulieu and Donoghue  detected six shifts in Campanulidae (35,000 spp.), Koenen et al.  detected nine shifts in Fabaceae (19,500 spp.), Xi et al.  found five significant shifts in Malpighiales (16,000 spp.), and Arakaki et al.  found seven shifts in Cactaceae (1850 spp.). Within monocots, Escudero and Hipp  supported three shifts in the family Cyperaceae (5480 spp.) while Baker and Couvreur  identified up to 13 shifts in palms (2500 spp.). Despite the dependence of the number of shifts detected on the taxonomic level of tips in the chronograms used (MEDUSA cannot detect shifts within supraspecific terminal compartments), it seems that Magnoliidae experienced a high variation in diversification rates. Because of their longer evolutionary history in comparison with these other clades and because their fossil record supports a global distribution during the Cretaceous and the Paleogene , Magnoliidae as a whole may have been affected by more events potentially influencing the speciation and extinction rates (e.g., climatic shifts ; variation in available area ; variation in geographic distribution ; dynamic of diversification of associated pollinators ).
Establishing the causes affecting diversification-rate variation can be difficult. There is a large diversity of intrinsic and extrinsic traits known to influence speciation and extinction , all of them being potentially correlated to each other . The influence of a trait on the birth and death of species is dependent on other taxa, other traits of the same organism, and its physical environment , implying that a trait could have different effects on different clades. Finally, the actual shift of diversification may happen several nodes after the evolution of an influencing trait. One explanation could be that an isolated trait will influence diversification in combination with other characters, the effect appearing when the entire set of traits needed is present . For these reasons, the investigation of the correlation between traits and diversification rates will need to take into account as many potential factors as possible .
The older estimates of the ages in the “angio-200” hypothesis led to slightly lower absolute rates of diversification than in the “angio-140” hypothesis, as expected. Because the general pattern (number and position of shifts) is the same in both analyses, in the present paragraph we refer only to rates of the angio-140 analysis. The background net diversification rate obtained in the tree of Magnoliidae (0.0401 ± 0.0138 sp.myr−1) is comparable to those found in angiosperms [14,24]. In the present study, the lowest net diversification rates, resulting from shifts 3 and 4 (Figure 2), were found in the crown group of Calycanthaceae (0.0175 ± 0.0066 sp.myr−1), and in the branches of Degeneriaceae and Himantandraceae (0.0057 ± 0.0088 sp.myr−1). These three families share distinctive floral features not or seldom found in other Magnoliidae groups, such as the presence of inner staminodes that cover the stigma at the end of the female phase , and spiral phyllotaxy of all floral organs [65,66]. However, it is difficult to link these two characters with mechanisms sustaining a low diversification rate (influence on speciation and/or extinction). Within Calycanthaceae, the low rate of net diversification is associated with a low speciation rate for the entire clade (0.0182 sp.myr−1) and near-zero extinction (0.0007 sp.myr−1). The absence of extinction would seem contradictory with the presence of several fossils in the Cretaceous  and could be an artefact of the method which does not take into account the extinct diversity (see below). However, because the phylogenetic position of these fossils within the family has not been tested (except for Virginianthus calycanthoides ), it is difficult to draw conclusions about extinction within the family. The pollination systems of Calycanthaceae involve trapping of a large variety of beetles and thrips [68,69]. In comparison with specific pollination systems involving a unique relationship for the two partners, this generalist interaction is not favorable for genetic isolation promoting diversification . In addition, Calycanthaceae have a temperate distribution, except for the monospecific genus Idiospermum growing in a restricted area in north-east Australia. The low ecological limits on species richness in temperate areas  could provide an explanation for the low diversification of the lineage leading to the genera Calycanthus and Chimonanthus. On the other hand, the apparent lack of diversification of the lineage leading to Idiospermum could be explained by different causes. Individuals of this species are long-living trees  involving low fixation rates in diverging populations. The genus is restricted to the very wet humid tropical lowland rainforests of Australia . Idiospermum presents one of the heaviest seeds among Australian plants and is therefore probably not dispersed by animals, limiting its capacity of dispersion [71,73]. If this low capacity of dispersion is ancestral in the lineage leading to the extant species of the genus, it restrained the number of available areas promoting diversification . This latter explanation might be generalized for all Calycanthaceae as the seeds of the family contain secondary metabolites that are toxic at least for mammals .
The highest diversification rate inferred in this study was found on the branch leading to the terminal compartment Miliuseae in subfamily Malmeoideae of Annonaceae (0.1527 ± 0.0548 sp.myr−1). A shift of +0.1231 ± 0.0306 sp.myr−1 at the node sustaining the branch of this clade was present in 71% of the trees tested and the associated rate of speciation on the branch is on average equal to 0.4558 sp.myr−1. Couvreur et al.  found that major clades in Annonaceae have undergone different rates of diversification, with subfamily Malmeoideae (previously the Short-Branch Clade) having significantly higher rates than subfamily Annonoideae (previously the Long-Branch Clade). However, analyses of family-level lineage-through-time plots (LTTs) did not detect significant shits in diversification rates during most of the evolutionary history of Annonaceae . Erkens et al. , using a tree with fewer taxa and two alternative approaches to detect rate shifts, identified up to three diversification rate shifts in Annonaceae, including the one we detected here at the crown of Miliuseae. Reasons for an increase in diversification within Miliuseae are not yet well known but could be related to founder effects after dispersal into South-East Asia or distinctive pollen characters .
Canellales and Piperales have, on average, higher speciation rates than Magnoliales and Laurales. Within Piperales, this may be explained by widespread herbaceous habits, resulting in faster life history (shorter generation time) and higher fixation rates in divergent populations . In addition, within this order, Aristolochiaceae typically have zygomorphic flowers, which has been shown to have a positive impact on diversification . However, Canellales and Piperales also share the highest extinction rates, making Laurales and Magnoliales the most productive orders in terms of net diversification. Differences in extinction rates are difficult to explain. For instance, there are no major differences in pollination systems among the four orders , and their current geographic distributions are similar . However, Doyle and Endress  found several synapomorphies for the clade of Laurales and Magnoliales, specifically the presence of more than two whorls of stamens and more than one whorl or series of carpels. These characters might have played a role in relation to pollination, but the underlying mechanisms remain unknown.
Finally, even though the MEDUSA approach represents a significant conceptual improvement over previous models, it still requires several important assumptions that may have influenced our results. First, MEDUSA does not implement diversification models in which the extinction rate can be higher than speciation rates. Second, the rates are assumed to be constant in each part of the tree, and only abrupt variations are evaluated, probably not in accordance with the biological reality in which rates could gradually fluctuate. Morlon et al.,  recently developed an approach specifically to relax these two assumptions, but it is not implemented yet to detect shifts automatically nor to analyze higher-level phylogenies with terminally unresolved clades. In addition, Rabosky  recently developed a new Bayesian framework to detect multiple diversification processes in a tree while allowing the speciation rate to vary continuously through time, but it is not available yet for a set of multiple trees and therefore cannot take phylogenetic and dating uncertainty into account. Last, because the fossil record is not taken into account in the diversification rate analyses, very low rates of speciation and extremely low extinction rates are always allocated to the old lineages with few species (long branches in chronograms; ), which have probably been more diverse in the past (e.g., Idiospermum). New research efforts are currently being made in this area [80,81]. Therefore, our understanding of Magnoliidae diversification may further improve in the future when more realistic methods and more complete phylogenies are available.
The present study suggests that Magnoliidae began to diversify somewhere between the Toarcian (Early Jurassic) and the Barremian (Early Cretaceous), or 178.78-126.82 Myr ago. Several key nodes within Magnoliidae are dated here for the first time. In general, our age estimates suggest a possibly older diversification of the group than previously inferred. The rich fossil record of Magnoliidae may eventually provide additional calibration points to refine the time scale proposed here. However, considerable work remains to be done to securely relate this extinct diversity to the extant one, a task for which an integrated morphological dataset will be essential. The tempo and mode of diversification within this clade appears to be characterized by several increases and decreases of diversification rates, suggesting that Magnoliidae have not undergone constant diversification but have instead been shaped by alternative, yet to be determined evolutionary processes. Our new dated trees provide a solid basis for future biogeographical studies and robust statistical tests of correlation among intrinsic traits, extrinsic factors, and diversification rates within Magnoliidae.
In order to conduct the molecular dating analyses we used the same 12-marker molecular dataset as . This matrix includes 12 coding and non-coding markers from the three genomes: atpB, matK, trnL intron, trnL-trnF spacer, ndhF, rbcL from the chloroplast; atp1, matR, mtSSU, mtLSU from the mitochondrion; and 18S rDNA and 26S rDNA from the nucleus. In this dataset, we used an exemplar approach, in which each genus was represented by one species. The problematic parasitic family Hydnoraceae was excluded because its exact position within Piperales remains difficult to assess [16,33]. All remaining 19 families of Magnoliidae are represented, and more than 75% of the genera were sampled (198 genera out of 262). In addition, we included 23 outgroup taxa sampled from early-diverging angiosperms (Amborellales, Nymphaeales, Austrobaileyales, Chloranthales) as well as eudicots and monocots [see Additional file 3].
Molecular dating analyses
Fossil species used to define the calibration scheme presented in Massoni et al. 
crown-group core Laurales
crown-group core Laurales
stem node of extant Saururus
Divergence time estimation
Uncorrelated Lognormal Clock (UCLN) analyses
To evaluate divergence times within Magnoliidae while taking into account phylogenetic uncertainty, we used BEAST v1.7.5  without fixing the tree. In all BEAST analyses, we partitioned our molecular dataset according to the 12 markers (as in Massoni et al. ). We used MrModeltest 2.3  to evaluate the best fitted model for each partition. For all partitions, GTR + Gamma was selected as the most appropriate model according to the Akaike Information Criterion, and the evaluation of invariant sites (I) was recommended for all, except for the plastid spacer trnL-trnF. A Birth-Death incomplete sampling prior was specified for the trees , and rate heterogeneity was modeled using the UCLN relaxed clock of Drummond et al.  with a uniform prior for the mean of the branch rate set between 0 and 1E100. All age constraints were applied using uniform priors (hard minimum and maximum ages) because of the difficulty of parameterizing non-uniform priors, which may involve subjectivity . These uniform distributions were bounded either by the age of the fossil (Table 2) and 1E100 Myr for minimum age constraints, or by the appropriate maximum age (see above) and 0 Myr for maximum age constraints. Because in BEAST the construction of the calibrated tree prior is a multiplicative construction involving the multiplication of the calibration densities (specified by the user) with the tree prior, the effective calibration priors used in analyses can differ from the specified priors . This can cause the effective prior to violate the desired paleontological constraints . We conducted all BEAST analyses without the data in order to estimate the mismatch between the specified and effective calibration priors, but also the effect of data on posterior distributions. Because of initial likelihood problems at the start of the analyses we used a starting tree obtained from the Bayesian analysis of Massoni et al.  and rendered ultrametric using Penalized Likelihood (PL) implemented in r8s v1.8 . Because of rounding to a limited number of decimals, branch times in this starting tree summed up to slightly older node ages than its original calibrations. Thus, in order to use this starting tree, we had to apply slightly older (+0.1 Myr) maximum age constraints (126.8, 130.1, 140.1, 150.1, 170.1, and 200.1 Myr). For all analyses, we ran four independent chains of a Markov Chain Monte Carlo procedure for 100 million generations each, sampling parameters and trees every 1000 generations [see Additional file 1]. All analyses were performed on the CIPRES cluster . We evaluated the size of the burnin phase for each run using Tracer v1.5 . The post-burnin posteriors of the four runs of each analysis were then combined using LogCombiner v1.7.5. Because of computational limitations in relation to the large size of generated files, we re-sampled trees and parameters every 10,000 or 15,000 generations [see Additional file 1]. We then used TreeAnnotator v1.7 to select the maximum clade credibility (MCC) tree of each analysis.
Penalized Likelihood (PL) analyses
We also conducted PL analyses for each set of age constraints using r8s v1.8 . The PL relaxed clock assumes some degree of autocorrelation of molecular substitution rates between a parent and its immediate descendants . Although autocorrelated relaxed clock models are also available in Bayesian frameworks, our rationale for using r8s was to test the sensitivity of our age estimates to a fundamentally different approach to molecular dating. In all analyses, we used the best scoring maximum likelihood (ML) phylogram of Massoni et al.  obtained with the 12-marker RAxML analysis excluding Hydnora. In order to determine the optimal level of autocorrelation across the tree (smoothing parameter), we conducted cross validation procedures testing 17 different values of smoothing in a range between 0.1 and 10,000,000 for each calibration scheme, using the penalty additive function. For angio-130, angio-140, angio-170 and angio-200 the optimal values were 10, and for angio-150 it was 3.2. In order to provide confidence intervals on age estimates with PL, we conducted non-parametric maximum likelihood bootstrapping of 1000 replicates on the original dataset of Massoni et al.  using RAxML v7.3.2  while fixing the tree topology to the best scoring ML tree obtained in this previous study. As RAxML requires the application of the same model to all partitions, we used the GTR + GAMMA + I model for each partition. As mentioned above, this model was selected for all partitions of the present dataset except for the plastid spacer trnL-trnF for which the invariant site was not recommended. Within the resulting 1000-tree collection, only branch lengths varied. We reconstructed the chronograms for both the best ML tree and the 1000 bootstrapped phylograms using the optimal smoothing parameter obtained for the best ML tree and the TN algorithm of r8s. All results were summarized using the software TreeAnnotator v1.7 with the 1000 PL trees as the input file and the ML-PL tree as the target tree.
We used the MEDUSA approach of Alfaro et al.  to test for significant diversification rate shifts in an incompletely sampled phylogeny where each terminal is assigned the total number of extant species it represents. MEDUSA detects, using a stepwise AIC approach, significant shifts of diversification rates across a given phylogeny by evaluating the fit of different pairwise birth-death models. In order to select the best model we used the AICc criterion, a birth-death model of diversification, and we allowed the placement of shifts either on stem or on crown nodes. To assess the sensitivity of MEDUSA to age and phylogenetic uncertainty, we conducted the analyses on two collections of 1000 posterior BEAST trees randomly sampled from the post-burnin phase, corresponding to calibration scenarios angio-140 and angio-200. Prior to the analyses, each randomly sampled chronogram was transformed in two ways, using functions in the ape package of R . First, all outgroups of Magnoliidae were pruned. Second, the chronograms were simplified by pruning selected terminal taxa so that each remaining terminal taxon would represent a monophyletic compartment with known extant diversity (see below). All analyses were performed in R using the multiMEDUSA procedure from package MEDUSA v0.93 4.33 available on the website authored by Joseph W. Brown (https://github.com/josephwb/turboMEDUSA).
Extant taxonomic richness data
In order to assign the entire species richness of Magnoliidae to the tips of the phylogeny, it was necessary in several instances to merge tips into a single terminal to represent larger monophyletic groups hereafter referred to as compartments. We proceeded in two steps. First, we maintained genera as tips or created larger compartments, depending on the monophyletic status of genera as tested in previous phylogenetic studies: if the literature supported the monophyly of a genus, this taxon was used as a compartment; if the monophyly of a genus had never been tested before, a compartment including this genus and its sister group was defined (or a larger compartment if the neighborhood relationships were not well supported); if the monophyly of a genus had been challenged by previous studies, we defined a larger compartment to include this genus and all other genera potentially involved in the paraphyly or polyphyly. After this first step, we modified our compartmentalization scheme to take into account the species numbers of genera not sampled in our trees (less than 25% of the total number of accepted genera). If there was enough information in the literature to support an accurate placement (phylogenetic or apomorphy-based) of a missing genus, we added its number of species to that of an earlier defined compartment (generic or supra-generic). In some cases, the missing genus could not be precisely placed within a larger clade of two or more previously defined compartments. As a guideline, we accepted to merge these compartments into a larger one if the number of species would represent more than three percent of the total number of species in the resulting compartment. Otherwise, the missing diversity was ignored in order to maintain enough compartments for conducting a meaningful analysis. In all cases, we defined these supra-generic compartments in such a way that they were present in all 1000 BEAST trees used to conduct the MEDUSA analyses (posterior probabilities [PP] equal to 100%). The resulting compartmentalization scheme consisted in 85 terminal taxa [see Additional file 5]. In total, 31 species (16 genera) were ignored due to unknown phylogenetic placement, representing ca. 0.3% of the 10,209 species of Magnoliidae. For 15 genera (out of 262) the number of species was not clearly indicated in previous publications in which case the Plant List (http://www.theplantlist.org/) was used to estimate the number of extant species currently accepted in the genus. The definition of supra-generic clades and the incorporation of the missing diversity are justified in detail in Additional file 5.
Availability of supporting data
Nexus files including the MCC trees of the BEAST analyses and the PL trees of the r8s analyses (annotated with the 95% credibility intervals of age estimates) are available in Dryad Digital repository, http://0-dx.doi.org.brum.beds.ac.uk/10.5061/dryad.ct231 .
We thank Luna Sánchez Reyes, Elisabeth Reyes, Hélène Morlon, Susana Magallón, James Doyle, and Michael Pirie for discussions on various aspects of this study. Two anonymous reviewers are also gratefully acknowledged for their helpful comments on this paper. This work was supported by the Agence Nationale de la Recherche (grant number ANR-12-JVS7-0015-01 to HS).
- Alfaro ME, Santini F, Brock C, Alamillo H, Dornburg A, Rabosky DL, et al. Nine exceptional radiations plus high turnover explain species diversity in jawed vertebrates. Proc Natl Acad Sci U S A. 2009;106:13410–4.View ArticlePubMed CentralPubMedGoogle Scholar
- Cantino PD, Doyle JA, Graham SW, Judd WS, Olmstead RG, Soltis DE, et al. Towards a phylogenetic nomenclature of Tracheophyta. Taxon. 2007;56:E1–44.View ArticleGoogle Scholar
- Gottsberger G. How diverse are Annonaceae with regard to pollination? Bot J Linn Soc. 2012;169:245–61.View ArticleGoogle Scholar
- Condamine FL, Sperling FAH, Wahlberg N, Rasplus J-Y, Kergoat GJ. What causes latitudinal gradients in species diversity? Evolutionary processes and ecological constraints on swallowtail biodiversity. Ecol Lett. 2012;15:267–77.View ArticlePubMedGoogle Scholar
- Ayyappan N, Parthasarathy N. Biodiversity inventory of trees in a large-scale permanent plot of tropical evergreen forest at Varagalaiar, Anamalais, Western Ghats, India. Biodivers Conserv. 1999;8:1533–54.View ArticleGoogle Scholar
- Pitman NCA, Terborgh JW, Silman MR, Núñez VP, Neill DA, Cerón CE, et al. Dominance and distribution of tree species in upper Amazonian Terra Firme forests. Ecology. 2001;82:2101–17.View ArticleGoogle Scholar
- McKenna DD, Sequeira AS, Marvaldi AE, Farrell BD. Temporal lags and overlap in the diversification of weevils and flowering plants. Proc Natl Acad Sci U S A. 2009;106:7083–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Arakaki M, Christin P-A, Nyffeler R, Lendel A, Eggli U, Ogburn RM, et al. Contemporaneous and recent radiations of the world’s major succulent plant lineages. Proc Natl Acad Sci U S A. 2011;108:8379–84.View ArticlePubMed CentralPubMedGoogle Scholar
- Van der Niet T, Johnson SD. Phylogenetic evidence for pollinator-driven diversification of angiosperms. Trends Ecol Evol. 2012;27:353–61.View ArticlePubMedGoogle Scholar
- Vamosi JC, Vamosi SM. Factors influencing diversification in angiosperms: at the crossroads of intrinsic and extrinsic traits. Am J Bot. 2011;98:460–71.View ArticlePubMedGoogle Scholar
- Thien LB, Azuma H, Kawano S. New perspectives on the pollination biology of basal angiosperms. Curr Perspect Basal Angiosperms. 2000;161 Suppl 6:225–35.Google Scholar
- Bell CD, Soltis DE, Soltis PS. The age of the angiosperms: a molecular timescale without a clock. Evolution. 2005;59:1245–58.View ArticlePubMedGoogle Scholar
- Magallón SA. Using fossils to break long branches in molecular dating: a comparison of relaxed clocks applied to the origin of angiosperms. Syst Biol. 2010;59:384–99.View ArticlePubMedGoogle Scholar
- Magallón SA, Castillo A. Angiosperm diversification through time. Am J Bot. 2009;96:349–65.View ArticlePubMedGoogle Scholar
- Wikström N, Savolainen V, Chase MW. Evolution of the angiosperms: calibrating the family tree. Proc R Soc B. 2001;268:2211–20.View ArticlePubMed CentralPubMedGoogle Scholar
- Naumann J, Salomo K, Der JP, Wafula EK, Bolin JF, Maass E, et al. Single-copy nuclear genes place haustorial Hydnoraceae within Piperales and reveal a Cretaceous origin of multiple parasitic angiosperm lineages. PLoS One. 2013. doi:10.1371/journal.pone.0079204.Google Scholar
- Zhou S, Renner SS, Wen J. Molecular phylogeny and intra- and intercontinental biogeography of Calycanthaceae. Mol Phylogenet Evol. 2006;39:1–15.View ArticlePubMedGoogle Scholar
- Renner SS. Variation in diversity among laurales, early cretaceous to present. Biol Skr Det K Danske Vidensk Selsk. 2004;55:441–58.Google Scholar
- Doyle JA, Sauquet H, Scharaschkin T, Le Thomas A. Phylogeny, molecular and fossil dating, and biogeographic history of Annonaceae and Myristicaceae (Magnoliales). Int J Plant Sci. 2004;165:55–67.View ArticleGoogle Scholar
- Pirie MD, Doyle JA. Dating clades with fossils and molecules: the case of Annonaceae. Bot J Linn Soc. 2012;169:84–116.View ArticleGoogle Scholar
- Smith SA, Beaulieu JM, Donoghue MJ. An uncorrelated relaxed-clock analysis suggests an earlier origin for flowering plants. Proc Natl Acad Sci U S A. 2010;107:5897–902.View ArticlePubMed CentralPubMedGoogle Scholar
- Magallón SA, Hilu KW, Quandt D. Land plant evolutionary timeline: gene effects are secondary to fossil constraints in relaxed clock estimation of age and substitution rates. Am J Bot. 2013;100:556–73.View ArticlePubMedGoogle Scholar
- Forest F, Chase MW. Magnoliids. In: Hedges SB, Kumar S, editors. Time tree life. Oxford, UK: Oxford University Press; 2009. p. 166–8.Google Scholar
- Magallón SA, Sanderson MJ. Absolute diversification rates in angiosperm clades. Evolution. 2001;55:1762–80.View ArticlePubMedGoogle Scholar
- Erkens RHJ, Chatrou LW, Couvreur TLP. Radiations and key innovations in an early branching angiosperm lineage (Annonaceae; Magnoliales). Bot J Linn Soc. 2012;169:117–34.View ArticleGoogle Scholar
- Couvreur TLP, Pirie MD, Chatrou LW, Saunders RMK, Su YCF, Richardson JE, et al. Early evolutionary history of the flowering plant family Annonaceae: steady diversification and boreotropical geodispersal. J Biogeogr. 2011;38:664–80.View ArticleGoogle Scholar
- Qiu Y-L, Dombrovska O, Lee J, Li L, Whitlock BA, Bernasconi-Quadroni F, et al. Phylogenetic analyses of basal angiosperms based on nine plastid, mitochondrial, and nuclear genes. Int J Plant Sci. 2005;166:815–42.View ArticleGoogle Scholar
- Graham SW, Olmstead RG. Utility of 17 chloroplast genes for inferring the phylogeny of the basal angiosperms. Am J Bot. 2000;87:1712–30.View ArticlePubMedGoogle Scholar
- Soltis DE, Smith SA, Cellinese N, Wurdack KJ, Tank DC, Brockington SF, et al. Angiosperm phylogeny: 17 genes, 640 taxa. Am J Bot. 2011;98:704–30.View ArticlePubMedGoogle Scholar
- Zanis MJ, Soltis DE, Soltis PS, Mathews S, Donoghue MJ. The root of the angiosperms revisited. Proc Natl Acad Sci U S A. 2002;99:6848–53.View ArticlePubMed CentralPubMedGoogle Scholar
- Hilu KW, Borsch T, Müller K, Soltis DE, Soltis PS, Savolainen V, et al. Angiosperm phylogeny based on matK sequence information. Am J Bot. 2003;90:1758–76.View ArticlePubMedGoogle Scholar
- Moore MJ, Soltis PS, Bell CD, Burleigh JG, Soltis DE. Phylogenetic analysis of 83 plastid genes further resolves the early diversification of eudicots. Proc Natl Acad Sci U S A. 2010;107:4623–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Massoni J, Forest F, Sauquet H. Increased sampling of both genes and taxa improves resolution of phylogenetic relationships within Magnoliidae, a large and early-diverging clade of angiosperms. Mol Phylogenet Evol. 2014;70:84–93.View ArticlePubMedGoogle Scholar
- Sauquet H, Doyle JA, Scharaschkin T, Borsch T, Hilu KW, Chatrou LW, et al. Phylogenetic analysis of Magnoliales and Myristicaceae based on multiple data sets: implications for character evolution. Bot J Linn Soc. 2003;142:125–86.View ArticleGoogle Scholar
- Renner SS. Circumscription and phylogeny of the Laurales: evidence from molecular and morphological data. Am J Bot. 1999;86:1301–15.View ArticlePubMedGoogle Scholar
- Massoni J, Doyle JA, Sauquet H. Fossil calibration of Magnoliidae, an ancient lineage of angiosperms. Palaeontol Electron. 2015;18.1.2FC:1-25Google Scholar
- Bell CD, Soltis DE, Soltis PS. The age and diversification of the angiosperms re-revisited. Am J Bot. 2010;97:1296–303.View ArticlePubMedGoogle Scholar
- Friis EM, Crane PR, Pedersen KR. Early flowers and angiosperms evolution. New York: Cambridge University Press; 2011. p. 585.View ArticleGoogle Scholar
- Gradstein FM, Ogg JG, Schmitz M, Ogg G. The geologic time scale 2012. Amsterdam: Elsvier; 2012.Google Scholar
- Moore MJ, Bell CD, Soltis PS, Soltis DE. Using plastid genome-scale data to resolve enigmatic relationships among basal angiosperms. Proc Natl Acad Sci U S A. 2007;104:19363–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Soltis DE, Bell CD, Kim S, Soltis PS. Origin and early evolution of angiosperms. Ann N Y Acad Sci. 2008;1133:3–25.View ArticlePubMedGoogle Scholar
- Zanne AE, Tank DC, Cornwell WK, Eastman JM, Smith SA, FitzJohn RG, et al. Three keys to the radiation of angiosperms into freezing environments. Nature. 2014;506:89–92.View ArticlePubMedGoogle Scholar
- Schneider H, Schuettpelz E, Pryer KM, Cranfill R, Magallón SA, Lupia R. Ferns diversified in the shadow of angiosperms. Nature. 2004;428:553–7.View ArticlePubMedGoogle Scholar
- Su YCF, Saunders RMK. Evolutionary divergence times in the Annonaceae: evidence of a late Miocene origin of Pseuduvaria in Sundaland with subsequent diversification in New Guinea. BMC Evol Biol. 2009;9:153.View ArticlePubMed CentralPubMedGoogle Scholar
- Renner SS, Foreman DB, Murray D. Timing transantarctic disjunctions in the Atherospermataceae (Laurales): evidence from coding and noncoding chloroplast sequences. Syst Biol. 2000;49:579–91.View ArticlePubMedGoogle Scholar
- Sanderson MJ, Doyle JA. Sources of error and confidence intervals in estimating the age of angiosperms from rbcL and 18S rDNA data. Am J Bot. 2001;88:1499–516.View ArticlePubMedGoogle Scholar
- Sauquet H, Ho SYW, Gandolfo MA, Jordan GJ, Wilf P, Cantrill DJ, et al. Testing the impact of calibration on molecular divergence times using a fossil-rich group: the case of Nothofagus (Fagales). Syst Biol. 2012;61:289–313.View ArticlePubMedGoogle Scholar
- Richardson JE, Chatrou LW, Mols JB, Erkens RHJ, Pirie MD. Historical biogeography of two cosmopolitan families of flowering plants: annonaceae and rhamnaceae. Philos Trans R Soc Lond B Biol Sci. 2004;359:1495–508.View ArticlePubMed CentralPubMedGoogle Scholar
- Pirie MD, Chatrou LW, Mols JB, Erkens RHJ, Oosterhof J. “Andean-centred” genera in the short-branch clade of annonaceae: testing biogeographical hypotheses using phylogeny reconstruction and molecular dating. J Biogeogr. 2006;33:31–46.View ArticleGoogle Scholar
- Michalak I, Zhang L-B, Renner SS. Trans-Atlantic, trans-Pacific and trans-Indian Ocean dispersal in the small Gondwanan Laurales family Hernandiaceae. J Biogeogr. 2010;37:1214–26.View ArticleGoogle Scholar
- Chanderbali AS, van der Werff H, Renner SS. Phylogeny and historical biogeography of Lauraceae: evidence from the chloroplast and nuclear genomes. Ann Missouri Bot Gard. 2001;88:104–34.View ArticleGoogle Scholar
- Renner SS, Strijk JS, Strasberg D, Thébaud C. Biogeography of the Monimiaceae (Laurales): a role for East Gondwana and long-distance dispersal, but not West Gondwana. J Biogeogr. 2010;37:1227–38.View ArticleGoogle Scholar
- Beaulieu JM, Donoghue MJ. Fruit evolution and diversification in campanulid angiosperms. Evolution. 2013;67:3132–44.View ArticlePubMedGoogle Scholar
- Koenen EJM, de Vos JM, Atchison GW, Simon MF, Schrire BD, de Souza ER, et al. Exploring the tempo of species diversification in legumes. South African J Bot. 2013;89:19–30.View ArticleGoogle Scholar
- Xi Z, Ruhfel BR, Schaefer H, Amorim AM, Sugumaran M, Wurdack KJ, et al. Phylogenomics and a posteriori data partitioning resolve the Cretaceous angiosperm radiation Malpighiales. Proc Natl Acad Sci U S A. 2012;109:17519–24.View ArticlePubMed CentralPubMedGoogle Scholar
- Escudero M, Hipp A. Shifts in diversification rates and clade ages explain species richness in higher-level sedge taxa (Cyperaceae). Am J Bot. 2013;100:2403–11.View ArticlePubMedGoogle Scholar
- Baker WJ, Couvreur TLP. Global biogeography and diversification of palms sheds light on the evolution of tropical lineages. II. Diversification history and origin of regional assemblages. J Biogeogr. 2013;40:286–98.View ArticleGoogle Scholar
- Jansson R, Dynesius M. The fate of clades in a world of recurrent climatic change: Milankovitch oscillations and evolution. Annu Rev Ecol Syst. 2002;33:741–77.View ArticleGoogle Scholar
- Payne JL, Finnegan S. The effect of geographic range on extinction risk during background and mass extinction. Proc Natl Acad Sci U S A. 2007;104:10506–11.View ArticlePubMed CentralPubMedGoogle Scholar
- Cardinal S, Danforth BN. Bees diversified in the age of eudicots. Proc R Soc B. 2013;280:20122686.View ArticlePubMed CentralPubMedGoogle Scholar
- Vamosi JC, Otto SP, Barrett SCH. Phylogenetic analysis of the ecological correlates of dioecy in angiosperms. J Evol Biol. 2003;16:1006–18.View ArticlePubMedGoogle Scholar
- De Queiroz A. Contingent predictability in evolution: Key traits and diversification. Syst Biol. 2002;51:917–29.View ArticlePubMedGoogle Scholar
- Leslie AB, Beaulieu JM, Crane PR, Donoghue MJ. Explaining the distribution of breeding and dispersal syndromes in conifers. Proc R Soc B. 2013;280:20131812.View ArticlePubMed CentralPubMedGoogle Scholar
- Endress PK. The evolution of floral biology in basal angiosperms. Philos Trans R Soc Lond B Biol Sci. 2010;365:411–21.View ArticlePubMed CentralPubMedGoogle Scholar
- Endress PK, Doyle JA. Floral phyllotaxis in basal angiosperms: development and evolution. Curr Opin Plant Biol. 2007;10:52–7.View ArticlePubMedGoogle Scholar
- Endress PK. Himantandraceae. In: Kubitzki K, Rohwer JG, Bittrich V, editors. Fam genera vasc plants. Berlin Heidelberg New-York: Springer-Verlag; 1993. p. 338–41.Google Scholar
- Crepet WL, Nixon KC, Gandolfo MA. An extinct calycanthoid taxon, Jerseyanthus calycanthoides, from the Late Cretaceous of New Jersey. Am J Bot. 2005;92:1475–85.View ArticlePubMedGoogle Scholar
- Grant V. The pollination of Calycanthus occidentalis. Am J Bot. 1950;37:294–7.View ArticleGoogle Scholar
- Worboys SJ, Jackes BR. Pollination processes in idiospermum australiense (Calycanthaceae), an arborescent basal angiosperm of Australia’s tropical rain forests. Plant Syst Evol. 2005;251:107–17.View ArticleGoogle Scholar
- Kay KM, Sargent RD. The role of animal pollination in plant speciation: integrating ecology, geography, and genetics. Annu Rev Ecol Evol Syst. 2009;40:637–56.View ArticleGoogle Scholar
- Jones LM, Gadek PA, Harrington MG. Population genetic structuring in a rare tropical plant: Idiospermum australiense (Diels) S.T. Blake. Plant Syst Evol. 2010;286:133–9.View ArticleGoogle Scholar
- Goosem S. Update of original Wet tropics of Queensland nomination dossier. Australia: Wet Tropics Management Authority Cairns, Queensland; 2002.Google Scholar
- Edwards W, Gadek P, Weber E, Worboys S. Idiosyncratic phenomenon of regeneration from cotyledons in the idiot fruit tree, Idiospermum australiense. Austral Ecol. 2001;26:254–8.View ArticleGoogle Scholar
- Endress PK. Calycanthaceae. In: Kubitzki K, editor. The evolution of floral biolology of basal angiosperms. Berlin Heidelberg New-York: Springer-Verlag; 1993. p. 197–200.Google Scholar
- Smith SA, Donoghue MJ. Rates of molecular evolution are linked to life history in flowering plants. Science. 2008;322:86–9.View ArticlePubMedGoogle Scholar
- Heywood VH, Brummitt RK, Culham A, Seberg O. Flowering plant families of the world. Richmond: Kew Publishing; 2007.Google Scholar
- Doyle JA, Endress PK. Tracing the early evolutionary diversification of the angiosperm flower. In: Wanntorp L, Ronse De Craene LP, editors. Flowers on the tree of life. Cambridge, UK: Cambridge University Press; 2011. p. 88–119.View ArticleGoogle Scholar
- Morlon H, Parsons TL, Plotkin JB. Reconciling molecular phylogenies with the fossil record. Proc Natl Acad Sci U S A. 2011;108:16327–32.View ArticlePubMed CentralPubMedGoogle Scholar
- Rabosky DL. Automatic detection of key innovations, rate shifts, and diversity-dependence on phylogenetic trees. PLoS One. 2014;9:e89543.View ArticlePubMed CentralPubMedGoogle Scholar
- Heath TA, Huelsenbeck JP, Stadler T. The fossilized birth-death process for coherent calibration of divergence-time estimates. Proc Natl Acad Sci U S A. 2014;111:2957–66.View ArticleGoogle Scholar
- Silvestro D, Schnitzler J, Liow LH, Antonelli A, Salamin N. Bayesian estimation of speciation and extinction from incomplete fossil occurrence data. Syst Biol. 2014;63:349–67.View ArticlePubMed CentralPubMedGoogle Scholar
- Friis EM, Pedersen KR, Crane PR. Cretaceous angiosperm flowers: innovation and evolution in plant reproduction. Palaeogeogr Palaeoclimatol Palaeoecol. 2006;232:251–93.View ArticleGoogle Scholar
- Ogg JG, Hinnov LA. Cretaceous. In: Gradstein FM, Ogg JG, Schmitz M, Ogg G, editors. Geological time scale 2012. Amsterdam: Elseiver; 2012. p. 793–853.View ArticleGoogle Scholar
- Doyle JA. Molecular and fossil evidence on the origin of angiosperms. Annu Rev Earth Planet Sci. 2012;40:301–26.View ArticleGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.View ArticlePubMed CentralPubMedGoogle Scholar
- Nylander JAA. MrModeltest v2. Evolutionary Biology Centre, Uppsala Univ. 2004.Google Scholar
- Stadler T. On incomplete sampling under birth–death models and connections to the sampling-based coalescent. J Theor Biol. 2009;261:58–66.View ArticlePubMedGoogle Scholar
- Drummond AJ, Ho SYW, Phillips MJ, Rambaut A. Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006;4:e88.View ArticlePubMed CentralPubMedGoogle Scholar
- Heled J, Drummond AJ. Calibrated tree priors for relaxed phylogenetics and divergence time estimation. Syst Biol. 2011;61:138–49.View ArticlePubMed CentralPubMedGoogle Scholar
- Warnock RCM, Parham JF, Joyce WG, Lyson TR, Donoghue PCJ. Calibration uncertainty in molecular dating analyses: there is no substitute for the prior evaluation of time priors. Proc R Soc B. 2014;282:20141013.View ArticleGoogle Scholar
- Sanderson MJ. r8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics. 2003;19:301–2.View ArticlePubMedGoogle Scholar
- Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. Proc Gatew Comput Environ Work. 2010;1:1–8.Google Scholar
- Tracer, version 1.5. [http://tree.bio.ed.ac.uk/software/tracer/].
- Sanderson MJ. Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach. Mol Biol Evol. 2002;19:101–9.View ArticlePubMedGoogle Scholar
- Stamatakis A, Hoover P, Rougemont J. A rapid bootstrap algorithm for the RAxML Web servers. Syst Biol. 2008;57:758–71.View ArticlePubMedGoogle Scholar
- Paradis E, Claude J, Strimmer K. APE: Analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–90.View ArticlePubMedGoogle Scholar
- Massoni J, Couvreur TLP, Sauquet H. Chronograms of the BEAST and r8s analyses. Dryad Digital Repository. 2015. http://0-dx.doi.org.brum.beds.ac.uk/10.5061/dryad.ct231.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.