Reconciling fossils with phylogenies reveals the origin and macroevolutionary processes explaining the global cycad biodiversity

Accounting for sampling heterogeneity suggests a low paleolatitude origin for dinosaurs

Dinosaurs dominated Mesozoic terrestrial ecosystems for ∼160 million years, but their biogeographic origin remains poorly understood. The earliest unequivocal dinosaur fossils appear in the Carnian (∼230 Ma) of southern South America and Africa, leading most authors to propose southwestern Gondwana as the likely center of origin. However, the high taxonomic and morphological diversity of these earliest assemblages suggests a more ancient evolutionary history that is currently unsampled. Phylogenetic uncertainty at the base of Dinosauria, combined with the subsequent appearance of dinosaurs throughout Laurasia in their early evolutionary history, further complicates this picture. Here, we estimate the distribution of early dinosaurs and their archosaurian relatives under a phylogenetic maximum likelihood framework, testing alternative topological arrangements and incorporating potential abiotic barriers to dispersal into our biogeographic models. For the first time, we include spatiotemporal sampling heterogeneity in these models, which frequently supports a low-latitude Gondwanan origin for dinosaurs. These results are best supported when silesaurids are constrained as early-diverging ornithischians, which is likely because this topology accounts for the otherwise substantial ornithischian ghost lineage, explaining the group's absence from the fossil record prior to the Early Jurassic. Our results suggest that the archosaur radiation also took place within low-latitude Gondwana following the end-Permian extinction before lineages dispersed across Pangaea into ecologically and climatically distinct provinces during the Late Triassic. Mesozoic terrestrial vertebrates are under-sampled at low paleolatitudes, and our findings suggest that heterogeneous sampling has hitherto obscured the true paleobiogeographic origin of dinosaurs and their kin.

Bayesian tip-dated timeline for diversification and major biogeographic events in Muroidea (Rodentia), the largest mammalian radiation

BACKGROUND

Extinct organisms provide vital information about the time of origination and biogeography of extant groups. The development of phylogenetic methods to study evolutionary processes through time has revolutionized the field of evolutionary biology and led to an unprecedented expansion of our knowledge of the tree of life. Recent developments applying Bayesian approaches, using fossil taxa as tips to be included alongside their living relatives, have revitalized the use of morphological data in evolutionary tree inferences. Eumuroida rodents represent the largest group of mammals including more than a quarter of all extant mammals and have a rich fossil record spanning the last ~ 45 million years. Despite this wealth of data, our current understanding of the classification, major biogeographic patterns, and divergence times for this group comes from molecular phylogenies that use fossils only as a source of node calibrations. However, node calibrations impose several constraints on tree topology and must necessarily make a priori assumptions about the placement of fossil taxa without testing their placement in the tree.

RESULTS

We present the first morphological dataset with extensive fossil sampling for Muroidea. By applying Bayesian morphological clocks with tip dating and process-based biogeographic models, we provide a novel hypothesis for muroid relationships and revised divergence times for the clade that incorporates uncertainty in the placement of all fossil species. Even under strong violation of the clock model, we found strong congruence between results for divergence times, providing a robust timeline for muroid diversification. This new timeline was used for biogeographic analyses, which revealed a dynamic scenario mostly explained by dispersal events between and within the Palearctic and North African regions.

CONCLUSIONS

Our results provide important insights into the evolution of Muroidea rodents and clarify the evolutionary pathways of their main lineages. We exploited the advantage of tip dating Bayesian approaches in morphology-based datasets and provided a classification of the largest superfamily of mammals resulting from robust phylogenetic inference, inferring the biogeographical history, diversification, and divergence times of its major lineages.

JC. A possible case of adaptive radiation in cycads. A commentary on 'Transcriptome sequencing data provide a solid base to understand the phylogenetic relationships, biogeography and reticulated evolution of the genus Zamia L. (Cycadales: Zamiaceae)'

This article comments on:

Anders Lindstrom, Sadaf Habib, Shanshan Dong, Yiqing Gong, Jian Liu, Michael Calonje, Dennis Stevenson and Shouzhou Zhang, Transcriptome sequencing data provide a solid base to understand the phylogenetic relationships, biogeography and reticulated evolution of the genus Zamia L. (Cycadales: Zamiaceae), Annals of Botany, Volume 134, Issue 5, 1 November 2024, Pages 747–768 https://doi.org/10.1093/aob/mcae065

Transcriptome sequencing data provide a solid base to understand the phylogenetic relationships, biogeography and reticulated evolution of the genus Zamia L. (Cycadales: Zamiaceae)

BACKGROUND AND AIMS

Cycads are a key lineage to understand the early evolution of seed plants and their response to past environmental changes. However, tracing the evolutionary trajectory of cycad species is challenging when the robust relationships at inter- or infrageneric level are not well resolved.

METHODS

Here, using 2901 single-copy nuclear genes, we explored the species relationships and gene flow within the second largest genus of cycads, i.e. Zamia, based on phylotranscriptomic analyses of 90 % extant Zamia species. Based on a well-resolved phylogenetic framework, we performed gene flow analyses, molecular dating and biogeographical reconstruction to examine the spatiotemporal evolution of Zamia. We also performed ancestral state reconstruction of a total of 62 traits of the genus to comprehensively investigate its morphological evolution.

KEY RESULTS

Zamia comprises seven major clades corresponding to seven distinct distribution areas in the Americas, with at least three reticulation nodes revealed in this genus. Extant lineages of Zamia initially diversified around 18.4-32.6 (29.14) million years ago in Mega-Mexico, and then expanded eastward into the Caribbean and southward into Central and South America. Ancestral state reconstruction revealed homoplasy in most of the morphological characters.

CONCLUSIONS

This study revealed congruent phylogenetic relationships from comparative methods/datasets, with some conflicts being the result of incomplete lineage sorting and ancient/recent hybridization events. The strong association between the clades and the biogeographic areas suggested that ancient dispersal events shaped the modern distribution pattern, and regional climatic factors may have resulted in the following in situ diversification. Climate cooling starting during the mid-Miocene is associated with the global expansion of Zamia to tropical South America that has dramatically driven lineage diversification in the New World flora, as well as the extinction of cycad species in the nowadays cooler regions of both hemispheres, as indicated by the fossil records.

Evolutionary trends of reproductive phenotype in Cycadales: an analysis of morphological evolution in Ceratozamia

BACKGROUND AND AIMS

The size and shape of reproductive structures is especially relevant in evolution because these characters are directly related to the capacity for pollination and seed dispersal, a process that plays a basic role in evolutionary patterns. The evolutionary trajectories of reproductive phenotypes in gymnosperms have received special attention in terms of pollination and innovations related to the emergence of the spermatophytes. However, variability of reproductive structures, evolutionary trends and the role of environment in the evolution of cycad species have not been well documented and explored. This study considered this topic under an explicitly phylogenetic and evolutionary approach that included a broad sampling of reproductive structures in the genus Ceratozamia.

METHODS

We sampled 1400 individuals of 36 Ceratozamia species to explore the evolutionary pattern and identify and evaluate factors that potentially drove their evolution. We analysed characters for both pollen and ovulate strobili within a phylogenetic framework using different methods and characters (i.e. molecular and both quantitative and qualitative morphological) to infer phylogenetic relationships. Using this phylogenetic framework, evolutionary models of trait evolution for strobilar size were evaluated. In addition, quantitative morphological variation and its relation to environmental variables across species were analysed.

KEY RESULTS

We found contrasting phylogenetic signals between characters of pollen and ovulate strobili. These structures exhibited high morphological disparity in several characters related to size. Results of analyses of evolutionary trajectories suggested a stabilizing selection model. With regard to phenotype-environment, the analysis produced mixed results and differences for groups in the vegetation type where the species occur; however, a positive relationship with climatic variables was found.

CONCLUSIONS

The integrated approach synthesized reproductive phenotypic variation with current phylogenetic hypotheses and provided explicit statements of character evolution. The characters of volume for ovulate strobili were the most informative, and could provide a reference for further study of the evolutionary complexity in Ceratozamia. Finally, heterogeneous environments, which are under changing weather conditions, promote variability of reproductive structures.

The impact of thermogenesis on the origin of insect pollination

Thermogenesis in plants is the ability to raise their temperature above that of the surrounding air through metabolic processes, and is especially detected in reproductive organs. Warming benefits plants by facilitating the transmission of odours and compounds that attract insects. As a result, these plants increase their odds of being pollinated by the attracted insect. Modern thermogenesis has been reported in extant cycads and a small number of angiosperm lineages. Although thermogenesis is not directly preserved in the fossil record, it can be inferred by examining extant thermogenic plant lineages and comparing their features with those of the fossil record. We suggest that thermogenesis has probably occurred in seed plants for at least the past 200 million years, long before the origin of angiosperms. Thermogenesis in plants is an important factor that facilitated entomophilous pollination by enhancing the attraction of insects, complementary to other factors, thereby participating in the success of the two groups of organisms and providing many facets of past and recent reproductive biology for future exploration.

Fossils indicate marine dispersal in osteoglossid fishes, a classic example of continental vicariance

The separation of closely related terrestrial or freshwater species by vast marine barriers represents a biogeographical riddle. Such cases can provide evidence for vicariance, a process whereby ancient geological events like continental rifting divided ancestral geographical ranges. With an evolutionary history extending tens of millions of years, freshwater ecology, and distribution encompassing widely separated southern landmasses, osteoglossid bonytongue fishes are a textbook case of vicariance attributed to Mesozoic fragmentation of the Gondwanan supercontinent. Largely overlooked fossils complicate the clean narrative invoked for extant species by recording occurrences on additional continents and in marine settings. Here, we present a new total-evidence phylogenetic hypothesis for bonytongue fishes combined with quantitative models of range evolution and show that the last common ancestor of extant osteoglossids was likely marine, and that the group colonized freshwater settings at least four times when both extant and extinct lineages are considered. The correspondence between extant osteoglossid relationships and patterns of continental fragmentation therefore represents a striking example of biogeographical pseudocongruence. Contrary to arguments against vicariance hypotheses that rely only on temporal or phylogenetic evidence, these results provide direct palaeontological support for enhanced dispersal ability early in the history of a group with widely separated distributions in the modern day.

Limited effects of xylem anatomy on embolism resistance in cycad leaves

Drought-induced xylem embolism is a primary cause of plant mortality. Although c. 70% of cycads are threatened by extinction and extant cycads diversified during a period of increasing aridification, the vulnerability of cycads to embolism spread has been overlooked. We quantified the vulnerability to drought-induced embolism, pressure-volume curves, in situ water potentials, and a suite of xylem anatomical traits of leaf pinnae and rachises for 20 cycad species. We tested whether anatomical traits were linked to hydraulic safety in cycads. Compared with other major vascular plant clades, cycads exhibited similar embolism resistance to angiosperms and pteridophytes but were more vulnerable to embolism than noncycad gymnosperms. All 20 cycads had both tracheids and vessels, the proportions of which were unrelated to embolism resistance. Only vessel pit membrane fraction was positively correlated to embolism resistance, contrary to angiosperms. Water potential at turgor loss was significantly correlated to embolism resistance among cycads. Our results show that cycads exhibit low resistance to xylem embolism and that xylem anatomical traits - particularly vessels - may influence embolism resistance together with tracheids. This study highlights the importance of understanding the mechanisms of drought resistance in evolutionarily unique and threatened lineages like the cycads.

The meaning of mass extinctions and what the fossil record tells us about angiosperm survival at K-Pg: a reply to Hagen (2024)

Last year, we published research using phylogenetic comparative methods (PCMs) to reveal no phylogenetic evidence for elevated lineage-level extinction rates in angiosperms across K-Pg (Thompson JB, Ramírez-Barahona S. 2023 No phylogenetic evidence for angiosperm mass extinction at the Cretaceous-Palaeogene (K-Pg) boundary. Biol. Lett. 19, 20230314. (doi:10.1098/rsbl.2023.0314)), results that are in step with the global angiosperm fossil record. In a critique of our paper (Hagen ER. 2024 A critique of Thompson and Ramírez-Barahona (2023) or: how I learned to stop worrying and love the fossil record. EcoEvoRxiv. (doi:10.32942/X2631W)), simulation work is presented to argue we erred in our methodological choices and interpretations, and that we should have deferred to fossil evidence. In our opinion, underlying this critique are poor methodological choices on simulations and philosophical problems surrounding the definition of a mass extinction event, which leads to incorrect interpretations of both the fossil record and PCMs. We further argue that deferring to one source of evidence in favour of the other shuts the door to important evolutionary and philosophical questions.

Parallel evolution of angiosperm-like venation in Peltaspermales: a reinvestigation of Furcula

Leaf venation is a pivotal trait in the success of vascular plants. Whereas gymnosperms have single or sparsely branched parallel veins, angiosperms developed a hierarchical structure of veins that form a complex reticulum. Its physiological consequences are considered to have enabled angiosperms to dominate terrestrial ecosystems in the Late Cretaceous and Cenozoic. Although a hierarchical-reticulate venation also occurs in some groups of extinct seed plants, it is unclear whether these are stem relatives of angiosperms or have evolved these traits in parallel. Here, we re-examine the morphology of the enigmatic foliage taxon Furcula, a potential early Mesozoic angiosperm relative, and argue that its hierarchical vein network represents convergent evolution (in the Late Triassic) with flowering plants (which developed in the Early Cretaceous) based on details of vein architecture and the absence of angiosperm-like stomata and guard cells. We suggest that its nearest relatives are Peltaspermales similar to Scytophyllum and Vittaephyllum, the latter being a genus that originated during the Late Triassic (Carnian) and shares a hierarchical vein system with Furcula. We further suggest that the evolution of hierarchical venation systems in the early Permian, the Late Triassic, and the Early Cretaceous represent 'natural experiments' that might help resolve the selective pressures enabling this trait to evolve.

Eco-evolutionary evidence for the global diversity pattern of Cycas (Cycadaceae)

The evolution of the latitudinal diversity gradient (LDG), characterized by a peak in diversity toward the tropics, has captured significant attention in evolutionary biology and ecology. However, the inverse LDG (i-LDG) mechanism, wherein species richness increases toward the poles, remains inadequately explored. Cycads are among one of the oldest lineages of extant seed plants and have undergone extensive diversification in the tropics. Intriguingly, the extant cycad abundance exhibits an i-LDG pattern, and the underlying causes for this phenomenon remain largely elusive. Here, using 1,843 nuclear genes from a nearly complete sampling, we conducted comprehensive phylogenomic analyses to establish a robust species-level phylogeny for Cycas, the largest genus within cycads. We then reconstructed the spatial-temporal dynamics and integrated global environmental data to evaluate the roles of species ages, diversification rates, contemporary environment, and conservatism to ancestral niches in shaping the i-LDG pattern. We found Cycas experienced decreased diversification rates, coupled with the cooling temperature since its origin in the Eocene from continental Asia. Different regions have distinctively contributed to the formation of i-LDG for Cycas, with the northern hemisphere acting as evolutionary museums and the southern hemisphere serving as cradles. Moreover, water-related climate variables, specifically precipitation seasonality and potential evapotranspiration, were identified as paramount factors constraining Cycas species richness in the rainforest biome near the equator. Notably, the adherence to ancestral monsoonal climates emerges as a critical factor in sustaining the diversity pattern. This study underscores the imperative of integrating both evolutionary and ecological approaches to comprehensively unravel the mechanisms underpinning global biodiversity patterns.

Ceratonotha, a new erotylid genus (Coleoptera, Erotylidae) from late Eocene amber

Ceratonotha danica gen. nov., sp. nov., a new genus and species of Erotylidae (Coleoptera), is described and illustrated from Eocene Danish amber. Ceratonotha is morphologically similar to fossil erotylids from European amber, differing from them mainly in the length of the 4th tarsomere. Cycadophila mumia Alekseev & Bukejs, 2017, described from Baltic amber, belongs in this new genus and is here transferred, becoming Ceratonotha mumia (Alekseev & Bukejs) comb. nov..

Disparity of cycad leaves dispels the living fossil metaphor

The living fossil metaphor is tightly linked with the cycads. This group of gymnosperms is supposed to be characterised by long-term morphological stasis, particularly after their peak of diversity and disparity in the Jurassic. However, no formal test of this hypothesis exists. Here, we use a recent phylogenetic framework and an improved character matrix to reconstruct the Disparity Through Time for cycad leaves using a Principal Coordinate Analysis and employing Pre-Ordination Ancestral State Reconstruction to test the impact of sampling on the results. Our analysis shows that the cycad leaf morsphospace expanded up to the present, with numerous shifts in its general positioning, independently of sampling biases. Moreover, they also show that Zamiaceae expanded rapidly in the Early Cretaceous and continued to expand up to the present, while now-extinct clades experienced a slow contraction from their peak in the Triassic. We also show that rates of evolution were constantly high up to the Early Cretaceous, and then experienced a slight decrease in the Paleogene, followed by a Neogene acceleration. These results show a much more dynamic history for cycads, and suggest that the 'living fossil' metaphor is actually a hindrance to our understanding of their macroevolution.

Phylometagenomics of cycad coralloid roots reveals shared symbiotic signals

Cycads are known to host symbiotic cyanobacteria, including Nostocales species, as well as other sympatric bacterial taxa within their specialized coralloid roots. Yet, it is unknown if these bacteria share a phylogenetic origin and/or common genomic functions that allow them to engage in facultative symbiosis with cycad roots. To address this, we obtained metagenomic sequences from 39 coralloid roots sampled from diverse cycad species and origins in Australia and Mexico. Culture-independent shotgun metagenomic sequencing was used to validate sub-community co-cultures as an efficient approach for functional and taxonomic analysis. Our metanalysis shows a host-independent microbiome core consisting of seven bacterial orders with high species diversity within the identified taxa. Moreover, we recovered 43 cyanobacterial metagenome-assembled genomes, and in addition to Nostoc spp., symbiotic cyanobacteria of the genus Aulosira were identified for the first time. Using this robust dataset, we used phylometagenomic analysis to reveal three monophyletic cyanobiont clades, two host-generalist and one cycad-specific that includes Aulosira spp. Although the symbiotic clades have independently arisen, they are enriched in certain functional genes, such as those related to secondary metabolism. Furthermore, the taxonomic composition of associated sympatric bacterial taxa remained constant. Our research quadruples the number of cycad cyanobiont genomes and provides a robust framework to decipher cyanobacterial symbioses, with the potential of improving our understanding of symbiotic communities. This study lays a solid foundation to harness cyanobionts for agriculture and bioprospection, and assist in conservation of critically endangered cycads.

Embracing uncertainty: The way forward in plant fossil phylogenetics

Although molecular phylogenetics remains the most widely used method of inferring the evolutionary history of living groups, the last decade has seen a renewed interest in morphological phylogenetics, mostly driven by the promises that integrating the fossil record in phylogenetic trees offers to our understanding of macroevolutionary processes and dynamics and the possibility that the inclusion of fossil taxa could lead to more accurate phylogenetic hypotheses. The plant fossil record presents some challenges to its integration in a phylogenetic framework. Phylogenies including plant fossils often retrieve uncertain relationships with low support, or lack of resolution. This low support is due to the pervasiveness of morphological convergence among plant organs and the fragmentary nature of many plant fossils, and it is often perceived as a fundamental weakness reducing the utility of plant fossils in phylogenetics. Here I discuss the importance of uncertainty in morphological phylogenetics and how we can identify important information from different patterns and types of uncertainty. I also review a set of methodologies that can allow us to understand the causes underpinning uncertainty and how these practices can help us to further our knowledge of plant fossils. I also propose that a new visual language, including the use of networks instead of trees, represents an improvement on the old visualization based on consensus trees and more adequately serves phylogeneticists working with plant fossils. This set of methods and visualization tools represents an important way forward in a fundamental field for our understanding of the evolutionary history of plants.

Nitrogen isotopes reveal independent origins of N2-fixing symbiosis in extant cycad lineages

Cycads are ancient seed plants (gymnosperms) that emerged by the early Permian. Although they were common understory flora and food for dinosaurs in the Mesozoic, their abundance declined markedly in the Cenozoic. Extant cycads persist in restricted populations in tropical and subtropical habitats and, with their conserved morphology, are often called 'living fossils.' All surviving taxa receive nitrogen from symbiotic N2-fixing cyanobacteria living in modified roots, suggesting an ancestral origin of this symbiosis. However, such an ancient acquisition is discordant with the abundance of cycads in Mesozoic fossil assemblages, as modern N2-fixing symbioses typically occur only in nutrient-poor habitats where advantageous for survival. Here, we use foliar nitrogen isotope ratios-a proxy for N2 fixation in modern plants-to probe the antiquity of the cycad-cyanobacterial symbiosis. We find that fossilized cycad leaves from two Cenozoic representatives of extant genera have nitrogen isotopic compositions consistent with microbial N2 fixation. In contrast, all extinct cycad genera have nitrogen isotope ratios that are indistinguishable from co-existing non-cycad plants and generally inconsistent with microbial N2 fixation, pointing to nitrogen assimilation from soils and not through symbiosis. This pattern indicates that, rather than being ancestral within cycads, N2-fixing symbiosis arose independently in the lineages leading to living cycads during or after the Jurassic. The preferential survival of these lineages may therefore reflect the effects of competition with angiosperms and Cenozoic climatic change.

Geologically calibrated mammalian tree and its correlation with global events, including the emergence of humans

A robust timetree for Mammalia was constructed using the time calibration function of BEAST v1.10.4 and MEGA 11. The analysis involved the application of times of the most recent common ancestors, including a total of 19 mammalian fossil calibration ages following Benton et al. (Palaeontologia Electronica, 2015, 1-106) for their minimum ages. Additionally, fossil calibration ages for Gorilla, Pan, and a geologic event calibration age for otters were incorporated. Using these calibration ages, I constructed a geologically calibrated tree that estimates the age of the Homo and Pan splitting to be 5.69 Ma. The tree carries several significant implications. First, after the initial rifting at 120 Ma, the Atlantic Ocean expanded by over 500 km around Chron 34 (84 Ma), and vicariant speciation between Afrotheria (Africa) and Xenarthra (South America) appears to have commenced around 70 Ma. Additionally, ordinal level differentiations began immediately following the K-Pg boundary (66.0 Ma), supporting previous hypothesis that mammalian radiation rapidly filled ecological niches left vacant by non-avian dinosaurs. I constructed a diagram depicting the relationship between base substitution rate and age using an additional function in BEAST v1.10.4. The diagram reveals an exponential increase in the base substitution rate approaching recent times. This increased base substitution rate during the Neogene period may have contributed to the expansion of biodiversity, including the extensive adaptive radiation that led to the evolution of Homo sapiens. One significant driving factor behind this radiation could be attributed to the emergence and proliferation of C4 grasses since 20 Ma. These grasses have played a role in increasing carbon fixation, reducing atmospheric CO2 concentration, inducing global cooling, and initiating Quaternary glacial-interglacial cycles, thereby causing significant climatic changes.

Why plant diversity and distribution matter

This article is the Editorial for the Special Collection ‘Global plant diversity and distribution’. See https://www.newphytologist.org/global-plant-diversity for more details.