Eco-evolutionary evidence for the global diversity pattern of Cycas (Cycadaceae).

The evolution of the latitudinal diversity gradient (LDG), characterized by a peak in diversity toward the tropics, has captured significant attention in evolutionary biology and ecology. However, the inverse LDG (i-LDG) mechanism, wherein species richness increases toward the poles, remains inadequately explored. Cycads are among one of the oldest lineages of extant seed plants and have undergone extensive diversification in the tropics. Intriguingly, the extant cycad abundance exhibits an i-LDG pattern, and the underlying causes for this phenomenon remain largely elusive. Here, using 1,843 nuclear genes from a nearly complete sampling, we conducted comprehensive phylogenomic analyses to establish a robust species-level phylogeny for Cycas, the largest genus within cycads. We then reconstructed the spatial-temporal dynamics and integrated global environmental data to evaluate the roles of species ages, diversification rates, contemporary environment, and conservatism to ancestral niches in shaping the i-LDG pattern. We found Cycas experienced decreased diversification rates, coupled with the cooling temperature since its origin in the Eocene from continental Asia. Different regions have distinctively contributed to the formation of i-LDG for Cycas, with the northern hemisphere acting as evolutionary museums and the southern hemisphere serving as cradles. Moreover, water-related climate variables, specifically precipitation seasonality and potential evapotranspiration, were identified as paramount factors constraining Cycas species richness in the rainforest biome near the equator. Notably, the adherence to ancestral monsoonal climates emerges as a critical factor in sustaining the diversity pattern. This study underscores the imperative of integrating both evolutionary and ecological approaches to comprehensively unravel the mechanisms underpinning global biodiversity patterns.

Genomes tell tales of spores versus seeds.

Genomes from ferns and a cycad reveal deep roots of plant reproduction.

The Cycas genome and the early evolution of seed plants.

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

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

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

The morpho-genetic and ecological niche analyses reveal the existence of climatically restricted Cycas zeylanica complex in Sri Lanka.

Taxonomy and phylogenesis of Sri Lankan cycad species of the subsection Rumphiae has not been fully resolved and therefore, we conducted an island-wide survey of cycads of the subsection to assess their morphological or genetic variations while exploring the phylogenetic relationship between Sri Lankan Rumphiae and other world cycad species. Further, we assessed the possible distribution of the species in the region through climatic profiling, using maximum entropy modeling approach. We analyzed 21 variable morphological features in collected specimens and used the polymorphism of trnH-psbA locus to understand the phylogeny. The distance tree drawn from the principal component analysis revealed a significant variation in female reproductive structures. The maximum likelihood tree separated Sri Lankan Cycas zeylanica to a well-supported unigeneric clade (bootstrap = 96, posterior probability = 100) with shallow divergence. Ecological niche modeling supported the existence of Cycas zeylanica in South East Asia and in southern Western Ghats in India in addition to the Wet Zone of Sri Lanka. We rename the taxa as Cycas zeylanica complex based on the observed high morphological diversity of female reproductive structures which might have ascended due to multiple introductions of South East Asian cycads by long distance dispersal of seeds through sea currents.

RADseq as a valuable tool for plants with large genomes-A case study in cycads.

Full genome sequencing of organisms with large and complex genomes is intractable and cost ineffective under most research budgets. Cycads (Cycadales) represent one of the oldest lineages of the extant seed plants and, partly due to their age, have incredibly large genomes up to ~60 Gbp. Restriction site-associated DNA sequencing (RADseq) offers an approach to find genome-wide informative markers and has proven to be effective with both model and nonmodel organisms. We tested the application of RADseq using ezRAD across all 10 genera of the Cycadales including an example data set of Cycas calcicola representing 72 samples from natural populations. Using previously available plastid and mitochondrial genomes as references, reads were mapped recovering plastid and mitochondrial genome regions and nuclear markers for all of the genera. De novo assembly generated up to 138,407 high-depth clusters and up to 1,705 phylogenetically informative loci for the genera, and 4,421 loci for the example assembly of C. calcicola. The number of loci recovered by de novo assembly was lower than previous RADseq studies, yet still sufficient for downstream analysis. However, the number of markers could be increased by relaxing our assembly parameters, especially for the C. calcicola data set. Our results demonstrate the successful application of RADseq across the Cycadales to generate a large number of markers for all genomic compartments, despite the large number of plastids present in a typical plant cell. Our modified protocol was adapted to be applied to cycads and other organisms with large genomes to yield many informative genome-wide markers.

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.

Middle-Upper Pleistocene climate changes shaped the divergence and demography of Cycas guizhouensis (Cycadaceae): Evidence from DNA sequences and microsatellite markers.

Climatic oscillations in the Pleistocene have had profound effects on the demography and genetic diversity of many extant species. Cycas guizhouensis Lan &R.F. Zou is an endemic and endangered species in Southwest China that is primarily distributed along the valleys of the Nanpan River. In this study, we used four chloroplast DNAs (cpDNA), three nuclear genes (nDNA) and 13 microsatellite (SSR) loci to investigate the genetic structure, divergence time and demographic history of 11 populations of C. guizhouensis. High genetic diversity and high levels of genetic differentiation among the populations were observed. Two evolutionary units were revealed based on network and Structure analysis. The divergence time estimations suggested that haplotypes of C. guizhouensis were diverged during the Middle-Upper Pleistocene. Additionally, the demographic histories deduced from different DNA sequences were discordant, but overall indicated that C. guizhouensis had experienced a recent population expansion during the post-glacial period. Microsatellite data revealed that there was a contraction in effective population size in the past. These genetic features allow conservation measures to be taken to ensure the protection of this endangered species from extinction.

Selection and Validation of Reference Genes for qRT-PCR in Cycas elongata.

Quantitative reverse transcription PCR (qRT-PCR) is a sensitive technique used in gene expression studies. To achieve a reliable quantification of transcripts, appropriate reference genes are required for comparison of transcripts in different samples. However, few reference genes are available for non-model taxa, and to date, reliable reference genes in Cycas elongata have not been well characterized. In this study, 13 reference genes (ACT7, TUB, UBQ, EIF4, EF1, CLATHRIN1, PP2A, RPB2, GAPC2, TIP41, MAPK, SAMDC and CYP) were chosen from the transcriptome database of C. elongata, and these genes were evaluated in 8 different organ samples. Three software programs, NormFinder, GeNorm and BestKeeper, were used to validate the stability of the potential reference genes. Results obtained from these three programs suggested that CeGAPC2 and CeRPB2 are the most stable reference genes, while CeACT7 is the least stable one among the 13 tested genes. Further confirmation of the identified reference genes was established by the relative expression of AGAMOUSE gene of C. elongata (CeAG). While our stable reference genes generated consistent expression patterns in eight tissues, we note that our results indicate that an inappropriate reference gene might cause erroneous results. Our systematic analysis for stable reference genes of C. elongata facilitates further gene expression studies and functional analyses of this species.

Astonishing 35S rDNA diversity in the gymnosperm species Cycas revoluta Thunb.

In all eukaryotes, the highly repeated 35S ribosomal DNA (rDNA) sequences encoding 18S-5.8S-26S ribosomal RNA (rRNA) typically show high levels of intragenomic uniformity due to homogenisation processes, leading to concerted evolution of 35S rDNA repeats. Here, we compared 35S rDNA divergence in several seed plants using next generation sequencing and a range of molecular and cytogenetic approaches. Most species showed similar 35S rDNA homogeneity indicating concerted evolution. However, Cycas revoluta exhibits an extraordinary diversity of rDNA repeats (nucleotide sequence divergence of different copies averaging 12 %), influencing both the coding and non-coding rDNA regions nearly equally. In contrast, its rRNA transcriptome was highly homogeneous suggesting that only a minority of genes (<20 %) encode functional rRNA. The most common SNPs were C > T substitutions located in symmetrical CG and CHG contexts which were also highly methylated. Both functional genes and pseudogenes appear to cluster on chromosomes. The extraordinary high levels of 35S rDNA diversity in C. revoluta, and probably other species of cycads, indicate that the frequency of repeat homogenisation has been much lower in this lineage, compared with all other land plant lineages studied. This has led to the accumulation of methylation-driven mutations and pseudogenisation. Potentially, the reduced homology between paralogs prevented their elimination by homologous recombination, resulting in long-term retention of rDNA pseudogenes in the genome.

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.

Genetic diversity, genetic structure and demographic history of Cycas simplicipinna (Cycadaceae) assessed by DNA sequences and SSR markers.

BACKGROUND: Cycas simplicipinna (T. Smitinand) K. Hill. (Cycadaceae) is an endangered species in China. There were seven populations and 118 individuals that we could collect were genotyped in this study. Here, we assessed the genetic diversity, genetic structure and demographic history of this species. RESULTS: Analyses of data of DNA sequences (two maternally inherited intergenic spacers of chloroplast, cpDNA and one biparentally inherited internal transcribed spacer region ITS4-ITS5, nrDNA) and sixteen microsatellite loci (SSR) were conducted in the species. Of the 118 samples, 86 individuals from the seven populations were used for DNA sequencing and 115 individuals from six populations were used for the microsatellite study. We found high genetic diversity at the species level, low genetic diversity within each of the seven populations and high genetic differentiation among the populations. There was a clear genetic structure within populations of C. simplicipinna. A demographic history inferred from DNA sequencing data indicates that C. simplicipinna experienced a recent population contraction without retreating to a common refugium during the last glacial period. The results derived from SSR data also showed that C. simplicipinna underwent past effective population contraction, likely during the Pleistocene. CONCLUSIONS: Some genetic features of C. simplicipinna such as having high genetic differentiation among the populations, a clear genetic structure and a recent population contraction could provide guidelines for protecting this endangered species from extinction. Furthermore, the genetic features with population dynamics of the species in our study would help provide insights and guidelines for protecting other endangered species effectively.

Asymmetric introgression in the horticultural living fossil cycas sect. Asiorientales using a genome-wide scanning approach.

The Asian cycads are mostly allopatric, distributed in small population sizes. Hybridization between allopatric species provides clues in determining the mechanism of species divergence. Horticultural introduction provides the chance of interspecific gene flow between allopatric species. Two allopatrically eastern Asian Cycas sect. Asiorientales species, C. revoluta and C. taitungensis, which are widely distributed in Ryukyus and Fujian Province and endemic to Taiwan, respectively, were planted in eastern Taiwan for horticultural reason. Higher degrees of genetic admixture in cultivated samples than wild populations in both cycad species were detected based on multilocus scans by neutral AFLP markers. Furthermore, bidirectional but asymmetric introgression by horticultural introduction of C. revoluta is evidenced by the reanalyses of species associated loci, which are assumed to be diverged after species divergence. Partial loci introgressed from native cycad to the invaders were also detected at the loci of strong species association. Consistent results tested by all neutral loci, and the species-associated loci, specify the recent introgression from the paradox of sharing of ancestral polymorphisms. Phenomenon of introgression of cultivated cycads implies niche conservation among two geographic-isolated cycads, even though the habitats of the extant wild populations of two species are distinct.

RNA editing sites exist in protein-coding genes in the chloroplast genome of Cycas taitungensis.

RNA editing is a post-transcriptional process that results in modifications of ribonucleotides at specific locations. In land plants editing can occur in both mitochondria and chloroplasts and most commonly involves C-to-U changes, especially in seed plants. Using prediction and experimental determination, we investigated RNA editing in 40 protein-coding genes from the chloroplast genome of Cycas taitungensis. A total of 85 editing sites were identified in 25 transcripts. Comparison analysis of the published editotypes of these 25 transcripts in eight species showed that RNA editing events gradually disappear during plant evolution. The editing in the first and third codon position disappeared quicker than that in the second codon position. ndh genes have the highest editing frequency while serine and proline codons were more frequently edited than the codons of other amino acids. These results imply that retained RNA editing sites have imbalanced distribution in genes and most of them may function by changing protein structure or interaction. Mitochondrion protein-coding genes have three times the editing sites compared with chloroplast genes of Cycas, most likely due to slower evolution speed.

Maternal inheritance of plastids and mitochondria in Cycas L. (Cycadaceae).

Cycas is often considered a living fossil, thereby providing a unique model for revealing the evolution of spermatophytes. To date, the genetic inheritance of these archaic plants is not fully understood. The present study seeks to document the process of organelle inheritance in an interspecific cross of Cycas species. Extranuclear organelle DNA from chloroplasts and mitochondria was analyzed using both polymerase chain reaction-restriction fragment length polymorphism analysis and microscopy. Here, we show that the chloroplasts and mitochondria in the progeny of interspecific crosses between Cycas taitungensis and Cycas ferruginea were exclusively inherited from the female parent. Epifluorescence microscopic analyses of the pollen cells from Cycas elongata indicated that there was a significant degradation of organelle DNA in male reproductive cells following maturation; the DNA fluorescent signals were only seen after pollen mitosis two, but not detectable at mature stage. Lack of organelle DNA fluorescent signal in prothallial cells was confirmed by the absence of plastids and mitochondria in electronic microscopic images. In conclusion, these data suggest that the maternal plastid and mitochondrial inheritance in Cycas, native to the old world, are the same as seen in seed plants.

Microsatellite primers in the native perennial cycad Cycas taitungensis (Cycadaceae).

PREMISE OF THE STUDY: Microsatellite primers were developed for the native perennial cycad Cycas taitungensis to evaluate the genetic variation of this endangered insular species. METHODS AND RESULTS: Using a magnetic bead enrichment method and EST data, 16 primer sets were developed and identified for the native Taiwan cycad C. taitungensis. The primers amplified dinucleotide, trinucleotide, and complex repeats with 1-9 alleles per locus. Most primers also amplified DNA from C. revoluta and C. debaoensis. CONCLUSIONS: These results indicate the utility of primers for future studies of the genetic structure of C. taitungensis. In addition, the primers are useful for further phylogeographic studies between C. taitungensis and C. revoluta, which is a closely related species.

Conservation and canalization of gene expression during angiosperm diversification accompany the origin and evolution of the flower.

The origin and rapid diversification of the angiosperms (Darwin's "Abominable Mystery") has engaged generations of researchers. Here, we examine the floral genetic programs of phylogenetically pivotal angiosperms (water lily, avocado, California poppy, and Arabidopsis) and a nonflowering seed plant (a cycad) to obtain insight into the origin and subsequent evolution of the flower. Transcriptional cascades with broadly overlapping spatial domains, resembling the hypothesized ancestral gymnosperm program, are deployed across morphologically intergrading organs in water lily and avocado flowers. In contrast, spatially discrete transcriptional programs in distinct floral organs characterize the more recently derived angiosperm lineages represented by California poppy and Arabidopsis. Deep evolutionary conservation in the genetic programs of putatively homologous floral organs traces to those operating in gymnosperm reproductive cones. Female gymnosperm cones and angiosperm carpels share conserved genetic features, which may be associated with the ovule developmental program common to both organs. However, male gymnosperm cones share genetic features with both perianth (sterile attractive and protective) organs and stamens, supporting the evolutionary origin of the floral perianth from the male genetic program of seed plants.

Editing site analysis in a gymnosperm mitochondrial genome reveals similarities with angiosperm mitochondrial genomes.

Sequence analysis of organelle genomes and comprehensive analysis of C-to-U editing sites from flowering and non-flowering plants have provided extensive sequence information from diverse taxa. This study includes the first comprehensive analysis of RNA editing sites from a gymnosperm mitochondrial genome, and utilizes informatics analyses to determine conserved features in the RNA sequence context around editing sites. We have identified 565 editing sites in 21 full-length and 4 partial cDNAs of the 39 protein-coding genes identified from the mitochondrial genome of Cycas taitungensis. The information profiles and RNA sequence context of C-to-U editing sites in the Cycas genome exhibit similarity in the immediate flanking nucleotides. Relative entropy analyses indicate that similar regions in the 5' flanking 20 nucleotides have information content compared to angiosperm mitochondrial genomes. These results suggest that evolutionary constraints exist on the nucleotide sequences immediately adjacent to C-to-U editing sites, and similar regions are utilized in editing site recognition.

When North and South don't mix: genetic connectivity of a recently endangered oceanic cycad, Cycas micronesica, in Guam using EST-microsatellites.

Subject to environmental changes and recurrent isolation in the last ca. 250 Ma, cycads are often described as relicts of a previously common lineage, with populations characterized by low genetic variation and restricted gene flow. We found that on the island of Guam, the endemic Cycas micronesica has most of the genetic variation of 14 EST-microsatellites distributed within each of 18 genetic populations, from 24 original sampling sites. There were high levels of genetic variation in terms of total number of alleles and private alleles, and moderate levels of inbreeding. Restricted but ongoing gene flow among populations within Guam reveals a genetic mosaic, probably more typical of cycads than previously assumed. Contiguous cycad populations in the north of Guam had higher self-recruitment rates compared to fragmented populations in the south, with no substantial connection between them except for one population. Guam's genetic mosaic may be explained by the influence of forest continuity, seed size, edaphic differences, and human transport of cycads. Also important are the extent of synchrony among flushes of reproductive female seed-bearing sporophylls and restricted pollen movement by an obligate mutualist and generalist insects. An NADH EST-locus under positive selection may reflect pressure from edaphic differences across Guam. This and three other loci are ideal candidates for ecological genomic studies. Given this species' vulnerability due to the recent introduction of the cycad aulacaspis scale, we also identify priority populations for ex situ conservation, and provide a genetic baseline for understanding the effects of invasive species on cycads in the Western Pacific, and islands in general.