Comparative Analysis of the Complete Chloroplast Genomes of Six Endangered Cycas Species: Genomic Features, Comparative Analysis, and Phylogenetic Implications

Cycas (family Cycadaceae), which spread throughout tropical and subtropical regions, is crucial in conservation biology. Due to subtle morphological variations between species, a solid species-level phylogeny for Cycas is lacking. In the present study, we assembled and analyzed the chloroplast genomes of six Cycas plants, including their genome structure, GC content, and nucleotide diversity. The Cycas chloroplast genome spans from 162,038 to 162,159 bp and contains 131 genes, including 86 protein-coding genes, 37 transfer RNA (tRNA) genes, and 8 ribosomal RNA (rRNA) genes. Through a comparative analysis, we found that the chloroplast genome of Cycas was highly conserved, as indicated by the contraction and expansion of the inverted repeat (IR) regions and sequence polymorphisms. In addition, several non-coding sites (psbK-psbI, petN-psbM, trnE-UUC-psbD, ndhC-trnM-CAU, and rpl32-trnP-GGG) showed significant variation. The utilization of phylogenetic analysis relying on protein-coding genes has substantiated the division of Cycas primarily into four groups. The application of these findings will prove valuable in evaluating genetic diversity and the phylogenetic connections among closely related species. Moreover, it will provide essential support for the advancement of wild germplasm resources.