Origin and diversification of living cycads: a cautionary tale on the impact of the branching process prior in Bayesian molecular dating

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

This article comments on:

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

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

BACKGROUND AND AIMS

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

METHODS

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

KEY RESULTS

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

CONCLUSIONS

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

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

BACKGROUND AND AIMS

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

METHODS

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

KEY RESULTS

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

CONCLUSIONS

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

The non-flowering plants of a near-polar forest in East Gondwana, Tasmania, Australia, during the Early Eocene Climatic Optimum

PREMISE

The Cenozoic Macquarie Harbour Formation (MHF) hosts one of the oldest and southernmost post-Cretaceous fossil plant assemblages in Australia. Coinciding with the Early Eocene Climatic Optimum (EECO) and predating the breakup of Australia from Antarctica, it offers critical data to study the diversity and extent of the Austral Polar Forest Biome, and the floristic divergence between Australasia and South America resulting from the Gondwana breakup.

METHODS

The micromorphology and macromorphology of new fossil plant compressions from the MHF were described and systematically analyzed. Previously published non-flowering plant records were reviewed and revised. Macrofossil abundance data were provided. The flora was compared with other early Paleogene assemblages from across the Southern Hemisphere.

RESULTS

Twelve species of non-flowering plants were identified from the macrofossil record. Conifers include Araucariaceae (Araucaria macrophylla, A. readiae, A. timkarikensis sp. nov., and Araucarioides linearis), Podocarpaceae (Acmopyle glabra, Dacrycarpus mucronatus, Podocarpus paralungatikensis sp. nov., and Retrophyllum sp.), and Cupressaceae (Libocedrus microformis). Dacrycarpus linifolius was designated a junior synonym of D. mucronatus. Further components include a cycad (Bowenia johnsonii, Zamiaceae), a pteridosperm (Komlopteris cenozoicus, Umkomasiaceae), and a fern (Lygodium dinmorphyllum, Schizaeaceae).

CONCLUSIONS

The fossil assemblage represents a mixed near-polar forest with a high diversity of conifers. The morphology and preservation of several species indicate adaptations to life at high latitudes. The coexistence of large- and small-leaved conifers implies complex, possibly open forest structures. Comparisons with contemporaneous assemblages from Argentina support a circumpolar biome during the EECO, reaching from southern Australia across Antarctica to southern South America.

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.

Sequential trait evolution did not drive deep-time diversification in sharks

Estimating how traits evolved and impacted diversification across the tree of life represents a critical topic in ecology and evolution. Although there has been considerable research in comparative biology, large parts of the tree of life remain underexplored. Sharks are an iconic clade of marine vertebrates, and key components of marine ecosystems since the early Mesozoic. However, few studies have addressed how traits evolved or whether they impacted their extant diversity patterns. Our study aimed to fill this gap by reconstructing the largest time-calibrated species-level phylogeny of sharks and compiling an exhaustive database for ecological (diet, habitat) and biological (reproduction, maximum body length) traits. Using state-of-the-art models of evolution and diversification, we outlined the major character shifts and modes of trait evolution across shark species. We found support for sequential models of trait evolution and estimated a small to medium-sized lecithotrophic and coastal-dwelling most recent common ancestor for extant sharks. However, our exhaustive hidden traits analyses do not support trait-dependent diversification for any examined traits, challenging previous works. This suggests that the role of traits in shaping sharks' diversification dynamics might have been previously overestimated and should motivate future macroevolutionary studies to investigate other drivers of diversification in this clade.

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 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.

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.

Gut microbiomes of cycad-feeding insects tolerant to β-methylamino-L-alanine (BMAA) are rich in siderophore biosynthesis

Ingestion of the cycad toxins β-methylamino-L-alanine (BMAA) and azoxyglycosides is harmful to diverse organisms. However, some insects are specialized to feed on toxin-rich cycads with apparent immunity. Some cycad-feeding insects possess a common set of gut bacteria, which might play a role in detoxifying cycad toxins. Here, we investigated the composition of gut microbiota from a worldwide sample of cycadivorous insects and characterized the biosynthetic potential of selected bacteria. Cycadivorous insects shared a core gut microbiome consisting of six bacterial taxa, mainly belonging to the Proteobacteria, which we were able to isolate. To further investigate selected taxa from diverging lineages, we performed shotgun metagenomic sequencing of co-cultured bacterial sub-communities. We characterized the biosynthetic potential of four bacteria from Serratia, Pantoea, and two different Stenotrophomonas lineages, and discovered a suite of biosynthetic gene clusters notably rich in siderophores. Siderophore semi-untargeted metabolomics revealed a broad range of chemically related yet diverse iron-chelating metabolites, including desferrioxamine B, suggesting the occurrence of an unprecedented desferrioxamine-like biosynthetic pathway that remains to be identified. These results provide a foundation for future investigations into how cycadivorous insects tolerate diets rich in azoxyglycosides, BMAA, and other cycad toxins, including a possible role for bacterial siderophores.

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.

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.

DNA barcodes provide insights into the diversity and biogeography of the non-biting midge Polypedilum (Diptera, Chironomidae) in South America

South America, particularly within its tropical belt, is renowned for its unparalleled high levels of species richness, surpassing other major biomes. Certain neotropical areas harbor fragmented knowledge of insect diversity and face imminent threats from biodiversity loss and climate change. Hence, there is an urgent need for rapid estimation methods to complement slower traditional taxonomic approaches. A variety of algorithms for delimiting species through single-locus DNA barcodes have been developed and applied for rapid species diversity estimates across diverse taxa. However, tree-based and distance-based methods may yield different group assignments, leading to potential overestimation or underestimation of putative species. Here, we investigate the performance of different DNA-based species delimitation approaches to rapidly estimate the diversity of Polypedilum (Chironomidae, Diptera) in South America. Additionally, we test the hypothesis that significant differences exist in the community structure of Polypedilum fauna between South America and its neighboring regions, particularly the Nearctic. Our analysis encompasses a dataset of 1492 specimens from 598 locations worldwide, with a specific focus on South America. Within this region, we analyzed a subset of 247 specimens reported from 37 locations. Using various methods including the Barcode Index Number (BIN), Bayesian Poisson tree processes (bPTP), multi-rate Poisson tree processes (mPTP), single-rate Poisson tree processes (sPTP), and generalized mixed Yule coalescent (sGMYC), we identify molecular operational taxonomic units (MOTUs) ranging from 267 to 520. Our results indicate that the sGMYC method is the most suitable for estimating putative species in our dataset, resulting in the identification of 75 species in the Neotropical region, particularly in South America. Notably, this region exhibited higher species richness in comparison to the Palearctic and Oriental realms. Additionally, our findings suggest potential differences in species composition of Polypedilum fauna between the Neotropical and the adjacent Nearctic realms, highlighting high levels of endemism and species richness in the first. These results support our hypothesis that there are substantial differences exist in species composition between the Polypedilum fauna in South America and the neighboring regions.

Increasing Information Content and Diagnosability in Family-Level Classifications

Higher-level classifications often must account for monotypic taxa representing depauperate evolutionary lineages and lacking synapomorphies of their better-known, well-defined sister clades. In a ranked (Linnean) or unranked (phylogenetic) classification system, discovering such a depauperate taxon does not necessarily invalidate the rank classification of sister clades. Named higher taxa must be monophyletic to be phylogenetically valid. Ranked taxa above the species level should also maximize information content, diagnosability, and utility (e.g., in biodiversity conservation). In spider classification, families are the highest rank that is systematically catalogued, and incertae sedis is not allowed. Consequently, it is important that family-level taxa be well defined and informative. We revisit the classification problem of Orbipurae, an unranked suprafamilial clade containing the spider families Nephilidae, Phonognathidae, and Araneidae sensu stricto. We argue that, to maximize diagnosability, information content, conservation utility, and practical taxonomic considerations, this "splitting" scheme is superior to its recently proposed alternative, which lumps these families together as Araneidae sensu lato. We propose to redefine Araneidae and recognize a monogeneric spider family, Paraplectanoididae fam. nov. to accommodate the depauperate lineage Paraplectanoides. We present new subgenomic data to stabilize Orbipurae topology which also supports our proposed family-level classification. Our example from spiders demonstrates why classifications must be able to accommodate depauperate evolutionary lineages, for example, Paraplectanoides. Finally, although clade age should not be a criterion to determine rank, other things being equal, comparable ages of similarly ranked taxa do benefit comparative biology. [Classification, family rank, phylogenomics, systematics, monophyly, spider phylogeny.].

Ghost introgression facilitates genomic divergence of a sympatric cryptic lineage in Cycas revoluta

A cryptic lineage is a genetically diverged but morphologically unrecognized variant of a known species. Clarifying cryptic lineage evolution is essential for quantifying species diversity. In sympatric cryptic lineage divergence compared with allopatric divergence, the forces of divergent selection and mating patterns override geographical isolation. Introgression, by supplying preadapted or neutral standing genetic variations, can promote sympatric cryptic lineage divergence via selection. However, most studies concentrated on extant species introgression, ignoring the genetic legacy of introgression from extinct or unsampled lineages ("ghost introgression"). Cycads are an ideal plant for studying the influence of ghost introgression because of their common interspecific gene flow and past high extinction rate. Here, we utilized reference-based ddRADseq to clarify the role of ghost introgression in the evolution of a previously identified sympatric cryptic lineage in Cycas revoluta. After re-evaluating the evolutionary independency of cryptic lineages, the group-wise diverged single-nucleotide polymorphisms among sympatric and allopatric lineages were compared and functionally annotated. Next, we employed an approximate Bayesian computation method for hypothesis testing to clarify the cryptic lineage evolution and ghost introgression effect. SNPs with the genomic signatures of ghost introgression were further annotated. Our results reconfirmed the evolutionary independency of cryptic lineage among C. revoluta and demonstrated that ghost introgression to the noncryptic lineage facilitated their divergence. Gene function related to heat stress and disease resistance implied ecological adaptation of the main extant populations of C. revoluta.

Diversification Models Conflate Likelihood and Prior, and Cannot be Compared Using Conventional Model-Comparison Tools

Time-calibrated phylogenetic trees are a tremendously powerful tool for studying evolutionary, ecological, and epidemiological phenomena. Such trees are predominantly inferred in a Bayesian framework, with the phylogeny itself treated as a parameter with a prior distribution (a "tree prior"). However, we show that the tree "parameter" consists, in part, of data, in the form of taxon samples. Treating the tree as a parameter fails to account for these data and compromises our ability to compare among models using standard techniques (e.g., marginal likelihoods estimated using path-sampling and stepping-stone sampling algorithms). Since accuracy of the inferred phylogeny strongly depends on how well the tree prior approximates the true diversification process that gave rise to the tree, the inability to accurately compare competing tree priors has broad implications for applications based on time-calibrated trees. We outline potential remedies to this problem, and provide guidance for researchers interested in assessing the fit of tree models. [Bayes factors; Bayesian model comparison; birth-death models; divergence-time estimation; lineage diversification].

Formation of the Holarctic Fauna: Dated molecular phylogenetic and biogeographic insights from the Quedius-lineage of Ground-Dwelling Rove Beetles (Coleoptera, Staphylinidae)

Although the Holarctic fauna has been explored for centuries, many questions on its formation are still unanswered. For example, i) what was the impact of the uplift of the Himalaya and Tibetan Plateau?, ii) what were the timings and climate of the faunal bridges connecting the Nearctic and Palearctic regions?, and iii) how did insect lineages respond to the late Paleogene global cooling and regional aridification? To answer these, we developed a phylogenetic dataset of 1229 nuclear loci for a total of 222 species of rove beetles (Staphylinidae) with emphasis in the tribe Quediini, especially Quedius-lineage and its subclade Quedius sensu stricto. Using eight fossils for calibrating molecular clock, we estimated divergence times and then analysed in BioGeoBEARS paleodistributions of the most recent common ancestor for each target lineage. For each species we generated climatic envelopes of the temperature and precipitation and mapped them across the phylogeny to explore evolutionary shifts. Our results suggest that the warm and humid Himalaya and Tibetan Plateau acted as an evolutionary cradle for the Quedius-lineage originating during the Oligocene from where, in the Early Miocene, the ancestor of the Quedius s. str. dispersed into the West Palearctic. With the climate cooling from the Mid Miocene onwards, new lineages within Quedius s. str. emerged and gradually expanded distributions across the Palearctic. In Late Miocene, a member of the group dispersed to the Nearctic region via Beringia before the closure of this land bridge 5.3 Ma. Paleogene global cooling and regional aridification largely shaped the current biogeographic pattern for Quedius s. str. species, many of them originating during the Pliocene and shifting or contracting their ranges during Pleistocene.

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.

Modulation of T-cell function by myeloid-derived suppressor cells in hematological malignancies

Myeloid-derived suppressor cells (MDSCs) are pathologically activated neutrophils and monocytes that negatively regulate the immune response to cancer and chronic infections. Abnormal myelopoiesis and pathological activation of myeloid cells generate this heterogeneous population of myeloid-derived suppressor cells. They are characterized by their distinct transcription, phenotypic, biochemical, and functional features. In the tumor microenvironment (TME), myeloid-derived suppressor cells represent an important class of immunosuppressive cells that correlate with tumor burden, stage, and a poor prognosis. Myeloid-derived suppressor cells exert a strong immunosuppressive effect on T-cells (and a broad range of other immune cells), by blocking lymphocyte homing, increasing production of reactive oxygen and nitrogen species, promoting secretion of various cytokines, chemokines, and immune regulatory molecules, stimulation of other immunosuppressive cells, depletion of various metabolites, and upregulation of immune checkpoint molecules. Additionally, the heterogeneity of myeloid-derived suppressor cells in cancer makes their identification challenging. Overall, they serve as a major obstacle for many cancer immunotherapies and targeting them could be a favorable strategy to improve the effectiveness of immunotherapeutic interventions. However, in hematological malignancies, particularly B-cell malignancies, the clinical outcomes of targeting these myeloid-derived suppressor cells is a field that is still to be explored. This review summarizes the complex biology of myeloid-derived suppressor cells with an emphasis on the immunosuppressive pathways used by myeloid-derived suppressor cells to modulate T-cell function in hematological malignancies. In addition, we describe the challenges, therapeutic strategies, and clinical relevance of targeting myeloid-derived suppressor cells in these diseases.

A comprehensive phylogeny and revised taxonomy illuminate the origin and diversification of the global radiation of Papilio (Lepidoptera: Papilionidae)

The swallowtail genus Papilio (Lepidoptera: Papilionidae) is species rich, distributed worldwide, and has broad morphological habits and ecological niches. Because of its elevated species richness, it has been historically difficult to reconstruct a densely sampled phylogeny for this clade. Here we provide a taxonomic working list for the genus, resulting in 235 Papilio species, and assemble a molecular dataset of seven gene fragments representing ca. 80% of the currently described diversity. Phylogenetic analyses reconstructed a robust tree with highly supported relationships within subgenera, although a few nodes in the early history of the Old World Papilio remain unresolved. Contrasting with previous results, we found that Papilio alexanor is sister to all Old World Papilio and that the subgenus Eleppone is no longer monotypic. The latter includes the recently described Fijian Papilio natewa with the Australian Papilio anactus and is sister to subgenus Araminta (formerly included in subgenus Menelaides) occurring in Southeast Asia. Our phylogeny also includes rarely studied (P. antimachus, P. benguetana) or endangered species (P. buddha, P. chikae). Taxonomic changes resulting from this study are elucidated. Molecular dating and biogeographic analyses indicate that Papilio originated ca. 30 million years ago (Oligocene), in a northern region centered on Beringia. A rapid early Miocene radiation in the Paleotropics is revealed within Old World Papilio, potentially explaining their low early branch support. Most subgenera originated in the early to middle Miocene followed by synchronous southward biogeographic dispersals and repeated local extirpations in northern latitudes. This study provides a comprehensive phylogenetic framework for Papilio with clarification of subgeneric systematics and species taxonomic changes enumerated, which will facilitate further studies to address questions on their ecology and evolutionary biology using this model clade.

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.

Total evidence phylogeny of platyrrhine primates and a comparison of undated and tip-dating approaches

There have been multiple published phylogenetic analyses of platyrrhine primates (New World monkeys) using both morphological and molecular data, but relatively few that have integrated both types of data into a total evidence approach. Here, we present phylogenetic analyses of recent and fossil platyrrhines, based on a total evidence data set of 418 morphological characters and 10.2 kilobases of DNA sequence data from 17 nuclear genes taken from previous studies, using undated and tip-dating approaches in a Bayesian framework. We compare the results of these analyses with molecular scaffold analyses using maximum parsimony and Bayesian approaches, and we use a formal information theoretic approach to identify unstable taxa. After a posteriori pruning of unstable taxa, the undated and tip-dating topologies appear congruent with recent molecular analyses and support largely similar relationships, with strong support for Stirtonia as a stem alouattine, Neosaimiri as a stem saimirine, Cebupithecia as a stem pitheciine, and Lagonimico as a stem callitrichid. Both analyses find three Greater Antillean subfossil platyrrhines (Xenothrix, Antillothrix, and Paralouatta) to form a clade that is related to Callicebus, congruent with a single dispersal event by the ancestor of this clade to the Greater Antilles. They also suggest that the fossil Proteropithecia may not be closely related to pitheciines, and that all known platyrrhines older than the Middle Miocene are stem taxa. Notably, the undated analysis found the Early Miocene Panamacebus (currently recognized as the oldest known cebid) to be unstable, and the tip-dating analysis placed it outside crown Platyrrhini. Our tip-dating analysis supports a late Oligocene or earliest Miocene (20.8-27.0 Ma) age for crown Platyrrhini, congruent with recent molecular clock analyses.

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.

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.

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.

Seed predation and potential seed dispersers of the narrow endemic Ceratozamianorstogii (Zamiaceae)

In this study, we report the observation of potential seed dispersers of the endemic to Mexico and narrowly distributed Ceratozamianorstogii (Zamiaceae). Camera traps were installed in front of two plants of Ceratozamianorstogii and cone phenology until their maturity and disintegration was determined. The female cone of Ceratozamianorstogii has a development of ten months, from the time it emerges until it disintegrates. We were able to identify three stages of cone development: 1) Pre-pollination phase, 2) Pollination phase and 3) Seed maturation phase. Our results support an animal-dispersal hypothesis in Ceratozamia. Three mammals [a mouse (Pteromiscus sp.), a southern spotted skunk (Spilogaleangustifrons) and a kinkajou (Potusflavus)] were recorded biting, carrying or removing seeds of Ceratozamianorstogii. The camera traps recorded no evidence of birds or other mammals coming to the cones to feed. Thus, interaction of frugivores with seeds occurs at night. The most frequent visitor was the mouse, followed by the southern spotted skunk and the kinkajou. Significant differences (GLM, p< 0.05) in visitor frequency and time for interaction were found between species. We believe that the mouse is probably the most effective seed disperser for Ceratozamianorstogii. The results presented here have evolutionary implications that can be scaled to the entire genus Ceratozamia. Specifically, short-distance dispersal promotes allopatric speciation in this group of plants.

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.

Himalayan orogeny and monsoon intensification explain species diversification in an endemic ginger (Hedychium: Zingiberaceae) from the Indo-Malayan Realm

The Indo-Malayan Realm is a biogeographic realm that extends from the Indian Subcontinent to the islands of Southeast Asia (Malay Archipelago). Despite being megadiverse, evolutionary hypotheses explaining taxonomic diversity in this region have been rare. Here, we investigate the role of geoclimatic events such as Himalayan orogeny and monsoon intensification in the diversification of the ginger-lilies (Hedychium J.Koenig: Zingiberaceae). We first built a comprehensive, time-calibrated phylogeny of Hedychium with 75% taxonomic and geographic sampling. We found that Hedychium is a very young lineage that originated in Northern Indo-Burma, in the Late Miocene (c. 10.6 Ma). This was followed by a late Neogene and early Quaternary diversification, with multiple dispersal events to Southern Indo-Burma, Himalayas, Peninsular India, and the Malay Archipelago. The most speciose clade IV i.e., the predominantly Indo-Burmese clade also showed a higher diversification rate, suggesting its recent rapid radiation. Our divergence dating and GeoHiSSE results demonstrate that the diversification of Hedychium was shaped by both the intensifications in the Himalayan uplift as well as the Asian monsoon. Ancestral character-state reconstructions identified the occurrence of vegetative dormancy in both clades I and II, whereas the strictly epiphytic growth behavior, island dwarfism, lack of dormancy, and a distinct environmental niche were observed only in the predominantly island clade i.e., clade III. Finally, we show that the occurrence of epiphytism in clade III corresponds with submergence due to sea-level changes, suggesting it to be an adaptive trait. Our study highlights the role of recent geoclimatic events and environmental factors in the diversification of plants within the Indo-Malayan Realm and the need for collaborative work to understand biogeographic patterns within this understudied region. This study opens new perspectives for future biogeographic studies in this region and provides a framework to explain the taxonomic hyperdiversity of the Indo-Malayan Realm.

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.

Limitations of phylogenomic data can drive inferred speciation rate shifts

Biodiversity analyses of phylogenomic timetrees have produced many high-profile examples of shifts in the rate of speciation across the tree of life. Temporally correlated events in ecology, climate, and biogeography are frequently invoked to explain these rate shifts. In a re-examination of 15 genomic timetrees and 25 major published studies of the pattern of speciation through time, we observed an unexpected correlation between the timing of reported rate shifts and the information content of sequence alignments. Here, we show that the paucity of sequence variation and insufficient species sampling in phylogenomic datasets are the likely drivers of many inferred speciation rate shifts, rather than the proposed biological explanations. Therefore, data limitations can produce predictable but spurious signals of rate shifts even when speciation rates may be similar across taxa and time. Our results suggest that the reliable detection of speciation rate shifts requires the acquisition and assembly of long phylogenomic alignments with near-complete species sampling and accurate estimates of species richness for the clades of study.

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.

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.

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.].

Pulled Diversification Rates, Lineages-Through-Time Plots and Modern Macroevolutionary Modelling

Estimating time-dependent rates of speciation and extinction from dated phylogenetic trees of extant species (timetrees), and determining how and why they vary, is key to understanding how ecological and evolutionary processes shape biodiversity. Due to an increasing availability of phylogenetic trees, a growing number of process-based methods relying on the birth–death model have been developed in the last decade to address a variety of questions in macroevolution. However, this methodological progress has regularly been criticized such that one may wonder how reliable the estimations of speciation and extinction rates are. In particular, using lineages-through-time (LTT) plots, a recent study has shown that there are an infinite number of equally likely diversification scenarios that can generate any timetree. This has led to questioning whether or not diversification rates should be estimated at all. Here, we summarize, clarify, and highlight technical considerations on recent findings regarding the capacity of models to disentangle diversification histories. Using simulations, we illustrate the characteristics of newly proposed “pulled rates” and their utility. We recognize that the recent findings are a step forward in understanding the behavior of macroevolutionary modeling, but they in no way suggest we should abandon diversification modeling altogether. On the contrary, the study of macroevolution using phylogenetic trees has never been more exciting and promising than today. We still face important limitations in regard to data availability and methods, but by acknowledging them we can better target our joint efforts as a scientific community. [Birth–death models; extinction; phylogenetics; speciation.]

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.

Phylogenetic relationships amongst the snake-eyed lizards of the genus Ablepharus Fitzinger, 1823 (Sauria, Scincidae) in the Iranian Plateau based on mtDNA sequences

We recovered molecular phylogenetic relationships amongst species of the genus Ablepharus in Iran and Iraq. Partial sequences of three mitochondrial genes (cytochrome C oxidase subunit I – COI, 12S rRNA and 16S rRNA) were analysed. In addition, phylogenetic relationships and taxonomic evaluation of Ablepharus species in Cyprus, India, Greece, Turkey and Syria were performed using partial sequences of the 16S rRNA gene. Phylogenetic trees and estimated genetic distances showed that the Ablepharus populations of Iran and Iraq clustered into three distinct clades. One is found in northwest Iran (A. bivittatus in Ardabil, East and West Azerbaijan and Hamedan Provinces). The second clade, formed by A. chernovi, is found only in Uromia. The third and most heterogeneous clade is divided into two subclades, the first includes two lineages of Ablepharus in Khorasan Razavi and Semnan Provinces (A. pannonicus) and in eastern and south-eastern Iran (A. grayanus); the second subclade is distributed in the eastern part of Iraq and west and south-western Iran (Ablepharus sp.). Our analyses indicated that splitting of A. chernovi within the genus occurred in the early Miocene [about 22.5 million years ago (Mya)]. Ablepharus bivittatus diverged 15.2 Mya, in the middle Miocene. Ablepharus pannonicus diverged in the late Miocene (8.4 Mya) and A. grayanus separated in the late Miocene (6.7 Mya). The lineages of eastern Iraq and south-western Iran (Ablepharus sp.) diverged also in the late Miocene (7.0 Mya).

Sympatric and independently evolving lineages in the Thoropa miliaris - T. taophora species complex (Anura: Cycloramphidae)

Species delimitation can be challenging and affected by subjectivity. Sibling lineages that occur in sympatry constitute good candidates for species delimitation regardless of the adopted species concept. The Thoropa miliaris + T. taophora species complex exhibits high genetic diversity distributed in several lineages that occur sympatrically in the southeastern Atlantic Forest of Brazil. We used 414 loci obtained by anchored hybrid enrichment to characterize genetic variation in the Thoropa miliaris species group (T. saxatilis, T megatympanum, T. miliaris, and T. taophora), combining assignment analyses with traditional and coalescent phylogeny reconstruction. We also investigated evolutionary independence in co-occurring lineages by estimating gene flow, and validated lineages under the multispecies coalescent. We recovered most previously described lineages as unique populations in assignment analyses; exceptions include two lineages within T. miliaris that are further substructured, and the merging of all T. taophora lineages. We found very low probabilities of gene flow between sympatric lineages, suggesting independent evolution. Species tree inferences and species delimitation yielded resolved relationships and indicate that all lineages constitute putative species that diverged during the Pliocene and Pleistocene, later than previously estimated.

Understanding Diversity and Systematics in Australian Fabaceae Tribe Mirbelieae

Australia has a very diverse pea-flowered legume flora with 1715 native and naturalised species currently recognised. Tribe Mirbelieae s.l. includes 44% of Australia’s peas in 24 genera with 756 recognised species. However, several genera within the Pultenaea alliance in tribe Mirbelieae are considered to be non-monophyletic and two main options have been proposed: option one is to merge ca. 18 genera containing ca. 540 species (the largest genus, Pultenaea has nomenclatural priority); and option two is to re-circumscribe some genera and describe new genera as required to form monophyletic groups. At the species level, option one would require 76% of names to be changed; whereas based on available data, option two is likely to require, at most, 8.3% of names to change. Option two therefore provides the least nomenclatural disruption but cannot be implemented without a robust phylogenetic framework to define new generic limits. Here we present novel analyses of available plastid DNA data (trnL-F) which suggest that option two would be feasible once sufficient data are generated to resolve relationships. However, the reticulate evolutionary histories or past rapid speciation suggested for this group may prevent the resolution of all nodes. We propose targeted use of Next-Generation Sequencing technology as the best way to resolve relationships between the key clades in the tribe and present a framework for such a study. An overview of current taxonomy in the tribe is presented, along with the state of taxonomic knowledge and availability of published descriptions for electronic flora treatments. Several new combinations and typifications are published in an appendix.

Phylogenomic and Macroevolutionary Evidence for an Explosive Radiation of a Plant Genus in the Miocene

Mountain systems harbor a substantial fraction of global biodiversity and, thus, provide excellent opportunities to study rapid diversification and to understand the historical processes underlying the assembly of biodiversity hotspots. The rich biodiversity in mountains is widely regarded as having arisen under the influence of geological and climatic processes as well as the complex interactions among them. However, the relative contribution of geology and climate in driving species radiation is seldom explored. Here, we studied the evolutionary radiation of Oreocharis (Gesneriaceae), which has diversified extensively throughout East Asia, especially within the Hengduan Mountains (HDM), using transcriptomic data and a time calibrated phylogeny for 88% (111/126) of all species of the genus. In particular, we applied phylogenetic reconstructions to evaluate the extent of incomplete lineage sorting accompanying the early and rapid radiation in the genus. We then fit macroevolutionary models to explore its spatial and diversification dynamics in Oreocharis and applied explicit birth–death models to investigate the effects of past environmental changes on its diversification. Evidence from 574 orthologous loci suggest that Oreocharis underwent an impressive early burst of speciation starting ca. 12 Ma in the Miocene, followed by a drastic decline in speciation toward the present. Although we found no evidence for a shift in diversification rate across the phylogeny of Oreocharis, we showed a difference in diversification dynamics between the HDM and non-HDM lineages, with higher diversification rates in the HDM. The diversification dynamic of Oreocharis is most likely positively associated with temperature-dependent speciation and dependency on the Asian monsoons. We suggest that the warm and humid climate of the mid-Miocene was probably the primary driver of the rapid diversification in Oreocharis, while mountain building of the HDM might have indirectly affected species diversification of the HDM lineage. This study highlights the importance of past climatic changes, combined with mountain building, in creating strong environmental heterogeneity and driving diversification of mountain plants, and suggests that the biodiversity in the HDM cannot directly be attributed to mountain uplift, contrary to many recent speculations.[East Asian monsoons; environmental heterogeneity; Hengduan Mountains; incomplete lineage sorting; Oreocharis; past climate change; rapid diversification; transcriptome.]

Fast and accurate estimation of species-specific diversification rates using data augmentation

Diversification rates vary across species as a response to various factors, including environmental conditions and species-specific features. Phylogenetic models that allow accounting for and quantifying this heterogeneity in diversification rates have proven particularly useful for understanding clades diversification. Recently, we introduced the cladogenetic diversification rate shift model (ClaDS), which allows inferring multiple rate changes of small magnitude across lineages. Here we present a new inference technique for this model that considerably reduces computation time through the use of data augmentation and provide an implementation of this method in Julia. In addition to drastically reducing computation time, this new inference approach provides a posterior distribution of the augmented data, that is the tree with extinct and unsampled lineages as well as associated diversification rates. In particular, this allows extracting the distribution through time of both the mean rate and the number of lineages. We assess the statistical performances of our approach using simulations and illustrate its application on the entire bird radiation.