Recent Synchronous Radiation of a Living Fossil

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.

The evolution of ectomycorrhizal symbiosis and host-plant switches are the main drivers for diversification of Amanitaceae (Agaricales, Basidiomycota)

BACKGROUND

Evolutionary radiation is widely recognized as a mode of species diversification, but the drivers of the rapid diversification of fungi remain largely unknown. Here, we used Amanitaceae, one of the most diverse families of macro-fungi, to investigate the mechanism underlying its diversification.

RESULTS

The ancestral state of the nutritional modes was assessed based on phylogenies obtained from fragments of 36 single-copy genes and stable isotope analyses of carbon and nitrogen. Moreover, a number of time-, trait-, and paleotemperature-dependent models were employed to investigate if the acquisition of ectomycorrhizal (ECM) symbiosis and climate changes promoted the diversification of Amanitaceae. The results indicate that the evolution of ECM symbiosis has a single evolutionary origin in Amanitaceae. The earliest increase in diversification coincided with the acquisition of the ECM symbiosis with angiosperms in the middle Cretaceous. The recent explosive diversification was primarily triggered by the host-plant switches from angiosperms to the mixed forests dominated by Fagaceae, Salicaceae, and Pinaceae or to Pinaceae.

CONCLUSIONS

Our study provides a good example of integrating phylogeny, nutritional mode evolution, and ecological analyses for deciphering the mechanisms underlying fungal evolutionary diversification. This study also provides new insights into how the transition to ECM symbiosis has driven the diversification of fungi.

Phylogenomic analysis reveals the evolutionary history of Paleartic needle-leaved junipers

Needle-leaved junipers (Juniperus sect. Juniperus, Cupressaceae) are coniferous trees and shrubs with red or blue fleshy cones. They are distributed across Asia, Macaronesia and the Mediterranean Basin, with one species (J. communis) having a circumboreal distribution. Here we aim to resolve the phylogeny of this clade to infer its intricate evolutionary history. To do so, we built a comprehensive, time-calibrated phylogeny using genotyping-by-sequencing (GBS) and combine it with species occurrence using phylogeographic tools. Our results provide solid phylogenetic resolution to propose a new taxonomic classification and a biogeographical history of the section. Specifically, we confirm the monophyly of two groups within J. sect. Juniperus: the Asian (blue-cone) species including the circumboreal J. communis, and the Mediterranean-Macaronesian (red-cone) species. In addition, we provide strong phylogenetic evidence for three distinct species (J. badia, J. conferta, J. lutchuensis) previously considered subspecies or varieties, as well as for the differentiation between the eastern and western Mediterranean lineages of J. macrocarpa. Our findings suggest that the Mediterranean basin was the primary center of diversification for Juniperus sect. Juniperus, followed by an East Asian-Tethyan disjunction resulting from uplifts of the Qinghai-Tibetan Plateau and climatic shifts. The colonization history of Macaronesia by red-cone junipers from the western Mediterranean appears to have taken place independently in two different geological periods: the Miocene (Azores) and the Pliocene (Madeira-Canary Islands). Overall, genomic data and phylogenetic analysis are key to consider a new taxonomic proposal and reconstruct the biogeographical history of the iconic needle-leaved junipers across the Paleartic.

Genetically- and environmentally-dependent processes drive interspecific and intraspecific divergence in the Chinese relict endemic genus Dipteronia

China is a hotspot of relict plant species that were once widespread throughout the Northern Hemisphere. Recent research has demonstrated that the occurrence of long-term stable refugia in the mountainous regions of central and south-western China allowed their persistence through the late Neogene climate fluctuations. One of these relict lineages is Dipteronia, an oligotypic tree genus with a fossil record extending to the Paleocene. Here, we investigated the genetic variability, demographic dynamics and diversification patterns of the two currently recognized Dipteronia species (D ipteronia sinensis and D . dyeriana). Molecular data were obtained from 45 populations of Dipteronia by genotyping three cpDNA regions, two single copy nuclear genes and 15 simple sequence repeat loci. The genetic study was combined with niche comparison analyses on the environmental space, ecological niche modeling, and landscape connectivity analysis. We found that the two Dipteronia species have highly diverged both in genetic and ecological terms. Despite the incipient speciation processes that can be observed in D. sinensis, the occurrence of long-term stable refugia and, particularly, a dispersal corridor along Daba Shan-west Qinling, likely ensured its genetic and ecological integrity to date. Our study will not only help us to understand how populations of Dipteronia species responded to the tectonic and climatic changes of the Cenozoic, but also provide insight into how Arcto-Tertiary relict plants in East Asia survived, evolved, and diversified.

Phylotranscriptomic and ecological analyses reveal the evolution and morphological adaptation of Abies

Coniferous forests are under severe threat of the rapid anthropogenic climate warming. Abies (firs), the fourth-largest conifer genus, is a keystone component of the boreal and temperate dark-coniferous forests and harbors a remarkably large number of relict taxa. However, the uncertainty of the phylogenetic and biogeographic history of Abies significantly impedes our prediction of future dynamics and efficient conservation of firs. In this study, using 1,533 nuclear genes generated from transcriptome sequencing and a complete sampling of all widely recognized species, we have successfully reconstructed a robust phylogeny of global firs, in which four clades are strongly supported and all intersectional relationships are resolved, although phylogenetic discordance caused mainly by incomplete lineage sorting and hybridization was detected. Molecular dating and ancestral area reconstruction suggest a Northern Hemisphere high-latitude origin of Abies during the Late Cretaceous, but all extant firs diversified during the Miocene to the Pleistocene, and multiple continental and intercontinental dispersals took place in response to the late Neogene climate cooling and orogenic movements. Notably, four critically endangered firs endemic to subtropical mountains of China, including A. beshanzuensis, A. ziyuanensis, A. fanjingshanensis and A. yuanbaoshanensis from east to west, have different origins and evolutionary histories. Moreover, three hotspots of species richness, including western North America, central Japan, and the Hengduan Mountains, were identified in Abies. Elevation and precipitation, particularly precipitation of the coldest quarter, are the most significant environmental factors driving the global distribution pattern of fir species diversity. Some morphological traits are evolutionarily constrained, and those linked to elevational variation (e.g., purple cone) and cold resistance (e.g., pubescent branch and resinous bud) may have contributed to the diversification of global firs. Our study sheds new light on the spatiotemporal evolution of global firs, which will be of great help to forest management and species conservation in a warming world.

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.

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.

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.

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.

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.

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.

A novel seed cone of Pinus from the Miocene of coastal Southeast China indicates kinship with Southeast Asian pines

Pinus is an economically and ecologically important genus whose members are dominant components globally in low-latitude mountainous and mid-latitude temperate forests. Pinus species richness is currently concentrated in subtropical mid-low latitudes of the Northern Hemisphere, differing from the latitudinal diversity gradient mostly recognized in woody angiosperms. How the present pattern was developing in Earth's past is still poorly studied, particularly in eastern Asia. Here, a new fossil species, Pinus shengxianica sp. nov. is described based on a fossil seed cone from the Late Miocene Shengxian Formation in Zhejiang, southeast China. A co-occurring cone is recognized as a known fossil species, Pinus speciosa Li. Extensive comparison of extant and fossil members of Pinus suggests P. shengxianica shares a striking cone similarity to Pinus merkusii and Pinus latteri (subsection Pinus) from tropical Southeast Asia in having annular bulges around the umbo on the apophysis. The morphological resemblance indicates these two extant low-latitude pines probably possess a close affinity with the present newly-discovered P. shengxianica and originated from East Asian mid-low latitude ancestors during this generic re-diversification in the Miocene. This scenario is consistent with the evolutionary trajectory reflected by the pine fossil history and molecular data, marking the Miocene as a key period for the origin and evolution of most extant pines globally. The co-occurrences of diverse conifers and broadleaved angiosperms preferring diverse niches demonstrate Late Miocene eastern Zhejiang was one of the hot spots for coniferophyte diversity and hosted a needled-broadleaved mixed forest with complex vegetation structure and an altitudinal zonation.

Museomics unveil systematics, diversity and evolution of Australian cycad-pollinating weevils

Weevils have been shown to play significant roles in the obligate pollination of Australian cycads. In this study, we apply museomics to produce a first molecular phylogeny estimate of the Australian cycad weevils, allowing an assessment of their monophyly, placement and relationships. Divergence dating suggests that the Australian cycad weevils originated from the Late Oligocene to the Middle Miocene and that the main radiation of the cycad-pollinating groups occurred from the Middle to the Late Miocene, which is congruent with the diversification of the Australian cycads, thus refuting any notion of an ancient ciophilous system in Australia. Taxonomic studies reveal the existence of 19 Australian cycad weevil species and that their associations with their hosts are mostly non-species-specific. Co-speciation analysis shows no extensive co-speciation events having occurred in the ciophilous system of Australian cycads. The distribution pattern suggests that geographical factors, rather than diversifying coevolution, constitute the overriding process shaping the Australian cycad weevil diversity. The synchronous radiation of cycads and weevil pollinators is suggested to be a result of the post-Oligocene diversification common in Australian organisms.

Cretaceous pollen cone with three-dimensional preservation sheds light on the morphological evolution of cycads in deep time

The Cycadales are an ancient and charismatic group of seed plants. However, their morphological evolution in deep time is poorly understood. While molecular divergence time analyses estimate a Cretaceous origin for most major living cycad clades, much of the extant diversity is inferred to be a result of Neogene diversifications. This leads to long branches throughout the cycadalean phylogeny that, with few exceptions, have yet to be rectified by unequivocal fossil cycads. We report a permineralized pollen cone from the Campanian Holz Shale located in Silverado Canyon, CA, USA (c. 80 million yr ago). This fossil was studied via serial sectioning, SEM, 3D reconstruction and phylogenetic analyses. Microsporophyll and pollen morphology indicate this cone is assignable to Skyttegaardia, a recently described genus based on disarticulated lignitized microsporophylls from the Early Cretaceous of Denmark. Data from this new species, including a simple cone architecture, anatomical details and vasculature organization, indicate cycadalean affinities for Skyttegaardia. Phylogenetic analyses support this assignment and recover Skyttegaardia as crown-group Cycadales, nested within Zamiaceae. Our findings support a Cretaceous diversification for crown-group Zamiaceae, which included the evolution of morphological divergent extinct taxa with unique traits that have yet to be widely identified in the fossil record.

Time-calibrated phylogenies reveal mediterranean and pre-mediterranean origin of the thermophilous vegetation of the Canary Islands

BACKGROUND AND AIMS

The Canary Islands have strong floristic affinities with the Mediterranean Basin. One of the most characteristic and diverse vegetation belts of the archipelago is the thermophilous woodland (between 200 and 900 m.a.s.l.). This thermophilous plant community consists of many non-endemic species shared with the Mediterranean Floristic Region together with Canarian endemic species. Consequently, phytogeographic studies have historically proposed the hypothesis of an origin of the Canarian thermophilous species following the establishment of the summer-dry mediterranean climate in the Mediterranean Basin around 2.8 million years ago.

METHODS

Time-calibrated phylogenies for 39 plant groups including Canarian thermophilous species were primarily analysed to infer colonization times. In particular, we used 26 previously published phylogenies together with 13 new time-calibrated phylogenies (including newly generated plastid and nuclear DNA sequence data) to assess whether the time interval between stem and crown ages of Canarian thermophilous lineages postdates 2.8 Ma. For lineages postdating this time threshold, we additionally conducted ancestral area reconstructions to infer the potential source area for colonization.

KEY RESULTS

A total of 43 Canarian thermophilous lineages were identified from 39 plant groups. Both mediterranean (16) and pre-mediterranean (9) plant lineages were found. However, we failed to determine the temporal origin for 18 lineages because a stem-crown time interval overlaps with the 2.8-Ma threshold. The spatial origin of thermophilous lineages was also heterogeneous, including ancestral areas from the Mediterranean Basin (nine) and other regions (six).

CONCLUSIONS

Our findings reveal an unexpectedly heterogeneous origin of the Canarian thermophilous species in terms of colonization times and mainland source areas. A substantial proportion of the lineages arrived in the Canaries before the summer-dry climate was established in the Mediterranean Basin. The complex temporal and geographic origin of Canarian thermophilous species challenges the view of the Canary Islands (and Madeira) as a subregion within the Mediterranean Floristic Region.

Cycad phylogeny predicts host plant use of Eumaeus butterflies

Eumaeus butterflies are obligate herbivores of Zamia, the most diverse neotropical genus of cycads. Eumaeus-Zamia interactions have been characterized mainly for species distributed in North and Central America. However, larval host plant use by the southern Eumaeus clade remains largely unknown, precluding a comprehensive study of co-evolution between the genera. Here, we combine fieldwork with museum and literature surveys to expand herbivory records for Eumaeus from 21 to 38 Zamia species. We inferred a time-calibrated phylogeny of Eumaeus to test for distinct macroevolutionary scenarios of larval host plant conservatism and co-evolution. We found a remarkable coincidence between Eumaeus and Zamia diversification, with the butterfly stem group diverging at the same time as the most recent radiation of Zamia in the Miocene. Cophylogenetic reconciliation analyses show a strong cophylogenetic signal between cycads and their butterfly herbivores. Bipartite model-based approaches indicate that this is because closely related Zamia species are used by the same Eumaeus species, suggesting larval host plant resource tracking by the butterfly herbivores. Our results highlight a case of tight evolution between Eumaeus butterflies and cycads, pointing to the generality of correlated evolution and phylogenetic tracking in plant-herbivore interactions across seed plants.

Assessment of population, habitat, and threats to Cycas pectinata Buch.-Ham. (Cycadaceae), a vulnerable cycad in Bhutan

Cycas pectinata Buch.-Ham. is an ancient gymnosperm that is now confined to pockets of habitats in the Indian subcontinent but was once widely distributed. Its decline is attributable to habitat loss, and has reached the point where C. pectinata is listed as ‘Vulnerable’ in the IUCN Red List. C. pectinata is the only species of Cycas found in Bhutan, and in this biodiversity-rich area it is present as a relic of great scientific and conservation value because of its rarity and long evolutionary history. Although it is well known in India, it has not been studied in detail in Bhutan. This study assessed populations and threats to C. pectinata in two places in Bhutan. Field visits were made to document the distribution, habitats, and associated threats to the populations. Plants were observed growing in steep rugged terrain in the open Chir Pine forest. Populations are significantly threatened due to human activities such as habitat destruction and over collection as ornamental plant. Possible expansion of populations is naturally threatened by low seed production and by predators. 

Phytophthora pseudocryptogea, P. nicotianae and P. multivora Associated to Cycas revoluta: First Report Worldwide

A dieback was observed on three-year-old pot-grown plants of Cycas revoluta in Sicily (Italy). Symptoms, including stunting, yellowing and blight of the leaf crown, root rot and internal browning and decay of the basal stem, closely resembled the Phytophthora root and crown rot syndrome, common in other ornamentals. Isolations from rotten stem and roots, using a selective medium, and from rhizosphere soil of symptomatic plants, using leaf baiting, yielded three Phytophthora species, P. multivora, P. nicotianae and P. pseudocryptogea, were obtained. Isolates were identified by both morphological characters and DNA barcoding analysis, using three gene regions: ITS, β-tub and COI. Phytophthora pseudocryptogea was the sole species isolated directly from the stem and roots. The pathogenicity of the isolates of the three Phytophthora species was tested on one-year-old potted plants of C. revoluta, using both stem inoculation by wounding, and root inoculation through infested soil. Phytophthora pseudocryptogea was the most virulent and, like P. nicotianae, reproduced all the symptoms of natural infections, while P. multivora was the least virulent and induced solely very mild symptoms. Phytophthora pseudocryptogea was identified as the causal agent of the decline of C. revoluta, as it was re-isolated from both the roots and stems of artificially infected symptomatic plants, thus fulfilling Koch's postulates.

Phylotranscriptomics Shed Light on Intrageneric Relationships and Historical Biogeography of Ceratozamia (Cycadales)

Ceratozamia Brongn. is one of the species-rich genera of Cycadales comprising 38 species that are mainly distributed in Mexico, with a few species reported from neighboring regions. Phylogenetic relationships within the genus need detailed investigation based on extensive datasets and reliable systematic approaches. Therefore, we used 30 of the known 38 species to reconstruct the phylogeny based on transcriptome data of 3954 single-copy nuclear genes (SCGs) via coalescent and concatenated approaches and three comparative datasets (nt/nt12/aa). Based on all these methods, Ceratozamia is divided into six phylogenetic subclades within three major clades. There were a few discrepancies regarding phylogenetic position of some species within these subclades. Using these phylogenetic trees, biogeographic history and morphological diversity of the genus are explored. Ceratozamia originated from ancestors in southern Mexico since the mid-Miocene. There is a distinct distribution pattern of species through the Trans-Mexican Volcanic Belt (TMVB), that act as a barrier for the species dispersal at TMVB and its southern and northern part. Limited dispersal events occurred during the late Miocene, and maximum diversification happened during the Pliocene epoch. Our study provides a new insight into phylogenetic relationships, the origin and dispersal routes, and morphological diversity of the genus Ceratozamia. We also explain how past climatic changes affected the diversification of this Mesoamerica-native genus.

Arrested diversification? The phylogenetic distribution of poorly-diversifying lineages

Rapidly-diversifying lineages have been a major focus of modern evolutionary biology, with many hypotheses seeking to explain how they contribute to the uneven distribution of species in space and among taxa. However, an alternative view that is rarely explored is that some lineages evolve to become depauperate clades and show disproportionately low diversification, in a phenomenon we here call arrested diversification. In this study, we analyse several large-scale datasets including amphibian, squamate, mammal, and seed plant species to assess the extent to which poorly-diversifying lineages show distinct phylogenetic and spatial distributions in relation to other lineages. We found significant evidence that clades with low diversification rates tend to be more phylogenetically overdispersed than expected and show more idiosyncratic spatial distributions. These results suggest that arrested diversification is a real phenomenon that might play an important (yet largely overlooked) role in explaining asymmetries in the distribution of species across lineages.

Geographic factors and climatic fluctuation drive the genetic structure and demographic history of Cycas taiwaniana (Cycadaceae), an endemic endangered species to Hainan Island in China

Hainan Island had experienced several cold-warm and dry-humid fluctuations since the Late Pleistocene period, resulting in separating and connecting from the mainland several times with the cyclic rise and fall of sea level. The fluctuations can change the biota and ecological environment in the island. Cycas taiwaniana Carruthers is endemic to Hainan Island and is classified as endangered by the International Union for Conservation of Nature (IUCN). To comprehensively understand the genetic dynamics of C. taiwaniana, we sampled 12 wild populations in Hainan Island and one cultivated population in Fujian province, and analyzed the genetic diversity, genetic structure, and demographic history based on the molecular data. Results revealed that C. taiwaniana had relatively low genetic diversity and high genetic differentiation. Haplotypes of C. taiwaniana diversified during the Pleistocene based on the chloroplast DNA (cpDNA) and the concatenated nuclear DNA (nDNA) data. Genetic cluster analyses based on the microsatellite (SSR) data showed that the 12 wild populations were separated into three clusters which could be three evolutionary significant units (ESUs), indicating three basic units of protection were identified. Moreover, we also confirmed the cultivated population FJ derived from the DLS1-GSL clade. Demographic inference from different data was discordant, but overall, it uncovered that C. taiwaniana had experienced population contraction events twice during the Pleistocene and Holocene, and then expanded recently. Our study elucidated the population genetic characteristics of C. taiwaniana, and guided us to develop targeted conservation and management strategies for this endangered species.

Mega-sized pericentromeric blocks of simple telomeric repeats and their variants reveal patterns of chromosome evolution in ancient Cycadales genomes

Simple telomeric repeats composed of six to seven iterating nucleotide units are important sequences typically found at the ends of chromosomes. Here we analyzed their abundance and homogeneity in 42 gymnosperm (29 newly sequenced), 29 angiosperm (one newly sequenced), and eight bryophytes using bioinformatics, conventional cytogenetic and molecular biology approaches to explore their diversity across land plants. We found more than 10 000-fold variation in the amounts of telomeric repeats among the investigated taxa. Repeat abundance was positively correlated with increasing intragenomic sequence heterogeneity and occurrence at non-telomeric positions, but there was no correlation with genome size. The highest abundance/heterogeneity was found in the gymnosperm genus Cycas (Cycadaceae), in which megabase-sized blocks of telomeric repeats (i.e., billions of copies) were identified. Fluorescent in situ hybridization experiments using variant-specific probes revealed canonical Arabidopsis-type telomeric TTTAGGG repeats at chromosome ends, while pericentromeric blocks comprised at least four major telomeric variants with decreasing abundance: TTTAGGG>TTCAGGG >TTTAAGG>TTCAAGG. Such a diversity of repeats was not found in the sister cycad family Zamiaceae or in any other species analyzed. Using immunocytochemistry, we showed that the pericentromeric blocks of telomeric repeats overlapped with histone H3 serine 10 phosphorylation signals. We show that species of Cycas have amplified their telomeric repeats in centromeric and telomeric positions on telocentric chromosomes to extraordinary high levels. The ancestral chromosome number reconstruction suggests their occurrence is unlikely to be the product of ancient Robertsonian chromosome fusions. We speculate as to how the observed chromosome dynamics may be associated with the diversification of cycads.

Hidden species diversity in an iconic living fossil vertebrate

Ancient, species-poor lineages persistently occur across the Tree of life. These lineages are likely to contain unrecognized species diversity masked by the low rates of morphological evolution that characterize living fossils. Halecomorphi is a lineage of ray-finned fishes that diverged from its closest relatives before 200 Ma and is represented by only one living species in eastern North America, the bowfin, Amia calva Linnaeus. Here, we use double digest restriction-site-associated DNA sequencing and morphology to illuminate recent speciation in bowfins. Our results support the delimitation of a second living species of Amia, with the timing of diversification dating to the Plio-Pleistocene. This delimitation expands the species diversity of an ancient lineage that is integral to studies of vertebrate genomics and development, yet is facing growing conservation threats driven by the caviar fishery.

Monograph of Ceratozamia (Zamiaceae, Cycadales): an endangered genus

Ceratozamia (Zamiaceae, Cycadales), is a member one of the most endangered seed plant groups. Species of Ceratozamia grow from near sea level up to 2,100 meters in Mexico and Mesoamerica. We present a modern taxonomic treatment of Ceratozamia, based on fieldwork combined with detailed study of herbarium specimens in and from Mexico and Central America. This new revision is based on incorporation of morphological, molecular and biogeographic data that have been previously published in circumscriptions of species complexes by our group. Detailed morphological descriptions of the 36 species of Ceratozamia are provided and relevant characters for the genus are discussed and described. A majority are endemic to Mexico and are concentrated at high elevations in mountainous areas. Synonymies, lectotypifications, etymologies, taxonomic notes, distribution maps, illustrations and detailed species-level comparisons are included, as well as a dichotomous key for identification of all species. Data on distributional ranges and habitats of all species are summarized. Ceratozamiaosbornei D.W.Stev., Mart.-Domínguez & Nic.-Mor., sp. nov. is described from evergreen tropical forests of Belize and we highlight new populations and distributional ranges for C.subroseophylla Mart.-Domínguez & Nic.-Mor. and C.vovidesii Pérez-Farr. & Iglesias in the Mesoamerican region.

Correlations between leaf economics, mechanical resistance and drought tolerance across 41 cycad species

BACKGROUND AND AIMS

We conducted a comprehensive analysis of the functional traits of leaves (leaflets) of cycads. The aim of this study was to clarify the functional divergence between the earlier origin Cycadaceae and the later differentiated Zamiaceae, and the differences in trait associations between cycads and angiosperms.

METHODS

We selected 20 Cycadaceae species and 21 Zamiaceae species from the same cycad garden in South China, and measured their leaf structure, economic traits, mechanical resistance (Fp) and leaf water potential at the turgor loss point (πtlp). In addition, we compiled a dataset of geographical distribution along with climatic variables for these cycad species, and some leaf traits of tropical-sub-tropical angiosperm woody species from the literature for comparison.

KEY RESULTS

The results showed significantly contrasting leaf trait syndromes between the two families, with Zamiaceae species exhibiting thicker leaves, higher carbon investments and greater Fp than Cycadaceae species. Leaf thickness (LT) and πtlp were correlated with mean climatic variables in their native distribution ranges, indicating their evolutionary adaptation to environmental conditions. Compared with the leaves of angiosperms, the cycad leaves were thicker and tougher, and more tolerant to desiccation. Greater Fp was associated with a higher structural investment in both angiosperms and cycads; however, cycads showed lower Fp at a given leaf mass per area or LT than angiosperms. Enhancement of Fp led to more negative πtlp in angiosperms, but the opposite trend was observed in cycads.

CONCLUSIONS

Our results reveal that variations in leaf traits of cycads are mainly influenced by taxonomy and the environment of their native range. We also demonstrate similar leaf functional associations in terms of economics, but different relationships with regard to mechanics and drought tolerance between cycads and angiosperms. This study expands our understanding of the ecological strategies and likely responses of cycads to future climate change.

Early diversifications of angiosperms and their insect pollinators: were they unlinked?

The present-day ubiquity of angiosperm-insect pollination has led to the hypothesis that these two groups coevolved early in their evolutionary history. However, recent fossil discoveries and fossil-calibrated molecular dating analyses challenge the notion that early diversifications of angiosperms and insects were inextricably linked. In this article, we examine (i) the discrepancies between dates of emergence for angiosperms and major clades of insects; (ii) the long history of gymnosperm-insect pollination modes, which likely shaped early angiosperm-insect pollination mutualisms; and (iii) how the K-Pg (Cretaceous-Paleogene) mass extinction event was vital in propelling modern angiosperm-insect mutualisms. We posit that the early diversifications of angiosperms and their insect pollinators were largely decoupled until the end of the Cretaceous.

Evolutionary history of the Asplenium scolopendrium complex (Aspleniaceae), a relictual fern with a northern pan-temperate disjunct distribution

Asplenium scolopendrium is distributed in northern temperate forests with many global biogeographic disjunctions. The species complex of A. scolopendrium has been generated by spatial segregation coupled with divergent evolution. We elucidated the biogeographic history of the A. scolopendrium complex by exploring its origin, dispersal and evolution, thus providing insights into the evolutionary history of the Tertiary floras with northern pan-temperate disjunct distributions. The results revealed that all infraspecific taxa descended from a widely distributed common ancestor in the Northern Hemisphere. This pan-temperate ancestral population formed by unidirectional westward dispersal from European origins primarily during the Early Eocene when the Earth’s climate was much warmer than today. The splitting of European, American and East Asian lineages occurred during the Early Miocene due to geo-climatic vicariances. Polyploidy events in the American ancestral populations created additional reproductive barriers. The star-shaped haplotypes in each continent indicated that local disjunctions also led to derived genotypes with potential to diverge into different taxa. This intracontinental lineage splitting is likely related to latitudinal range shift and habitat fragmentation caused by glacial cycles and climate change during the Pleistocene. The evolutionary history of the A. scolopendrium complex supported the Boreotropical hypothesis exhibiting range expansion during the Early Eocene Climatic Optimum.

Are legumes different? Origins and consequences of evolving nitrogen fixing symbioses

Nitrogen fixing symbioses between plants and bacteria are ancient and, while not numerous, are formed in diverse lineages of plants ranging from microalgae to angiosperms. One symbiosis stands out as the most widespread one is that between legumes and rhizobia, leading to the formation of nitrogen-fixing nodules. The legume family is one of the largest and most diverse group of plants and legumes have been used by humans since the beginning of agriculture, both as high nitrogen food, as well as pastures and rotation crops. One open question is whether their ability to form a nitrogen-fixing symbiosis has contributed to legumes' success, and whether legumes have any unique characteristics that have made them more diverse and widespread than other groups of plants. This review examines the evolutionary journey that has led to the diversification of legumes, in particular its nitrogen-fixing symbiosis, and asks four questions to investigate which legume traits might have contributed to their success: 1. In what ways do legumes differ from other plant groups that have evolved nitrogen-fixing symbioses? In order to answer this question, the characteristics of the symbioses, and efficiencies of nitrogen fixation are compared between different groups of nitrogen fixing plants. 2. Could certain unique features of legumes be a reason for their success? This section examines the manifestations and possible benefits of a nitrogen-rich 'lifestyle' in legumes. 3. If nitrogen fixation was a reason for such a success, why have some species lost the symbiosis? Formation of symbioses has trade-offs, and while these are less well known for non-legumes, there are known energetic and ecological reasons for loss of symbiotic potential in legumes. 4. What can we learn from the unique traits of legumes for future crop improvements? While exploiting some of the physiological properties of legumes could be used to improve legume breeding, our increasing molecular understanding of the essential regulators of root nodule symbioses raise hope of creating new nitrogen fixing symbioses in other crop species.

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.

Taxonomic Revision of the Genus Miltotranes Zimmerman, 1994 (Coleoptera: Curculionidae: Molytinae), the Bowenia-Pollinating Cycad Weevils in Australia, with Description of a New Species and Implications for the Systematics of Bowenia

The Australian endemic weevils of the genus Miltotranes Zimmerman, 1994 (Curculionidae: Molytinae: Tranes group), comprising two species, M. prosternalis (Lea, 1929) and M. subopacus (Lea, 1929), are highly host-specific and the only known pollinators of Bowenia cycads, which comprise two CITES-protected species restricted to Tropical Queensland in Australia. In the present study, the taxonomy of Miltotranes is reviewed, a lectotype for the name Tranes prosternalis Lea, 1929 is designated and a new species associated with the Bowenia population in the McIlwraith Range is described as M. wilsoni sp. n. The descriptions and diagnoses of all species are supplemented with illustrations of their habitus and salient structures, and an identification key to all species and a distribution map are provided. Potential implications of the new species and of the taxonomy and biogeography of Miltotranes overall on the systematics and conservation of Bowenia are discussed.

Divergence With Gene Flow and Contrasting Population Size Blur the Species Boundary in Cycas Sect. Asiorientales, as Inferred From Morphology and RAD-Seq Data

The divergence process of incipient species is fascinating but elusive by incomplete lineage sorting or gene flow. Species delimitation is also challenging among those morphologically similar allopatric species, especially when lacking comprehensive data. Cycas sect. Asiorientales, comprised of C. taitungensis and C. revoluta in the Ryukyu Archipelago and Taiwan, diverged recently with continuous gene flow, resulting in a reciprocal paraphyletic relationship. Their previous evolutionary inferences are questioned from few genetic markers, incomplete sampling, and incomprehensive morphological comparison by a long-term taxonomic misconception. By whole range sampling, this study tests the geographic mode of speciation in the two species of Asiorientales by approximate Bayesian computation (ABC) using genome-wide single nucleotide polymorphisms (SNPs). The individual tree was reconstructed to delimit the species and track the gene-flow trajectory. With the comparison of diagnostic morphological traits and genetic data, the allopatric speciation was rejected. Alternatively, continuous but spatially heterogeneous gene flow driven by transoceanic vegetative dispersal and pollen flow with contrasting population sizes blurred their species boundary. On the basis of morphological, genetic, and evolutionary evidence, we synonymized these two Cycas species. This study highlights not only the importance of the Kuroshio Current to species evolution but also the disadvantage of using species with geographically structured genealogies as conservation units.

The Cycas genome and the early evolution of seed plants

Cycads represent one of the most ancient lineages of living seed plants. Identifying genomic features uniquely shared by cycads and other extant seed plants, but not non-seed-producing plants, may shed light on the origin of key innovations, as well as the early diversification of seed plants. Here, we report the 10.5-Gb reference genome of Cycas panzhihuaensis, complemented by the transcriptomes of 339 cycad species. Nuclear and plastid phylogenomic analyses strongly suggest that cycads and Ginkgo form a clade sister to all other living gymnosperms, in contrast to mitochondrial data, which place cycads alone in this position. We found evidence for an ancient whole-genome duplication in the common ancestor of extant gymnosperms. The Cycas genome contains four homologues of the fitD gene family that were likely acquired via horizontal gene transfer from fungi, and these genes confer herbivore resistance in cycads. The male-specific region of the Y chromosome of C. panzhihuaensis contains a MADS-box transcription factor expressed exclusively in male cones that is similar to a system reported in Ginkgo, suggesting that a sex determination mechanism controlled by MADS-box genes may have originated in the common ancestor of cycads and Ginkgo. The C. panzhihuaensis genome provides an important new resource of broad utility for biologists.

Towards the plastome evolution and phylogeny of Cycas L. (Cycadaceae): molecular-morphology discordance and gene tree space analysis

BACKGROUND

Plastid genomes (plastomes) present great potential in resolving multiscale phylogenetic relationship but few studies have focused on the influence of genetic characteristics of plastid genes, such as genetic variation and phylogenetic discordance, in resolving the phylogeny within a lineage. Here we examine plastome characteristics of Cycas L., the most diverse genus among extant cycads, and investigate the deep phylogenetic relationships within Cycas by sampling 47 plastomes representing all major clades from six sections.

RESULTS

All Cycas plastomes shared consistent gene content and structure with only one gene loss detected in Philippine species C. wadei. Three novel plastome regions (psbA-matK, trnN-ndhF, chlL-trnN) were identified as containing the highest nucleotide variability. Molecular evolutionary analysis showed most of the plastid protein-coding genes have been under purifying selection except ndhB. Phylogenomic analyses that alternatively included concatenated and coalescent methods, both identified four clades but with conflicting topologies at shallow nodes. Specifically, we found three species-rich Cycas sections, namely Stangerioides, Indosinenses and Cycas, were not or only weakly supported as monophyly based on plastomic phylogeny. Tree space analyses based on different tree-inference methods both revealed three gene clusters, of which the cluster with moderate genetic properties showed the best congruence with the favored phylogeny.

CONCLUSIONS

Our exploration in plastomic data for Cycas supports the idea that plastid protein-coding genes may exhibit discordance in phylogenetic signals. The incongruence between molecular phylogeny and morphological classification reported here may largely be attributed to the uniparental attribute of plastid, which cannot offer sufficient information to resolve the phylogeny. Contrasting to a previous consensus that genes with longer sequences and a higher proportion of variances are superior for phylogeny reconstruction, our result implies that the most effective phylogenetic signals could come from loci that own moderate variation, GC content, sequence length, and underwent modest selection.

Not that young: combining plastid phylogenomic, plate tectonic and fossil evidence indicates a Palaeogene diversification of Cycadaceae

BACKGROUND AND AIMS

Previous molecular dating studies revealed historical mass extinctions and recent radiations of extant cycads, but debates still exist between palaeobotanists and evolutionary biologists regarding the origin and evolution of Cycadaceae.

METHODS

Using whole plastomic data, we revisited the phylogeny of this family and found the Palawan endemic Cycas clade was strongly related to all lineages from Southeast Eurasia, coinciding with a plate drift event occurring in the Early Oligocene. By integrating fossil and biogeographical calibrations as well as molecular data from protein-coding genes, we established different calibration schemes and tested competing evolutionary timelines of Cycadaceae.

KEY RESULTS

We found recent dispersal cannot explain the distribution of Palawan Cycas, yet the scenario including the tectonic calibration yielded a mean crown age of extant Cycadaceae of ~69-43 million years ago by different tree priors, consistent with multiple Palaeogene fossils assigned to this family. Biogeographical analyses incorporating fossil distributions revealed East Asia as the ancestral area of Cycadaceae.

CONCLUSIONS

Our findings challenge the previously proposed Middle-Late Miocene diversification of cycads and an Indochina origin for Cycadaceae and highlight the importance of combining phylogenetic clades, tectonic events and fossils for rebuilding the evolutionary history of lineages that have undergone massive extinctions.

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.

Evolutionary legacies in contemporary tetrapod imperilment

The Tree of Life will be irrevocably reshaped as anthropogenic extinctions continue to unfold. Theory suggests that lineage evolutionary dynamics, such as age since origination, historical extinction filters and speciation rates, have influenced ancient extinction patterns - but whether these factors also contribute to modern extinction risk is largely unknown. We examine evolutionary legacies in contemporary extinction risk for over 4000 genera, representing ~30,000 species, from the major tetrapod groups: amphibians, birds, turtles and crocodiles, squamate reptiles and mammals. We find consistent support for the hypothesis that extinction risk is elevated in lineages with higher recent speciation rates. We subsequently test, and find modest support for, a primary mechanism driving this pattern: that rapidly diversifying clades predominantly comprise range-restricted, and extinction-prone, species. These evolutionary patterns in current imperilment may have important consequences for how we manage the erosion of biological diversity across the Tree of Life.

Congruence between ocean-dispersal modelling and phylogeography explains recent evolutionary history of Cycas species with buoyant seeds

Ocean currents play a significant role in driving the long-distance dispersal (LDD), spatial distribution and phylogeographic patterns of many organisms. Integrating phylogeographic analyses and mechanistic ocean current modelling can provide novel insights into the evolutionary history of terrestrial littoral species but has been rarely applied in this context. We focused on a group of Cycas that have buoyant seeds and occupy coastal habitats. By integrating evidence from mechanistic simulations and whole plastomic data, we examined the role of ocean circulation in shaping the phylogeography of these Cycas species. Plastomes of the studied Cycas species showed extreme conservatism, following a post-Pleistocene divergence. Phylogenies revealed three subclades, corresponding to the Pacific Ocean, Sunda Shelf and Indian Ocean. The ocean modelling results indicate that hotspots of seed stranding coincide well with the contemporary distribution of the Cycas species and that drifting trajectories from the three subclades are largely confined to separate regions. These findings suggest that ocean current systems, by driving long-distance dispersal, have shaped the distribution and phylogeography for Cycas with buoyant seeds. This study highlights how the combination of genomic data and ocean drift modelling can help explain phylogeographic patterns and diversity in terrestrial littoral ecosystems.

Spatiotemporal evolution of the global species diversity of Rhododendron

Evolutionary radiation is a widely recognized mode of species diversification, but its underlying mechanisms have not been unambiguously resolved for species-rich cosmopolitan plant genera. In particular, it remains largely unknown how biological and environmental factors have jointly driven its occurrence in specific regions. Here, we use Rhododendron, the largest genus of woody plants in the Northern Hemisphere, to investigate how geographic and climatic factors, as well as functional traits, worked together to trigger plant evolutionary radiations and shape the global patterns of species richness based on a solid species phylogeny. Using 3,437 orthologous nuclear genes, we reconstructed the first highly supported and dated backbone phylogeny of Rhododendron comprising 200 species that represent all subgenera, sections, and nearly all multispecies subsections, and found that most extant species originated by evolutionary radiations when the genus migrated southward from circumboreal areas to tropical/subtropical mountains, showing rapid increases of both net diversification rate and evolutionary rate of environmental factors in the Miocene. We also found that the geographically uneven diversification of Rhododendron led to a much higher diversity in Asia than in other continents, which was mainly driven by two environmental variables, that is, elevation range and annual precipitation, and were further strengthened by the adaptation of leaf functional traits. Our study provides a good example of integrating phylogenomic and ecological analyses in deciphering the mechanisms of plant evolutionary radiations, and sheds new light on how the intensification of the Asian monsoon has driven evolutionary radiations in large plant genera of the Himalaya-Hengduan Mountains.

Diversity and Genetic Structure of Dioon holmgrenii (Cycadales: Zamiaceae) in the Mexican Pacific Coast Biogeographic Province: Implications for Conservation

Dioon holmgrenii De Luca, Sabato et Vázq.Torres is an endangered species; it is endemic and its distribution is restricted to the biogeographic province of the Mexican Pacific Coast. The aim of this work was to determine the diversity and genetic structure of nine populations. The genetic diversity parameters and Wright’s F statistics were determined with six microsatellite loci. The genetic structure was determined by using the Structure software and by a discriminant analysis. The genetic diversity of the populations was high. The proportion of polymorphic loci was 0.89, the observed heterogeneity was higher (Ho = 0.62 to 0.98) than expected (He = 0.48 to 0.78), and the fixation index was negative (IF = −0.091 to −0.601). Heterozygous deficiency (FIT = 0.071) was found at the species level and heterozygotes excess (FIS = −0.287) at the population level. The genetic differentiation between populations was high (FST = 0.287), with the number of migrants less than one. Three groups of populations were differentiated, and the variation within populations, between populations, and between groups was: 65.5, 26.3, and 8.2%, respectively. Multiple factors explain the high genetic diversity, while the genetic structure is due to geographic barriers. Community reserves are urgent in at least one most diverse population of each group.

Inferring the Total-Evidence Timescale of Marattialean Fern Evolution in the Face of Model Sensitivity

Phylogenetic divergence-time estimation has been revolutionized by two recent developments: 1) total-evidence dating (or "tip-dating") approaches that allow for the incorporation of fossils as tips in the analysis, with their phylogenetic and temporal relationships to the extant taxa inferred from the data and 2) the fossilized birth-death (FBD) class of tree models that capture the processes that produce the tree (speciation, extinction, and fossilization) and thus provide a coherent and biologically interpretable tree prior. To explore the behavior of these methods, we apply them to marattialean ferns, a group that was dominant in Carboniferous landscapes prior to declining to its modest extant diversity of slightly over 100 species. We show that tree models have a dramatic influence on estimates of both divergence times and topological relationships. This influence is driven by the strong, counter-intuitive informativeness of the uniform tree prior, and the inherent nonidentifiability of divergence-time models. In contrast to the strong influence of the tree models, we find minor effects of differing the morphological transition model or the morphological clock model. We compare the performance of a large pool of candidate models using a combination of posterior-predictive simulation and Bayes factors. Notably, an FBD model with epoch-specific speciation and extinction rates was strongly favored by Bayes factors. Our best-fitting model infers stem and crown divergences for the Marattiales in the mid-Devonian and Late Cretaceous, respectively, with elevated speciation rates in the Mississippian and elevated extinction rates in the Cisuralian leading to a peak diversity of ${\sim}$2800 species at the end of the Carboniferous, representing the heyday of the Psaroniaceae. This peak is followed by the rapid decline and ultimate extinction of the Psaroniaceae, with their descendants, the Marattiaceae, persisting at approximately stable levels of diversity until the present. This general diversification pattern appears to be insensitive to potential biases in the fossil record; despite the preponderance of available fossils being from Pennsylvanian coal balls, incorporating fossilization-rate variation does not improve model fit. In addition, by incorporating temporal data directly within the model and allowing for the inference of the phylogenetic position of the fossils, our study makes the surprising inference that the clade of extant Marattiales is relatively young, younger than any of the fossils historically thought to be congeneric with extant species. This result is a dramatic demonstration of the dangers of node-based approaches to divergence-time estimation, where the assignment of fossils to particular clades is made a priori (earlier node-based studies that constrained the minimum ages of extant genera based on these fossils resulted in much older age estimates than in our study) and of the utility of explicit models of morphological evolution and lineage diversification. [Bayesian model comparison; Carboniferous; divergence-time estimation; fossil record; fossilized birth-death; lineage diversification; Marattiales; models of morphological evolution; Psaronius; RevBayes.].

Evolutionary patterns of nucleotide substitution rates in plastid genomes of Quercus

Molecular evolution, including nucleotide substitutions, plays an important role in understanding the dynamics and mechanisms of species evolution. Here, we sequenced whole plastid genomes (plastomes) of Quercus fabri, Quercus semecarpifolia, Quercus engleriana, and Quercus phellos and compared them with 14 other Quercus plastomes to explore their evolutionary relationships using 67 shared protein-coding sequences. While many previously identified evolutionary relationships were found, our findings do not support previous research which retrieve Quercus subg. Cerris sect. Ilex as a monophyletic group, with sect. Ilex found to be polyphyletic and composed of three strongly supported lineages inserted between sections Cerris and Cyclobalanposis. Compared with gymnosperms, Quercus plastomes showed higher evolutionary rates (Dn/Ds = 0.3793). Most protein-coding genes experienced relaxed purifying selection, and the high Dn value (0.1927) indicated that gene functions adjusted to environmental changes effectively. Our findings suggest that gene interval regions play an important role in Quercus evolution. We detected greater variation in the intergenic regions (trnH-psbA, trnK_UUU-rps16, trnfM_CAU-rps14, trnS_GCU-trnG_GCC, and atpF-atpH), intron losses (petB and petD), and pseudogene loss and degradation (ycf15). Additionally, the loss of some genes suggested the existence of gene exchanges between plastid and nuclear genomes, which affects the evolutionary rate of the former. However, the connective mechanism between these two genomes is still unclear.

Genome-Wide Evidence for Complex Hybridization and Demographic History in a Group of Cycas From China

Cycads represent one of the most ancestral living seed plants as well as one of the most threatened plant groups in the world. South China is a major center and potential origin of Cycas, the most rapidly diversified lineage of cycads. However, genomic-wide diversity of Cycas remains poorly understood due to the challenge of generating genomic markers associated with their inherent large genomes. Here, we perform a comprehensive conservation genomic study based on restriction-site associated DNA sequencing (RADseq) data in six representative species of Cycas in South China. Consistently low genetic diversity and strong genetic differentiation were detected across species. Both phylogenetic inference and genetic structure analysis via several methods revealed generally congruent groups among the six Cycas species. The analysis with ADMIXTURE showed low mixing of genetic composition among species, while individuals of C. dolichophylla exhibited substantial genetic admixture with C. bifida, C. changjiangensis, and C. balansae. Furthermore, the results from Treemix, f₄-statistic, and ABBA-BABA test were generally consistent and revealed the complex patterns of interspecific gene flow. Relatively strong signals of hybridization were detected between C. dolichophylla and C. szechuanensis, and the ancestor of C. taiwaniana and C. changjiangensis. Distinct patterns of demographic history were inferred for these species by Stairway Plot, and our results suggested that both climate fluctuation and frequent geological activities during the late Pleistocene exerted deep impacts on the population dynamics of these species in South China. Finally, we explore the practical implications of our findings for the development of conservation strategies in Cycas. The present study demonstrates the efficiency of RADseq for conservation genomic studies on non-model species with large and complex genomes. Given the great significance of cycads as a radical transition in the evolution of plant biodiversity, our study provides important insights into the mechanisms of diversification in such recently radiated living fossil taxa.

Gene duplications and phylogenomic conflict underlie major pulses of phenotypic evolution in gymnosperms

Inferring the intrinsic and extrinsic drivers of species diversification and phenotypic disparity across the tree of life is a major challenge in evolutionary biology. In green plants, polyploidy (or whole-genome duplication, WGD) is known to play a major role in microevolution and speciation, but the extent to which WGD has shaped macroevolutionary patterns of diversification and phenotypic innovation across plant phylogeny remains an open question. Here, we examine the relationship of various facets of genomic evolution-including gene and genome duplication, genome size, and chromosome number-with macroevolutionary patterns of phenotypic innovation, species diversification, and climatic occupancy in gymnosperms. We show that genomic changes, such as WGD and genome-size shifts, underlie the origins of most major extant gymnosperm clades, and notably, our results support an ancestral WGD in the gymnosperm lineage. Spikes of gene duplication typically coincide with major spikes of phenotypic innovation, while increased rates of phenotypic evolution are typically found at nodes with high gene-tree conflict, representing historic population-level dynamics during speciation. Most shifts in gymnosperm diversification since the rise of angiosperms are decoupled from putative WGDs and instead are associated with increased rates of climatic occupancy evolution, particularly in cooler and/or more arid climatic conditions, suggesting that ecological opportunity, especially in the later Cenozoic, and environmental heterogeneity have driven a resurgence of gymnosperm diversification. Our study provides critical insight on the processes underlying diversification and phenotypic evolution in gymnosperms, with important broader implications for the major drivers of both micro- and macroevolution in plants.

The complete mitochondrial genome of Cycas debaoensis revealed unexpected static evolution in gymnosperm species

Mitochondrial genomes of vascular plants are well known for their liability in architecture evolution. However, the evolutionary features of mitogenomes at intra-generic level are seldom studied in vascular plants, especially among gymnosperms. Here we present the complete mitogenome of Cycas debaoensis, an endemic cycad species to the Guangxi region in southern China. In addition to assemblage of draft mitochondrial genome, we test the conservation of gene content and mitogenomic stability by comparing it to the previously published mitogenome of Cycas taitungensis. Furthermore, we explored the factors such as structural rearrangements and nuclear surveillance of double-strand break repair (DSBR) proteins in Cycas in comparison to other vascular plant groups. The C. debaoensis mitogenome is 413,715 bp in size and encodes 69 unique genes, including 40 protein coding genes, 26 tRNAs, and 3 rRNA genes, similar to that of C. taitungensis. Cycas mitogenomes maintained the ancestral intron content of seed plants (26 introns), which is reduced in other lineages of gymnosperms, such as Ginkgo biloba, Taxus cuspidata and Welwitschia mirabilis due to selective pressure or retroprocessing events. C. debaoensis mitogenome holds 1,569 repeated sequences (> 50 bp), which partially account for fairly large intron size (1200 bp in average) of Cycas mitogenome. The comparison of RNA-editing sites revealed 267 shared non-silent editing site among predicted vs. empirically observed editing events. Another 33 silent editing sites from empirical data increase the total number of editing sites in Cycas debaoensis mitochondrial protein coding genes to 300. Our study revealed unexpected conserved evolution between the two Cycas species. Furthermore, we found strict collinearity of the gene order along with the identical set of genomic content in Cycas mt genomes. The stability of Cycas mt genomes is surprising despite the existence of large number of repeats. This structural stability may be related to the relative expansion of three DSBR protein families (i.e., RecA, OSB, and RecG) in Cycas nuclear genome, which inhibit the homologous recombinations, by monitoring the accuracy of mitochondrial chromosome repair.

Out of the Qinghai-Tibetan Plateau and rapid radiation across Eurasia for Allium section Daghestanica (Amaryllidaceae)

The disjunctive distribution (Europe-Caucasus-Asia) and species diversification across Eurasia for the genus Allium sect. Daghestanica has fascinating attractions for researchers aiming to understanding the development and history of modern Eurasia flora. However, no any studies have been carried out to address the evolutionary history of this section. Based on the nrITS and cpDNA fragments (trnL-trnF and rpl32-trnL), the evolutionary history of the third evolutionary line (EL3) of the genus Allium was reconstructed and we further elucidated the evolutionary line of sect. Daghestanica under this background. Our molecular phylogeny recovered two highly supported clades in sect. Daghestanica: the Clade I includes Caucasian-European species and Asian A. maowenense, A. xinlongense and A. carolinianum collected in Qinghai; the Clade II comprises Asian yellowish tepal species, A. chrysanthum, A. chrysocephalum, A. herderianum, A. rude and A. xichuanense. The divergence time estimation and biogeography inference indicated that Asian ancestor located in the Qinghai-Tibetan Plateau (QTP) and the adjacent region could have migrated to Caucasus and Europe distributions around the Late Miocene and resulted in further divergence and speciation; Asian ancestor underwent the rapid radiation in the QTP and the adjacent region most likely due to the heterogeneous ecology of the QTP resulted from the orogeneses around 4-3 million years ago (Mya). Our study provides a picture to understand the origin and species diversification across Eurasia for sect. Daghestanica.