- Research article
- Open Access
Diversification of land plants: insights from a family-level phylogenetic analysis
BMC Evolutionary Biologyvolume 11, Article number: 341 (2011)
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 [1–3]. During the history of life, fluctuations of the world's climate have most likely caused major extinctions  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. ), mosses (112 families, ca. 12,000 spp. [8, 9]), hornworts (five families, ca. 150 spp. ) and tracheophytes. The latter include ferns (45 families, ca. 9,000 spp. ), lycophytes (three families, ca. 1,200 spp. ), and seed plants, which in turn are separated into gymnosperms (14 families, ca. 1,000 spp. ) and angiosperms (456 families, ca. 260,000 spp. ).
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 . 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 . 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 [15–17]. 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  (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).
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 . 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.
Using a topology-based test of diversification , 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  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  (Additional file 3).
Using another diversification test that take into account branch lengths (i.e., LASER ), 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).
Finally, diversification test incorporating multiple birth and death models as implemented in MEDUSA  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.
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 , or to the origin of major groups of angiosperms (asterids, rosids) as found on our unconstrained analysis and as suggested by previous studies . 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 . 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 . 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 . 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  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 .
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  and Cox et al. . 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  using 200 bootstrap replicates and GTR+GAMMA model, as selected by ModelTest . Divergence times were calculated using penalized likelihood in r8s  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 . We used a further sixteen calibration points as minimum constrains, plus a maximum age of 725 Ma  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 . The dating procedure was repeated constraining the age of angiosperms to a maximum of 130 Ma following Brenner , i.e. "constrained tree".
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 .
Secondly, to take into account extinction rates we used the approach of Magallón and Sanderson . 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 .
Thirdly, we applied a topological-based test of diversification. Diversification rate shifts were calculated using the Δ1 statistics of Moore et al.  as implemented by Bouchenak-Khelladi et al.  in ApTreeshape  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 . 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 . 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. , Buck and Goffinet , Stevens , Heinrichs et al.  and Smith et al.  (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 ; 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. , as well as the global temperature model of Scotese .
Pounds JA, Fogden MPL, Campbell JH: Biological response to climate change on a tropical mountain. Nature. 1999, 398: 611-615. 10.1038/19297.
Parmesan C, Yohe G: A globally coherent fingerprint of climate change impacts across natural systems. Nature. 2003, 421: 37-42. 10.1038/nature01286.
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.
Wignall PB: Causes of mass extinctions. Extinctions in the history of life. Edited by: Taylor PD. 2004, Cambridge University Press, 151-174.
Wellman CH, Osterhoff PL, Mohiuddin U: Fragments of the earliest land plants. Nature. 2003, 425: 282-285. 10.1038/nature01884.
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.
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.
Crosby MR, Magill RE, Allen B, He S: A Checklist of the Mosses. 1999, [http://www.mobot.org/MOBOT/tropicos/most/checklist.shtml]
Buck WR, Goffinet B: Morphology and classification of mosses. Bryophyte Biology. Edited by: Shaw AJ, Goffinet B. 2000, Cambridge: Cambridge University Press, 71-123.
Shaw J, Renzaglia K: Phylogeny and diversification of bryophytes. Am J Bot. 2004, 91: 1557-1581. 10.3732/ajb.91.10.1557.
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.
Starr A, Evers CA, Starr L: Biology: Concepts and Applications Without Physiology. 2008, New York: Brooks/Cole Publishing
Stevens PF: Angiosperm phylogeny website, v8. 2001, [http://www.mobot.org/MOBOT/research/APweb]
Kenrick P, Crane PR: The origin and early evolution of plants on land. Nature. 1997, 389: 33-39. 10.1038/37918.
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.
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.
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.
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.
Magallón S, Sanderson MJ: Absolute diversification rates in angiosperm clades. Evolution. 2001, 55: 1762-1780.
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.
Frakes L, Francis JE, Syktus JI: Climate modes of the phanerozoic. 1992, Cambridge: Cambridge Univ. Press
Scotese CR: PALEOMAP Project. 2001, [http://www.scotese.com/Default.htm]
Rabosky DL: LASER: A Maximum Likelihood Toolkit for Detecting Temporal Shifts in Diversification Rates From Molecular Phylogenies. Evol Bioinf Online. 2006, 2: 257-260.
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.
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.
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.
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.
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.
Burger WC: Why are there so many kinds of flowering plants?. Bioscience. 1981, 31: 572-581. 10.2307/1308218.
Magallón S, Castillo A: Angiosperm diversification through time. Am J Bot. 2009, 96: 349-365. 10.3732/ajb.0800060.
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.
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.
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.
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.
Stamatakis A, Hoover P, Rougemont J: A Rapid Bootstrap Algorithm for the RAxML Web-Servers. Syst Biol. 2008, 75: 758-771.
Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.
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.
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.
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.
Paradis E, Claude J, Strimmer K: APE: Analyses of Phylogenetics and Evolution in R language. Bioinformatics. 2004, 20: 289-290. 10.1093/bioinformatics/btg412.
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.
Bortolussi N, Durand E, Blum MGB, Francois O: ApTreeshape: a statistical analysis of phylogenetic treeshape. Bioinformatics. 2006, 22: 363-364. 10.1093/bioinformatics/bti798.
Stotler RE, Crandall-Stotler BJ: Bryophytes. 2011, [http://bryophytes.plant.siu.edu/]
Chesters KIM: Some plant remains from the Upper Cretaceous and Tertiary of West Africa. Ann Mag Nat Hist. 1955, 12: 498-504.
Crane PR, Friis EM, Pedersen KR: Palaeobotanical evidence on the early radiation of magnoliid angiosperms. Plant Syst Evol. 1994, 8: 51-72.
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.
Magallon S, Crane PR, Herendeen PS: Phylogenetic pattern, diversity, and diversification of eudicots. Ann Mo Bot Gard. 1999, 86: 297-372. 10.2307/2666180.
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.
Friis EM, Pedersen KR, Crane PR: Angiosperm floral structures from the Early Cretaceous of Portugal. Plant Syst Evol. 1994, 8: 31-49.
Friis EM, Pedersen KR, Crane PR: Fossil evidence of water lilies (Nymphaeales) in the Early Cretaceous. Nature. 2001, 410: 356-360.
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.
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.
Doyle JA: Early Cretaceous tetrads, zonasulcate pollen, and Winteraceae. II. Cladistic analysis and implications. Am J Bot. 1990, 77: 1558-1568. 10.2307/2444488.
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.
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.
Miller CN: Implications of fossil conifers for the phylogenetic relationships of living families. Bot Rev. 1999, 30: 239-277.
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.
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).
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.