Assembly and phylogenetic analysis of the mitochondrial genome of endangered medicinal plant Huperzia crispata.

Huperzia crispata is a traditional Chinese herb plant and has attracted special attention in recent years for its products Hup A can serve as an acetylcholinesterase inhibitor (AChEI). Although the chloroplast (cp) genome of H. crispata has been studied, there are no reports regarding the Huperzia mitochondrial (mt) genome since the previously reported H. squarrosa has been revised as Phlegmariurus squarrosus. The mt genome of H. crispata was sequenced using a combination of long-read nanopore and Illumina sequencing platforms. The entire H. crispata mt genome was assembled in a circular with a length of 412,594 bp and a total of 91 genes, including 45 tRNAs, 6 rRNAs, 37 protein-coding genes (PCGs), and 3 pseudogenes. Notably, the rps8 gene was present in P. squarrosus and a pseudogene rps8 was presented in H. crispata, which was lacking in most of Pteridophyta and Gymnospermae. Intron-encoded maturase (mat-atp9i85 and mat-cobi787) genes were present in H. crispata and P. squarrosus, but lost in other examined lycophytes, ferns, and Gymnospermae plants. Collinearity analysis showed that the mt genome of H. crispata and P. squarrossus is highly conservative compared to other ferns. Relative synonymous codon usage (RSCU) analysis showed that the amino acids most frequently found were phenylalanine (Phe) (4.77%), isoleucine (Ile) (4.71%), lysine (Lys) (4.26%), while arginine (Arg) (0.32%), and histidine (His) (0.42%) were rarely found. Simple sequence repeats (SSR) analysis revealed that a total of 114 SSRs were identified in the mt genome of H. crispata and account for 0.35% of the whole mt genome. Monomer repeats were the majority types of SSRs and represent 91.89% of the total SSRs. In addition, a total of 1948 interspersed repeats (158 forward, 147 palindromic, and 5 reverse repeats) with a length ranging from 30 bp to 14,945 bp were identified in the H. crispata mt genome and the 30-39-bp repeats were the most abundant type. Gene transfer analysis indicated that a total of 12 homologous fragments were discovered between the cp and mt genomes of H. crispata, accounting for 0.93% and 2.48% of the total cp and mt genomes, respectively. The phylogenetic trees revealed that H. crispata was the sister of P. squarrosus. The Ka/Ks analysis results suggested that most PCGs, except atp6 gene, were subject to purification selection during evolution. Our study provides extensive information on the features of the H. crispata mt genome and will help unravel evolutionary relationships, and molecular identification within lycophytes.

The complete chloroplast genome of the medical plant Huperzia crispata from the Huperziaceae family: structure, comparative analysis, and phylogenetic relationships.

BACKGROUND: Huperzia crispata, belonging to the Huperziaceae family, is one of the most essential resources of huperzine A for candidate drugs to treat Alzheimer's diseases. However, there is very limited information about H. crispat, and its taxonomic status and interspecific relationships between Huperzia species are still unclear. To investigate the taxonomic classification of Huperzia species and identify species discrimination markers, the complete chloroplast (cp) genome of H. crispata was sequenced and characterized for the first time. METHODS AND RESULTS: Total genomic DNA was isolated and sequenced using the next-generation Illumina NovaSeq 6000 platform. The data were filtered, assembled and annotated by a series software and web service. The results were as follows: the cp genome of H. crispata was 154,320 bp long with a large single-copy (LSC) region of 104,023 bp, a small single-copy (SSC) region of 19,671 bp, and a pair of inverted repeat (IRa and IRb) regions of 15,313 bp. A total of 131 genes, including 87 protein-coding genes, 36 transfer RNA genes (tRNAs), and eight ribosome RNA genes (rRNAs), were annotated in the cp genome. The contraction and expansion of the inverted repeat (IR) regions were relatively conserved in the Huperzia genus. Codon usage bias analysis showed that the encoding rate at the 3-end of codon A/T (74.34%) was significantly higher than that of C/G (25.66%). A total of 8 hotspot loci with high Pi values (> 0.06) were identified in the four Huperzia species based on nucleic acid diversity analysis. Ka/Ks selective pressure analysis demonstrated that the cemA gene is the most common gene undergoing positive selection among Huperzia. In addition, a total of 261 simple sequence repeats and 179 interspersed repeats were identified in the cp genome. Phylogenetic tree analysis based on the complete protein sequences of 23 related species of H. crispata indicated that H. serrata f. longipetiolata is a sister of H. crispata, suggesting that H. serrata f. longipetiolata and H. crispata are more closely related than H. serrata and H. lucidula. CONCLUSIONS: The results strongly supported that H. crispata was more closely related to H. serrata f. longipetiolata than to H. serrata and H. lucidula within the Huperzia genus. The outcome provided important information for the phylogenetic analysis of the subsequent specific molecular species identification in Huperzia. The present results will provide valuable information for further research into the classification, phylogeny and species identification of Huperzia plants.