Cenozoic extinctions account for the low diversity of extant gymnosperms compared with angiosperms

Phylogenetics of the rust fungi (Pucciniales) of South Africa, with notes on their life histories and possible origins

South Africa has an indigenous rust (Pucciniales) funga of approximately 460 species. This funga was sampled with species from as many genera as possible. The nuclear ribosomal large subunit (28S) region was amplified from samples representing 110 indigenous species, as well as the small subunit (18S) region and the cytochrome c oxidase subunit 3 (CO3) in some cases, and these were used in phylogenetic analyses. One new species is described, 12 new combinations made, six names reinstated, and two life history connections made. The life histories of this funga were summarized; it is dominated by species with contracted life histories. The majority of species are autoecious, with a small proportion being heteroecious. Of the autoecious species, many will likely be homothallic with no spermagonia. A shortened life history with homothallism allows for a single basidiospore infection to initiate a local population buildup under the prevailing unpredictable climatic conditions. Suggestions are made as to the possible origin of this funga based on the development of the modern South African flora. It is postulated that the rusts of South Africa are of relatively recent origin, consisting of three groups. Firstly, there is an African tropical element with members of the Mikronegerineae (Hemileia), the Sphaerophragmiaceae (Puccorchidium, Sphaerophragmium), and certain Uredinineae (Stomatisora). Their immediate ancestors likely occurred in the tropical forests of Africa during the Paleogene. Secondly, there is a pantropical element including the Raveneliaceae (e.g., Diorchidium, Maravalia, Ravenelia sensu lato, Uropyxis). This likely diversified during the Neogene, when the mimosoids became the dominant trees of the developing savannas. Thirdly, the Pucciniaceae invaded Africa as this continent pushed northward closing the Tethys Sea. They diversified with the development of the savannas as these become the dominant habitat in most of Africa, and are by far the largest component of the South African rust funga.

A revised phylogeny of Boletaceae using whole genome sequences

The porcini mushroom family Boletaceae is a diverse, widespread group of ectomycorrhizal (ECM) mushroom-forming fungi that so far has eluded intrafamilial phylogenetic resolution based on morphology and multilocus data sets. In this study, we present a genome-wide molecular data set of 1764 single-copy gene families from a global sampling of 418 Boletaceae specimens. The resulting phylogenetic analysis has strong statistical support for most branches of the tree, including the first statistically robust backbone. The enigmatic Phylloboletellus chloephorus from non-ECM Argentinian subtropical forests was recovered as a new subfamily sister to the core Boletaceae. Time-calibrated branch lengths estimate that the family first arose in the early to mid-Cretaceous and underwent a rapid radiation in the Eocene, possibly when the ECM nutritional mode arose with the emergence and diversification of ECM angiosperms. Biogeographic reconstructions reveal a complex history of vicariance and episodic long-distance dispersal correlated with historical geologic events, including Gondwanan origins and inferred vicariance associated with its disarticulation. Together, this study represents the most comprehensively sampled, data-rich molecular phylogeny of the Boletaceae to date, establishing a foundation for future robust inferences of biogeography in the group.

The angiosperm radiation played a dual role in the diversification of insects and insect pollinators

Interactions with angiosperms have been hypothesised to play a crucial role in driving diversification among insects, with a particular emphasis on pollinator insects. However, support for coevolutionary diversification in insect-plant interactions is weak. Macroevolutionary studies of insect and plant diversities support the hypothesis that angiosperms diversified after a peak in insect diversity in the Early Cretaceous. Here, we used the family-level fossil record of insects as a whole, and insect pollinator families in particular, to estimate diversification rates and the role of angiosperms on insect macroevolutionary history using a Bayesian process-based approach. We found that angiosperms played a dual role that changed through time, mitigating insect extinction in the Cretaceous and promoting insect origination in the Cenozoic, which is also recovered for insect pollinator families only. Although insects pollinated gymnosperms before the angiosperm radiation, a radiation of new pollinator lineages began as angiosperm lineages increased, particularly significant after 50 Ma. We also found that global temperature, increases in insect diversity, and spore plants were strongly correlated with origination and extinction rates, suggesting that multiple drivers influenced insect diversification and arguing for the investigation of different explanatory variables in further studies.

The phylogenetic reconstruction of the Neotropical cycad genus Ceratozamia (Zamiaceae) reveals disparate patterns of niche evolution

The cycad genus Ceratozamia comprises 40 species from Mexico, Guatemala, Belize, and Honduras, where cycads occur throughout climatically varied montane habitats. Ceratozamia has the potential to reveal the history and processes of species diversification across diverse Neotropical habitats in this region. However, the species relationships within Ceratozamia and the ecological trends during its evolution remain unclear. Here, we aimed to clarify the phylogenetic relationships, the timing of clade and species divergences, and the niche evolution throughout the phylogenetic history of Ceratozamia. Genome-wide DNA sequences were obtained with MIG-seq, and multiple data-filtering steps were used to optimize the dataset used to construct an ultrametric species tree. Divergence times among branches and ancestral niches were estimated. The niche variation among species was evaluated, summarized into two principal components, and their ancestral states were reconstructed to test whether niche shifts among branches can be explained by random processes, under a Brownian Motion model. Ceratozamia comprises three main clades, and most species relationships within the clades were resolved. Ceratozamia has diversified since the Oligocene, with major branching events occurring during the Miocene. This timing is consistent with fossil evidence, the timing estimated for other Neotropical plant groups, and the major geological events that shaped the topographic and climatic variation in Mexico. Patterns of niche evolution in the genus do not accord with the Brownian Motion model. Rather, non-random evolution with shifts towards more seasonal environments at high latitudes, or shifts towards humid or dry environments at low latitudes explain the diversification of Ceratozamia. We present a comprehensive phylogenetic reconstruction for Ceratozamia and identify for the first time the environmental factors involved in clade and species diversification within the genus. This study alleviates the controversies regarding the species relationships in the genus and provides the first evidence that latitude-associated environmental factors may influence processes of niche evolution in cycads.

Diploid species phylogeny and evolutionary reticulation indicate early radiation of Ephedra in the Tethys coast

Reconstructing a robust species phylogeny and disentangling the evolutionary and biogeographic history of the gymnosperm genus Ephedra, which has a large genome and rich polyploids, remain a big challenge. Here we reconstructed a transcriptome-based phylogeny of 19 diploid Ephedra species, and explored evolutionary reticulations in this genus represented by 50 diploid and polyploid species, using four low-copy nuclear and nine plastid genes. The diploid species phylogeny indicates that the Mediterranean species diverged first, and the remaining species split into three clades, including the American species (Clade A), E. rhytidosperma, and all other Asian species (Clade B). The single-gene trees placed E. rhytidosperma sister to Clade A, Clade B, or Clades A + B in similar proportions, suggesting that radiation and gene flow likely occurred in the early evolution of Ephedra. In addition, reticulate evolution occurred not only among the deep nodes, but also in the recently evolved South American species, which further caused difficulty in phylogenetic reconstruction. Moreover, we found that allopolyploid speciation was pervasive in Ephedra. Our study also suggests that Ephedra very likely originated in the Tethys coast during the late Cretaceous, and the South American Ephedra species have a single origin by dispersal from Mexico or North America.

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

The determinants of biodiversity patterns can be understood using macroevolutionary analyses. The integration of fossils into phylogenies offers a deeper understanding of processes underlying biodiversity patterns in deep time. Cycadales are considered a relict of a once more diverse and globally distributed group but are restricted to low latitudes today. We still know little about their origin and geographic range evolution. Combining molecular data for extant species and leaf morphological data for extant and fossil species, we study the origin of cycad global biodiversity patterns through Bayesian total-evidence dating analyses. We assess the ancestral geographic origin and trace the historical biogeography of cycads with a time-stratified process-based model. Cycads originated in the Carboniferous on the Laurasian landmass and expanded in Gondwana in the Jurassic. Through now-vanished continental connections, Antarctica and Greenland were crucial biogeographic crossroads for cycad biogeography. Vicariance is an essential speciation mode in the deep and recent past. Their latitudinal span increased in the Jurassic and restrained toward subtropical latitudes in the Neogene in line with biogeographic inferences of high-latitude extirpations. We show the benefits of integrating fossils into phylogenies to estimate ancestral areas of origin and to study evolutionary processes explaining the global distribution of present-day relict groups.

Phylotranscriptomics reveal the spatio-temporal distribution and morphological evolution of Macrozamia, an Australian endemic genus of Cycadales

BACKGROUND AND AIMS

Cycads are regarded as an ancient lineage of living seed plants, and hold important clues to understand the early evolutionary trends of seed plants. The molecular phylogeny and spatio-temporal diversification of one of the species-rich genera of cycads, Macrozamia, have not been well reconstructed.

METHODS

We analysed a transcriptome dataset of 4740 single-copy nuclear genes (SCGs) of 39 Macrozamia species and two outgroup taxa. Based on concatenated (maximum parsimony, maximum likelihood) and multispecies coalescent analyses, we first establish a well-resolved phylogenetic tree of Macrozamia. To identify cyto-nuclear incongruence, the plastid protein coding genes (PCGs) from transcriptome data are extracted using the software HybPiper. Furthermore, we explore the biogeographical history of the genus and shed light on the pattern of floristic exchange between three distinct areas of Australia. Six key diagnostic characters are traced on the phylogenetic framework using two comparative methods, and infra-generic classification is investigated.

KEY RESULTS

The tree topologies of concatenated and multi-species coalescent analyses of SCGs are mostly congruent with a few conflicting nodes, while those from plastid PCGs show poorly supported relationships. The genus contains three major clades that correspond to their distinct distributional areas in Australia. The crown group of Macrozamia is estimated to around 11.80 Ma, with a major expansion in the last 5-6 Myr. Six morphological characters show homoplasy, and the traditional phenetic sectional division of the genus is inconsistent with this current phylogeny.

CONCLUSIONS

This first detailed phylogenetic investigation of Macrozamia demonstrates promising prospects of SCGs in resolving phylogenetic relationships within cycads. Our study suggests that Macrozamia, once widely distributed in Australia, underwent major extinctions because of fluctuating climatic conditions such as cooling and mesic biome disappearance in the past. The current close placement of morphologically distinct species in the phylogenetic tree may be related to neotenic events that occurred in the genus.

Anatomical and biochemical evolutionary ancient traits of Araucaria araucana (Molina) K. Koch and their effects on carbon assimilation

The study of ancient species provides valuable information concerning the evolution of specific adaptations to past and current environmental conditions. Araucaria araucana (Molina) K. Koch belongs to one of the oldest families of conifers in the world, but despite this, there are few studies focused on its physiology and responses to changes in environmental conditions. We used an integrated approach aimed at comprehensively characterizing the ecophysiology of this poorly known species, focusing in its stomatal, mesophyll and biochemical traits, hypothesizing that these traits govern the carbon assimilation of A. araucana under past and present levels of atmospheric CO2. Results indicated that A. araucana presents the typical traits of an ancient species, such as large stomata and low stomatal density, which trigger low stomatal conductance and slow stomatal responsiveness to changing environmental conditions. Interestingly, the quantitative analysis showed that photosynthetic rates were equally limited by both diffusive and biochemical components. The Rubisco catalytic properties proved to have a low Rubisco affinity for CO2 and O2, similar to other ancient species. This affinity for CO2, together with the low carboxylation turnover rate, are responsible for the low Rubisco catalytic efficiency of carboxylation. These traits could be the result of the diverse environmental selective pressures that A. araucana was exposed during its diversification. The increase in measured temperatures induced an increase in stomatal and biochemical limitations, which together with a lower Rubisco affinity for CO2 could explain the low photosynthetic capacity of A. araucana in warmer conditions.

Deep evaluation of the evolutionary history of the Heat Shock Factor (HSF) gene family and its expansion pattern in seed plants

Heat shock factor (HSF) genes are essential in some of the basic developmental pathways in plants. Despite extensive studies on the structure, functional diversification, and evolution of HSF genes, their divergence history and gene duplication pattern remain unknown. To further illustrate the probable divergence patterns in these subfamilies, we analyzed the evolutionary history of HSF genes using phylogenetic reconstruction and genomic syntenic analyses, taking advantage of the increased sampling of genomic data from pteridophytes, gymnosperms and basal angiosperms. We identified a novel clade that includes HSFA2, HSFA6, HSFA7, and HSFA9 with a complex relationship, which is very likely due to orthologous or paralogous genes retained after frequent gene duplication events. We hypothesized that HSFA9 derives from HSFA2 through gene duplication in eudicots at the ancestral state, and then expanded in a lineage-specific way. Our findings indicate that HSFB3 and HSFB5 emerged before the divergence of ancestral angiosperms, but were lost in the most recent common ancestors of monocots. We also presumed that HSFC2 derives from HSFC1 in ancestral monocots. This work proposes that during the radiation of flowering plants, an era during which there was a differentiation of angiosperms, the size of the HSF gene family was also being adjusted with considerable sub- or neo-functionalization. The independent evolution of HSFs in eudicots and monocots, including lineage-specific gene duplication, gave rise to a new gene in ancestral eudicots and monocots, and lineage-specific gene loss in ancestral monocots. Our analyses provide essential insights for studying the evolutionary history of this multigene family.

Recent advances on phylogenomics of gymnosperms and a new classification

Living gymnosperms comprise four major groups: cycads, Ginkgo, conifers, and gnetophytes. Relationships among/within these lineages have not been fully resolved. Next generation sequencing has made available a large number of sequences, including both plastomes and single-copy nuclear genes, for reconstruction of solid phylogenetic trees. Recent advances in gymnosperm phylogenomic studies have updated our knowledge of gymnosperm systematics. Here, we review major advances of gymnosperm phylogeny over the past 10 years and propose an updated classification of extant gymnosperms. This new classification includes three classes (Cycadopsida, Ginkgoopsida, and Pinopsida), five subclasses (Cycadidae, Ginkgoidae, Cupressidae, Pinidae, and Gnetidae), eight orders (Cycadales, Ginkgoales, Araucariales, Cupressales, Pinales, Ephedrales, Gnetales, and Welwitschiales), 13 families, and 86 genera. We also described six new tribes including Acmopyleae Y. Yang, Austrocedreae Y. Yang, Chamaecyparideae Y. Yang, Microcachrydeae Y. Yang, Papuacedreae Y. Yang, and Prumnopityeae Y. Yang, and made 27 new combinations in the genus Sabina.

Limits to reproduction and seed size-number trade-offs that shape forest dominance and future recovery

The relationships that control seed production in trees are fundamental to understanding the evolution of forest species and their capacity to recover from increasing losses to drought, fire, and harvest. A synthesis of fecundity data from 714 species worldwide allowed us to examine hypotheses that are central to quantifying reproduction, a foundation for assessing fitness in forest trees. Four major findings emerged. First, seed production is not constrained by a strict trade-off between seed size and numbers. Instead, seed numbers vary over ten orders of magnitude, with species that invest in large seeds producing more seeds than expected from the 1:1 trade-off. Second, gymnosperms have lower seed production than angiosperms, potentially due to their extra investments in protective woody cones. Third, nutrient-demanding species, indicated by high foliar phosphorus concentrations, have low seed production. Finally, sensitivity of individual species to soil fertility varies widely, limiting the response of community seed production to fertility gradients. In combination, these findings can inform models of forest response that need to incorporate reproductive potential.

Dissecting the difference in tree species richness between Africa and South America

Our full-scale comparison of Africa and South America’s lowland tropical tree floras shows that both Africa and South America’s moist and dry tree floras are organized similarly: plant families that are rich in tree species on one continent are also rich in tree species on the other continent, and these patterns hold across moist and dry environments. Moreover, we confirm that there is an important difference in tree species richness between the two continents, which is linked to a few families that are exceptionally diverse in South American moist forests, although dry formations also contribute to this difference. Plant families only present on one of the two continents do not contribute substantially to differences in tree species richness.

Differences in species diversity over continental scales represent imprints of evolutionary, ecological, and biogeographic events. Here, we investigate whether the higher tree species richness in South America relative to Africa is due to higher richness in certain taxonomic clades, irrespective of vegetation type, or instead due to higher richness in specific biomes across all taxonomic clades. We used tree species inventory data to address this topic and began by clustering inventories from each continent based on species composition to derive comparable vegetation units. We found that moist forests in South America hold approximately four times more tree species than do moist forests in Africa, supporting previous studies. We also show that dry vegetation types in South America, such as tropical dry forests and savannas, hold twice as many tree species as do those in Africa, even though they cover a much larger area in Africa, at present and over geological time. Overall, we show that the marked species richness difference between South America and Africa is due primarily to a key group of families in the South American Amazon and Atlantic moist forests, which while present and speciose in Africa, are markedly less diverse there. Moreover, we demonstrate that both South American and African tree floras are organized similarly and that speciose families on one continent are likely speciose on the other. Future phylogenetic and functional trait work focusing on these key families should provide further insight into the processes leading to South America’s exceptional plant species diversity.

The Plastome Sequences of Triticum sphaerococcum (ABD) and Triticum turgidum subsp. durum (AB) Exhibit Evolutionary Changes, Structural Characterization, Comparative Analysis, Phylogenomics and Time Divergence

The mechanism and course of Triticum plastome evolution is currently unknown; thus, it remains unclear how Triticum plastomes evolved during recent polyploidization. Here, we report the complete plastomes of two polyploid wheat species, Triticum sphaerococcum (AABBDD) and Triticum turgidum subsp. durum (AABB), and compare them with 19 available and complete Triticum plastomes to create the first map of genomic structural variation. Both T. sphaerococcum and T. turgidum subsp. durum plastomes were found to have a quadripartite structure, with plastome lengths of 134,531 bp and 134,015 bp, respectively. Furthermore, diploid (AA), tetraploid (AB, AG) and hexaploid (ABD, AGA^m) Triticum species plastomes displayed a conserved gene content and commonly harbored an identical set of annotated unique genes. Overall, there was a positive correlation between the number of repeats and plastome size. In all plastomes, the number of tandem repeats was higher than the number of palindromic and forward repeats. We constructed a Triticum phylogeny based on the complete plastomes and 42 shared genes from 71 plastomes. We estimated the divergence of Hordeum vulgare from wheat around 11.04-11.9 million years ago (mya) using a well-resolved plastome tree. Similarly, Sitopsis species diverged 2.8-2.9 mya before Triticum urartu (AA) and Triticum monococcum (AA). Aegilops speltoides was shown to be the maternal donor of polyploid wheat genomes and diverged ~0.2-0.9 mya. The phylogeny and divergence time estimates presented here can act as a reference framework for future studies of Triticum evolution.

Prediction of Suitable Distribution of a Critically Endangered Plant Glyptostrobus pensilis

Glyptostrobus pensilis is a critically endangered living fossil plant species of the Mesozoic era, with high scientific research and economic value. The aim of this study was to assess the impact of climate change on the potential habitat area of G. pensilis in East Asia. The MaxEnt (maximum entropy) model optimized by the ENMeval data package was used to simulate the potential distribution habitats of G. pensilis since the last interglacial period (LIG, 120–140 ka). The results showed that the optimized MaxEnt model has a high prediction accuracy with the area under the receiver operating characteristic curve (AUC) of 0.9843 ± 0.005. The Current highly suitable habitats were found in the Northeast Jiangxi, Eastern Fujian and Eastern Guangdong; the main climatic factors affecting the geographic distribution of G. pensilis are temperature and precipitation, with precipitation as the temperature factor. The minimum temperature of coldest month (Bio6) may be the key factor restricting the northward distribution of G. pensilis; during the LIG, it contracted greatly in the highly suitable habitat area. Mean Diurnal Range (Bio2), Minimum Temperature of Coldest Month (Bio6), Annual Precipitation (Bio12) and Mean Temperature of Driest Quarter (Bio9) may be important climatic factors causing the changes in geographic distribution. In the next four periods, the suitable areas all migrated southward. Except for the RCP2.6-2070s, the highly suitable areas in the other three periods showed varying degrees of shrinkage. The results will provide a theoretical basis for the management and resource protection of G. pensilis.

The number of tree species on Earth

One of the most fundamental questions in ecology is how many species inhabit the Earth. However, due to massive logistical and financial challenges and taxonomic difficulties connected to the species concept definition, the global numbers of species, including those of important and well-studied life forms such as trees, still remain largely unknown. Here, based on global ground-sourced data, we estimate the total tree species richness at global, continental, and biome levels. Our results indicate that there are ∼73,000 tree species globally, among which ∼9,000 tree species are yet to be discovered. Roughly 40% of undiscovered tree species are in South America. Moreover, almost one-third of all tree species to be discovered may be rare, with very low populations and limited spatial distribution (likely in remote tropical lowlands and mountains). These findings highlight the vulnerability of global forest biodiversity to anthropogenic changes in land use and climate, which disproportionately threaten rare species and thus, global tree richness.

Insights from an ancient gymnosperm lineage: ambient temperature and light and the timing of thermogenesis in cycad cones

PREMISE

Although maintaining the appropriate mid-day timing of the diel thermogenic events of cones of the dioecious cycads Macrozamia lucida and M. macleayi is central to the survival of both plant and pollinator in this obligate pollination mutualism, the nature of the underlying mechanism remains obscure. We investigated whether it is under circadian control. Circadian mechanisms control the timing of many ecologically important processes in angiosperms, yet only a few gymnosperms have been studied in this regard.

METHODS

We subjected cones to different ambient temperature and lighting regimens (constant temperature and darkness; stepwise cool/warm ambient temperatures in constant darkness; stepwise dark/light exposures at constant temperature) to determine whether the resulting timing of their thermogenic events was consistent with circadian control.

RESULTS

Cones exposed to constant ambient temperature and darkness generated multiple temperature peaks endogenously, with an average interpeak-temperature period of 20.7 (±0.20) h that is temperature-compensated (Q₁₀  = 1.02). Exposure to 24-h ambient temperature cycles (12 h cool/12 h warm, constant darkness) yielded an interpeak-temperature period of 24.0 (±0.05) h, accurately and precisely replicating the ambient temperature period. Exposure to 24-h photo-cycles (12 h light/12 h dark, constant ambient temperature) yielded a shorter, more variable interpeak-temperature period of 23 (±0.23) h.

CONCLUSIONS

Our results indicate that cycad cone thermogenesis is under circadian clock control and differentially affected by ambient temperature and light cycles. Our data from cycads (an ancient gymnosperm lineage) adds to what little is known about circadian timing in gymnosperms, which have rarely been studied from the circadian perspective.

The dynamic history of gymnosperm plastomes: Insights from structural characterization, comparative analysis, phylogenomics, and time divergence

Gymnosperms are among the most endangered groups of plant species; they include ginkgo, pines (Conifers I), cupressophytes (Conifers II), cycads, and gnetophytes. The relationships among the five extant gymnosperm groups remain equivocal. We analyzed 167 available gymnosperm plastomes and investigated their diversity and phylogeny. We found that plastome size, structure, and gene order were highly variable in the five gymnosperm groups, of which Parasitaxus usta (Vieill.) de Laub. and Macrozamia mountperriensis F.M.Bailey had the smallest and largest plastomes, respectively. The inverted repeats (IRs) of the five groups were shown to have evolved through distinctive evolutionary scenarios. The IRs have been lost in all conifers but retained in cycads and gnetophytes. A positive association between simple sequence repeat (SSR) abundance and plastome size was observed, and the SSRs with the most variation were found in Pinaceae. Furthermore, the number of repeats was negatively correlated with IR length; thus, the highest number of repeats was detected in Conifers I and II, in which the IRs had been lost. We constructed a phylogeny based on 29 shared genes from 167 plastomes. With the plastome tree and 13 calibrations, we estimated the tree height between present-day angiosperms and gymnosperms to be ∼380 million years ago (mya). The placement of Gnetales in the tree agreed with the Gnetales-other gymnosperms hypothesis. The divergence between Ginkgo and cycads was estimated as ∼284 mya; the crown age of the cycads was 251 mya. Our time-calibrated plastid-based phylogenomic tree provides a framework for comparative studies of gymnosperm evolution.

Stomatal development in the cycad family Zamiaceae

BACKGROUND AND AIMS

The gymnosperm order Cycadales is pivotal to our understanding of seed-plant phylogeny because of its phylogenetic placement close to the root node of extant spermatophytes and its combination of both derived and plesiomorphic character states. Although widely considered a 'living fossil' group, extant cycads display a high degree of morphological and anatomical variation. We investigate stomatal development in Zamiaceae to evaluate variation within the order and homologies between cycads and other seed plants.

METHODS

Leaflets of seven species across five genera representing all major clades of Zamiaceae were examined at various stages of development using light microscopy and confocal microscopy.

KEY RESULTS

All genera examined have lateral subsidiary cells of perigenous origin that differ from other pavement cells in mature leaflets and could have a role in stomatal physiology. Early epidermal patterning in a 'quartet' arrangement occurs in Ceratozamia, Zamia and Stangeria. Distal encircling cells, which are sclerified at maturity, are present in all genera except Bowenia, which shows relatively rapid elongation and differentiation of the pavement cells during leaflet development.

CONCLUSIONS

Stomatal structure and development in Zamiaceae highlights some traits that are plesiomorphic in seed plants, including the presence of perigenous encircling subsidiary cells, and reveals a clear difference between the developmental trajectories of cycads and Bennettitales. Our study also shows an unexpected degree of variation among subclades in the family, potentially linked to differences in leaflet development and suggesting convergent evolution in cycads.

Non-adaptive evolutionary processes governed the diversification of a temperate conifer lineage after its migration into the tropics

Constructing phylogenetic relationships among closely related species is a recurrent challenge in evolutionary biology, particularly for long-lived taxa with large effective population sizes and uncomplete reproductive isolation, like conifers. Conifers further have slow evolutionary rates, which raises the question of whether adaptive or non/adaptive processes were predominantly involved when they rapidly diversified after migrating from temperate regions into the tropical mountains. Indeed, fine-scale phylogenetic relationships within several conifer genus remain under debate. Here, we studied the phylogenetic relationships of endemic firs (Abies, Pinaceae) discontinuously distributed in the montane forests from the Southwestern United States to Guatemala, and addressed several hypotheses related to adaptive and non-adaptive radiations. We derived over 80 K SNPs from genotyping by sequencing (GBS) for 45 individuals of nine Mesoamerican species to perform phylogenetic analyses. Both Maximum Likelihood and quartets-inference phylogenies resulted in a well-resolved topology, showing a single fir lineage divided in four subgroups that coincided with the main mountain ranges of Mesoamerica; thus having important taxonomic implications. Such subdivision fitted a North-South isolation by distance framework, in which non-adaptive allopatric processes seemed the rule. Interestingly, several reticulations were observed within subgroups, especially in the central-south region, which may explain past difficulties for generating infrageneric phylogenies. Further evidence for non-adaptive processes was obtained from analyses of 21 candidate-gene regions, which exhibited diminishing values of π_a/π_s and K_a/K_s with latitude, thus indicating reduced efficiency of purifying selection towards the Equator. Our study indicates that non-adaptive allopatric processes may be key generators of species diversity and endemism in the tropics.

Evolutionary analysis of chloroplast tRNA of Gymnosperm revealed the novel structural variation and evolutionary aspect

Gymnosperms such as ginkgo, conifers, cycads, and gnetophytes are vital components of land ecosystems, and they have significant economic and ecologic value, as well as important roles as forest vegetation. In this study, we investigated the structural variation and evolution of chloroplast transfer RNAs (tRNAs) in gymnosperms. Chloroplasts are important organelles in photosynthetic plants. tRNAs are key participants in translation where they act as adapter molecules between the information level of nucleic acids and functional level of proteins. The basic structures of gymnosperm chloroplast tRNAs were found to have family-specific conserved sequences. The tRNAΨ -loop was observed to contain a conforming sequence, i.e., U-U-C-N-A-N2. In gymnosperms, tRNAIle was found to encode a "CAU" anticodon, which is usually encoded by tRNAMet. Phylogenetic analysis suggested that plastid tRNAs have a common polyphyletic evolutionary pattern, i.e., rooted in abundant common ancestors. Analyses of duplication and loss events in chloroplast tRNAs showed that gymnosperm tRNAs have experienced little more gene loss than gene duplication. Transition and transversion analysis showed that the tRNAs are iso-acceptor specific and they have experienced unequal evolutionary rates. These results provide new insights into the structural variation and evolution of gymnosperm chloroplast tRNAs, which may improve our comprehensive understanding of the biological characteristics of the tRNA family.

The rise of angiosperms pushed conifers to decline during global cooling

Competition among species and entire clades can impact species diversification and extinction, which can shape macroevolutionary patterns. The fossil record shows successive biotic turnovers such that a dominant group is replaced by another. One striking example involves the decline of gymnosperms and the rapid diversification and ecological dominance of angiosperms in the Cretaceous. It is generally believed that angiosperms outcompeted gymnosperms, but the macroevolutionary processes and alternative drivers explaining this pattern remain elusive. Using extant time trees and vetted fossil occurrences for conifers, we tested the hypotheses that clade competition or climate change led to the decline of conifers at the expense of angiosperms. Here, we find that both fossil and molecular data show high congruence in revealing 1) low diversification rates, punctuated by speciation pulses, during warming events throughout the Phanerozoic and 2) that conifer extinction increased significantly in the Mid-Cretaceous (100 to 110 Ma) and remained high ever since. Their extinction rates are best explained by the rise of angiosperms, rejecting alternative models based on either climate change or time alone. Our results support the hypothesis of an active clade replacement, implying that direct competition with angiosperms increased the extinction of conifers by pushing their remaining species diversity and dominance out of the warm tropics. This study illustrates how entire branches on the Tree of Life may actively compete for ecological dominance under changing climates.

Niche conservatism promotes speciation in cycads: the case of Dioon merolae (Zamiaceae) in Mexico

Niche conservatism is the tendency of lineages to retain the same niche as their ancestors. It constrains biological groups and prevents ecological divergence. However, theory predicts that niche conservatism can hinder gene flow, strengthen drift and increase local adaptation: does it mean that it also can facilitate speciation? Why does this happen? We aim to answer these questions. We examined the variation of chloroplast DNA, genome-wide single nucleotide polymorphisms, morphological traits and environmental variables across the Dioon merolae cycad populations. We tested geographical structure, scenarios of demographic history, and niche conservatism between population groups. Lineage divergence is associated with the presence of a geographical barrier consisting of unsuitable habitats for cycads. There is a clear genetic and morphological distinction between the geographical groups, suggesting allopatric divergence. However, even in contrasting available environmental conditions, groups retain their ancestral niche, supporting niche conservatism. Niche conservatism is a process that can promote speciation. In D. merolae, lineage divergence occurred because unsuitable habitats represented a barrier against gene flow, incurring populations to experience isolated demographic histories and disparate environmental conditions. This study explains why cycads, despite their ancient lineage origin and biological stasis, have been able to diversify into modern ecosystems worldwide.

Eocene Araucaria Sect. Eutacta from Patagonia and floristic turnover during the initial isolation of South America

PREMISE

Eocene floras of Patagonia document biotic response to the final separation of Gondwana. The conifer genus Araucaria, distributed worldwide during the Mesozoic, has a disjunct extant distribution between South America and Australasia. Fossils assigned to Australasian Araucaria Sect. Eutacta usually are represented by isolated organs, making diagnosis difficult. Araucaria pichileufensis E.W. Berry, from the middle Eocene Río Pichileufú (RP) site in Argentine Patagonia, was originally placed in Sect. Eutacta and later reported from the early Eocene Laguna del Hunco (LH) locality. However, the relationship of A. pichileufensis to Sect. Eutacta and the conspecificity of the Araucaria material among these Patagonian floras have not been tested using modern methods.

METHODS

We review the type material of A. pichileufensis alongside large (n = 192) new fossil collections of Araucaria from LH and RP, including multi-organ preservation of leafy branches, ovuliferous complexes, and pollen cones. We use a total evidence phylogenetic analysis to analyze relationships of the fossils to Sect. Eutacta.

RESULTS

We describe Araucaria huncoensis sp. nov. from LH and improve the whole-plant concept for Araucaria pichileufensis from RP. The two species respectively resolve in the crown and stem of Sect. Eutacta.

CONCLUSIONS

Our results confirm the presence and indicate the survival of Sect. Eutacta in South America during early Antarctic separation. The exceptionally complete fossils significantly predate several molecular age estimates for crown Eutacta. The differentiation of two Araucaria species demonstrates conifer turnover during climate change and initial South American isolation from the early to middle Eocene.

Risk and Uncertainty of Losing Suitable Habitat Areas Under Climate Change Scenarios: A Case Study for 109 Gymnosperm Species in China

Taking 109 gymnosperm species in China as a case, the uncertainty and risk of losing habitat areas of gymnosperm species under future climate conditions were investigated via representative concentration pathways climate change scenarios, fuzzy set classifications and Monte Carlo techniques. Under nonrandom climate change scenarios, the richness of 109 species increased in the partial locations of northwestern and northeastern China and declined in the partial locations of eastern and central and southeastern China; the numbers of species that losing <20%, 20-40%, 40-60%, 60-80%, and over 80% of their current habitat areas were ~33-49, 36-40, 11-24, 7-9, and 2-8, respectively; ~99-105 species occupied over 80% of their total suitable areas and ~4-9 species occupied 60-80% their total suitable areas. Under random climate change scenarios, the number of species that losing various level of the habitat areas declined with enhancing probability; with a probabilities of over 0.6, the numbers of species that losing <20%, 20-40%, 40-60%, 60-80% and over 80% of their current habitat areas were ~19-28, 3-19, 0-3, 1-2, and 9-14, respectively, and the numbers of species that occupying ~20%, 20-40%, 40-60%, 60-80%, and over 80% of their total suitable areas were ~9-14, 4-11, 2-6, 1-3, and 34-45, respectively. Approximately 41% of 109 species will face extinction risks from climate change; the losing habitat areas in future climate condition will cause the varying of coniferous forest composition and the losing of ecosystem service related to the species; the uncertainty of losing distribution areas for species should not be ignored.

Early diversification and permeable species boundaries in the Mediterranean firs

BACKGROUND AND AIMS

Inferring the evolutionary relationships of species and their boundaries is critical in order to understand patterns of diversification and their historical drivers. Despite Abies (Pinaceae) being the second most diverse group of conifers, the evolutionary history of Circum-Mediterranean firs (CMFs) remains under debate.

METHODS

We used restriction site-associated DNA sequencing (RAD-seq) on all proposed CMF taxa to investigate their phylogenetic relationships and taxonomic status.

KEY RESULTS

Based on thousands of genome-wide single nucleotide polymorphisms (SNPs), we present here the first formal test of species delimitation, and the first fully resolved, complete species tree for CMFs. We discovered that all previously recognized taxa in the Mediterranean should be treated as independent species, with the exception of Abies tazaotana and Abies marocana. An unexpectedly early pulse of speciation in the Oligocene-Miocene boundary is here documented for the group, pre-dating previous hypotheses by millions of years, revealing a complex evolutionary history encompassing both ancient and recent gene flow between distant lineages.

CONCLUSIONS

Our phylogenomic results contribute to shed light on conifers' diversification. Our efforts to resolve the CMF phylogenetic relationships help refine their taxonomy and our knowledge of their evolution.

Araucarioid wood from the late Oligocene-early Miocene of Hainan Island: first fossil evidence for the genus Agathis in the Northern Hemisphere

Although many fossil and molecular data suggest migrations from Malesia and Asia to Australia appear to dominate floristic exchange between Australian and Asian rainforests, evidence is emerging that demonstrate dispersal of plant groups from Australia to Asia. In this paper, a new species Agathis ledongensis sp. nov. is described on the basis of silicified wood from the late Oligocene-early Miocene of the Qiutangling Formation in Ledong, Hainan Island, South China. It is the first fossil record of Agathis in the Northern Hemisphere, and the only known fossil evidence of its dispersal outside of Gondwana. The close affinity of the fossil wood from Ledong with the genus Agathis was confirmed by comparing quantitative traits in 31 wood samples of 20 species representing all three extant genera of the Araucariaceae. The percentage of tracheids with uniseriate pitting on radial walls is shown as an additional diagnostic trait for separating Agathis and Wollemia from Araucaria. The wood of Agathis ledongensis provides evidence for the dispersal of this important plant group from Australia, or another Gondwanan terrane, to eastern Asia based on reliable fossil data. It records the occurrence of this genus in Hainan Island by the early Miocene, i.e. at the beginning of the formation of the island chains between Australia and the South-East Asia and thus the provision of a land migration route. As the land routes between these continents were restricted at that time, the migration of Agathis to Malesia and Asia was presumably facilitated by long-distance dispersal of its winged seeds by wind.

Steep topography buffers threatened gymnosperm species against anthropogenic pressures in China

China is one of the most species-rich countries in the world, harboring many rare gymnosperms. Following recent human-led loss of forests, China is now experiencing increases in forest cover resulting from efforts of reforestation schemes. As anthropogenic activities have previously been found to interact with topography in shaping forest cover in China and considering the large human population and the ongoing population increase of the country, it is important to understand the role of anthropogenic pressures relative to environmental drivers for shaping species distributions here. Based on the well-established relationship between human population density and topography, we propose a hypothesis for explaining species distributions in a country dominated by human activities, predicting that species are more likely to occur in areas of steep topography under medium human population densities compared to low and high human population densities. Using species occurrence data from the Chinese Vascular Plant Distribution Database along with a common SDM method (maximum entropy modeling), we tested this hypothesis. Our results show that steep topography has the highest importance for predicting Chinese gymnosperm species occurrences in general, and threatened species specifically, in areas of medium human population densities. Consequently, these species are more often found in areas of steep terrain, supporting the proposed hypothesis. Results from this study highlight the need to include topographically heterogeneous habitats when planning new protected areas for species conservation.

Functional and morphological evolution in gymnosperms: A portrait of implicated gene families

Gymnosperms diverged from their sister plant clade of flowering plants 300 Mya. Morphological and functional divergence between the two major seed plant clades involved significant changes in their reproductive biology, water-conducting systems, secondary metabolism, stress defense mechanisms, and small RNA-mediated epigenetic silencing. The relatively recent sequencing of several gymnosperm genomes and the development of new genomic resources have enabled whole-genome comparisons within gymnosperms, and between angiosperms and gymnosperms. In this paper, we aim to understand how genes and gene families have contributed to the major functional and morphological differences in gymnosperms, and how this information can be used for applied breeding and biotechnology. In addition, we have analyzed the angiosperm versus gymnosperm evolution of the pleiotropic drug resistance (PDR) gene family with a wide range of functionalities in plants' interaction with their environment including defense mechanisms. Some of the genes reviewed here are newly studied members of gene families that hold potential for biotechnological applications related to commercial and pharmacological value. Some members of conifer gene families can also be exploited for their potential in phytoremediation applications.

Climate and evolutionary history define the phylogenetic diversity of vegetation types in the central region of South America

In South America the biogeographic history has produced different biomes with different vegetation types and distinct floras. As these vegetation types may diverge in evolutionary histories, we analysed how alpha and beta phylogenetic diversity vary across them and determine the main drivers of variation in phylogenetic diversity. To this end, we compiled a list of 205 sites and 1222 tree species spread over four biomes and eight vegetation types in central South America. For each site we evaluated six measures of evolutionary alpha diversity (species richness, phylogenetic diversity sensu stricto and the standardized effect size of phylogenetic diversity, mean phylogenetic distance and mean nearest taxon distance) and beta diversity (phylogenetic Sorensen's similarity). We checked the influence of spatial and environmental variables using generalized least squares models. The greatest phylogenetic differentiation was found between west and east of central South America, mainly between the Chaco communities and the other vegetation types, suggesting that species found in this biome come from different lineages, comparing with the others vegetation types. Our results also showed a clustered phylogenetic structure for the Dry Chaco woodlands, which may be associated with harsh environmental conditions. In addition to historical process, climatic conditions are the main drivers shaping phylogenetic patterns among the distinct vegetation types. Understanding patterns of phylogenetic diversity and distribution can greatly improve conservation planning and management since it allows the conservation of unique biome characteristics.

A General and Efficient Algorithm for the Likelihood of Diversification and Discrete-Trait Evolutionary Models

As the size of phylogenetic trees and comparative data continue to grow and more complex models are developed to investigate the processes that gave rise to them, macroevolutionary analyses are becoming increasingly limited by computational requirements. Here, we introduce a novel algorithm, based on the “flow” of the differential equations that describe likelihoods along tree edges in backward time, to reduce redundancy in calculations and efficiently compute the likelihood of various macroevolutionary models. Our algorithm applies to several diversification models, including birth–death models and models that account for state- or time-dependent rates, as well as many commonly used models of discrete-trait evolution, and provides an alternative way to describe macroevolutionary model likelihoods. As a demonstration of our algorithm’s utility, we implemented it for a popular class of state-dependent diversification models—BiSSE, MuSSE, and their extensions to hidden-states. Our implementation is available through the R package $\texttt{castor}$. We show that, for these models, our algorithm is one or more orders of magnitude faster than existing implementations when applied to large phylogenies. Our algorithm thus enables the fitting of state-dependent diversification models to modern massive phylogenies with millions of tips and may lead to potentially similar computational improvements for many other macroevolutionary models.

Resequencing 545 ginkgo genomes across the world reveals the evolutionary history of the living fossil

As Charles Darwin anticipated, living fossils provide excellent opportunities to study evolutionary questions related to extinction, competition, and adaptation. Ginkgo (Ginkgo biloba L.) is one of the oldest living plants and a fascinating example of how people have saved a species from extinction and assisted its resurgence. By resequencing 545 genomes of ginkgo trees sampled from 51 populations across the world, we identify three refugia in China and detect multiple cycles of population expansion and reduction along with glacial admixture between relict populations in the southwestern and southern refugia. We demonstrate multiple anthropogenic introductions of ginkgo from eastern China into different continents. Further analyses reveal bioclimatic variables that have affected the geographic distribution of ginkgo and the role of natural selection in ginkgo’s adaptation and resilience. These investigations provide insights into the evolutionary history of ginkgo trees and valuable genomic resources for further addressing various questions involving living fossil species.

Ginkgo is one of the living fossils from the plant kingdom. Here, authors conduct population genomics analyses to reveal its refugia and demographic history, and provide evidence of multiple anthropogenic introductions of ginkgo from eastern China into different continents.

Reconciling the concepts and measures of diversity, rarity and originality in ecology and evolution

The concept of biological diversity, or biodiversity, is at the core of evolutionary and ecological studies. Many indices of biodiversity have been developed in the last four decades, with species being one of the central units of these indices. However, evolutionary and ecological studies need a precise description of species' characteristics to best quantify inter-species diversity, as species are not equivalent and exchangeable. One of the first concepts characterizing species in biodiversity studies was abundance-based rarity. Abundance-based rarity was then complemented by trait- and phylo-based rarity, called species' trait-based and phylogenetic originalities, respectively. Originality, which is a property of an individual species, represents a species' contribution to the overall diversity of a reference set of species. Originality can also be defined as the rarity of a species' characteristics such as the state of a functional trait, which is often assumed to be represented by the position of the species on a phylogenetic tree. We review and compare various approaches for measuring originality, rarity and diversity and demonstrate that (i) even if attempts to bridge these concepts do exist, only a few ecological and evolutionary studies have tried to combine them all in the past two decades; (ii) phylo- and trait-based diversity indices can be written as a function of species rarity and originality measures in several ways; and (iii) there is a need for the joint use of these three types of indices to understand community assembly processes and species' roles in ecosystem functioning in order to protect biodiversity efficiently.

The early wasp plucks the flower: disparate extant diversity of sawfly superfamilies (Hymenoptera: ‘Symphyta’) may reflect asynchronous switching to angiosperm hosts

The insect order Hymenoptera originated during the Permian nearly 300 Mya. Ancestrally herbivorous hymenopteran lineages today make up the paraphyletic suborder ‘Symphyta’, which encompasses c. 8200 species with very diverse host-plant associations. We use phylogeny-based statistical analyses to explore the drivers of diversity dynamics within the ‘Symphyta’, with a particular focus on the hypothesis that diversification of herbivorous insects has been driven by the explosive radiation of angiosperms during and after the Cretaceous. Our ancestral-state estimates reveal that the first symphytans fed on gymnosperms, and that shifts onto angiosperms and pteridophytes – and back – have occurred at different time intervals in different groups. Trait-dependent analyses indicate that average net diversification rates do not differ between symphytan lineages feeding on angiosperms, gymnosperms or pteridophytes, but trait-independent models show that the highest diversification rates are found in a few angiosperm-feeding lineages that may have been favoured by the radiations of their host taxa during the Cenozoic. Intriguingly, lineages-through-time plots show signs of an early Cretaceous mass extinction, with a recovery starting first in angiosperm-associated clades. Hence, the oft-invoked assumption of herbivore diversification driven by the rise of flowering plants may overlook a Cretaceous global turnover in insect herbivore communities during the rapid displacement of gymnosperm- and pteridophyte-dominated floras by angiosperms.

A trait-based approach to the conservation of threatened plant species

Traditionally the vulnerability of threatened species to extinction has been assessed by studying their environment, genetics and population dynamics. A more comprehensive understanding of the factors promoting or limiting the long-term persistence of threatened species could be achieved by conducting an analysis of their functional responses to changing environments, their ecological interactions, and their role in ecosystem functioning. These less traditional research areas can be unified in a trait-based approach, a recent methodological advance in ecology that is being used to link individual-level functions to species, community and ecosystem processes to provide mechanistic explanations of observed patterns, particularly in changing environments. We illustrate how trait-based information can be translated into well-defined conservation strategies, using the example of Dioon sonorense, an Endangered cycad endemic to north-western Mexico. Scientific information yielded by trait-based research, coupled with existing knowledge derived from well-established traditional approaches, could facilitate the development of more integrative conservation strategies to promote the long-term persistence of individual threatened species. A comprehensive database of functional traits of threatened species would be of value in assisting the implementation of the trait-based approach.

The macroecology and macroevolution of plant species at risk

Contents Summary 708 I. Introduction 708 II. Fossils and phylogenies: learning from our past 708 III. Threatened species and Red Lists 710 IV. The geography of threat 710 V. The taxonomy of threat 710 VI. Predicting species at risk 711 VII. Conclusion 711 Acknowledgements 712 References 712 SUMMARY: Current rates of extinction are unprecedented in human history. The fossil record and newer molecular phylogenies illuminate historical patterns of speciation and extinction. They reveal both spectacular radiations and the characteristic features of mass extinction events in our geological past. The IUCN Red List provides insight into present-day species declines. There is emerging synthesis that species at risk are nonrandomly distributed across space and phylogeny. This pattern may be explained by geographical variation in driver intensity and species differential sensitivities. However, traits that confer resistance to one global change driver may increase susceptibility to a different driver. A complete understanding of extinction risk requires consideration of the interaction between extinction drivers, ecological traits, and species' evolutionary histories.

Turnover of southern cypresses in the post-Gondwanan world: extinction, transoceanic dispersal, adaptation and rediversification

Cupressaceae subfamily Callitroideae has been an important exemplar for vicariance biogeography, but its history is more than just disjunctions resulting from continental drift. We combine fossil and molecular data to better assess its extinction and, sometimes, rediversification after past global change. Key fossils were reassessed and their phylogenetic placement for calibration was determined using trait mapping and Bayes Factors. Five vicariance hypotheses were tested by comparing molecular divergence times with the timing of tectonic rifting. The role of adaptation to fire (serotiny) in its spread across a drying Australia was tested for Callitris. Our findings suggest that three transoceanic disjunctions within the Callitroideae probably arose from long-distance dispersal. A signature of extinction, centred on the end-Eocene global climatic chilling and drying, is evident in lineages-through-time plots and in the fossil record. Callitris, the most diverse extant callitroid genus, suffered extinctions but surviving lineages adapted and re-radiated into dry, fire-prone biomes that expanded in the Neogene. Serotiny, a key adaptation to fire, likely evolved in Callitris coincident with the biome shift. Both extinction and adaptive shifts have probably played major roles in this chronicle of turnover and renewal, but better understanding of biogeographical history requires improved taxonomy of fossils.

The dimensionality of niche space allows bounded and unbounded processes to jointly influence diversification

There are two prominent and competing hypotheses that disagree about the effect of competition on diversification processes. The first, the bounded hypothesis, suggests that species diversity is limited (bounded) by competition between species for finite ecological niche space. The second, the unbounded hypothesis, proposes that innovations associated with evolution render competition unimportant over macroevolutionary timescales. Here we use phylogenetically structured niche modelling to show that processes consistent with both of these diversification models drive species accumulation in conifers. In agreement with the bounded hypothesis, niche competition constrained diversification, and in line with the unbounded hypothesis, niche evolution and partitioning promoted diversification. We then analyse niche traits to show that these diversification enhancing and inhibiting processes can occur simultaneously on different niche dimensions. Together these results suggest a new hypothesis for lineage diversification based on the multi-dimensional nature of ecological niches that can accommodate both bounded and unbounded evolutionary processes.

Pleistocene diversification in an ancient lineage: a role for glacial cycles in the evolutionary history of Dioon Lindl. (Zamiaceae)

PREMISE OF THE STUDY

Recent estimates of crown ages for cycad genera (Late Miocene) challenge us to consider what processes have produced the extant diversity of this ancient group in such relatively little time. Pleistocene climate change has driven major shifts in species distributions in Mexico and may have led to speciation in the genus Dioon by forcing populations to migrate up in elevation, thereby becoming separated by topography.

METHODS

We inferred orthologs from transcriptomes of five species and sequenced these in 42 individuals representing all Dioon species. From these data and published plastid sequences, we inferred dated species trees and lineage-specific diversification rates.

KEY RESULTS

Analyses of 84 newly sequenced nuclear orthologs and published plastid data confirm four major clades within Dioon, all of Pleistocene age. Gene tree analysis, divergence dates, and an increase in diversification rate support very recent and rapid divergence of extant taxa.

CONCLUSIONS

This study confirms the Pleistocene age of Dioon species and implicates Pleistocene climate change and established topography in lineage spitting. These results add to our understanding of the cycads as evolutionarily dynamic lineages, not relicts or evolutionary dead ends. We also find that well-supported secondary calibration points can be reliable in the absence of fossils. Our hypothesis of lineage splitting mediated by habitat shifts may be applicable to other taxa that are restricted to elevation specific ecotones.

An overview of extant conifer evolution from the perspective of the fossil record

PREMISE OF THE STUDY

Conifers are an important living seed plant lineage with an extensive fossil record spanning more than 300 million years. The group therefore provides an excellent opportunity to explore congruence and conflict between dated molecular phylogenies and the fossil record.

METHODS

We surveyed the current state of knowledge in conifer phylogenetics to present a new time-calibrated molecular tree that samples ~90% of extant species diversity. We compared phylogenetic relationships and estimated divergence ages in this new phylogeny with the paleobotanical record, focusing on clades that are species-rich and well known from fossils.

KEY RESULTS

Molecular topologies and estimated divergence ages largely agree with the fossil record in Cupressaceae, conflict with it in Araucariaceae, and are ambiguous in Pinaceae and Podocarpaceae. Molecular phylogenies provide insights into some fundamental questions in conifer evolution, such as the origin of their seed cones, but using them to reconstruct the evolutionary history of specific traits can be challenging.

CONCLUSIONS

Molecular phylogenies are useful for answering deep questions in conifer evolution if they depend on understanding relationships among extant lineages. Because of extinction, however, molecular datasets poorly sample diversity from periods much earlier than the Late Cretaceous. This fundamentally limits their utility for understanding deep patterns of character evolution and resolving the overall pattern of conifer phylogeny.

Rethinking Living Fossils

Biologists would be mistaken if they relegated living fossils to paleontological inquiry or assumed that the concept is dead. It is now used to describe entities ranging from viruses to higher taxa, despite recent warnings of misleading inferences. Current work on character evolution illustrates how analyzing living fossils and stasis in terms of parts (characters) and wholes (e.g., organisms and lineages) advances our understanding of prolonged stasis at many hierarchical levels. Instead of viewing the concept's task as categorizing living fossils, we show how its primary role is to mark out what is in need of explanation, accounting for the persistence of both molecular and morphological traits. Rethinking different conceptions of living fossils as specific hypotheses reveals novel avenues for research that integrate phylogenetics, ecological and evolutionary modeling, and evo-devo to produce a more unified theoretical outlook.

Agathis trees of Patagonia's Cretaceous-Paleogene death landscapes and their evolutionary significance

PREMISE OF THE STUDY

The fossil record of Agathis historically has been restricted to Australasia. Recently described fossils from the Eocene of Patagonian Argentina showed a broader distribution than found previously, which is reinforced here with a new early Paleocene Agathis species from Patagonia. No previous phylogenetic analyses have included fossil Agathis species.

METHODS

We describe macrofossils from Patagonia of Agathis vegetative and reproductive organs from the early Danian, as well as leaves with Agathis affinities from the latest Maastrichtian. A total evidence phylogenetic analysis is performed, including the new Danian species together with other fossil species having agathioid affinities.

KEY RESULTS

Early Danian Agathis immortalis sp. nov. is the oldest definite occurrence of Agathis and one of the most complete Agathis species in the fossil record. Leafy twigs, leaves, pollen cones, pollen, ovuliferous complexes, and seeds show features that are extremely similar to the living genus. Dilwynites pollen grains, associated today with both Wollemia and Agathis and known since the Turonian, were found in situ within the pollen cones.

CONCLUSIONS

Agathis was present in Patagonia ca. 2 million years after the K-Pg boundary, and the putative latest Cretaceous fossils suggest that the genus survived the K-Pg extinction. Agathis immortalis sp nov. is recovered in a stem position for the genus, while A. zamunerae (Eocene, Patagonia) is recovered as part of the crown. A Mesozoic divergence for the Araucariaceae crown group, previously challenged by molecular divergence estimates, is supported by the combined phylogenetic analyses including the fossil taxa.

Araucaria lefipanensis (Araucariaceae), a new species with dimorphic leaves from the Late Cretaceous of Patagonia, Argentina

PREMISE OF THE STUDY

We describe a new araucarian species, Araucaria lefipanensis, from the Late Cretaceous flora of the Lefipán Formation, in Patagonia (Argentina) based on reproductive and vegetative remains, with a combination of characters that suggest mosaic evolution in the Araucaria lineage.

METHODS

The studied fossils were found at the Cañadón del Loro locality. Specimens were separated into two leaf morphotypes, and their morphological differences were tested with MANOVA.

KEY RESULTS

The new species Araucaria lefipanensis is erected based on the association of dimorphic leaves with cuticle remains and isolated cone scale complexes. The reproductive morphology is characteristic of the extant section Eutacta, whereas the vegetative organs resemble those of the sections Intermedia, Bunya, and Araucaria (the broad-leaved clade).

CONCLUSIONS

The leaf dimorphism of A. lefipanensis is similar to that of extant A. bidwillii, where dimorphism is considered to be related to seasonal growth. The leaf dimorphism in A. lefipanensis is consistent with the paleoclimatic and paleoenvironmental reconstructions previously suggested for the Lefipán Formation, which is thought to have been a seasonal subtropical forest. The new species shows evidence of mosaic evolution, with cone scale complexes morphologically similar to section Eutacta and leaves similar to the sections of the broad-leaved clade, constituting a possible transitional form between these two well-defined lineages. More complete plant concepts, especially those including both reproductive and vegetative remains are necessary to understand the evolution of ancient plant lineages. This work contributes to this aim by documenting a new species that may add to the understanding of the early evolution of the sections of Araucaria.