Skip to main content
  • Research article
  • Open access
  • Published:

Diversification of land plants: insights from a family-level phylogenetic analysis



Some of the evolutionary history of land plants has been documented based on the fossil record and a few broad-scale phylogenetic analyses, especially focusing on angiosperms and ferns. Here, we reconstructed phylogenetic relationships among all 706 families of land plants using molecular data. We dated the phylogeny using multiple fossils and a molecular clock technique. Applying various tests of diversification that take into account topology, branch length, numbers of extant species as well as extinction, we evaluated diversification rates through time. We also compared these diversification profiles against the distribution of the climate modes of the Phanerozoic.


We found evidence for the radiations of ferns and mosses in the shadow of angiosperms coinciding with the rather warm Cretaceous global climate. In contrast, gymnosperms and liverworts show a signature of declining diversification rates during geological time periods of cool global climate.


This broad-scale phylogenetic analysis helps to reveal the successive waves of diversification that made up the diversity of land plants we see today. Both warm temperatures and wet climate may have been necessary for the rise of the diversity under a successive lineage replacement scenario.


It is believed that climate change is one of the main factors affecting global biodiversity [13]. During the history of life, fluctuations of the world's climate have most likely caused major extinctions [4] and led to the development of new ecosystems, promoting new biotic interactions and the evolution of novel adaptive traits. The dynamics of such diversification events can be studied based on phylogenetic trees dated with fossils. Here we focus on land plants. The origin and diversification of land plants has intrigued biologists for centuries. According to the fossil record, land plants diverged from green algae before 475 million years ago (Ma; first land plant fossil) and led to the major clades found today [5, 6]. These are liverworts (74 families, ca. 6,000 spp. [7]), mosses (112 families, ca. 12,000 spp. [8, 9]), hornworts (five families, ca. 150 spp. [10]) and tracheophytes. The latter include ferns (45 families, ca. 9,000 spp. [11]), lycophytes (three families, ca. 1,200 spp. [12]), and seed plants, which in turn are separated into gymnosperms (14 families, ca. 1,000 spp. [13]) and angiosperms (456 families, ca. 260,000 spp. [13]).

There are various possible scenarios to describe the processes that influenced land plant diversification throughout geological time. One frequently proposed scenario is based on a successive replacement of ancestral lineages by more derived lineages, which in turn evolved similar habits (e.g., tree-like structure for forested ecosystems), and diversified to fill up the niches left empty after the extinction of the 'previous' taxon. In this kind of scenario, extant taxa of liverworts, mosses, and ferns, are considered to be relicts of previous radiations [14]. An alternative scenario suggests a coincidence between diversification events in each of the extant land plant lineages instead of a 'continuous replacement' idea. In this case, the majority of extant diversity is either the result of recent radiation events or of a long accumulation of species diversity throughout a taxon's history [14]. External factors, such as the break-up of continents and climate fluctuations, are prominent factors influencing the branching of the tree of life.

In this study we ask two questions: (1) Do we find evidence for non-constant rate of diversification in land plants? (2) Are major shifts of diversification rates, if any, correlated with some major external factors such as global climate warming or cooling?


We inferred the divergence times of over 98% of all families of land plants in a single phylogenetic analysis based on multiple genes from two genomes (Additional file 1; TreeBase study ID S11106). The topology and divergence times retrieved from the various analyses are broadly congruent with previous studies with limited sampling [1517]. All major lineages of land plants as well as the relationships among them were supported (bootstrap support > 74%) being mosses (62%), lycophytes as sister to seed plants (68%) and hornworts as sister to mosses (47%) the clades with lowest bootstrap values (TreeBase study ID S11106). This topology is congruent with the analysis using the three most complete markers (18S, rbcL and atpB; 6.3% missing data).

The tree was calibrated using multiple fossils. In one of the calibration procedures, we also constrained the age of angiosperms to a maximum of 130 Ma following Brenner [18] (hereafter the constrained tree). The estimated crown age of land plants was 544.7 Ma (confidence interval [C.I.] = 563.1-536.5) and that of angiosperms was 267.6 Ma (C.I. = 289.9-263.2; Additional file 2; TreeBase study ID 11106) whereas for the constrained tree we obtained a crown age for land plants of 510.8 Ma (C.I. = 512.9-475.5).

We produced lineage through time (LTT) plots for both time estimations, presented in Figure 1. These show a roughly constant rate of lineage increase (at least for the family-level studied here), although for angiosperms, ferns and mosses some acceleration is apparent since the Cretaceous, while for liverworts and gymnosperms a slowdown is observed (Figure 1).

Figure 1
figure 1

Lineage Through Time plot. Lineage Through Time (LTT) plot for liverworts (blue), mosses (green), ferns (purple), gymnosperms (yellow) and angiosperms (red) with indication of average global temperature [22] and cool climate modes (blue bars [21]). Triangles pointing up or down indicate diversification rate shifts as detected with LASER (increasing or decreasing, respectively, see Methods). (A) unconstrained tree; (B) constrained tree (i.e., angiosperms not older than 130 Ma). The y-axis indicates the number of lineages N on a logarithmic scale.

A congruent pattern is obtained when we explore the data applying a high level of background extinction using a methodology developed by Magallón & Sanderson [19]. Figure 2 shows sizes of the major clades against a 95% confidence interval of background diversification through time for land plants as a whole. In recent times, most clade sizes for mosses (Figure 2C), ferns (Figure 2D) and angiosperms (Figures 2E & 2F, the former being the tree from a constrained analysis) fall above these confidence intervals.

Figure 2
figure 2

Clade-size through time plot. Clade-size through time plot with indication of the 95% C.I. depicted from the rate of diversification of land plants as a whole and assuming some background extinction (see Methods): (A) liverworts, (B) gymnosperms, (C) mosses, (D) ferns (E) angiosperms unconstrained tree (F) angiosperms constrained tree.

Using a topology-based test of diversification [20], a total of 135 significant rate shifts were identified, with a similar figure found for the constrained tree (139; Table 1). The inclusion of 11 families with no DNA data resulted in the identification of just one more shift in diversification, i.e. on the branch leading to Balanophoraceae. We then explored the concordance of these shifts with the major cool and warm climatic modes [21] and we found some striking correlations. The majority of shifts in diversification rates in angiosperms, ferns, and mosses coincide with the last warm climate mode (Table 1). For liverworts, the highest number of shifts (5) took place in another warm climate mode (184-252 Ma; Table 1). For gymnosperms, only one shift in net diversification occurred, but in this case, during a cool period. This pattern appears to also hold if we compare the timing of shifts in diversification rates with the more continuous global temperature change presented in Scotese [22] (Additional file 3).

Table 1 Number of significant shifts in net diversification rate (topology-based method of Moore et al.[20]) with indication of cool (black) and warm (bold) climate modes of the Phanerozoic [21].

Using another diversification test that take into account branch lengths (i.e., LASER [23]), constant rates of lineage diversification were rejected for all major subclades. In gymnoperms, the best model was one with a rate shift occurring ca. 154 Ma, corresponding to a decrease in diversification during a cool climate mode (Table 2, Figure 1A). In the other subclades, two-variable rates were favoured (Table 2, Figure 1A). In angiosperms, two consecutive slowdowns in diversification were identified for the current cool climate mode. In liverworts, a similar pattern was encountered but decreases in rates of diversification occurred firstly in a cool climate mode (ca. 184 Ma) and secondly during a warm climate mode (ca. 99 Ma). In ferns and mosses, we first observe two increases in diversification during a cool climate mode (ca. 106 and 133 Ma, respectively, Table 2, Figure 1A). Subsequently, two decreases took place 60 Ma (warm mode) for ferns and 35 Ma (cool mode) for mosses (Table 2, Figure 1A). With the constrained tree, the pattern is similar for gymnosperms (Figure 1B, Additional file 4). In the case of angiosperms and liverworts, only one decrease was retrieved about 34 (cool mode) and 99 Ma (warm mode) ago, respectively (Figure 1B, Additional file 4). This pattern is similar to that obtained in the unconstrained tree for mosses and ferns (Figure 1B, Additional file 4), although this time fern diversification increases during a warm mode (ca. 93 Ma) and decreases during a cool mode (ca. 52 Ma; Figure 1B, Additional file 4).

Table 2 LASER analysis using a constant-rate birth-death model with no extinction (a = 0) against variable-rates models with 2 and 3 rates (r) and 1 or 2 time shifts given for best fitting model (ts; time unit is million years ago).

Finally, diversification test incorporating multiple birth and death models as implemented in MEDUSA [24] located 69 diversification rate shifts being the highest overall net diversification rates for different clades within angiosperms (Figure 3, Additional file 5). Among land plants we also found rate shifts leading to high clade-size in mosses and ferns for individual families and clades (Figure 3, Additional file 5). On other hand rates among liverworts and gymnosperms were among the lowest: their background rate were similar to the overall background rate, and their highest rates were lower than most rates found for mosses and ferns (see Additional file 5). Results using the constrained tree were widely congruent (Additional file 5) and a new rate shift for gymnosperms and higher net diversification rates across monocots were recovered.

Figure 3
figure 3

MEDUSA chronogram. Diversification chronogram with rate shifts located using MEDUSA [24] for different groups of land plants. Numbers correspond to the rate shifts located by MEDUSA being the numbers in increasing order from the highest to lowest net diversification rate (see Additional file 5). Different colours indicate different net diversification rates found in the tree. Boostrap support for the main nodes are indicated with one asterisk (> 70%) and two asterisk (50%-70%). Angiosperms (grey asterisk) have been simplified for this figure.

Discussion and Conclusion

By combining data for all families of land plants we are now able to clarify the picture of their evolution through geological times. Lineages of extant gymnosperms radiated in the Permian and experienced a decrease in diversification rate towards the end of the Jurassic (analysis with unconstrained tree) or early Cretaceous (analysis with constrained tree), during a cool climate mode. Although their early history may have involved various lineage replacements associated with the evolution of new ecosystems [25, 26], we found that the slowdown in diversification of gymnosperms took place in the same period as liverworts while mosses were diversifying intensely, pointing towards a role of climate in determining such patterns.

In this study we were also able to evaluate the diversification dynamics of all families of mosses within a phylogenetic framework for the first time. Our analyses converge to show that the diversification rate of this group experienced an important acceleration in the Cretaceous, potentially 'replacing' the diversity of gymnosperms and liverworts. This occurred during a warm climate mode when tropical habitats were undergoing expansion. Significantly, it also corresponds to the origin of the angiosperms according to Brenner [18], or to the origin of major groups of angiosperms (asterids, rosids) as found on our unconstrained analysis and as suggested by previous studies [27]. In this sense, mosses have diversified at the same period as the one reported for fern as "the shadow of angiosperms". It is important to note that the Cretaceous could be divided in three main intervals with regards to vegetation and climate: i) "Early" Cretaceous (ca. Berriasian--Barremian) with few angiosperms, probably no closed canopy angiosperm forests, largely dry climates at low palaeolatitudes; ii) mid-Cretaceous (ca. Aptian-Santonian/Campanian), where we observe the rapid diversification of angiosperms, with presence of some angiosperm-dominated forests but still no tropical everwet forests at low palaeolatitudes; and iii) "Late" Cretaceous (ca. Campanian-Maastrichtian), where we see an early development of angiosperm dominated forests, possibly with everwet forests in low palaeolatitudes of the Old World, and perhaps also in the New World [14]. This is then followed by iv) the Early Cenozoic, when temperatures were warm and climate wet--and where there is strong evidence of widespread tropical-sub-tropical warm wet forests [14]. We found that all six shifts in diversification for mosses (Table 1) fall within this last two intervals (i.e. 75, 69, 64, 57, 43 and 37 mya; see details in Additional files 3), pointing to the importance of both warm temperatures and wet climate for the rise of moss diversity.

According to our analyses, mosses were not the only group to have diversified in the shadow of angiosperms: ferns have also radiated in a period that coincides with the rise of angiosperms. Such a pattern, had previously been reported [28]. Here, we find further support for this hypothesis of diversification in the shadow of angiosperms, identifying a significant increase in diversification during the warmest period of the Cretaceous, and decrease during the coolest period of the Tertiary (see Figure 1). More specifically, three out of the five rate shifts (Table 1) fall within interval iii) above of the Cretaceous, when climate was warm and wet.

Finally, we found that angiosperm diversity has accumulated sharply in recent time (as shown by the LTT plots), but diversification decreased in the coolest period of the Tertiary (Figure 1). This is in agreement with the idea that angiosperms have outcompeted and outnumber gymnosperms and free-sporing plants [29, 30]. Subsequently, ferns (especially polypods [28, 31]) and mosses [32] opportunistically diversified in the ecological niches provided by the angiosperms as the climate became warmer and more humid. In this sense, our study favours the "successive replacement" of ancestral lineages [14].



We put together phylogenetic data for at least one representative of each of the 706 currently accepted families of land plants (Additional file 1). Our dataset was assembled using plastid rbcL, atpB and rps4 genes, as well as 18S and 26S nuclear ribosomal regions (hereafter 18S rDNA and 26S rDNA). We downloaded sequences from GenBank when available and filled some of the gaps by sequencing missing taxa when we were able to obtain suitable material (Additional file 1). DNA extraction and PCR amplification used standard protocols and primers for nuclear and plastid genes from Nickrent and Starr [33] and Cox et al. [34]. We sequenced the 18S rDNA for 22 angiosperms and 13 mosses, rbcL for two angiosperms, 10 mosses and one liverwort, and atpB for 39 mosses, two liverworts, one hornwort, and 18 angiosperms (Additional file 1). In total we produced a 6,950 base pairs data matrix consisting of 699 families (including four outgroups) with 65% of data presence. Only one gene could be obtained for 55 of these 699 families (Additional file 1). Streptophytes and Chlorophytes were used as the outgroup.

Due to the large size of the matrix, maximum likelihood analyses were performed in RAxML [35] using 200 bootstrap replicates and GTR+GAMMA model, as selected by ModelTest [36]. Divergence times were calculated using penalized likelihood in r8s [37] and the smoothing parameter (smooth = 1000) was calculated by cross-validation. We calibrated the chronogram with the age of eudicots at 121 mya, corresponding to the appearance of the tricolpate pollen grain typical of this clade [38]. We used a further sixteen calibration points as minimum constrains, plus a maximum age of 725 Ma [39] for the root of the tree (Marchantiopsida, Monilophytes, Mosses, Seed plants, Annonaceae, Calycanthaceae, Hedyosmum, Lauraceae, Magnoliaceae, Meliosma, Menispermaceae, Nelumbaceae, Nymphaceae, Platanaceae, Trochodendron, Winteraceae, Additional file 6). Confidence intervals (C.I.) for divergences times were calculated by repeating the dating procedures in r8s using 100 bootstrapped matrices produced in RAxML [35]. The dating procedure was repeated constraining the age of angiosperms to a maximum of 130 Ma following Brenner [18], i.e. "constrained tree".

Diversification tests

We examined diversification through time using several methods.

Firstly, we plotted the number of lineages through time (hereafter LTT plots) for each major subclade of land plants using the APE 1.8 package [40].

Secondly, to take into account extinction rates we used the approach of Magallón and Sanderson [19]. For time intervals of one million year, we calculated net diversification rates under a relative high level of background extinction (0.9 using equation 10 of Magallón and Sanderson [19].

Thirdly, we applied a topological-based test of diversification. Diversification rate shifts were calculated using the Δ1 statistics of Moore et al. [20] as implemented by Bouchenak-Khelladi et al. [41] in ApTreeshape [42] using 0.05 significance level as the cut-off point. This test uses the tree but also takes into account the total number of species per family (Table 1).

Fourthly, we used a test of diversification that takes into account branch lengths, i.e. the elapsed time between the nodes of the family-level tree, LASER [23]. Using the Akaike Information Criterion (AIC), LASER can compare models with various rates of diversifications (yule model with rates r) against the null expectation of a constant rate (birth-death model with no extinction). LASER also allows to identify at which points in time a given rate shift occurred (ts). LASER was applied to all major subclades.

Fifthly, we tested for multiple shifts in birth and death rates using a stepwise approach implemented in MEDUSA until improvement in AIC score was < 4 [24]. Net diversification rates together with relative extinction rates and AIC improvements were retrieved (Additional file 5).

Also, to comply with other phylogenetic analyses that have combined more genes but for fewer taxa, we also re-ran the analyses above with the following two modifications. First 11 families for which we could not obtain any DNA data (i.e., five families of liverworts, three of mosses, and three of angiosperms; Additional file 1) were placed in the DNA-based phylogenetic tree using taxonomic information following Crosby et al. [8], Buck and Goffinet [9], Stevens [13], Heinrichs et al. [7] and Smith et al. [11] (see Additional file 2). Although this procedure is suboptimal, it allowed us to perform diversification tests on a complete-family level tree. Second we enforced hornworts and lycophytes to be sister to vascular plants [15, 16] plus we set the maximum age for angiosperms to 130 Ma (following Brenner [18]; our "constrained tree"). Results were compared for the constrained vs. unconstrained topologies. Finally, we compared these diversification profiles and metrics against the distribution of the climate modes of the Panerozoic following Frakes et al. [21], as well as the global temperature model of Scotese [22].


  1. Pounds JA, Fogden MPL, Campbell JH: Biological response to climate change on a tropical mountain. Nature. 1999, 398: 611-615. 10.1038/19297.

    Article  CAS  Google Scholar 

  2. Parmesan C, Yohe G: A globally coherent fingerprint of climate change impacts across natural systems. Nature. 2003, 421: 37-42. 10.1038/nature01286.

    Article  CAS  PubMed  Google Scholar 

  3. Root TL, Price JT, Hall KR, Schneider H, Rosenzweig C, Pounds JA: Fingerprints of global warming on wild animals and plants. Nature. 2003, 421: 57-60. 10.1038/nature01333.

    Article  CAS  PubMed  Google Scholar 

  4. Wignall PB: Causes of mass extinctions. Extinctions in the history of life. Edited by: Taylor PD. 2004, Cambridge University Press, 151-174.

    Google Scholar 

  5. Wellman CH, Osterhoff PL, Mohiuddin U: Fragments of the earliest land plants. Nature. 2003, 425: 282-285. 10.1038/nature01884.

    Article  CAS  PubMed  Google Scholar 

  6. Steemans P, Le Herisse A, Melvin J, Miller MA, Paris F, Verniers J, Wellmann CH: Origin and radiation of the earliest vascular land plants. Science. 2009, 324: 353-10.1126/science.1169659.

    Article  CAS  PubMed  Google Scholar 

  7. Heinrichs J, Gradstein SR, Wilson R, Schneider H: Towards a natural classification of liverworts (Marchantiophyta) based on the chloroplast gene rbcL. Crypt Bryol. 2005, 26: 131-150.

    Google Scholar 

  8. Crosby MR, Magill RE, Allen B, He S: A Checklist of the Mosses. 1999, []

    Google Scholar 

  9. Buck WR, Goffinet B: Morphology and classification of mosses. Bryophyte Biology. Edited by: Shaw AJ, Goffinet B. 2000, Cambridge: Cambridge University Press, 71-123.

    Chapter  Google Scholar 

  10. Shaw J, Renzaglia K: Phylogeny and diversification of bryophytes. Am J Bot. 2004, 91: 1557-1581. 10.3732/ajb.91.10.1557.

    Article  PubMed  Google Scholar 

  11. Smith AR, Pryer KM, Schuettpelz E, Korall P, Schneider H, Wolf PG: A classification of extant ferns. Taxon. 2006, 55: 705-731. 10.2307/25065646.

    Article  Google Scholar 

  12. Starr A, Evers CA, Starr L: Biology: Concepts and Applications Without Physiology. 2008, New York: Brooks/Cole Publishing

    Google Scholar 

  13. Stevens PF: Angiosperm phylogeny website, v8. 2001, []

    Google Scholar 

  14. Kenrick P, Crane PR: The origin and early evolution of plants on land. Nature. 1997, 389: 33-39. 10.1038/37918.

    Article  CAS  Google Scholar 

  15. Pryer KM, Schneider H, Smith AR, Cranfill R, Wolf PG, Hunt JS, Sipes SD: Horsetails and ferns are a monophyletic group and the closest living relatives to seed plants. Nature. 2001, 409: 618-622. 10.1038/35054555.

    Article  CAS  PubMed  Google Scholar 

  16. Qiu Y, et al: The deepest divergences in land plants inferred from phylogenomic evidence. Proc Natl Acad Sci USA. 2006, 103: 15511-15516. 10.1073/pnas.0603335103.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Wang H, Moore MJ, Soltis PS, Bell CD, Brockington SF, Alexandre R, Davis CC, Latvis M, Manchester SR, Soltis DE: Rosid radiation and the rapid rise of angiosperm-dominated forests. Proc Natl Acad Sci USA. 2009, 106: 3853-3858. 10.1073/pnas.0813376106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Brenner GJ: Evidence for the earliest stage of angiosperm pollen evolution: A paleoequatorial section from Israel. Flowering plant origin, evolution and phylogeny. Edited by: Taylor DW, Hickey LJ. 1996, New York: Chapman & Hall, 91-115.

    Chapter  Google Scholar 

  19. Magallón S, Sanderson MJ: Absolute diversification rates in angiosperm clades. Evolution. 2001, 55: 1762-1780.

    Article  PubMed  Google Scholar 

  20. Moore BR, Chan KMA, Donoghue MJ: Detecting diversification rate variation in supertrees. Phylogenetic Supertrees: Combining Information to Reveal the Tree of Life. Edited by: Bininda-Emonds ORP. 2004, Dordrecht: Netherlands Kluwer Academic, 487-533.

    Chapter  Google Scholar 

  21. Frakes L, Francis JE, Syktus JI: Climate modes of the phanerozoic. 1992, Cambridge: Cambridge Univ. Press

    Book  Google Scholar 

  22. Scotese CR: PALEOMAP Project. 2001, []

    Google Scholar 

  23. Rabosky DL: LASER: A Maximum Likelihood Toolkit for Detecting Temporal Shifts in Diversification Rates From Molecular Phylogenies. Evol Bioinf Online. 2006, 2: 257-260.

    CAS  Google Scholar 

  24. Alfaro ME, Santini F, Brock C, Alamillo H, Dornburg A, Rabosky DL, Carnelale G, Harmon LJ: Nine exceptional radiations plus high turnover explain species diversity in jawed vertebrates. Proc Natl Acad Sci USA. 2009, 106: 13410-13414. 10.1073/pnas.0811087106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Looy CV, Brugman WA, Dilcher DL, Visscher H: The delayed resurgence of equatorial forests after the Permian-Triassic ecologic crisis. Proc Natl Acad Sci USA. 1999, 96: 13857-13862. 10.1073/pnas.96.24.13857.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Looy CV, Twitchett RJ, Dilcher DL, Van Konijnenburg-Van JHA, Visscher H: Life in the end-Permian dead zone. Proc Natl Acad Sci USA. 2001, 98: 7879-7883. 10.1073/pnas.131218098.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Bell CD, Soltis DE, Soltis PS: The age and diversification of the angiosperms re-revisited. Am J Bot. 2010, 97: 1296-1303. 10.3732/ajb.0900346.

    Article  PubMed  Google Scholar 

  28. Schneider H, Schuettpelz E, Pryer KM, Cranfill R, Magallón S, Lupia R: Ferns diversified in the shadow of the angiosperms. Nature. 2004, 428: 553-557. 10.1038/nature02361.

    Article  CAS  PubMed  Google Scholar 

  29. Burger WC: Why are there so many kinds of flowering plants?. Bioscience. 1981, 31: 572-581. 10.2307/1308218.

    Article  Google Scholar 

  30. Magallón S, Castillo A: Angiosperm diversification through time. Am J Bot. 2009, 96: 349-365. 10.3732/ajb.0800060.

    Article  PubMed  Google Scholar 

  31. Schuettpelz E, Pryer KM: Evidence for a Cenozoic radiation of ferns in an angiosperm-dominated canopy. Proc Natl Acad Sci USA. 2009, 106: 11200-11205. 10.1073/pnas.0811136106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Newton AE, Wikstrom N, Bell N, Forrest NL, Ignatov MS: Dating the Diversification of the Pleurocarpous Mosses. Pleurocarpous Mosses: Systematics and Evolution. Edited by: Newton AE, Tangney RS. 2007, CRC Press, 337-366.

    Chapter  Google Scholar 

  33. Nickrent D, Starr EM: High rates of nucleotide substitution in nuclear small-subunit (18S) rDNA from holoparasitic flowering plants. J Mol Evol. 1994, 39: 62-70.

    Article  CAS  PubMed  Google Scholar 

  34. Cox C, Goffinet B, Newton AE, Shaw AJ, Hedderson TAJ: Phylogenetic relationships among the diplolepidous-alternate mosses (Bryidae) inferred from nuclear and chloroplast DNA sequences. Bryologist. 2000, 103: 224-241. 10.1639/0007-2745(2000)103[0224:PRATDA]2.0.CO;2.

    Article  CAS  Google Scholar 

  35. Stamatakis A, Hoover P, Rougemont J: A Rapid Bootstrap Algorithm for the RAxML Web-Servers. Syst Biol. 2008, 75: 758-771.

    Article  Google Scholar 

  36. Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.

    Article  CAS  PubMed  Google Scholar 

  37. Sanderson MJ: r8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics. 2003, 19: 301-303. 10.1093/bioinformatics/19.2.301.

    Article  CAS  PubMed  Google Scholar 

  38. Doyle JA, Hotton CL: Diversity of early angiosperm pollen in a cladistic context. Pollen and spores. Edited by: Blackmore S, Barnes A. 1991, Oxford: Clarendon Press, 169-195.

    Google Scholar 

  39. Zimmer A, Lang D, Richardt S, Frank W, Reski R, Rensing R: Dating the early evolution of plants: detection and molecular clock analyses of orthologs. Mol Genet Genom. 2007, 278: 393-402. 10.1007/s00438-007-0257-6.

    Article  CAS  Google Scholar 

  40. Paradis E, Claude J, Strimmer K: APE: Analyses of Phylogenetics and Evolution in R language. Bioinformatics. 2004, 20: 289-290. 10.1093/bioinformatics/btg412.

    Article  CAS  PubMed  Google Scholar 

  41. Bouchenak-Khelladi Y, Verboom GA, Hodkinson TR, Salamin N, Francois O, Chonghaile GN, Savolainen V: The origins and diversification of C4 grasses and savanna-adapted ungulates. Global Change Biol. 2009, 1-21.

    Google Scholar 

  42. Bortolussi N, Durand E, Blum MGB, Francois O: ApTreeshape: a statistical analysis of phylogenetic treeshape. Bioinformatics. 2006, 22: 363-364. 10.1093/bioinformatics/bti798.

    Article  CAS  PubMed  Google Scholar 

  43. Stotler RE, Crandall-Stotler BJ: Bryophytes. 2011, []

    Google Scholar 

  44. Chesters KIM: Some plant remains from the Upper Cretaceous and Tertiary of West Africa. Ann Mag Nat Hist. 1955, 12: 498-504.

    Article  Google Scholar 

  45. Crane PR, Friis EM, Pedersen KR: Palaeobotanical evidence on the early radiation of magnoliid angiosperms. Plant Syst Evol. 1994, 8: 51-72.

    Google Scholar 

  46. Dilcher DL, Crane PR: Archaeanthus: an early angiosperm from the Cenomanian of the western interior of North America. Ann Mo Bot Gard. 1984, 71: 351-383. 10.2307/2399030.

    Article  Google Scholar 

  47. Magallon S, Crane PR, Herendeen PS: Phylogenetic pattern, diversity, and diversification of eudicots. Ann Mo Bot Gard. 1999, 86: 297-372. 10.2307/2666180.

    Article  Google Scholar 

  48. Friis EM, Eklund H, Pedersen KR, Crane PR: Virgininanthus calycanthoides gen. et sp. nov.: a calycanthaceous flower from the Potomac Group (Early Cretaceous) of eastern North America. Int J Plant Sci. 1994, 155: 772-785. 10.1086/297217.

    Article  Google Scholar 

  49. Friis EM, Pedersen KR, Crane PR: Angiosperm floral structures from the Early Cretaceous of Portugal. Plant Syst Evol. 1994, 8: 31-49.

    Google Scholar 

  50. Friis EM, Pedersen KR, Crane PR: Fossil evidence of water lilies (Nymphaeales) in the Early Cretaceous. Nature. 2001, 410: 356-360.

    Article  Google Scholar 

  51. Manchester SR, Crane PR, Dilcher DL: Nordenskioldia and Trochodendron (Trochodendraceae) from the Miocene of Northwestern North America. Bot Gazette. 1991, 152: 357-368. 10.1086/337898.

    Article  Google Scholar 

  52. Doyle JA, Hotton CL, Ward JV: Early Cretaceous tetrads, zonasulcate pollen, and Winteraceae. I. Taxonomy, morphology, and ultrastructure. Am J Bot. 1990, 77: 1544-1557. 10.2307/2444487.

    Article  Google Scholar 

  53. Doyle JA: Early Cretaceous tetrads, zonasulcate pollen, and Winteraceae. II. Cladistic analysis and implications. Am J Bot. 1990, 77: 1558-1568. 10.2307/2444488.

    Article  Google Scholar 

  54. Heinrichs J, Hentschel J, Wilson R, Feldberg K, Schneider H: Evolution of leafy liverworts (Jungermanniidae, Marchantiophyta): estimating divergence times from chloroplast DNA sequences using penalized likelihood with integrated fossil evidence. Taxon. 2007, 56: 31-44.

    Google Scholar 

  55. Bateman RM: Palaeobiology and phylogenetic implications of anatomically preserved Archaeocalamites from the Dinantian of Oxroad Bay and Loch Humphrey Bun, Scotland. Palaeontographica Abt B Palaephytol. 1991, 223: 1-59.

    Google Scholar 

  56. Miller CN: Implications of fossil conifers for the phylogenetic relationships of living families. Bot Rev. 1999, 30: 239-277.

    Article  Google Scholar 

  57. Newton AE, Wikström N, Bell N, Forrest LL, Ignatov MS: Dating the diversification of the pleurocarpous mosses. Pleurocarpous mosses: Systematics and Evolution. Edited by: Newton A, Tangney N. 2006, CRC Press, Boca Raton, Systematics Association, 337-366.

    Google Scholar 

Download references


We thank P. Crane for early review of the manuscript and helpful information, A. Newton for plant material, A. Stamatakis, Y. Bouchenak-Khelladi and P. Ajawatanawong for their help with bioinformatics, C. Cox and N. Bell for providing 18S sequences of mosses, M. Grundmann, A. Navarro-Gomez and S. Russel for laboratory assistance and A. Papadoupulos, J. Schnitzler, P. Rymer, L. Valente, J. Hortal, M. Romeralo and the editor and three anonymous reviewers for helpful comments. This study was supported by the Spanish Ministerio de Ciencia and the Bentham-Moxon Trust (to O.F.P), and Leverhulme Trust, NERC, Royal Society and European Research Council (to V.S).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Omar Fiz-Palacios.

Additional information

Authors' contributions

OFP carried out the lab work, assembled the data and performed the analyses. OFP, HS and VS designed the study. JH provided crucial material and information. OFP, HS and VS coordinated the study and wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Taxa and GenBank accession numbers with new sequences generated for this study in bold. Families for which DNA data could not be obtained are indicated in italics. (XLS 200 KB)


Additional file 2: Chronogram of the unconstrained tree. Numbers after family name are species number considered for our analysis following Stevens [13] for angiosperms and gymnosperms, Crosby et al. [8] and Buck and Goffinet [9] for mosses, Smith et al. [11] for ferns and Stotler et al. [43] for liverworts. The x axis indicate time in million years. The placement of 11 families in which no molecular data could be collected are indicated; they were connected to node numbers as follows: 1 - Monocarpaceae, 2 - Sandeothallaceae, 3 - Chonecoleaceae, 4 - Grolleaceae, 5 - Trichotemnomaceae, 6 - Viridivelleraceae, 7 - Microtheliaceae, 8 - Sorapillaceae, 9 - Hapthantaceae, 10 - Balanophoraceae and 11 - Rafflesiaceae. (PDF 6 MB)


Additional file 3: Number of significant shift in net diversification rate (topology-based method of Moore et al. [20]) with indication of cool (grey) and warm (black) temperatures of the Phanerozoic [22]. Bold number columns correspond to periods of rate increases. Time unit is Mya (million years ago). For each lineage the bottom row (italics) corresponds to the constrained tree. (DOC 50 KB)


Additional file 4: LASER analysis for the constrained tree using a constant-rate birth-death model with no extinction (a = 0) against variable-rates models with 2 and 3 rates (r) and 1 or 2 time shifts given for best fitting model (ts; time unit is million years ago). (DOC 61 KB)


Additional file 5: Diversification rate shift retrieved from MEDUSA [24]. Numbers on the first column correspond to the net diversification rate from highest to lowest and are depicted on Figure 3. Non-angiosperms cases are highlighted in bold, followed by the name of the group they belong to in brackets. "r" are the estimates for net diversification rate, "e" are estimates for relative extinction rate, "ΔAIC" is the increase on the stepwise AIC procedure and "ΔAICc" is the increase when corrected for small sample size [24]. Results for the constrained tree are presented at the bottom. (XLS 40 KB)


Additional file 6: Minimum-age calibration points used in divergence time reconstructions [4457]. (DOC 58 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Fiz-Palacios, O., Schneider, H., Heinrichs, J. et al. Diversification of land plants: insights from a family-level phylogenetic analysis. BMC Evol Biol 11, 341 (2011).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: