Phylogeny of the cycads based on multiple single-copy nuclear genes: congruence of concatenated parsimony, likelihood and species tree inference methods

The Systematics and Evolution of Gymnosperms with an Emphasis on a Few Problematic Taxa

Gymnosperms originated in the Middle Devonian and have experienced a long evolutionary history with pulses of speciation and extinction, which resulted in the four morphologically distinct extant groups, i.e., cycads, Ginkgo, conifers and gnetophytes. For over a century, the systematic relationships within the extant gymnosperms have been debated because different authors emphasized different characters. Recent phylogenomic studies of gymnosperms have given a consistent topology, which aligns well with extant gymnosperms classified into three classes, five subclasses, eight orders, and 13 families. Here, we review the historical opinions of systematics of gymnosperms with special reference to several problematic taxa and reconsider the evolution of some key morphological characters previously emphasized by taxonomists within a phylogenomic context. We conclude that (1) cycads contain two families, i.e., the Cycadaceae and the Zamiaceae; (2) Ginkgo is sister to cycads but not to conifers, with the similarities between Ginkgo and conifers being the result of parallel evolution including a monopodial growth pattern, pycnoxylic wood in long shoots, and the compound female cones, and the reproductive similarities between Ginkgo and cycads are either synapomorphic or plesiomorphic, e.g., the boat-shaped pollen, the branched pollen tube, and the flagellate sperms; (3) conifers are paraphyletic with gnetophytes nested within them, thus gnetophytes are derived conifers, and our newly delimited coniferophytes are equivalent to the Pinopsida and include three subclasses, i.e., Pinidae, Gnetidae, and Cupressidae; (4) fleshy cones of conifers originated multiple times, the Podocarpaceae are sister to the Araucariaceae, the Cephalotaxaceae and the Taxaceae comprise a small clade, which is sister to the Cupressaceae; (5) the Cephalotaxaceae are distinct from the Taxaceae, because the former family possesses typical female cones and the fleshy part of the seed is derived from the fleshiness of integument, while the latter family has reduced female cones and preserves no traces of the seed scale complexes.

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.

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.

Phylogenomics illuminates the phylogeny of flower weevils (Curculioninae) and reveals ten independent origins of brood-site pollination mutualism in true weevils

Weevils are an unusually species-rich group of phytophagous insects for which there is increasing evidence of frequent involvement in brood-site pollination. This study examines phylogenetic patterns in the emergence of brood-site pollination mutualism among one of the most speciose beetle groups, the flower weevils (subfamily Curculioninae). We analysed a novel phylogenomic dataset consisting of 214 nuclear loci for 202 weevil species, with a sampling that mainly includes flower weevils as well as representatives of all major lineages of true weevils (Curculionidae). Our phylogenomic analyses establish a uniquely comprehensive phylogenetic framework for Curculioninae and provide new insights into the relationships among lineages of true weevils. Based on this phylogeny, statistical reconstruction of ancestral character states revealed at least 10 independent origins of brood-site pollination in higher weevils through transitions from ancestral associations with reproductive structures in the larval stage. Broadly, our results illuminate the unexpected frequency with which true weevils-typically specialized phytophages and hence antagonists of plants-have evolved mutualistic interactions of ecological significance that are key to both weevil and plant evolutionary fitness and thus a component of their deeply intertwined macroevolutionary success.

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.

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.

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.

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.

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.

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.

Reassessment of Annamocarya sinesis (Carya sinensis) Taxonomy through Concatenation and Coalescence Phylogenetic Analysis

Due to its peculiar morphological characteristics, there is dispute as to whether the genus of Annamocarya sinensis, a species of Juglandaceae, is Annamocarya or Carya. Most morphologists believe it should be distinguished from the Carya genus while genomicists suggest that A. sinensis belongs to the Carya genus. To explore the taxonomic status of A. sinensis using chloroplast genes, we collected chloroplast genomes of 16 plant species and assembled chloroplast genomes of 10 unpublished Carya species. We analyzed all 26 species' chloroplast genomes through two analytical approaches (concatenation and coalescence), using the entire and unique chloroplast coding sequence (CDS) and entire and protein sequences. Our results indicate that the analysis of the CDS and protein sequences or unique CDS and unique protein sequence of chloroplast genomes shows that A. sinensis indeed belongs to the Carya genus. In addition, our analysis shows that, compared to single chloroplast genes, the phylogeny trees constructed using numerous genes showed higher consistency. Moreover, the phylogenetic analysis calculated with the coalescence method and unique gene sequences was more robust than that done with the concatenation method, particularly for analyzing phylogenetically controversial species. Through the analysis, our results concluded that A. sinensis should be called C. sinensis.

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

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

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.

Stomatal development in the cycad family Zamiaceae

BACKGROUND AND AIMS

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

METHODS

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

KEY RESULTS

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

CONCLUSIONS

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

Comparative transcriptomics provides a strategy for phylogenetic analysis and SSR marker development in Chaenomeles

The genus Chaenomeles has long been considered an important ornamental, herbal and cash crop and is widely cultivated in East Asia. Traditional studies of Chaenomeles mainly focus on evolutionary relationships at the phenotypic level. In this study, we conducted RNA-seq on 10 Chaenomeles germplasms supplemented with one outgroup species, Docynia delavayi (D. delavayi), on the Illumina HiSeq2500 platform. After de novo assemblies, we generated from 40,084 to 49,571 unigenes for each germplasm. After pairwise comparison of the orthologous sequences, 9,659 orthologues within the 11 germplasms were obtained, with 6,154 orthologous genes identified as single-copy genes. The phylogenetic tree was visualized to reveal evolutionary relationships for these 11 germplasms. GO and KEGG analyses were performed for these common single-copy genes to compare their functional similarities and differences. Selective pressure analysis based on 6,154 common single-copy genes revealed that 45 genes were under positive selection. Most of these genes are involved in building the plant disease defence system. A total of 292 genes containing simple sequence repeats (SSRs) were used to develop SSR markers and compare their functions in secondary metabolism pathways. Finally, 10 primers were chosen as SSR marker candidates for Chaenomeles germplasms by comprehensive standards. Our research provides a new methodology and reference for future related research in Chaenomeles and is also useful for improvement, breeding and selection projects in other related species.

Phylogenomics of the Andean Tetraploid Clade of the American Amaryllidaceae (Subfamily Amaryllidoideae): Unlocking a Polyploid Generic Radiation Abetted by Continental Geodynamics

One of the two major clades of the endemic American Amaryllidaceae subfam. Amaryllidoideae constitutes the tetraploid-derived (n = 23) Andean-centered tribes, most of which have 46 chromosomes. Despite progress in resolving phylogenetic relationships of the group with plastid and nrDNA, certain subclades were poorly resolved or weakly supported in those previous studies. Sequence capture using anchored hybrid enrichment was employed across 95 species of the clade along with five outgroups and generated sequences of 524 nuclear genes and a partial plastome. Maximum likelihood phylogenetic analyses were conducted on concatenated supermatrices, and coalescent-based species tree analyses were run on the gene trees, followed by hybridization network, age diversification and biogeographic analyses. The four tribes Clinantheae, Eucharideae, Eustephieae, and Hymenocallideae (sister to Clinantheae) are resolved in all analyses with > 90 and mostly 100% support, as are almost all genera within them. Nuclear gene supermatrix and species tree results were largely in concordance; however, some instances of cytonuclear discordance were evident. Hybridization network analysis identified significant reticulation in Clinanthus, Hymenocallis, Stenomesson and the subclade of Eucharideae comprising Eucharis, Caliphruria, and Urceolina. Our data support a previous treatment of the latter as a single genus, Urceolina, with the addition of Eucrosia dodsonii. Biogeographic analysis and penalized likelihood age estimation suggests an origin in the Cauca, Desert and Puna Neotropical bioprovinces for the complex in the mid-Oligocene, with more dispersals than vicariances in its history, but no extinctions. Hymenocallis represents the only instance of long-distance vicariance from the tropical Andean origin of its tribe Hymenocallideae. The absence of extinctions correlates with the lack of diversification rate shifts within the clade. The Eucharideae experienced a sudden lineage radiation ca. 10 Mya. We tie much of the divergences in the Andean-centered lineages to the rise of the Andes, and suggest that the Amotape-Huancabamba Zone functioned as both a corridor (dispersal) and a barrier to migration (vicariance). Several taxonomic changes are made. This is the largest DNA sequence data set to be applied within Amaryllidaceae to date.

Modeling the Distribution of the Invasive Alien Cycad Aulacaspis Scale in Africa Under Current and Future Climate Scenarios

The cycad aulacaspis scale, Aulacaspis yasumatsui Takagi (Hemiptera: Coccoidea: Diaspididae), is native to Southeast Asia but an invasive pest of the gymnosperm order Cycadales in many parts of the world. Aulacaspis yasumatsui was recently reported on the cycad genus Encephalartos in South Africa and is currently categorized as a 'prohibited terrestrial invertebrate' in the invasive species legislation, National Environmental Management: Biodiversity Act, 2004 (NEM:BA). Encephalartos is endemic to Africa, and 11 species are listed as critically endangered and four species as endangered. Seeing the limited distribution of A. yasumatsui in South Africa and only one unconfirmed record from the Ivory Coast, understanding the potential distribution range is essential for control and management. Here we model the potential distribution of A. yasumatsui under current and future climate scenarios in Africa, with a focus on South Africa. Future climatic scenarios were simulated using a bio-climatic software, CLIMEX. The model indicates that, under the current climatic scenario, all 17 African countries possessing Encephalartos are susceptible to A. yasumatsui establishment. However, under climatic change, the suitability decreases for large parts of Africa. In South Africa, 93% of the winter rainfall areas, and 90% of the temperate, summer rainfall areas are suitable for A. yasumatsui establishment. In this study, we highlight the urgent need for regulation, management, and research on A. yasumatsui in African countries with native cycads.

Stem Branching of Cycad Plants Informs Horticulture and Conservation Decisions

The number of branches in male and female plants of Cycas micronesica K.D. Hill, Cycas edentata de Laub., Cycas wadei Merr., and Zamia encephalartoides D.W. Stev. were counted in Guam, Philippines, and Colombia, to confirm earlier reports that female plants develop fewer branches than males. Cycas plants produce determinate male strobili and indeterminate female strobili, but Zamia plants produce determinate strobili for both sexes. More than 80% of the female trees for each of the Cycas species were unbranched with a single stem, but more than 80% of the male trees exhibited two or more branches. The mean number of branches on male plants was more than double that of female plants. The number of branches of the Zamia male plants was almost triple that of female plants. Moreover, the Zamia plants produced 2.8-fold greater numbers of branches than the mean of the Cycas plants. Most of Guam’s unsexed C. micronesica trees in 2004 were unbranched, but after 15 years of damage from non-native insect herbivores, most of the remaining live trees in 2020 contained three or more branches. The results confirm that male Cycas and Zamia plants produce more branches than female plants and suggest cycad species with determinate female strobili produce more branches on female plants than species with indeterminate female strobili. Our results indicate that the years of plant mortality on Guam due to non-native insect herbivores have selectively killed more female C. micronesica trees. Horticulture and conservation decisions may be improved with this sexual dimorphism knowledge.

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

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

An ancient push-pull pollination mechanism in cycads

Most cycads engage in brood-site pollination mutualisms, yet the mechanism by which the Cycadales entice pollination services from diverse insect mutualists remains unknown. Here, we characterize a push-pull pollination mechanism between a New World cycad and its weevil pollinators that mirrors the mechanism between a distantly related Old World cycad and its thrips pollinators. The behavioral convergence between weevils and thrips, combined with molecular phylogenetic dating and a meta-analysis of thermogenesis and coordinated patterns of volatile attraction and repulsion suggest that a push-pull pollination mutualism strategy is ancestral in this ancient, dioecious plant group. Hence, it may represent one of the earliest insect/plant pollination mechanisms, arising long before the evolution of visual floral signaling commonly used by flowering plants.

Genetic characterization of the entire range of Cycas panzhihuaensis (Cycadaceae)

Cycas panzhihuaensis L. Zhou & S. Y. Yang (Cycadaceae) is an endangered gymnosperm species endemic to the dry-hot valley of the Jinsha River basin in southwest China. Although the wild C. panzhihuaensis population from Panzhihua Cycad Natural Reserve is well protected and its genetic diversity has been well assessed, the genetic characteristics of populations outside the nature reserve, which face larger risks of extinction, remain unknown. Furthermore, the population genetics and historical dynamics of this endemic and endangered species have not been examined across its entire range. In this study, to analyze the genetic diversity, phylogeographical structure and demographic history of C. panzhihuaensis from all its seven known locations, we sequenced and compared molecular data from chloroplastic DNA (psbA-trnH, psbM-trnD, and trnS-trnG), single-copy nuclear genes (PHYP, AC5, HSP70, and AAT) from 61 individuals, as well as 11 nuclear microsatellite loci (SSR) from 102 individuals. We found relatively high genetic diversity within populations and high genetic differentiation among populations of C. panzhihuaensis, which is consistent with the patterns of other Asian inland cycads. Although no significant phylogeographical structure was detected, we found that small and unprotected populations possess higher genetic diversity and more unique haplotypes, which revises our understanding of diversity within this species and deserves due attention. Analysis of demographic dynamics suggest that human activity might be the key threat to C. panzhihuaensis. Based on the genetic characterization of C. panzhihuaensis, we propose several practical guidelines for the conservation of this species, especially for the populations with small sizes.

Target Nuclear and Off-Target Plastid Hybrid Enrichment Data Inform a Range of Evolutionary Depths in the Orchid Genus Epidendrum

Universal angiosperm enrichment probe sets designed to enrich hundreds of putatively orthologous nuclear single-copy loci are increasingly being applied to infer phylogenetic relationships of different lineages of angiosperms at a range of evolutionary depths. Studies applying such probe sets have focused on testing the universality and performance of the target nuclear loci, but they have not taken advantage of off-target data from other genome compartments generated alongside the nuclear loci. Here we do so to infer phylogenetic relationships in the orchid genus Epidendrum and closely related genera of subtribe Laeliinae. Our aims are to: 1) test the technical viability of applying the plant anchored hybrid enrichment (AHE) method (Angiosperm v.1 probe kit) to our focal group, 2) mine plastid protein coding genes from off-target reads; and 3) evaluate the performance of the target nuclear and off-target plastid loci in resolving and supporting phylogenetic relationships along a range of taxonomical depths. Phylogenetic relationships were inferred from the nuclear data set through coalescent summary and site-based methods, whereas plastid loci were analyzed in a concatenated partitioned matrix under maximum likelihood. The usefulness of target and flanking non-target nuclear regions and plastid loci was assessed through the estimation of their phylogenetic informativeness. Our study successfully applied the plant AHE probe kit to Epidendrum, supporting the universality of this kit in angiosperms. Moreover, it demonstrated the feasibility of mining plastome loci from off-target reads generated with the Angiosperm v.1 probe kit to obtain additional, uniparentally inherited sequence data at no extra sequencing cost. Our analyses detected some strongly supported incongruences between nuclear and plastid data sets at shallow divergences, an indication of potential lineage sorting, hybridization, or introgression events in the group. Lastly, we found that the per site phylogenetic informativeness of the ycf1 plastid gene surpasses that of all other plastid genes and several nuclear loci, making it an excellent candidate for assessing phylogenetic relationships at medium to low taxonomic levels in orchids.

The complete chloroplast genome of Microcycas calocoma (Miq.) A. DC. (Zamiaceae, Cycadales) and evolution in Cycadales

Cycadales is an extant group of seed plants occurring in subtropical and tropical regions comprising putatively three families and 10 genera. At least one complete plastid genome sequence has been reported for all of the 10 genera except Microcycas, making it an ideal plant group to conduct comprehensive plastome comparisons at the genus level. This article reports for the first time the plastid genome of Microcycas calocoma. The plastid genome has a length of 165,688 bp with 134 annotated genes including 86 protein-coding genes, 47 non-coding RNA genes (39 tRNA and eight rRNA) and one pseudogene. Using global sequence variation analysis, the results showed that all cycad genomes share highly similar genomic profiles indicating significant slow evolution and little variation. However, identity matrices coinciding with the inverted repeat regions showed fewer similarities indicating that higher polymorphic events occur at those sites. Conserved non-coding regions also appear to be more divergent whereas variations in the exons were less discernible indicating that the latter comprises more conserved sequences. Phylogenetic analysis using 81 concatenated protein-coding genes of chloroplast (cp) genomes, obtained using maximum likelihood and Bayesian inference with high support values (>70% ML and = 1.0 BPP), confirms that Microcycas is closest to Zamia and forms a monophyletic clade with Ceratozamia and Stangeria. While Stangeria joined the Neotropical cycads Ceratozamia, Zamia and Microcyas, Bowenia grouped with the Southern Hemisphere cycads Encephalartos, Lepidozamia and Macrozamia. All Cycas species formed a distinct clade separated from the other genera. Dioon, on the other hand, was outlying from the rest of Zamiaceae encompassing two major clades—the Southern Hemisphere cycads and the Neotropical cycads. Analysis of the whole cp genomes in phylogeny also supports that the previously recognized family—Stangeriaceae—which contained Bowenia and Stangeria, is not monophyletic. Thus, the cp genome topology obtained in our study is congruent with other molecular phylogenies recognizing only a two-family classification (Cycadaceae and Zamiaceae) within extant Cycadales.

Population Differentiation and Demographic History of the Cycas taiwaniana Complex (Cycadaceae) Endemic to South China as Indicated by DNA Sequences and Microsatellite Markers

Historical geology, climatic oscillations, and seed dispersal capabilities are thought to influence the population dynamics and genetics of plants, especially for distribution-restricted and threatened species. Investigating the genetic resources within and among taxa is a prerequisite for conservation management. The Cycas taiwaniana complex consists of six endangered species that are endemic to South China. In this study, we investigated the relationship between phylogeographic history and the genetic structure of the C. taiwaniana complex. To estimate the phylogeographic history of the complex, we assessed the genetic structure and divergence time, and performed phylogenetic and demographic historical analyses. Two chloroplast DNA intergenic regions (cpDNA), two single-copy nuclear genes (SCNGs), and six microsatellite loci (SSR) were sequenced for 18 populations. The SCNG data indicated a high genetic diversity within populations, a low genetic diversity among populations, and significant genetic differentiation among populations. Significant phylogeographical structure was detected. Structure and phylogenetic analyses both revealed that the 18 populations of the C. taiwaniana complex have two main lineages, which were estimated to diverge in the Middle Pleistocene. We propose that Cycas fairylakea was incorporated into Cycas szechuanensis and that the other populations, which are mainly located on Hainan Island, merged into one lineage. Bayesian skyline plot analyses revealed that the C. taiwaniana complex experienced a recent decline, suggesting that the complex probably experienced a bottleneck event. We infer that the genetic structure of the C. taiwaniana complex has been affected by Pleistocene climate shifts, sea-level oscillations, and human activities. In addition to providing new insights into the evolutionary legacy of the genus, the genetic characterizations will be useful for the conservation of Cycas species.

Characterization of the complete chloroplast genome of Microcycas calocoma (Zamiaceae), an Endangered monotypic cycad species from Cuba

Microcycas calocoma is the monotypic species from the critical endangered and endemic cycad genus Microcycas in Cuba, an important taxon to study the evolution of extant gymnosperms. Here we report the complete chloroplast sequences of M. calocoma and characterize the genome structure of this species. The genome size of M. calocoma is 165,667 bp in length which contains 135 genes, including 88 protein-coding genes, eight ribosomal RNA genes and 39 transfer RNA genes. The GC content of this genome is 39.6%. Phylogenomic study shows M. calocoma is mostly closely related to the cycad genus Zamia, which corresponds to previous studies based on multiple nuclear genes and combined plastid and nuclear evidence. The plastome information of M. calocoma offered by this study can contribute to further comparative chloroplast genome in cycads/gymnosperms as well as conservation genetic studies of M. calocoma.

Transcriptome analysis of two radiated Cycas species and the subsequent species delimitation of the Cycas taiwaniana complex

PREMISE

Cycas is an important gymnosperm genus, and the most diverse of all cycad genera. The C. taiwaniana complex of species are morphologically similar and difficult to distinguish due to a lack of genomic resources.

METHODS

We characterized the transcriptomes of two closely related and endangered Cycas species endemic to Hainan, China: C. hainanensis and C. changjiangensis. Three single-copy nuclear genes in the C. taiwaniana complex were sequenced based on these transcriptomes, enabling us to evaluate the species boundaries using the multispecies coalescent method implemented in the Bayesian Phylogenetics and Phylogeography program.

RESULTS

We obtained 68,184 and 81,561 unigenes for C. changjiangensis and C. hainanensis, respectively. We identified six positively selected genes that are mainly involved in stimulus responses, suggesting that environmental adaptation may have played an important role in the relatively recent divergence of these species. The similar K S distribution peaks at 1.0 observed for the paralogs in the two species indicate a common whole-genome duplication event. Our species delimitation analysis indicated that the C. taiwaniana complex consists of three distinct species, which correspond to the previously reported morphological differences.

DISCUSSION

Our study provides valuable genetic resources for Cycas species and guidance for the taxonomic treatment of the C. taiwaniana complex, as well as new insights into evolution of species within Cycas.

First cycad seedling foliage from the fossil record and inferences for the Cenozoic evolution of cycads

The morphology of the early ontogenetic stages of cycad foliage may help resolve the relationships between extinct to extant cycad lineages. However, prior to this study, fossil evidence of cycad seedlings was not known. We describe a compression fossil of cycad eophylls with co-occurring fully developed leaves of adult specimens from the early Palaeocene (ca 63.8 Ma) Castle Rock flora from the Denver Basin, CO, USA and assign it to the fossil genus Dioonopsis (Cycadales) based on leaf morphology and anatomy. The new fossil seedling foliage is particularly important because fully differentiated pinnate leaves of adult plants and the eophylls belong to the same species based on shared epidermal micromorphology, therefore, increasing the number of morphological characteristics that can be used to place Dioonopsis phylogenetically. Significantly, the seedling fossil has a basic foliage structure that is very similar to seedlings of extant cycads, which is consistent with a cycadalean affinity of Dioonopsis. Nevertheless, the set of morphological characters in the seedling and adult specimens of Dioonopsis suggests a distant relationship between Dioonopsis and extant Dioon. This indicates that extinct lineages of cycads were present and widespread during the early Cenozoic (Palaeogene) coupled with the subordinate role of extant genera in the Palaeogene fossil record of cycads.

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

PREMISE OF THE STUDY

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

METHODS

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

KEY RESULTS

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

CONCLUSIONS

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

Beetle Pollination of Cycads in the Mesozoic

Cycads, unlike modern wind-pollinated conifers and Ginkgo, are unusual in that they are an ancient group of gymnosperms pollinated by insects [1-3]. Although it is well documented that cycads were diverse and abundant during the mid-Mesozoic, little is known about their biogeography and pollination before the rise of angiosperms. Direct fossil evidence illuminating the evolutionary history of cycads is extremely rare [4, 5]. Here we report a specialized beetle-mediated pollination mode from the mid-Cretaceous of Myanmar, wherein a new boganiid beetle, Cretoparacucujus cycadophilus, with specialized pollen-feeding adaptations in its mouthparts and legs, was associated with many pollen grains of Cycadopites. Phylogenetic analyses indicate Cretoparacucujus as a sister group to the extant Australian Paracucujus, which pollinate the cycad Macrozamia riedlei. Our discovery, along with the current disjunct distribution of related beetle-herbivore (tribe Paracucujini) and cycad-host (tribe Encephalarteae) pairs in South Africa and Australia, indicate a probable ancient origin of beetle pollination of cycads at least in the Early Jurassic, long before angiosperm dominance and the radiation of flowering-plant pollinators later in the Cretaceous.

Comparative Biology of Cycad Pollen, Seed and Tissue - A Plant Conservation Perspective

Cycads are the most endangered of plant groups based on IUCN Red List assessments; all are in Appendix I or II of CITES, about 40% are within biodiversity 'hotspots,' and the call for action to improve their protection is long-standing. We contend that progress in this direction will not be made until there is better understanding of cycad pollen, seed and tissue biology, which at the moment is limited to relatively few (<10%) species. We review what is known about germplasm (seed and pollen) storage and germination, together with recent developments in the application of contemporary technologies to tissues, such as isotype labelling, biomolecular markers and tissue culture. Whilst progress is being made, we conclude that an acceleration of comparative studies is needed to facilitate the integration of in situ and ex situ conservation programmes to better safeguard endangered cycads.

Predicting the potential distributions of the invasive cycad scale Aulacaspis yasumatsui (Hemiptera: Diaspididae) under different climate change scenarios and the implications for management

Cycads are an ancient group of gymnosperms that are popular as landscaping plants, though nearly all of them are threatened or endangered in the wild. The cycad aulacaspis scale (CAS), Aulacaspis yasumatsui Takagi (Hemiptera: Diaspididae), has become one of the most serious pests of cycads in recent years; however, the potential distribution range and the management approach for this pest are unclear. A potential risk map of cycad aulacaspis scale was created based on occurrence data under different climatic conditions and topology factors in this study. Furthermore, the future potential distributions of CAS were projected for the periods 2050s and 2070s under three different climate change scenarios (GFDL-CM3, HADGEM2-AO and MIROC5) described in the Special Report on Emissions Scenarios of the IPCC (Intergovernmental Panel on Climate Change). The model suggested high environmental suitability for the continents of Asia and North America, where the species has already been recorded. The potential distribution expansions or reductions were also predicted under different climate change conditions. Temperature of Driest Quarter (Bio9) was the most important factor, explaining 48.1% of the distribution of the species. The results also suggested that highly suitable habitat for CAS would exist in the study area if the mean temperature of 15–20 °C in the driest quarter and a mean temperature of 25–28 °C the wettest quarter. This research provides a theoretical reference framework for developing policy to manage and control this invasive pest.

Phylogeny of the gymnosperm genus Cycas L. (Cycadaceae) as inferred from plastid and nuclear loci based on a large-scale sampling: Evolutionary relationships and taxonomical implications

The gymnosperm genus Cycas is the sole member of Cycadaceae, and is the largest genus of extant cycads. There are about 115 accepted Cycas species mainly distributed in the paleotropics. Based on morphology, the genus has been divided into six sections and eight subsections, but this taxonomy has not yet been tested in a molecular phylogenetic framework. Although the monophyly of Cycas is broadly accepted, the intrageneric relationships inferred from previous molecular phylogenetic analyses are unclear due to insufficient sampling or uninformative DNA sequence data. In this study, we reconstructed a phylogeny of Cycas using four chloroplast intergenic spacers and seven low-copy nuclear genes and sampling 90% of extant Cycas species. The maximum likelihood and Bayesian inference phylogenies suggest: (1) matrices of either concatenated cpDNA markers or of concatenated nDNA lack sufficient informative sites to resolve the phylogeny alone, however, the phylogeny from the combined cpDNA-nDNA dataset suggests the genus can be roughly divided into 13 clades and six sections that are in agreement with the current classification of the genus; (2) although with partial support, a clade combining sections Panzhihuaenses + Asiorientales is resolved as the earliest diverging branch; (3) section Stangerioides is not monophyletic because the species resolve as a grade; (4) section Indosinenses is not monophyletic as it includes Cycas macrocarpa and C. pranburiensis from section Cycas; (5) section Cycas is the most derived group and its subgroups correspond with geography.

Pollination of the Australian cycad Cycas ophiolitica (Cycadaceae): the limited role of wind pollination in a cycad with beetle pollinator mutualists, and its ecological significance

Cycads in the Zamiaceae are well known for their host-specific insect pollination mutualisms. Pollination of Cycas in the sister family Cycadaceae is less well-documented, with beetle pollination possibly coexisting with a limited potential for wind pollination, a hypothesis we tested for C. ophiolitica in Central Queensland, Australia. Cones were associated with three species of beetle: an undescribed weevil (Curculionidae), Hapalips sp. (Erotylidae) and Ulomoides sp. (Tenebrionidae). Pollination-vector exclusion experiments compared the pollination success (quantified as % ovules pollinated per cone) of control cones against bagged or netted cones that excluded wind or insects respectively (n = 10 for all treatments). Insects do pollinate C. ophiolitica in the absence of wind, the median (first quartile-third quartile) pollination success of control plants being 83.7% (60.8–87.2%) while bagged cones, from which wind, but not insects, were excluded, pollinated at 52.9% (19.5–74.8%). For netted cones, (excluding insects but not wind), pollination fell to 12.6% (10.9–45.9%). Airborne pollen (as quantified by capture on a series of adhesive pollen traps) decreased rapidly with distance from male cones, potentially becoming ineffective for wind pollination at ~5 m. Airborne pollen load in the vicinity of female cones, and distance of females from neighbouring males, suggests wind pollination may occur sporadically, but only at high spatial densities. Although Cycas appears to be primarily insect pollinated, this limited potential for ambophily may be significant given the history of dispersal and pollinator host shifts among these cycads.

Cracking the omega code: hydraulic architecture of the cycad leaf axis

Background and Aims

The leaf axis of members of the order Cycadales ('cycads') has long been recognized by its configuration of independent vascular bundles that, in transverse section, resemble the Greek letter omega (hence the 'omega pattern'). This provides a useful diagnostic character for the order, especially when applied to paleobotany. The function of this pattern has never been elucidated. Here we provide a three-dimensional analysis and explain the pattern in terms of the hydraulic architecture of the pinnately compound cycad leaf.

Methods

The genus Cycas was used as a simple model, because each leaflet is supplied by a single vascular bundle. Sequential sectioning was conducted throughout the leaf axis and photographed with a digital camera. Photographs were registered and converted to a cinematic format, which provided an objective method of analysis.

Key Results

The omega pattern in the petiole can be sub-divided into three vascular components, an abaxial 'circle', a central 'column' and two adaxial 'wings', the last being the only direct source of vascular supply to the leaflets. Each leaflet is supplied by a vascular bundle that has divided or migrated directly from the closest wing bundle. There is neither multiplication nor anastomoses of vascular bundles in the other two components. Thus, as one proceeds from base to apex along the leaf axis, the number of vascular bundles in circle and column components is reduced distally by their uniform migration throughout all components. Consequently, the distal leaflets are irrigated by the more abaxial bundles, guaranteeing uniform water supply along the length of the axis.

Conclusions

The omega pattern exemplifies one of the many solutions plants have achieved in supplying distal appendages of an axis with a uniform water supply. Our method presents a model that can be applied to other genera of cycads with more complex vascular organization.

The phylogeography of the cycad genus Dioon (Zamiaceae) clarifies its Cenozoic expansion and diversification in the Mexican transition zone

Background and Aims

Biogeographic transition zones are promising areas to study processes of biogeographic evolution and its influence on biological groups. The Mexican transition zone originated due to the overlap of Nearctic and Neotropical biota, which promoted great biological diversification. However, since most previous studies in this area were focused on revealing the phylogeography of Nearctic plants, how historical biogeographic configuration influenced the expansion and diversification of the Neotropical flora remains almost unknown. Using the cycad genus Dioon (Zamiaceae), this study aimed to test whether the biogeographic provinciality of the Mexican transition zone reflects the history of diversification of Neotropical plants.

Methods

Two chloroplast DNA (cpDNA) regions were analysed from 101 specimens of 15 Dioon species to reveal the distribution of haplogroups. In addition, genome-wide single nucleotide polymorphisms (SNPs) from 84 specimens were used to test the concordance between phylogenetic clusters and the biogeographic provinces. An ultrametric tree was constructed from the sequences containing SNPs to reconstruct the biogeographic events of vicariance and dispersal of Dioon across the Neotropical biogeographic provinces.

Key Results

Four Dioon lineages with strong phylogeographic structures were recognized using both cpDNA and SNP data. The lineages correspond to two clades that originated from a common ancestor in Eastern Mexico. One clade expanded and diversified in South-east Mexico and Central America. Another clade diversified into three lineages that dispersed to North-east, South and North-west Mexico. Each lineage was biogeographically delimitated. Biogeographic provinces might have provided disparate ecological conditions that facilitated speciation in Dioon since the Miocene.

Conclusions

The current genetic structure and species diversity of Dioon depict the history of expansion and diversification of the northernmost Neotropical provinces. Past biogeographic connectivities were favoured by elevated topographies, since mountain systems served as corridors for the migration of Dioon and as refugia of tropical communities that diversified during the formation of modern Neotropical forests.

Aridification as a driver of biodiversity: a case study for the cycad genus Dioon (Zamiaceae)

Background and Aims

Aridification is considered a selective pressure that might have influenced plant diversification. It is suggested that plants adapted to aridity diversified during the Miocene, an epoch of global aridification (≈15 million years ago). However, evidence supporting diversification being a direct response to aridity is scarce, and multidisciplinary evidence, besides just phylogenetic estimations, is necessary to support the idea that aridification has driven diversification. The cycad genus Dioon (Zamiaceae), a tropical group including species occurring from humid forests to arid zones, was investigated as a promising study system to understand the associations among habitat shifts, diversification times, the evolution of leaf epidermal adaptations, and aridification of Mexico.

Methods

A phylogenetic tree was constructed from seven chloroplast DNA sequences and the ITS2 spacer to reveal the relationships among 14 Dioon species from habitats ranging from humid forests to deserts. Divergence times were estimated and the habitat shifts throughout Dioon phylogeny were detected. The epidermal anatomy among Dioon species was compared and correlation tests were performed to associate the epidermal variations with habitat parameters.

Key Results

Events of habitat shifts towards arid zones happened exclusively in one of the two main clades of Dioon. Such habitat shifts happened during the species diversification of Dioon, mainly during the Miocene. Comparative anatomy showed epidermal differences between species from arid and mesic habitats. The variation of epidermal structures was found to be correlated with habitat parameters. Also, most of the analysed epidermal traits showed significant phylogenetic signals.

Conclusions

The diversification of Dioon has been driven by the aridification of Mexico. The Miocene timing corresponds to the expansion of arid zones that embedded the ancestral Dioon populations. As response, species in arid zones evolved epidermal traits to counteract aridity stress. This case study provides a robust body of evidence supporting the idea that aridification is an important driver of biodiversity.

Historical Biogeography of endemic seed plant genera in the Caribbean: Did GAARlandia play a role?

The Caribbean archipelago is a region with an extremely complex geological history and an outstanding plant diversity with high levels of endemism. The aim of this study was to better understand the historical assembly and evolution of endemic seed plant genera in the Caribbean, by first determining divergence times of endemic genera to test whether the hypothesized Greater Antilles and Aves Ridge (GAARlandia) land bridge played a role in the archipelago colonization and second by testing South America as the main colonization source as expected by the position of landmasses and recent evidence of an asymmetrical biotic interchange. We reconstructed a dated molecular phylogenetic tree for 625 seed plants including 32 Caribbean endemic genera using Bayesian inference and ten calibrations. To estimate the geographic range of the ancestors of endemic genera, we performed a model selection between a null and two complex biogeographic models that included timeframes based on geological information, dispersal probabilities, and directionality among regions. Crown ages for endemic genera ranged from Early Eocene (53.1 Ma) to Late Pliocene (3.4 Ma). Confidence intervals for divergence times (crown and/or stem ages) of 22 endemic genera occurred within the GAARlandia time frame. Contrary to expectations, the Antilles appears as the main ancestral area for endemic seed plant genera and only five genera had a South American origin. In contrast to patterns shown for vertebrates and other organisms and based on our sampling, we conclude that GAARlandia did not act as a colonization route for plants between South America and the Antilles. Further studies on Caribbean plant dispersal at the species and population levels will be required to reveal finer-scale biogeographic patterns and mechanisms.

Shifting Quaternary migration patterns in the Bahamian archipelago: Evidence from the Zamia pumila complex at the northern limits of the Caribbean island biodiversity hotspot

PREMISE OF THE STUDY

The Bahamas archipelago is formed by young, tectonically stable carbonate banks that harbor direct geological evidence of global ice-volume changes. We sought to detect signatures of major changes on gene flow patterns and reconstruct the phylogeographic history of the monophyletic Zamia pumila complex across the Bahamas.

METHODS

Nuclear molecular markers with both high and low mutation rates were used to capture two different time scale signatures and test several gene flow and demographic hypotheses.

KEY RESULTS

Single-copy nuclear genes unveiled apparent ancestral admixture on Andros, suggesting a significant role of this island as main hub of diversity of the archipelago. We detected demographic and spatial expansion of the Zamia pumila complex on both paleo-provinces around the Piacenzian (Pliocene)/Gelasian (Pleistocene). Populations evidenced signatures of different migration models that have occurred at two different times. Populations on Long Island (Z. lucayana) may either represent a secondary colonization of the Bahamas by Zamia or a rapid and early-divergence event of at least one population on the Bahamas.

CONCLUSIONS

Despite changes in migration patterns with global climate, expected heterozygosity with both marker systems remains within the range reported for cycads, but with significant levels of increased inbreeding detected by the microsatellites. This finding is likely associated with reduced gene flow between and within paleo-provinces, accompanied by genetic drift, as rising seas enforced isolation. Our study highlights the importance of the maintenance of the predominant direction of genetic exchange and the role of overseas dispersion among the islands during climate oscillations.

Eobowenia gen. nov. from the Early Cretaceous of Patagonia: indication for an early divergence of Bowenia?

BACKGROUND

Even if they are considered the quintessential "living fossils", the fossil record of the extant genera of the Cycadales is quite poor, and only extends as far back as the Cenozoic. This lack of data represents a huge hindrance for the reconstruction of the recent history of this important group. Among extant genera, Bowenia (or cuticles resembling those of extant Bowenia) has been recorded in sediments from the Late Cretaceous and the Eocene of Australia, but its phylogenetic placement and the inference from molecular dating still imply a long ghost lineage for this genus.

RESULTS

We re-examine the fossil foliage Almargemia incrassata from the Lower Cretaceous Anfiteatro de Ticó Formation in Patagonia, Argentina, in the light of a comparative cuticular analysis of extant Zamiaceae. We identify important differences with the other member of the genus, viz. A. dentata, and bring to light some interesting characters shared exclusively between A. incrassata and extant Bowenia. We interpret our results to necessitate the erection of the new genus Eobowenia to accommodate the fossil leaf earlier assigned as Almargemia incrassata. We then perfom phylogenetic analyses, including the first combined morphological and molecular analysis of the Cycadales, that indicate that the newly erected genus could be related to extant Bowenia.

CONCLUSION

Eobowenia incrassata could represent an important clue for the understanding of evolution and biogeography of the extant genus Bowenia, as the presence of Eobowenia in Patagonia is yet another piece of the biogeographic puzzle that links southern South America with Australasia.

New insights on the origin of the woody flora of the Chihuahuan Desert: The case of Lindleya

PREMISE OF THE STUDY

We reconstructed aspects of the evolutionary history of Lindleya, a shrubby element in the rose family, as a model for identifying the ancestral area of the biota of the Chihuahuan Desert and for understanding the effect that pre-Pleistocene or Pleistocene climate fluctuations had on the evolution of the plants of this desert.

METHODS

We conducted phylogenetic, biogeographic, and phylogeographic analyses based on one nuclear DNA and two plastid markers sequenced for 20 populations comprising 153 accessions and representing the species' entire distribution range. We carried out ecological niche modeling to forecast Lindleya's potential distribution under Current, Last Interglacial, and Last Glacial Maximum conditions.

KEY RESULTS

The ancestral area of Lindleya populations was probably in the central and southern parts of the Chihuahuan Desert. Levels of genetic and haplotype variation were the highest in the northernmost and southernmost areas and low levels of genetic variation in the central part. Extended Bayesian skyline plots including all populations identified a slight demographic expansion ∼1 Ma. The predicted potential distribution during the Last Glacial Maximum was very restricted.

CONCLUSIONS

We suggest that pre-Pleistocene changes had an influence on the divergence of Lindleya populations. We also suggest that the ancestral area for Lindleya was the Chihuahuan Desert and that the uplift of the Trans-Mexican Volcanic Belt isolated the Tehuacán Valley populations. Genetic and demographic analyses as well as and ecological niche modeling indicate that populations of Lindleya experienced genetic bottlenecks and that they have expanded since the Last Glacial Maximum.