Patterns of chloroplast DNA variation in Cycas debaoensis (Cycadaceae): conservation implications

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

Full-Length 16S rRNA and ITS Gene Sequencing Revealed Rich Microbial Flora in Roots of Cycas spp. in China

Cycads have developed a complex root system categorized either as normal or coralloid roots. Past literatures revealed that a great diversity of key microbes is associated with these roots. This recent study aims to comprehensively determine the diversity and community structure of bacteria and fungi associated with the roots of two Cycas spp. endemic to China, Cycas debaoensis Zhong & Chen and Cycas fairylakea D.Y. Wang using high-throughput amplicon sequencing of the full-length 16S rRNA (V1-V9 hypervariable) and short fragment ITS region. The total DNA from 12 root samples were extracted, amplified, sequenced, and analyzed. Resulting sequences were clustered into 61 bacteria and 2128 fungal OTUs. Analysis of community structure revealed that the coralloid roots were dominated mostly by the nitrogen-fixer Nostocaceae but also contain other non-diazotrophic bacteria. The sequencing of entire 16S rRNA gene identified four different strains of cyanobacteria under the heterocystous genera Nostoc and Desmonostoc. Meanwhile, the top bacterial families in normal roots were Xanthobacteraceae, Burkholderiaceae, and Bacillaceae. Moreover, a diverse fungal community was also found in the roots of cycads and the predominating families were Ophiocordycipitaceae, Nectriaceae, Bionectriaceae, and Trichocomaceae. Our results demonstrated that bacterial diversity in normal roots of C. fairylakea is higher in richness and abundance than C. debaoensis. On the other hand, a slight difference, albeit insignificant, was noted for the diversity of fungi among root types and host species as the number of shared taxa is relatively high (67%). Our results suggested that diverse microbes are present in roots of cycads which potentially interact together to support cycads survival. Our study provided additional knowledge on the microbial diversity and composition in cycads and thus expanding our current knowledge on cycad-microbe association. Our study also considered the possible impact of ex situ conservation on cyanobiont community of cycads.

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.

Which Species Should We Focus On? Umbrella Species Assessment in Southwest China

In conservation biology, umbrella species are often used as agents for a broader set of species, or as representatives of an ecosystem, and their conservation is expected to benefit a large number of naturally co-occurring species. Southwest China is home to not only global biodiversity hotspots, but also rapid economic and population growth and extensive changes in land use. However, because of the large regional span, the diverse species distributions, and the difficulty of field investigations, traditional methods used to assess umbrella species are not suitable for implementation in Southwest China. In the current study, we assessed 810 key protected species from seven taxa by indicator value analysis, correlation analysis, and factor analysis. We selected 32 species as umbrella species, whose habitats overlapped the habitats of 97% of the total species. Furthermore, the selected species were significantly correlated with 70% of all species in the study area. A total of 16 out of 19 selected animal species have been previously mentioned as umbrella species, compared with only 3 out of 13 plants species; this is despite plants accounting for a large proportion of the total species in Southwest China. We discuss the roles of indicator species and co-occurring species, and provide suggestions for species protection in Southwest China based on the current results. Our research provides valuable scientific information for research on umbrella conservation species over large geographical scales, and related fields of biodiversity conservation.

Investigating the Genetic Diversity, Population Differentiation and Population Dynamics of Cycas segmentifida (Cycadaceae) Endemic to Southwest China by Multiple Molecular Markers

Climate change, species dispersal ability and habitat fragmentation are major factors influencing species distribution and genetic diversity, especially for the range-restricted and threatened taxa. Here, using four sequences of chloroplast DNAs (cpDNAs), three nuclear genes (nDNAs) and 12 nuclear microsatellites (SSRs), we investigated the genetic diversity, genetic structure, divergence time and population dynamics of Cycas segmentifida D. Y. Wang and C. Y. Deng, a threatened cycad species endemic to Southwest China. High levels of genetic diversity and genetic differentiation were revealed in C. segmentifida. Haplotypes of networks showed two evolutionary units in C. segmentifida, with the exception of the nuclear gene GTP network. Meanwhile, the UPGMA tree, structure and PCoA analyses suggested that 14 populations of C. segmentifida were divided into two clades. There was significant effect of isolation by distance (IBD) in this species. However, this species did not display a significant phylogeographic structure. The divergence time estimation suggested that its haplotypes diverged during the Middle Pleistocene. Additionally, the population dynamics inferred from different DNA sequences analyses were discordant. Bottleneck analysis showed that populations of C. segmentifida did not experience any recent bottleneck effect, but rather pointed to a contraction of its effective population size over time. Furthermore, our results suggested that the population BM which held an intact population structure and occupied undisturbed habitat was at the Hardy-Weinberg equilibrium, implying that this population is a free-mating system. These genetic features provide important information for the sustainable management of C. segmentifida.

The distribution, diversity, and conservation status of Cycas in China

As ancient gymnosperm and woody plants, cycads have survived through dramatic tectonic activities, climate fluctuation, and environmental variations making them of great significance in studying the origin and evolution of flora biodiversity. However, they are among the most threatened plant groups in the world. The principal aim of this review is to outline the distribution, diversity, and conservation status of Cycas in China and provide suggestions for conservation practices. In this review, we describe the taxonomy, distribution, and conservation status of Cycas in China. By comparing Chinese Cycas species with its relatives worldwide, we then discuss the current genetic diversity, genetic differentiation of Cycas, and try to disentangle the potential effects of Quaternary climate changes and topographical events on Cycas. We review conservation practices from both researchers and practitioners for these rare and endangered species. High genetic diversity at the species level and strong genetic differentiation within Cycas have been observed. Most Cycas species in southwest China have experienced population retreats in contrast to the coastal Cycas's expansion during the Quaternary glaciation. Additionally, human activities and habitat fragmentation have pushed these endangered taxa to the brink of extinction. Although numerous efforts have been made to mitigate threats to Cycas survival, implementation and compliance monitoring in protection zones are currently inadequate. We outline six proposals to strengthen conservation measures for Cycas in China and anticipate that these measures will provide guidelines for further research on population genetics as well as conservation biology of not only cycads but also other endangered species worldwide.

Genetic structure and demographic history of Cycas chenii (Cycadaceae), an endangered species with extremely small populations

Geological activities and climate oscillations during the Quaternary period profoundly impacted the distribution of species in Southwest China. Some plant species may be harbored in refugia, such as the dry-hot valleys of Southwest China. Cycas chenii X. Gong & W. Zhou, a critically endangered cycad species, which grows under the canopy in subtropical evergreen broad-leaved forests along the upstream drainage area of the Red River, is endemic to this refugium. In this study, 60 individuals of C. chenii collected from six populations were analyzed by sequencing two chloroplast intergenic spacers (cpDNA: psbA-trnH and trnL-trnF) and two nuclear genes (PHYP and RBP-1). Results showed high genetic diversity at the species level, but low within-population genetic diversity and high interpopulation genetic differentiation. A Bayesian phylogenetic tree based on cpDNA showed that five chloroplast haplotypes were clustered into two clades, which corresponds to the division of the western and eastern bank of the Red River. These data indicate a possible role for the Red River as a geographic barrier to gene flow in C. chenii. Based on our findings, we propose appropriate in situ and ex situ conservation strategies for C. chenii.

Comparison of intraspecific, interspecific and intergeneric chloroplast diversity in Cycads

Cycads are among the most threatened plant species. Increasing the availability of genomic information by adding whole chloroplast data is a fundamental step in supporting phylogenetic studies and conservation efforts. Here, we assemble a dataset encompassing three taxonomic levels in cycads, including ten genera, three species in the genus Cycas and two individuals of C. debaoensis. Repeated sequences, SSRs and variations of the chloroplast were analyzed at the intraspecific, interspecific and intergeneric scale, and using our sequence data, we reconstruct a phylogenomic tree for cycads. The chloroplast was 162,094 bp in length, with 133 genes annotated, including 87 protein-coding, 37 tRNA and 8 rRNA genes. We found 7 repeated sequences and 39 SSRs. Seven loci showed promising levels of variations for application in DNA-barcoding. The chloroplast phylogeny confirmed the division of Cycadales in two suborders, each of them being monophyletic, revealing a contradiction with the current family circumscription and its evolution. Finally, 10 intraspecific SNPs were found. Our results showed that despite the extremely restricted distribution range of C. debaoensis, using complete chloroplast data is useful not only in intraspecific studies, but also to improve our understanding of cycad evolution and in defining conservation strategies for this emblematic group.

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.

Species delimitation, genetic diversity and population historical dynamics of Cycas diannanensis (Cycadaceae) occurring sympatrically in the Red River region of China

Delimitating species boundaries could be of critical importance when evaluating the species' evolving process and providing guidelines for conservation genetics. Here, species delimitation was carried out on three endemic and endangered Cycas species with resembling morphology and overlapped distribution range along the Red River (Yuanjiang) in China: Cycas diananensis Z. T. Guan et G. D. Tao, Cycas parvula S. L. Yang and Cycas multiovula D. Y. Wang. A total of 137 individuals from 15 populations were genotyped by using three chloroplastic (psbA-trnH, atpI-atpH, and trnL-rps4) and two single copy nuclear (RPB1 and SmHP) DNA sequences. Basing on the carefully morphological comparison and cladistic haplotype aggregation (CHA) analysis, we propose all the populations as one species, with the rest two incorporated into C. diannanensis. Genetic diversity and structure analysis of the conflated C. diannanensis revealed this species possessed a relative lower genetic diversity than estimates of other Cycas species. The higher genetic diversity among populations and relative lower genetic diversity within populations, as well as obvious genetic differentiation among populations inferred from chloroplastic DNA (cpDNA) suggested a recent genetic loss within this protected species. Additionally, a clear genetic structure of C. diannanensis corresponding with geography was detected based on cpDNA, dividing its population ranges into "Yuanjiang-Nanhun" basin and "Ejia-Jiepai" basin groups. Demographical history analyses based on combined cpDNA and one nuclear DNA (nDNA) SmHP both showed the population size of C. diannanensis began to decrease in Quaternary glaciation with no subsequent expansion, while another nDNA RPB1 revealed a more recent sudden expansion after long-term population size contraction, suggesting its probable bottleneck events in history. Our findings offer grounded views for clarifying species boundaries of C. diannanensis when determining the conservation objectives. For operational guidelines, the downstream populations which occupy high and peculiar haplotypes should be given prior in-situ conservation. In addition, ex-situ conservation and reintroduction measures for decades of generations are supplemented for improving the population size and genetic diversity of the endemic and endangered species.

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.

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.