Comprehensive comparative analysis and development of molecular markers for Lasianthus species based on complete chloroplast genome sequences

Lasianthus species are widely used in traditional Chinese folk medicine with high medicinal value. However, source materials and herbarium specimens are often misidentified due to morphological characteristics and commonly used DNA barcode fragments are not sufficient for accurately identifying Lasianthus species. To improve the molecular methods for distinguishing among Lasianthus species, we report the complete chloroplast (CP) genomes of Lasianthus attenuatus, Lasianthus henryi, Lasianthus hookeri, Lasianthus sikkimensis, obtained via high-throughput Illumina sequencing. These showed CP genomes size of 160164-160246 bp and a typical quadripartite structure, including a large single-copy region (86675–86848 bp), a small single-copy region (17177–17326 bp), and a pair of inverted repeats (28089–28135 bp). As a whole, the gene order, GC content and IR/SC boundary structure were remarkably similar among of the four Lasianthus CP genomes, the partial gene length and IR, LSC and SSC regions length are still different. The average GC content of the CP genomes was 36.71–36.75%, and a total of 129 genes were detected, including 83 different protein-coding genes, 8 different rRNA genes and 38 different tRNA genes. Furthermore, we compared our 4 complete CP genomes data with publicly available CP genome data from six other Lasianthus species, and we initially screened eleven highly variable region fragments were initially screened. We then evaluated the identification efficiency of eleven highly variable region fragments and 5 regular barcode fragments. Ultimately, we found that the optimal combination fragment' ITS2 + psaI-ycf4' could authenticated the Lasianthus species well. Additionally, the results of genome comparison of Rubiaceae species showed that the coding region is more conservative than the non-coding region, and the ycf1 gene shows the most significant variation. Finally, 49 species of CP genome sequences belonging to 16 genera of the Rubiaceae family were used to construct phylogenetic trees. Our research is the first to analyze the chloroplast genomes of four species of Lasianthus in detail and we ultimately determined that the combination fragment' ITS2 + psaI-ycf4' is the optimal barcode combination for identifying the genus of Lasianthus. Meanwhile, we gathered the available CP genome sequences from the Rubiaceae and used them to construct the most comprehensive phylogenetic tree for the Rubiaceae family. These investigations provide an important reference point for further studies in the species identification, genetic diversity, and phylogenetic analyses of Rubiaceae species.