A phylogeny of cycads (Cycadales) inferred from chloroplast matK gene, trnK intron, and nuclear rDNA ITS region

Cycads: A comprehensive review of its botany, traditional uses, phytochemistry, pharmacology and toxicology

Cycads, which primarily consist of the families Cycadaceae and Zamiaceae, possess intrinsic therapeutic attributes that are prominently expressed across their morphological spectrum, including roots, leaves, flowers, and seeds. In Chinese traditional medicine, the leaves of cycads are particularly revered for their profound healing capabilities. This meticulous review engages with existing literature on cycads and presents insightful avenues for future research. Over 210 phytoconstituents have been isolated and identified from various cycad tissues, including flavonoids, azoxy metabolites, sterols, lignans, non-proteogenic amino acids, terpenoids, and other organic constituents. The contemporary pharmacological discourse highlights the antineoplastic, antimicrobial, and antidiabetic activities inherent in these ancient plants, which are of particular importance to the field of oncology. Despite the prevalent focus on crude extracts and total flavonoid content, our understanding of the nuanced pharmacodynamics of cycads lags considerably behind. The notoriety of cycads derived toxicity, notably within the context of Guam's neurological disease cluster, has precipitated an established emphasis on toxicological research within this field. As such, this critical review emphasizes nascent domains deserving of academic and clinical pursuit, whilst nested within the broader matrix of current scientific understanding. The systematic taxonomy, traditional applications, phytochemical composition, therapeutic potential, and safety profile of cycads are holistically interrogated, assimilating an indispensable repository for future scholarly inquiries. In conclusion, cycads stand as a veritable treasure trove of pharmacological virtue, displaying remarkable therapeutic prowess and holding vast promise for ongoing scientific discovery and clinical utilization.

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.

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.

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.

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.

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.

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.

Sex-specific differences in functional traits and resource acquisition in five cycad species

Selective pressures acting on plant life histories can drive extreme specialization. One example of such specialization is the evolution of dioecious breeding systems. Evolutionary and ecological theory posits that dioecy may subject male and female individuals to different selective pressures and result in unique sex-mediated adaptive traits related to resource allocation and ecophysiology. Cycads are the earliest diverging lineage of seed plants with strict dioecy, yet we know almost nothing about the ecology and physiology of this group. Especially limited is our understanding of potential sex-specific differences and how such differences may influence species ecology. Here we examine the ecophysiology of male and female cycads to understand better, the role that dioecy plays in this group. We evaluated sex-specific differences in ecophysiological traits and resource acquisition in five species. Specifically, we compared photosynthetic physiology, nitrogen and carbon content, isotope discrimination (δ15N and δ13C), and stomatal density. In some cycads, (i) males and females have similar investments in leaf nitrogen but females exhibit greater incorporation of nitrogen from nitrogen-fixing soil bacteria, (ii) males display higher photosynthetic capacity but females show decreased [corrected] water-use efficiency, and (iii) males have higher stomatal conductance but similar stomatal density to females. This study is the first to examine the ecophysiological differences that have evolved in the oldest dioecious lineage of seed-bearing plants. Our results show unexpected differences in photosynthetic physiology and highlight the co-evolution with nitrogen fixing soil bacteria as a potential new key player in an old lineage.

Exposing the illegal trade in cycad species (Cycadophyta:Encephalartos) at two traditional medicine markets in South Africa using DNA barcoding

Species in the cycad genus Encephalartos are listed in CITES Appendix I and as Threatened or Protected Species in terms of South Africa’s National Environmental Management: Biodiversity Act (NEM:BA) of 2004. Despite regulations, illegal plant harvesting for medicinal trade has continued in South Africa and resulted in declines in cycad populations and even complete loss of sub-populations. Encephalartos is traded at traditional medicine markets in South Africa in the form of bark strips and stem sections; thus, determining the species traded presents a major challenge due to a lack of characteristic plant parts. Here, a case study is presented on the use of DNA barcoding to identify cycads sold at the Faraday and Warwick traditional medicine markets in Johannesburg and Durban, respectively. Market samples were sequenced for the core DNA barcodes (rbcLa and matK) as well as two additional regions: nrITS and trnH-psbA. The barcoding database for cycads at the University of Johannesburg was utilized to assign query samples to known species. Three approaches were followed: tree-based, similarity-based, and character-based (BRONX) methods. Market samples identified were Encephalartos ferox (Near Threatened), Encephalartos lebomboensis (Endangered), Encephalartos natalensis (Near Threatened), Encephalartos senticosus (Vulnerable), and Encephalartos villosus (Least Concern). Results from this study are crucial for making appropriate assessments and decisions on how to manage these markets.

Evolutionary Stasis in Cycad Plastomes and the First Case of Plastome GC-Biased Gene Conversion

In angiosperms, gene conversion has been known to reduce the mutational load of plastid genomes (the plastomes). Particularly, more frequent gene conversions in inverted repeat (IR) than in single copy (SC) regions result in contrasting substitution rates between these two regions. However, little has been known about the effect of gene conversion in the evolution of gymnosperm plastomes. Cycads (Cycadophyta) are the second largest gymnosperm group. Evolutionary study of their plastomes is limited to the basal cycad genus, Cycas. In this study, we addressed three questions. 1) Do the plastomes of other cycad genera evolve slowly as previously observed in the plastome of Cycas taitungensis? 2) Do substitution rates differ between their SC and IR regions? And 3) Does gene conversion occur in the cycad plastomes? If yes, is it AT-biased or GC-biased? Plastomes of eight species from other eight genera of cycads were sequenced. These plastomes are highly conserved in genome organization. Excluding ginkgo, cycad plastomes have significantly lower synonymous and nonsynonymous substitution rates than other gymnosperms, reflecting their evolutionary stasis in nucleotide mutations. In the IRs of cycad plastomes, the reduced substitution rates and GC-biased mutations are associated with a GC-biased gene conversion (gBGC) mechanism. Further investigations suggest that in cycads, gBGC is able to rectify plastome-wide mutations. Therefore, this study is the first to uncover the plastomic gBGC in seed plants. We also propose a gBGC model to interpret the dissimilar evolutionary patterns as well as the compositionally biased mutations in the SC and IR regions of cycad plastomes.

Complete Chloroplast Genome of the Wollemi Pine (Wollemia nobilis): Structure and Evolution

The Wollemi pine (Wollemia nobilis) is a rare Southern conifer with striking morphological similarity to fossil pines. A small population of W. nobilis was discovered in 1994 in a remote canyon system in the Wollemi National Park (near Sydney, Australia). This population contains fewer than 100 individuals and is critically endangered. Previous genetic studies of the Wollemi pine have investigated its evolutionary relationship with other pines in the family Araucariaceae, and have suggested that the Wollemi pine genome contains little or no variation. However, these studies were performed prior to the widespread use of genome sequencing, and their conclusions were based on a limited fraction of the Wollemi pine genome. In this study, we address this problem by determining the entire sequence of the W. nobilis chloroplast genome. A detailed analysis of the structure of the genome is presented, and the evolution of the genome is inferred by comparison with the chloroplast sequences of other members of the Araucariaceae and the related family Podocarpaceae. Pairwise alignments of whole genome sequences, and the presence of unique pseudogenes, gene duplications and insertions in W. nobilis and Araucariaceae, indicate that the W. nobilis chloroplast genome is most similar to that of its sister taxon Agathis. However, the W. nobilis genome contains an unusually high number of repetitive sequences, and these could be used in future studies to investigate and conserve any remnant genetic diversity in the Wollemi pine.

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

Background

Bayesian relaxed-clock dating has significantly influenced our understanding of the timeline of biotic evolution. This approach requires the use of priors on the branching process, yet little is known about their impact on divergence time estimates. We investigated the effect of branching priors using the iconic cycads. We conducted phylogenetic estimations for 237 cycad species using three genes and two calibration strategies incorporating up to six fossil constraints to (i) test the impact of two different branching process priors on age estimates, (ii) assess which branching prior better fits the data, (iii) investigate branching prior impacts on diversification analyses, and (iv) provide insights into the diversification history of cycads.

Results

Using Bayes factors, we compared divergence time estimates and the inferred dynamics of diversification when using Yule versus birth-death priors. Bayes factors were calculated with marginal likelihood estimated with stepping-stone sampling. We found striking differences in age estimates and diversification dynamics depending on prior choice. Dating with the Yule prior suggested that extant cycad genera diversified in the Paleogene and with two diversification rate shifts. In contrast, dating with the birth-death prior yielded Neogene diversifications, and four rate shifts, one for each of the four richest genera. Nonetheless, dating with the two priors provided similar age estimates for the divergence of cycads from Ginkgo (Carboniferous) and their crown age (Permian). Of these, Bayes factors clearly supported the birth-death prior.

Conclusions

These results suggest the choice of the branching process prior can have a drastic influence on our understanding of evolutionary radiations. Therefore, all dating analyses must involve a model selection process using Bayes factors to select between a Yule or birth-death prior, in particular on ancient clades with a potential pattern of high extinction. We also provide new insights into the history of cycad diversification because we found (i) periods of extinction along the long branches of the genera consistent with fossil data, and (ii) high diversification rates within the Miocene genus radiations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0347-8) contains supplementary material, which is available to authorized users.

Extending the generality of leaf economic design principles in the cycads, an ancient lineage

Cycads are the most ancient lineage of living seed plants, but the design of their leaves has received little study. We tested whether cycad leaves are governed by the same fundamental design principles previously established for ferns, conifers and angiosperms, and characterized the uniqueness of this relict lineage in foliar trait relationships. Leaf structure, photosynthesis, hydraulics and nutrient composition were studied in 33 cycad species from nine genera and three families growing in two botanical gardens. Cycads varied greatly in leaf structure and physiology. Similarly to other lineages, light-saturated photosynthetic rate per mass (Am ) was related negatively to leaf mass per area and positively to foliar concentrations of chlorophyll, nitrogen (N), phosphorus and iron, but unlike angiosperms, leaf photosynthetic rate was not associated with leaf hydraulic conductance. Cycads had lower photosynthetic N use efficiency and higher photosynthetic performance relative to hydraulic capacity compared with other lineages. These findings extend the relationships shown for foliar traits in angiosperms to the cycads. This functional convergence supports the modern synthetic understanding of leaf design, with common constraints operating across lineages, even as they highlight exceptional aspects of the biology of this key relict lineage.

The historical demography and genetic variation of the endangered Cycas multipinnata (Cycadaceae) in the red river region, examined by chloroplast DNA sequences and microsatellite markers

Cycas multipinnata C.J. Chen & S.Y. Yang is a cycad endemic to the Red River drainage region that occurs under evergreen forest on steep limestone slopes in Southwest China and northern Vietnam. It is listed as endangered due to habitat loss and over-collecting for the ornamental plant trade, and only several populations remain. In this study, we assess the genetic variation, population structure, and phylogeography of C. multipinnata populations to help develop strategies for the conservation of the species. 60 individuals from six populations were used for chloroplast DNA (cpDNA) sequencing and 100 individuals from five populations were genotyped using 17 nuclear microsatellites. High genetic differentiation among populations was detected, suggesting that pollen or seed dispersal was restricted within populations. Two main genetic clusters were observed in both the cpDNA and microsatellite loci, corresponding to Yunnan China and northern Vietnam. These clusters indicated low levels of gene flow between the regions since their divergence in the late Pleistocene, which was inferred from both Bayesian and coalescent analysis. In addition, the result of a Bayesian skyline plot based on cpDNA portrayed a long history of constant population size followed by a decline in the last 50,000 years of C. multipinnata that was perhaps affected by the Quaternary glaciations, a finding that was also supported by the Garza-Williamson index calculated from the microsatellite data. The genetic consequences produced by climatic oscillations and anthropogenic disturbances are considered key pressures on C. multipinnata. To establish a conservation management plan, each population of C. multipinnata should be recognized as a Management Unit (MU). In situ and ex situ actions, such as controlling overexploitation and creating a germplasm bank with high genetic diversity, should be urgently implemented to preserve this species.

Nuclear ribosomal ITS functional paralogs resolve the phylogenetic relationships of a late-Miocene radiation cycad Cycas (Cycadaceae)

Cycas is the most widespread and diverse genus among the ancient cycads, but the extant species could be the product of late Miocene rapid radiations. Taxonomic treatments to date for this genus are quite controversial, which makes it difficult to elucidate its evolutionary history. We cloned 161 genomic ITS sequences from 31 species representing all sections of Cycas. The divergent ITS paralogs were examined within each species and identified as putative pseudogenes, recombinants and functional paralogs. Functional paralogs were used to reconstruct phylogenetic relationships with pseudogene sequences as molecular outgroups, since an unambiguous ITS sequence alignment with their closest relatives, the Zamiaceae, is unachievable. A fully resolved and highly supported tree topology was obtained at the section level, with two major clades including six minor clades. The results fully supported the classification scheme proposed by Hill (2004) at the section level, with the minor clades representing his six sections. The two major clades could be recognised as two subgenera. The obtained pattern of phylogenetic relationships, combined with the different seed dispersal capabilities and paleogeography, allowed us to propose a late Miocene rapid radiation of Cycas that might have been promoted by vicariant events associated with the complex topography and orogeny of South China and adjacent regions. In contrast, transoceanic dispersals might have played an important role in the rapid diversification of sect. Cycas, whose members have evolved a spongy layer in their seeds aiding water dispersals.

Phylotranscriptomics: saturated third codon positions radically influence the estimation of trees based on next-gen data

Recent advancements in molecular sequencing techniques have led to a surge in the number of phylogenetic studies that incorporate large amounts of genetic data. We test the assumption that analyzing large number of genes will lead to improvements in tree resolution and branch support using moths in the superfamily Bombycoidea, a group with some interfamilial relationships that have been difficult to resolve. Specifically, we use a next-gen data set that included 19 taxa and 938 genes (∼1.2M bp) to examine how codon position and saturation might influence resolution and node support among three key families. Maximum likelihood, parsimony, and species tree analysis using gene tree parsimony, on different nucleotide and amino acid data sets, resulted in largely congruent topologies with high bootstrap support compared with prior studies that included fewer loci. However, for a few shallow nodes, nucleotide and amino acid data provided high support for conflicting relationships. The third codon position was saturated and phylogenetic analysis of this position alone supported a completely different, potentially misleading sister group relationship. We used the program RADICAL to assess the number of genes needed to fix some of these difficult nodes. One such node originally needed a total of 850 genes but only required 250 when synonymous signal was removed. Our study shows that, in order to effectively use next-gen data to correctly resolve difficult phylogenetic relationships, it is necessary to assess the effects of synonymous substitutions and third codon positions.

Phylogeny and divergence times of gymnosperms inferred from single-copy nuclear genes

Phylogenetic reconstruction is fundamental to study evolutionary biology and historical biogeography. However, there was not a molecular phylogeny of gymnosperms represented by extensive sampling at the genus level, and most published phylogenies of this group were constructed based on cytoplasmic DNA markers and/or the multi-copy nuclear ribosomal DNA. In this study, we use LFY and NLY, two single-copy nuclear genes that originated from an ancient gene duplication in the ancestor of seed plants, to reconstruct the phylogeny and estimate divergence times of gymnosperms based on a complete sampling of extant genera. The results indicate that the combined LFY and NLY coding sequences can resolve interfamilial relationships of gymnosperms and intergeneric relationships of most families. Moreover, the addition of intron sequences can improve the resolution in Podocarpaceae but not in cycads, although divergence times of the cycad genera are similar to or longer than those of the Podocarpaceae genera. Our study strongly supports cycads as the basal-most lineage of gymnosperms rather than sister to Ginkgoaceae, and a sister relationship between Podocarpaceae and Araucariaceae and between Cephalotaxaceae-Taxaceae and Cupressaceae. In addition, intergeneric relationships of some families that were controversial, and the relationships between Taxaceae and Cephalotaxaceae and between conifers and Gnetales are discussed based on the nuclear gene evidence. The molecular dating analysis suggests that drastic extinctions occurred in the early evolution of gymnosperms, and extant coniferous genera in the Northern Hemisphere are older than those in the Southern Hemisphere on average. This study provides an evolutionary framework for future studies on gymnosperms.

An ancient satellite DNA has maintained repetitive units of the original structure in most species of the living fossil plant genus Zamia

ZpS1 satellite DNA is specific to the genus Zamia and presents repetitive units organized as long arrays and also as very short arrays dispersed in the genome. We have characterized the structure of the ZpS1 repeats in 12 species representative of the whole geographic distribution of the genus. In most species, the clone most common sequences (cMCS) were so similar that a general most common sequence (GMCS) of the ZpS1 repetitive unit in the genus could be obtained. The few partial variations from the GMCS found in cMCS of some species correspond to variable positions present in most other species, as indicated by the clone consensus sequences (cCS). Two species have an additional species-specific variety of ZpS1 satellite. The dispersed repeats were found to contain more mutations than repeats from long arrays. Our results indicate that all or most species of Zamia inherited the ZpS1 satellite from a common ancestor in Miocene and have maintained repetitive units of the original structure till present. The features of ZpS1 satellite in the genus Zamia are poorly compatible with the model of concerted evolution, but they are perfectly consistent with a new model of satellite evolution based on experimental evidences indicating that a specific amplification-substitution repair mechanism maintains the homogeneity and stability of the repeats structure in each satellite DNA originally present in a species as long as the species exists.

Conserved genetic regions across angiosperms as tools to develop single-copy nuclear markers in gymnosperms: an example using cycads

Several individuals of the Caribbean Zamia clade and other cycad genera were used to identify single-copy nuclear genes for phylogeographic and phylogenetic studies in Cycadales. Two strategies were employed to select target loci: (i) a tblastX search of Arabidopsis conserved ortholog sequence (COS) set and (ii) a tblastX search of Arabidopsis-Populus-Vitis-Oryza Shared Single-Copy genes (APVO SSC) against the EST Zamia databases in GenBank. From the first strategy, 30 loci were selected, and from the second, 16 loci. In both cases, the matching GenBank accessions of Zamia were used as a query for retrieving highly similar sequences from Cycas, Picea, Pinus species or Ginkgo biloba. After retrieving and aligning all the sequences in each locus, intron predictions were completed to assist in primer design. PCR was carried out in three rounds to detect paralogous loci. A total of 29 loci were successfully amplified as a single band of which 20 were likely single-copy loci. These loci showed different diversity and divergence levels. A preliminary screening allowed us to select 8 promising loci (40S, ATG2, BG, GroES, GTP, LiSH, PEX4 and TR) for the Zamia pumila complex and 4 loci (COS26, GroES, GTP and HTS) for all other cycad genera.

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

BACKGROUND AND AIMS

Despite a recent new classification, a stable phylogeny for the cycads has been elusive, particularly regarding resolution of Bowenia, Stangeria and Dioon. In this study, five single-copy nuclear genes (SCNGs) are applied to the phylogeny of the order Cycadales. The specific aim is to evaluate several gene tree-species tree reconciliation approaches for developing an accurate phylogeny of the order, to contrast them with concatenated parsimony analysis and to resolve the erstwhile problematic phylogenetic position of these three genera.

METHODS

DNA sequences of five SCNGs were obtained for 20 cycad species representing all ten genera of Cycadales. These were analysed with parsimony, maximum likelihood (ML) and three Bayesian methods of gene tree-species tree reconciliation, using Cycas as the outgroup. A calibrated date estimation was developed with Bayesian methods, and biogeographic analysis was also conducted.

KEY RESULTS

Concatenated parsimony, ML and three species tree inference methods resolve exactly the same tree topology with high support at most nodes. Dioon and Bowenia are the first and second branches of Cycadales after Cycas, respectively, followed by an encephalartoid clade (Macrozamia-Lepidozamia-Encephalartos), which is sister to a zamioid clade, of which Ceratozamia is the first branch, and in which Stangeria is sister to Microcycas and Zamia.

CONCLUSIONS

A single, well-supported phylogenetic hypothesis of the generic relationships of the Cycadales is presented. However, massive extinction events inferred from the fossil record that eliminated broader ancestral distributions within Zamiaceae compromise accurate optimization of ancestral biogeographical areas for that hypothesis. While major lineages of Cycadales are ancient, crown ages of all modern genera are no older than 12 million years, supporting a recent hypothesis of mostly Miocene radiations. This phylogeny can contribute to an accurate infrafamilial classification of Zamiaceae.

Chloroplast phylogenomics indicates that Ginkgo biloba is sister to cycads

Molecular phylogenetic studies have not yet reached a consensus on the placement of Ginkgoales, which is represented by the only living species, Ginkgo biloba (common name: ginkgo). At least six discrepant placements of ginkgo have been proposed. This study aimed to use the chloroplast phylogenomic approach to examine possible factors that lead to such disagreeing placements. We found the sequence types used in the analyses as the most critical factor in the conflicting placements of ginkgo. In addition, the placement of ginkgo varied in the trees inferred from nucleotide (NU) sequences, which notably depended on breadth of taxon sampling, tree-building methods, codon positions, positions of Gnetopsida (common name: gnetophytes), and including or excluding gnetophytes in data sets. In contrast, the trees inferred from amino acid (AA) sequences congruently supported the monophyly of a ginkgo and Cycadales (common name: cycads) clade, regardless of which factors were examined. Our site-stripping analysis further revealed that the high substitution saturation of NU sequences mainly derived from the third codon positions and contributed to the variable placements of ginkgo. In summary, the factors we surveyed did not affect results inferred from analyses of AA sequences. Congruent topologies in our AA trees give more confidence in supporting the ginkgo-cycad sister-group hypothesis.

Thrips pollination of Mesozoic gymnosperms

Within modern gymnosperms, conifers and Ginkgo are exclusively wind pollinated whereas many gnetaleans and cycads are insect pollinated. For cycads, thrips are specialized pollinators. We report such a specialized pollination mode from Early Cretaceous amber of Spain, wherein four female thrips representing a genus and two species in the family Melanthripidae were covered by abundant Cycadopites pollen grains. These females bear unique ring setae interpreted as specialized structures for pollen grain collection, functionally equivalent to the hook-tipped sensilla and plumose setae on the bodies of bees. The most parsimonious explanation for this structure is parental food provisioning for larvae, indicating subsociality. This association provides direct evidence of specialized collection and transportation of pollen grains and likely gymnosperm pollination by 110-105 million years ago, possibly considerably earlier.

Genome sizes for all genera of Cycadales

Nuclear DNA content (2C) is reported for all genera of the Cycadales, using flow cytometry with propidium iodide. Nuclear DNA content ranges from 24 to 64 pg in cycads. This implies that the largest genome contains roughly 40 × 10(9) more base pairs than the smallest genome. The narrow range in nuclear DNA content within a genus is remarkable for such an old group. Furthermore, 42 of the 58 plants measured, covering five genera, have 18 chromosomes. They vary from 36.1 to 64.7 pg, covering the whole range of genome sizes (excluding the genome of Cycas). Hence, their does not seem to be a correlation between genome size and the number of chromosomes.

Defining the limits of flowers: the challenge of distinguishing between the evolutionary products of simple versus compound strobili

Recent phylogenetic reconstructions suggest that axially condensed flower-like structures evolved iteratively in seed plants from either simple or compound strobili. The simple-strobilus model of flower evolution, widely applied to the angiosperm flower, interprets the inflorescence as a compound strobilus. The conifer cone and the gnetalean 'flower' are commonly interpreted as having evolved from a compound strobilus by extreme condensation and (at least in the case of male conifer cones) elimination of some structures present in the presumed ancestral compound strobilus. These two hypotheses have profoundly different implications for reconstructing the evolution of developmental genetic mechanisms in seed plants. If different flower-like structures evolved independently, there should intuitively be little commonality of patterning genes. However, reproductive units of some early-divergent angiosperms, including the extant genus Trithuria (Hydatellaceae) and the extinct genus Archaefructus (Archaefructaceae), apparently combine features considered typical of flowers and inflorescences. We re-evaluate several disparate strands of comparative data to explore whether flower-like structures could have arisen by co-option of flower-expressed patterning genes into independently evolved condensed inflorescences, or vice versa. We discuss the evolution of the inflorescence in both gymnosperms and angiosperms, emphasising the roles of heterotopy in dictating gender expression and heterochrony in permitting internodal compression.

Molecular evolution and positive Darwinian selection of the chloroplast maturase matK

It is not clear whether matK evolves under Darwinian selection. In this study, 70 plant groups, representing 2,279 species at various evolutionary levels, were used to illustrate the molecular adaptation and evolutionary dynamics of gene divergence in matKs. Selective influences were investigated using standard dN/dS ratio methods. Analyses revealed the presence of positive selection in matKs of 32 plant groups. More positively selected sites were detected in the N-terminal region than in the RT domain and domain X of matK. Moreover, removing amino acid sites that are under positive selection has a significant effect on the bootstrap values of phylogenetic reconstruction. Our results suggest that the rapidly evolving matK evolves under positive selection in some lineages of land plants. Several regions of matK have experienced molecular adaptation, which fine-tunes maturase performance.

Beetle pollination of the fruit-scented cones of the South African cycad Stangeria eriopus

There has been considerable uncertainty about the importance of wind vs. insects in cycad pollination, but recent studies in several cycad genera have indicated that these are pollinated primarily, if not exclusively, by insects. Stangeria represents an isolated southern African cycad lineage previously thought to be wind-pollinated. Unlike in most other cycads, there is no evidence of cone thermogenesis in Stangeria. We found that the scent of both male and female Stangeria cones mimics that of fermented fruit, the main volatiles being esters of acetic acid, ketones, and aldehydes. We found a large variety of insect visitors on the cones, the most common ones being sap and rove beetles (Coleoptera: Nitidulidae, Staphylinidae) and fruit flies (Diptera: Drosophilidae). Of these, only sap beetles (Nitidulidae) were able to effect pollination under experimental conditions. Because sap beetles are also pollinators of Cycas and members of several ancient angiosperm families, their role in the pollination of Stangeria adds interesting details to the role this group of insects has played in the history of plant-pollinator interactions.

Phylogenetic signal in matK vs. trnK: a case study in early diverging eudicots (angiosperms)

The matK gene has been among the most useful loci for resolving plant phylogenetic relationships at different evolutionary time-scales, but much less is known about the phylogenetic utility of the flanking trnK intron, especially for deep level phylogenetics. We compared the relative performance of matK and trnK intron regions for resolving the relationships of the early diverging eudicots (angiosperms). The two regions display similar nucleotide compositions and distributions of rate variation among sites. The trnK intron sequences also provide similar levels of phylogenetic information per-site as matK. Combining the trnK intron sequences with matK increases overall bootstrap support for the early diverging eudicots compared to analyses of matK alone. MP, ML and Bayesian analyses provide strong support for eudicots, the sister group relationship of Ranunculales to remaining eudicots, and a Buxales+Trochodendraceae+core eudicots clade. matK and the trnK intron support conflicting positions for Buxales and Trochodendrales in relation to the core eudicots.

DNA barcoding in the cycadales: testing the potential of proposed barcoding markers for species identification of cycads

Barcodes are short segments of DNA that can be used to uniquely identify an unknown specimen to species, particularly when diagnostic morphological features are absent. These sequences could offer a new forensic tool in plant and animal conservation-especially for endangered species such as members of the Cycadales. Ideally, barcodes could be used to positively identify illegally obtained material even in cases where diagnostic features have been purposefully removed or to release confiscated organisms into the proper breeding population. In order to be useful, a DNA barcode sequence must not only easily PCR amplify with universal or near-universal reaction conditions and primers, but also contain enough variation to generate unique identifiers at either the species or population levels. Chloroplast regions suggested by the Plant Working Group of the Consortium for the Barcode of Life (CBoL), and two alternatives, the chloroplast psbA-trnH intergenic spacer and the nuclear ribosomal internal transcribed spacer (nrITS), were tested for their utility in generating unique identifiers for members of the Cycadales. Ease of amplification and sequence generation with universal primers and reaction conditions was determined for each of the seven proposed markers. While none of the proposed markers provided unique identifiers for all species tested, nrITS showed the most promise in terms of variability, although sequencing difficulties remain a drawback. We suggest a workflow for DNA barcoding, including database generation and management, which will ultimately be necessary if we are to succeed in establishing a universal DNA barcode for plants.

Developmental bases for key innovations in the seed-plant microgametophyte

The early evolution of seed plants was marked by a sequence of putative key innovations associated with the microgametophyte, including reversal of microspore polarity, loss of sperm motility and co-option of the pollen tube to a new role, siphonogamy. Data from several interfacing fields of plant development are reviewed here in the context of recent palaeobotanical discoveries and improved understanding of seed-plant relationships. The proximal-distal microspore polarity transition that occurred early in seed-plant evolution represents a fossil fingerprint for an underlying series of radical developmental shifts involving quadripartite partitioning at meiosis and an asymmetric primary mitosis. It had important downstream effects, not only on aperture location and site of germination but also on microgametophyte polarity, and, perhaps, indirectly on sperm motility.

Molecular phylogenetics and asexuality in the brine shrimp Artemia

Explaining cases of long-term persistence of parthenogenesis has proven an arduous task for evolutionary biologists. Interpreting sexual-asexual interactions though has recently advanced owing to methodological design, increased taxon sampling and choice of model organisms. We inferred the phylogeny of Artemia, a halophilic branchiopod genus of sexual and parthenogenetic forms with cosmopolitan distribution, marked geographic patterns and ecological partitioning. Joint analysis of newly derived ITS1 sequences and 16S RFLP markers from global isolates indicates significant interspecific divergence as well as pronounced diversity for parthenogens, matching that of sexual ancestors. Maximum parsimony, maximum likelihood, and Bayesian methods were largely congruent in reconstructing the phylogeny of the genus. Given the current sampling, at least four independent origins of parthenogenesis are deduced. Molecular clock calibrations based on biogeographic landmarks indicate that the lineage leading to A. persimilis diverged from the common ancestor of all Artemia species between 80 and 90 MYA at the time of separation of Africa from South America, whereas parthenogenesis first appeared at least 3 MYA. Common mitochondrial DNA haplotypes delineate A. urmiana and A. tibetiana as possible maternal parents of several clonal lineages. A novel topological placement of A. franciscana as a sister clade to all Asian Artemia and parthenogenetic forms is proposed and also supported by ITS1 length and other existing data.

Cataphylls of the Middle Triassic cycad Antarcticycas schopfii and new insights into cycad evolution

Cataphylls associated with the Middle Triassic stem genus Antarcticycas are described, and their impact on understanding cycad evolution is discussed. The cataphylls of Antarcticycas are triangular in outline and flattened adaxially with lateral flanges. The outer surfaces are covered with a ramentum of filamentous hairs, the epidermis is a single cell layer thick, and the ground tissue is parenchymatous with mucilage canals and sclereids. Vascular bundles form a distinct inverted omega-shaped pattern characteristic of the Cycadales observed in petioles of extant species. The structures in Antarcticycas are interpreted as cataphylls based on overall morphology, presence of straight vascular strands in the cortex of the associated stem, and lack of fascicular cambia in the vascular bundles. Because much of the overall diversity of Cycadales is represented by fossils, integrating fossil taxa into explicit phylogenetic hypotheses is important for understanding cycad evolution. Therefore, character and minimum age mapping were performed on a phylogeny of extant and fossil taxa including Antarcticycas. The results suggest that major extant lineages of Cycadales had diverged by the Permian to Triassic and that certain synapomorphies for Cycadales had evolved by the Permian. Evidence of insect feeding on Antarcticycas suggests that associations between cycads and insects are ancient.

The chloroplast genome of Phalaenopsis aphrodite (Orchidaceae): comparative analysis of evolutionary rate with that of grasses and its phylogenetic implications

Whether the Amborella/Amborella-Nymphaeales or the grass lineage diverged first within the angiosperms has recently been debated. Central to this issue has been focused on the artifacts that might result from sampling only grasses within the monocots. We therefore sequenced the entire chloroplast genome (cpDNA) of Phalaenopsis aphrodite, Taiwan moth orchid. The cpDNA is a circular molecule of 148,964 bp with a comparatively short single-copy region (11,543 bp) due to the unusual loss and truncation/scattered deletion of certain ndh subunits. An open reading frame, orf91, located in the complementary strand of the rrn23 was reported for the first time. A comparison of nucleotide substitutions between P. aphrodite and the grasses indicates that only the plastid expression genes have a strong positive correlation between nonsynonymous (Ka) and synonymous (Ks) substitutions per site, providing evidence for a generation time effect, mainly across these genes. Among the intron-containing protein-coding genes of the sampled monocots, the Ks of the genes are significantly correlated to transitional substitutions of their introns. We compiled a concatenated 61 protein-coding gene alignment for the available 20 cpDNAs of vascular plants and analyzed the data set using Bayesian inference, maximum parsimony, and neighbor-joining (NJ) methods. The analyses yielded robust support for the Amborella/Amborella-Nymphaeales-basal hypothesis and for the orchid and grasses together being a monophyletic group nested within the remaining angiosperms. However, the NJ analysis using Ka, the first two codon positions, or amino acid sequences, respectively, supports the monocots-basal hypothesis. We demonstrated that these conflicting angiosperm phylogenies are most probably linked to the transitional sites at all codon positions, especially at the third one where the strong base-composition bias and saturation effect take place.